Source author record

D. Harsono

D. Harsono appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2026arXiv

Characterizing the physical and chemical properties of the Class I protostellar system Oph-IRS 44. Binarity, infalling streamers, and accretion shocks

(Abridged) In the low-mass star formation process, theoretical models predict that material from the infalling envelope could be shocked as it encounters the outer regions of the disk. Nevertheless, only a few protostars show evidence of these shocks at the disk-envelope interface, and the main formation path of shocked-related species is still unclear. We present new ALMA observations of IRS 44, a Class I source that has previously been associated with accretion shocks, taken at high angular resolution (0.1"). We target multiple molecular transitions of CO, H2CO, and simple sulfur-bearing species. In continuum emission, the binary nature of IRS 44 is observed for the first time at sub-millimeter wavelengths. Infalling signatures are seen for the CO line and the emission peaks at the edges of the continuum emission around IRS 44 B, the same region where bright SO and SO2 emission is seen. Weak CS and H2CO emission is observed, while OCS, H2S, and H2CS transitions are not detected. IRS 44 B seems to be more embedded than IRS 44 A, indicating a non-coeval formation scenario or the rejuvenation of source B due to late infall. CO emission is tracing the outflow component at large scales, infalling envelope material at intermediate scales, and two infalling streamer candidates are identified at disk scales. Infalling streamers might produce accretion shocks when they encounter the outer regions of the infalling-rotating envelope. These shocks heat the dust and release S-bearing species as well as promoting a lukewarm chemistry in the gas phase. With the majority of carbon locked in CO, there is little free C available to form CS and H2CS in the gas, leaving an oxygen-rich environment. The high column densities of SO and SO2 might be a consequence of two processes: direct thermal desorption from dust grains and gas-phase formation due to the availability of O and S.

preprint2025arXiv

The factors that influence protostellar multiplicity II. Gas temperature and mass in Perseus with APEX

Protostellar multiplicity is a common outcome of the star formation process. To fully understand the formation and evolution of these systems, the physical parameters of the molecular gas together with the dust must be systematically characterized. Using observations of molecular gas tracers, we characterize the physical properties of cloud cores in the Perseus molecular cloud (average distance of 295 pc) at envelope scales (5000-8000 AU). We used Atacama Pathfinder EXperiment (APEX) and Nobeyama 45m Radio Observatory (NRO) observations of DCO$^+$, H$_2$CO and c-C$_3$H$_2$ in several transitions to derive the physical parameters of the gas toward 31 protostellar systems in Perseus. Gas kinetic temperature was obtained from DCO$^+$, H$_2$CO and c-C$_3$H$_2$ line ratios. Column densities and gas masses were then calculated for each species and transition. Gas kinetic temperature and gas masses were compared with bolometric luminosity, envelope dust mass, and multiplicity to search for statistically significant correlations. Gas kinetic temperature derived from DCO$^+$, H$_2$CO and c-C$_3$H$_2$ line ratios have average values of 14 K, 26 and 16 K, respectively, with a range of 10-26 K for DCO$^+$ and c-C$_3$H$_2$. The gas kinetic temperature obtained from H$_2$CO line ratios have a range of 13-82 K. Column densities of all three molecular species are on the order of 10$^{11}$ to 10$^{14}$ cm$^{-2}$, resulting in gas masses of 10$^{-11}$ to 10$^{-9}$ M$_{\odot}$. Statistical analysis of the physical parameters finds: i) similar envelope gas and dust masses for single and binary protostellar systems; ii) multiple (>2 components) protostellar systems tend to have slightly higher gas and dust masses than binaries and single protostars; iii) a continuous distribution of gas and dust masses is observed regardless of separation between components in protostellar systems.

preprint2023arXiv

Constraints on the non-thermal desorption of methanol in the cold core LDN 429-C

Cold cores are an early step of star formation, characterized by densities > 10$^4$ cm$^{-3}$, low temperatures (< 15 K), and very low external UV radiation. We investigate the physico-chemical processes at play to tracing the origin of molecules that are predominantly formed via reactions on dust grain surfaces. We observed the cold core LDN 429-C with the NOEMA interferometer and the IRAM 30m single dish telescope in order to obtain the gas-phase abundances of key species, including CO and CH$_3$OH. Comparing the observed gas phase of methanol to its solid phase previously observed with Spitzer allows us to put quantitative constraints on the efficiency of the non-thermal desorption of this species. With physical parameters determined from available Herschel data, we computed abundance maps of 11 detected molecules with a non-local thermal equilibrium radiative transfer model. These observations allowed us to probe the molecular abundances as a function of density and visual extinction, with the variation in temperature being restrained between 12 and 18 K. We then compared the observed abundances to the predictions of the Nautilus astrochemical model. We find that all molecules have lower abundances at high densities and visual extinctions with respect to lower density regions, except for methanol. Comparing these observations with a grid of chemical models based on the local physical conditions, we were able to reproduce these observations, allowing only the parameter time to vary. Comparing the observed gas-phase abundance of methanol with previous measurements of the methanol ice, we estimate a non-thermal desorption efficiency between 0.002% and 0.09%, increasing with density. The apparent increase in the desorption efficiency cannot be reproduced by our model unless the yield of cosmic-ray sputtering is altered due to the ice composition varying as a function of density.

preprint2022arXiv

Observational constraints on disc sizes in protoplanetary discs in multiple systems in the Taurus region. II. Gas disc sizes

The formation of multiple stellar systems is a natural by-product of the star-formation process, and its impact on the properties of protoplanetary discs and on the formation of planets is still to be fully understood. To date, no detailed uniform study of the gas emission from a sample of protoplanetary discs around multiple stellar systems has been performed. Here we analyse new ALMA observations at a $\sim$21 au resolution of the molecular CO gas emission targeting discs in eight multiple stellar systems in the Taurus star-forming regions. $^{12}$CO gas emission is detected around all primaries and in seven companions. With these data, we estimate the inclination and the position angle for all primary discs and for five secondary or tertiary discs, and measure the gas disc radii of these objects with a cumulative flux technique on the spatially resolved zeroth moment images. When considering the radius including 95\% of the flux as a metric, the estimated gas disc size in multiple stellar systems is found to be on average $\sim 4.2$ times larger than the dust disc size. This ratio is higher than what was recently found in a population of more isolated and single systems. On the contrary, when considering the radius including 68\% of the flux, no difference between multiple and single discs is found in the distribution of ratios. This discrepancy is due to the sharp truncation of the outer dusty disc observed in multiple stellar systems. The measured gas disc sizes are consistent with tidal truncation models in multiple stellar systems assuming eccentricities of $\sim0.15$-$0.5$, as expected in typical binary systems.

preprint2022arXiv

Tracing pebble drift and trapping using radial carbon depletion profiles in protoplanetary disks

The composition of planets may be largely determined by the chemical processing and accretion of icy pebbles in protoplanetary disks. Recent observations of protoplanetary disks hint at wide-spread depletion of gaseous carbon. The missing volatile carbon is likely frozen in CO and/or CO$_2$ ice on grains and locked into the disk through pebble trapping in pressure bumps or planetesimals. We present the results of the first successful ACA (Atacama Compact Array) [C I] $J$ = 1-0 mini-survey of seven protoplanetary disks. Using tailored azimuthally symmetric DALI (Dust And LInes) thermo-chemical disk models, supported by the [C I] $J$ = 1-0 and resolved CO isotopologue data, we determine the system-averaged elemental volatile carbon abundance in the outer disk of three sources. Six out of seven sources are detected in [C I] $J$ = 1-0 with ACA, four of which show a distinct disk component. Based on the modeling we find severe cold gaseous carbon depletion in the outer disk of DL Tau and moderate depletion in the outer disks of DR Tau and DO Tau. Combining the outer and inner disk carbon abundances, we demonstrate definitive evidence for radial drift in the disk of DL Tau, where the existence of multiple dust rings points to either short lived or leaky dust traps. We find dust locking in the compact and smooth disks of DO Tau and DR Tau, hinting at unresolved dust substructure. Comparing our results with stars of different ages and luminosities, we identify an observational evolutionary trend in gaseous carbon depletion that is consistent with dynamical models of CO depletion processes. Transport efficiency of solids in protoplanetary disks can significantly differ from what we expect based on the current resolved substructure in the continuum observations. This has important implications for our understanding of the impact of radial drift and pebble accretion on planetary compositions.

preprint2020arXiv

APEX-SEPIA660 Early Science: Gas at densities above $10^7$ cm$^{-3}$ towards OMC-1

Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1 region located in the vicinity of the Orion Nebula Cluster (ONC). Methods. We mapped the molecular emission at 652 GHz in OMC-1 as part of the APEX-SEPIA660 Early Science. Results. We detect bright and extended N$_2$H$^+$ (J=7-6) line emission along the entire OMC-1 region. Comparisons with previous ALMA data of the (J=1-0) transition and radiative transfer models indicate that the line intensities observed in this N$_2$H$^+$ (7-6) line are produced by large mass reservoirs of gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Conclusions. The first detection of this N$_2$H$^+$ (7-6) line at parsec-scales demonstrates the extreme density conditions of the star-forming gas in young massive clusters such as the ONC. Our results highlight the unique combination of sensitivity and mapping capabilities of the new SEPIA660 receiver for the study of the ISM properties at high frequencies.

preprint2020arXiv

Spirals inside the millimeter cavity of transition disk SR 21

Hydrodynamical simulations of planet-disk interactions suggest that planets may be responsible for a number of the sub-structures frequently observed in disks in both scattered light and dust thermal emission. Despite the ubiquity of these features, direct evidence of planets embedded in disks and of the specific interaction features like spiral arms within planetary gaps still remain rare. In this study we discuss recent observational results in the context of hydrodynamical simulations in order to infer the properties of a putative embedded planet in the cavity of a transition disk. We imaged the transition disk SR 21 in H-band in scattered light with SPHERE/IRDIS and in thermal dust emission with ALMA band 3 (3mm) observations at a spatial resolution of 0.1". We combine these datasets with existing band 9 (430um) and band 7 (870um) ALMA continuum data. The Band 3 continuum data reveals a large cavity and a bright ring peaking at 53 au strongly suggestive of dust trapping.The ring shows a pronounced azimuthal asymmetry, with a bright region in the north-west that we interpret as a dust over-density. A similarly-asymmetric ring is revealed at the same location in polarized scattered light, in addition to a set of bright spirals inside the mm cavity and a fainter spiral bridging the gap to the outer ring. These features are consistent with a number of previous hydrodynamical models of planet-disk interactions, and suggest the presence of a ~1 MJup planet at 44 au and PA=11°. This makes SR21 the first disk showing spiral arms inside the mm cavity, as well as one for which the location of a putative planet can be precisely inferred. With the location of a possible planet being well-constrained by observations, it is an ideal candidate for follow-up observations to search for direct evidence of a planetary companion still embedded in its disk.

preprint2016arXiv

A multi-wavelength characterization of proto-brown dwarf candidates in Serpens

We present results from a deep sub-millimeter survey in the Serpens Main and Serpens/G3-G6 clusters, conducted with the Submillimetre Common-User Bolometer Array (SCUBA-2) at the James Clerk Maxwell Telescope. We have combined Herschel PACS far-infrared photometry, sub-millimeter continuum and molecular gas line observations, with the aim to conduct a detailed multi-wavelength characterization of `proto-brown dwarf' candidates in Serpens. We have performed continuum and line radiative transfer modeling, and have considered various classification schemes to understand the structure and the evolutionary stage of the system. We have identified four proto-brown dwarf candidates, of which the lowest luminosity source has an Lbol ~0.05 Lsun. Two of these candidates show characteristics consistent with Stage 0/I systems, while the other two are Stage I-T/Class Flat systems with tenuous envelopes. Our work has also revealed a ~20% fraction of mis-identified Class 0/I/Flat sources that show characteristics consistent with Class II edge-on disk systems. We have set constraints on the mass of the central object using the measured bolometric luminosities and numerical simulations of stellar evolution. Considering the available gas+dust mass reservoir and the current mass of the central source, three of these candidates are likely to evolve into brown dwarfs.

preprint2016arXiv

Constraining the physical structure of the inner few 100 AU scales of deeply embedded low-mass protostars

(Abridged) The physical structure of deeply-embedded low-mass protostars (Class 0) on scales of less than 300 AU is still poorly constrained. Determining this is crucial for understanding the physical and chemical evolution from cores to disks. In this study two models of the emission, a Gaussian disk intensity distribution and a parametrized power-law disk model, are fitted to sub-arcsecond resolution interferometric continuum observations of five Class 0 sources, including one source with a confirmed Keplerian disk. For reference, a spherically symmetric single power-law envelope is fitted to the larger scale ($\sim$1000 AU) emission and investigated further for one of the sources on smaller scales. A thin disk model can approximate the emission and physical structure in the inner few 100 AU scales of the studied deeply-embedded low-mass protostars and paves the way for analysis of a larger sample with ALMA. While the disk radii agree with previous estimates the masses are different for some of the sources studied. Assuming a typical temperature distribution, the fractional amount of mass in the disk above 100 K varies in between 7% to 30%. Kinematic data are needed to determine the presence of any Keplerian disk. Using previous observations of p-H$_2^{18}$O, we estimate the relative gas phase water abundances roughly an order of magnitude higher than previously inferred when both warm and cold H$_2$ was used as reference. A spherically symmetric single power-law envelope model fails to simultaneously reproduce both the small and large scale emission.

preprint2016arXiv

Water around IRAS15398-3359 observed with ALMA

How protostars accrete mass is one of the fundamental problems of star formation. High column densities and complex kinematical structures make direct observations challenging and they only provide a snapshot. Chemical tracers provide an interesting alternative to characterise the infall histories of protostars. Previous observations of H13CO+ towards the low-mass protostar IRAS15398-3359 showed a depression in the abundance. This is a sign of destruction of HCO+ by an enhanced presence of gaseous water in an extended region, possibly related to a recent burst in the accretion. Direct observations of water vapour can determine the exact extent of the emission and confirm the hypothesis that HCO+ is indeed a good tracer of the water snow-line. IRAS15398 was observed using ALMA at 0.5" resolution. Maps of HDO(101-000) and H218O(414-321) were taken simultaneously with observations of the CS(8-7) and N2H+(5-4) lines and continuum at 0.65 and 0.75 mm. The maps were interpreted using dust radiative transfer calculations of the protostellar infalling envelope with an outflow cavity. HDO is clearly detected and extended over the scales of the H13CO+ depression, although it is displaced by ~500 AU in the direction of the outflow. H218O is tentatively detected towards the red-shifted outflow lobe, but otherwise it is absent from the mapped region, which suggests that temperatures are low. Based on the temperature structure obtained from dust radiative transfer models, we conclude that the water was most likely released from the grains in an extended hour-glass configuration during a recent accretion burst. HDO is only detected in the region closest to the protostar, at distances of up to 500 AU. These signatures can only be explained if the luminosity has recently been increased by orders of magnitudes. Additionally, the densities in the outflow cones must be sufficiently low.

preprint2015arXiv

Outflow forces in intermediate mass star formation

Intermediate mass protostarsprovide a bridge between theories of low- and high-mass star formation. Emerging molecular outflows can be used to determine the influence of fragmentation and multiplicity on protostellar evolution through the correlation of outflow forces of intermediate mass protostars with the luminosity. The aim of this paper is to derive outflow forces from outflows of six intermediate mass protostellar regions and validate the apparent correlation between total luminosity and outflow force seen in earlier work, as well as remove uncertainties caused by different methodology. By comparing CO 6--5 observations obtained with APEX with non-LTE radiative transfer model predictions, optical depths, temperatures, densities of the gas of the molecular outflows are derived. Outflow forces, dynamical timescales and kinetic luminosities are subsequently calculated. Outflow parameters, including the forces, were derived for all sources. Temperatures in excess of 50 K were found for all flows, in line with recent low-mass results. However, comparison with other studies could not corroborate conclusions from earlier work on intermediate mass protostars which hypothesized that fragmentation enhances outflow forces in clustered intermediate mass star formation. Any enhancement in comparison with the classical relation between outflow force and luminosity can be attributed the use of a higher excitation line and improvement in methods; They are in line with results from low-mass protostars using similar techniques. The role of fragmentation on outflows is an important ingredient to understand clustered star formation and the link between low and high-mass star formation. However, detailed information on spatial scales of a few 100 AU, covering all individual members is needed to make the necessary progress.

preprint2014arXiv

The deuterium fractionation of water on solar-system scales in deeply-embedded low-mass protostars

(Abridged) The water deuterium fractionation (HDO/H$_2$O abundance ratio) has traditionally been used to infer the amount of water brought to Earth by comets. Measuring this ratio in deeply-embedded low-mass protostars makes it possible to probe the critical stage when water is transported from clouds to disks in which icy bodies are formed. We present sub-arcsecond resolution observations of HDO in combination with H$_2^{18}$O from the PdBI toward the three low-mass protostars NGC 1333-IRAS 2A, IRAS 4A-NW, and IRAS 4B. The resulting HDO/H$_2$O ratio is $7.4\pm2.1\times10^{-4}$ for IRAS 2A, $19.1\pm5.4\times10^{-4}$ for IRAS 4A-NW, and $5.9\pm1.7\times10^{-4}$ for IRAS 4B. Derived ratios agree with radiative transfer models within a factor of 2-4 depending on the source. Our HDO/H$_2$O ratios for the inner regions (where $T>100$ K) of four young protostars are only a factor of 2 higher than those found for pristine, solar system comets. These small differences suggest that little processing of water occurs between the deeply embedded stage and the formation of planetesimals and comets.

preprint2013arXiv

ALMA imaging of the CO snowline of the HD 163296 disk with DCO+

The high spatial and line sensitivity of ALMA opens the possibility of resolving emission from molecules in circumstellar disks. With an understanding of physical conditions under which molecules have high abundance, they can be used as direct tracers of distinct physical regions. In particular, DCO+ is expected to have an enhanced abundance within a few Kelvin of the CO freezeout temperature of 19 K, making it a useful probe of the cold disk midplane. We compare ALMA line observations of HD 163296 to a grid of models. We vary the upper- and lower-limit temperatures of the region in which DCO+ is present as well as the abundance of DCO+ in order to fit channel maps of the DCO+ J=5-4 line. To determine the abundance enhancement compared to the general interstellar medium, we carry out similar fitting to HCO+ J=4-3 and H13CO+ J=4-3 observations. ALMA images show centrally peaked extended emission from HCO+ and H13CO+. DCO+ emission lies in a resolved ring from ~110 to 160 AU. The outer radius approximately corresponds to the size of the CO snowline as measured by previous lower resolution observations of CO lines in this disk. The ALMA DCO+ data now resolve and image the CO snowline directly. In the best fitting models, HCO+ exists in a region extending from the 19 K isotherm to the photodissociation layer with an abundance of 3x10^-10 relative to H2. DCO+ exists within the 19-21 K region of the disk with an abundance ratio [DCO+] / [HCO+] = 0.3. This represents a factor of 10^4 enhancement of the DCO+ abundance within this narrow region of the HD 163296 disk. Such a high enhancement has only previously been seen in prestellar cores. The inferred abundances provide a lower limit to the ionization fraction in the midplane of the cold outer disk (approximately greater than 4x10^-10), and suggest the utility of DCO+ as a tracer of its parent molecule H2D+. Abridged

preprint2012arXiv

Water in star-forming regions with Herschel (WISH): II. Evolution of 557 GHz 110-101 emission in low-mass protostars

(Abridged) Water is a key tracer of dynamics and chemistry in low-mass protostars, but spectrally resolved observations have so far been limited in sensitivity and angular resolution. In this first systematic survey of spectrally resolved water emission in low-mass protostellar objects, H2O was observed in the ground-state transition at 557 GHz with HIFI on Herschel in 29 embedded Class 0 and I protostars. Complementary far-IR and sub-mm continuum data (including PACS data from our program) are used to constrain the spectral energy distribution of each source. H2O intensities are compared to inferred envelope and outflow properties and CO 3-2 emission. H2O emission is detected in all objects except one. The line profiles are complex and consist of several kinematic components. The profiles are typically dominated by a broad Gaussian emission feature, indicating that the bulk of the water emission arises in outflows, not the quiescent envelope. Several sources show multiple shock components in either emission or absorption, thus constraining the internal geometry of the system. Furthermore, the components include inverse P-Cygni profiles in 7 sources (6 Class 0, 1 Class I) indicative of infalling envelopes, and regular P-Cygni profiles in 4 sources (3 Class I, 1 Class 0) indicative of expanding envelopes. "Bullets" moving at >50 km/s are seen in 4 Class 0 sources; 3 of these are new detections. In the outflow, the H2O/CO abundance ratio as a function of velocity is nearly the same for all sources, increasing from 10^-3 at <5 km/s to >10^-1 at >10 km/s. The H2O abundance in the outer envelope is low, ~10^-10. The different H2O profile components show a clear evolutionary trend: in the Class 0 sources, emission is dominated by outflow components originating inside an infalling envelope. When the infall diminishes during the Class I phase, the outflow weakens and H2O emission disappears.

preprint2011arXiv

Modelling Herschel observations of hot molecular gas emission from embedded low-mass protostars

Aims. Young stars interact vigorously with their surroundings, as evident from the highly rotationally excited CO (up to Eup=4000 K) and H2O emission (up to 600 K) detected by the Herschel Space Observatory in embedded low-mass protostars. Our aim is to construct a model that reproduces the observations quantitatively, to investigate the origin of the emission, and to use the lines as probes of the various heating mechanisms. Methods. The model consists of a spherical envelope with a bipolar outflow cavity. Three heating mechanisms are considered: passive heating by the protostellar luminosity, UV irradiation of the outflow cavity walls, and C-type shocks along the cavity walls. Line fluxes are calculated for CO and H2O and compared to Herschel data and complementary ground-based data for the protostars NGC1333 IRAS2A, HH 46 and DK Cha. The three sources are selected to span a range of evolutionary phases and physical characteristics. Results. The passively heated gas in the envelope accounts for 3-10% of the CO luminosity summed over all rotational lines up to J=40-39; it is best probed by low-J CO isotopologue lines such as C18O 2-1 and 3-2. The UV-heated gas and the C-type shocks, probed by 12CO 10-9 and higher-J lines, contribute 20-80% each. The model fits show a tentative evolutionary trend: the CO emission is dominated by shocks in the youngest source and by UV-heated gas in the oldest one. This trend is mainly driven by the lower envelope density in more evolved sources. The total H2O line luminosity in all cases is dominated by shocks (>99%). The exact percentages for both species are uncertain by at least a factor of 2 due to uncertainties in the gas temperature as function of the incident UV flux. However, on a qualitative level, both UV-heated gas and C-type shocks are needed to reproduce the emission in far-infrared rotational lines of CO and H2O.

preprint2010arXiv

Global Gravitational Instabilities in Discs with Infall

Gravitational instability plays an important role in driving gas accretion in massive protostellar discs. Particularly strong is the global gravitational instability, which arises when the disc mass is of order 0.1 of the mass of the central star and has a characteristic spatial scale much greater than the disc's vertical scale-height. In this paper we use three-dimensional numerical hydrodynamics to study the development of gravitational instabilities in a disc which is embedded in a dense, gaseous envelope. We find that global gravitational instabilities are the dominant mode of angular momentum transport in the disc with infall, in contrast to otherwise identical isolated discs. The accretion torques created by low-order, global modes of the gravitational instability in a disc subject to infall are larger by a factor of several than an isolated disc of the same mass. We show that this global gravitational instability is driven by the strong vertical shear at the interface between the disc and the envelope, and suggest that this process may be an important means of driving accretion on to young stars.