-
HCN as a probe of the inner disk in a candidate proto-brown dwarf
Authors:
B. Riaz,
W. -F. Thi,
M. N. Machida
Abstract:
The detection of Keplerian rotation is rare among Class 0 protostellar systems. We have investigated the high-density tracer HCN as a probe of the inner disk in a Class 0 proto-brown dwarf candidate. Our ALMA high angular resolution observations show the peak in the HCN (3-2) line emission arises from a compact component near the proto-brown dwarf with a small bar-like structure and a deconvolved…
▽ More
The detection of Keplerian rotation is rare among Class 0 protostellar systems. We have investigated the high-density tracer HCN as a probe of the inner disk in a Class 0 proto-brown dwarf candidate. Our ALMA high angular resolution observations show the peak in the HCN (3-2) line emission arises from a compact component near the proto-brown dwarf with a small bar-like structure and a deconvolved size of $\sim$50 au. Radiative transfer modelling indicates that this HCN feature is tracing the innermost, dense regions in the proto-brown dwarf where a small Keplerian disk is expected to be present. The limited velocity resolution of the observations, however, makes it difficult to confirm the rotational kinematics of this feature. A brightening in the HCN emission towards the core center suggests that HCN can survive in the gas phase in the inner, dense regions of the proto-brown dwarf. In contrast, modelling of the HCO$^{+}$ (3-2) line emission indicates that it originates from the outer pseudo-disk/envelope region and is centrally depleted. HCN line emission can reveal the small-scale structures and can be an efficient observational tool to study the inner disk properties in such faint compact objects where spatially resolving the disk is nearly impossible.
△ Less
Submitted 18 May, 2024;
originally announced May 2024.
-
Cloud Dissipation and Disk Wind in the Late Phase of Star Formation
Authors:
Masahiro N. Machida,
Shantanu Basu
Abstract:
We perform a long-term simulation of star and disk formation using three-dimensional non-ideal magnetohydrodynamics. The simulation starts from a prestellar cloud and proceeds through the long-term evolution of the circumstellar disk until $\sim 1.5\times10^5$ yr after protostar formation. The disk has size $\lesssim 50$ au and little substructure in the main accretion phase because of the action…
▽ More
We perform a long-term simulation of star and disk formation using three-dimensional non-ideal magnetohydrodynamics. The simulation starts from a prestellar cloud and proceeds through the long-term evolution of the circumstellar disk until $\sim 1.5\times10^5$ yr after protostar formation. The disk has size $\lesssim 50$ au and little substructure in the main accretion phase because of the action of magnetic braking and the magnetically-driven outflow to remove angular momentum. The main accretion phase ends when the outflow breaks out of the cloud, causing the envelope mass to decrease rapidly. The outflow subsequently weakens as the mass accretion rate also weakens. While the envelope-to-disk accretion continues, the disk grows gradually and develops transient spiral structures due to gravitational instability. When the envelope-to-disk accretion ends, the disk becomes stable and reaches a size $\gtrsim 300$ au. In addition, about 30% of the initial cloud mass has been ejected by the outflow. A significant finding of this work is that after the envelope dissipates, a revitalization of the wind occurs, and there is mass ejection from the disk surface that lasts until the end of the simulation. This mass ejection (or disk wind) is generated since the magnetic pressure significantly dominates both the ram pressure and thermal pressure above and below the disk at this stage. Using the angular momentum flux and mass loss rate estimated from the disk wind, the disk dissipation timescale is estimated to be $\sim10^6$ yr.
△ Less
Submitted 13 May, 2024;
originally announced May 2024.
-
Delivery of Dust Particles from Protoplanetary Disks onto Circumplanetary Disks of Giant Planets
Authors:
Natsuho Maeda,
Keiji Ohtsuki,
Ryo Suetsugu,
Yuhito Shibaike,
Takayuki Tanigawa,
Masahiro N. Machida
Abstract:
Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a…
▽ More
Principal regular satellites of gas giants are thought to be formed by the accumulation of solid materials in circumplanetary disks (CPDs). While there has been significant progress in the study of satellite formation in CPDs, details of the supply of satellite building blocks to CPDs remain unclear. We performed orbital integration of solid particles in the protoplanetary disk (PPD) approaching a planet, considering the gas drag force using the results of three-dimensional hydrodynamical simulations of a local region around the planet. We investigated planetary-mass dependence of the capture positions and capture rates of dust particles accreting onto the CPD. We also examined the degree of dust retention in accreting gas onto the CPD, which is important for determining the ratio of dust-to-gas inflow rates, a key parameter in satellite formation. We found that the degree of dust retention increases with increasing planetary mass for a given dust scale height in the PPD. In the case of a small planet ($M_{\rm p}=0.2M_{\rm Jup}$), most particles with insufficient initial altitudes in the PPD are isolated from the gas in the accreting region. On the other hand, in the case of a massive planet ($M_{\rm p}=1M_{\rm Jup}$), dust particles can be coupled to the vertically accreting gas, even when the dust scale height is about $10-30$\% of the gas scale height. The results of this study can be used for models of dust delivery and satellite formation in the CPDs of gas giants of various masses, including exoplanets.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Internal 1000 AU-scale Structures of the R CrA Cluster-forming Cloud -- I: Filamentary Structures
Authors:
Kengo Tachihara,
Naofumi Fukaya,
Kazuki Tokuda,
Yasumasa Yamasaki,
Takeru Nishioka,
Daisei Abe,
Tsuyoshi Inoue,
Naoto Harada,
Ayumu Shoshi,
Shingo Nozaki,
Asako Sato,
Mitsuki Omura,
Kakeru Fujishiro,
Misato Fukagawa,
Masahiro N. Machida,
Takahiro Kanai,
Yumiko Oasa,
Toshikazu Onishi,
Kazuya Saigo,
Yasuo Fukui
Abstract:
We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected b…
▽ More
We report on ALMA ACA observations of a high-density region of the Corona Australis cloud forming a young star cluster, and the results of resolving internal structures. In addition to embedded Class 0/I protostars in continuum, a number of complex dense filamentary structures are detected in the C18O and SO lines by the 7m array. These are sub-structures of the molecular clump that are detected by the TP array as the extended emission. We identify 101 and 37 filamentary structures with a few thousand AU widths in C18O and SO, respectively, called as feathers. The typical column density of the feathers in C18O is about 10^{22} cm^{-2}, and the volume density and line mass are ~ 10^5 cm^{-3}, and a few times M_{sun} pc^{-1}, respectively. This line mass is significantly smaller than the critical line mass expected for cold and dense gas. These structures have complex velocity fields, indicating a turbulent internal property. The number of feathers associated with Class 0/I protostars is only ~ 10, indicating that most of them do not form stars but rather being transient structures. The formation of feathers can be interpreted as a result of colliding gas flow as the morphology well reproduced by MHD simulations, supported by the the presence of HI shells in the vicinity. The colliding gas flows may accumulate gas and form filaments and feathers, and trigger the active star formation of the R CrA cluster.
△ Less
Submitted 17 April, 2024;
originally announced April 2024.
-
Multiple Outflows around a Single Protostar IRAS 15398$-$3359
Authors:
Jinshi Sai,
Hsi-Wei Yen,
Masahiro N. Machida,
Nagayoshi Ohashi,
Yusuke Aso,
Ana�lle J. Maury,
S�bastien Maret
Abstract:
We present the results of our mosaic observations of a single Class 0 protostar IRAS 15398$-$3359 with Atacama Compact Array (ACA) in the CO $J=2\mbox{-}1$ line. The new observations covering a $\sim\!2'$ square region revealed elongated redshifted and blueshifted components, which are located at distances of $\sim\!30''\mbox{-}75''$ on the northern and southern sides of the protostar, respectivel…
▽ More
We present the results of our mosaic observations of a single Class 0 protostar IRAS 15398$-$3359 with Atacama Compact Array (ACA) in the CO $J=2\mbox{-}1$ line. The new observations covering a $\sim\!2'$ square region revealed elongated redshifted and blueshifted components, which are located at distances of $\sim\!30''\mbox{-}75''$ on the northern and southern sides of the protostar, respectively, in addition to the previously observed primary and secondary outflows. These elongated components exhibit Hubble-law like velocity structures, i.e., an increase of velocity with increasing distance from the protostar, suggesting that it is the third outflow associated with the protostar. Besides, a new redshifted component is detected at radii of $\sim\!40''\mbox{-}75''$ on the northwestern side of the protostar. This redshifted component also exhibits a Hubble-law like velocity profile, which could be the counterpart of the secondary outflow mostly detected at blueshifted velocities in a previous study. The three outflows are all misaligned by $\sim\!20\mbox{-}90^\circ$, and the dynamical timescale of the primary outflow is shorter than those of the other outflows approximately by an order of magnitude. These facts hint that the outflow launch direction has significantly changed with time. The outflow direction may change if the rotational axis and the magnetic field are misaligned, or if the dense core is turbulent. We favor the second scenario as the origin of the multiple outflows in IRAS 15398$-$3359 based on a comparison between the observational results and numerical simulations.
△ Less
Submitted 16 March, 2024;
originally announced March 2024.
-
Observations of spiral and streamer on a candidate proto-brown dwarf
Authors:
B. Riaz,
D. Stamatellos,
M. Machida
Abstract:
Spirals and streamers are the hallmarks of mass accretion during the early stages of star formation. We present the first observations of a large-scale spiral and a streamer towards a very young brown dwarf candidate in its early formation stages. These observations show, for the first time, the influence of external environment that results in asymmetric mass accretion via feeding filaments onto…
▽ More
Spirals and streamers are the hallmarks of mass accretion during the early stages of star formation. We present the first observations of a large-scale spiral and a streamer towards a very young brown dwarf candidate in its early formation stages. These observations show, for the first time, the influence of external environment that results in asymmetric mass accretion via feeding filaments onto a candidate proto-brown dwarf in the making. The impact of the streamer has produced emission in warm carbon-chain species close to the candidate proto-brown dwarf. Two contrasting scenarios, a pseudo-disk twisted by core rotation and the collision of dense cores, can both explain these structures. The former argues for the presence of a strong magnetic field in brown dwarf formation while the latter suggests that a minimal magnetic field allows large-scale spirals and clumps to form far from the candidate proto-brown dwarf.
△ Less
Submitted 12 March, 2024;
originally announced March 2024.
-
Discovery of Asymmetric Spike-like Structures of the 10 au Disk around the Very Low-luminosity Protostar Embedded in the Taurus Dense Core MC 27/L1521F with ALMA
Authors:
Kazuki Tokuda,
Naoto Harada,
Mitsuki Omura,
Tomoaki Matsumoto,
Toshikazu Onishi,
Kazuya Saigo,
Ayumu Shoshi,
Shingo Nozaki,
Kengo Tachihara,
Naofumi Fukaya,
Yasuo Fukui,
Shu-ichiro Inutsuka,
Masahiro N. Machida
Abstract:
Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations have revealed an increasing number of compact protostellar disks with radii of less than a few tens of astronomical units and that young Class 0/I objects have an intrinsic size diversity. To deepen our understanding of the origin of such tiny disks, we performed the highest-resolution configuration observations with ALMA at a…
▽ More
Recent Atacama Large Millimeter/submillimeter Array (ALMA) observations have revealed an increasing number of compact protostellar disks with radii of less than a few tens of astronomical units and that young Class 0/I objects have an intrinsic size diversity. To deepen our understanding of the origin of such tiny disks, we performed the highest-resolution configuration observations with ALMA at a beam size of $\sim$0$''$03 (4 au) on the very low-luminosity Class 0 protostar embedded in the Taurus dense core MC 27/L1521F. The 1.3 mm continuum measurement successfully resolved a tiny, faint ($\sim$1 mJy) disk with a major axis length of $\sim$10 au, one of the smallest examples in the ALMA protostellar studies. In addition, we detected spike-like components in the northeastern direction at the disk edge. Gravitational instability or other fragmentation mechanisms cannot explain the structures, given the central stellar mass of $\sim$0.2 $M_{\odot}$ and the disk mass of $\gtrsim$10$^{-4}$ $M_{\odot}$. Instead, we propose that these small spike structures were formed by a recent dynamic magnetic flux transport event due to interchange instability that would be favorable to occur if the parental core has a strong magnetic field. The presence of complex arc-like structures on a larger ($\sim$2000 au) scale in the same direction as the spike structures suggests that the event was not single. Such episodic, dynamical events may play an important role in maintaining the compact nature of the protostellar disk in the complex gas envelope during the main accretion phase.
△ Less
Submitted 3 April, 2024; v1 submitted 1 March, 2024;
originally announced March 2024.
-
An Extremely Young Protostellar Core, MMS 1/ OMC-3: Episodic Mass Ejection History Traced by the Micro SiO Jet
Authors:
Satoko Takahashi,
Masahiro N. Machida,
Mitsuki Omura,
Doug Johnstone,
Kazuya Saigo,
Naoto Harada,
Kohji Tomisaka,
Paul T. P. Ho,
Luis A. Zapata,
Steve Mairs,
Gregory J. Herczeg,
Kotomi Taniguchi,
Yuhua Liu,
Asako Sato
Abstract:
We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, havin…
▽ More
We present ${\sim}0.2$ arcsec ($\sim$80 au) resolution observations of the CO (2-1) and SiO (5-4) lines made with the Atacama large millimeter/submillimeter array toward an extremely young intermediate-mass protostellar source (t$_{\rm dyn}<$1000 years), MMS 1 located in the Orion Molecular Cloud-3 region. We have successfully imaged a very compact CO molecular outflow associated with MMS 1, having deprojected lobe sizes of $\sim$18000 au (red-shifted lobe) and $\sim$35000 au (blue-shifted lobe). We have also detected an extremely compact ($\lesssim$1000 au) and collimated SiO protostellar jet within the CO outflow. The maximum deprojected jet speed is measured to be as high as 93 km s$^{-1}$. The SiO jet wiggles and displays a chain of knots. Our detection of the molecular outflow and jet is the first direct evidence that MMS 1 already hosts a protostar. The position-velocity diagram obtained from the SiO emission shows two distinct structures: (i) bow-shocks associated with the tips of the outflow, and (ii) a collimated jet, showing the jet velocities linearly increasing with the distance from the driving source. Comparisons between the observations and numerical simulations quantitatively share similarities such as multiple-mass ejection events within the jet and Hubble-like flow associated with each mass ejection event. Finally, while there is a weak flux decline seen in the 850 $μ$m light curve obtained with JCMT/SCUBA 2 toward MMS 1, no dramatic flux change events are detected. This suggests that there has not been a clear burst event within the last 8 years.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Outflows Driven from a Magnetic Pseudodisk
Authors:
Shantanu Basu,
Mahmoud Sharkawi,
Masahiro N. Machida
Abstract:
Outflows play a pivotal role in star formation as one of its most visible markers and a means of transporting mass, momentum, and angular momentum from the infalling gas into the surrounding molecular cloud. Their wide reach (at least thousands of au) is a contrast to typical disk sizes ($\sim 10-100$ au). We employ high-resolution three-dimensional nested-grid nonideal magnetohydrodynamic (MHD) s…
▽ More
Outflows play a pivotal role in star formation as one of its most visible markers and a means of transporting mass, momentum, and angular momentum from the infalling gas into the surrounding molecular cloud. Their wide reach (at least thousands of au) is a contrast to typical disk sizes ($\sim 10-100$ au). We employ high-resolution three-dimensional nested-grid nonideal magnetohydrodynamic (MHD) simulations to study outflow properties in the Class 0 phase. We find that no disk wind is driven from the extended centrifugal disk that has weak magnetic coupling. The low-velocity winds emerge instead from the infalling magnetic pseudodisk. Much of the disk actually experiences an infall of matter rather than outflowing gas. Some of the pseudodisk wind (PD-wind) moves inward to regions above the disk and either falls onto the disk or proceeds upward. The upward flow gives the impression of a disk wind above a certain height even if the gas is originally emerging from the pseudodisk. The PD-wind has the strongest flow coming from a disk interaction zone that lies just outside the disk and is an interface between the inwardly advected magnetic field of the pseudodisk and the outwardly diffusing magnetic field of the disk. The low-velocity wind exhibits the features of a flow driven by the magnetic pressure gradient force in some regions and those of a magnetocentrifugal wind in other regions. We interpret the structure and dynamics of the outflow zone in terms of the basic physics of gravity, angular momentum, magnetic fields, and nonideal MHD.
△ Less
Submitted 14 March, 2024; v1 submitted 8 January, 2024;
originally announced January 2024.
-
Revealing multiple nested molecular outflows with rotating signatures in HH270mms1-A with ALMA
Authors:
Mitsuki Omura,
Kazuki Tokuda,
Masahiro N. Machida
Abstract:
We present molecular line observations of the protostellar outflow associated with HH270mms1 in the Orion B molecular cloud with ALMA. The 12CO(J = 3 - 2) emissions show that the outflow velocity structure consists of four distinct components of low ($\gtrsim$ 10 km s-1), intermediate (~ 10 - 25 km s-1) and high ($\gtrsim$ 40 km s-1) velocities in addition to the entrained gas velocity (~ 25 - 40…
▽ More
We present molecular line observations of the protostellar outflow associated with HH270mms1 in the Orion B molecular cloud with ALMA. The 12CO(J = 3 - 2) emissions show that the outflow velocity structure consists of four distinct components of low ($\gtrsim$ 10 km s-1), intermediate (~ 10 - 25 km s-1) and high ($\gtrsim$ 40 km s-1) velocities in addition to the entrained gas velocity (~ 25 - 40 km s-1). The high- and intermediate-velocity flows have well-collimated structures surrounded by the low-velocity flow. The chain of knots is embedded in the high-velocity flow or jet, which is the evidence of episodic mass ejections induced by time-variable mass accretion. We could detect the velocity gradients perpendicular to the outflow axis in both the low- and intermediate-velocity flows. We confirmed the rotation of the envelope and disk in the 13CO and C17O emission and found that their velocity gradients are the same as those of the outflow. Thus, we concluded that the velocity gradients in the low- and intermediate-velocity flows are due to the outflow rotation. Using observational outflow properties, we estimated the outflow launching radii to be 67.1 - 77.1 au for the low-velocity flow and 13.3 - 20.8 au for the intermediate-velocity flow. Although we could not detect the rotation in the jets due to the limited spatial resolution, we estimated the jet launching radii to be (2.36 - 3.14) x 10^-2 au using the observed velocity of each knots. Thus, the jet is driven from the inner disk region. We could identify the launching radii of distinct velocity components within a single outflow with all the prototypical characteristics expected from recent theoretical works.
△ Less
Submitted 5 January, 2024;
originally announced January 2024.
-
Dust Polarization of Prestellar and Protostellar Sources in OMC-3
Authors:
Yuhua Liu,
Satoko Takahashi,
Masahiro Machida,
Kohji Tomisaka,
Josep Miquel Girart,
Paul T. P. Ho,
Kouichiro Nakanishi,
Asako Sato
Abstract:
We present the Atacama Large Millimeter/submillimeter Array (ALMA) observations of linearly polarized 1.1 mm continuum emission at $\sim$0.14" (55 au) resolution and CO ($J$=2$-$1) emission at $\sim$1.5" (590 au) resolution towards one prestellar (MMS 4), four Class 0 (MMS$\,$1, MMS$\,$3, MMS$\,$5, and MMS$\,$6), one Class I (MMS$\,$7), and one flat-spectrum (MMS$\,$2) sources in the Orion Molecul…
▽ More
We present the Atacama Large Millimeter/submillimeter Array (ALMA) observations of linearly polarized 1.1 mm continuum emission at $\sim$0.14" (55 au) resolution and CO ($J$=2$-$1) emission at $\sim$1.5" (590 au) resolution towards one prestellar (MMS 4), four Class 0 (MMS$\,$1, MMS$\,$3, MMS$\,$5, and MMS$\,$6), one Class I (MMS$\,$7), and one flat-spectrum (MMS$\,$2) sources in the Orion Molecular Cloud$\,$3 region. The dust disk-like structures and clear CO outflows are detected towards all sources except for MMS$\,$4. The diameters of these disk-like structures, ranging from 16 au to 97 au, are estimated based on the deconvolved full width half maximum (FWHM) values obtained from the multi-Gaussian fitting. Polarized emissions are detected towards MMS$\,$2, MMS$\,$5, MMS$\,$6, and MMS$\,$7, while no polarized emission is detected towards MMS$\,$1, MMS$\,$3, and MMS$\,$4. MMS$\,$2, MMS$\,$5, and MMS$\,$7 show organized polarization vectors aligned with the minor axes of the disk-like structures, with mean polarization fractions ranging from 0.6$\%$ to 1.2$\%$. The strongest millimeter source, MMS$\,$6, exhibits complex polarization orientations and a remarkably high polarization fraction of $\sim$10$\%$ around the Stokes $I$ peak, and 15$-$20$\%$ on the arm-like structure, as reported by Takahashi et al. (2019). The origins of the polarized emission, such as self-scattering and dust alignment due to the magnetic field or radiative torque, are discussed for individual sources. Some disk-like sources exhibit a polarized intensity peak shift towards the nearside of the disk, which supports that the polarized emission originates from self-scattering.
△ Less
Submitted 20 December, 2023;
originally announced December 2023.
-
Ring Gap Structure around Class I Protostar WL 17
Authors:
Ayumu Shoshi,
Naoto Harada,
Kazuki Tokuda,
Yoshihiro Kawasaki,
Hayao Yamasaki,
Asako Sato,
Mitsuki Omura,
Masayuki Yamaguchi,
Kengo Tachihara,
Masahiro N. Machida
Abstract:
WL 17 is a Class I object and was considered to have a ring-hole structure. We analyzed the structure around WL 17 to investigate the detailed properties of WL 17. We used ALMA archival data, which have a higher angular resolution than previous observations. We investigated the WL 17 system with the 1.3 mm dust continuum and 12CO and C18O (J = 2-1) line emissions. The dust continuum emission showe…
▽ More
WL 17 is a Class I object and was considered to have a ring-hole structure. We analyzed the structure around WL 17 to investigate the detailed properties of WL 17. We used ALMA archival data, which have a higher angular resolution than previous observations. We investigated the WL 17 system with the 1.3 mm dust continuum and 12CO and C18O (J = 2-1) line emissions. The dust continuum emission showed a clear ring structure with inner and outer edges of ~11 and ~21 au, respectively. In addition, we detected an inner disk of < 5 au radius enclosing the central star within the ring, the first observation of this structure. Thus, WL 17 has a ring-gap structure, not a ring-hole structure. We did not detect any marked emission in either the gap or inner disk, indicating that there is no sign of a planet, circumplanetary disk, or binary companion. We identified the base of both blue-shifted and red-shifted outflows based on the 12CO emission, which is clearly associated with the disk around WL 17. The outflow mass ejection rate is ~3.6x10^-7 Msun yr-1 and the dynamical timescale is as short as ~ 10^4 yr. The C18O emission showed that an inhomogeneous infalling envelope, which can induce episodic mass accretion, is distributed in the region within ~1000 au from the central protostar. With these new findings, we can constrain the planet formation and dust growth scenarios in the accretion phase of star formation.
△ Less
Submitted 5 December, 2023;
originally announced December 2023.
-
Secondary outflow driven by the protostar Ser-emb 15 in Serpens
Authors:
Asako Sato,
Kazuki Tokuda,
Masahiro N. Machida,
Kengo Tachihara,
Naoto Harada,
Hayao Yamasaki,
Shingo Hirano,
Toshikazu Onishi,
Yuko Matsushita
Abstract:
We present the detection of a secondary outflow associated with a Class I source, Ser-emb 15, in the Serpens Molecular Cloud. We reveal two pairs of molecular outflows consisting of three lobes, namely primary and secondary outflows, using ALMA 12CO and SiO line observations at a resolution of 318 au. The secondary outflow is elongated approximately perpendicular to the axis of the primary outflow…
▽ More
We present the detection of a secondary outflow associated with a Class I source, Ser-emb 15, in the Serpens Molecular Cloud. We reveal two pairs of molecular outflows consisting of three lobes, namely primary and secondary outflows, using ALMA 12CO and SiO line observations at a resolution of 318 au. The secondary outflow is elongated approximately perpendicular to the axis of the primary outflow in the plane of the sky. We also identify two compact structures, Sources A and B, within an extended structure associated with Ser-emb 15 in the 1.3 mm continuum emission at a resolution of 40 au. The projected sizes of Sources A and B are 137 au and 60 au, respectively. Assuming a dust temperature of 20 K, we estimate the dust mass to be 0.0024 Msun for Source A and 0.00033 Msun for Source B. C18O line data imply the existence of rotational motion around the extended structure, however, cannot resolve rotational motion in Source A and/or B, due to insufficient angular and frequency resolutions. Therefore, we cannot conclude whether Ser-emb 15 is a single or binary system. Thus, either Source A or B could drive the secondary outflow. We discuss two scenarios to explain the driving mechanism of the primary and secondary outflows: the Ser-emb 15 system is (1) a binary system composed of Source A and B or (2) a single star system composed of only Source A. In either case, the system could be a suitable target for investigating the disk and/or binary formation processes in complicated environments. Detecting these outflows should contribute to understanding complex star-forming environments, which may be common in the star-formation processes.
△ Less
Submitted 9 October, 2023;
originally announced October 2023.
-
Simulations of two-temperature jets in galaxy clusters: I. Effect of jet magnetization on dynamics and electron heating
Authors:
Takumi Ohmura,
Mami Machida
Abstract:
Non-radiating protons in the radio lobes have an essential role to form the jet morphology which is shown by recent radio and X-ray observations. However, since protons and electrons are not always in energy equilibrium due to weak Coulomb coupling, it is difficult to estimate the energy contribution of protons for inflation of radio lobes. The main focus of this study is to examine the effect of…
▽ More
Non-radiating protons in the radio lobes have an essential role to form the jet morphology which is shown by recent radio and X-ray observations. However, since protons and electrons are not always in energy equilibrium due to weak Coulomb coupling, it is difficult to estimate the energy contribution of protons for inflation of radio lobes. The main focus of this study is to examine the effect of the variable model for electron heating by turbulence and shock waves on the thermal energy distribution of electron and proton. We performed two-temperature three-dimensional magnetohydrodynamic simulations of sub-relativistic jets in the galaxy cluster while varying jet magnetization parameters. Because the energy partition rate between electrons and protons in shock and turbulence is determined by plasma kinetic scale physics, our global simulations include electron instantaneous heating sub-grid models for shock waves and turbulence. We find that most of the bulk kinetic energy of the jet is converted into thermal energy of protons through both shocks and turbulence. Thus, protons are energetically dominant. Meanwhile, thermal electrons stored in the lobe evolve toward energy equipartition with magnetic energy through turbulent dissipation. We further estimated the radio power and the mechanical jet power of radio lobes following the same method as used for radio and X-ray observations, and compared these powers with that of the observed radio jets. The two-temperature model quantitatively explains the radiatively inefficient radio cavities, but cannot reproduce the radiatively efficient cavity, even for strong magnetized jets. This implies that a significant population of pair-plasma is needed in the radiatively efficient radio cavities.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
Simulations of two-temperature jets in galaxy clusters: II. X-ray property of forward shock
Authors:
Takumi Ohmura,
Mami Machida,
Hiroki Akamatsu
Abstract:
Forward shocks by radio jets, driven into the intracluster medium, are one of the indicators that can be used to evaluate the power of the jet. Meanwhile high-angular-resolution X-ray observations show the Mach numbers of powerful radio jets are smaller compared to that of theoretical and numerical studies, $\mathcal{M_{\rm obs}} < 2$. Our aim is to systematically investigate various factors, such…
▽ More
Forward shocks by radio jets, driven into the intracluster medium, are one of the indicators that can be used to evaluate the power of the jet. Meanwhile high-angular-resolution X-ray observations show the Mach numbers of powerful radio jets are smaller compared to that of theoretical and numerical studies, $\mathcal{M_{\rm obs}} < 2$. Our aim is to systematically investigate various factors, such as projection effects and temperature non-equilibration between protons and electrons, that influence the Mach number estimate in a powerful jet. Using two-temperature magnetohydrodynamic simulation data for the Cygnus A radio jets, whose Mach number is approximately 6, we construct mock X-ray maps of simulated jets from various viewing angles. Further, we evaluate the shock Mach number from density/temperature jump using the same method of X-ray observations. Our results demonstrate that measurements from density jump significantly underestimate the Mach numbers, $\mathcal{M} < 2$, around the jet head at a low viewing angle, $\lessapprox 50^{\circ}$. The observed post-shock temperature is strongly reduced by the projection effect, as our jet is in the cluster center where the gas density is high. On the other hand, the temperature jump is almost unity, even if thermal electrons are in instant equilibration with protons. Upon comparison, we find that shock property of our model at viewing angle of $<$ $55^{\circ}$ is in a good agreement with that of Cygnus A observations.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
An ALMA-resolved view of 7000 au Protostellar Gas Ring around the Class I source CrA-IRS 2 as a possible sign of magnetic flux advection
Authors:
Kazuki Tokuda,
Naofumi Fukaya,
Kengo Tachihara,
Mitsuki Omura,
Naoto Harada,
Shingo Nozaki,
Ayumu Shoshi,
Masahiro N. Machida
Abstract:
Transferring a significant fraction of the magnetic flux from a dense cloud core is essential in the star formation process. A ring-like structure produced by magnetic flux loss has been predicted theoretically, but no observational identification has been presented. We have performed ALMA observations of the Class I protostar IRS 2 in the Corona Australis star-forming region and resolved a distin…
▽ More
Transferring a significant fraction of the magnetic flux from a dense cloud core is essential in the star formation process. A ring-like structure produced by magnetic flux loss has been predicted theoretically, but no observational identification has been presented. We have performed ALMA observations of the Class I protostar IRS 2 in the Corona Australis star-forming region and resolved a distinctive gas ring in the C$^{18}$O ($J$ = 2-1) line emission. The center of this gas ring is $\sim$5,000 au away from the protostar, with a diameter of $\sim$7,000 au. The radial velocity of the gas is $\lesssim1$ km s$^{-1}$ blueshifted from that of the protostar, with a possible expanding feature judged from the velocity-field (moment 1) map and position-velocity diagram. These features are either observationally new or have been discovered but not discussed in depth because they are difficult to explain by well-studied protostellar phenomena such as molecular outflows and accretion streamers. A plausible interpretation is a magnetic wall created by the advection of magnetic flux which is theoretically expected in the Class 0/I phase during star formation as a removal mechanism of magnetic flux. Similar structures reported in the other young stellar sources could likely be candidates formed by the same mechanism, encouraging us to revisit the issue of magnetic flux transport in the early stages of star formation from an observational perspective.
△ Less
Submitted 15 October, 2023; v1 submitted 24 September, 2023;
originally announced September 2023.
-
Exact Calculation of Nonideal Fields Demonstrates Dominance of Injection in Relativistic Reconnection
Authors:
Samuel R. Totorica,
Seiji Zenitani,
Shuichi Matsukiyo,
Mami Machida,
Kazuhiro Sekiguchi,
Amitava Bhattacharjee
Abstract:
Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-…
▽ More
Magnetic reconnection is an important source of energetic particles in systems ranging from astrophysics to the laboratory. The large separation of spatiotemporal scales involved makes it critical to determine the minimum physical model containing the necessary physics for modeling particle acceleration. By resolving the energy gain from ideal and nonideal magnetohydrodynamic electric fields self-consistently in kinetic particle-in-cell simulations of reconnection, we conclusively show the dominant role of the nonideal field for the early stage of energization known as injection. The importance of the nonideal field increases with magnetization, guide field, and in three-dimensions, indicating its general importance for reconnection in natural astrophysical systems. We obtain the statistical properties of the injection process from the simulations, paving the way for the development of extended MHD models capable of accurately modeling particle acceleration in large-scale systems. The novel analysis method developed in this study can be applied broadly to give new insight into a wide range of processes in plasma physics.
△ Less
Submitted 15 September, 2023;
originally announced September 2023.
-
An ALMA Glimpse of Dense Molecular Filaments Associated with High-mass Protostellar Systems in the Large Magellanic Cloud
Authors:
Kazuki Tokuda,
Naoto Harada,
Kei E. I. Tanaka,
Tsuyoshi Inoue,
Takashi Shimonishi,
Yichen Zhang,
Marta Sewiło,
Yuri Kunitoshi,
Ayu Konishi,
Yasuo Fukui,
Akiko Kawamura,
Toshikazu Onishi,
Masahiro N. Machida
Abstract:
Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3)…
▽ More
Recent millimeter/sub-millimeter facilities have revealed the physical properties of filamentary molecular clouds in relation to high-mass star formation. A uniform survey of the nearest, face-on star-forming galaxy, the Large Magellanic Cloud (LMC), complements the Galactic knowledge. We present ALMA survey data with a spatial resolution of $\sim$0.1 pc in the 0.87 mm continuum and HCO$^{+}$(4-3) emission toward 30 protostellar objects with luminosities of 10$^4$-10$^{5.5}$ $L_{\odot}$ in the LMC. The spatial distributions of the HCO$^{+}$(4-3) line and thermal dust emission are well correlated, indicating that the line effectively traces dense, filamentary gas with an H$_2$ volume density of $\gtrsim$10$^5$ cm$^{-3}$ and a line mass of $\sim$10$^3$-10$^{4}$ $M_{\odot}$ pc$^{-1}$. Furthermore, we obtain an increase in the velocity linewidths of filamentary clouds, which follows a power-law dependence on their H$_2$ column densities with an exponent of $\sim$0.5. This trend is consistent with observations toward filamentary clouds in nearby star-forming regions withiin $ \lesssim$1 kpc from us and suggests enhanced internal turbulence within the filaments owing to surrounding gas accretion. Among the 30 sources, we find that 14 are associated with hub-filamentary structures, and these complex structures predominantly appear in protostellar luminosities exceeding $\sim$5 $\times$10$^4$ $L_{\odot}$. The hub-filament systems tend to appear in the latest stages of their natal cloud evolution, often linked to prominent H$\;${\sc ii} regions and numerous stellar clusters. Our preliminary statistics suggest that the massive filaments accompanied by hub-type complex features may be a necessary intermediate product in forming extremely luminous high-mass stellar systems capable of ultimately dispersing the parent cloud.
△ Less
Submitted 10 August, 2023;
originally announced August 2023.
-
First observations of warm and cold methanol in Class 0/I proto-brown dwarfs
Authors:
B. Riaz,
W. -F. Thi,
M. N. Machida
Abstract:
We present results from the first molecular line survey to search for the fundamental complex organic molecule, methanol (CH$_{3}$OH), in 14 Class 0/I proto-brown dwarfs (proto-BDs). IRAM 30-m observations over the frequency range of 92-116 GHz and 213-280 GHz have revealed emission in 14 CH$_{3}$OH transition lines, at upper state energy level, E$_{upper}\sim$7-49 K, and critical densities,…
▽ More
We present results from the first molecular line survey to search for the fundamental complex organic molecule, methanol (CH$_{3}$OH), in 14 Class 0/I proto-brown dwarfs (proto-BDs). IRAM 30-m observations over the frequency range of 92-116 GHz and 213-280 GHz have revealed emission in 14 CH$_{3}$OH transition lines, at upper state energy level, E$_{upper}\sim$7-49 K, and critical densities, $n_{crit}$ of 10$^{5}$ to 10$^{9}$ cm$^{-3}$. The most commonly detected lines are at E$_{upper} <$ 20 K, while 11 proto-BDs also show emission in the higher excitation lines at E$_{upper}\sim$21-49 K and $n_{crit}\sim$10$^{5}$ to 10$^{8}$ cm$^{-3}$. In comparison with the brown dwarf formation models, the high excitation lines likely probe the warm ($\sim$25-50 K) corino region at $\sim$10-50 au in the proto-BDs, while the low-excitation lines trace the cold ($<$ 20 K) gas at $\sim$50-150 au. The column density for the cold component is an order of magnitude higher than the warm component. The CH$_{3}$OH ortho-to-para ratios range between $\sim$0.3-2.3. The volume-averaged CH$_{3}$OH column densities show a rise with decreasing bolometric luminosity among the proto-BDs, with the median column density higher by a factor of $\sim$3 compared to low-mass protostars. Emission in high-excitation (E$_{upper}>$ 25 K) CH$_{3}$OH lines together with the model predictions suggest that a warm corino is present in $\sim$78\% of the proto-BDs in our sample. The remaining show evidence of only the cold component, possibly due to the absence of a strong, high-velocity jet that can stir up the warm gas around it.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
Non-thermal emissions from a head-tail radio galaxy in 3D magnetohydrodynamic simulations
Authors:
Takumi Ohmura,
Katsuaki Asano,
Kosuke Nishiwaki,
Mami Machida,
Haruka Sakemi
Abstract:
We present magnetohydrodynamic simulations of a jet-wind interaction in a galaxy cluster and the radio to gamma-ray and the neutrino emissions from this "head-tail galaxy". Our simulation follows the evolution of cosmic-ray (CR) particle spectra with energy losses and the stochastic turbulence acceleration. We find that the reacceleration is essential to explain the observed radio properties of he…
▽ More
We present magnetohydrodynamic simulations of a jet-wind interaction in a galaxy cluster and the radio to gamma-ray and the neutrino emissions from this "head-tail galaxy". Our simulation follows the evolution of cosmic-ray (CR) particle spectra with energy losses and the stochastic turbulence acceleration. We find that the reacceleration is essential to explain the observed radio properties of head-tail galaxies, in which the radio flux and spectral index do not drastically change. Our models suggest that hard X-ray emissions can be detected around the head-tail galaxy in the Perseus cluster by the hard X-ray satellites, such as FORCE, and it will potentially constrain the acceleration efficiency. We also explore the origin of the collimated synchrotron threads, which are found in some head-tail galaxies by recent high-quality radio observations. Thin and elongated flux tubes, connecting the two tails, are formed by strong backflows at an early phase. We find that these threads advect with the wind for over 300 Myr without disrupting. The radio flux from the flux tubes is much lower than the typical observed flux. An efficient CR diffusion process along the flux tubes, however, may solve this discrepancy.
△ Less
Submitted 8 May, 2023;
originally announced May 2023.
-
Impact of turbulence intensity and fragmentation velocity on dust particle size evolution and non-ideal magnetohydrodynamics effects
Authors:
Yoshihiro Kawasaki,
Masahiro N. Machida
Abstract:
We investigate the influence of dust particle size evolution on non-ideal magnetohydrodynamic effects during the collapsing phase of star-forming cores, taking both the turbulence intensity in the collapsing cloud core and the fragmentation velocity of dust particles as parameters. When the turbulence intensity is small, the dust particles do not grow significantly, and the non-ideal MHD effects w…
▽ More
We investigate the influence of dust particle size evolution on non-ideal magnetohydrodynamic effects during the collapsing phase of star-forming cores, taking both the turbulence intensity in the collapsing cloud core and the fragmentation velocity of dust particles as parameters. When the turbulence intensity is small, the dust particles do not grow significantly, and the non-ideal MHD effects work efficiently in high-density regions. The dust particles rapidly grow in a strongly turbulent environment, while the efficiency of non-ideal MHD effects in such an environment depends on the fragmentation velocity of the dust particles. When the fragmentation velocity is small, turbulence promotes coagulation growth and collisional fragmentation of dust particles, producing small dust particles. In this case, the adsorption of charged particles on the dust particle surfaces becomes efficient and the abundance of charged particles decreases, making non-ideal MHD effects effective at high densities. On the other hand, when the fragmentation velocity is high, dust particles are less likely to fragment, even if the turbulence is strong. In this case, the production of small dust particles become inefficient and non-ideal MHD effects become less effective. We also investigate the effect of the dust composition on the star and disk formation processes. We constrain the turbulence intensity of a collapsing core and the fragmentation velocity of dust for circumstellar disk formation due to the dissipation of the magnetic field.
△ Less
Submitted 25 April, 2023;
originally announced April 2023.
-
Co-evolution of dust grains and protoplanetary disks
Authors:
Yusuke Tsukamoto,
Masahiro N. Machida,
Shu-ichiro Inutsuka
Abstract:
We propose a new evolutionary process of protoplanetary disks "co-evolution of dust grains and protoplanetary disks", revealed by dust-gas two-fluid non-ideal magnetohydrodynamics simulations considering the growth of dust and associated changes in magnetic resistivity. We found that the dust growth significantly affects disk evolution by changing the coupling between the gas and magnetic field. M…
▽ More
We propose a new evolutionary process of protoplanetary disks "co-evolution of dust grains and protoplanetary disks", revealed by dust-gas two-fluid non-ideal magnetohydrodynamics simulations considering the growth of dust and associated changes in magnetic resistivity. We found that the dust growth significantly affects disk evolution by changing the coupling between the gas and magnetic field. Moreover, once the dust grains sufficiently grow and the adsorption of charged particles on dust grains becomes negligible, the physical quantities (e.g., density and magnetic field) of the disk are well described by characteristic power laws. In this disk structure, the radial profile of density is steeper and the disk mass is smaller than those of the model ignoring dust growth. We analytically derive these power laws from the basic equations of non-ideal magnetohydrodynamics. The analytical power laws are determined only by observable physical quantities, e.g., central stellar mass and mass accretion rate, and do not include difficult-to-determine parameters e.g., viscous parameter $α$. Therefore, our model is observationally testable and this disk structure is expected to provide a new perspective for future studies on protostar and disk evolution.
△ Less
Submitted 14 June, 2023; v1 submitted 18 March, 2023;
originally announced March 2023.
-
Pseudo-observation of spiral galaxies in the radio band to verify depolarization models
Authors:
Yuta Tashima,
Takumi Ohmura,
Mami Machida
Abstract:
Studies of the three-dimensional (3D) structures of galactic magnetic fields are now entering a new era, with broadband, highly sensitive radio observations and new analysis methods. To reveal the magnetic field configuration from the observed value integrated along the line of sight, it is necessary to derive an appropriate model that can reproduce the observational characteristics. We aim to cla…
▽ More
Studies of the three-dimensional (3D) structures of galactic magnetic fields are now entering a new era, with broadband, highly sensitive radio observations and new analysis methods. To reveal the magnetic field configuration from the observed value integrated along the line of sight, it is necessary to derive an appropriate model that can reproduce the observational characteristics. We aim to clarify the relationship between the radiation field and the spatial distribution of physical quantities through pseudo-observations using global 3D magnetohydrodynamics (MHD) simulation results. In particular, we focus on using the depolarization effect to verify the polarization model and to identify the emission region. First, we show that wavelength-independent depolarization, which takes into account anisotropic turbulence, does not work efficiently because the polarized emission is stronger in regions of ordered spiral fields than in regions dominated by isotropic turbulent fields. Beam depolarization becomes more effective below 1 GHz. Although close to the equatorial plane there will be strong depolarization which increases with observing wavelength, this effect is less in the halo, making halo magnetic fields detectable through their polarized emission at meter-wavelength bands. Although polarized emission from the halo is below the detection limit of current facilities, it will be detectable within the SKA era. In addition, we found that the spiral polarization projected on a screen is produced by overlapping magnetic flux tubes extending to different heights from the equatorial plane. This suggests that the traditional classification of global magnetic fields has difficulty reproducing the global structure of the magnetic fields. Finally, we demonstrate the method to separate magnetic flux tubes at different heights by using peak frequencies which cause the decreasing of polarized flux.
△ Less
Submitted 5 March, 2023;
originally announced March 2023.
-
Direct Imaging Explorations for Companions around Mid-Late M Stars from the Subaru/IRD Strategic Program
Authors:
Taichi Uyama,
Charles Beichman,
Masayuki Kuzuhara,
Markus Janson,
Takayuki Kotani,
Dimitri Mawet,
Bun'ei Sato,
Motohide Tamura,
Hiroyuki Tako Ishikawa,
Bryson Cale,
Thayne Currie,
Hiroki Harakawa,
Thomas Henning,
Teruyuki Hirano,
Klaus Hodapp,
Yasunori Hori,
Masato Ishizuka,
Shane Jacobson,
Yui Kasagi,
Eiichiro Kokubo,
Mihoko Konishi,
Tomoyuki Kudo,
Takashi Kurokawa,
Nobuhiko Kusakabe,
Jungmi Kwon
, et al. (10 additional authors not shown)
Abstract:
The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this…
▽ More
The Subaru telescope is currently performing a strategic program (SSP) using the high-precision near-infrared (NIR) spectrometer IRD to search for exoplanets around nearby mid/late-M~dwarfs via radial velocity (RV) monitoring. As part of the observing strategy for the exoplanet survey, signatures of massive companions such as RV trends are used to reduce the priority of those stars. However, this RV information remains useful for studying the stellar multiplicity of nearby M~dwarfs. To search for companions around such ``deprioritized" M~dwarfs, we observed 14 IRD-SSP targets using Keck/NIRC2 observations with pyramid wavefront sensing at NIR wavelengths, leading to high sensitivity to substellar-mass companions within a few arcseconds. We detected two new companions (LSPM~J1002+1459~B and LSPM~J2204+1505~B) and two new candidates that are likely companions (LSPM~J0825+6902~B and LSPM~J1645+0444~B) as well as one known companion. Including two known companions resolved by the IRD fiber injection module camera, we detected seven (four new) companions at projected separations between $\sim2-20$~au in total. A comparison of the colors with the spectral library suggests that LSPM~J2204+1505~B and LSPM~J0825+6902~B are located at the boundary between late-M and early-L spectral types. Our deep high-contrast imaging for targets where no bright companions were resolved did not reveal any additional companion candidates. The NIRC2 detection limits could constrain potential substellar-mass companions ($\sim10-75\ M_{\rm Jup}$) at 10~au or further. The failure with Keck/NIRC2 around the IRD-SSP stars having significant RV trends makes these objects promising targets for further RV monitoring or deeper imaging with JWST to search for smaller-mass companions below the NIRC2 detection limits.
△ Less
Submitted 26 February, 2023;
originally announced February 2023.
-
Crescent-Shaped Molecular Outflow from the Intermediate-mass Protostar DK Cha Revealed by ALMA
Authors:
Naoto Harada,
Kazuki Tokuda,
Hayao Yamasaki,
Asako Sato,
Mitsuki Omura,
Shingo Hirano,
Toshikazu Onishi,
Kengo Tachihara,
Masahiro N. Machida
Abstract:
We report on an Atacama Large Millimeter/submillimeter Array (ALMA) study of the Class I or II intermediate-mass protostar DK Cha in the Chamaeleon II region. The 12CO (J=2-1) images have an angular resolution of ~1'' (~250 au) and show high-velocity blueshifted (>70 km s-1) and redshifted (>50 km s-1) emissions which have 3000 au scale crescent-shaped structures around the protostellar disk trace…
▽ More
We report on an Atacama Large Millimeter/submillimeter Array (ALMA) study of the Class I or II intermediate-mass protostar DK Cha in the Chamaeleon II region. The 12CO (J=2-1) images have an angular resolution of ~1'' (~250 au) and show high-velocity blueshifted (>70 km s-1) and redshifted (>50 km s-1) emissions which have 3000 au scale crescent-shaped structures around the protostellar disk traced in the 1.3mm continuum. Because the high-velocity components of the CO emission are associated with the protostar, we concluded that the emission traces the pole-on outflow. The blueshifted outflow lobe has a clear layered velocity gradient with a higher velocity component located on the inner side of the crescent shape, which can be explained by a model of an outflow with a higher velocity in the inner radii. Based on the directly driven outflow scenario, we estimated the driving radii from the observed outflow velocities and found that the driving region extends over two orders of magnitude. The 13CO emission traces a complex envelope structure with arc-like substructures with lengths of ~1000au. We identified the arc-like structures as streamers because they appear to be connected to a rotating infalling envelope. DK Cha is useful for understanding characteristics that are visible by looking at nearly face-on configurations of young protostellar systems, providing an alternative perspective for studying the star-formation process.
△ Less
Submitted 3 February, 2023;
originally announced February 2023.
-
Molecular clouds at the eastern edge of radio nebula W50
Authors:
Haruka Sakemi,
Mami Machida,
Hiroaki Yamamoto,
Kengo Tachihara
Abstract:
Microquasar SS 433 located at the geometric center of radio nebula W50 is a suitable source for investigating the physical process of how galactic jets affect the surrounding interstellar medium (ISM). Previous studies have searched for evidence of the interaction between the SS 433 jet and ISM, such as neutral hydrogen gas and molecular clouds; however, it is still unclear which ISM interacts wit…
▽ More
Microquasar SS 433 located at the geometric center of radio nebula W50 is a suitable source for investigating the physical process of how galactic jets affect the surrounding interstellar medium (ISM). Previous studies have searched for evidence of the interaction between the SS 433 jet and ISM, such as neutral hydrogen gas and molecular clouds; however, it is still unclear which ISM interacts with the jet. We looked for new molecular clouds that possibly interact at the terminal of the SS 433 eastern jet using the Nobeyama 45-m telescope and the Atacama Submillimeter Telescope Experiment (ASTE). We identified two molecular clouds, comprising many small clumps, in the velocity range of 30.1--36.5 km s$^{-1}$ for the first time. These clouds have complex velocity structures, and one of them has a density gradient toward SS 433. Although it is difficult to conclude the relation between the molecular clouds and the SS 433/W50 system, there is a possibility that the eastern structure of W50 constructed by the SS 433 jet swept up tiny molecular clumps drifting in the surroundings and formed the molecular clouds that we identified in this study.
△ Less
Submitted 30 January, 2023;
originally announced January 2023.
-
Environmental Effects of Star-Forming Cores on Mass Accretion Rate
Authors:
Shingo Nozaki,
Masahiro N. Machida
Abstract:
We calculate the evolution of cloud cores embedded in different envelopes to investigate environmental effects on the mass accretion rate onto protostars. As the initial state, we neglect the magnetic field and cloud rotation, and adopt star-forming cores composed of two parts: a centrally condensed core and an outer envelope. The inner core has a critical Bonnor-Ebert density profile and is enclo…
▽ More
We calculate the evolution of cloud cores embedded in different envelopes to investigate environmental effects on the mass accretion rate onto protostars. As the initial state, we neglect the magnetic field and cloud rotation, and adopt star-forming cores composed of two parts: a centrally condensed core and an outer envelope. The inner core has a critical Bonnor-Ebert density profile and is enclosed by the outer envelope. We prepare 15 star-forming cores with different outer envelope densities and gravitational radii, within which the gas flows into the collapsing core, and calculate their evolution until $\sim 2 \times10^5$ yr after protostar formation. The mass accretion rate decreases as the core is depleted when the outer envelope density is low. In contrast, the mass accretion rate is temporarily enhanced when the outer envelope density is high and the resultant protostellar mass exceeds the initial mass of the centrally condensed core. Some recent observations indicate that the mass of prestellar cores is too small to reproduce the stellar mass distribution. Our simulations show that the mass inflow from outside the core contributes greatly to protostellar mass growth when the core is embedded in a high-density envelope, which could explain the recent observations.
△ Less
Submitted 26 December, 2022;
originally announced December 2022.
-
Dust Motion and Possibility of Dust Growth in a Growing Circumstellar Disk
Authors:
Shunta Koga,
Masahiro N. Machida
Abstract:
We calculate the evolution of a star-forming cloud core using a three-dimensional resistive magnetohydrodynamics simulation, treating dust grains as Lagrangian particles, to investigate the dust motion in the early star formation stage. We prepare six different-sized set of dust particles in the range $a_{\rm d}=0.01$--$1000\,μ$m, where $a_{\rm d}$ is the dust grain size. In a gravitationally coll…
▽ More
We calculate the evolution of a star-forming cloud core using a three-dimensional resistive magnetohydrodynamics simulation, treating dust grains as Lagrangian particles, to investigate the dust motion in the early star formation stage. We prepare six different-sized set of dust particles in the range $a_{\rm d}=0.01$--$1000\,μ$m, where $a_{\rm d}$ is the dust grain size. In a gravitationally collapsing cloud, a circumstellar disk forms around a protostar and drives a protostellar outflow. Almost all the small dust grains ($a_{\rm d} \lesssim 10$--$100\,μ$m) initially distributed in the region $θ_0 \lesssim 45^\circ$ are ejected from the center by the outflow, where $θ_0$ is the initial zenith angle relative to the rotation axis, whereas only a small number of the large dust grains ($a_{\rm d} \gtrsim 100\,μ$m) distributed in the region are ejected. All other grains fall onto either the protostar or disk without being ejected by the outflow. Regardless of the dust grain size, the behavior of the dust motion is divided into two trends after dust particles settle into the circumstellar disk. The dust grains reaching the inner disk region from the upper envelope preferentially fall onto the protostar, while those reaching the outer disk region or disk outer edge from the envelope can survive without an inward radial drift. These surviving grains can induce dust growth. Thus, we expect that the outer disk regions could be a favored place of planet formation.
△ Less
Submitted 28 November, 2022;
originally announced November 2022.
-
ALMA Fragmented Source Catalogue in Orion (FraSCO) I. Outflow interaction within an embedded cluster in OMC-2/FIR3, FIR4, and FIR5
Authors:
Asako Sato,
Satoko Takahashi,
Shun Ishii,
Paul T. P. Ho,
Masahiro N. Machida,
John Carpenter,
Luis A. Zapata,
Paula Stella Teixeira,
Sümeyye Suri
Abstract:
We present a high angular resolution ($\sim1"$) and wide-field ($2'.9 \times 1'.9$) image of the 1.3-mm continuum, CO($J$ = 2--1) line, and SiO($J$ = 5--4) line emissions toward an embedded protocluster, FIR3, FIR4, and FIR5, in the Orion Molecular Cloud 2 obtained from the Atacama Large Millimeter/submillimeter Array (ALMA). We identify 51 continuum sources, 36 of which are newly identified in th…
▽ More
We present a high angular resolution ($\sim1"$) and wide-field ($2'.9 \times 1'.9$) image of the 1.3-mm continuum, CO($J$ = 2--1) line, and SiO($J$ = 5--4) line emissions toward an embedded protocluster, FIR3, FIR4, and FIR5, in the Orion Molecular Cloud 2 obtained from the Atacama Large Millimeter/submillimeter Array (ALMA). We identify 51 continuum sources, 36 of which are newly identified in this study. Their dust masses, projected sizes, and $\mathrm{H_2}$ gas number densities are estimated to be $3.8 \times 10^{-5}$--$ 1.1 \times 10^{-2} \mathrm{M_{\odot}}$, 290--2000 au, and $6.4 \times 10^{6}$--$3.3 \times 10^{8}\,\mathrm{cm^{-3}}$, respectively. The results of a Jeans analysis show that $\sim80\,\%$ of the protostellar sources and $\sim15\,\%$ of the prestellar sources are gravitationally bound. We identify 12 molecular outflows traced in the CO($J$ = 2--1) emission, six of which are newly detected. We spatially resolve shocked gas structures traced by the SiO($J$ = 5--4) emission in this region for the first time. We identify shocked gas originating from outflows and other shocked regions. These results provide direct evidence of an interaction between a dust condensation, FIR4, and an energetic outflow driven by HOPS-370 located within FIR3. A comparison of the outflow dynamical timescales, fragmentation timescales, and protostellar ages shows that the previously proposed triggered star-formation scenario in FIR4 is not strongly supported. We also discuss the spatial distribution of filaments identified in our continuum image by comparing it with a previously identified hub-fiber system in the $\mathrm{N_2H^+}$ line.
△ Less
Submitted 23 November, 2022; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Twisted magnetic field in star formation processes of L1521 F revealed by submillimeter dual band polarimetry using James Clerk Maxwell Telescope
Authors:
Sakiko Fukaya,
Hiroko Shinnaga,
Ray S. Furuya,
Kohji Tomisaka,
Masahiro N. Machida,
Naoto Harada
Abstract:
Understanding the initial conditions of star formation requires both observational studies and theoretical works taking into account the magnetic field, which plays an important role in star formation processes. Herein, we study the young nearby dense cloud core L1521 F ($n$(H$_2$) $\sim 10^{4-6}$ cm$^{-3}$) in the Taurus Molecular Cloud. This dense core hosts a 0.2 $M_\odot$ protostar, categorize…
▽ More
Understanding the initial conditions of star formation requires both observational studies and theoretical works taking into account the magnetic field, which plays an important role in star formation processes. Herein, we study the young nearby dense cloud core L1521 F ($n$(H$_2$) $\sim 10^{4-6}$ cm$^{-3}$) in the Taurus Molecular Cloud. This dense core hosts a 0.2 $M_\odot$ protostar, categorized as a Very Low Luminosity Objects with complex velocity structures, particularly in the vicinity of the protostar. To trace the magnetic field within the dense core, we conducted high sensitivity submillimeter polarimetry of the dust continuum at $λ$= 850 $μ$m and 450 $μ$m using the POL-2 polarimeter situated in front of the SCUBA-2 submillimeter bolometer camera on James Clerk Maxwell Tetescope. This was compared with millimeter polarimetry taken at $λ$= 3.3 mm with ALMA. The magnetic field was detected at $λ$= 850 $μ$m in the peripheral region, which is threaded in a north-south direction, while the central region traced at $λ$= 450 $μ$m shows a magnetic field with an east-west direction, i.e., orthogonal to that of the peripheral region. Magnetic field strengths are estimated to be $\sim$70 $μ$G and 200 $μ$G in the peripheral- and central-regions, respectively, using the Davis-Chandrasekhar-Fermi method. The resulting mass-to-flux ratio of 3 times larger than that of magnetically critical state for both regions indicates that L1521 F is magnetically supercritical, i.e., gravitational forces dominate over magnetic turbulence forces. Combining observational data with MHD simulations, detailed parameters of the morphological properties of this puzzling object are derived for the first time.
△ Less
Submitted 16 November, 2022;
originally announced November 2022.
-
The role of magnetic fields in the formation of protostars, disks, and outflows
Authors:
Yusuke Tsukamoto,
Anaëlle Maury,
Benoît Commerçon,
Felipe O. Alves,
Erin G. Cox,
Nami Sakai,
Tom Ray,
Bo Zhao,
Masahiro N. Machida
Abstract:
We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. A…
▽ More
We present our current understanding of the formation and early evolution of protostars, protoplanetary disks, and the driving of outflows as dictated by the interplay of magnetic fields and partially ionized gas in molecular cloud cores. In recent years, the field has witnessed enormous development through sub-millimeter observations which in turn have constrained models of protostar formation. As a result of these observations % that the observations provided, the state-of-the-art theoretical understanding of the formation and evolution of young stellar objects is described. In particular, we emphasize the importance of the coupling, decoupling, and re-coupling between weakly ionized gas and the magnetic field on appropriate scales. This highlights the complex and intimate relationship between gravitational collapse and magnetic fields in young protostars.
△ Less
Submitted 27 September, 2022;
originally announced September 2022.
-
Magnetic Effects Promote Supermassive Star Formation in Metal-enriched Atomic-cooling Halos
Authors:
Shingo Hirano,
Masahiro N. Machida,
Shantanu Basu
Abstract:
Intermediate-mass black holes (with $\geq\!10^5\,M_\odot$) are promising candidates for the origin of supermassive black holes (with $\sim\!10^9\,M_\odot$) in the early universe (redshift $z\sim6$). Chon & Omukai (2020) firstly pointed out the direct collapse black hole (DCBH) formation in metal-enriched atomic-cooling halos (ACHs), which relaxes the DCBH formation criterion. On the other hand, Hi…
▽ More
Intermediate-mass black holes (with $\geq\!10^5\,M_\odot$) are promising candidates for the origin of supermassive black holes (with $\sim\!10^9\,M_\odot$) in the early universe (redshift $z\sim6$). Chon & Omukai (2020) firstly pointed out the direct collapse black hole (DCBH) formation in metal-enriched atomic-cooling halos (ACHs), which relaxes the DCBH formation criterion. On the other hand, Hirano et al. (2021) showed that the magnetic effects promote the DCBH formation in metal-free ACHs. We perform a set of magnetohydrodynamical simulations to investigate star formation in the magnetized ACHs with metallicities $Z/Z_\odot = 0$, $10^{-5}$, and $10^{-4}$. Our simulations show that the mass accretion rate onto the protostars becomes lower in metal-enriched ACHs than that of metal-free ACHs. However, many protostars form from gravitationally and thermally unstable metal-enriched gas clouds. Under such circumstances, the magnetic field rapidly increases as the magnetic field lines wind up due to the spin of protostars. The region with the amplified magnetic field expands outwards due to the orbital motion of protostars and the rotation of the accreting gas. The amplified magnetic field extracts the angular momentum from the accreting gas, promotes the coalescence of the low-mass protostars, and increases the mass growth rate of the primary protostar. We conclude that the magnetic field amplification is always realized in the metal-enriched ACHs regardless of the initial magnetic field strength, which affects the DCBH formation criterion. In addition, we find a qualitatively different trend from the previous unmagnetized simulations in that the mass growth rate is maximal for the extremely metal-poor ACHs with $Z/Z_\odot = 10^{-5}$.
△ Less
Submitted 13 June, 2023; v1 submitted 8 September, 2022;
originally announced September 2022.
-
Exponentially amplified magnetic field eliminates disk fragmentation around the Population III protostar
Authors:
Shingo Hirano,
Masahiro N. Machida
Abstract:
One critical remaining issue to unclear the initial mass function of the first (Population III) stars is the final fate of secondary protostars formed in the accretion disk, specifically whether they merge or survive. We focus on the magnetic effects on the first star formation under the cosmological magnetic field. We perform a suite of ideal magnetohydrodynamic simulations until 1000 years after…
▽ More
One critical remaining issue to unclear the initial mass function of the first (Population III) stars is the final fate of secondary protostars formed in the accretion disk, specifically whether they merge or survive. We focus on the magnetic effects on the first star formation under the cosmological magnetic field. We perform a suite of ideal magnetohydrodynamic simulations until 1000 years after the first protostar formation. Instead of the sink particle technique, we employ a stiff equation of state approach to represent the magnetic field structure connecting to protostars. Ten years after the first protostar formation in the cloud initialized with $B_0 = 10^{-20}$ G at $n_0 = 10^4\,{\rm cm^{-3}}$, the magnetic field strength around protostars amplifies from pico- to kilo-gauss, which is the same strength as the present-day star. The magnetic field rapidly winds up since the gas in the vicinity of the protostar ($\leq\!10$ au) has undergone several tens orbital rotations in the first decade after protostar formation. As the mass accretion progresses, the vital magnetic field region extends outward, and the magnetic braking eliminates fragmentation of the disk that would form in the unmagnetized model. On the other hand, assuming a gas cloud with small angular momentum, this amplification might not work because the rotation would be slower. However, disk fragmentation would not occur in that case. We conclude that the exponential amplification of the cosmological magnetic field strength, about $10^{-18}$ G, eliminates disk fragmentation around the Population III protostars.
△ Less
Submitted 1 August, 2022;
originally announced August 2022.
-
Implementation of dust particles in three-dimensional magnetohydrodynamics simulation: Dust dynamics in a collapsing cloud core
Authors:
Shunta Koga,
Yoshihiro Kawasaki,
Masahiro N. Machida
Abstract:
The aim of this study is to examine dust dynamics on a large scale and investigate the coupling of dust with gas fluid in the star formation process. We propose a method for calculating the dust trajectory in a gravitationally collapsing cloud, where the dust grains are treated as Lagrangian particles and are assumed to be neutral. We perform the dust trajectory calculations in combination with no…
▽ More
The aim of this study is to examine dust dynamics on a large scale and investigate the coupling of dust with gas fluid in the star formation process. We propose a method for calculating the dust trajectory in a gravitationally collapsing cloud, where the dust grains are treated as Lagrangian particles and are assumed to be neutral. We perform the dust trajectory calculations in combination with non-ideal magnetohydrodynamics simulation. Our simulation shows that dust particles with a size of $\le 10\,{\rm μm}$ are coupled with gas in a star-forming cloud core. We investigate the time evolution of the dust-to-gas mass ratio and the Stokes number, which is defined as the stopping time normalized by the freefall time-scale, and show that large dust grains ($\gtrsim 100\,{\rm μm}$) have a large Stokes number (close to unity) and tend to concentrate in the central region (i.e., protostar and rotationally supported disk) faster than do small grains ($\lesssim 10\,{\rm μm}$). Thus, large grains significantly increase the dust-to-gas mass ratio around and inside the disk. We also confirm that the dust trajectory calculations, which trace the physical quantities of each dust particle, reproduce previously reported results obtained using the Eulerian approach.
△ Less
Submitted 26 July, 2022;
originally announced July 2022.
-
The First Detection of a Protostellar CO Outflow in the Small Magellanic Cloud with ALMA
Authors:
Kazuki Tokuda,
Sarolta Zahorecz,
Yuri Kunitoshi,
Kosuke Higashino,
Kei E. I. Tanaka,
Ayu Konishi,
Taisei Suzuki,
Naoya Kitano,
Naoto Harada,
Takashi Shimonishi,
Naslim Neelamkodan,
Yasuo Fukui,
Akiko Kawamura,
Toshikazu Onishi,
Masahiro N. Machida
Abstract:
Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularl…
▽ More
Protostellar outflows are one of the most outstanding features of star formation. Observational studies over the last several decades have successfully demonstrated that outflows are ubiquitously associated with low- and high-mass protostars in the solar-metallicity Galactic conditions. However, the environmental dependence of protostellar outflow properties is still poorly understood, particularly in the low-metallicity regime. Here we report the first detection of a molecular outflow in the Small Magellanic Cloud with 0.2 $Z_{\odot}$, using Atacama Large Millimeter/submillimeter Array observations at a spatial resolution of 0.1 pc toward the massive protostar Y246. The bipolar outflow is nicely illustrated by high-velocity wings of CO(3-2) emission at $\gtrsim$15 km s$^{-1}$. The evaluated properties of the outflow (momentum, mechanical force, etc.) are consistent with those of the Galactic counterparts. Our results suggest that the molecular outflows, i.e., the guidepost of the disk accretion at the small scale, might be universally associated with protostars across the metallicity range of $\sim$0.2-1 $Z_{\odot}$.
△ Less
Submitted 7 August, 2022; v1 submitted 18 July, 2022;
originally announced July 2022.
-
Effects of Depolarizing Intervening Galaxies on Background Radio Emission I. Global Disk Magnetic Field
Authors:
Rikuto Omae,
Takuya Akahori,
Mami Machida
Abstract:
External galaxies often intervene in front of background radio sources such as quasars and radio galaxies. Linear polarization of the background emission is depolarized by Faraday rotation of inhomogeneous magnetized plasma of the intervening galaxies. Exploring the depolarizing intervening galaxies (DINGs) can be a powerful tool to investigate the cosmological evolution of the galactic magnetic f…
▽ More
External galaxies often intervene in front of background radio sources such as quasars and radio galaxies. Linear polarization of the background emission is depolarized by Faraday rotation of inhomogeneous magnetized plasma of the intervening galaxies. Exploring the depolarizing intervening galaxies (DINGs) can be a powerful tool to investigate the cosmological evolution of the galactic magnetic field. In this paper, we investigate the effects of DINGs on background radio emission using theoretical DING models. We find that complex structures of galaxy result in complicated depolarization features and the Faraday dispersion functions (FDFs), but, for the features of depolarizations and FDFs, the global component of magnetic fields is important. We show the simplest results with ring magnetic field in the galactic disk. We find that the degree of depolarization significantly depends on the inclination angle and the impact parameter of the DING. We found that the larger the standard deviation, the more likely it is that depolarization will occur. The FDF represents the RM structure within the beam. The FDF exhibits multi-components due mainly to the RM structure within the beam and the fraction of the DING that covers the background emission (the filling factor). The peak Faraday depth of the FDF is different from the beam-averaged RM of the DING. The Monte-Carlo simulations indicate that DING's contribution to the standard deviation of observed RMs follows $σ_{\rm RM} \propto 1/{(1+z)^k}$ with $k \sim 2.7$ and exhibits a steeper redshift dependence than the wavelength squared. DINGs will have a significant impact on RM catalogs created by future survey projects such as the SKA and SKA Precursor/Pathfinder.
△ Less
Submitted 15 July, 2022;
originally announced July 2022.
-
Delivery of gas onto the circumplanetary disk of giant planets: Planetary-mass dependence of the source region of accreting gas and mass accretion rate
Authors:
Natsuho Maeda,
Keiji Ohtsuki,
Takayuki Tanigawa,
Masahiro N. Machida,
Ryo Suetsugu
Abstract:
Gas accretion onto the circumplanetary disks and the source region of accreting gas are important to reveal dust accretion that leads to satellite formation around giant planets. We performed local three-dimensional high-resolution hydrodynamic simulations of isothermal and inviscid gas flow around a planet to investigate planetary-mass dependence of gas accretion band width and gas accretion rate…
▽ More
Gas accretion onto the circumplanetary disks and the source region of accreting gas are important to reveal dust accretion that leads to satellite formation around giant planets. We performed local three-dimensional high-resolution hydrodynamic simulations of isothermal and inviscid gas flow around a planet to investigate planetary-mass dependence of gas accretion band width and gas accretion rate onto circumplanetary disks. We examined cases with various planetary masses corresponding to M_p=0.05-1M_{Jup} at 5.2 au, where M_{Jup} is the current Jovian mass. We found that the radial width of the gas accretion band is proportional to M_p^{1/6} for the low-mass regime with M_p < 0.2 M_{Jup} while it is proportional to M_p for the high-mass regime with M_p > 0.2M_{Jup}. We found that the ratio of the mass accretion rate onto the circumplanetary disk to that into the Hill sphere is about 0.4 regardless of planetary mass for the cases we examined. Combining our results with the gap model obtained from global hydrodynamic simulations, we derive semi-analytical formulae of mass accretion rate onto circumplanetary disks. We found that the mass dependence of our three-dimensional accretion rates is the same as the previously-obtained two-dimensional case, although the qualitative behavior of accretion flow onto the CPD is quite different between the two cases.
△ Less
Submitted 7 July, 2022;
originally announced July 2022.
-
Dust coagulation and fragmentation in a collapsing cloud core and their influence on non-ideal magnetohydrodynamic effects
Authors:
Yoshihiro Kawasaki,
Shunta Koga,
Masahiro N. Machida
Abstract:
We determine the time evolution of the dust particle size distribution during the collapse of a cloud core, accounting for both dust coagulation and dust fragmentation, to investigate the influence of dust growth on non-ideal magnetohydrodynamic effects.The density evolution of the collapsing core is given by a one-zone model. We assume two types of dust model: dust composed only of silicate (sili…
▽ More
We determine the time evolution of the dust particle size distribution during the collapse of a cloud core, accounting for both dust coagulation and dust fragmentation, to investigate the influence of dust growth on non-ideal magnetohydrodynamic effects.The density evolution of the collapsing core is given by a one-zone model. We assume two types of dust model: dust composed only of silicate (silicate dust) and dust with a surface covered by $\mathrm{H_{2}O}$ ice ($\mathrm{H_{2}O}$ ice dust). When only considering collisional coagulation, the non-ideal magnetohydrodynamic effects are not effective in the high-density region for both the silicate and $\mathrm{H_{2}O}$ ice dust cases. This is because dust coagulation reduces the abundance of small dust particles, resulting in less efficient adsorption of charged particles on the dust surface. For the silicate dust case, when collisional fragmentation is included, the non-ideal magnetohydrodynamic effects do apply at a high density of $n_{\mathrm{H}}>10^{12} \ \mathrm{cm^{-3}}$ because of the abundant production of small dust particles. On the other hand, for the $\mathrm{H_{2}O}$ ice dust case, the production of small dust particles due to fragmentation is not efficient. Therefore, for the $\mathrm{H_{2}O}$ ice dust case, non-ideal magnetohydrodynamic effects apply only in the range $n_{\mathrm{H}}\gtrsim 10^{14} \ \mathrm{cm^{-3}}$, even when collisional fragmentation is considered. Our results suggest that it is necessary to consider both dust collisional coagulation and fragmentation to activate non-ideal magnetohydrodynamic effects, which should play a significant role in the star and disk formation processes.
△ Less
Submitted 6 July, 2022;
originally announced July 2022.
-
Synthetic Polarization Maps of an Outflow Zone from Magnetohydrodynamic Simulations
Authors:
Gianfranco Bino,
Shantanu Basu,
Masahiro N Machida,
Aris Tritsis,
Mahmoud Sharkawi,
Kundan Kadam,
Indrani Das
Abstract:
The canonical theory of star formation in a magnetized environment predicts the formation of hourglass-shaped magnetic fields during the prestellar collapse phase. In protostellar cores, recent observations reveal complex and strongly distorted magnetic fields in the inner regions that are sculpted by rotation and outflows. We conduct resistive, nonideal magnetohydrodynamic (MHD) simulations of a…
▽ More
The canonical theory of star formation in a magnetized environment predicts the formation of hourglass-shaped magnetic fields during the prestellar collapse phase. In protostellar cores, recent observations reveal complex and strongly distorted magnetic fields in the inner regions that are sculpted by rotation and outflows. We conduct resistive, nonideal magnetohydrodynamic (MHD) simulations of a protostellar core and employ the radiative transfer code POLARIS to produce synthetic polarization segment maps. Comparison of our mock-polarization maps based on the toroidal-dominated magnetic field in the outflow zone with the observed polarization vectors of SiO lines in Orion Source I shows a reasonable agreement when the magnetic axis is tilted at an angle $θ= 15^{\circ}$ with respect to the plane-of-sky and if the SiO lines have a net polarization parallel to the local magnetic field. Although the observed polarization is from SiO lines and our synthetic maps are due to polarized dust emission, a comparison is useful and allows us to resolve the ambiguity of whether the line polarization is parallel or perpendicular to the local magnetic field direction.
△ Less
Submitted 22 August, 2022; v1 submitted 4 July, 2022;
originally announced July 2022.
-
Elemental abundances of nearby M dwarfs based on high-resolution near-infrared spectra obtained by the Subaru/IRD survey: Proof of concept
Authors:
Hiroyuki Tako Ishikawa,
Wako Aoki,
Teruyuki Hirano,
Takayuki Kotani,
Masayuki Kuzuhara,
Masashi Omiya,
Yasunori Hori,
Eiichiro Kokubo,
Tomoyuki Kudo,
Takashi Kurokawa,
Nobuhiko Kusakabe,
Norio Narita,
Jun Nishikawa,
Masahiro Ogihara,
Akitoshi Ueda,
Thayne Currie,
Thomas Henning,
Yui Kasagi,
Jared R. Kolecki,
Jungmi Kwon,
Masahiro N. Machida,
Michael W. McElwain,
Takao Nakagawa,
Sebastien Vievard,
Ji Wang
, et al. (2 additional authors not shown)
Abstract:
Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolu…
▽ More
Detailed chemical analyses of M dwarfs are scarce but necessary to constrain the formation environment and internal structure of planets being found around them. We present elemental abundances of 13 M dwarfs (2900 < Teff < 3500 K) observed in the Subaru/IRD planet search project. They are mid-to-late M dwarfs whose abundance of individual elements has not been well studied. We use the high-resolution (~70,000) near-infrared (970-1750 nm) spectra to measure the abundances of Na, Mg, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Sr by the line-by-line analysis based on model atmospheres, with typical errors ranging from 0.2 dex for [Fe/H] to 0.3-0.4 dex for other [X/H]. We measure radial velocities from the spectra and combine them with Gaia astrometry to calculate the Galactocentric space velocities UVW. The resulting [Fe/H] values agree with previous estimates based on medium-resolution K-band spectroscopy, showing a wide distribution of metallicity (-0.6 < [Fe/H] < +0.4). The abundance ratios of individual elements [X/Fe] are generally aligned with the solar values in all targets. While the [X/Fe] distributions are comparable to those of nearby FGK stars, most of which belong to the thin disk population, the most metal-poor object, GJ 699, could be a thick disk star. The UVW velocities also support this. The results raise the prospect that near-infrared spectra of M dwarfs obtained in the planet search projects can be used to grasp the trend of elemental abundances and Galactic stellar population of nearby M dwarfs.
△ Less
Submitted 30 November, 2021;
originally announced December 2021.
-
"Ash-fall" induced by molecular outflow in protostar evolution
Authors:
Yusuke Tsukamoto,
Masahiro N. Machida,
Shu-ichiro Inutsuka
Abstract:
Dust growth and its associated dynamics play key roles in the first phase of planet formation in young stellar objects (YSOs). Observations have detected signs of dust growth in very young protoplanetary disks. Furthermore, signs of planet formation, gaps in the disk at a distance of several 10 astronomical units (AU) from the central protostar are also reported. From a theoretical point of view,…
▽ More
Dust growth and its associated dynamics play key roles in the first phase of planet formation in young stellar objects (YSOs). Observations have detected signs of dust growth in very young protoplanetary disks. Furthermore, signs of planet formation, gaps in the disk at a distance of several 10 astronomical units (AU) from the central protostar are also reported. From a theoretical point of view, however, it is not clear how planet form at the outer region of a disk despite the difficulty due to rapid inward drift of dust so called radial drift barrier. Here, on the basis of three-dimensional magneto-hydrodynamical simulations of disk evolution with the dust growth, we propose a mechanism named "ash-fall" phenomenon induced by powerful molecular outflow driven by magnetic field which may circumvent the radial drift barrier. We found that the large dust which grows to a size of $\sim \cm$ in the inner region of a disk is entrained by an outflow from the disk. Then large dust decoupled from gas is ejected from the outflow due to centrifugal force, enriching the grown dust in the envelope and is eventually fall onto the outer edge of the disk. The overall process is similar to behaviour of ash-fall from volcanic eruptions. In the ash-fall phenomenon, the Stokes number of dust increases by reaccreting to the less dense disk outer edge. This may make the dust grains overcome the radial drift barrier. Consequently, the ash-fall phenomenon can provide a crucial assist for making the formation of the planetesimals in outer region of the disk possible, and hence the formation of wide-orbit planets and the formation of the gaps.
△ Less
Submitted 19 October, 2021; v1 submitted 26 September, 2021;
originally announced September 2021.
-
Impact of Magnetic Braking on High-mass Close Binary Formation
Authors:
Naoto Harada,
Shingo Hirano,
Masahiro N. Machida,
Takashi Hosokawa
Abstract:
Combining numerical simulations and analytical modeling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the…
▽ More
Combining numerical simulations and analytical modeling, we investigate whether close binary systems form by the effect of magnetic braking. Using magnetohydrodynamics simulations, we calculate the cloud evolution with a sink, for which we do not resolve the binary system or binary orbital motion to realize long-term time integration. Then, we analytically estimate the binary separation using the accreted mass and angular momentum obtained from the simulation. In unmagnetized clouds, wide binary systems with separations of >100 au form, in which the binary separation continues to increase during the main accretion phase. In contrast, close binary systems with separations of <100 au can form in magnetized clouds. Since the efficiency of magnetic braking strongly depends on both the strength and configuration of the magnetic field, they also affect the formation conditions of a close binary. In addition, the protostellar outflow has a negative impact on close binary formation, especially when the rotation axis of the prestellar cloud is aligned with the global magnetic field. The outflow interrupts the accretion of gas with small angular momentum, which is expelled from the cloud, while gas with large angular momentum preferentially falls from the side of the outflow onto the binary system and widens the binary separation. This study shows that a cloud with a magnetic field that is not parallel to the rotation axis is a favorable environment for the formation of close binary systems.
△ Less
Submitted 24 September, 2021;
originally announced September 2021.
-
Can High-velocity Protostellar Jets Help to Drive Low-velocity Outflow?
Authors:
Masahiro N. Machida
Abstract:
Using three-dimensional magnetohydrodynamics simulations, the driving of protostellar jets is investigated in different star-forming cores with the parameters of magnetic field strength and mass accretion rate. Powerful high-velocity jets appear in strongly magnetized clouds when the mass accretion rate onto the protostellar system is lower than $\dot{M} \lesssim 10^{-3}\,{\rm M}_\odot$ yr$^{-1}$.…
▽ More
Using three-dimensional magnetohydrodynamics simulations, the driving of protostellar jets is investigated in different star-forming cores with the parameters of magnetic field strength and mass accretion rate. Powerful high-velocity jets appear in strongly magnetized clouds when the mass accretion rate onto the protostellar system is lower than $\dot{M} \lesssim 10^{-3}\,{\rm M}_\odot$ yr$^{-1}$. On the other hand, even at this mass accretion rate range, no jets appear for magnetic fields of prestellar clouds as weak as $μ_0 \gtrsim 5$--$10$, where $μ_0$ is the mass-to-flux ratio normalized by the critical value $(2πG^{1/2})^{-1}$. For $\dot{M}\gtrsim 10^{-3}\,{\rm M}_\odot$ yr$^{-1}$, although jets usually appear just after protostar formation independent of the magnetic field strength, they soon weaken and finally disappear. Thus, they cannot help drive the low-velocity outflow when there is no low-velocity flow just before protostar formation. As a result, no significant mass ejection occurs during the early mass accretion phase either when the prestellar cloud is weaky magnetized or when the mass accretion rate is very high. Thus, protostars formed in such environments would trace different evolutionary paths from the normal star formation process.
△ Less
Submitted 13 September, 2021;
originally announced September 2021.
-
Super-fast Rotation in the OMC 2/FIR 6b Jet
Authors:
Yuko Matsushita,
Satoko Takahashi,
Shun Ishii,
Kohji Tomisaka,
Paul T. P. Ho,
John M. Carpenter,
Masahiro N. Machida
Abstract:
We present ALMA CO ($J$=2--1) and 1.3 mm continuum observations of the high-velocity jet associated with the FIR 6b protostar located in the Orion Molecular Cloud-2. We detect a velocity gradient along the short axis of the jet in both the red- and blue-shifted components. The position-velocity diagrams along the short axis of the red-shifted jet show a typical characteristic of a rotating cylinde…
▽ More
We present ALMA CO ($J$=2--1) and 1.3 mm continuum observations of the high-velocity jet associated with the FIR 6b protostar located in the Orion Molecular Cloud-2. We detect a velocity gradient along the short axis of the jet in both the red- and blue-shifted components. The position-velocity diagrams along the short axis of the red-shifted jet show a typical characteristic of a rotating cylinder. We attribute the velocity gradient in the red-shifted component to rotation of the jet. The rotation velocity ($>20\,\ \rm{km s^{-1}}$) and specific angular momentum ($>10^{22}\, \rm{cm^{2}\, s^{-1}}$) of the jet around FIR 6b are the largest among all jets in which rotation has been observed. By combining disk wind theory with our observations, the jet launching radius is estimated to be in the range of $2.18-2.96$\,au. The rapid rotation, large specific angular momentum, and a launching radius far from the central protostar can be explained by a magnetohydrodynamic disk wind that contributes to the angular momentum transfer in the late stages of protostellar accretion.
△ Less
Submitted 26 July, 2021;
originally announced July 2021.
-
Jets from MRC 0600-399 bent by magnetic fields in the cluster Abell 3376
Authors:
James O. Chibueze,
Haruka Sakemi,
Takumi Ohmura,
Mami Machida,
Hiroki Akamatsu,
Takuya Akahori,
Hiroyuki Nakanishi,
Viral Parekh,
Ruby van Rooyen,
Tsutomu T. Takeuchi
Abstract:
Galaxy clusters are known to harbour magnetic fields. The nature of the intra-cluster magnetic fields remains an unresolved question. Intra-cluster magnetic field can be observed at the density contact discontinuity formed by cool and dense plasma running into hot ambient plasma, and the discontinuity exists near the 2nd BCG MRC 0600-399 of a merging galaxy cluster Abell 3376 (z=0.0461, hereafter…
▽ More
Galaxy clusters are known to harbour magnetic fields. The nature of the intra-cluster magnetic fields remains an unresolved question. Intra-cluster magnetic field can be observed at the density contact discontinuity formed by cool and dense plasma running into hot ambient plasma, and the discontinuity exists near the 2nd BCG MRC 0600-399 of a merging galaxy cluster Abell 3376 (z=0.0461, hereafter as A3376). Elongated X-ray image in the east-west direction with a comet-like structure reaches a Mpc-scale (Fig1.(a)). Previous radio observations detected the bent jets from MRC 0600-399, moving in same direction as the sub-cluster's motion against ram pressure.Here we report a new radio observation of a radio galaxy MRC 0600-399 which is 3.4 times and 11 times higher resolution and sensitivity than the previous results. Contrary to typical jets, the MRC 0600-399 shows a 90deg bend at the contact discontinuity and the collimated jets further extend over 100 kpc from the bend point. Diffuse, elongated emission named "double-scythe" structures were detected for the first time. The spectral index flattens downstream of the bend point, indicating cosmic-ray re-acceleration. High-resolution numerical simulations reveal that the ordered magnetic field along the discontinuity plays a significant role in the change in the jet direction. The morphology of the "double-scythe" bear remarkable similarities with the simulations, which strengthens our understanding of the interaction between relativistic electrons and intra-cluster magnetic field.
△ Less
Submitted 24 June, 2021;
originally announced June 2021.
-
Conditions for justifying single-fluid approximation for charged and neutral dust fluids and a smoothed particle magnetohydrodynamics method for dust-gas mixture
Authors:
Y. Tsukamoto,
M. N. Machida,
S. Inutsuka
Abstract:
We describe a numerical scheme for magnetohydrodynamics simulations of dust-gas mixture by extending smoothed particle magnetohydrodynamics. We employ the single-species particle approach to describe dust-gas mixture with several modifications from the previous studies. We assume that the charged and neutral dusts can be treated as single-fluid and the electro-magnetic force acts on the gas and th…
▽ More
We describe a numerical scheme for magnetohydrodynamics simulations of dust-gas mixture by extending smoothed particle magnetohydrodynamics. We employ the single-species particle approach to describe dust-gas mixture with several modifications from the previous studies. We assume that the charged and neutral dusts can be treated as single-fluid and the electro-magnetic force acts on the gas and that on the charged dust is negligible. The validity of these assumption in the context of protostar formation is not obvious and is extensively evaluated. By investigating the electromagnetic force and electric current with terminal velocity approximation, it is found that as the dust size increases, the contribution of dust to them becomes smaller and negligible. We conclude that our assumptions of the electro-magnetic force on the dusts is negligible are valid for the dust size with a d & 10μm. On the other hand, they do not produce the numerical artifact for the dust a d . 10μm in envelope and disk where the perfect coupling between gas and dusts realizes. However, we also found that our assumptions may break down in outflow (or under environment with very strong magnetic field and low density) for the dust a d . 10μm. We conclude that our assumptions are valid in almost all cases where macroscopic dust dynamics is important in the context of protostar formation. We conduct numerical tests of dusty wave, dusty magnetohydrodynamics shock, and gravitational collapse of magnetized cloud core with our simulation code. The results show that our numerical scheme well reproduces the dust dynamics in the magnetized medium.
△ Less
Submitted 16 June, 2021;
originally announced June 2021.
-
Supermassive Star Formation in Magnetized Atomic-Cooling Gas Clouds: Enhanced Accretion, Intermittent Fragmentation, and Continuous Mergers
Authors:
Shingo Hirano,
Masahiro N. Machida,
Shantanu Basu
Abstract:
The origin of supermassive black holes (with $\gtrsim\!10^9\,M_{\odot}$) in the early universe (redshift $z \sim 7$) remains poorly understood. Gravitational collapse of a massive primordial gas cloud is a promising initial process, but theoretical studies have difficulty growing the black hole fast enough. We focus on the magnetic effects on star formation that occurs in an atomic-cooling gas clo…
▽ More
The origin of supermassive black holes (with $\gtrsim\!10^9\,M_{\odot}$) in the early universe (redshift $z \sim 7$) remains poorly understood. Gravitational collapse of a massive primordial gas cloud is a promising initial process, but theoretical studies have difficulty growing the black hole fast enough. We focus on the magnetic effects on star formation that occurs in an atomic-cooling gas cloud. Using a set of three-dimensional magnetohydrodynamic (MHD) simulations, we investigate the star formation process in the magnetized atomic-cooling gas cloud with different initial magnetic field strengths. Our simulations show that the primordial magnetic seed field can be quickly amplified during the early accretion phase after the first protostar formation. The strong magnetic field efficiently extracts angular momentum from accreting gas and increases the accretion rate, which results in the high fragmentation rate in the gravitationally unstable disk region. On the other hand, the coalescence rate of fragments is also enhanced by the angular momentum transfer due to the magnetic effects. Almost all the fragments coalesce to the primary star, so the mass growth rate of the massive star increases due to the magnetic effects. We conclude that the magnetic effects support the direct collapse scenario of supermassive star formation.
△ Less
Submitted 7 June, 2021;
originally announced June 2021.
-
ALMA Observations toward the S-shaped Outflow and the Envelope around NGC1333 IRAS 4A2
Authors:
Chen-Yu Chuang,
Yusuke Aso,
Naomi Hirano,
Shingo Hirano,
Masahiro N. Machida
Abstract:
We have analyzed the ALMA archival data of the SO ($J_N=6_5-5_4$ and $J_N=7_6-6_5$), CO ($J=2-1$), and CCH ($N=3-2, J=7/2-5/2, F=4-3$) lines from the class 0 protobinary system, NGC1333 IRAS 4A. The images of SO ($J_N = 6_5-5_4$) and CO ($J=2-1$) successfully separate two northern outflow lobes connected to each protostar, IRAS 4A1 and IRAS 4A2. The outflow from IRAS 4A2 shows an S-shaped morpholo…
▽ More
We have analyzed the ALMA archival data of the SO ($J_N=6_5-5_4$ and $J_N=7_6-6_5$), CO ($J=2-1$), and CCH ($N=3-2, J=7/2-5/2, F=4-3$) lines from the class 0 protobinary system, NGC1333 IRAS 4A. The images of SO ($J_N = 6_5-5_4$) and CO ($J=2-1$) successfully separate two northern outflow lobes connected to each protostar, IRAS 4A1 and IRAS 4A2. The outflow from IRAS 4A2 shows an S-shaped morphology, consisting of a flattened envelope around IRAS 4A2 with two outflow lobes connected to both edges of the envelope. The flattened envelope surrounding IRAS 4A2 has an opposite velocity gradient to that of the circumbinary envelope. The observed features are reproduced by the magnetohydrodynamic simulation of the collapsing core whose magnetic field direction is misaligned to the rotational axis. Our simulation shows that the intensity of the outflow lobes is enhanced on one side, resulting in the formation of S-shaped morphology. The S-shaped outflow can also be explained by the precessing outflow launched from an unresolved binary with a separation larger than 12 au (0.04arcsec). Additionally, we discovered a previously unknown extremely high velocity component at $\sim$45-90 km/s near IRAS 4A2 with CO. CCH ($J_{N,F}=7/2_{3,4}-5/2_{2,3}$) emission shows two pairs of blobs attaching to the bottom of shell like feature, and the morphology is significantly different from those of SO and CO lines. Toward IRAS 4A2, the S-shaped outflow shown in SO is overlapped with the edges of CCH shells, while CCH shells have the velocity gradients opposite to the flattened structure around IRAS 4A2.
△ Less
Submitted 10 May, 2021;
originally announced May 2021.
-
Growth of Magnetorotational Instability in Circumstellar Disks around Class 0 Protostars
Authors:
Yoshihiro Kawasaki,
Shunta Koga,
Masahiro N. Machida
Abstract:
We investigate the possibility of the growth of magnetorotational instability (MRI) in disks around Class 0 protostars. We construct a disk model and calculate the chemical reactions of neutral and charged atoms, molecules and dust grains to derive the abundance of each species and the ionization degree of the disk. Then, we estimate the diffusion coefficients of non-ideal magnetohydrodynamics eff…
▽ More
We investigate the possibility of the growth of magnetorotational instability (MRI) in disks around Class 0 protostars. We construct a disk model and calculate the chemical reactions of neutral and charged atoms, molecules and dust grains to derive the abundance of each species and the ionization degree of the disk. Then, we estimate the diffusion coefficients of non-ideal magnetohydrodynamics effects such as ohmic dissipation, ambipolar diffusion and the Hall effect. Finally, we evaluate the linear growth rate of MRI in each area of the disk. We investigate the effect of changes in the strength and direction of the magnetic field in our disk model and we adopt four different dust models to investigate the effect of dust size distribution on the diffusion coefficients. Our results indicate that an MRI active region possibly exists with a weak magnetic field in a region far from the protostar where the Hall effect plays a role in the growth of MRI. On the other hand, in all models the disk is stable against MRI in the region within $<20$ au from the protostar on the equatorial plane. Since the size of the disks in the early stage of star formation is limited to $\lesssim 10-$$20$ au, it is difficult to develop MRI-driven turbulence in such disks.
△ Less
Submitted 27 April, 2021;
originally announced April 2021.
-
Gravitational Wave Physics and Astronomy in the nascent era
Authors:
Makoto Arimoto,
Hideki Asada,
Michael L. Cherry,
Michiko S. Fujii,
Yasushi Fukazawa,
Akira Harada,
Kazuhiro Hayama,
Takashi Hosokawa,
Kunihito Ioka,
Yoichi Itoh,
Nobuyuki Kanda,
Koji S. Kawabata,
Kyohei Kawaguchi,
Nobuyuki Kawai,
Tsutomu Kobayashi,
Kazunori Kohri,
Yusuke Koshio,
Kei Kotake,
Jun Kumamoto,
Masahiro N. Machida,
Hideo Matsufuru,
Tatehiro Mihara,
Masaki Mori,
Tomoki Morokuma,
Shinji Mukohyama
, et al. (28 additional authors not shown)
Abstract:
The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion,…
▽ More
The detections of gravitational waves (GW) by LIGO/Virgo collaborations provide various possibilities to physics and astronomy. We are quite sure that GW observations will develop a lot both in precision and in number owing to the continuous works for the improvement of detectors, including the expectation to the newly joined detector, KAGRA, and the planned detector, LIGO-India. In this occasion, we review the fundamental outcomes and prospects of gravitational wave physics and astronomy. We survey the development focusing on representative sources of gravitational waves: binary black holes, binary neutron stars, and supernovae. We also summarize the role of gravitational wave observations as a probe of new physics.
△ Less
Submitted 6 April, 2021;
originally announced April 2021.