Search

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

We present a study of the molecular gas in early-mid stage major-mergers, with a sample of 43 major-merger galaxy pairs selected from the Mapping Nearby Galaxies at Apache Point Observatory (MaNGA) survey and a control sample of 195 isolated galaxies selected from the xCOLD GASS survey. Adopting kinematic asymmetry as a new effective indicator to describe the merger stage, we aim to study the role of molecular gas in the merger-induced star formation enhancement along the merger sequence of galaxy pairs. We obtain the molecular gas properties from CO observations with the James Clerk Maxwell Telescope (JCMT), Institut de Radioastronomie Milimetrique (IRAM) 30-m telescope, and the MASCOT survey. Using these data, we investigate the differences in molecular gas fraction ($f_{\rm H_{2}}$), star formation rate (SFR), star formation efficiency (SFE), molecular-to-atomic gas ratio ($M_{\rm H_{2}}/M_{\rm HI}$), total gas fraction ($f_{\rm gas}$), and the star formation efficiency of total gas (${\rm SFE_{gas}}$) between the pair and control samples. In the full pair sample, our results suggest the $f_{\rm H_{2}}$ of paired galaxies is significantly enhanced, while the SFE is comparable to that of isolated galaxies. We detect significantly increased $f_{\rm H_{2}}$ and $M_{\rm H_{2}}/M_{\rm HI}$ in paired galaxies at the pericenter stage, indicating an accelerated transition from atomic gas to molecular gas due to interactions. Our results indicate that the elevation of $f_{\rm H_{2}}$ plays a major role in the enhancement of global SFR in paired galaxies at the pericenter stage, while the contribution of enhanced SFE in specific regions requires further explorations through spatially resolved observations of a larger sample spanning a wide range of merger stages., Comment: 25 pages, 12 figures, 5 tables, accepted for publication in ApJS

Filamentary structures are ubiquitously found in high-mass star-forming clouds. To investigate the relationship between filaments and star formation, we carry out the INFANT (INvestigations of massive Filaments ANd sTar formation) survey, a multi-scale, multi-wavelength survey of massive filamentary clouds with ALMA band 3/band 6 and VLA K band. In this first paper, we present the ALMA band 6 continuum observations toward a sample of 8 high-mass star forming filaments. We covered each target with approximately rectangular mosaic field of view with two 12-m array configurations, achieving an angular resolution of $\sim$0.6" (2700 AU at 4.5 kpc) and a continuum rms of $\sim$0.1 mJy/beam ($\sim$0.06 Msun in gas mass assuming 15 K). We identify cores using the getsf and astrodendro and find the former is more robust in terms of both identification and measuring flux densities. We identify in total 183 dense cores (15--36 cores in each cloud) and classify their star formation states via outflow and warm gas tracers. The protostellar cores are statistically more massive than the prestellar cores, possibly indicating further accretion onto cores after formation of protostars. For the high-mass end ($M_\text{core}$ $>$ 1.5 Msun) of the core mass function (CMF) we derive a power-law index of $-$1.15 $\pm$ 0.12 for the whole sample, and $-$1.70 $\pm$ 0.25 for the prestellar population. We also find a steepening trend in CMF with cloud evolution ($-$0.89 $\pm$ 0.15 for the young group v.s. $-$1.44 $\pm$ 0.25 for the evolved group) and discuss its implication for cluster formation., Comment: 25 pages, 8 figures, accepted for ApJ

The variations of elemental abundance and their ratios along the Galactocentric radius result from the chemical evolution of the Milky Way disks. The $\rm ^{12}C/^{13}C$ ratio in particular is often used as a proxy to determine other isotopic ratios, such as $\rm ^{16}O/^{18}O$ and $\rm ^{14}N/^{15}N$. Measurements of $\rm ^{12}CN$ and $\rm ^{13}CN$ (or $\rm C^{15}N$) -- with their optical depths corrected via their hyper-fine structure lines -- have traditionally been exploited to constrain the Galactocentric gradients of the CNO isotopic ratios. Such methods typically make several simplifying assumptions (e.g. a filling factor of unity, the Rayleigh-Jeans approximation, and the neglect of the cosmic microwave background) while adopting a single average gas phase. However, these simplifications introduce significant biases to the measured $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$. We demonstrate that exploiting the optically thin satellite lines of $\rm ^{12}CN$ constitutes a more reliable new method to derive $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$ from CN isotopologues. We apply this satellite-line method to new IRAM 30-m observations of $\rm ^{12}CN$, $\rm ^{13}CN$, and $\rm C^{15}N$ $N=1\to0$ towards 15 metal-poor molecular clouds in the Galactic outer disk ($R_{\rm gc} > $ 12 kpc), supplemented by data from the literature. After updating their Galactocentric distances, we find that $\rm ^{12}C/^{13}C$ and $\rm ^{14}N/^{15}N$ gradients are in good agreement with those derived using independent optically thin molecular tracers, even in regions with the lowest metallicities. We therefore recommend using optically thin tracers for Galactic and extragalactic CNO isotopic measurements, which avoids the biases associated with the traditional method., Comment: 41 pages, 29 figures, accepted by MNRAS. Meeting materials related to this work at https://box.nju.edu.cn/d/5035a574e236408eab94/

We studied the unique kinematic properties in massive filament G352.63-1.07 at $10^3$-AU spatial scale with the dense molecular tracers observed with the Atacama Large Millimeter/submillimeter Array (ALMA). We find the central massive core M1 (12 $M_\odot$) being separated from the surrounding filament with a velocity difference of $v- {v}_{sys}=-2$ km/s and a transverse separation within 3 arcsec. Meanwhile, as shown in multiple dense-gas tracers, M1 has a spatial extension closely aligned with the main filament and is connected to the filament towards its both ends. M1 thus represents a very beginning state for a massive young star-forming core escaping from the parental filament, within a time scale of $\sim 4000$ years. Based on its kinetic energy ($3.5\times10^{44}$ erg), the core escape is unlikely solely due to the original filament motion or magnetic field, but requires more energetic events such as a rapid intense anisotropic collapse. The released energy also seems to noticeably increase the environmental turbulence. This may help the filament to become stabilized again., Comment: 15 pages, 6 figures, accepted for publication in the Astrophysical Journal

Studying the physical and chemical properties of cold and dense molecular clouds is crucial for the understanding of how stars form. Under the typical conditions of infrared dark clouds, CO is removed from the gas phase and trapped on to the surface of dust grains by the so-called depletion process. This suggests that the CO depletion factor ($f_{\rm D}$) can be a useful chemical indicator for identifying cold and dense regions (i.e., prestellar cores). We have used the 1.3 mm continuum and C$^{18}$O(2-1) data observed at the resolution of $\sim$5000 au in the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES) to construct averaged maps of $f_{\rm D}$ in twelve clumps to characterise the earliest stages of the high-mass star formation process. The average $f_{\rm D}$ determined for 277 of the 294 ASHES cores follows an unexpected increase from the prestellar to the protostellar stage. If we exclude the temperature effect due to the slight variations in the NH$_3$ kinetic temperature among different cores, we explain this result as a dependence primarily on the average gas density, which increases in cores where protostellar conditions prevail. This shows that $f_{\rm D}$ determined in high-mass star-forming regions at the core scale is insufficient to distinguish among prestellar and protostellar conditions for the individual cores, and should be complemented by information provided by additional tracers. However, we confirm that the clump-averaged $f_{\rm D}$ values correlates with the luminosity-to-mass ratio of each source, which is known to trace the evolution of the star formation process., Comment: 20 pages, 10 figures, 2 tables. Accepted for publication in The Astrophysical Journal (ApJ)

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed., Comment: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

We present a study of narrow filaments toward a massive infrared dark cloud, NGC 6334S, using the Atacama Large Millimeter/submillimeter Array (ALMA). Thirteen gas filaments are identified using the H$^{13}$CO$^{+}$ line, while a single continuum filament is revealed by the continuum emission. The filaments present a compact radial distribution with a median filament width of $\sim$0.04 pc narrower than the previously proposed `quasi-universal' 0.1~pc filament width. The higher spatial resolution observations and higher-density gas tracer tend to identify even narrower and lower mass filaments. The filament widths are roughly twice the size of embedded cores. The gas filaments are largely supported by thermal motions. The nonthermal motions are predominantly subsonic and transonic in both identified gas filaments and embedded cores, which may imply that stars are likely born in environments of low turbulence. A fraction of embedded objects show a narrower velocity dispersion compared with their corresponding natal filaments, which may indicate that the turbulent dissipation is taking place in these embedded cores. The physical properties (mass, mass per unit length, gas kinematics, and width) of gas filaments are analogous to those of narrow filaments found in low- to high-mass star-forming regions. The more evolved sources are found to be farther away from the filaments, a situation that may have resulted from the relative motions between the YSOs and their natal filaments., Comment: 28 pages, 12 figures, 1 table. Accepted for publication in ApJ

With a mass of $\sim$1000 $M_\odot$ and a surface density of $\sim$0.5 g cm$^{-2}$, G023.477+0.114 also known as IRDC 18310-4 is an infrared dark cloud (IRDC) that has the potential to form high-mass stars and has been recognized as a promising prestellar clump candidate. To characterize the early stages of high-mass star formation, we have observed G023.477+0.114 as part of the ALMA Survey of 70 $\mu$m Dark High-mass Clumps in Early Stages (ASHES). We have conducted $\sim$1."2 resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm in dust continuum and molecular line emission. We identified 11 cores, whose masses range from 1.1 $M_\odot$ to 19.0 $M_\odot$. Ignoring magnetic fields, the virial parameters of the cores are below unity, implying that the cores are gravitationally bound. However, when magnetic fields are included, the prestellar cores are close to virial equilibrium, while the protostellar cores remain sub-virialized. Star formation activity has already started in this clump. Four collimated outflows are detected in CO and SiO. H$_2$CO and CH$_3$OH emission coincide with the high-velocity components seen in the CO and SiO emission. The outflows are randomly oriented for the natal filament and the magnetic field. The position-velocity diagrams suggest that episodic mass ejection has already begun even in this very early phase of protostellar formation. The masses of the identified cores are comparable to the expected maximum stellar mass that this IRDC could form (8-19 $M_\odot$). We explore two possibilities on how IRDC G023.477+0.114 could eventually form high-mass stars in the context of theoretical scenarios., Comment: Accepted for publication in ApJ (September 1, 2021). 33 pages, 20 figures, and 5 tables

We present the results of on-the-fly mapping observations of 44 fields containing 107 SCUBA-2 cores in the emission lines of molecules, N$_2$H$^+$, HC$_3$N, and CCS at 82$-$94 GHz using the Nobeyama 45-m telescope. This study aimed at investigating the physical properties of cores that show high deuterium fractions and might be close to the onset of star formation. We found that the distributions of the N$_2$H$^+$ and HC$_3$N line emissions are approximately similar to that of 850-$\mu$m dust continuum emission, whereas the CCS line emission is often undetected or is distributed in a clumpy structure surrounding the peak position of the 850-$\mu$m dust continuum emission. Occasionally (12%), we observe the CCS emission which is an early-type gas tracer toward the young stellar object, probably due to local high excitation. Evolution toward star formation does not immediately affect nonthermal velocity dispersion., Comment: 45 pages, 15 figures, ApJS, in press

Pre-stellar cores represent the initial conditions of star formation. Although these initial conditions in nearby low-mass star-forming regions have been investigated in detail, such initial conditions remain vastly unexplored for massive star-forming regions. We report the detection of a cluster of low-mass starless and pre-stellar core candidates in a massive star protocluster forming cloud, NGC6334S. With the ALMA observations at a $\sim$0.02 pc spatial resolution, we identified 17 low-mass starless core candidates that do not show any evidence of protostellar activity. These candidates present small velocity dispersions, high fractional abundances of NH$_{2}$D, high NH$_{3}$ deuterium fractionations, and are completely dark in the infrared wavelengths from 3.6 up to 70~$\mu$m. Turbulence is significantly dissipated and the gas kinematics are dominated by thermal motions toward these candidates. Nine out of the 17 cores are gravitationally bound, and therefore are identified as pre-stellar core candidates. The embedded cores of NGC6334S show a wide diversity in masses and evolutionary stages., Comment: 13 pages, 4 figures, accepted for publication in ApJ Letters

We observe 1.3~mm spectral lines at 2000~AU resolution toward four massive molecular clouds in the Central Molecular Zone of the Galaxy to investigate their star formation activities. We focus on several potential shock tracers that are usually abundant in protostellar outflows, including SiO, SO, CH$_3$OH, H$_2$CO, HC$_3$N, and HNCO. We identify 43 protostellar outflows, including 37 highly likely ones and 6 candidates. The outflows are found toward both known high-mass star forming cores and less massive, seemingly quiescent cores, while 791 out of the 834 cores identified based on the continuum do not have detected outflows. The outflow masses range from less than 1~$M_\odot$ to a few tens of $M_\odot$, with typical uncertainties of a factor of 70. We do not find evidence of disagreement between relative molecular abundances in these outflows and in nearby analogs such as the well-studied L1157 and NGC7538S outflows. The results suggest that i) protostellar accretion disks driving outflows ubiquitously exist in the CMZ environment, ii) the large fraction of candidate starless cores is expected if these clouds are at very early evolutionary phases, with a caveat on the potential incompleteness of the outflows, iii) high-mass and low-mass star formation is ongoing simultaneously in these clouds, and iv) current data do not show evidence of difference between the shock chemistry in the outflows that determines the molecular abundances in the CMZ environment and in nearby clouds., Comment: 43 pages, 25 figures, 3 tables, accepted for publication in ApJ

We have observed the very low-mass Class 0 protostar IRAS 15398-3359 at scales ranging from 50 au to 1800 au, as part of the ALMA Large Program FAUST. We uncover a linear feature, visible in H2CO, SO, and C18O line emission, which extends from the source along a direction almost perpendicular to the known active outflow. Molecular line emission from H2CO, SO, SiO, and CH3OH further reveals an arc-like structure connected to the outer end of the linear feature and separated from the protostar, IRAS 15398-3359, by 1200 au. The arc-like structure is blue-shifted with respect to the systemic velocity. A velocity gradient of 1.2 km/s over 1200 au along the linear feature seen in the H2CO emission connects the protostar and the arc-like structure kinematically. SO, SiO, and CH3OH are known to trace shocks, and we interpret the arc-like structure as a relic shock region produced by an outflow previously launched by IRAS 15398-3359. The velocity gradient along the linear structure can be explained as relic outflow motion. The origins of the newly observed arc-like structure and extended linear feature are discussed in relation to turbulent motions within the protostellar core and episodic accretion events during the earliest stage of protostellar evolution.

We present the results of a single-pointing survey of 207 dense cores embedded in Planck Galactic Cold Clumps distributed in five different environments ($\lambda$ Orionis, Orion A, B, Galactic plane, and high latitudes) to identify dense cores on the verge of star formation for the study of the initial conditions of star formation. We observed these cores in eight molecular lines at 76-94 GHz using the Nobeyama 45-m telescope. We find that early-type molecules (e.g., CCS) have low detection rates and that late-type molecules (e.g., N$_2$H$^+$, c-C$_3$H$_2$) and deuterated molecules (e.g., N$_2$D$^+$, DNC) have high detection rates, suggesting that most of the cores are chemically evolved. The deuterium fraction (D/H) is found to decrease with increasing distance, indicating that it suffers from differential beam dilution between the D/H pair of lines for distant cores ($>$1 kpc). For $\lambda$ Orionis, Orion A, and B located at similar distances, D/H is not significantly different, suggesting that there is no systematic difference in the observed chemical properties among these three regions. We identify at least eight high D/H cores in the Orion region and two at high latitudes, which are most likely to be close to the onset of star formation. There is no clear evidence of the evolutionary change in turbulence during the starless phase, suggesting that the dissipation of turbulence is not a major mechanism for the beginning of star formation as judged from observations with a beam size of 0.04 pc., Comment: 55 pages, 12 Figures, 8 Tables, accepted for publication in ApJS

We mapped two molecular cloud cores in the Orion A cloud with the ALMA ACA 7-m Array and with the Nobeyama 45-m radio telescope. These cores have bright N$_2$D$^+$ emission in single-pointing observations with the Nobeyama 45-m radio telescope, have relatively high deuterium fraction, and are thought to be close to the onset of star formation. One is a star-forming core, and the other is starless. These cores are located along filaments observed in N$_2$H$^+$, and show narrow linewidths of 0.41 km s$^{-1}$ and 0.45 km s$^{-1}$ in N$_2$D$^+$, respectively, with the Nobeyama 45-m telescope. Both cores were detected with the ALMA ACA 7m Array in the continuum and molecular lines at Band 6. The starless core G211 shows clumpy structure with several sub-cores, which in turn show chemical differences. Also, the sub-cores in G211 have internal motions that are almost purely thermal. The starless sub-core G211D, in particular, shows a hint of the inverse P Cygni profile, suggesting infall motion. The star-forming core G210 shows an interesting spatial feature of two N$_2$D$^+$ peaks of similar intensity and radial velocity located symmetrically with respect to the single dust continuum peak. One interpretation is that the two N$_2$D$^+$ peaks represent an edge-on pseudo-disk. The CO outflow lobes, however, are not directed perpendicular to the line connecting both N$_2$D$^+$ peaks., Comment: 31 pages, 17 figures, 9 tables, to appear in ApJ (accepted April 23, 2020)

We present an observational study of the sulfur (S)-bearing species towards Orion KL at 1.3 mm by combining ALMA and IRAM-30\,m single-dish data. At a linear resolution of $\sim$800 au and a velocity resolution of 1 $\mathrm{km\, s^{-1}\, }$, we have identified 79 molecular lines from 6 S-bearing species. In these S-bearing species, we found a clear dichotomy between carbon-sulfur compounds and carbon-free S-bearing species in various characteristics, e.g., line profiles, spatial morphology, and molecular abundances with respect to $\rm H_2$. Lines from the carbon-sulfur compounds (i.e., OCS, $^{13}$CS, H$_2$CS) exhibit spatial distributions concentrated around the continuum peaks and extended to the south ridge. The full width at half maximum (FWHM) linewidth of these molecular lines is in the range of 2 $\sim$ 11 $\mathrm{km\, s^{-1}\, }$. The molecular abundances of OCS and H$_2$CS decrease slightly from the cold ($\sim$68 K) to the hot ($\sim$176 K) regions. In contrast, lines from the carbon-free S-bearing species (i.e., SO$_2$, $^{34}$SO, H$_2$S) are spatially more extended to the northeast of mm4, exhibiting broader FWHM linewidths (15 $\sim$ 26 $\mathrm{km\, s^{-1}\, }$). The molecular abundances of carbon-free S-bearing species increase by over an order of magnitude as the temperature increase from 50 K to 100 K. In particular, $\mathrm{^{34}SO/^{34}SO_2}$ and $\mathrm{OCS/SO_2}$ are enhanced from the warmer regions ($>$100 K) to the colder regions ($\sim$50 K). Such enhancements are consistent with the transformation of SO$_2$ at warmer regions and the influence of shocks.

Towards the pre-stellar core L1544, the methanol (CH$_3$OH) emission forms an asymmetric ring around the core centre, where CH$_3$OH is mostly in solid form, with a clear peak 4000~au to the north-east of the dust continuum peak. As part of the NOEMA Large Project SOLIS (Seeds of Life in Space), the CH$_3$OH peak has been spatially resolved to study its kinematics and physical structure and to investigate the cause behind the local enhancement. We find that methanol emission is distributed in a ridge parallel to the main axis of the dense core. The centroid velocity increases by about 0.2~km~s$^{-1}$ and the velocity dispersion increases from subsonic to transonic towards the central zone of the core, where the velocity field also shows complex structure. This could be indication of gentle accretion of material onto the core or interaction of two filaments, producing a slow shock. We measure the rotational temperature and show that methanol is in local thermodynamic equilibrium (LTE) only close to the dust peak, where it is significantly depleted. The CH$_3$OH column density, $N_{tot}({\rm CH_3OH})$, profile has been derived with non-LTE radiative transfer modelling and compared with chemical models of a static core. The measured $N_{tot}({\rm CH_3OH})$ profile is consistent with model predictions, but the total column densities are one order of magnitude lower than those predicted by models, suggesting that the efficiency of reactive desorption or atomic hydrogen tunnelling adopted in the model may be overestimated; or that an evolutionary model is needed to better reproduce methanol abundance., Comment: 19 pages, 18 figures. Accepted by ApJ

Filamentary structures are ubiquitous in high-mass star-forming molecular clouds. Their relation with high-mass star formation is still to be understood. Here we report interferometric observations toward 8 filamentary high-mass star-forming clouds. A total of 50 dense cores are identified in these clouds, most of which present signatures of high-mass star formation. Five of them are not associated with any star formation indicators, hence are prestellar core candidates. Evolutionary phases of these cores and their linewidths, temperatures, NH$_3$ abundances, and virial parameters are found to be correlated. In a sub-sample of 4 morphologically well-defined filaments, we find that their fragmentation can not be solely explained by thermal or turbulence pressure support. We also investigate distributions of gas temperatures and non-thermal motions along the filaments, and find a spatial correlation between non-thermal linewidths and star formation activities. We find evidence of gas flows along these filaments, and derive an accretion rate along filaments of $\sim$10$^{-4}$ M$_\odot$ yr$^{-1}$. These results suggest a strong relationship between massive filaments and high-mass star formation, through i) filamentary fragmentation in very early evolutionary phases to form dense cores, ii) accretion flows along filaments that are important for the growth of dense cores and protostars, and iii) enhancement of non-thermal motion in the filaments by the feedback or accretion during star formation., Comment: 31 pages, 11 figures, 8 tables, ApJ Accepted. 3-D interactive figures to visualize outflows (see Figure 7) will be published online, and are available at https://xinglunju.github.io/outflows.html