THESTAR

FORMATION

NEWSLETTER

An electronic publication dedicated to early stellar/planetary evolution and molecular clouds
No. 319 — 17 July 2019

Editor: Bo Reipurth (reipurth@ifa.hawaii.edu)

List of Contents
The Star Formation Newsletter
Editor: Bo Reipurthreipurth@ifa.hawaii.edu

Interview . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
Perspective . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
Abstracts of Newly Accepted Papers . . . . . . . . . . 10

Associate Editor: Anna McLeod
anna.mcleod@berkeley.edu

Dissertation Abstracts . . . . . . . . . . . . . . . . . . . . . . . . 35

Technical Editor: Hsi-Wei Yen
hwyen@asiaa.sinica.edu.tw

Meetings . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 38

Editorial Board

Short Announcements . . . . . . . . . . . . . . . . . . . . . . . . 42

Joao Alves
Alan Boss
Jerome Bouvier
Lee Hartmann
Thomas Henning
Paul Ho
Jes Jorgensen
Charles J. Lada
Thijs Kouwenhoven
Michael R. Meyer
Ralph Pudritz
Luis Felipe Rodrı́guez
Ewine van Dishoeck
Hans Zinnecker
The Star Formation Newsletter is a vehicle forfast distribution of information of interest for astronomers working on star and planet formationand molecular clouds. You can submit material
for the following sections: Abstracts of recently
accepted papers (only for papers sent to refereed
journals), Abstracts of recently accepted major reviews (not standard conference contributions), Dissertation Abstracts (presenting abstracts of new
Ph.D dissertations), Meetings (announcing meetings broadly of interest to the star and planet formation and early solar system community), New
Jobs (advertising jobs specifically aimed towards
persons within the areas of the Newsletter), andShort Announcements (where you can inform or request information from the community). Additionally, the Newsletter brings short overview articles
on objects of special interest, physical processes ortheoretical results, the early solar system, as wellas occasional interviews.

Newsletter Archive
www.ifa.hawaii.edu/users/reipurth/newsletter.htm

New Jobs . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37

Summary of Upcoming Meetings . . . . . . . . . . . . . 41

Cover Picture
IC 2944 and IC 2948 form a large HII region in Centaurus, energized by a group of OB stars dominatedby the O7 III star HD 101205. Projected towards
the OB stars is a group of compact dark globules,
discovered by A.D. Thackeray in 1955. The globules are splintered and have convulsed shapes, aswell illustrated in this image taken with the ESO
VLT. The globule complex is likely in the processof being destroyed by the powerful UV radiation.
Image courtesy ESO.

Submitting your abstracts
Latex
macrosfor
submitting
abstractsand dissertation abstracts (by e-mail toreipurth@ifa.hawaii.edu)
areappendedtoeach Call for Abstracts.
You can also
submitvia
theNewsletterwebinterfaceat
http://www2.ifa.hawaii.edu/starformation/index.cfm

ject of your 1971 Annual Reviews article. What were theissues back then?

Carl Heiles

A: My thesis work covered an area about 160 square degrees, which happened to contain some dust clouds withsizes of a degree or so. We had the idea that the dust/gas
ratio was roughly constant, which meant that these clouds
should have contained enormous amounts of interstellar
gas. But the 21-cm line showed no excess intensity towards these clouds; if anything, it tended to be weaker!
This planted the idea that the gas in the clouds might bemolecular. When I got to Berkeley immediately after receiving my PhD, I was in the amazingly fortunate position
that Berkeley’s Hat Creek telescope was tuned to the OH
line. So it was easy pickings to see if OH resided in thoseclouds. And it did! Using the dust content as a surrogate
for the gas provided an estimate for the strength of selfgravity. And self-gravity turned out to be important! Sothese were likely areas of star formation. Connecting all
those dots was an exciting prospect for an Annual Reviews
article. History repeats itself: these days, the presence of‘dark gas’—large-density regions with no CO emission—isinferred in the same way.

in conversation with Bo Reipurth

Q: Your PhD thesis from 1966 dealt with observations of
the spatial structure of interstellar hydrogen. What motivated you to choose that topic?

Q: Your discovery of H I shells and supershells at highGalactic latitudes appeared in 1979, with a follow-up in
1984, and these papers have been very widely cited. Pleasetell how this discovery came about.

A: As a kid I was interested in technical things and becamea ‘ham’ radio operator. This began my interest in electronics and sparked my desire to be an electrical engineer,which was my original major at Cornell University in 1958.
Cornell had a program of public lectures on various science
topics. In my third year it focused on astronomy and had
four amazing lecturers: Tommy Gold, Fred Hoyle, Hermann Bondi, and Geoffrey Burbidge. They talked about
the great astronomical questions, ranging from Solar magnetohydrodynamics (MHD) to cosmology. I was hooked: Ichanged my major to Engineering Physics and took every
course related to those areas that I could. Marshall Cohentaught Plasma physics; he went on to CalTech and specialized in VLBI. I took two courses on stellar structure: theconventional version with the differential equations, beautifully taught by John Cox; in contrast, Ed Salpeter eschewed the equations and emphasized phenomenology andorder-of-magnitude estimates. Those two courses opened
my eyes: you achieve deeper understanding by looking at
problems from multiple viewpoints.

A: It was totally obvious that supershells would be detectable in the Galactic plane! Harm Habing and I had
recently finished the Hat Creek 85-foot telescope survey ofhigh-latitude 21-cm line emission—half a million H i line
profiles! We pioneered using grayscale and color imaging
to look at them. What jumped out of those first images
were the Orion-Eridanus superbubble and the Sco-Oph supershell. They were tens of degrees in diameter and were
nearby. Given the Copernican Principle, the Galaxy had
to be peppered with such objects. At 1 Kpc distance, theangular diameter would be about 6 degrees. So I lookedat Harold Weaver’s Hat Creek H i survey of the Galactic
plane—and there they were.
Q: A few years later you calculated the rate at which clustered Type II supernovae produce bubbles that break throughthe dense, thin layer of disk gas. How is the correspondence with current observations?

All this, plus being a grad student at Princeton underGeorge Field and Lyman Spitzer, made the interstellar
medium a hugely attractive research area. When the timecame for my thesis work, the 300-foot telescope had justbeen finished and was superbly equipped for 21-cm line
work with a 100-channel 1-bit digital correlator. There I
was, a lowly graduate student, empowered with the besttelescope/electronics in the world for this project. Whowouldn’t jump in that direction?

A: That’s an excellent question because such events are responsible for transferring mass from the disk to the halo.
This was explicitly shown by Colin Norman and Satoru
Ikeuchi, who in 1989 wrote my favorite paper on thissubject and invented the compelling moniker ‘interstellar
chimneys’. Also, infalling pristine intergalactic gas clouds
have to transit through the ambient halo gas. These days
there’s a lot of activity pursuing these issues at high redQ: Shortly after your PhD you became interested in the shift, because they directly affect galactic evolution, butphysical conditions inside dust clouds, which was the sub3

it’s my impression that more work on synthesizing multi- A: The agglomeration was an update of the original one
ple data sets for local galaxies would reveal lots of impor- by Mathewson and Ford; their spectacular result was thetant details.
amazingly coherent pattern of linear polarization aligned
Q: You were early on interested in magnetic fields in clouds, with Radio Loop I (the North Polar Spur). The new agand summarized what was known in an Annual Reviews glomeration, together with the then-new all-sky H i maps,
article in 1976. You have maintained this interest all your showed that the morphological similarity was essentiallylife, so you must have a unique perspective on how the sub- universal and applied to structure on all scales. Thesedays, Susan Clark has used polarized dust emission fromject has developed over the years.
Planck together with the GALFA H i survey data to find
A: My interest in magnetic fields stems from those public highly elongated H i ‘fibers’ that are aligned with the neighlectures at Cornell, where they were reinforced by a spec- boring magnetic field. This is proving to be a rich subject!
tacular display of aurorae—rare at Cornell’s magnetic latitude. In many astronomical systems, magnetic pressure Q: In the early 2000s you embarked, together with Thomas
is comparable to other pressures, so the field has to be Troland, on the Millennium Arecibo 21 cm Absorptionimportant. But specifying the influence of the magnetic Line Survey. What were the goals and achievements?
field as a single number misses the fundamental point that A: Tom was the only grad student crazy enough to work
the magnetic forces are anisotropic, meaning that not only with me on interstellar magnetic fields, specifically meapressure but also magnetic tension plays a role, particu- suring Zeeman splitting of the 21-cm line seen in emission.
larly because of the coupling of cosmic rays to the field and Why ‘crazy’ ? The detected signatures were weak, requirits frozen-in gas. It’s these aspects that produce unique ing huge integration times, and the instrumental effects
and important dynamics, such as sunspots and angular (especially polarized beam response) comparable to themomentum transport.
actual splitting. Fortunately, Tom is very methodical andRegarding my ‘unique perspective’, there’s the scientific careful, compensating for my deficiencies in those areas.
one and the human one. I think the more important, right He got results for his successful thesis and moved on tonow, is the human one, because technology has advanced study fields in molecular clouds, often with Dick Crutcher.
to enable people to study magnetic phenomenae, both Meanwhile, we improved the Hat Creek instrumentation
observationally (think optical CCDs and radio wideband and measured H i Zeeman splitting for hundreds of poFaraday rotation maps) and theoretically (think large- sitions. Those measurements were controversial because
scale computers). So we need smart and clever people of the instrumental effects, with criticisms especially by
to exploit these opportunities. There’s a big future here! Gerrit Verschuur (who made the first real measurementsof Zeeman splitting in the ISM).
Q: The tiny-scale atomic structure of the ISM was the
subject of a study you published in 1997. What techniques The obvious way around polarized beam response problems was to forgo emission lines and instead use absorption
did you use, and what were your conclusions?
lines. The best telescope for that was Arecibo. Tom and I
A: It is straightforward to show that tiny-scale atomic had kept in collaborative and social contact, and this idea
structure (TSAS) cannot exist: where the H i column den- popped up. It was a good project with useful results, notsity changes over small angles, the volume density must be only for magnetic fields, but also statistics on nonthermallarge (scales of tens of AU), and this means high pressure, line broadening (turbulence), cloud temperature, ISM fillso the region will explode—and annihilate itself. The only ing factors. . . Some highly-motivated crazy person needs
problem: TSAS is commonly observed! Two ways out: (1) to expand this work and attain a decent statistical sample!
the structures are anisotropic, so lengths along the lineof-sight are longer than those across and volume densities Q: What are your current interests?
don’t need to be so high; (2) the rapid changes with an- A: With retirement, I resolved to become more consciengle are simply a manifestation of the spatial spectrum of tious with my collaborative partners. But I’m a sucker forinterstellar turbulence (proposed by Avinash Deshpande). new and interesting projects, and I’m busier than ever on
Almost certainly both occur, but we don’t know which pre- a number of seemingly disparate projects. These include
dominates. Similar statements apply to the ionized struc- time variations in OH/IR star masers (an outgrowth oftures, which produce interstellar scintillation and scatter- radio results on Solar-orbiting spacecraft(!)), linear polaring of pulsars and radio-loud quasars. These fascinating ization of the 21-cm line, ‘dark gas’, magnetic fields, H itiny structures are not receiving the attention they de- morphology and dynamics. I’m also working with Chinese
serve, making this a field ripe for future studies!
radio astronomers on FAST—both instrumentally, with
Q: In 2000 you published the results of the herculean task calibration, and scientifically, working with grad students.
to agglomerate all existing stellar polarization catalogs, to- All of this activity continues to make my scientific life
stimulating, productive, and fulfilling.
taling 9286 stars. What did you learn?

4

!"

#"

Perspective
Planet formation in clusters
Susanne Pfalzner

Figure 1: Environmental effects on protoplanetary discs.
a) HST image of an externally photo-evaporated disc(credit: NASA/ESA) and b) simulation result of the effectof a stellar flyby.
For both questions still no definite answer exists, but there
has been considerable progress during past years.
The central parameter is the “stellar density” since it determines both the frequency and the strength of the external influence. This applies equally to radiation as togravitational interactions. Sparse stellar groups like Cha I
consist of just a few dozens of stars, whereas dense groups
like NGC 3603 and Trumpler 14 can contain tens of thousands of stars packed within less than 1 pc3 . Even in
high-mass stellar groups (Mc > 103 M⊙ ) the average stelvary over many orders of magnitude from <In the context of planet formation this question arises to: lar densities
−3
0.1
M
pc
to > 105 M⊙ pc−3 (Wolff et al. 2007).
⊙is it sufficient to study the nascent planetary system in isolation? Stars usually do not form in isolation but as part Even without large simulations it is intuitively obviousof a stellar group (Lada & Lada 2003, Porras 2003). The that environmental effects will play a much larger role infirst important question in this field is then: How impor- Trumpler 14 than in Cha I. This is why investigations totant is the influence of the surrounding stars on circum- date have tended to concentrate on fairly dense systems
stellar discs and forming planetary systems? Distance to like the Orion Nebula Cluster (ONC) due to its relatively
other stars in the same system is the key factor here, so proximity to us (383 ± 3 pc; Kounkel et al. 2017). For
the main parameter is obviously the local stellar density. this reason the ONC is often taken as model cluster for
There are basically two main channels of influence exerted theoretical investigations, too (for example, Adams 2010,by the surrounding stars on discs: their radiation and Wijnen et al. 2017, Winter et al. 2018, Cuello et al. 2019,gravitational interactions. Both effects might truncate or Portegies Zwart 2019). The general outcome of these inpossibly even completely destroy a planet-forming disc. vestigations is that the environment is only responsible
Famous examples of the influence of radiation on discs are for the destruction of 5-10% of all discs. Some small discthe so-called proplyds in Orion, which show the destruc- sizes can be attributed to the environment, but typicallytive power of external photo-evaporation on young discs the disc size in the ONC is only reduced to about 100 AU.
(Fig. 1a). The signs of gravitational truncation of discs are This is only slightly smaller than the typical disc sizes ofmore difficult to track because close flybys mainly occur in 100-200 AU observed in sparse associations.
One well-tested method in science is to separate the object of interest from its surroundings and look at it inisolation. The advantage is that unimportant informationis removed and the true properties of the object are seen
more clearly. However, sometimes the influences of thesurroundings actually determine the properties of an object. In this case, not taking the environment into accountcan lead to incomplete or even false conclusions.

the strongly embedded, and therefore difficult to observe,
phase of young clusters. Tell-tale signs are spiral arms
(Fig. 1b). However, they can also be caused by other processes and tend to evolve into ring-like structures withinjust a few thousand years (Cuello et al. 2019). Thus, we
can pose a second important question: Does radiation or
gravitational interaction dominate?

However, the density in a stellar group is far from constant, and changes by several orders of magnitude withinjust a few Myr. In my opinion, the key here lies in understanding early cluster dynamics. This means that theresult obtained for the ONC is far from universal: it is
only valid for this particular stellar group at this specificphase of its temporal development. If one is interested in5

1
Upp Cen-Crux
I Lac 1
Low Cen-Crux
Ori 1a1

RSGC01

Westerlund 1
Westerlund 2

0,4

NGC 3603
Arches

0

1

10

0

2

¡ Cha

CrA

0,2

MBM 12

Taurusclusters

Trumpler 14

TWA Hya

χ Per
h Per
DSB2003
Quintuplett

0,6

d Cha

RSGC02dense clusters
associations
sparse associations
Haisch (2003)
Mamajek(2009)

Trumpler 16

reff [pc]

RSGC03

TTC

Ori 1cIC 1805
NGC 7380
NGC 6611
Ori 1bCygOB2
NGC 2244

Trumpler 14
Arches0,8Disc fraction

10

Upp ScoNGC 3603

associations

4

6

8

10

12

Cluster age [Myr]Age [Myr]

Figure 2: Two types of stellar groups — compact clus- Figure 3: Disc fraction in compact clusters (red), associters (red) and loose associations (blue, figure adapted from ations (blue) and co-moving groups (green), Pfalzner et
al. 2014.
Pfalzner 2009).
Again, observations are snapshots in time, yet young stellar groups are highly dynamical objects whose density isconstantly changing (Concha-Ramı́rez et al. 2019). Hence,
understanding the dynamics of young clusters is the key
Fortunately, we are spared the mammoth task of consid- to answering the question of how strong the influence ofering every observed stellar group individually, because it the surroundings actually is.
turns out there seem to be only two types of stellar groups:
For a long time it was unclear whether stellar groups
those that develop into long-lived clusters and those thatformed by merging of subgroups, or just appeared thelargely dissolve within about 5-10 Myr. These two types
way we observe them, and whether they expand as soon
are often referred to as clusters and associations, respecas they have formed. Recently the advent of Gaia results
1tively . Clusters are much more efficient in converting gas
clarified this situation. At ages 1–3 Myr, stellar groups
and dust into stars which allows them to form much more
typically expand with a velocity of 0.5–1 pc/Myr (Kuhncompact entities: at any given age their stellar density is
et al. 2019) and remain in the very dense phase 3 Myr
at least 100× higher than that of an association of theat most, and probably more likely only 1-2 Myr. Thissame age and mass. However, it should be stressed thatmeans that the stellar groups with little-to-no gas have inabout 80-90% of the field population was formed in shortfact been much denser in the past, when the influence of
lived associations and only 10-20% in long-lived clusters
the environment has also been stronger. Our simulations(Adamo et al. 2011, Bastian 2013).
(Vincke & Pfalzner 2018) have shown that as soon as the
Associations and clusters develop along different, but well- gas is expelled and the stellar groups start to expand thedefined evolutionary trajectories in terms of cluster size as influence of the environment drops considerably.
a function of time (see Fig. 2). This means that simulaGiven the difference in density, environmental effects aretions conforming to either of these temporal paths are notgenerally much stronger in clusters than in associations.
only valid for one specific stellar group at a given point of
Two properties can be tested in this respect by simulatime, but for an entire class of stellar groups over the fulltions and observations alike: disc frequency and disc size.
time period considered — a huge reduction in the possible
Roughly speaking it can be said that simulations show thatparameter space.
complete disc destruction rarely happens in associations
1 Historically it was thought that clusters formed as bound andthat contain only a few hundred stars; only in clusters
associations as unbound systems. Nowadays we know that most that contain a few thousand or ten thousand stars does
stellar groups are bound as long as they are still embedded in gas.
disc destruction become an issue, at least near the cluster
The historical misinterpretation leads to the confusing nomenclaturethat many embedded stellar groups are labelled clusters, even though centre. However, even in massive associations probably
they will largely dissolve as soon as they lose their gas like most of less than 20% of discs can be completely destroyed by thethe stellar groups in the solar neighbourhood.
environment (Olzcak et al. 2006, Vincke & Pfalzner 2016).
the environmental influence on discs in general, one has tolook at different stellar groups and also consider the timedependence of the stellar density.

6

onto eccentric orbits, which is not only the case for discmatter, but also for already formed planets. As most flybys will have some angle to the plane of the planets, thisleads to inclined orbits. Pr 0211 in M44 is possibly a firstexample for such a system (Pfalzner et al. 2018).
The values given above are averages over entire stellar
groups. However, it is the most central areas where thedensity is highest, and thus where disc destruction happens most frequently and disc fractions are lower thanelsewhere. As these stellar groups expand so do these innerareas. As observations mainly focus on the similar-sized
regions independent of cluster age, the cluster expansion
might lead to a faster drop in disc fraction than really
Figure 4: Average disc sizes expected from simulations of happens (Pfalzner et al. 2014)
different environments at ages 3 Myr and older. Summary There has been a long debate about whether externalof the results by Vincke & Pfalzner 2016, 2018
photo-evaporation or gravitational interactions dominate

Total disc destruction is of course an extreme and rather
negative outcome as far as forming planetary systems areconcerned! What is much more common is that the environment leads to disc truncation or redistribution of discmaterial. Just considering gravitational interactions, thehigh flyby frequency in stellar groups like the ONC could
be the reason that discs are rarely larger than 100-200 AU
in such environments (Vincke 2016).

in star forming regions (recently, Champion et al. 2017,
Haworth et al. 2018, Winter et al. 2018). To me, the answer is just that “it depends”. Fig. 5 illustrates this in aqualitative way. Both effects require high stellar densities
to be efficient. However, for external photo-evaporation towork it is also necessary for the radiation to penetrate efficiently through the cluster/association. This means thatthe gas/dust density in the stellar group has to be low
enough not to absorb the radiation. The deeply embedded phase lasts less than 1 Myr and so far very few observational constraints exist. However, during this earlyphase it is gravitational interactions that set the scene.
Afterwards the efficiency of close stellar flybys steadilydecreases, because as soon as most discs are truncated tothe disc size that is typical for this stellar group, events
that truncate the disc even further become seldom. The
frequency of events that change the disc size is reduced further when the stellar density eventually drops due to cluster expansion. The situation might be different in massive
long-lived clusters, but very few observational constraints
exist so far.

Unfortunately most clusters are further away from us thanassociations, so that there is little information about thedisc sizes in such environments. Fig. 4 gives a schematic
illustration of the average disc sizes that are expected indifferent environments from simulations. In the typicalstellar groups found in the solar neighbourhood containing
just a few hundred stars, the disc size is probably hardly
affected by the environment. By contrast, in dense longlived clusters, the mean disc size should be of the orderof just 10-20 AU in environments like Westerlund 2 andTrumpler 14. This means that planetary systems in clusters should be much more compact than those that formed
around field stars (Fujii & Hori 2019). In addition, it canbe expected that the outermost planets of these planetarysystems are often on eccentric and inclined orbits. As illustrated in Fig. 1b, any close stellar flybys excites matter

In associations external photo-evaporation (Fig. 5, blackline) basically inherits the disc frequency and disc size distribution that close flybys have created and potentiallyleads to additional disc truncation or even destruction.
The proplyds in the ONC are a prime example of such asituation. They are so clearly visible because the stellar
group is just at the stage when it has become gas-free inthe centre, while at the same time the stellar group has not
yet started to expand. This is the moment when external
photo-evaporation is at its strongest. However, this phaseis relatively short-lived, because stellar groups are extremely dynamical. How short depends somewhat on themass of the association: heavier ones develop faster thanlighter ones. Most of the environmental influence happens approximately during the first 2-3 Myr. This means
the time-window for efficient external photo-evaporation

Conversely, in clusters complete disc destruction by the environment is a real factor. A strong observational indicator is the much lower disc frequency in clusters compared
to associations (see Fig. 3) at a given age. This seems toindicate that for very dense massive clusters environmental effects might be the dominant disc destruction process.
However, one has to be careful with such a judgement as
we have no observations covering the formation phase ofsuch dense clusters. Therefore we do not know whether
these clusters start out with basically all stars being surrounded by discs or not.

7

If there are really 1015 of such objects floating around ininterstellar space, as anticipated, some of them could become part of a newly formed circumstellar disc. Due totheir size they could help to jump the meter-sized barrier. This way it is possible that these small objects arevital agents in seeding or accelerating planetary growth(Pfalzner & Bannister 2019).

Figure 5: Schematics of the relative importance of gravitational interactions versus external photo-evaporation during the different phases of cluster development.
is only 1-2 Myr in associations. Afterwards, the massive
stars dominating the radiation are just too distant to makeexternal photo-evaporation work efficiently.
By contrast, although the frequency of close flybys that
actually have an effect on the discs decreases dramatically, they never completely vanish. This means even in
the planetary system phase such interactions can also influence the newly formed planetary systems (Malmberg etal. 2007, Hao et al. 2013, Wang et al. 2015, Cai et al. 2017).
This is especially true in the most massive clusters like forexample Westerlund 1. For a long time it was thoughtthat such dense environments would hinder the formationof planets or at least reduce it considerably. However,in recent years over twenty planets have been found inlong-lived clusters (Pfalzner et al. 2018), so that at least
for short-period planets the environment does not prevent
their formation. However, these systems will have likely
more planets on highly eccentric orbits and fewer distantplanets; something that future exoplanet surveys will test.

[]

The influence of the environment is not restricted to distant stellar groups: our own solar system was once shapedby its environment. Several properties like the sharp outeredge at 30 AU or the composition of meteorites (Adams etal. 2006, Adams 2010, Parker 2017, Portegies Zwart 2019)
seem to bear the hallmarks of the past external influences.
Even the orbits of the hundreds of thousands of small objects moving outside Neptune’s orbit might be due to aclose stellar flyby in the Sun’s past (Pfalzner et al. 2018).

[]

Figure 6: Range of clustered star formation, stretchingfrom (a) low-mass associations like σ Ori, where the environment plays a minor role, to (b) high-mass clusters like
Westerlund 2 where disc sizes are dominated by the enviThe discovery of Oumuamua in 2017 (Meech et al. 2017) ronment. Image credits: ESO and Digitized Sky Survey 2;
has brought a dramatic new aspect to the interplay be- NASA, ESA, the Hubble Heritage Team (STScI/AURA),
tween discs and their environment. Namely, what hap- A. Nota (ESA/STScI), and the Westerlund 2 Science
pens to all the material that becomes unbound due to the Team.
interactions with the stellar cluster (Hands et al. 2019)?

8

In summary, I think it is important to see planetary systems not necessarily as isolated objects, but also to takethe interactions with their environment into account. Star
formation in the Milky Way mainly takes place in the
short-lived associations, like the ones surrounding us. Herethe influence on discs and forming planetary systems is
non-negligible but moderate (Fig. 6a). By contrast, in
long-lived clusters the environment plays a dominant role(Fig. 6b). Another important point is that even if specific clusters do not show signs of environmental influence
going on at the moment, the situation might have beenvery different in their past. The relative importance ofgravitational- vs. radiation-driven impact is a question ofdevelopmental phase. Gravitational interactions are mostimportant during the strongly embedded phase, they set
the initial conditions in terms of disc sizes and frequencies which can be modified by external photo-evaporation.
However, there is only a relatively short time available forradiation to act.

References:
Adamo, A., Östlin, G., & Zackrisson, E. 2011, MNRAS, 417, 1904
Adams, F. C., Proszkow, E. M., et al. 2006, ApJ, 641, 504
Adams, F. C. 2010, ARAA, 48, 47
Bastian, N. 2013, 370 Years of Astronomy in Utrecht, 287
Cai, M. X., et al. 2017, MNRAS, 470, 4337
Champion, J., Berné, O., Vicente, S., et al. 2017, A&A, 604, A69
Concha-Ramı́rez, F., Vaher, E., & Portegies Zwart, S. 2019,
MNRAS, 482, 732
Cuello, N., Dipierro, G., et al. 2019, MNRAS, 483, 4114
Facchini, S., Clarke, C. J., & Bisbas, T. G. 2016, MNRAS, 457,
3593
Farias, J. P., Tan, J. C., & Chatterjee, S. 2019, MNRAS, 483, 4999
Fujii, M. S., & Hori, Y. 2019, A&A, 624, A110
Hao, W., Kouwenhoven, M. B. N., & Spurzem, R. 2013, MNRAS,
433, 867
Hands, T. O., Dehnen, W., et al. 2019, MNRAS, 1064
Haworth, T. J., Clarke, C. J., et al. 2018, MNRAS, 481, 452
Kounkel, M., Hartmann, L., Loinard, L., et al. 2017, ApJ, 834, 142
Kuhn, M. A., Hillenbrand, L. A., et al. 2019, ApJ, 870, 32
Lada, C. J., & Lada, E. A. 2003, ARAA, 41, 57
Malmberg, D., de Angeli, F., Davies, M. B., et al. 2007, MNRAS,
378, 1207
Meech, K. J., Weryk, R., Micheli, M., et al. 2017, Nature, 552, 378
Olczak, C., Kaczmarek, T., et al. 2012, ApJ, 756, 123
Olczak, C., Pfalzner, S., & Spurzem, R. 2006, ApJ, 642, 1140
Parker, R. J. 2017, Handbook of Supernovae, 2313
Pfalzner, S. 2009, A&A, 498, L37
Pfalzner, S., Steinhausen, M., & Menten, K. 2014, ApJ, 793, L34
Pfalzner, S., & Bannister, M. T. 2019, ApJ, 874, L34
Pfalzner, S., Bhandare, A., et al. 2018, ApJ, 863, 45
Pfalzner, S., Bhandare, A., & Vincke, K. 2018, A&A, 610, A33
Porras, A., Christopher, M., et al. 2003, AJ, 126, 1916
Portegies Zwart, S. 2019, A&A, 622, A69
Richert, A. J. W., Getman, K. V., et al. 2018, MNRAS, 477, 5191
Scally, A., & Clarke, C. 2001, MNRAS, 325, 449
van Elteren, A., Portegies Zwart, S., et al. 2019, arXiv:1902.04652
Vincke, K., & Pfalzner, S. 2016, ApJ, 828, 48
Vincke, K., & Pfalzner, S. 2018, ApJ, 868, 1
Wang, L., Kouwenhoven, M. B. N., et al. 2015, MNRAS, 449, 3543
Wijnen, T. P. G., Pols, O. R., et al. 2017, A&A, 604, A91
Winter, A. J., Clarke, C. J., et al. 2018, MNRAS, 478, 2700
Wolff, S. C., Strom, S. E., Dror, D., et al. 2007, AJ, 133, 1092

What next? As pointed out above, cluster dynamics is a
key factor here. The wealth of the Gaia data will hopefully
enable us to put even tighter constraints on the dynamics of young clusters and associations alike. Equally, it
would be important to obtain better information on discs
in the embedded phase and in dense clusters, which is admittedly an observational challenge. On the theoretical
front, it is important to take the newly found observational constraints into account and eventually bridge thestill existing gap between cluster formation models and thefollowing phase of cluster expansion (Farias et al. 2019).
Moreover, so far it is rarely taken into account that suchclusters contain a large fraction of binaries and multiple
systems, and future studies should include that. Given theconsiderable progress during the last few years, I think allthese aims are reachable in the near future.

9

Abstracts of recently accepted papers

1

The structure of the Orion Nebula in the direction of Theta 1 Ori C
Nicholas Abel1 , Gary Ferland2 and Bob O’Dell3
1

University of Cincinnati, Clermont College, 4200 Clermont College Drive, Batavia, OH, 45103, USA; 2 University of
Kentucky, 505 Rose St, Lexington, KY 40508, USA; 3 Vanderbilt University, Box 1807-B, Nashville, TN 37235, USA
E-mail contact: npabel2 at gmail.com
We have used existing optical emission and absorption lines, [C II] emission lines, and H I absorption lines to create
a new model for a Central Column of material near the Trapezium region of the Orion Nebula. This was necessary
because recent high spectral resolution spectra of optical emission lines and imaging spectra in the [C II] 158 micronline have shown that there are new velocity systems associated with the foreground Veil and the material lying between
Theta 1 Ori C and the Main Ionization Front of the nebula. When a family of models generated with the spectral
synthesis code Cloudy were compared with the surface brightness of the emission lines and strengths of the Veilabsorption lines seen in the Trapezium stars, distances from Theta 1 Ori C, were derived, with the closest, highest
ionization layer being 1.3 pc. The line of sight distance of this layer is comparable with the size of the inner Huygens
Region in the plane of the sky. These layers are all blueshifted with respect to the Orion Nebula Cluster of stars,
probably because of the pressure of a hot central bubble created by Theta 1 Ori C’s stellar wind. We find velocity
components that are ascribed to both sides of this bubble. Our analysis shows that the foreground [C II] 158 micron
emission is part of a previously identified layer that forms a portion of a recently discovered expanding shell of material
covering most of the larger Extended Orion Nebula.
Accepted by ApJ
http://arxiv.org/pdf/1906.07779

2

Revealing the chemical structure of the Class I disc Oph-IRS 67
E. Artur de la Villarmois1, L. E. Kristensen1 and J. K. Jørgensen1
1

Niels Bohr Institute & Centre for Star and Planet Formation, University of Copenhagen, Øster Voldgade 5–7, 1350
Copenhagen K., DenmarkE-mail contact: elizabeth.artur at nbi.ku.dk

Recent results suggest that the first steps towards planet formation may be already taking place in protoplanetary discs
during the first 100,000 years after stars form. It is therefore crucial to unravel the physical and chemical structuresof such discs in their earliest stages while they are still embedded in their natal envelopes and compare them withmore evolved systems. The purpose of this paper is to explore the structure of a line-rich Class I protobinary source,
Oph-IRS 67, and analyse the differences and similarities with Class 0 and Class II sources. We present a systematic
molecular line study of IRS 67 with the Submillimeter Array (SMA) on 1–2 arcsec (150–300 AU) scales. The wide
instantaneous band-width of the SMA observations (∼30 GHz) provide detections of a range of molecular transitions
that trace different physics, such as CO isotopologues, sulphur-bearing species, deuterated species, and carbon-chain
molecules. We see significant differences between different groups of species. For example, the CO isotopologues andsulphur-bearing species show a rotational profile and are tracing the larger-scale circumbinary disc structure, while
CN, DCN, and carbon-chain molecules peak at the southern edge of the disc at blue-shifted velocities. In addition,
the cold gas tracer DCO+ is seen beyond the extent of the circumbinary disc. The detected molecular transitions canbe grouped into three main components: cold regions far from the system, the circumbinary disc, and a UV-irradiatedregion likely associated with the surface layers of the disc that are reached by the UV radiation from the sources. The
different components are consistent with the temperature structure derived from the ratio of two H2 CO transitions,that is, warm temperatures are seen towards the outflow direction, lukewarm temperatures are associated with the
UV-radiated region, and cold temperatures are related with the circumbinary disc structure. The chemistry towards
10

IRS 67 shares similarities with both Class 0 and Class II sources, possibly due to the high gas column density andthe strong UV radiation arising from the binary system. IRS 67 is, therefore, highlighting the intermediate chemistry
between deeply embedded sources and T-Tauri discs.
Accepted by Astronomy & Astrophysics
http://arxiv.org/pdf/1906.00685

3

The Effect of Binarity on Circumstellar Disk Evolution
Scott A. Barenfeld1 , John M. Carpenter2 , Anneila I. Sargent1 , Aaron C. Rizzuto3 , Adam L. Kraus3 ,Tiffany Meshkat4,5 , Rachel L. Akeson6 , Eric L. N. Jensen7 and Sasha Hinkley8
1

California Institute of Technology, Department of Astronomy, MC 249-17, Pasadena, CA 91125, USA; 2 Joint ALMA
Observatory, Av. Alonso de Cordova 3107, Vitacura, Santiago, Chile; 3 The University of Texas at Austin, Departmentof Astronomy, Austin, TX 78712, USA; 4 IPAC, Caltech, M/C 100-22, 1200 East California Boulevard, Pasadena,
CA 91125, USA; 5 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena,
CA 91109, USA; 6 IPAC-NExScI, Caltech, Pasadena, CA 91125, USA; 7 Swarthmore College, Dept. of Physics &Astronomy, 500 College Ave., Swarthmore, PA 19081-1390, USA; 8 University of Exeter, Physics Department, Stocker
Road, Exeter, EX4 4QL, UK
E-mail contact: scott.barenfeld at gmail.com
We present new results on how the presence of stellar companions affects disk evolution based on a study of the 5-11
Myr old Upper Scorpius OB Association. Of the 50 G0-M3 Upper Sco members with disks in our sample, only sevenhost a stellar companion within 2” and brighter than K = 15, compared to 35 of 75 members without disks. This
matches a trend seen in the 1-2 Myr old Taurus region, where systems with a stellar companion within 40 au havea lower fraction of infrared-identified disks than those without such companions, indicating shorter disk lifetimes inclose multiple systems. However, the fractions of disk systems with a stellar companion within 40 au match in Upper
Sco and Taurus. Additionally, we see no difference in the millimeter brightnesses of disks in Upper Sco systems withand without companions, in contrast to Taurus where systems with a companion within 300 au are significantly fainterthan wider and single systems. These results suggest that the effects of stellar companions on disk lifetimes occur
within the first 1-2 Myr of disk evolution, after which companions play little further role. By contrast, disks aroundsingle stars lose the millimeter-sized dust grains in their outer regions between ages of 1-2 Myr and 5-11 Myr. The
end result of small dust disk sizes and faint millimeter luminosities is the same whether the disk has been truncatedby a companion or has evolved through internal processes.
Accepted by Astrophysical Journal

http://arxiv.org/pdf/1906.04878

4

Variability of young stellar objects in the star-forming region Pelican Nebula
A. Bhardwaj1 , N. Panwar2 , G.J. Herczeg1 , W.P. Chen3 , and H.P. Singh4
1

Kavli Institute for Astronomy and Astrophysics, Peking University, Yi He Yuan Lu 5, Hai Dian District, Beijing
100871, China; 2 Aryabhatta Research Institute of Observational Sciences, Manora Peak, Nainital-263002, Uttarakhand, India; 3 Graduate Institute of Astronomy, National Central University, Jhongli 32001, Taiwan; 4 Departmentof Physics and Astrophysics, University of Delhi, Delhi-110007, India
E-mail contact: anupam.bhardwajj at gmail.com
We observed a field of 16′ × 16′ in the star-forming region Pelican Nebula (IC 5070) at BVRI wavelengths for 90 nights
spread over one year in 2012–2013. More than 250 epochs in VRI-bands are used to identify and classify variables
up to V ∼ 21 mag. We present a catalogue of optical time-series photometry with periods, mean-magnitudes andclassifications for 95 variable stars including 67 pre-main-sequence variables towards star-forming region IC 5070. The
pre-main-sequence variables are further classified as candidate classical T Tauri and weak-line T Tauri stars basedon their light curve variations and the locations on the color-color and color-magnitude diagrams using optical andinfrared data together with Gaia DR2 astrometry. Classical T Tauri stars display variability amplitudes up to threetimes the maximum fluctuation in disk-free weak-line T Tauri stars, which show strong periodic variations. Short-term
variability is missed in our photometry within single nights. Several classical T Tauri stars display long-lasting (≥10
11

days) single or multiple fading and brightening events up to a couple of magnitudes at optical wavelengths. The
typical mass and age of the pre-main-sequence variables from the isochrone-fitting and spectral energy distributions
are estimated to be ≤1 M⊙ and ∼2 Myr, respectively. We do not find any correlation between the optical amplitudes
or periods with the physical parameters (mass and age) of pre-main-sequence stars.
Accepted by A&A
http://arxiv.org/pdf/1906.00256

5

The complex chemistry of hot cores in Sagittarius B2(N): Influence of cosmic-ray ionization and thermal history
M. Bonfand1 , A. Belloche1 , R. T. Garrod2,3 , K. M. Menten1 , E. Willis3, G. Stphan3 and H. S. P. Mller4
1

Max-Planck-Institut fr Radioastronomie, Auf dem Hgel 69, Bonn, Germany; 2 Departments of Chemistry, Universityof Virginia, Charlottesville, VA 22904, USA; 3 Departments of Astronomy, University of Virginia, Charlottesville, VA
22904, USA; 4 I. Physikalisches Institut, Universitt zu Kln, Zlpicher Str. 77, 50937 Kln, Germany
E-mail contact: bonfand at mpifr-bonn.mpg.de
As the number of complex organic molecules (COMs) detected in the interstellar medium increases, it becomes
important to place meaningful constraints on the formation pathways of these species. The molecular cloud Sagittarius B2(N)-orth is host to several hot molecular cores in the early stage of star formation, where a great variety ofCOMs are detected in the gas phase. Because of its exposure to the extreme conditions of the Galactic center region,
Sgr B2(N) is one of the best targets to study the impact of environmental conditions on the production of COMs. Our
main goal is to characterize the physico-chemical evolution of Sgr B2(N)’s sources in order to explain their chemicaldifferences and constrain their environmental conditions. The chemical composition of Sgr B2(N)’s hot cores, N2,
N3, N4, and N5 is derived by modeling their 3 mm emission spectra extracted from the EMoCA imaging spectralline survey performed with ALMA. We derive the density distribution in the envelope of the sources based on themasses computed from the ALMA dust continuum emission maps. We use the radiative transfer code RADMC-3D
to compute temperature profiles based on the COM rotational temperatures derived from population diagrams. We
use published results of 3D RMHD simulations of high-mass star formation to estimate the time evolution of thesources properties. We employ the chemical code MAGICKAL to compute time-dependent chemical abundances inthe sources and investigate how physical properties and environmental conditions influence the production of COMs.
We find that chemical models with a cosmic-ray ionization rate of 7×10−16 s−1 best reproduce the abundances withrespect to methanol of ten COMs observed toward Sgr B2(N2-N5). We also show that COMs still form efficiently on
dust grains with minimum dust temperatures in the prestellar phase as high as 15 K, but that minimum temperatureshigher than 25 K are excluded.
Accepted by A&A

https://arxiv.org/pdf/1906.04695

6

GG Tau A: Dark shadows on the ringworld
Robert Brauer1,2 , Eric Pantin1 , Emmanuel Di Folco3,4 , Emilie Habart2 , Anne Dutrey3,4 and Stphane
Guilloteau3,4
1

CEA Saclay - Service d’Astrophysique, Orme des Merisiers, Bt 709, 91191 Gif sur Yvette, France; 2 Institutd’astrophysique spatiale, CNRS UMR 8617, Universit Paris-Sud 11, Bt 121, 91405, Orsay, France; 3 Univ. Bordeaux, Laboratoire d’Astrophysique de Bordeaux, UMR 5804, 33270 Floirac, France; 4 CNRS, LAB, UMR 5804,
33270 Floirac, France
E-mail contact: robert.brauer at cea.fr

Context. With its high complexity, large size, and close distance, the ringworld around GG Tau A is an appealing caseto study the formation and evolution of protoplanetary disks around multiple star systems. However, investigations
with radiative transfer models are usually neglecting the influence of the circumstellar dust around the individual
stars.
Aims. We investigate how circumstellar disks around the stars of GG Tau A are influencing the emission that is
scattered at the circumbinary disk and if constraints on these circumstellar disks can be derived.
12

Methods. We perform radiative transfer simulations with the code POLARIS to obtain spectral energy distributions
and emission maps in the H-Band (near-infrared). Subsequently, we compare them with observations to achieve ouraims.
Results. We studied the ratio of polarized intensity at different locations in the circumbinary disk and conclude thatthe observed scattered-light near-infrared emission is best reproduced, if the circumbinary disk lies in the shadow
of at least two co-planar circumstellar disks surrounding the central stars. This implies that the inner wall of thecircumbinary disk is strongly obscured around the midplane, while the observed emission is actually dominated by
the most upper disk layers. In addition, the inclined dark lane (”gap”) on the western side of the circumbinary disk,
which is a stable (non rotating) feature since ∼20 yr, can only be explained by the self-shadowing of a misaligned
circumstellar disk surrounding one of the two components of the secondary close-binary star GG Tau Ab.
Accepted by Astronomy & Astrophysics
https://arxiv.org/pdf/1906.11582

7

Sub-arcsecond (sub)millimeter imaging of the massive protocluster G358.93-0.03: Discovery of 14 new methanol maser lines associated with a hot core
C. L. Brogan1 , T. R. Hunter1 , A. P. M. Towner1,2 , B. A. McGuire1,3 , G. C. MacLeod4,5 , M. A. Gurwell3 ,C. J. Cyganowski6 , J. Brand7,8 , R. A. Burns9,10 , A. Caratti o Garatti11 , X. Chen12 , J. O. Chibueze13,14 ,N. Hirano15 , T. Hirota9,16 , K.-T. Kim10 , B. H. Kramer17,18 , H. Linz19 , K. M. Menten17 , A. Remijan1 ,A. Sanna17 , A. M. Sobolev20 , T. K. Sridharan3 , B. Stecklum21 , K. Sugiyama22 , G. Surcis23 , J. Van der
Walt24 , A. E. Volvach25,26 and L. N. Volvach25
1

National Radio Astronomy Observatory (NRAO), 520 Edgemont Rd, Charlottesville, VA 22903-2475, USA; 2 Department of Astronomy, University of Virginia, P.O. Box 3818, Charlottesville, VA 22904-4325, USA; 3 Center for
Astrophysics | Harvard & Smithsonian, Cambridge, MA 02138, USA; 4 Hartebeesthoek Radio Astronomy Observatory,
PO Box 443, Krugersdorp 1740, South Africa; 5 The University of Western Ontario, 1151 Richmond Street, London,
ON N6A 3K7, Canada; 6 SUPA, School of Physics and Astronomy, University of St. Andrews, North Haugh, St.
Andrews KY16 9SS, UK; 7 INAF-Istituto di Radioastronomia, via P. Gobetti 101, 40129 Bologna, Italy; 8 Italian
ALMA Regional Centre, via P. Gobetti 101, 40129 Bologna, Italy; 9 Mizusawa VLBI Observatory, National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan; 10 Korea Astronomy and Space Science
Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055, Republic of Korea; 11 Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, D02 XF86, Dublin, Ireland; 12 Center for Astrophysics, GuangZhou University, Guangzhou
510006, China; 13 Space Research Unit, Physics Department, North-West University, Potchefstroom, 2520, South
Africa; 14 Department of Physics and Astronomy, University of Nigeria, Carver Building, 1 University Road, Nsukka,
410001, Nigeria; 15 Institute of Astronomy and Astrophysics, Academia Sinica, 11F of Astronomy-Mathematics Building, AS/NTU, No.1, Sec. 4, Roosevelt Rd, Taipei 10617, Taiwan, R.O.C.; 16 Department of Astronomical Sciences,
SOKENDAI (The Graduate University for Advanced Studies), Osawa 2-21-1, Mitaka-shi, Tokyo 181-8588, Japan; 17
Max-Planck-Institut für Radioastronomie, Auf dem Hügel 69, D-53121 Bonn, Germany; 18 National Astronomical
Research Institute of Thailand, 260 Moo 4, T. Donkaew, Amphur Maerim, Chiang Mai, 50180, Thailand; 19 MaxPlanck-Institut für Astronomie, Königstuhl 17, 69117 Heidelberg, Germany; 20 Ural Federal University, 19 Mira street,
620002 Ekaterinburg, Russia; 21 Thüringer Landessternwarte, Sternwarte 5, 07778 Tautenburg, Germany; 22 National
Astronomical Observatory of Japan, Osawa 2-21-1, Mitaka-shi, Tokyo 181-8588, Japan; 23 INAF-Osservatorio Astronomico di Cagliari, Via della Scienza 5, I-09047, Selargius, Italy; 24 Centre for Space Research, North-West University,
Potchefstroom, South Africa; 25 Radio Astronomy Laboratory of Crimean Astrophysical Observatory, Katsively, RT-22
Crimea; 26 Institute of Applied Astronomy, Russian Academy of Sciences, St. Petersburg, Russia
E-mail contact: cbrogan at nrao.edu
We present (sub)millimeter imaging at 0.′′ 5 resolution of the massive star-forming region G358.93−0.03 acquired in
multiple epochs at 2 and 3 months following the recent flaring of its 6.7 GHz CH3 OH maser emission. Using SMA
and ALMA, we have discovered 14 new Class II CH3 OH maser lines ranging in frequency from 199 to 361 GHz, which
originate mostly from vt =1 torsionally-excited transitions and include one vt =2 transition. The latter detection
provides the first observational evidence that Class II maser pumping involves levels in the vt =2 state. The masers
are associated with the brightest continuum source (MM1), which hosts a line-rich hot core. The masers present a
consistent curvilinear spatial velocity pattern that wraps around MM1, suggestive of a coherent physical structure13

1200 au in extent. In contrast, the thermal lines exhibit a linear pattern that crosses MM1 but at progressive position
angles that appear to be a function of either increasing temperature or decreasing optical depth. The maser spectral
profiles evolved significantly over one month, and the intensities dropped by factors of 3.0 to 7.2, with the vt =2 line
showing the largest change. A small area of maser emission from only the highest excitation lines closest to MM1
has disappeared. There are seven additional dust continuum sources in the protocluster, including another hot core(MM3). We do not find evidence for a significant change in (sub)millimeter continuum emission from any of thesources during the one month interval, and the total protocluster emission remains comparable to prior single dishmeasurements.
Accepted by The Astrophysical Journal Letters
http://arxiv.org/pdf/1907.02470

8

Modeling C-Shock Chemistry in Isolated Molecular Outflows
Andrew M. Burkhardt1,2 , Christopher N. Shingledecker3,4,5 , Romane Le Gal1 , Brett A. McGuire1,6 ,Anthony J. Remijan6 and Eric Herbst2,5
1

Harvard-Smithsonian Center for Astrophysics, Cambridge, MA 02138, USA; 2 Astronomy Department, Universityof Virginia, Charlottesville, VA 22904, USA; 3 Max-Planck-Institut für extraterrestrische Physik, D-85748 Garching,
Germany; 4 Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569, Germany; 5 Chemistry Department, University of Virginia, Charlottesville, VA 22904, USA; 6 National Radio Astronomy Observatory,
Charlottesville, VA 22904, USA
E-mail contact: andrew.burkhardt at cfa.harvard.edu
Shocks are a crucial probe for understanding the ongoing chemistry within ices on interstellar dust grains where
many complex organic molecules (COMs) are believed to be formed. However, previous work has been limited to theinitial liberation into the gas phase through non-thermal desorption processes such as sputtering. Here, we present
results from the adapted three-phase gas-grain chemical network code NAUTILUS, with the inclusion of additional
high-temperature reactions, non-thermal desorption, collisional dust heating, and shock-physics parameters. Thisenhanced model is capable of reproducing many of the molecular distributions and abundance ratios seen in ourprior observations of the prototypical shocked-outflow L1157. In addition, we find that, among others, NH2 CHO,HCOOCH3 , and CH3 CHO have significant post-shock chemistry formation routes that differ from those of manyother COMs observed in shocks. Finally, a number of selected species and phenomena are studied here with respect
to their usefulness as shock tracers in various astrophysical sources.
Accepted by ApJ
http://arxiv.org/pdf/1906.08774

9

Molecular complexity on disk-scales uncovered by ALMA: The chemical composition of
the high-mass protostar AFGL 4176
Eva G. Bøgelund1 , Andrew G. Barr1 , Vianney Taquet1,2 , Niels F.W. Ligterink3 , Magnus V. Persson4 ,Michiel R. Hogerheijde1,5 , and Ewine F. van Dishoeck1,6
1

Leiden Observatory, Leiden University, PO Box 9513, 2300 RA Leiden, The Netherlands; 2 INAF, Osservatorio
Astrofisico di Arcetri, Largo E. Fermi 5, 50125 Firenze, Italy; 3 Center for Space and Habitability (CSH), University ofBern, Sidlerstrasse 5, 3012 Bern, Switzerland; 4 Department of Space, Earth and Environment, Chalmers Universityof Technology, Onsala Space Observatory, 439 92, Onsala, Sweden; 5 Anton Pannekoek Institute for Astronomy,
University of Amsterdam, Science Park 904, 1098 XH Amsterdam, The Netherlands; 6 Max-Planck Institut für
Extraterrestrische Physik, Giessenbachstr. 1, 85748 Garching, Germany
E-mail contact: bogelund at strw.leidenuniv.nl
The chemical composition of high-mass protostars reflects the physical evolution associated with different stages of starformation. In this study, the molecular inventory of the forming high-mass star AFGL 4176 is studied in detail at highangular resolution (∼0.′′ 35) using ALMA. This high resolution makes it possible to separate the emission associatedwith the inner hot envelope and disk around the forming star from that of its cool outer envelope. In addition, thehigh sensitivity of ALMA makes it possible to identify weak and optically thin lines and allows for many isotopologues
14to be detected, providing a more complete and accurate inventory of the source. A total of 23 different molecularspecies and their isotopologues are detected in the spectrum towards AFGL 4176. The most abundant species ismethanol (CH3 OH), remaining species are present at levels between 0.003% and 15% with respect to CH3 OH. Hints
that N-bearing species peak slightly closer to the location of the peak continuum emission than the O-bearing species
are seen. AFGL 4176 comprises a rich chemical inventory including many complex species present on disk-scales. On
average, the derived column density ratios with respect to methanol of O-bearing species are higher than those derived
for N-bearing species by a factor of three. This may indicate that AFGL 4176 is a relatively young source since nitrogenchemistry generally takes longer to evolve in the gas-phase. Taking methanol as a reference, the composition of AFGL
4176 more closely resembles that of the low-mass protostar IRAS 16293−2422B than that of high-mass star-forming
regions located near the Galactic centre. This similarity hints that the chemical composition of complex species isalready set in the cold cloud stage and implies that AFGL 4176 is a young source whose chemical composition has
not yet been strongly processed by the central protostar.
Accepted by A&A
http://arxiv.org/pdf/1906.06156

Upper limits on the water vapour content of the β Pictoris debris disk

10 M. Cavallius1, G. Cataldi1,2,3,4, A. Brandeker1 , G. Olofsson1, B. Larsson1 , and R. Liseau5
1

Department of Astronomy, AlbaNova University Center, Stockholm University, 106 91 Stockholm, Sweden; 2 Konkoly
Observatory, Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, Konkoly-Thege
Miklós út 15–17, 1121 Budapest, Hungary; 3 Subaru Telescope, National Astronomical Observatory of Japan, 650
North A’ohoku Place, Hilo, HI 96720, USA; 4 Department of Astronomy and Astrophysics, University of Toronto,
Toronto, ON M5S 3H4, Canada; 5 Department of Space, Earth and Environment, Chalmers University of Technology,
Onsala Space Observatory, 439 92 Onsala, Sweden
E-mail contact: maria.cavallius at astro.su.se
The debris disk surrounding β Pictoris has been observed with ALMA to contain a belt of CO gas with a distinct
peak at ∼85 au. This CO clump is thought to be the result of a region of enhanced density of solids that collide
and release CO through vaporisation. The parent bodies are thought to be comparable to solar system comets, inwhich CO is trapped inside a water ice matrix. Since H2 O should be released along with CO, we aim to put an upper
limit on the H2 O gas mass in the disk of β Pictoris. We use archival data from the Heterodyne Instrument for the
Far-Infrared (HIFI) aboard the Herschel Space Observatory to study the ortho-H2 O 110 –101 emission line. The line
is undetected. Using a python implementation of the radiative transfer code RADEX, we convert upper limits on theline flux to H2 O gas masses. The resulting lower limits on the CO/H2 O mass ratio are compared to the composition ofsolar system comets. Depending on the assumed gas spatial distribution, we find a 95% upper limit on the ortho-H2 Oline flux of 7.5 × 10−20 W m−2 or 1.2 × 10−19 W m−2 . These translate into an upper limit on the H2 O mass of
7.4 × 1016 –1.1 × 1018 kg depending on both the electron density and gas kinetic temperature. The range of derived
gas-phase CO/H2 O ratios is marginally consistent with low-ratio solar system comets.
Accepted by A&A
http://arxiv.org/pdf/1906.11106

Silicate features in the circumstellar envelopes of the Class I binary driving source of

11 HH 250

F. Comerón1 , B. Merı́n2 , B. Reipurth3 and H.-W. Yen1
1

European Southern Observatory, Karl-Schwarzschild-Strasse 2, D-85748 Garching bei München, Germany; 2 European Space Astronomy Centre (ESAC), European Space Agency (ESA), E-28691 Villanueva de la Cañada, Spain; 3
Institute for Astronomy, University of Hawaii at Manoa, 640 N. Aohoku Place, Hilo, HI 96720, USA
E-mail contact: fcomeron at eso.org
The silicate feature near 10µm one of the main tools available to study the mineralogy of circumstellar disks andenvelopes, providing information on the thermal processing, growth, location, and circulation of dust grains. We
investigate the silicate feature of the two Class I components of HH 250-IRS, a resolved binary system with a separation15

of 0”53 driving a Herbig-Haro flow. Each component has its own circumstellar envelope, and the system is surrounded
by a circumbinary disk. We have carried out low resolution spectroscopy in the 8-13µm range using VISIR, thethermal infrared imager and spectrograph at ESOs Very Large Telescope. The silicate features of both sources areclearly different. The NW component has a broad, smooth absorption profile lacking structure. We attribute most of it
to foreground interstellar dust absorption, but estimate that additional absorption by amorphous silicates takes placein the circumstellar envelope of the young stellar object. The SE component shows the silicate feature in emission,
with structure longwards of 9.5µm indicating the presence of crystalline dust in the dominant form of forsterite. The
apparent lack of an absorption feature caused by foreground dust is probably due to the filling of the band withemission by amorphous silicates in the envelope of the object. Despite their virtually certain coevality, the differencesin the components of the HH 250-IRS binary are most likely due to markedly different circumstellar environments.
The NW component displays an unevolved envelope, whereas dust growth and crystallization has taken place in the
SE component. The weak or absent signatures of enstatite in the latter are fairly unusual among envelopes withcrystalline dust, and we tentatively relate it to a possible wide gap or an inner truncation of the disk already hinted
in previous observations by a drop in the L’-band flux, which might indicate that the SE component could actually be
a very close binary. We speculate that the clear differences between the silicate feature spectra of both components of
HH 250-IRS may be due either to disk evolution sped up by multiplicity, or by accretion variability leading to episodesof crystal formation. Different inclinations with respect to the line of sight may play a role as well, although it is very
unlikely that they are the sole responsible for the differences between both objects.
Accepted by Astronomy and Astrophysics
http://www.eso.org/~fcomeron/hh250_spec_final.pdf

Planet formation and stability in polar circumbinary discs

12 Nicolás Cuello1,2 and Cristian A. Giuppone3
1

Instituto de Astrofı́sica, Pontificia Universidad Católica de Chile, Santiago, Chile; 2 Núcleo Milenio de FormaciónPlanetaria (NPF), Chile; 3 Universidad Nacional de Córdoba, Observatorio Astronómico - IATE, Laprida 854,
X5000BGR Córdoba, Argentina
E-mail contact: cuellonicolas at gmail.com
Context. Dynamical studies suggest that most of the circumbinary discs (CBDs) should be coplanar, i.e. the rotation
vectors of the binary and the disc are aligned. However, some theoretical works show that under certain conditions a
CBD can become polar, which means that its rotation vector is orthogonal with respect to the binary orbital plane.
Interestingly, very recent observations showed that polar CBDs exist in Nature (e.g. HD 98800).
Aims. We test the predictions of CBD alignment around eccentric binaries based on linear theory. In particular, we
compare prograde and retrograde CBD configurations. Then, assuming planets form in these systems, we thoroughly
characterise the orbital behaviour and stability of misaligned (P-type) particles. This is done for massless and massive
particles.
Methods. The evolution of the CBD alignment for various configurations is modelled through three-dimensionalhydrodynamical simulations. For the orbital characterisation and the analysis stability, we relied on long-term N-body
integrations and structure and chaos indicators, such as ∆e and MEGNO.
Results. We confirm previous analytical predictions on CBD alignment, but find an unexpected symmetry breaking
between prograde and retrograde configurations. More specifically, we observe polar alignment for a retrograde misaligned CBD that was expected to become coplanar with respect to the binary disc plane. Therefore, the likelihood of
becoming polar for a highly misaligned CBD is higher than previously thought. Regarding the stability of circumbinary P-type planets (aka Tatooines), polar orbits are stable over a wide range of binary parameters. In particular,
for binary eccentricities below 0.4, the orbits are stable for any value of the binary mass ratio. In the absence ofgas, planets with masses below 10−5 M⊙ have negligible effects on the binary orbit. Finally, we predict that Polar
Tatooines should be searched around mildly eccentric equal-mass binaries.
Accepted by Astronomy & Astrophysics
https://arxiv.org/pdf/1906.10579

16

13

The Degree of Alignment Between Circumbinary Disks and Their Binary Hosts
Ian Czekala1 , Eugene Chiang1,2 , Sean M. Andrews3 , Eric L.N. Jensen4 , Guillermo Torres3 , David J.
Wilner3 , Keivan G. Stassun5,6 and Bruce Macintosh7
1

Department of Astronomy, University of California at Berkeley, Campbell Hall, CA 94720-3411, USA; 2 Departmentof Earth and Planetary Science, University of California at Berkeley, McCone Hall, Berkeley, CA 94720-4767, USA;
3
Center for Astrophysics — Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA; 4 Departmentof Physics and Astronomy, Swarthmore College, 500 College Avenue, Swarthmore, PA 19081, USA; 5 Departmentof Physics and Astronomy, Vanderbilt University, 6301 Stevenson Center, Nashville, TN 37235, USA; 6 Departmentof Physics, Fisk University, Nashville, TN 37208, USA; 7 Kavli Institute for Particle Astrophysics and Cosmology,
Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA
E-mail contact: iczekala at berkeley.edu
All four circumbinary (CB) protoplanetary disks orbiting short-period (P < 20 day) double-lined spectroscopic binaries(SB2s)—a group that includes UZ Tau E, for which we present new ALMA data—exhibit sky-plane inclinations idisk
which match, to within a few degrees, the sky-plane inclinations i⋆ of their stellar hosts. Although for these systems
the true mutual inclinations θ between disk and binary cannot be directly measured because relative nodal angles are
unknown, the near-coincidence of idisk and i⋆ suggests that θ is small for these most compact of systems. We confirm
this hypothesis using a hierarchical Bayesian analysis, showing that 68% of CB disks around short-period SB2s have
θ < 3.0◦ . Near co-planarity of CB disks implies near co-planarity of CB planets discovered by Kepler, which in turnimplies that the occurrence rate of close-in CB planets is similar to that around single stars. By contrast, at longer
periods ranging from 30–105 days (where the nodal degeneracy can be broken via, e.g., binary astrometry), CB disks
exhibit a wide range of mutual inclinations, from co-planar to polar. Many of these long-period binaries are eccentric,
as their component stars are too far separated to be tidally circularized. We discuss how theories of binary formationand disk-binary gravitational interactions can accommodate all these observations.
Accepted by ApJ
http://arxiv.org/pdf/1906.03269

14

Dynamics of magnetic flux tubes in accretion discs of T Tauri stars
A.E. Dudorov1 , S.A. Khaibrakhmanov1,2 , A.M. Sobolev2
1

Chelyabinsk state university, 129 Br. Kashirinykh str., Chelyabinsk 454001, Russia; 2 Ural Federal University, 51
Lenin str., Ekaterinburg 620000, Russia
E-mail contact: dudorov at csu.ru
Dynamics of slender magnetic flux tubes (MFT) in the accretion discs of T Tauri stars is investigated. We perform
simulations taking into account buoyant, aerodynamic and turbulent drag forces, radiative heat exchange between
MFT and ambient gas, magnetic field of the disc. The equations of MFT dynamics are solved using Runge-Kutta
method of the fourth order. The simulations show that there are two regimes of MFT motion in absence of external
magnetic field. In the region r < 0.2 au, the MFT of radii 0.05 ≤ a0 ≤ 0.16H (H is the scale height of the disc) withinitial plasma β of 1 experience thermal oscillations above the disc. The oscillations decay over some time, and MFTcontinue upward motion afterwards. Thinner or thicker MFT do not oscillate. MFT velocity increases with initial
radius and magnetic field strength. MFT rise periodically with velocities up to 5–15 km s−1 and periods of 0.5–10
yr determined by the toroidal magnetic field generation time. Approximately 20% of disc mass and magnetic flux
can escape to disc atmosphere via the magnetic buoyancy over characteristic time of disc evolution. MFT dispersal
forms expanding magnetized corona of the disc. External magnetic field causes MFT oscillations near the disc surface.
These magnetic oscillations have periods from several days to 1–3 months at r < 0.6 au. The magnetic oscillations
decay over few periods. We simulate MFT dynamics in accretion discs in the Chameleon I cluster. The simulations
demonstrate that MFT oscillations can produce observed IR-variability of T Tauri stars.
Accepted by MNRAS
http://arxiv.org/pdf/1906.05133

17

15

Double-peaked OI profile: a likely signature of the gaseous ring around KH 15D
Min Fang1,2 , Ilaria Pascucci3,2 , Jinyoung Serena Kim1,2 and Suzan Edwards4
1

Department of Astronomy, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA; 2 Earths in
Other Solar Systems Team, NASA Nexus for Exoplanet System Science; 3 Department of Planetary Sciences, Universityof Arizona, 1629 East University Boulevard, Tucson, AZ 85721, USA; 4 Five College Astronomy Department, Smith
College, Northampton, MA 01063, USA
E-mail contact: mfang.cn at gmail.com

KH 15D is a well-known spectroscopic binary because of its unique and dramatic photometric variability. The variability is explained by a circumbinary dust ring but the ring itself was never directly detected. We present a newinterpretation of the double-peaked [O i] λ6300 profiles as originating from the hot disk surface of KH 15D. By modeling these profiles, we measure emitting radii between ∼0.5–5 au, basically a gaseous ring very similar in radial extent
to the dust ring inferred from modeling the system’s photometric variability. We also discuss the possibility thatexternal photoevaporation driven by UV photons from the nearby massive star HD 47887 has truncated the outeredge of the disk to the observed value.
Accepted by ApJL
http://arxiv.org/pdf/1906.04723

16

CO Detected in CI Tau b: Hot Start Implied by Planet Mass and M K
Laura Flagg1 , Christopher M. Johns-Krull1,2 , Larissa Nofi3,4 , Joe Llama3 , L. Prato2,3 , Kendall Sullivan5 ,D. T. Jaffe5 and Gregory Mace5
1Department of Physics and Astronomy, Rice University, 6100 Main St. MS-108, Houston, TX 77005, USA; 2 Visiting
Astronomer at McDonald Observatory of The University of Texas at Austin; 3 Lowell Observatory, 1400 W Mars Hill
Road, Flagsta?, AZ 86001, USA; 4 Institute for Astronomy, University of Hawaii Manoa, 2680 Woodlawn Drive,
Honolulu, HI 96822, USA; 5 McDonald Observatory and the Department of Astronomy, The University of Texas atAustin, Austin, TX 78712, USA
E-mail contact: laura.flagg at rice.edu
We acquired high resolution IR spectra of CI Tau, the host star of one of the few young planet candidates amenable
to direct spectroscopic detection. We confirm the planet’s existence with a direct detection of CO in the planet’s
atmosphere. We also calculate a mass of 11.6 MJ based on the amplitude of its radial velocity variations. We
estimate its flux contrast with its host star to get an absolute magnitude estimate for the planet of 8.17 in the Kband. This magnitude implies the planet formed via a “hot start” formation mechanism. This makes CI Tau b
the youngest confirmed exoplanet as well as the first exoplanet around a T Tauri star with a directly determined,
model-independent, dynamical mass.
Accepted by ApJ letters
https://arxiv.org/pdf/1906.02860

17

The SPHERE view of the jet and the envelope of RY Tau
A. Garufi1 , L. Podio1 , F. Bacciotti1 , S. Antoniucci2 , A. Boccaletti3 , C. Codella1,4 , C. Dougados4 , F.
Menard4 , D. Mesa5 , M. Meyer6 , B. Nisini2, H.M. Schmid7 , T. Stolker7 , J.L. Baudino8 , B. Biller9 , M.
Bonavita9 , M. Bonnefoy4 , F. Cantalloube10 , G. Chauvin4 , A. Cheetham11 , S. Desidera5 , V. D’Orazi5 ,M. Feldt10 , R. Galicher3 , A. Grandjean4 , R. Gratton5 , J. Hagelberg7, A.M. Lagrange4 , M. Langlois12 ,J. Lannier4 , C. Lazzoni13 , A.L. Maire14 , C. Perrot3,15,16 , E. Rickman11 , T. Schmidt3 , A. Vigan17 and A.
Zurlo17,18,19
1

INAF, Osservatorio Astrofisico di Arcetri, Largo Enrico Fermi 5, I-50125 Firenze, Italy; 2 INAF - Osservatorio
Astronomico di Roma, Via Frascati 33, 00078, Monte Porzio Catone (RM), Italy; 3 LESIA, Observatoire de Paris,
Université PSL, CNRS, Sorbonne Université, Univ. Paris Diderot, Sorbonne Paris Cité, 5 place Jules Janssen, 92195
Meudon, France; 4 Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France; 5 INAF - Osservatorio Astronomico
di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy; 6 Department of Astronomy, University of Michigan, 1085
18

S. University Ave, Ann Arbor, MI 48109-1107, USA; 7 Institute for Particle Physics and Astrophysics, ETH Zurich,
Wolfgang-Pauli-Strasse 27, 8093 Zurich, Switzerland; 8 Department of Physics, University of Oxford, Oxford, UK; 9
Institute for Astronomy, University of Edinburgh, EH9 3HJ, Edinburgh, UK; 10 Max Planck Institute for Astronomy,Königstuhl 17, 69117 Heidelberg, Germany; 11 Geneva Observatory, University of Geneva, Ch. des Maillettes 51,
1290, Versoix, Switzerland; 12 CRAL, CNRS, Université Lyon 1, 9 avenue Charles André, 69561 Saint Genis Laval
Cedex, France; 13 Università degli Studi di Padova, dipartimento di Fisica e Astronomia, vicolo dell’osservatorio 3,35122 Padova; 14 STAR Institute, University of Liège, Allée du Six Août 19c, B-4000 Liège, Belgium; 15 Instituto deFı́sica y Astronomı́a, Facultad de Ciencias, Universidad de Valparaı́so, Av. Gran Bretaña 1111, Valparaı́so, Chile; 16
Núcleo Milenio Formación Planetaria - NPF, Universidad de Valparaı́so, Av. Gran Bretaña 1111, Valparaı́so, Chile;
17
Aix Marseille Université, CNRS, LAM - Laboratoire d’Astrophysique de Marseille, UMR 7326, 13388, Marseille,
France; 18 Núcleo de Astronomı́a, Facultad de Ingenierı́a y Ciencias, Universidad Diego Portales, Av. Ejercito 441,
Santiago, Chile; 19 Escuela de Ingenierı́a Industrial, Facultad de Ingenierı́a y Ciencias, Universidad Diego Portales,
Av. Ejercito 441, Santiago, Chile
E-mail contact: agarufi at arcetri.astro.it
Jets are rarely associated with pre-main-sequence intermediate-mass stars. Optical and near-IR observations of jetdriving sources are often hindered by the presence of a natal envelope. Jets around partly embedded sources are auseful diagnostic to constrain the geometry of the concealed protoplanetary disk. In fact, the jet-driving mechanisms
are affected by both spatial anisotropies and episodic variations at the (sub-)au scale from the star. We obtained a rich
set of high-contrast VLT/SPHERE observations from 0.6 micron to 2.2 micron of the young intermediate-mass star
RY Tau. Given the proximity to the Sun of this source, our images have the highest spatial resolution ever obtained
for an atomic jet. Optical observations in polarized light show no sign of the protoplanetary disk detected by ALMA.
Instead, we observed a diffuse signal resembling a remnant envelope with an outflow cavity. The jet is detected infour spectral lines. The jet appears to be wiggling and its radial width increasing with the distance is complementary
to the shape of the outflow cavity suggesting a strong jet/envelope interaction. Through the estimated tangential
velocity, we revealed a possible connection between the launching time of the jet sub-structures and the stellar activityof RY Tau. RY Tau is at an intermediate stage toward the dispersal of the natal envelope. This source shows episodic
increases of mass accretion/ejection similarly to other known intermediate-mass stars. The amount of observed jet
wiggle is consistent with the presence of a precessing disk warp or misaligned inner disk that would be induced byan unseen planetary/sub-stellar companion at sub-/few-au scales. The high disk mass of RY Tau and of two otherjet-driving intermediate-mass stars, HD163296 and MWC480, suggests that massive, full disks are more efficient atlaunching prominent jets.
Accepted by A&A
https://arxiv.org/pdf/1906.06910

18

Warm gas in protostellar outflows. II. EHV jet and outflows from OMC-2/3
A.I. Gómez-Ruiz1 , A. Gusdorf2 , S. Leurini3 , K.M. Menten4 , S. Takahashi5,6,7 , F. Wyrowski4 , R. Güsten4
1

CONACYT–Instituto Nacionalde Astrofı́sica, Ópticay Electrónica, Luis E. Erro1, C.P. 72840, Tonantzintla, Puebla,México; 2 LERMA, UMR 8112 du CNRS, Observatoire de Paris, École Normale Supérieure, 61 Av. de l’Observatoire,
75014, Paris, France; 3 INAF-Osservatorio Astronomico di Cagliari, Via della Scienza 5, I-09047, Selargius (CA),
Italy; 4 Max-Planck-Institut für Radioastronomie (MPIfR), Auf dem Hügel 69, 53121 Bonn, Germany; 5 Joint ALMA
Observatory, Alonso de Córdova 3107, Vitacura, Santiago, Chile; 6 NAOJ Chile Observatory, Alonso de Córdova
3788, Oficina 61B, Vitacura, Santiago, Chile; 7 Department of Astronomical Science, School of Physical Sciences,
SOKENDAI (The Graduate University for Advanced Studies), Mitaka, Tokyo 181-8588, Japan
E-mail contact: aigomez at inaoep.mx
OMC-2/3 is one of the nearest embedded cluster-forming region that includes intermediate-mass protostars at earlystages of evolution. A previous CO (3–2) mapping survey towards this region revealed outflow activity related tosources at different evolutionary phases. The present work aims to study the warm gas in the high-velocity emission
from several outflows found in CO (3–2) emission by previous observations, determine their physical conditions, andmake comparison with previous results in low-mass star-forming regions. We use the CHAMP+ heterodyne array on
the APEX telescope to map the CO (6–5) and CO (7–6) emission in the OMC-2 FIR 6 and OMC-3 MMS 1-6 regions,and to observe 13 CO (6–5) at selected positions. We analyze these data together with previous CO (3–2) observations.
19

In addition, we mapped the SiO (5–4) emission in OMC-2 FIR 6. The CO (6–5) emission was detected in most of theoutflow lobes in the mapped regions, while the CO (7–6) was found mostly from the OMC-3 outflows. In the OMC-3
MMS 5 outflow, a previously undetected extremely high velocity gas was found in CO (6–5). This extremely highvelocity emission arises from the regions close to the central object MMS 5. Radiative transfer models revealed thatthe high-velocity gas from MMS 5 outflow consists of gas with nH2 = 104 –105 cm−3 and T > 200 K, similar to what
is observed in young Class 0 low-mass protostars. For the other outflows, values of nH2 > 104 cm−3 were found. The
physical conditions and kinematic properties of the young intermediate-mass outflows presented here are similar to
those found in outflows from Class 0 low-mass objects. Due to their excitation requirements, mid-J CO lines are good
tracers of extremely high velocity gas in young outflows likely related to jets.
Accepted by A&A
http://arxiv.org/pdf/1906.03522

Two accreting protoplanets around the young star PDS 70

19

S.Y. Haffert1 , A.J. Bohn1 , J. de Boer1 , I.A.G. Snellen1 , J. Brinchmann1,2 , J.H. Girard3,4 , C.U. Keller1 ,R. Bacon5
1

Leiden Observatory, Leiden University, Postbus 9513, 2300 RA Leiden, The Netherlands; 2 Instituto de Astrofı́sica eCiências do Espaço, Universidade do Porto, CAUP, Rua das Estrelas, PT4150- 762 Porto, Portugal; 3 Space Telescope
Science Institute, 3700 San Martin Drive, Baltimore, MD 21218, USA; 4 Université Grenoble Alpes, CNRS, IPAG,
38000 Grenoble, France; 5 CRAL, Observatoire de Lyon, CNRS, Université Lyon 1, 9 Avenue Ch. André, F-69561
Saint Genis Laval Cedex, France
E-mail contact: haffert at strw.leidenuniv.nl
Newly forming proto-planets are expected to create cavities and substructures in young, gas-rich proto-planetarydisks, but they are difficult to detect as they could be confused with disk features affected by advanced image-analysis
techniques. Recently, a planet was discovered inside the gap of the transitional disk of the T-Tauri star PDS 70. Here
we report on the detection of strong Hα emission from two distinct locations in the PDS 70 system, one corresponding
to the previously discovered planet PDS 70 b, which confirms the earlier Hα detection, and another located close to
the outer-edge of the gap, coinciding with a previously identified bright dust spot in the disk and with a small opening
in a ring of molecular emission. We identify this second Hα peak as a second proto-planet in the PDS 70 system. The
Hα emission spectra of both proto-planets indicate ongoing accretion onto the proto-planets, which appear to be neara 2:1 mean motion resonance. Our observations show that adaptive-optics-assisted, medium-resolution, integral-field
spectroscopy with MUSE targeting accretion signatures will be a powerful way to trace ongoing planet formationin transitional disks at different stages of their evolution. Finding more young planetary systems in mean motion
resonance would give credibility to the Grand Tack hypothesis in which Jupiter and Saturn migrated in a resonanceorbit during the early formation period of our Solar System.
Accepted by Nature Astronomy
http://arxiv.org/pdf/1906.01486

Submillimeter emission associated with candidate protoplanets

20

Andrea Isella1 , Myriam Benisty2,3,4 , Richard Teague5 , Jaehan Bae6 , Miriam Keppler7 , Stefano Facchini8and Laura Pérez2
1

Department of Physics and Astronomy, Rice University, 6100 Main Street, MS-108, Houston, TX 77005, USA; 2
Departamento de Astronomı́a, Universidad de Chile, Camino El Observatorio 1515, Las Condes, Santiago, Chile;
3
Unidad Mixta Internacional Franco-Chilena de Astronomı́a, CNRS, UMI 3386; 4 Univ. Grenoble Alpes, CNRS,
IPAG, 38000 Grenoble, France; 5 Department of Astronomy, University of Michigan, 311 West Hall, 1085 S. University
Avenue, Ann Arbor, MI 48109, USA; 6 Department of Terrestrial Magnetism, Carnegie Institution for Science, 5241
Broad Branch Road, NW, Washington, DC 20015, USA; 7 Max Planck Institute for Astronomy, Königstuhl 17, 69117,
Heidelberg, Germany; 8 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching, Germany
E-mail contact: isella at rice.edu
We present the discovery of a spatially unresolved source of sub-millimeter continuum emission (λ = 855µm) associated20

with a young planet, PDS 70 c, recently detected in Hα emission around the 5 Myr old T Tauri star PDS 70. We
interpret the emission as originating from a dusty circumplanetary disk with a dust mass between 2 × 10−3 and4.2 × 10−3 Earth masses. Assuming a standard gas-to-dust ratio of 100, the ratio between the total mass of thecircumplanetary disk and the mass of the central planet would be between 10−4 –10−5 . Furthermore, we report thediscovery of another compact continuum source located 0.′′ 074±0.′′ 013 South-West of a second known planet in thissystem, PDS 70 b, that was previously detected in near-infrared images. We speculate that the latter source might
trace dust orbiting in proximity of the planet, but more sensitive observations are required to unveil its nature.
Accepted by ApJL
http://arxiv.org/pdf/1906.06156

21Intrinsic polarisation of elongated porous dust grains
Florian Kirchschlager1 , Gesa H.-M. Bertrang2 and Mario Flock2
1

University College London, Gower Street, WC1E 6BT London, United Kingdom; 2 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany
E-mail contact: f.kirchschlager at ucl.ac.uk
ALMA observations revealed recently polarised radiation of several protoplanetary disks in the (sub-)millimetre wavelength range. Besides self-scattering of large particles, thermal emission by elongated grains is a potential source
for the detected polarisation signal. We calculate the wavelength dependent absorption and intrinsic polarisation of
spheroidally shaped, micrometre and sub-millimetre sized dust grains using the discrete dipole approximation. In particular, we analyse the impact of dust grain porosity which appears to be present in disks when small grains coagulate
to form larger aggregates. For the first time our results show that (a) the intrinsic polarisation decreases for increasing
grain porosity and (b) the polarisation orientation flips by 90 degree for certain ratios of wavelength to grain size. We
present a new method to constrain grain porosity and the grain size in protoplanetary disks using multi-wavelengthpolarisation observations in the far-infrared to millimetre wavelengths. Finally, we find that moderate grain porosities(P ≤ 0.7) potentially explain the observed polarisation fraction in the system HD 142527 while highly porous grains(P > 0.7) fail unless the grain’s axis ratio is extraordinarily large.
Accepted by MNRAS
https://arxiv.org/pdf/1906.10699

22

The Bridge: a transient phenomenon of forming stellar multiples
Michael Kuffmeier1 , Hannah Calcutt2,3 and Lars E. Kristensen2
1
Zentrum für Astronomie der Universität Heidelberg, Institut für Theoretische Astrophysik, Albert-Ueberle-Str. 2,69120 Heidelberg; 2 Centre for Star and Planet Formation (StarPlan) & Niels Bohr Institute, University of Copenhagen,Øster Voldgade 5–7, DK-1350 Copenhagen, Denmark; 3 Department of Space, Earth and Environment, Chalmers
University of Technology, Hörsalsvägen 11, 412 96, Gothenburg

E-mail contact: kuffmeier at uni-heidelberg.de
Observations with modern instruments such as Herschel reveal that stars form clustered inside filamentary arms of ∼1
pc length embedded in Giant Molecular Clouds. On smaller scales (∼1000 au), observations of, e.g., IRAS 16293–2422
show signs of filamentary ‘bridge’ structures connecting young protostars in their birth environment. We aim to findthe origin of bridges associated with deeply embedded protostars, by characterizing their connection to the filamentary
structure present on GMC scales and to the formation of protostellar multiples. Using the magnetohydrodynamicalcode ramses, we carry out zoom-in simulations of low-mass star formation starting from Giant-Molecular-Cloud-scales.
We analyze the morphology and dynamics involved in the formation process of a triple system. Colliding flows ofgas in the filamentary arms induce the formation of two protostellar companions at distances of ∼1000 au fromthe primary. After their birth, the stellar companions quickly (∆t ∼ 10 kyr) approach and orbit the primary on
eccentric orbits with separations of ∼100 au. The colliding flows induce transient structures lasting for up to a few
10 kyr connecting two forming protostellar objects that are kinematically quiescent along the line-of-sight. Colliding
flows compress the gas and trigger the formation of stellar companions via turbulent fragmentation. Our results
suggest that protostellar companions initially form with a wide separation of ∼1000 au. Smaller separations (a ∼ 10021

au or less) are a consequence of subsequent migration and capturing. Associated with the formation phase of thecompanion, the turbulent environment induces the formation of arc- and bridge-like structures. These bridges canbecome kinematically quiescent, when the velocity components of the colliding flows eliminate each other. However,
the gas in bridges still contributes to stellar accretion later. Our results demonstrate: bridge-like structures are atransient phenomenon of stellar multiple formation.
Accepted by A&A
http://arxiv.org/pdf/1907.02083

Compact Disks in a High Resolution ALMA Survey of Dust Structures in the Taurus

23 Molecular Cloud

Feng Long1,2 , Gregory J. Herczeg1 , Daniel Harsono3 , Paola Pinilla4,5 , Marco Tazzari6 , Carlo F.
Manara7 , Ilaria Pascucci8,9 , Sylvie Cabrit10,11 , Brunella Nisini12, Doug Johnstone13,14 Suzan Edwards15 ,Colette Salyk16 , Francois Menard11 , Giuseppe Lodato17 , Yann Boehler11,18 , Gregory N. Mace19 , Yao
Liu20,21 , Gijs D. Mulders22,9 , Nathanial Hendler8,23 , Enrico Ragusa24 , William J. Fischer25 , Andrea
Banzatti8 , Elisabetta Rigliaco26 , Gerrit van de Plas11 , Giovanni Dipierro24 , Michael Gully-Santiago27
and Ricardo Lopez-Valdivia19
1

Kavli Institute for Astronomy and Astrophysics, Peking University, Beijing 100871, China; 2 Department of Astronomy, School of Physics, Peking University, Beijing 100871, China; 3 Leiden Observatory, Leiden University, P.O.
box 9513, 2300 RA Leiden, The Netherlands; 4 Department of Astronomy/Steward Observatory, The University of
Arizona, 933 North Cherry Avenue, Tucson, AZ 85721, USA; 5 Max-Planck-Institut für Astronomie, Königstuhl 17,69117, Heidelberg, Germany; 6 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB30HA, UK; 7 European Southern Observatory, Karl-Schwarzschild-Str. 2, D-85748 Garching bei München, Germany;
8
Lunar and Planetary Laboratory, University of Arizona, Tucson, AZ 85721, USA; 9 Earths in Other Solar Systems
Team, NASA Nexus for Exoplanet System Science, USA; 10 Sorbonne Université, Observatoire de Paris, Université
PSL, CNRS, LERMA, F-75014 Paris, France; 11 Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France; 12
INAF–Osservatorio Astronomico di Roma, via di Frascati 33, 00040 Monte Porzio Catone, Italy; 13 NRC Herzberg
Astronomy and Astrophysics, 5071 West Saanich Road, Victoria, BC, V9E 2E7, Canada; 14 Department of Physics
and Astronomy, University of Victoria, Victoria, BC, V8P 5C2, Canada; 15 Five College Astronomy Department,
Smith College, Northampton, MA 01063, USA; 16 Vassar College Physics and Astronomy Department, 124 Raymond
Avenue, Poughkeepsie, NY 12604, USA; 17 Dipartimento di Fisica, Universita Degli Studi di Milano, Via Celoria,
16, I-20133 Milano, Italy; 18 Rice University, Department of Physics and Astronomy, Main Street, 77005 Houston,
USA; 19 McDonald Observatory and Department of Astronomy, University of Texas at Austin, 2515 Speedway, Stop
C1400, Austin, TX 78712-1205, USA; 20 Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstrasse 1,
85748, Garching, Germany; 21 Purple Mountain Observatory, Chinese Academy of Sciences, 2 West Beijing Road,
Nanjing 210008, China; 22 Department of the Geophysical Sciences, The University of Chicago, 5734 South Ellis
Avenue, Chicago, IL 60637, USA; 23 LSSTC Data Science Fellow; 24 Department of Physics and Astronomy, University of Leicester, Leicester LE1 7RH, UK; 25 Space Telescope Science Institute Baltimore, MD 21218, USA; 26
INAF-Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, 35122 Padova, Italy; 27 NASA Ames ResearchCenter and Bay Area Environmental Research Institute, Moffett Field, CA 94035, USA
E-mail contact: longfeng at pku.edu.cn
We present a high-resolution (∼0.′′ 12, ∼16 au, mean sensitivity of 50 µJy beam−1 at 225 GHz) snapshot survey of 32
protoplanetary disks around young stars with spectral type earlier than M3 in the Taurus star-forming region using
Atacama Large Millimeter Array (ALMA). This sample includes most mid-infrared excess members that were notpreviously imaged at high spatial resolution, excluding close binaries and highly extincted objects, thereby providing amore representative look at disk properties at 1–2 Myr. Our 1.3 mm continuum maps reveal 12 disks with prominent
dust gaps and rings, 2 of which are around primary stars in wide binaries, and 20 disks with no resolved features atthe observed resolution (hereafter smooth disks), 8 of which are around the primary star in wide binaries. The smoothdisks were classified based on their lack of resolved substructures, but their most prominent property is that they are
all compact with small effective emission radii (Reff,95% <∼ 50 au). In contrast, all disks with Reff,95% of at least 55 au
in our sample show detectable substructures. Nevertheless, their inner emission cores (inside the resolved gaps) have
similar peak brightness, power law profiles, and transition radii to the compact smooth disks, so the primary difference22between these two categories is the lack of outer substructures in the latter. These compact disks may lose their outerdisk through fast radial drift without dust trapping, or they might be born with small sizes. The compact dust disks,
as well as the inner disk cores of extended ring disks, that look smooth at the current resolution will likely show
small-scale or low-contrast substructures at higher resolution. The correlation between disk size and disk luminosity
correlation demonstrates that some of the compact disks are optically thick at millimeter wavelengths.
Accepted by ApJ
http://arxiv.org/pdf/1906.10809

24

Angular momentum in bipolar outflows: dynamical evolutionary model
J.A. López-Vázquez1 , J. Cantó2 , and S. Lizano1
1

Instituto de Radioastronomı́a y Astrofı́sica, Universidad Nacional Autónoma de México, Apartado Postal 3-72, 58089
Morelia, Michoacán, México; 2 Instituto de Astronomı́a, Universidad Nacional Autónoma de México, Apartado Postal
70-264, 04510, CDMX, México
E-mail contact: j.lopez at irya.unam.mx
We model molecular outflows produced by the time dependent interaction between a stellar wind and a rotating cloudenvelope in gravitational collapse, studied by Ulrich. We consider spherical and anisotropic stellar winds. We assume
that the bipolar outflow is a thin shocked shell, with axial symmetry around the cloud rotation axis and obtain themass and momentum fluxes into the shell. We solve numerically a set of partial differential equations in space andtime, and obtain the shape of the shell, the mass surface density, the velocity field, and the angular momentum of
the material in the shell. We find that there is a critical value of the ratio between the wind and the accretion flow
momentum rates β that allows the shell to expand. As expected, the elongation of the shells increase with the stellar
wind anisotropy. In our models, the rotation velocity of the shell is the order to 0.1–0.2 km s−1 , a factor of 5–10
lower than the values measured in several sources. We compare our models with those of Wilkin and Stahler forearly evolutionary times and find that our shells have the same sizes at the pole, although we use different boundary
conditions at the equator.
Accepted by ApJ

http://arxiv.org/pdf/1906.02166

25

Observational signatures of outbursting protostars - I: From hydrodynamic simulations
to observations
Benjamin MacFarlane1 , Dimitris Stamatellos1 , Doug Johnstone2,3 , Gregory Herczeg4 , Giseon Baek5 ,Huei-Ru Vivien Chen6 , Sung-Ju Kang7 and Jeong-Eun Lee5
1

Jeremiah Horrocks Institute for Mathematics, Physics and Astronomy, University of Central Lancashire, Preston,PR1 2HE, UK; 2 NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Rd, Victoria, BC, V9E 2E7, Canada;
3
Department of Physics and Astronomy, University of Victoria, Victoria, BC, V8P 1A1, Canada; 4 Kavli Institute for
Astronomy and Astrophysics, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China; 5 School of Space
Research and Institute of Natural Sciences, Kyung Hee University, 1732 Deogyeong-daero, Giheung-gu, Yongin-si,
Gyeonggi-do 446-701, Republic of Korea; 6 Department of Physics and Institute of Astronomy, National Tsing Hua
University, Taiwan; 7 Korea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055,
Republic of Korea
E-mail contact: dstamatellos at uclan.ac.uk

Accretion onto protostars may occur in sharp bursts. Accretion bursts during the embedded phase of young protostars
are probably most intense, but can only be inferred indirectly through long-wavelength observations. We perform
radiative transfer calculations for young stellar objects (YSOs) formed in hydrodynamic simulations to predict thelong wavelength, sub-mm and mm, flux responses to episodic accretion events, taking into account heating from theyoung protostar and from the interstellar radiation field. We find that the flux increase due to episodic accretion events
is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of ∼ 570 increase in the luminosityof the young protostar leads to a flux increase of a factor of 47 at 250 µm but only a factor of 10 at 1.3 mm. Heatingfrom the interstellar radiation field may reduce further the flux increase observed at longer wavelengths. We find that23

during FU Ori-type outbursts the bolometric temperature and luminosity may incorrectly classify a source as a more
evolved YSO, due to a larger fraction of the radiation of the object being emitted at shorter wavelengths.
Accepted by MNRAS
https://arxiv.org/pdf/1906.01960

26

Observational signatures of outbursting protostars - II: Exploring a wide range of eruptive protostars
Benjamin MacFarlane1 , Dimitris Stamatellos1 , Doug Johnstone2,3 , Gregory Herczeg4 , Giseon Baek5 ,Huei-Ru Vivien Chen6 , Sung-Ju Kang7 and Jeong-Eun Lee5
1

Jeremiah Horrocks Institute for Mathematics, Physics and Astronomy, University of Central Lancashire, Preston,PR1 2HE, UK; 2 NRC Herzberg Astronomy and Astrophysics, 5071 West Saanich Rd, Victoria, BC, V9E 2E7, Canada;
3
Department of Physics and Astronomy, University of Victoria, Victoria, BC, V8P 1A1, Canada; 4 Kavli Institute for
Astronomy and Astrophysics, Peking University, Yiheyuan 5, Haidian Qu, 100871 Beijing, China; 5 School of Space
Research and Institute of Natural Sciences, Kyung Hee University, 1732 Deogyeong-daero, Giheung-gu, Yongin-si,
Gyeonggi-do 446-701, Republic of Korea; 6 Department of Physics and Institute of Astronomy, National Tsing Hua
University, Taiwan; 7 Korea Astronomy and Space Science Institute, 776 Daedeokdae-ro, Yuseong-gu, Daejeon 34055,
Republic of Korea
E-mail contact: dstamatellos at uclan.ac.uk
Young stars exhibit variability due to changes in the gas accretion rate onto them, an effect that should be quite
significant in the early stages of their formation. As protostars are embedded within their natal cloud, this variability
may only be inferred through long wavelength observations. We perform radiative transfer simulations of young stellar
objects (YSOs) formed in hydrodynamical simulations, varying the structure and luminosity properties in order toestimate the long-wavelength, sub-mm and mm, variations of their flux. We find that the flux increase due to an
outburst event depends on the protostellar structure and is more prominent at sub-mm wavelengths than at mm
wavelengths; e.g. a factor of 40 increase in the luminosity of the young protostar leads to a flux increase of a factor of
10 at 250 µm but only a factor of 2.5 at 1.3 mm. We find that the interstellar radiation field dilutes the flux increase
but that this effect may be avoided if resolution permits the monitoring of the inner regions of a YSO, where theheating is primarily due to protostellar radiation. We also confirm that the bolometric temperature and luminosity ofoutbursting protostars may result in an incorrect classification of their evolutionary stage.
Accepted by MNRAS
https://arxiv.org/pdf/1906.01966

27

Looking deep into the Rosette Nebula’s heart: the (sub)stellar content of the massive
young cluster NGC 2244
Koraljka Muzic1 , Alexander Scholz2 , Karla Peña Ramı́rez3 , Ray Jayawardhana4 , Rainer Schödel5 ,Vincent C. Geers6 , Lucas A. Cieza7 and Amelia Bayo8
1

CENTRA/SIM, Faculdade de Ciencias de Universidade de Lisboa, Portugal; 2 SUPA, School of Physics & Astronomy,
St. Andrews University, United Kingdom; 3 Centro de Astronomı́a (CITEVA), Universidad de Antofagasta, Chile;4
Department of Astronomy, Cornell University, USA; 5 Instituto de Astrofı́sica de Andalucı́a (CSIC), Spain; 6 UK
Astronomy Technology Centre, Royal Observatory Edinburgh, United Kingdom; 7 Núcleo de Astronomı́a, Facultad
de Ingenierı́a y Ciencias, Universidad Diego Portales, Santiago, Chile; 8 Instituto de Fı́sica y Astronomı́a, Universidad
de Valparaı́so, Chile
E-mail contact: kmuzic at sim.ul.pt
As part of the ongoing effort to characterize the low-mass (sub-)stellar population in a sample of massive young clusters,
we have targeted the ∼2 Myr old cluster NGC 2244. The distance to NGC 2244 from Gaia DR2 parallaxes is 1.59
kpc, with errors of 1% (statistical) and 11% (systematic). We used the FLAMINGOS-2 near-infrared camera at theGemini-South telescope for deep multi-band imaging of the central portion of the cluster (∼ 2.4 pc2 ). We determined
membership in a statistical manner, through a comparison of the cluster’s color-magnitude diagram to that of acontrol field. Masses and extinctions of the candidate members are then calculated with the help of evolutionary24

models, leading to the first Initial Mass Function (IMF) of the cluster extending into the substellar regime, with the90-percent completeness limit around 0.02 M⊙ . The IMF is well represented by a broken power-law (dN/dM ∝ M −α ),with a break at ∼0.4 M⊙ . The slope at the high mass side (0.4 − 7 M⊙ ) is α = 2.12 ± 0.08, close to the standard
Salpeter’s slope. In the low-mass range (0.02 − 0.4 M⊙), we find a slope α = 1.03 ± 0.02, which is on the high end of thetypical values obtained in nearby star-forming regions (α = 0.5 − 1.0), but still in agreement within the uncertainties.
Our results reveal no clear evidence for variations in the formation efficiency of brown dwarfs and very-low mass
stars due to the presence of OB stars, or a change in stellar densities. Our finding rules out photoevaporation andfragmentation of infalling filaments as substantial pathways for brown dwarf formation.
Accepted by ApJ
https://arxiv.org/pdf/1907.00617

28

Modeling of CoRoT and Spitzer lightcurves in NGC 2264 caused by an optically thick
warp
Erick Nagel1 and Jerome Bouvier2
1
2

Departamento de Astronomı́a, Universidad de Guanajuato,Callejón de Jalisco S/N, Guanajuato, Gto. 36240, México;
IPAG,Université Grenoble Alpes, 38000 Grenoble, France

E-mail contact: erick at astro.ugto.mx
Aims: We present an analysis of simultaneously observed CoRoT and Spitzer lightcurves for 4 systems in the stellar
forming region NGC 2264: Mon-660, Mon-811, Mon-1140 and Mon-1308. These objects share in common a highresemblance between the optical and infrared lightcurves, such that the mechanism responsible to produce them is thesame. The aim of this paper is to explain both lightcurves simultaneously with only one mechanism.
Methods: We modeled the infrared emission as coming from a warp composed of an optically thick wall and an opticallythick asymmetric disk beyond this location. We modeled the optical emission mainly by partial stellar occultation by
the warp.
Results: The magnitude amplitude of the CoRoT and Spitzer observations for all the objects can be described with theemission coming from the system components. The difference between them is the value of the disk flux compared withthe wall flux and the azimuthal variations of the former. This result points out the importance of the hydrodynamicalinteraction between the stellar magnetic field and the disk.
Conclusions: CoRoT and Spitzer lightcurves for the stellar systems Mon-660, Mon-811, Mon-1140 and Mon-1308 canbe simultaneously explained using the emission coming from an asymmetric disk and emission with stellar occultation
by an optically thick wall.
Accepted by Astronomy & Astrophysics
http://arxiv.org/pdf/1906.10110

29

Are the observed gaps in protoplanetary discs caused by growing planets?
N. Ndugu1,2 , B. Bitsch2 , E. Jurua1
1

Department of Physics, Mbarara University of Science and Technology, Mbarara, Uganda; 2 Max-Planck-Institut
for Astronomy, Königstuhl 17, D-69117 Heidelberg, Germany
E-mail contact: nndugu at must.ac.ug
Recent detailed observations of protoplanetary discs revealed a lot of sub-structures which are mostly ring-like. One
interpretation is that these rings are caused by growing planets. These potential planets are not yet opening very deep
gaps in their discs. These planets instead form small gaps in the discs to generate small pressure bumps exterior totheir orbits that stop the inflow of the largest dust particles. In the pebble accretion paradigm, this planetary mass
corresponds to the pebble isolation mass, where pebble accretion stops and efficient gas accretion starts. We perform
planet population synthesis via pebble and gas accretion including type-I and type-II migration. In the first stage of
our simulations, we investigate the conditions necessary for planets to reach the pebble isolation mass and compare
their position to the observed gaps. We find that in order to match the gap structures 2000 M⊕ in pebbles is needed,which would be only available for the most metal rich stars. We then follow the evolution of these planets for a few
My to compare the resulting population with the observed exoplanet populations. Planet formation in discs with25

these large amounts of pebbles result in mostly forming gas giants and only very little super-Earths, contradicting
observations. This leads to the conclusions that either (i) the observed discs are exceptions, (ii) not all gaps in observed
discs are caused by planets or (iii) that we miss some important ingredients in planet formation related to gas accretion
and/or planet migration.
Accepted by MNRAS
http://arxiv.org/pdf/1906.11491

30

Exploring the properties of warm and cold atomic hydrogen in the Taurus and Gemini
regions
Hiep Nguyen1,2 , J.R. Dawson1 , Min-Young Lee3,4 , Claire E. Murray5,6 , Snezana Stanimirović7 , Carl
Heiles8, M.-A. Miville-Deschênes9, Anita Petzler1
1

Department of Physics and Astronomy and MQ Research Centre in Astronomy, Astrophysics and Astrophotonics,
Macquarie University, NSW 2109, Australia; 2 Australia Telescope National Facility, CSIRO Astronomy and Space
Science, PO Box 76, Epping, NSW 1710, Australia; 3 Max-Planck-Institut fur Radioastronomie, Auf dem Hugel 69, D53121 Bonn, Germany; 4 Korea Astronomy and Space Science Institute, Daedeokdaero 776, 34055 Daejeon, Republicof Korea; 5 Department of Physics and Astronomy, Johns Hopkins University, Baltimore, MD 21218, USA; 6 NSF
Astronomy and Astrophysics Postdoctoral Fellow; 7 Department of Astronomy, University of Wisconsin-Madison, 475
North Charter Street, Madison, WI 53706, USA; 8 Department of Astronomy, University of California, Berkeley, 601
Campbell Hall 3411, Berkeley, CA 94720-3411; 9 Laboratoire AIM, Paris-Saclay, CEA/IRFU/DAp - CNRS - Université
Paris Diderot, 91191, Gif-sur-Yvette Cedex, France
E-mail contact: vanhiep.nguyen at hdr.mq.edu.au
We report Arecibo 21 cm absorption-emission observations to characterise the physical properties of neutral hydrogen(HI) in the proximity of five giant molecular clouds (GMCs): Taurus, California, Rosette, Mon OB1, and NGC 2264.
Strong HI absorption was detected toward all 79 background continuum sources in the ∼60×20 square degree region.
Gaussian decompositions were performed to estimate temperatures, optical depths and column densities of the cold andwarm neutral medium (CNM, WNM). The properties of individual CNM components are similar to those previously
observed along random Galactic sightlines and in the vicinity of GMCs, suggesting a universality of cold HI properties.
The CNM spin temperature (Ts ) histogram peaks at ∼50K. The turbulent Mach numbers of CNM vary widely, with atypical value of ∼4, indicating that their motions are supersonic. About 60% of the total HI gas is WNM, and nearly
40% of the WNM lies in thermally unstable regime 500–5000K. The observed CNM fraction is higher around GMCs
than in diffuse regions, and increases with increasing column density (NHI ) to a maximum of ∼75%. On average, the∗optically thin approximation (NHI
) underestimates the total NHI by ∼21%, but we find large regional differences in the∗relationship between NHI and the required correction factor, f = NHI /NHI
. We examine two different methods (linear∗
fit of f vs log10 NHI and uniform Ts ) to correct for opacity effects using emission data from the GALFA-HI survey.
We prefer the uniform Ts method, since the linear relationship does not produce convincing fits for all subregions.
Accepted by ApJ
http://arxiv.org/pdf/1906.06846

31

Sulfur-Bearing Species Tracing the Disk/Envelope System in the Class I Protostellar
Source Elias 29
Yoko Oya1 , Ana Lopez-Sepulcre2,3,4 , Nami Sakai5 , Yoshimasa Watanabe6,7,8 , Aya E. Higuchi5 , Tomoya
Hirota9 , Yuri Aikawa10 , Takeshi Sakai11 , Cecilia Ceccarelli3,4 , Bertrand Lefloch3,4 , Emmanuel Caux12,13 ,Charlotte Vastel12,13 , Claudine Kahane3,4 and Satoshi Yamamoto14
1

Department of Physics, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan; 2 Institut deRadioastronomie Millimétrique (IRAM), 38406, Saint Martin d’Héres, France; 3 Université Grenoble Alpes, IPAG, F38000 Grenoble, France; 4 Le Centre National de la Recherche Scientifique (CNRS), IPAG, F-38000 Grenoble, France;
5
RIKEN Cluster for Pioneering Research, 2-1, Hirosawa, Wako-shi, Saitama 351-0198, Japan; 6 Division of Physics,
Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki 305-8571, Japan; 7 Tomonaga Centerfor the History of the Universe, Faculty of Pure and Applied Sciences, University of Tsukuba, Tsukuba, Ibaraki

26

305-8571, Japan; 8 College of Engineering, Nihon University, 1 Nakagawara, Tokusada, Tamuramachi, Koriyama,
Fukushima 963-8642, Japan; 9 National Astronomical Observatory of Japan, Osawa, Mitaka, Tokyo 181-8588, Japan; 10
Department of Astronomy, The University of Tokyo, 7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan; 11 Departmentof Communication Engineering and Informatics, Graduate School of Informatics and Engineering, The University ofElectro-Communications, Chofugaoka, Chofu, Tokyo 182-8585, Japan; 12 Universite de Toulouse, UPS-OMP, F-31028
Toulouse Cedex 4, France; 13 Le Centre National de la Recherche Scientifique (CNRS), IRAP, 9 Av. Colonel Roche,
BP 44346, F-31028 Toulouse Cedex 4, France; 14 Research Center for the Early Universe, The University of Tokyo,
7-3-1, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan
E-mail contact: oya at taurus.phys.s.u-tokyo.ac.jp
We have observed the Class I protostellar source Elias 29 with Atacama Large Millimeter/submillimeter Array(ALMA). We have detected CS, SO, 34 SO, SO2 , and SiO line emissions in a compact component concentrated near
the protostar and a ridge component separated from the protostar by 4 arcsec (∼ 500 au). The former component is2
found to be abundant in SO and SO2 but deficient in CS. The abundance ratio SO/CS is as high as 3+13
−2 × 10 atthe protostar, which is even higher than that in the outflow-shocked region of L1157 B1. However, organic molecules(HCOOCH3 , CH3 OCH3 , CCH, and c-C3 H2 ) are deficient in Elias 29. We attribute the deficiency in organic molecules
and richness in SO and SO2 to the evolved nature of the source or the relatively high dust temperature ( >∼ 20 K)in the parent cloud of Elias 29. The SO and SO2 emissions trace rotation around the protostar. Assuming a highly
inclined configuration (i ≥ 65◦ ; 0◦ for a face-on configuration) and Keplerian motion for simplicity, the protostellar
mass is estimated to be (0.8 – 1.0) M⊙ . The 34 SO and SO2 emissions are asymmetric in their spectra; the blue-shiftedcomponents are weaker than the red-shifted ones. Although this may be attributed to the asymmetric molecular
distribution, other possibilities are also discussed.
Accepted by Astrophysical Journal
https://arxiv.org/pdf/1906.09724

32

JCMT BISTRO Survey observations of the Ophiuchus Molecular Cloud: Dust grainalignment properties inferred using a Ricean noise model
Kate Pattle1,2 , Shih-Ping Lai1,3 , Tetsuo Hasegawa2 , Jia-Wei Wang1 , Ray Furuya4 , Derek WardThompson5 , Pierre Bastien6 , Simon Coudé7 , Chakali Eswaraiah8 , Lapo Fanciullo3 , James diFrancesco9,10 , Thiem Hoang11,12 , Gwanjeong Kim13 , Woojin Kwon11,12 , Chang Won Lee11,12 , ShengYuan Liu3 , Tie Liu11,14 , Masafumi Matsumura15 , Takashi Onaka16 , Sarah Sadavoy17 and Archana Soam7
1

Institute for Astronomy and Department of Physics, National Tsing Hua University, Hsinchu 30013, Taiwan; 2National Astronomical Observatory of Japan, National Institutes of Natural Sciences, Osawa, Mitaka, Tokyo 181-8588,
Japan; 3 Academia Sinica Institute of Astronomy and Astrophysics, P.O. Box 23-141, Taipei 10617, Taiwan; 4 Instituteof Liberal Arts and Sciences, Tokushima University, Minami Jousanajima-machi 1-1, Tokushima 770-8502, Japan; 5Jeremiah Horrocks Institute, University of Central Lancashire, Preston PR1 2HE, UK; 6 Centre de recherche en
astrophysique du Québec & département de physique, Université de Montréal, C.P. 6128 Succ. Centre-ville, Montréal,
QC, H3C 3J7, Canada; 7 SOFIA Science Center, Universities Space Research Association, NASA Ames ResearchCenter, Moffett Field, California 94035, USA; 8 CAS Key Laboratory of FAST, National Astronomical Observatories,
Chinese Academy of Sciences, People?s Republic of China; 9 NRC Herzberg Astronomy and Astrophysics, 5071 West
Saanich Road, Victoria, BC V9E 2E7, Canada; 10 Department of Physics and Astronomy, University of Victoria,
Victoria, BC V8W 2Y2, Canada; 11 Korea Astronomy and Space Science Institute (KASI), 776 Daedeokdae-ro,
Yuseong-gu, Daejeon 34055, Republic of Korea; 12 University of Science and Technology, Korea (UST), 217 Gajeong-ro,
Yuseong-gu, Daejeon 34113, Republic of Korea; 13 Nobeyama Radio Observatory, National Astronomical Observatoryof Japan, National Institutes of Natural Sciences, Nobeyama, Minamimaki, Minamisaku, Nagano 384-1305, Japan; 14
East Asian Observatory, 660 N. A’ohōkū Place, University Park, Hilo, Hawaii 96720, USA; 15 Faculty of Education andCenter for Educational Development and Support, Kagawa University, Saiwai-cho 1-1, Takamatsu, Kagawa, 760-8522,
Japan; 16 Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku,
Tokyo 113-0033, Japan; 17 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138,
USA
E-mail contact: kpattle at gapp.nthu.edu.tw
The dependence of polarization fraction p on total intensity I in polarized submillimeter emission measurements is27

typically parameterized as p ∝ I −α (α ≤ 1), and used to infer dust grain alignment efficiency in star-forming regions,with an index α = 1 indicating near-total lack of alignment of grains with the magnetic field. In this work wedemonstrate that the non-Gaussian noise characteristics of polarization fraction may produce apparent measurementsof α ∼ 1 even in data with significant signal-to-noise in Stokes Q, U and I emission, and so with robust measurementsof polarization angle. We present a simple model demonstrating this behavior, and propose a criterion by which
well-characterized measurements of polarization fraction may be identified. We demonstrate that where our model isapplicable, α can be recovered by fitting the p−I relationship with the mean of the Rice distribution, without statisticaldebiasing of polarization fraction. We apply our model to JCMT BISTRO Survey POL-2 850µm observations of threeclumps in the Ophiuchus Molecular Cloud, finding that in the externally-illuminated Oph A region, α ≈ 0.34, whilein the more isolated Oph B and C, despite their differing star formation histories, α ∼ 0.6 − 0.7. Our results thus
suggest that dust grain alignment in dense gas is more strongly influenced by incident interstellar radiation field thanby star formation history. We further find that grains may remain aligned with the magnetic field at significantly
higher gas densities than has previously been believed, thus allowing investigation of magnetic field properties within
star-forming clumps and cores.
Accepted by ApJ
https://arxiv.org/pdf/1906.03391

33

Vortex instabilities triggered by low-mass planets in pebble-rich, inviscid protoplanetary
discs
A. Pierens1 , M.-K. Lin2 , S.N. Raymond1
1

Laboratoire d’Astrophysique de Bordeaux, CNRS and Université de Bordeaux, Allée Geoffroy St. Hilaire, 33165
Pessac, France; 2 Institute of Astronomy and Astrophysics, Academia Sinica, Taipei 10617, Taiwan
E-mail contact: arnaud.pierens at u-bordeaux.fr
In the innermost regions of protoplanerary discs, the solid-to-gas ratio can be increased considerably by a number ofprocesses, including photoevaporative and particle drift. MHD disc models also suggest the existence of a dead-zone atR<∼ 10 AU, where the regions close to the midplane remain laminar. In this context, we use two-fluid hydrodynamical
simulations to study the interaction between a low-mass planet (∼1.7 M⊕ ) on a fixed orbit and an inviscid pebblerich disc with solid-to-gas ratio ǫ ≥ 0.5. For pebbles with Stokes numbers St = 0.1, 0.5, multiple dusty vortices
are formed through the Rossby Wave Instability at the planet separatrix. Effects due to gas drag then lead to astrong enhancement in the solid-to-gas ratio, which can increase by a factor of ∼103 for marginally coupled particles
with St = 0.5. As in streaming instabilities, pebble clumps reorganize into filaments that may plausibly collapse toform planetesimals. When the planet is allowed to migrate in a MMSN disc, the vortex instability is delayed due
to migration but sets in once inward migration stops due a strong positive pebble torque. Again, particle filaments
evolving in a gap are formed in the disc while the planet undergoes an episode of outward migration. Our results
suggest that vortex instabilities triggered by low-mass planets could play an important role in forming planetesimals
in pebble-rich, inviscid discs, and may significantly modify the migration of low-mass planets. They also imply thatplanetary dust gaps may not necessarily contain planets if these migrated away.
Accepted by MNRAS
http://arxiv.org/pdf/1906.07959

34

The Duration of Star Formation in Galactic Giant Molecular Clouds. I. The Great
Nebula in Carina
Matthew S. Povich1 , Jessica T. Maldonado1 , Evan Haze Nuez1 and Thomas P. Robitaille2
1

Department of Physics and Astronomy, California State Polytechnic University, 3801 West Temple Ave, Pomona,
CA 91768 USA; 2 Headingley Enterprise and Arts Centre, Bennett Road Headingley, Leeds LS6 3HN, United Kingdom
E-mail contact: mspovich at cpp.edu
We present a novel infrared spectral energy distribution (SED) modeling methodology that uses likelihood-based
weighting of the model fitting results to construct probabilistic H–R diagrams (pHRD) for X-ray identified,intermediate-mass (2–8 M⊙ ), pre-main sequence young stellar populations. This methodology is designed specifi28

cally for application to young stellar populations suffering strong, differential extinction (∆AV > 10 mag), typical of
Galactic massive star-forming regions. We pilot this technique in the Carina Nebula Complex (CNC) by modeling the
1–8 µm SEDs of 2269 likely stellar members that exhibit no excess emission from circumstellar dust disks at 4.5 µm
or shorter wavelengths. A subset of ∼100 intermediate-mass stars in the lightly-obscured Trumpler 14 and 16 clusters
have available spectroscopic Teff , measured from the Gaia-ESO survey. We correctly identify the stellar temperature in
70% of cases, and the aggregate pHRD for all sources returns the same peak in the stellar age distribution as obtainedusing the spectroscopic Teff,S . The SED model parameter distributions of stellar mass and evolutionary age reveal
significant variation in the duration of star formation among four large-scale stellar overdensities within the CNC anda large distributed stellar population. Star formation began ∼10 Myr ago and continues to the present day, with thestar formation rate peaking < 3 Myr ago when the massive Trumpler 14 and 16 clusters formed. We make public theset of 100,000 SED models generated from standard pre-main sequence evolutionary tracks and our custom software
package for generating pHRDs and mass-age distributions from the SED fitting results.
Accepted by ApJ
https://arxiv.org/pdf/1906.01730

35

Gas accretion within the dust cavity in AB Aur
Pablo Riviére-Marichalar1 , Asunción Fuente1 , Clément Baruteau2 , Roberto Neri3 , Sandra P. TreviñoMorales4 , Andrés Carmona2 , Marcelino Agúndez5 and Rafael Bachiller1
1

Observatorio Astronómico Nacional (OAN, IGN), Calle Alfonso XII 3 E-28014 Madrid, Spain; 2 CNRS / Institut
de Recherche en Astrophysique et Planétologie, 14 avenue Edouard Belin, F-31400 Toulouse, France; 3 Institut deRadioastronomie Millimétrique, 300 rue de la Piscine, F-38406 Saint-Martin d’Héres, France; 4 Dept. of Space,
Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, SE-439 92 Onsala, Sweden;
5
Instituto de Fı́sica Fundamental, CSIC, Serrano 123, E-28006 Madrid, Spain
E-mail contact: p.riviere at oan.es
AB Aur is a Herbig Ae star hosting a well-known transitional disk. Because of its proximity and low inclination
angle, it is an excellent object to study planet formation. Our goal is to investigate the chemistry and dynamics ofthe molecular gas component in the AB Aur disk, and its relation with the prominent horseshoe shape observed incontinuum mm emission. We used the NOEMA interferometer to map with high angular resolution the J = 3–2 linesof HCO+ and HCN. By combining both, we can gain insight into the AB Aur disk structure. Chemical segregation isobserved in the AB Aur disk: HCO+ shows intense emission toward the star position, at least one bright molecularbridge within the dust cavity, and ring-like emission at larger radii, while HCN is only detected in an annular ringthat is coincident with the dust ring and presents an intense peak close to the dust trap. We use HCO+ to investigate
the gas dynamics inside the cavity. The observed bright HCO+ bridge connects the compact central source with theouter dusty ring. This bridge can be interpreted as an accretion flow from the outer ring to the inner disk/jet systemproving gas accretion through the cavity.
Accepted by MNRAS
http://arxiv.org/pdf/1906.11638

36

Pebble-driven planet formation for TRAPPIST-1 and other compact systems
Djoeke Schoonenberg1 , Beibei Liu2 , Chris W. Ormel1 , and Caroline Dorn3
1

Anton Pannekoek Institute for Astronomy, University of Amsterdam, Science Park 904, 1090 GE Amsterdam, The
Netherlands; 2 Lund Observatory, Department of Astronomy and Theoretical Physics, Lund University, Box 43, 221
00 Lund, Sweden; 3 University of Zurich, Institute of Computational Sciences, Winterthurerstrasse 190, CH-8057,
Zurich, Switzerland
E-mail contact: d.schoonenberg at uva.nl
Recently, seven Earth-sized planets were discovered around the M-dwarf star TRAPPIST-1. Thanks to transit-timing
variations, the masses and therefore the bulk densities of the planets have been constrained, suggesting that all
TRAPPIST-1 planets are consistent with water mass fractions on the order of 10%. These water fractions, as well as
the similar planet masses within the system, constitute strong constraints on the origins of the TRAPPIST-1 system.
29

In a previous work, we outlined a pebble-driven formation scenario. In this paper we investigate this formation
scenario in more detail. We used a Lagrangian smooth-particle method to model the growth and drift of pebbles andthe conversion of pebbles to planetesimals through the streaming instability. We used the N-body code MERCURY
to follow the composition of planetesimals as they grow into protoplanets by merging and accreting pebbles. This
code is adapted to account for pebble accretion, type-I migration, and gas drag. In this way, we modelled the entire
planet formation process (pertaining to planet masses and compositions, not dynamical configuration). We find thatplanetesimals form in a single, early phase of streaming instability. The initially narrow annulus of planetesimals
outside the snowline quickly broadens due to scattering. Our simulation results confirm that this formation pathway
indeed leads to similarly-sized planets and is highly efficient in turning pebbles into planets (∼50% solids-to-planets
conversion efficiency). The water content of planets resulting from our simulations is on the order of 10%, and our
results predict a ‘V-shaped’ trend in the planet water fraction with orbital distance: from relatively high (innermost
planets) to relatively low (intermediate planets) to relatively high (outermost planets).
Accepted by A&A
http://arxiv.org/pdf/1906.00669

37

The Physical and chemical structure of Sagittarius B2 – IV. Converging filaments in
the high-mass cluster forming region Sgr B2(N)A. Schwörer1 , Á. Sánchez-Monge1 , P. Schilke1 , T. Möller1 , A. Ginsburg2 , F. Meng1 , A. Schmiedeke3 ,H. S. P. Müller1, D. Lis4 and S. -L. Qin5
1
I. Physikalisches Institut, Universität zu Köln, Zülpicher Strasse 77, 50937 Köln, Germany; 2 National Radio Astronomy Observatory, 1003 Lopezville Rd., Socorro, NM 87801, USA; 3 Max-Planck-Institut für extraterrestrische Physik,
Giessenbachstrasse 1, 85748 Garching bei München, Germany; 4 Jet Propulsion Laboratory, California Institute ofTechnology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA; 5 Department of Astronomy, Yunnan University, andKey Laboratory of Astroparticle Physics of Yunnan Province, Kumming,650091, China

E-mail contact: schwoerer at ph1.uni-koeln.de
We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtainedwith ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2(north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure,
together with the presence of the massive central region, suggest that these filaments may be associated with accretion
processes. In order to derive the kinematic properties of the gas in a chemically line-rich source like SgrB2(N), we havedeveloped a new tool that stacks all the detected transition lines of any molecular species. This permits to increase
the signal-to-noise ratio of our observations and average out line blending effects, which are a common problem inline-rich regions. We derive velocity gradients along the filaments of about 20-100 km/s/pc, which are 10-100 times
larger than those typically found on larger scales (1 pc) in other star-forming regions. The mass accretion rates ofindividual filaments are about 0.05 Msun/yr, which result in a total accretion rate of 0.16 Msun/yr. Some filaments
harbor dense cores that are likely forming stars and stellar clusters. The stellar content of these dense cores is on theorder of 50% of the total mass. We conclude that the cores may merge in the center when already forming stellar
clusters but still containing a significant amount of gas, resulting in a ”damp” merger. The high density and massof the central region, combined with the presence of converging filaments with high mass, high accretion rates andembedded dense cores already forming stars, suggest that SgrB2(N) may have the potential to evolve into a superstellar cluster.
Accepted by Astronomy&Astrophysics
https://arxiv.org/pdf/1906.10979

38

A Radial Velocity Survey of Embedded Sources in the Rho Ophiuchi Cluster
Timothy Sullivan1 , Bruce A. Wilking1 , Thomas P. Greene2 , Lindsey Lisalda1 , Erika L. Gibb1 andChemeda Ejeta1
1

University of Missouri-St. Louis, 1 University Boulevard, St. Louis, MO 63121; 2 NASA Ames Research Center,
Space Science and Astrobiology Division, M.S. 245-6, Moffett Field, CA 94035

30

E-mail contact: bwilking at umsl.edu
We present the results of a radial velocity survey of young stellar objects (YSOs) in early stages of evolution in thecore of the L1688 molecular cloud. New and archival spectra obtained with four high-resolution infrared spectrographs
were analyzed using Markov chain Monte Carlo techniques that simultaneously fit for the radial velocity, Tef f , vsini,and veiling by comparison with synthetic spectra. The radial velocity distribution for 32 objects, most with Class I or
flat-spectrum spectral energy distributions, is marginally Gaussian, with a higher dispersion relative to optical surveys
at the 2σ level. When comparing the results from both proper-motion and radial velocity surveys in L1688, there isa trend for the 1D dispersions to be higher for samples of Class I/flat-spectrum YSOs that reside in the cloud corecompared to Class II/III dominated samples, which are located in the lower extinction periphery. In addition, there
is a velocity gradient along the major axis of the cloud core that appears more pronounced than that derived fromoptically visible objects at the cloud edges. If these higher dispersions for Class I/flat-spectrum objects are confirmed
by future surveys, this could imply a supervirial state for the less evolved objects in the cloud core and be a signatureof the initial collapse and rebound of the cluster as suggested by recent simulations of cluster evolution.
Accepted by The Astronomical Journal

39

Cyanopolyyne Chemistry around Massive Young Stellar Objects
Kotomi Taniguchi1,2 , Eric Herbst1,2 , Paola Caselli3 , Alec Paulive2 , Dominique M. Maffucci2 and Masao
Saito4,5
1

Department of Astronomy, University of Virginia, Charlottesville, VA 22904, USA; 2 Department of Chemistry,
University of Virginia, Charlottesville, VA 22903, USA; 3 Max-Planck-Institute for Extraterrestrial Physics (MPE),
Giessenbachstr, 1, D-85748 Garching, Germany; 4 National Astronomical Observatory of Japan (NAOJ), Osawa,
Mitaka, Tokyo 181-8588, Japan; 5 Department of Astronomical Science, School of Physical Science, SOKENDAI (The
Graduate University for Advanced Studies), Osawa, Mitaka, Tokyo 181-8588, Japan
E-mail contact: kt8pm at virginia.edu
Recent radio astronomical observations have revealed that HC5 N, the second shortest cyanopolyyne (HC2n+1 N), isabundant around some massive young stellar objects (MYSOs), which is not predicted by classical carbon-chain
chemistry. For example, the observed HC5 N abundance toward the G28.28−0.36 MYSO is higher than that in L1527,
which is one of the warm carbon chain chemistry (WCCC) sources, by more than one order of magnitude (Taniguchi
et al. 2017). In this paper, we present chemical simulations of hot-core models with a warm-up period using theastrochemical code Nautilus. We find that the cyanopolyynes are formed initially in the gas phase and accreted onto
the bulk and surface of granular ice mantles during the lukewarm phase, which occurs at 25 < T < 100 K. In slow warmup period models, the peak abundances occur as the cyanopolyynes desorb from dust grains after the temperature rises
above 100 K. The lower limits of the abundances of HC5 N, CH3 CCH, and CH3 OH observed in the G28.28−0.36 MYSO
can be reproduced in our hot-core models, after their desorption from dust grains. Moreover, previous observations
suggested chemical diversity in envelopes around different MYSOs. We discuss possible interpretations of relationships
between stages of the star-formation process and such chemical diversity, such as the different warm-up timescales.
This timescale depends not only on the mass of central stars but also on the relationship between the size of warm
regions and their infall velocity.
Accepted by The Astrophysical Journal
https://arxiv.org/pdf/1906.11296

40

Feedback from OB stars on their parent cloud: Gas exhaustion rather than gas ejection
E.J. Watkins1 , N. Peretto1 , K. Marsh2 , and G.A. Fuller3
1

School of Physics and Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK; 2 Infrared Processing
and Analysis Center, California Institute of Technology 100-22, Pasadena, CA 91125, USA; 3 Jodrell Bank Centre for
Astrophysics, School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13 9PL
E-mail contact: WatkinsEJ1 at cardiff.ac.uk
Feedback from high-mass stars shapes the ISM of galaxies and thereby impacts gas that will form future generationsof stars. However, due to our inability to track the time evolution of individual molecular clouds, quantifying the exact
31

role of feedback on their star formation history is an observationally challenging task. In the present study, we takeadvantage of the unique properties of the G316.75−00.00 high-mass star-forming ridge to determine how feedback fromO-stars impacts the dynamical stability of filaments. G316.75 is a 13.6 pc long ridge containing 18,900 M⊙ of H2 gas
which is half IR dark, half IR bright. The IR bright half has already formed 4 O-type stars over the past 2 Myr, while
the IR dark half is still quiescent. Therefore, by assuming the gas properties of the dark part represents an earlier
evolutionary stage of the bright part, we can quantify how feedback impacts these properties by contrasting them. We
use archive Herschel and molecular line data, tracing both dense (NH3 and N2 H+ ) and more diffuse (13 CO) gas, tomeasure the ratio of kinetic to gravitational energy per-unit-length, virial-ratio-per-line, across the ridge for a range ofgas volume densities. We show that despite the presence of 4 embedded O-stars, feedback cannot unbind the ridge’s
mass except for some small gas pockets near the O-stars. In fact, the virial-ratio-per-line is almost indistinguishable
for both parts of the ridge. These results are at odds with most simulations where O-star-forming clouds are dispersedby feedback within a few cloud free-fall time. We conclude that the discrepancy between such simulations and ourobservations originates from different cloud morphologies and average densities when the first O-star forms. These
results have important implications regarding, for instance, how feedback is implemented within cosmological andgalactic-scale simulations.
Accepted by A&A
http://arxiv.org/pdf/1906.09275

41

A Census Star Formation in the Outer Galaxy: the SMOG field
Elaine Winston1 , Joseph Hora1 , Robert Gutermuth2 and Volker Tolls1
1Center for Astrophysics Harvard & Smithsonian, 60 Garden St, Cambridge, MA 02138, USA; 2 Department of
Astronomy, University of Massachusetts Amherst, MA 01003, USA
E-mail contact: elaine.winston at cfa.harvard.edu

In this paper we undertake a study of the 21 square degree SMOG field, a Spitzer cryogenic mission Legacy program tomap a region of the outer Milky Way towards the Perseus and Outer spiral arms with the IRAC and MIPS instruments.
We identify 4648 YSOs across the field. Using the DBSCAN method we identify 68 clusters or aggregations of YSOs
in the region, having 8 or more members. We identify 1197 class Is, 2632 class IIs, 819 class IIIs, of which 45 are
candidate transition disk objects, utilizing the MIPS 24 photometry. The ratio of YSOs identified as members ofclusters was 2872/4648, or 62%. The ratios of class I to class II YSOs in the clusters are broadly consistent withthose found in the inner galactic and nearby Gould’s Belt young star formation regions. The clustering properties
indicate that the protostars may be more tightly bound to their natal sites than the class IIs, and the class IIIs are
generally widely distributed. We further perform an analysis of the WISE data of the SMOG field to determine how
the lower resolution and sensitivity of WISE affects the identification of YSOs as compared to Spitzer: we identify 931
YSOs using combined WISE and 2MASS photometry, 931/4648 or 20% of the total number identified with Spitzer.
Performing the same clustering analysis finds 31 clusters which reliably trace the larger associations identified with
the Spitzer data. Twelve of the clusters identified have previously measured distances from the WISE HII survey.
SEDFitter modeling of these YSOs is reported, leading to an estimation of the IMF in the aggregate of these clusters
which approximates that found in the inner galaxy, implying that the processes behind stellar mass distribution duringstar formation are not widely affected by the lower density and metallicity of the outer galaxy.
Accepted by ApJ
https://arxiv.org/pdf/1906.03267

42

HL Tau disk in HCO+ (3–2) and (1–0) with ALMA: gas density, temperature, gap, andone-arm spiral
Hsi-Wei Yen1 , Pin-Gao Gu1 , Naomi Hirano1 , Patrick M. Koch1 , Chin-Fei Lee1 , Hauyu Baobab Liu1
and Shigehisa Takakuwa2,1
1

Academia Sinica Institute of Astronomy and Astrophysics, 11F of Astro-Math Bldg, 1, Sec. 4, Roosevelt Rd, Taipei
10617, Taiwan; 2 Department of Physics and Astronomy, Graduate School of Science and Engineering, Kagoshima
University, 1-21-35 Korimoto, Kagoshima, Kagoshima 890-0065, Japan

32

E-mail contact: hwyen at asiaa.sinica.edu.tw
We present our observational results of the 1.1 mm continuum and the HCO+ (3–2) line in HL Tau at angularresolutions of 0.′′ 1 obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+
(3–2) and (1–0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+
(3–2) emission, and the stellar mass is estimated to be 2.1±0.2 M⊙ with a disk inclination angle of 47◦ . The radial
profiles of the HCO+ column density and excitation temperature are measured with the LTE analysis of the twotransitions of the HCO+ emission. An HCO+ gas gap at a radius of 30 au, where the column density drops by a factor
of 4–8, is found in the HCO+ column density profile, coincident with the dust gap traced by the continuum emission.
No other clear HCO+ gas gap is seen. This HCO+ gas gap can be opened by a planet with mass of 0.5–0.8 MJ , whichis comparable to the planet mass adopted in numerical simulations to form the dust gap at the same radius in the
HL Tau disk. In addition to the disk component, a one-arm spiral with a length of ∼3′′ (520 au) stretching out fromthe inner disk is observed in the HCO+ (3–2) emission. The observed velocity structures along the spiral suggest an
infalling and rotating gas stream toward the inner disk.
Accepted by ApJ
https://arxiv.org/pdf/1906.05535

43

Structure, kinematics, and ages of the young stellar populations in the Orion region
E. Zari1 , A.G.A. Brown1 , and P.T. de Zeeuw1
1

Leiden Observatory, Leiden University, Niels Bohrweg 2, 2333 CA Leiden, the Netherlands

E-mail contact: zariem at strw.leidenuniv.nl
We present a study of the three dimensional structure, kinematics, and age distribution of the Orion OB association,
based on the second data release of the Gaia satellite (Gaia DR2). Our goal is to obtain a complete picture of thestar formation history of the Orion complex and to relate our findings to theories of sequential and triggered starformation. We select the Orion population with simple photometric criteria, and we construct a three dimensional
map in galactic Cartesian coordinates to study the physical arrangement of the stellar clusters in the Orion region.
The map shows structures that extend for roughly 150 pc along the line of sight, divided in multiple sub-clusters.
We separate different groups by using the density based clustering algorithm DBSCAN. We study the kinematic
properties of all the groups found by DBSCAN first by inspecting their proper motion distribution, and then by
applying a kinematic modelling code based on an iterative maximum likelihood approach, which we use to derive theirmean velocity, velocity dispersion and isotropic expansion. By using an isochrone fitting procedure we provide ages
and extinction values for all the groups. We confirm the presence of an old population (∼15 Myr) towards the 25
Ori region, and we find that groups with ages of 12–15Myr are present also towards the Belt region. We notice thepresence of a population of ∼10 Myr also in front of the Orion A molecular cloud. Our findings suggest that starformation in Orion does not follow a simple sequential scenario, but instead consists of multiple events, which caused
kinematic and physical sub-structure. To fully explain the detailed sequence of events, specific simulations and further
radial velocity data are needed.
Accepted by A&A
http://arxiv.org/pdf/1906.07002

44

Probing the initial conditions of high-mass star formation.
triggered star formation

III. Fragmentation and

Chuan-Peng Zhang1,2,3 , Timea Csengeri3 , Friedrich Wyrowski3 , Guang-Xing Li4,5 , Thushara Pillai3 ,Karl M. Menten3 , Jennifer Hatchell6 , Mark A. Thompson7 , and Michele R. Pestalozzi8
1National Astronomical Observatories, Chinese Academy of Sciences, 100101 Beijing, P.R. China; 2 Max-PlanckInstitut für Astronomie, Königstuhl 17, D-69117 Heidelberg, Germany; 3 Max-Planck-Institut für Radioastronomie,
Auf dem Hügel 69, D-53121 Bonn, Germany; 4 University Observatory Munich, Scheinerstrasse 1, D-81679 München,
Germany; 5 South-Western Institute for Astronomy Research, Yunnan University, Kunming, 650500 Yunnan, P.R.
China; 6 School of Physics, University of Exeter, Stocker Rd, Exeter EX4 4QL, United Kingdom; 7 Centre for
Astrophysics Research, University of Hertfordshire, College Lane, Hatfield, Herts AL10 9AB, United Kingdom; 833

Istituto di Fisica dello Spazio Interplanetario - INAF, via Fosso del Cavaliere 100, 00133, Roma, Italy
E-mail contact: cpzhang at nao.cas.cn
Massive clumps tend to fragment into clusters of cores and condensations, some of which form high-mass stars. In thiswork, we study the structure of massive clumps at different scales, analyze the fragmentation process, and investigate
the possibility that star formation is triggered by nearby HII regions. We present a high angular resolution studyof a sample of 8 massive proto-cluster clumps. Combining infrared data, we use few-arcsecond resolution radio- andmillimeter interferometric data to study their fragmentation and evolution. Our sample is unique in the sense that
all the clumps have neighboring HII regions. Taking advantage of that, we test triggered star formation using a novel
method where we study the alignment of the centres of mass traced by dust emission at multiple scales. The eightmassive clumps have masses ranging from 228 to 2279 M⊙ . The brightest compact structures within infrared brightclumps are typically associated with embedded compact radio continuum sources. The smaller scale structures ofReff ∼ 0.02 pc observed within each clump are mostly gravitationally bound and massive enough to form at least a
B3-B0 type star. Many condensations have masses larger than 8 M⊙ at small scale of Reff ∼ 0.02 pc. Although theclumps are mostly infrared quiet, the dynamical movements are active at clump scale (∼1 pc). We studied the spatial
distribution of the gas conditions detected at different scales. For some sources we find hints of external triggering,whereas for others we find no significant pattern that indicates triggering is dynamically unimportant. This probably
indicates that the different clumps go through different evolutionary paths. In this respect, studies with larger samples
are highly desired.
Accepted by A&A
http://arxiv.org/pdf/1906.03216

34

Dissertation Abstracts

The formation, evolution, and survivability of discs around young
binary stars
Rajika Kuruwita

Australian National University
Mt Stromlo Observatory, Cotter Rd, Weston Creek, ACT 2611, Australia
Address as of April 2019: Centre for Star and Planet Formation, University of Copenhagen, Oster Voldgade 5-7
DK-1350, Copenhagen, Denmark
Electronic mail: rajikakuruwita at gmail.com
Ph.D dissertation directed by: Christoph FederrathPh.D degree awarded: July 2019
In the age of the Kepler space telescope and other exoplanet finding missions, a variety of exotic planets have been
discovered. Some of these planets have been found to be in binary star systems — systems which have historically been
overlooked in planet formation models. This is due to the single star scenario being simpler to model than binaries,
as well our anthropocentric bias towards single stars like our Sun. However, planet formation around binary stars is
an important topic because a large fraction (∼50%) of stars form in binary systems.
In this thesis I investigated the physics that influences the creation, stability, and survivability of discs around binarystars with the broad understanding that the longer the lifetime of a disc (around a single or binary star) the higher
the likelihood of producing planets.
The theoretical work of this thesis was conducted using the ideal magnetohydrodynamical numerical simulation program FLASH. I simulated the collapse of molecular cores until the formation of protostars and followed the early
evolution of these systems. For the first theoretical project I investigated the influence that binarity had on the global
evolution of a young stellar system. This included studying mechanisms such as accretion, jets and outflows, anddynamical interactions. I found that binary stars produce weaker outflows when considering the transport of mass,linear momentum, and angular momentum. For the second theoretical project I investigated the formation of discs inbinary stars with the inclusion of turbulence in the initial conditions. I found that the turbulence helped to build large
circumbinary discs which restructured the magnetic fields for efficient outflow launching, but too much turbulence may
also disrupt this organisation of magnetic fields. Given the environment where binary stars form (turbulent molecularcores), it appears that the formation of circumbinary discs should be common place, however circumstellar discs could
also be destroyed quickly in these same environments.
My observational work aimed to determine the typical survivability of discs around binary star systems. This work
was carried out by using the Wide Field Spectrograph (WiFeS) on the Australian National University 2.3 m Telescope35

to search for radial velocity variation in disc-bearing members of the 11 Myr and 17 Myr old star-forming regions
Upper Scorpius and Upper Centaurus-Lupus. I found that the binary fraction of disc-bearing stars in these regions
do not differ significantly from the field star binary fraction. I hypothesised that this is due to two competing factors:
circumstellar discs are disrupted by companions and are dispersed quickly, but circumbinary discs are more common
than equivalently sized discs around single stars. These results suggest that the typical lifetimes of discs in single andbinary stars are comparable.
Overall, I found that in some scenarios binary stars may produce hostile environments for planet formation via thedestruction of circumstellar discs, but the formation of large circumbinary discs is likely to be a common occurrence.
This suggests that planet formation is as likely around binary stars as single stars. Therefore, planet formation aroundbinary stars needs to be considered to understand overall planet formation.

Moving ... ??
If you move or your e-mail address changes, please
send the editor your new address. If the Newsletterbounces back from an address for three consecutive
months, the address is deleted from the mailing list.

Abstract submission deadline
The deadline for submitting abstracts and othersubmissions is the first day of the month. Abstracts
submitted after the deadline will appear in the following month’s issue.

36

New Jobs

Postdoctoral position in dusty turbulence
Applications are invited for a 2-year Postdoctoral Research Fellow to work with Dr. Benoı̂t Commerçon and Dr.
Guillaume Laibe at the Centre de Recherche Astrophysique de Lyon (http://cral.univ-lyon1.fr) in the AstroENS
team based at Ecole Normale Supérieure de Lyon on characterising and modelling astrophysical dusty turbulence inthe context of star and planet formation.
This includes any of the following aspects:- Numerical simulations of dusty turbulence at different scales with the Adaptative Mesh-Refinement grid-based or
Smoothed Particle Hydrodynamics methods.
- Modelling the statistics of interstellar dusty turbulence.
- Interpreting the experiments of dusty turbulence conducted at the Physics Department at ENS de Lyon (UMR5672).
The successful candidate should have a Ph.D. in Astrophysics, Physics or Applied mathematics in relation withhydrodynamics with experience in turbulence, high-performance computing or dust-gas flows.
By August 16th 2019, enthusiastic candidates should submit on the web platform on the web platform indicatedjust below i) a CV including a publication list, ii) a brief statement of research interests and iii) arrange for threeletters of reference to be provided separately (sent to benoit.commercon@ens-lyon.fr and guillaume.laibe@ens-lyon.fr).
Applications which will not be submitted on the web platform can not be considered.
Web platform for applications:
https://emploi.cnrs.fr/Offres/CDD/UMR5574-BENCOM-001/Default.aspx?lang=EN
Complementary information:
Informal enquiries: benoit.commercon@ens-lyon.fr guillaume.laibe@ens-lyon.fr
Starting date: Between Fall 2019 and Spring 2020
Salary: The initial monthly net salary is around 2000 and is comensurate with experience.
Included Benefits: National medical insurance, maternity/paternity leaves, family supplement for children, participation to public transport fees, pension contributions. School is free in France for children above 3. The Universit de
Lyon provides support for settling in Lyon and France.
Miscellaneous:- Since 2016, the École Normale Supérieure de Lyon appears in the world’s Top 5 small universities of the Times
Higher Education ranking.
- Lyon is a vibrant city with high quality of life (e.g. https://www.nytimes.com/2006/07/02/travel/02goingto.html)- The project is funded by the IDEXLyon project (contract n?ANR-16-IDEX-0005) under University of Lyon auspices.

37

Meetings

European Conference on Laboratory Astrophysics

Linking Dust, Ice and Gas in Space19 - 24 April 2020, Anacapri, Italy
The interplay between ubiquitous dust, ice and gas in space knits an interesting tale from collapsing interstellar clouds
to the formation of new stars, planets, moons and comets. Along this path, the formation of complex organic molecules
necessary to construct the building blocks of life brings us a step closer to the understanding of the evolution of life.
The advancement in the understanding of these vast intricacies of space lies in the development of varied laboratory
techniques and astrophysical modelling in close connection with space exploration and astronomical observations.
We are living a very exciting time at the interplay between current and future missions/facilities. Return of extraterrestrial samples is one of the drivers of these missions (e.g. Hayabusa 2, OSIRIS-Rex) together with a more andmore detailed exploration of the Solar System planets, moons and other bodies (e.g. Martian rovers ExoMars and Mars
2020, BepiColumbo for Mercury, the future JUICEand Europa Clipper missions for Jupiter and its moon Europa).
From the astronomical side, new perspectives arise from progress in high angular resolution observations from theoptical/infrared (VLT and the coming JWST and ELT) to the sub-mm/mm range (ALMA, NOEMA) to the midradio frequencies (SKA) together with a better coverage of the electromagnetic spectrum (e.g. ATHENA for X-rays)
and the discovery of new exo-worlds (e.g. PLATO, ARIEL).
The ECLA 2020 conference will provide the opportunity for interdisciplinary exchanges between scientists activein different research fields, which range from astronomy to geology and from chemistry to instrumentation. These
interactions will foster collaborations that will result in better scientific exploitation of missions/facilities in planetology
and astronomy in relation with the major questions on our origins.
http://www.frcongressi.it/ecla2020/”Cool Stars, Stellar Systems, and the Sun” (CS21) first announcement
22 - 26 June 2020, Toulouse, France
We are pleased to announce that the international conference ”Cool Stars, Stellar Systems, and the Sun (CS21)” will
be held in Toulouse, France in June 2020, 22nd to 26th. The pre-registration is open. Please go to:https://coolstars21.github.io (see website for details)
Please note the following important dates:
June 17th 2019 : opening of pre-registration
Early September 2019: Call for Splinter sessions
Late November 2019: Deadline for submission of Splinter Sessions
Mid December 2019: list of Splinter sessions selected
January 2020 : Early bird Registration fee and Abstract soumission to main conference and splinters
Early February 2020: Deadline for travel grant support request
Early March 2020 : Normal registration fee
Mid April 2020 : announcement of program and selected grantees
38

Early May 2020: Late registration fee
Please mark your calendar, looking forward to seeing you all in Toulouse, next Summer. You can contact us at:
cs21@irap.omp.eu
A.S. Brun, J. Bouvier and P. Petit for the SOC/LOC

The Physics of Star Formation: From Stellar Cores to Galactic Scales
June 29 - July 3 2020 - Lyon, France
http://staratlyon.univ-lyon1.fr
We are happy to announce the conference “The Physics of Star Formation: From Stellar Cores to Galactic Scales”
which will take place from June 29 to July 3 2020 in Lyon.
Pre-registration is open on the website of the conference.
Looking forward to seeing you in Lyon in early summer 2020!
Rationale
Stars are the visible building blocks of the cosmic structures, and as such are essential for our understanding of the
Universe and the physical processes that govern its evolution. Star formation lies at the centre of the interstellar
matter cycle that drives cosmic evolution. Although star formation is recognised to be central in a wide range ofastrophysical processes (galaxy formation and evolution, molecular cloud evolution and destruction, planet formation,and potentially the development of life), it is one of the least understood processes in the Universe.
The aim of this conference is to bring together the communities related to star formation through the physical scales
and cosmic times. This topic is particularly timely for defining major explorations of the star forming Universe to be
carried out in the next decade with new instruments designed for observations and exascale HPC. The state-of-the-artof the three pillars of modern astrophysics, theory, observation, and numerical simulations, will be represented. The
program will be divided into eight sessions. Each session will be covered by invited reviews, invited presentations aswell a large fraction of contributed talks.
Invited Reviews- Nate Bastian (Liverpool John Moores University, UK)- Avishai Dekel (Hebrew University of Jerusalem, Israel)- James Di Francesco (NRCC, Canada)- Reinhard Genzel (MPE Garching, Germany)- Patrick Hennebelle (CEA Saclay, France)- Philip Hopkins (CalTech, USA)
- Shu-ichiro Inutsuka (University of Nagoya, Japan)- Diederik Kruijssen (University of Heidelberg, Germany)- Mark Krumholz (Australian National University, Australia)- Adam Leroy (Ohio State University, USA)- Anaëlle Maury (CEA Saclay, France)- Raffaella Schneider (University of Rome, Italy)- Kengo Tomida (University of Osaka, Japan)- Yusuke Tsukamoto (University of Kagoshima, Japan)- Stefanie Walch (University of Cologne, Germany)
Topics
1/ Physic and numerics for star formation
2/ Isolated star formation: from brown dwarfs to massive stars
3/ Cluster formation
4/ Protostellar disc formation and evolution
5/ Initial Mass Function and Star Formation Rate39

6/ Stellar feedback on molecular clouds
7/ Galactic disc formation and evolution
8/ Star formation across redshift
Scientific Organization Committee
Philippe André, Avishai Dekel, Patrick Hennebelle, Shu-ichiro Inutsuka, Helen Kirk, Ralf Klessen, Steven Longmore,Frédérique Motte, Christopher McKee, Nami Sakai, Joe Silk, Simon White.
Local Organization Committee
Gilles Cabrier (co-chair), Benoı̂t Commercon (co-chair), Congresss Team of the University of Lyon 1
E-mail contact: staratlyon@univ-lyon1.fr

Planet Formation: From Dust Coagulation to Final Orbit Assembly
1 - 26 June 2020, Munich, Germany
Organizers:
Man Hoi Lee (Department of Earth Sciences, The University of Hong Kong)
Nader Haghighipour (Institute for Astronomy, University of Hawaii-Manoa)
Soko Matsumura (Department of Physics, University of Dundee)
Hilke Schlichting, (Department of Earth, Space, and Planetary Sciences, UCLA)
Dear colleagues,
It gives us great pleasure to announce the 2020 Planet Formation program at the Munich Institute for Astro- andParticle Physics (MIAPP) in Garching, Germany. The program will be held on June 1-26, 2020. This four-weekprogram aims to bring together experts with complementary expertise in observation and theory in solar systemand extrasolar planetary astronomy to assess the current status of planet formation models, highlight problems ineach formation stage, and explore the possibility of developing comprehensive models that can be applied to different
planetary systems. Please see the program website for more details: http://www.munich-iapp.de/planetformation
The structure of the program will be informal. The main goal of the program is to create an environment that will
facilitate collaborations and new initiatives. There will be daily gathering for discussing specific topics, new results,
and brain storming.
Each attendee will be provided with an office and Internet access. Financial support at the rate of EUR 80 per
day will be provided for accommodation and local expenses. Additional financial support for attendees with family and children and for graduate students is available. Details with regard to financial support can be found at
http://www.munich-iapp.de/home/information-for-participants/financial-support/
In order to attend, it is necessary to apply to the program.
Please go to the program website(http://www.munich-iapp.de/planetformation), scroll down to the bottom of the page, and click on ”Register
here” (that is actually the link to apply for the program). Please note that the application deadline is September 1,
2019. Please also note that MIAPP requires attendance for at least two weeks and that the program can accommodate
only 45 attendees in each week. So please apply early.
We invite and encourage you to apply. We also request that you kindly pass this information to your colleagues,
collaborators, students and postdocs who are involved and/or interested in planet formation.
Looking forward to receiving your applications.
Best regards,
Man Hoi, Nader, Soko, and Hilke
E-mail contact: nader@ifa.hawaii.edu
http://www.munich-iapp.de/planetformation

40

Summary of Upcoming Meetings

Great Barriers in Planet Formation
21 - 26 July 2019 Palm Cove, Australia
https://dustbusters.bitbucket.io/great-barriers-2019/
Summer School Protoplanetary Disks and Planet Formation
5 - 9 August 2019 Copenhagen, Denmark
http://nbia.nbi.ku.dk/nbia-school-2019
Orion Uncovered
26 - 30 August 2019 Leiden, The Netherlands
https://sites.google.com/view/OrionLeiden2019
Understanding the Nearby Star-forming Universe with JWST
26 - 30 Aug 2019 Courmayeur, Italy
http://www.stsci.edu/institute/conference/unsfjwst2019
Celebrating the first 40 Years of Alexander Tielens’ Contribution to Science: The Physics and Chemistry of the ISM
2 - 6 september 2019, Avignon, France https://tielens2019.sciencesconf.org
From Gas to Stars: The Links between Massive Star and Star Cluster Formation
16-20 September 2019 York, UK
https://starformmapper.org/final-conference/
Crete III - Through dark lanes to new stars Celebrating the career of Prof. Charles Lada
23 - 27 September 2019 Crete, Greece
http://crete3.org
The UX Ori type stars and related topics
30 September - 4 October 2019 St. Petersburg, Russia
http://uxors-2019.craocrimea.ru
First Stars VI
1 - 6 March 2020 Concepcion, Chile
http://www.astro.udec.cl/FirstStarsVI/
Linking Dust, Ice, and Gas in Space19 - 24 April 2020, Capri, Italy
http://frcongressi.net/ecla2020.meet
COOL STARS 21: Cambridge Workshop on Cool Stars, Stellar Systems and the Sun
22 -26 June 2020, Toulouse, France
https://coolstars21.github.io/
The Physics of Star Formation: From Stellar Cores to Galactic Scales
29 June - 3 July 2020 Lyon, France https://cral.univ-lyon1.fr

41

Short Announcements

GaussPy+: A fully automated Gaussian decomposition package for emission line
spectraOur understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity
structure of emission line observations. However, the accumulation of ever increasing datasets, such as large Galactic
plane surveys, requires tools to simplify the enormous quantity of information present in these observations andfacilitate the extraction of interesting physical quantities.
With this goal in mind, we have developed GaussPy+, a fully automated Gaussian decomposition package that canbe applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package
upon the existing GaussPy algorithm (Lindner et al., 2015, AJ, 149, 138), but added many new functionalities andconvenience functions and significantly improved its performance for low signal-to-noise data.
The GaussPy+ package was designed to deal with different populations of emission peaks, such as a mixture of narrowand broad line shapes that are expected from contributions of the cold and warm neutral medium in HI observations,
respectively. GaussPy+ also contains automated spatial refitting routines that can add spatial coherence to thedecomposition results by refitting spectra based on neighbouring fit solutions, which allows, for instance, a better
decision on how to decompose blended emission peaks.
GaussPy+ should be of assistance to anyone interested in extracting and analyzing the detailed velocity structure fromemission line observations that can be well fitted with Gaussian components to, for example, study the turbulence
properties of molecular clouds or gain information on the (integrated) intensity and velocity dispersion values present
in the dataset, whenever simpler techniques such as moment analysis cannot be applied straightforwardly.
The
GaussPy+
packageiswrittenhttps://github.com/mriener/gausspyplus

inPython3

andfreelyavailable

on

GitHub:For a technical description of the GaussPy+ package and a discussion about decomposition results for synthetic spectra
and a test field from the Galactic Ring Survey (Jackson et al., 2006, ApJS, 163, 145) please see:
Riener et al., 2019, accepted for publication in A&A https://arxiv.org/pdf/1906.10506
Manuel Riener (manuel.riener at mpia-hd.mpg.de)

SOFIA Cycle 8 Calls for Proposals / Funding opportunities
The
SOFIA
Observatory,the
only
aircraft-based
facility
offering
anaccesstothe
4-600µm
wavelength
range,
hasrecently
releasedcalls
for
proposals(https://www.sofia.usra.edu/science/proposing-and-observing/proposal-calls/cycle-8) for the Cycle 8 period, with a deadline of September 6, 2019 (9 p.m. PDT). Anyone in the astronomy scientific community iswelcome to apply for time, and US-based proposers are also eligible for funding opportunities listed below.
In this observing Cycle, we are offering:- up to 300 hours of observations for Regular Proposals and up to 400 hours for Legacy Proposals (1-4 proposals up
to 200 h of observations each, spread over 2 observing cycles)- one Southern Hemisphere deployment to New Zealand with two instruments- improved mapping modes on instruments FORCAST (spectral mapping) and HAWC+ (on-the-fly polarimetry mapping)- new filters for FIFI-LS improving the sensitivity at the [OIII] 52µm line and for HAWC+ making band B (63µm)
available42

Also note the following funding opportunities for US-based proposers:- Up to $3M for Regular Proposals- Up to $2M per year for Legacy Proposals- $300k available through our new archival research proposal program- for proposals which are central to a PhD thesis, additional funding can be requested through the Thesis-enabling
Program (up to two years of graduate student funding)In addition, the SOFIA Science Center provides specialist support for proposal preparation and data analysis via email, phone, and otherwise. Please feel free to contact the Science Center through our HelpDesk:
sofia help@sofia.usra.edu. We will also hold a webinar on August 9th to provide support for proposers on how todesign effective observing proposals and use SOFIA proposal tools (more information to come on the SOFIA website
https://www.sofia.usra.edu/science).
The selection results will be announced in December of 2019. Good luck!

43