-
Gamma-ray Bursts as Distance Indicators by a Statistical Learning Approach
Authors:
Maria Giovanna Dainotti,
Aditya Narendra,
Agnieszka Pollo,
Vahe Petrosian,
Malgorzata Bogdan,
Kazunari Iwasaki,
Jason Xavier Prochaska,
Enrico Rinaldi,
David Zhou
Abstract:
Gamma-ray bursts (GRBs) can be probes of the early universe, but currently, only 26% of GRBs observed by the Neil Gehrels Swift Observatory GRBs have known redshifts ($z$) due to observational limitations. To address this, we estimated the GRB redshift (distance) via a supervised statistical learning model that uses optical afterglow observed by Swift and ground-based telescopes. The inferred reds…
▽ More
Gamma-ray bursts (GRBs) can be probes of the early universe, but currently, only 26% of GRBs observed by the Neil Gehrels Swift Observatory GRBs have known redshifts ($z$) due to observational limitations. To address this, we estimated the GRB redshift (distance) via a supervised statistical learning model that uses optical afterglow observed by Swift and ground-based telescopes. The inferred redshifts are strongly correlated (a Pearson coefficient of 0.93) with the observed redshifts, thus proving the reliability of this method. The inferred and observed redshifts allow us to estimate the number of GRBs occurring at a given redshift (GRB rate) to be 8.47-9 $yr^{-1} Gpc^{-1}$ for $1.9<z<2.3$. Since GRBs come from the collapse of massive stars, we compared this rate with the star formation rate highlighting a discrepancy of a factor of 3 at $z<1$.
△ Less
Submitted 2 May, 2024; v1 submitted 6 February, 2024;
originally announced February 2024.
-
Filamentary Network and Magnetic Field Structures Revealed with BISTRO in the High-Mass Star-Forming Region NGC2264 : Global Properties and Local Magnetogravitational Configurations
Authors:
Jia-Wei Wang,
Patrick M. Koch,
Seamus D. Clarke,
Gary Fuller,
Nicolas Peretto,
Ya-Wen Tang,
Hsi-Wei Yen,
Shih-Ping Lai,
Nagayoshi Ohashi,
Doris Arzoumanian,
Doug Johnstone,
Ray Furuya,
Shu-ichiro Inutsuka,
Chang Won Lee,
Derek Ward-Thompson,
Valentin J. M. Le Gouellec,
Hong-Li Liu,
Lapo Fanciullo,
Jihye Hwang,
Kate Pattle,
Frédérick Poidevin,
Mehrnoosh Tahani,
Takashi Onaka,
Mark G. Rawlings,
Eun Jung Chung
, et al. (132 additional authors not shown)
Abstract:
We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from…
▽ More
We report 850 $μ$m continuum polarization observations toward the filamentary high-mass star-forming region NGC 2264, taken as part of the B-fields In STar forming Regions Observations (BISTRO) large program on the James Clerk Maxwell Telescope (JCMT). These data reveal a well-structured non-uniform magnetic field in the NGC 2264C and 2264D regions with a prevailing orientation around 30 deg from north to east. Field strengths estimates and a virial analysis for the major clumps indicate that NGC 2264C is globally dominated by gravity while in 2264D magnetic, gravitational, and kinetic energies are roughly balanced. We present an analysis scheme that utilizes the locally resolved magnetic field structures, together with the locally measured gravitational vector field and the extracted filamentary network. From this, we infer statistical trends showing that this network consists of two main groups of filaments oriented approximately perpendicular to one another. Additionally, gravity shows one dominating converging direction that is roughly perpendicular to one of the filament orientations, which is suggestive of mass accretion along this direction. Beyond these statistical trends, we identify two types of filaments. The type-I filament is perpendicular to the magnetic field with local gravity transitioning from parallel to perpendicular to the magnetic field from the outside to the filament ridge. The type-II filament is parallel to the magnetic field and local gravity. We interpret these two types of filaments as originating from the competition between radial collapsing, driven by filament self-gravity, and the longitudinal collapsing, driven by the region's global gravity.
△ Less
Submitted 23 January, 2024;
originally announced January 2024.
-
Dynamics Near the Inner Dead-Zone Edges in a Proprotoplanetary Disk
Authors:
Kazunari Iwasaki,
Kengo Tomida,
Shinsuke Takasao,
Satoshi Okuzumi,
Takeru K. Suzuki
Abstract:
We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Elsässer numbers of the Ohmic r…
▽ More
We perform three-dimensional global non-ideal magnetohydrodynamic simulations of a protoplanetary disk containing the inner dead-zone edge. We take into account realistic diffusion coefficients of the Ohmic resistivity and ambipolar diffusion based on detailed chemical reactions with single-size dust grains. We found that the conventional dead zone identified by the Elsässer numbers of the Ohmic resistivity and ambipolar diffusion is divided into two regions: "the transition zone" and "the coherent zone". The coherent zone has the same properties as the conventional dead zone, and extends outside of the transition zone in the radial direction. Between the active and coherent zones, we discover the transition zone, the inner edge of which is identical to that of the conventional dead zone. The transition zone extends out over the regions where thermal ionization determines diffusion coefficients. The transition zone has completely different physical properties than the conventional dead zone, the so-called undead zone, and the zombie zone. The combination of amplification of the radial magnetic field owing to the ambipolar diffusion and a steep radial gradient of the Ohmic diffusivity causes the efficient evacuation of the net vertical magnetic flux from the transition zone within several rotations. Surface gas accretion occurs in the coherent zone but not in the transition zone. The presence of the transition zone prohibits mass and magnetic flux transport from the coherent zone to the active zone. Mass accumulation occurs at both edges of the transition zone as a result of mass supply from the active and coherent zones.
△ Less
Submitted 8 June, 2024; v1 submitted 8 January, 2024;
originally announced January 2024.
-
A Study of the Long-Term Behavior of Hybrid Systems with Symmetries via Reduction and the Frobenius-Perron Operator
Authors:
Maria Oprea,
Aden Shaw,
Robi Huq,
Kaito Iwasaki,
Dora Kassabova,
W. Clark
Abstract:
Hybrid dynamical systems are systems which undergo both continuous and discrete transitions. As typical in dynamical analysis, an essential goal is to study the long-term behavior of these systems. In this work, we present two different novel approaches for studying these systems. The first approach is based on constructing an analog of the Frobenius-Perron (transport) operator for hybrid systems.…
▽ More
Hybrid dynamical systems are systems which undergo both continuous and discrete transitions. As typical in dynamical analysis, an essential goal is to study the long-term behavior of these systems. In this work, we present two different novel approaches for studying these systems. The first approach is based on constructing an analog of the Frobenius-Perron (transport) operator for hybrid systems. Rather than tracking the evolution of a single trajectory, this operator encodes the asymptotic nature of an ensemble of trajectories. The second approach presented applies to an important subclass of hybrid systems, mechanical impact systems. We develop an analog of Lie-Poisson(-Suslov) reduction for left-invariant impact systems on Lie groups. In addition to the Hamiltonian (and constraints) being left-invariant, the impact surface must also be a right coset of a normal subgroup. This procedure allows a reduction from a $2n$-dimensional system to an $(n+1)$-dimensional one. We conclude the paper by presenting numerical results on a diverse array of applications.
△ Less
Submitted 21 September, 2023;
originally announced September 2023.
-
Lorentz- and CPT-violating effects in Penning traps at linear boost
Authors:
Ariam J. Acevedo-Lopez,
Yunhua Ding,
Kaito Iwasaki,
Arnaldo J. Vargas
Abstract:
We present in this work an analysis of Lorentz- and CPT-violating signals at linear boost order in Penning-trap experiments. The theory of quantum electrodynamics with Penning traps is revisited and the dominant shifts in the cyclotron and anomaly frequencies of confined particles and antiparticles are reproduced. To study time variations of the experimental signals at linear boost order, we provi…
▽ More
We present in this work an analysis of Lorentz- and CPT-violating signals at linear boost order in Penning-trap experiments. The theory of quantum electrodynamics with Penning traps is revisited and the dominant shifts in the cyclotron and anomaly frequencies of confined particles and antiparticles are reproduced. To study time variations of the experimental signals at linear boost order, we provide a general discussion on transformations of coefficients for Lorentz violation between different frames, and derive the expressions of the cyclotron and anomaly frequency shifts in the Sun-centered frame. Relating these frequency shifts to the charge-to-mass ratios, the $g$ factors, and their comparisons between particles and antiparticles, we extract numerous new or improved bounds on coefficients for Lorentz violation from existing Penning-trap measurements.
△ Less
Submitted 27 August, 2023;
originally announced August 2023.
-
Simulation of Head-on Collisions Between Filamentary Molecular Clouds Threaded by a Lateral Magnetic Field and Subsequent Evolution
Authors:
Raiga Kashiwagi,
Kazunari Iwasaki,
Kohji Tomisaka
Abstract:
Filamentary molecular clouds are regarded as the place where newborn stars are formed. In particular, a hub region, a place where it appears as if several filaments are colliding, often indicates active star formation. To understand the star formation in filament structures, we investigate the collisions between two filaments using two-dimensional magnetohydrodynamical simulations. As a model of f…
▽ More
Filamentary molecular clouds are regarded as the place where newborn stars are formed. In particular, a hub region, a place where it appears as if several filaments are colliding, often indicates active star formation. To understand the star formation in filament structures, we investigate the collisions between two filaments using two-dimensional magnetohydrodynamical simulations. As a model of filaments, we assume that the filaments are in magnetohydrostatic equilibrium under a global magnetic field perpendicular to the filament axis. We set two identical filaments with an infinite length and collided them with a zero-impact parameter (head-on). When the two filaments collide while sharing the same magnetic flux, we found two types of evolution after a merged filament is formed: runaway radial collapse and stable oscillation with a finite amplitude. The condition for the radial collapse is independent of the collision velocity and is given by the total line mass of the two filaments exceeding the magnetically critical line mass for which no magnetohydrostatic solution exists. The radial collapse proceeds in a self-similar manner, resulting in a unique distribution irrespective of the various initial line masses of the filament, as the collapse progresses. When the total line mass is less massive than the magnetically critical line mass, the merged filament oscillates, and the density distribution is well-fitted by a magnetohydrostatic equilibrium solution. The condition necessary for the radial collapse is also applicable to the collision whose direction is perpendicular to the global magnetic field.
△ Less
Submitted 14 July, 2023;
originally announced July 2023.
-
Metallicity Dependence of Molecular Cloud Hierarchical Structure at Early Evolutionary Stages
Authors:
Masato I. N. Kobayashi,
Kazunari Iwasaki,
Kengo Tomida,
Tsuyoshi Inoue,
Kazuyuki Omukai,
Kazuki Tokuda
Abstract:
The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however,…
▽ More
The formation of molecular clouds out of HI gas is the first step toward star formation. Its metallicity dependence plays a key role to determine star formation through the cosmic history. Previous theoretical studies with detailed chemical networks calculate thermal equilibrium states and/or thermal evolution under one-zone collapsing background. The molecular cloud formation in reality, however, involves supersonic flows, and thus resolving the cloud internal turbulence/density structure in three dimension is still essential. We here perform magnetohydrodynamics simulations of 20 km s^-1 converging flows of Warm Neutral Medium (WNM) with 1 micro Gauss mean magnetic field in the metallicity range from the Solar (1.0 Zsun) to 0.2 Zsun environment. The Cold Neutral Medium (CNM) clumps form faster with higher metallicity due to more efficient cooling. Meanwhile, their mass functions commonly follow dn/dm proportional to m^-1.7 at three cooling times regardless of the metallicity. Their total turbulence power also commonly shows the Kolmogorov spectrum with its 80 percent in the solenoidal mode, while the CNM volume alone indicates the transition towards the Larson's law. These similarities measured at the same time in the unit of the cooling time suggest that the molecular cloud formation directly from the WNM alone requires a longer physical time in a lower metallicity environment in the 1.0--0.2 Zsun range. To explain the rapid formation of molecular clouds and subsequent massive star formation possibly within less than 10 Myr as observed in the Large/Small Magellanic Clouds (LMC/SMC), the HI gas already contains CNM volume instead of pure WNM.
△ Less
Submitted 23 July, 2023; v1 submitted 3 July, 2023;
originally announced July 2023.
-
Primordial or Secondary? Testing models of debris disk gas with ALMA
Authors:
Gianni Cataldi,
Yuri Aikawa,
Kazunari Iwasaki,
Sebastian Marino,
Alexis Brandeker,
Antonio Hales,
Thomas Henning,
Aya E. Higuchi,
A. Meredith Hughes,
Markus Janson,
Quentin Kral,
Luca Matrà,
Attila Moór,
Göran Olofsson,
Seth Redfield,
Aki Roberge
Abstract:
The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA obs…
▽ More
The origin and evolution of gas in debris disks is still not well understood. Secondary gas production from cometary material or a primordial origin have been proposed. So far, observations have mostly concentrated on CO, with only few C observations available. We create an overview of the C and CO content of debris disk gas and use it test state-of-the-art models. We use new and archival ALMA observations of CO and CI emission, complemented by CII data from Herschel, for a sample of 14 debris disks. This expands the number of disks with ALMA measurements of both CO and CI by ten disks. We present new detections of CI emission towards three disks: HD 21997, HD 121191 and HD 121617. We use a simple disk model to derive gas masses and column densities. We find that current state-of-the-art models of secondary gas production overpredict the neutral carbon content of debris disk gas. This does not rule out a secondary origin, but might indicate that the models require an additional C removal process. Alternatively, the gas might be produced in transient events rather than a steady-state collisional cascade. We also test a primordial gas origin by comparing our results to a simplified thermo-chemical model. This yields promising results, but more detailed work is required before a conclusion can be reached. Our work demonstrates that the combination of C and CO data is a powerful tool to advance our understanding of debris disk gas.
△ Less
Submitted 5 June, 2023; v1 submitted 20 May, 2023;
originally announced May 2023.
-
The JCMT BISTRO Survey: Studying the Complex Magnetic Field of L43
Authors:
Janik Karoly,
Derek Ward-Thompson,
Kate Pattle,
David Berry,
Anthony Whitworth,
Jason Kirk,
Pierre Bastien,
Tao-Chung Ching,
Simon Coude,
Jihye Hwang,
Woojin Kwon,
Archana Soam,
Jia-Wei Wang,
Tetsuo Hasegawa,
Shih-Ping Lai,
Keping Qiu,
Doris Arzoumanian,
Tyler L. Bourke,
Do-Young Byun,
Huei-Ru Vivien Chen,
Wen Ping Chen,
Mike Chen,
Zhiwei Chen,
Jungyeon Cho,
Minho Choi
, et al. (133 additional authors not shown)
Abstract:
We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar…
▽ More
We present observations of polarized dust emission at 850 $μ$m from the L43 molecular cloud which sits in the Ophiuchus cloud complex. The data were taken using SCUBA-2/POL-2 on the James Clerk Maxwell Telescope as a part of the BISTRO large program. L43 is a dense ($N_{\rm H_2}\sim 10^{22}$-10$^{23}$ cm$^{-2}$) complex molecular cloud with a submillimetre-bright starless core and two protostellar sources. There appears to be an evolutionary gradient along the isolated filament that L43 is embedded within, with the most evolved source closest to the Sco OB2 association. One of the protostars drives a CO outflow that has created a cavity to the southeast. We see a magnetic field that appears to be aligned with the cavity walls of the outflow, suggesting interaction with the outflow. We also find a magnetic field strength of up to $\sim$160$\pm$30 $μ$G in the main starless core and up to $\sim$90$\pm$40 $μ$G in the more diffuse, extended region. These field strengths give magnetically super- and sub-critical values respectively and both are found to be roughly trans-Alfvénic. We also present a new method of data reduction for these denser but fainter objects like starless cores.
△ Less
Submitted 22 May, 2023; v1 submitted 18 May, 2023;
originally announced May 2023.
-
Reducing the uncertainty on the Hubble constant up to 35\% with an improved statistical analysis: different best-fit likelihoods for Supernovae Ia, Baryon Acoustic Oscillations, Quasars, and Gamma-Ray Bursts
Authors:
Maria Giovanna Dainotti,
Giada Bargiacchi,
Małgorzata Bogdan,
Aleksander Łukasz Lenart,
Kazunari Iwasaki,
Salvatore Capozziello,
Bing Zhang,
Nissim Fraija
Abstract:
Cosmological models and their parameters are widely debated, especially about whether the current discrepancy between the values of the Hubble constant, $H_{0}$, obtained by type Ia supernovae (SNe Ia), and the Planck data from the Cosmic Microwave Background Radiation could be alleviated when alternative cosmological models are considered. Thus, combining high-redshift probes, such as Gamma-Ray B…
▽ More
Cosmological models and their parameters are widely debated, especially about whether the current discrepancy between the values of the Hubble constant, $H_{0}$, obtained by type Ia supernovae (SNe Ia), and the Planck data from the Cosmic Microwave Background Radiation could be alleviated when alternative cosmological models are considered. Thus, combining high-redshift probes, such as Gamma-Ray Bursts (GRBs) and Quasars (QSOs), together with Baryon Acoustic Oscillations (BAO) and SNe Ia is important to assess the viability of these alternative models and if they can cast further light on the Hubble tension. In this work, for GRBs, we use a 3-dimensional relation between the peak prompt luminosity, the rest-frame time at the end of the X-ray plateau, and its corresponding luminosity in X-rays: the 3D Dainotti fundamental plane relation. Regarding QSOs, we use the Risaliti-Lusso relation among the UV and X-ray luminosities for a sample of 2421 sources. We correct both the QSO and GRB relations by accounting for selection and evolutionary effects with a reliable statistical method. We here use both the traditional Gaussian likelihoods ($\cal L_G$) and the new best-fit likelihoods ($\cal L_N$) to infer cosmological parameters of a non-flat $Λ$CDM and flat $w$CDM models. We obtain for all the parameters reduced uncertainties, up to $35\%$ for $H_{0}$, when applying the new $\cal L_N$ likelihoods in place of the Gaussian ones. Our results remain consistent with a flat $Λ$CDM model, although with a shift of the dark energy parameter $w$ toward $w<-1$ and a curvature density parameter toward $Ω_k<0$.
△ Less
Submitted 17 May, 2023;
originally announced May 2023.
-
Insights on the Sun birth environment in the context of star-cluster formation in hub-filament systems
Authors:
Doris Arzoumanian,
Sota Arakawa,
Masato I. N. Kobayashi,
Kazunari Iwasaki,
Kohei Fukada,
Shoji Mori,
Yutaka Hirai,
Masanobu Kunitomo,
M. S. Nanda Kumar,
Eiichiro Kokubo
Abstract:
Cylindrical molecular filaments are observed to be the main sites of Sun-like star formation, while massive stars form in dense hubs, at the junction of multiple filaments. The role of hub-filament configurations has not been discussed yet in relation to the birth environment of the solar system and to infer the origin of isotopic ratios of Short-Lived Radionuclides (SLR, such as $^{26}$Al) of Cal…
▽ More
Cylindrical molecular filaments are observed to be the main sites of Sun-like star formation, while massive stars form in dense hubs, at the junction of multiple filaments. The role of hub-filament configurations has not been discussed yet in relation to the birth environment of the solar system and to infer the origin of isotopic ratios of Short-Lived Radionuclides (SLR, such as $^{26}$Al) of Calcium-Aluminum-rich Inclusions (CAIs) observed in meteorites. In this work, we present simple analytical estimates of the impact of stellar feedback on the young solar system forming along a filament of a hub-filament system. We find that the host filament can shield the young solar system from the stellar feedback, both during the formation and evolution of stars (stellar outflow, wind, and radiation) and at the end of their life (supernovae). We show that the young solar system formed along a dense filament can be enriched with supernova ejecta (e.g., $^{26}$Al) during the formation timescale of CAIs. We also propose that the streamers recently observed around protostars may be channeling the SLR-rich material onto the young solar system. We conclude that considering hub-filament configurations as the birth environment of the Sun is important when deriving theoretical models explaining the observed properties of the solar system.
△ Less
Submitted 27 March, 2023;
originally announced March 2023.
-
A Constraint on the Amount of Hydrogen from the CO Chemistry in Debris Disks
Authors:
Kazunari Iwasaki,
Hiroshi Kobayashi,
Aya E. Higuchi,
Yuri Aikawa
Abstract:
The faint CO gases in debris disks are easily dissolved into C by UV irradiation, while CO can be reformed via reactions with hydrogen. The abundance ratio of C/CO could thus be a probe of the amount of hydrogen in the debris disks. We conduct radiative transfer calculations with chemical reactions for debris disks. For a typical dust-to-gas mass ratio of debris disks, CO formation proceeds withou…
▽ More
The faint CO gases in debris disks are easily dissolved into C by UV irradiation, while CO can be reformed via reactions with hydrogen. The abundance ratio of C/CO could thus be a probe of the amount of hydrogen in the debris disks. We conduct radiative transfer calculations with chemical reactions for debris disks. For a typical dust-to-gas mass ratio of debris disks, CO formation proceeds without the involvement of H$_2$ because a small amount of dust grains makes H$_2$ formation inefficient. We find that the CO to C number density ratio depends on a combination of $n_\mathrm{H}Z^{0.4}χ^{-1.1}$, where $n_\mathrm{H}$ is the hydrogen nucleus number density, $Z$ is the metallicity, and $χ$ is the FUV flux normalized by the Habing flux. Using an analytic formula for the CO number density, we give constraints on the amount of hydrogen and metallicity for debris disks. CO formation is accelerated by excited H$_2$ either when the dust-to-gas mass ratio is increased or the energy barrier of chemisorption of hydrogen on the dust surface is decreased. This acceleration of CO formation occurs only when the shielding effects of CO are insignificant. In shielded regions, the CO fractions are almost independent of the parameters of dust grains.
△ Less
Submitted 20 March, 2023;
originally announced March 2023.
-
First BISTRO observations of the dark cloud Taurus L1495A-B10: the role of the magnetic field in the earliest stages of low-mass star formation
Authors:
Derek Ward-Thompson,
Janik Karoly,
Kate Pattle,
Anthony Whitworth,
Jason Kirk,
David Berry,
Pierre Bastien,
Tao-Chung Ching,
Simon Coude,
Jihye Hwang,
Woojin Kwon,
Archana Soam,
Jia-Wei Wang,
Tetsuo Hasegawa,
Shih-Ping Lai,
Keping Qiu,
Doris Arzoumanian,
Tyler L. Bourke,
Do-Young Byun,
Huei-Ru Vivien Chen,
Wen Ping Chen,
Mike Chen,
Zhiwei Chen,
Jungyeon Cho,
Minho Choi
, et al. (133 additional authors not shown)
Abstract:
We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to…
▽ More
We present BISTRO Survey 850 μm dust emission polarisation observations of the L1495A-B10 region of the Taurus molecular cloud, taken at the JCMT. We observe a roughly triangular network of dense filaments. We detect 9 of the dense starless cores embedded within these filaments in polarisation, finding that the plane-of-sky orientation of the core-scale magnetic field lies roughly perpendicular to the filaments in almost all cases. We also find that the large-scale magnetic field orientation measured by Planck is not correlated with any of the core or filament structures, except in the case of the lowest-density core. We propose a scenario for early prestellar evolution that is both an extension to, and consistent with, previous models, introducing an additional evolutionary transitional stage between field-dominated and matter-dominated evolution, observed here for the first time. In this scenario, the cloud collapses first to a sheet-like structure. Uniquely, we appear to be seeing this sheet almost face-on. The sheet fragments into filaments, which in turn form cores. However, the material must reach a certain critical density before the evolution changes from being field-dominated to being matter-dominated. We measure the sheet surface density and the magnetic field strength at that transition for the first time and show consistency with an analytical prediction that had previously gone untested for over 50 years (Mestel 1965).
△ Less
Submitted 23 February, 2023;
originally announced February 2023.
-
JCMT BISTRO Observations: Magnetic Field Morphology of Bubbles Associated with NGC 6334
Authors:
Mehrnoosh Tahani,
Pierre Bastien,
Ray S. Furuya,
Kate Pattle,
Doug Johnstone,
Doris Arzoumanian,
Yasuo Doi,
Tetsuo Hasegawa,
Shu-ichiro Inutsuka,
Simon Coudé,
Laura Fissel,
Michael Chun-Yuan Chen,
Frédérick Poidevin,
Sarah Sadavoy,
Rachel Friesen,
Patrick M. Koch,
James Di Francesco,
Gerald H. Moriarty-Schieven,
Zhiwei Chen,
Eun Jung Chung,
Chakali Eswaraiah,
Lapo Fanciullo,
Tim Gledhill,
Valentin J. M. Le Gouellec,
Thiem Hoang
, et al. (120 additional authors not shown)
Abstract:
We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles…
▽ More
We study the HII regions associated with the NGC 6334 molecular cloud observed in the sub-millimeter and taken as part of the B-fields In STar-forming Region Observations (BISTRO) Survey. In particular, we investigate the polarization patterns and magnetic field morphologies associated with these HII regions. Through polarization pattern and pressure calculation analyses, several of these bubbles indicate that the gas and magnetic field lines have been pushed away from the bubble, toward an almost tangential (to the bubble) magnetic field morphology. In the densest part of NGC 6334, where the magnetic field morphology is similar to an hourglass, the polarization observations do not exhibit observable impact from HII regions. We detect two nested radial polarization patterns in a bubble to the south of NGC 6334 that correspond to the previously observed bipolar structure in this bubble. Finally, using the results of this study, we present steps (incorporating computer vision; circular Hough Transform) that can be used in future studies to identify bubbles that have physically impacted magnetic field lines.
△ Less
Submitted 21 December, 2022;
originally announced December 2022.
-
The JCMT BISTRO-2 Survey: Magnetic Fields of the Massive DR21 Filament
Authors:
Tao-Chung Ching,
Keping Qiu,
Di Li,
Zhiyuan Ren,
Shih-Ping Lai,
David Berry,
Kate Pattle,
Ray Furuya,
Derek Ward-Thompson,
Doug Johnstone,
Patrick M. Koch,
Chang Won Lee,
Thiem Hoang,
Tetsuo Hasegawa,
Woojin Kwon,
Pierre Bastien,
Chakali Eswaraiah,
Jia-Wei Wang,
Kyoung Hee Kim,
Jihye Hwang,
Archana Soam,
A-Ran Lyo,
Junhao Liu,
Valentin J. M. Le Gouellec,
Doris Arzoumanian
, et al. (132 additional authors not shown)
Abstract:
We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parall…
▽ More
We present 850 $μ$m dust polarization observations of the massive DR21 filament from the B-fields In STar-forming Region Observations (BISTRO) survey, using the POL-2 polarimeter and the SCUBA-2 camera on the James Clerk Maxwell Telescope. We detect ordered magnetic fields perpendicular to the parsec-scale ridge of the DR21 main filament. In the sub-filaments, the magnetic fields are mainly parallel to the filamentary structures and smoothly connect to the magnetic fields of the main filament. We compare the POL-2 and Planck dust polarization observations to study the magnetic field structures of the DR21 filament on 0.1--10 pc scales. The magnetic fields revealed in the Planck data are well aligned with those of the POL-2 data, indicating a smooth variation of magnetic fields from large to small scales. The plane-of-sky magnetic field strengths derived from angular dispersion functions of dust polarization are 0.6--1.0 mG in the DR21 filament and $\sim$ 0.1 mG in the surrounding ambient gas. The mass-to-flux ratios are found to be magnetically supercritical in the filament and slightly subcritical to nearly critical in the ambient gas. The alignment between column density structures and magnetic fields changes from random alignment in the low-density ambient gas probed by Planck to mostly perpendicular in the high-density main filament probed by JCMT. The magnetic field structures of the DR21 filament are in agreement with MHD simulations of a strongly magnetized medium, suggesting that magnetic fields play an important role in shaping the DR21 main filament and sub-filaments.
△ Less
Submitted 4 December, 2022;
originally announced December 2022.
-
Three-dimensional Simulations of Magnetospheric Accretion in a T Tauri Star: Accretion and Wind Structures Just Around Star
Authors:
Shinsuke Takasao,
Kengo Tomida,
Kazunari Iwasaki,
Takeru K. Suzuki
Abstract:
We perform three-dimensional magnetohydrodynamic simulations of magnetospheric accretion in a T Tauri star to study the accretion and wind structures in the close vicinity of the star. The gas accreting onto the star consists of the gas from the magnetospheric boundary and the failed disk winds. The accreting gas is commonly found as a multi-column accretion, which is consistent with observations.…
▽ More
We perform three-dimensional magnetohydrodynamic simulations of magnetospheric accretion in a T Tauri star to study the accretion and wind structures in the close vicinity of the star. The gas accreting onto the star consists of the gas from the magnetospheric boundary and the failed disk winds. The accreting gas is commonly found as a multi-column accretion, which is consistent with observations. A significant fraction of the angular momentum of the accreting flows is removed by the magnetic fields of conical disk winds and turbulent failed winds inside and near the magnetosphere. As a result, the accretion torque is significantly reduced compared to the simple estimation based on the mass accretion rate. The stellar spin affects the time variability of the conical disk wind by changing the stability condition of the magnetospheric boundary. However, the time-averaged magnetospheric radius only weakly depends on the stellar spin, which is unlike the prediction of classical theories that the stellar spin controls the magnetospheric radius through the magnetic torque. The ratio of the toroidal to the poloidal field strengths at the magnetospheric boundary, which is a key parameter for the magnetic torque, is also insensitive to the spin; it is rather determined by the disk dynamics. Considering newly found three-dimensional effects, we obtain a scaling relation of the magnetospheric radius very similar to the Ghosh & Lamb relation from the steady angular momentum transport equation.
△ Less
Submitted 2 November, 2022;
originally announced November 2022.
-
The JCMT BISTRO Survey: A Spiral Magnetic Field in a Hub-filament Structure, Monoceros R2
Authors:
Jihye Hwang,
Jongsoo Kim,
Kate Pattle,
Chang Won Lee,
Patrick M. Koch,
Doug Johnstone,
Kohji Tomisaka,
Anthony Whitworth,
Ray S. Furuya,
Ji-hyun Kang,
A-Ran Lyo,
Eun Jung Chung,
Doris Arzoumanian,
Geumsook Park,
Woojin Kwon,
Shinyoung Kim,
Motohide Tamura,
Jungmi Kwon,
Archana Soam,
Ilseung Han,
Thiem Hoang,
Kyoung Hee Kim,
Takashi Onaka,
Eswaraiah Chakali,
Derek Ward-Thompson
, et al. (135 additional authors not shown)
Abstract:
We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R…
▽ More
We present and analyze observations of polarized dust emission at 850 $μ$m towards the central 1 pc $\times$ 1 pc hub-filament structure of Monoceros R2 (Mon R2). The data are obtained with SCUBA-2/POL-2 on the James Clerk Maxwell Telescope (JCMT) as part of the BISTRO (B-fields in Star-forming Region Observations) survey. The orientations of the magnetic field follow the spiral structure of Mon R2, which are well-described by an axisymmetric magnetic field model. We estimate the turbulent component of the magnetic field using the angle difference between our observations and the best-fit model of the underlying large-scale mean magnetic field. This estimate is used to calculate the magnetic field strength using the Davis-Chandrasekhar-Fermi method, for which we also obtain the distribution of volume density and velocity dispersion using a column density map derived from $Herschel$ data and the C$^{18}$O ($J$ = 3-2) data taken with HARP on the JCMT, respectively. We make maps of magnetic field strengths and mass-to-flux ratios, finding that magnetic field strengths vary from 0.02 to 3.64 mG with a mean value of 1.0 $\pm$ 0.06 mG, and the mean critical mass-to-flux ratio is 0.47 $\pm$ 0.02. Additionally, the mean Alfvén Mach number is 0.35 $\pm$ 0.01. This suggests that in Mon R2, magnetic fields provide resistance against large-scale gravitational collapse, and magnetic pressure exceeds turbulent pressure. We also investigate the properties of each filament in Mon R2. Most of the filaments are aligned along the magnetic field direction and are magnetically sub-critical.
△ Less
Submitted 13 December, 2022; v1 submitted 12 October, 2022;
originally announced October 2022.
-
Universal Properties of Dense Clumps in Magnetized Molecular Clouds Formed through Shock Compression of Two-phase Atomic Gases
Authors:
Kazunari Iwasaki,
Kengo Tomida
Abstract:
We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamical simulations,including non-equilibrium chemical reactions, heating/cooling processes, and self-gravity by changing the collision speed $V_0$ and the angle $θ$ between the magnetic field and colliding flow. We found that the efficiency of the dense gas formation depends on $θ$. For small…
▽ More
We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamical simulations,including non-equilibrium chemical reactions, heating/cooling processes, and self-gravity by changing the collision speed $V_0$ and the angle $θ$ between the magnetic field and colliding flow. We found that the efficiency of the dense gas formation depends on $θ$. For small $θ$, anisotropic super-Alfvénic turbulence delays the formation of gravitationally unstable clumps. An increase in $θ$ develops shock-amplified magnetic fields along which the gas is accumulated, making prominent filamentary structures. We further investigate the statistical properties of dense clumps identified with different density thresholds. The statistical properties of the dense clumps with lower densities depend on $V_0$ and $θ$ because their properties are inherited from the global turbulence structure of the molecular clouds. By contrast, denser clumps appear to have asymptotic universal statistical properties, which do not depend on the properties of the colliding flow significantly. The internal velocity dispersions approach subsonic and plasma $β$ becomes order of unity. We develop an analytic formula of the virial parameter which reproduces the simulation results reasonably well. This property may be one of the reasons for the universality of the initial mass function of stars.
△ Less
Submitted 13 June, 2022; v1 submitted 7 June, 2022;
originally announced June 2022.
-
Nature of supersonic turbulence and density distribution function in the multiphase interstellar medium
Authors:
Masato I. N. Kobayashi,
Tsuyoshi Inoue,
Kengo Tomida,
Kazunari Iwasaki,
Hiroki Nakatsugawa
Abstract:
Supersonic flows in the interstellar medium (ISM) are believed to be a key driver of the molecular cloud formation and evolution. Among molecular clouds' properties, the ratio between the solenoidal and compressive modes of turbulence plays important roles in determining the star formation efficiency. We use numerical simulations of supersonic converging flows of the warm neutral medium (WNM) reso…
▽ More
Supersonic flows in the interstellar medium (ISM) are believed to be a key driver of the molecular cloud formation and evolution. Among molecular clouds' properties, the ratio between the solenoidal and compressive modes of turbulence plays important roles in determining the star formation efficiency. We use numerical simulations of supersonic converging flows of the warm neutral medium (WNM) resolving the thermal instability to calculate the early phase of molecular cloud formation, and investigate the turbulence structure and the density probability distribution function (density PDF) of the multiphase ISM. We find that both the solenoidal and compressive modes have their power spectrum similar to the Kolmogorov spectrum. The solenoidal (compressive) modes account for >~80% (<~20%) of the total turbulence power. When we consider both the cold neutral medium (CNM) and the thermally unstable neutral medium (UNM) up to T <~ 400 K, the density PDF follows the log-normal distribution whose width sigma_s is well explained by the known relation from the isothermal turbulence as sigma_s = ln(1 + b^2 * M^2) (where b is the parameter representing the turbulence mode ratio and M is the turbulent Mach number). The density PDF of the CNM component alone (T <= 50 K), however, exhibits a narrower sigma_s by a factor of ~ 2. These results suggest that observational estimations of b based on the CNM density PDF requires the internal turbulence within each CNM clump but not the inter-clump relative velocity, the latter of which is instead powered by the WNM/UNM turbulence.
△ Less
Submitted 1 March, 2022;
originally announced March 2022.
-
B-fields in Star-Forming Region Observations (BISTRO): Magnetic Fields in the Filamentary Structures of Serpens Main
Authors:
Woojin Kwon,
Kate Pattle,
Sarah Sadavoy,
Charles L. H. Hull,
Doug Johnstone,
Derek Ward-Thompson,
James Di Francesco,
Patrick M. Koch,
Ray Furuya,
Yasuo Doi,
Valentin J. M. Le Gouellec,
Jihye Hwang,
A-Ran Lyo,
Archana Soam,
Xindi Tang,
Thiem Hoang,
Florian Kirchschlager,
Chakali Eswaraiah,
Lapo Fanciullo,
Kyoung Hee Kim,
Takashi Onaka,
Vera Könyves,
Ji-hyun Kang,
Chang Won Lee,
Motohide Tamura
, et al. (127 additional authors not shown)
Abstract:
We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filament…
▽ More
We present 850 $μ$m polarimetric observations toward the Serpens Main molecular cloud obtained using the POL-2 polarimeter on the James Clerk Maxwell Telescope (JCMT) as part of the B-fields In STar-forming Region Observations (BISTRO) survey. These observations probe the magnetic field morphology of the Serpens Main molecular cloud on about 6000 au scales, which consists of cores and six filaments with different physical properties such as density and star formation activity. Using the histogram of relative orientation (HRO) technique, we find that magnetic fields are parallel to filaments in less dense filamentary structures where $N_{H_2} < 0.93\times 10^{22}$ cm$^{-2}$ (magnetic fields perpendicular to density gradients), while being perpendicular to filaments (magnetic fields parallel to density gradients) in dense filamentary structures with star formation activity. Moreover, applying the HRO technique to denser core regions, we find that magnetic field orientations change to become perpendicular to density gradients again at $N_{H_2} \approx 4.6 \times 10^{22}$ cm$^{-2}$. This can be interpreted as a signature of core formation. At $N_{H_2} \approx 16 \times 10^{22}$ cm$^{-2}$ magnetic fields change back to being parallel to density gradients once again, which can be understood to be due to magnetic fields being dragged in by infalling material. In addition, we estimate the magnetic field strengths of the filaments ($B_{POS} = 60-300~μ$G)) using the Davis-Chandrasekhar-Fermi method and discuss whether the filaments are gravitationally unstable based on magnetic field and turbulence energy densities.
△ Less
Submitted 13 January, 2022;
originally announced January 2022.
-
K3 Surfaces, Picard Numbers and Siegel Disks
Authors:
Katsunori Iwasaki,
Yuta Takada
Abstract:
If a K3 surface admits an automorphism with a Siegel disk, then its Picard number is an even integer between $0$ and $18$. Conversely, using the method of hypergeometric groups, we are able to construct K3 surface automorphisms with Siegel disks that realize all possible Picard numbers. The constructions involve extensive computer searches for appropriate Salem numbers and computations of algebrai…
▽ More
If a K3 surface admits an automorphism with a Siegel disk, then its Picard number is an even integer between $0$ and $18$. Conversely, using the method of hypergeometric groups, we are able to construct K3 surface automorphisms with Siegel disks that realize all possible Picard numbers. The constructions involve extensive computer searches for appropriate Salem numbers and computations of algebraic numbers arising from holomorphic Lefschetz-type formulas and related Grothendieck residues.
△ Less
Submitted 5 November, 2021;
originally announced November 2021.
-
The JCMT BISTRO Survey: An 850/450$μ$m Polarization Study of NGC 2071IR in OrionB
Authors:
A-Ran Lyo,
Jongsoo Kim,
Sarah Sadavoy,
Doug Johnstone,
David Berry,
Kate Pattle,
Woojin Kwon,
Pierre Bastien,
Takashi Onaka,
James Di Francesco,
Ji-Hyun Kang,
Ray Furuya,
Charles L. H. Hull,
Motohide Tamura,
Patrick M. Koch,
Derek Ward-Thompson,
Tetsuo Hasegawa,
Thiem Hoang,
Doris Arzoumanian,
Chang Won Lee,
Chin-Fei Lee,
Do-Young Byun,
Florian Kirchschlager,
Yasuo Doi,
Kee-Tae Kim
, et al. (121 additional authors not shown)
Abstract:
We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could b…
▽ More
We present the results of simultaneous 450 $μ$m and 850 $μ$m polarization observations toward the massive star forming region NGC 2071IR, a target of the BISTRO (B-fields in Star-Forming Region Observations) Survey, using the POL-2 polarimeter and SCUBA-2 camera mounted on the James Clerk Maxwell Telescope. We find a pinched magnetic field morphology in the central dense core region, which could be due to a rotating toroidal disk-like structure and a bipolar outflow originating from the central young stellar object, IRS 3. Using the modified Davis-Chandrasekhar-Fermi method, we obtain a plane-of-sky magnetic field strength of 563$\pm$421 $μ$G in the central $\sim$0.12 pc region from 850 $μ$m polarization data. The corresponding magnetic energy density of 2.04$\times$10$^{-8}$ erg cm$^{-3}$ is comparable to the turbulent and gravitational energy densities in the region. We find that the magnetic field direction is very well aligned with the whole of the IRS 3 bipolar outflow structure. We find that the median value of polarization fractions, 3.0 \%, at 450 $μ$m in the central 3 arcminute region, which is larger than the median value of 1.2 \% at 850 $μ$m. The trend could be due to the better alignment of warmer dust in the strong radiation environment. We also find that polarization fractions decrease with intensity at both wavelengths, with slopes, determined by fitting a Rician noise model, of $0.59 \pm 0.03$ at 450 $μ$m and $0.36 \pm 0.04$ at 850 $μ$m, respectively. We think that the shallow slope at 850 $μ$m is due to grain alignment at the center being assisted by strong radiation from the central young stellar objects.
△ Less
Submitted 28 September, 2021;
originally announced September 2021.
-
Revealing the diverse magnetic field morphologies in Taurus dense cores with sensitive sub-millimeter polarimetry
Authors:
Chakali Eswaraiah,
Di Li,
Ray S. Furuya,
Tetsuo Hasegawa,
Derek Ward-Thompson,
Keping Qiu,
Nagayoshi Ohashi,
Kate Pattle,
Sarah Sadavoy,
Charles L. H. Hull,
David Berry,
Yasuo Doi,
Tao-Chung Ching,
Shih-Ping Lai,
Jia-Wei Wang,
Patrick M. Koch,
Jungmi Kwon,
Woojin Kwon,
Pierre Bastien,
Doris Arzoumanian,
Simon Coudé,
Archana Soam,
Lapo Fanciullo,
Hsi-Wei Yen,
Junhao Liu
, et al. (120 additional authors not shown)
Abstract:
We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protost…
▽ More
We have obtained sensitive dust continuum polarization observations at 850 $μ$m in the B213 region of Taurus using POL-2 on SCUBA-2 at the James Clerk Maxwell Telescope (JCMT), as part of the BISTRO (B-fields in STar-forming Region Observations) survey. These observations allow us to probe magnetic field (B-field) at high spatial resolution ($\sim$2000 au or $\sim$0.01 pc at 140 pc) in two protostellar cores (K04166 and K04169) and one prestellar core (Miz-8b) that lie within the B213 filament. Using the Davis-Chandrasekhar-Fermi method, we estimate the B-field strengths in K04166, K04169, and Miz-8b to be 38$\pm$14 $μ$G, 44$\pm$16 $μ$G, and 12$\pm$5 $μ$G, respectively. These cores show distinct mean B-field orientations. B-field in K04166 is well ordered and aligned parallel to the orientations of the core minor axis, outflows, core rotation axis, and large-scale uniform B-field, in accordance with magnetically regulated star formation via ambipolar diffusion taking place in K04166. B-field in K04169 is found to be ordered but oriented nearly perpendicular to the core minor axis and large-scale B-field, and not well-correlated with other axes. In contrast, Miz-8b exhibits disordered B-field which show no preferred alignment with the core minor axis or large-scale field. We found that only one core, K04166, retains a memory of the large-scale uniform B-field. The other two cores, K04169 and Miz-8b, are decoupled from the large-scale field. Such a complex B-field configuration could be caused by gas inflow onto the filament, even in the presence of a substantial magnetic flux.
△ Less
Submitted 3 March, 2021;
originally announced March 2021.
-
Dust polarized emission observations of NGC 6334; BISTRO reveals the details of the complex but organized magnetic field structure of the high-mass star-forming hub-filament network
Authors:
D. Arzoumanian,
R. Furuya,
T. Hasegawa,
M. Tahani,
S. Sadavoy,
C. L. H. Hull,
D. Johnstone,
P. M. Koch,
S. -i. Inutsuka,
Y. Doi,
T. Hoang,
T. Onaka,
K. Iwasaki,
Y. Shimajiri,
T. Inoue,
N. Peretto,
P. André,
P. Bastien,
D. Berry,
H. -R. V. Chen,
J. Di Francesco,
C. Eswaraiah,
L. Fanciullo,
L. M. Fissel,
J. Hwang
, et al. (123 additional authors not shown)
Abstract:
[Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We presen…
▽ More
[Abridged] Filaments and hubs have received special attention recently thanks to studies showing their role in star formation. While the column density and velocity structures of both filaments and hubs have been studied, their magnetic fields (B-field) are not yet characterized. We aim to understand the role of the B-field in the dynamical evolution of the NGC 6334 hub-filament network. We present new observations of the dust polarized emission at 850$μ$m towards NGC 6334 obtained with the JCMT/POL-2. We study the distribution and dispersion of the polarized intensity ($PI$), the polarization fraction ($PF$), and the B-field angle ($θ_{B}$). We derive the power spectrum of the intensity and $θ_{B}$ along the ridge crest. Our analyses show a complex B-field structure when observed over the whole region ($\sim10$ pc), however, at smaller scales ($\sim1$ pc), $θ_{B}$ varies coherently along the filaments. The observed power spectrum of $θ_{B}$ can be well represented with a power law function with a slope $-1.33\pm0.23$, which is $\sim20\%$ shallower than that of $I$. This result is compatible with the properties of simulated filaments and may indicate the processes at play in the formation of filaments. $θ_{B}$ rotates from being mostly perpendicular to the filament crests to mostly parallel as they merge with the hubs. This variation of $θ_{B}$ may be tracing local velocity flows of matter in-falling onto the hubs. Our analysis suggests a variation of the energy balance along the crests of these filaments, from magnetically critical/supercritical at their far ends to magnetically subcritical near the hubs. We detect an increase of $PF$ towards the high-column density star cluster-forming hubs that may result from the increase of grain alignment efficiency due to stellar radiation from the newborn stars.
△ Less
Submitted 23 December, 2020;
originally announced December 2020.
-
Observations of magnetic fields surrounding LkH$α$ 101 taken by the BISTRO survey with JCMT-POL-2
Authors:
Nguyen Bich Ngoc,
Pham Ngoc Diep,
Harriet Parsons,
Kate Pattle,
Thiem Hoang,
Derek Ward-Thompson,
Le Ngoc Tram,
Charles L. H. Hull,
Mehrnoosh Tahani,
Ray Furuya,
Pierre Bastien,
Keping Qiu,
Tetsuo Hasegawa,
Woojin Kwon,
Yasuo Doi,
Shih-Ping Lai,
Simon Coude,
David Berry,
Tao-Chung Ching,
Jihye Hwang,
Archana Soam,
Jia-Wei Wang,
Doris Arzoumanian,
Tyler L. Bourke,
Do-Young Byun
, et al. (124 additional authors not shown)
Abstract:
We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to…
▽ More
We report the first high spatial resolution measurement of magnetic fields surrounding LkH$α$ 101, a part of the Auriga-California molecular cloud. The observations were taken with the POL-2 polarimeter on the James Clerk Maxwell Telescope within the framework of the B-fields In Star-forming Region Observations (BISTRO) survey. Observed polarization of thermal dust emission at 850 $μ$m is found to be mostly associated with the red-shifted gas component of the cloud. The magnetic field displays a relatively complex morphology. Two variants of the Davis-Chandrasekhar-Fermi method, unsharp masking and structure function, are used to calculate the strength of magnetic fields in the plane of the sky, yielding a similar result of $B_{\rm POS}\sim 115$ $\mathrmμ$G. The mass-to-magnetic-flux ratio in critical value units, $λ\sim0.3$, is the smallest among the values obtained for other regions surveyed by POL-2. This implies that the LkH$α$ 101 region is sub-critical and the magnetic field is strong enough to prevent gravitational collapse. The inferred $δB/B_0\sim 0.3$ implies that the large scale component of the magnetic field dominates the turbulent one. The variation of the polarization fraction with total emission intensity can be fitted by a power-law with an index of $α=0.82\pm0.03$, which lies in the range previously reported for molecular clouds. We find that the polarization fraction decreases rapidly with proximity to the only early B star (LkH$α$ 101) in the region. The magnetic field tangling and the joint effect of grain alignment and rotational disruption by radiative torques are potential of explaining such a decreasing trend.
△ Less
Submitted 8 December, 2020;
originally announced December 2020.
-
Bimodal Behavior and Convergence Requirement in Macroscopic Properties of the Multiphase Interstellar Medium Formed by Atomic Converging Flows
Authors:
Masato I. N. Kobayashi,
Tsuyoshi Inoue,
Shu-Ichiro Inutsuka,
Kengo Tomida,
Kazunari Iwasaki,
Kei E. I. Tanaka
Abstract:
We systematically perform hydrodynamics simulations of 20 km s^-1 converging flows of the warm neutral medium (WNM) to calculate the formation of the cold neutral medium (CNM), especially focusing on the mean properties of the multiphase interstellar medium (ISM), such as the average shock front position and the mean density on a 10 pc scale. Our results show that the convergence in those mean pro…
▽ More
We systematically perform hydrodynamics simulations of 20 km s^-1 converging flows of the warm neutral medium (WNM) to calculate the formation of the cold neutral medium (CNM), especially focusing on the mean properties of the multiphase interstellar medium (ISM), such as the average shock front position and the mean density on a 10 pc scale. Our results show that the convergence in those mean properties requires 0.02 pc spatial resolution that resolves the cooling length of the thermally unstable neutral medium (UNM) to follow the dynamical condensation from the WNM to CNM. We also find that two distinct post-shock states appear in the mean properties depending on the amplitude of the upstream WNM density fluctuation (= sqrt(<drho^2>)/rho_0). When the amplitude > 10 %, the interaction between shocks and density inhomogeneity leads to a strong driving of the post-shock turbulence of > 3 km s^-1, which dominates the energy budget in the shock-compressed layer. The turbulence prevents the dynamical condensation by cooling and the following CNM formation, and the CNM mass fraction remains as ~ 45 %. In contrast, when the amplitude <= 10 %, the shock fronts maintain an almost straight geometry and CNM formation efficiently proceeds, resulting in the CNM mass fraction of ~ 70 %. The velocity dispersion is limited to the thermal-instability mediated level of 2 - 3 km s^-1 and the layer is supported by both turbulent and thermal energy equally. We also propose an effective equation of state that models the multiphase ISM formed by the WNM converging flow as a one-phase ISM.
△ Less
Submitted 23 October, 2020;
originally announced October 2020.
-
The JCMT BISTRO Survey: Magnetic Fields Associated with a Network of Filaments in NGC 1333
Authors:
Yasuo Doi,
Tetsuo Hasegawa,
Ray S. Furuya,
Simon Coudé,
Charles L. H. Hull,
Doris Arzoumanian,
Pierre Bastien,
Michael Chun-Yuan Chen,
James di Francesco,
Rachel Friesen,
Martin Houde,
Shu-ichiro Inutsuka,
Steve Mairs,
Masafumi Matsumura,
Takashi Onaka,
Sarah Sadavoy,
Yoshito Shimajiri,
Mehrnoosh Tahani,
Kohji Tomisaka,
Chakali Eswaraiah,
Patrick M. Koch,
Kate Pattle,
Chang Won Lee,
Motohide Tamura,
David Berry
, et al. (113 additional authors not shown)
Abstract:
We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary…
▽ More
We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multiscale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~1 pc and remains continuous from the scales of filaments (~0.1 pc) to that of protostellar envelopes (~0.005 pc or ~1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data toward the filament network.
△ Less
Submitted 21 July, 2020; v1 submitted 30 June, 2020;
originally announced July 2020.
-
Hypergeometric Groups and Dynamics on K3 Surfaces
Authors:
Katsunori Iwasaki,
Yuta Takada
Abstract:
A hypergeometric group is a matrix group modeled on the monodromy group of a generalized hypergeometric differential equation. This article presents a fruitful interaction between the theory of hypergeometric groups and dynamics on K3 surfaces by showing that a certain class of hypergeometric groups and related lattices lead to a lot of K3 surface automorphisms of positive entropy, especially such…
▽ More
A hypergeometric group is a matrix group modeled on the monodromy group of a generalized hypergeometric differential equation. This article presents a fruitful interaction between the theory of hypergeometric groups and dynamics on K3 surfaces by showing that a certain class of hypergeometric groups and related lattices lead to a lot of K3 surface automorphisms of positive entropy, especially such automorphisms with Siegel disks.
△ Less
Submitted 30 April, 2021; v1 submitted 30 March, 2020;
originally announced March 2020.
-
First sub-arcsecond submillimeter-wave [C I] image of 49 Ceti with ALMA
Authors:
Aya E. Higuchi,
Kazuya Saigo,
Hiroshi Kobayashi,
Kazunari Iwasaki,
Munetake Momose,
Kang Lou Soon,
Nami Sakai,
Masanobu Kunitomo,
Daisuke Ishihara,
Satoshi Yamamoto
Abstract:
We present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substan…
▽ More
We present the first sub-arcsecond images of 49 Ceti in the [C I]3P1-3P0 emission and the 614 micron dust continuum emission observed with ALMA, as well as that in the CO(J=3-2) emission prepared by using the ALMA archival data. The spatial distribution of the 614 micron dust continuum emission is found to have a broad-ring structure with a radius of about 100 au around the central star. A substantial amount of gas is also associated with 49 Ceti. The [C I] emission map shows two peaks inside the dust ring, and its overall extent is comparable to that of the dust continuum emission and the CO emission. We find that the [C I]/CO(J=3-2) intensity ratio significantly varies along the major axis. The ratio takes the minimum value of 1.8 around the dust peak position, and increases inwards and outwards. The enhanced ratio around the central star (~3) likely originates from the stellar UV radiation, while that in the outer disk (~10) from the interstellar UV radiation. Such complex distributions of the [C I] and CO(J=3-2) emission will be a key to understand the origin of the gas in 49 Ceti, and will also provide a stringent constraint on physical and chemical models of gaseous debris disks.
△ Less
Submitted 19 August, 2019;
originally announced August 2019.
-
The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region
Authors:
Simon Coudé,
Pierre Bastien,
Martin Houde,
Sarah Sadavoy,
Rachel Friesen,
James Di Francesco,
Doug Johnstone,
Steve Mairs,
Tetsuo Hasegawa,
Woojin Kwon,
Shih-Ping Lai,
Keping Qiu,
Derek Ward-Thompson,
David Berry,
Michael Chun-Yuan Chen,
Jason Fiege,
Erica Franzmann,
Jennifer Hatchell,
Kevin Lacaille,
Brenda C. Matthews,
Gerald H. Moriarty-Schieven,
Andy Pon,
Philippe André,
Doris Arzoumanian,
Yusuke Aso
, et al. (96 additional authors not shown)
Abstract:
We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We th…
▽ More
We present the POL-2 850 $μ$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $μ$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $μ$G and a criticality criterion $λ_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.
△ Less
Submitted 23 April, 2019; v1 submitted 15 April, 2019;
originally announced April 2019.
-
Discrete Laplace Method and Truncation Error of Gauss Continued Fraction
Authors:
Katsunori Iwasaki
Abstract:
The leading asymptotics of the truncation error for Gauss's continued fraction is determined exactly. Not only for this purpose but also for wider applicability elsewhere the discrete analogue of Laplace's method for hypergeometric series containing a large parameter, which was developed in a previous paper, is generalized in two directions.
The leading asymptotics of the truncation error for Gauss's continued fraction is determined exactly. Not only for this purpose but also for wider applicability elsewhere the discrete analogue of Laplace's method for hypergeometric series containing a large parameter, which was developed in a previous paper, is generalized in two directions.
△ Less
Submitted 5 April, 2019;
originally announced April 2019.
-
The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $ρ$ Ophiuchus C
Authors:
Junhao Liu,
Keping Qiu,
David Berry,
James Di Francesco,
Pierre Bastien,
Patrick M. Koch,
Ray S. Furuya,
Kee-Tae Kim,
Simon Coudé,
Chang Won Lee,
Archana Soam,
Chakali Eswaraiah,
Di Li,
Jihye Hwang,
A-Ran Lyo,
Kate Pattle,
Tetsuo Hasegawa,
Woojin Kwon,
Shih-Ping Lai,
Derek Ward-Thompson,
Tao-Chung Ching,
Zhiwei Chen,
Qilao Gu,
Dalei Li,
Hua-bai Li
, et al. (106 additional authors not shown)
Abstract:
We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic…
▽ More
We report 850~$μ$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $ρ$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $μ$G, 136 $\pm$ 69 $μ$G, and 213 $\pm$ 115 $μ$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.
△ Less
Submitted 20 February, 2019;
originally announced February 2019.
-
Giant protostellar flares: accretion-driven accumulation and reconnection-driven ejection of magnetic flux in protostars
Authors:
Shinsuke Takasao,
Kengo Tomida,
Kazunari Iwasaki,
Takeru K. Suzuki
Abstract:
Protostellar flares are rapid magnetic energy release events associated with formation of hot plasma in protostars. In the previous models of protostellar flares, the interaction between a protostellar magnetosphere with the surrounding disk plays crucial roles in building-up and releasing the magnetic energy. However, it remains unclear if protostars indeed have magnetospheres because vigorous di…
▽ More
Protostellar flares are rapid magnetic energy release events associated with formation of hot plasma in protostars. In the previous models of protostellar flares, the interaction between a protostellar magnetosphere with the surrounding disk plays crucial roles in building-up and releasing the magnetic energy. However, it remains unclear if protostars indeed have magnetospheres because vigorous disk accretion and strong disk magnetic fields in the protostellar phase may destroy the magnetosphere. Considering this possibility, we investigate the energy accumulation and release processes in the absence of a magnetosphere using a three-dimensional magnetohydrodynamic simulation. Our simulation reveals that protostellar flares are repeatedly produced even in such a case. Unlike in the magnetospheric models, the protostar accumulates magnetic energy by acquiring large-scale magnetic fields from the disk by accretion. Protostellar flares occur when a portion of the large-scale magnetic fields are removed from the protostar as a result of magnetic reconnection. Protostellar flares in the simulation are consistent with observations; the released magnetic energy (up to $\sim 3\times 10^{38}$ erg) is large enough to drive observed flares, and the flares produce hot ejecta. The expelled magnetic fields enhance accretion, and the energy build-up and release processes are repeated as a result. The magnetic flux removal via reconnection leads to redistribution of magnetic fields in the inner disk. We therefore consider that protostellar flares will play an important role in the evolution of the disk magnetic fields in the vicinity of protostars.
△ Less
Submitted 21 May, 2019; v1 submitted 5 February, 2019;
originally announced February 2019.
-
JCMT BISTRO survey: Magnetic Fields within the Hub-Filament Structure in IC 5146
Authors:
Jia-Wei Wang,
Shih-Ping Lai,
Chakali Eswaraiah,
Kate Pattle,
James Di Francesco,
Doug Johnstone,
Patrick M. Koch,
Tie Liu,
Motohide Tamura,
Ray S. Furuya,
Takashi Onaka,
Derek Ward-Thompson,
Archana Soam,
Kee-Tae Kim,
Chang Won Lee,
Chin-Fei Lee,
Steve Mairs,
Doris Arzoumanian,
Gwanjeong Kim,
Thiem Hoang,
Jihye Hwang,
Sheng-Yuan Liu,
David Berry,
Pierre Bastien,
Tetsuo Hasegawa
, et al. (108 additional authors not shown)
Abstract:
We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scal…
▽ More
We present the 850 $μ$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfvénic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.
△ Less
Submitted 27 March, 2019; v1 submitted 14 December, 2018;
originally announced December 2018.
-
Does misalignment between magnetic field and angular momentum enhance or suppress circumstellar disk formation?
Authors:
Y. Tsukamoto,
S. Okuzumi,
K. Iwasaki,
M. N. Machida,
S. Inutsuka
Abstract:
The effect of misalignment between the magnetic field $\magB$ and the angular momentum $\Jang$ of molecular cloud cores on the angular momentum evolution during the gravitational collapse is investigated by ideal and non-ideal MHD simulations. For the non-ideal effect, we consider the ohmic and ambipolar diffusion. Previous studies that considered the misalignment reported qualitatively contradict…
▽ More
The effect of misalignment between the magnetic field $\magB$ and the angular momentum $\Jang$ of molecular cloud cores on the angular momentum evolution during the gravitational collapse is investigated by ideal and non-ideal MHD simulations. For the non-ideal effect, we consider the ohmic and ambipolar diffusion. Previous studies that considered the misalignment reported qualitatively contradicting results. Magnetic braking was reported as being either strengthened or weakened by misalignment in different studies. We conducted simulations of cloud-core collapse by varying the stability parameter $α$ (the ratio of the thermal to gravitational energy of the core) with and without including magnetic diffusion The non-ideal MHD simulations show the central angular momentum of the core with $θ=0^\circ$ ($\Jang \parallel \magB$) being always greater than that with $θ=90^\circ$ ($\Jang \perp \magB$), independently of $α$, meaning that circumstellar disks form more easily form in a core with $θ=0^\circ$. The ideal MHD simulations, in contrast, show the the central angular momentum of the core with $θ=90^\circ$ being greater than with $θ=0^\circ$ for small $α$, and is smaller for large $α$. Inspection of the angular momentum evolution of the fluid elements reveals three mechanisms contributing to the evolution of the angular momentum: (i) magnetic braking in the isothermal collapse phase, (ii) selective accretion of the rapidly (for $θ=90^\circ$ ) or slowly (for $θ=0^\circ$) rotating fluid elements to the central region, and (iii) magnetic braking in the first-core and the disk. The difference between the ideal and non-ideal simulations arises from the different efficiencies of (iii).
△ Less
Submitted 16 November, 2018;
originally announced November 2018.
-
The Early Stage of Molecular Cloud Formation by Compression of Two-phase Atomic Gases
Authors:
Kazunari Iwasaki,
Kengo Tomida,
Tsuyoshi Inoue,
Shu-ichiro Inutsuka
Abstract:
We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamic simulations, including non-equilibrium chemical reactions and heating/cooling processes. We consider super-Alfvénic head-on colliding flows of atomic gas possessing the two-phase structure that consists of HI clouds and surrounding warm diffuse gas. We examine how the formation of molecu…
▽ More
We investigate the formation of molecular clouds from atomic gas by using three-dimensional magnetohydrodynamic simulations, including non-equilibrium chemical reactions and heating/cooling processes. We consider super-Alfvénic head-on colliding flows of atomic gas possessing the two-phase structure that consists of HI clouds and surrounding warm diffuse gas. We examine how the formation of molecular clouds depends on the angle $θ$ between the upstream flow and the mean magnetic field. We find that there is a critical angle $θ_\mathrm{cr}$ above which the shock-amplified magnetic field controls the post-shock gas dynamics. If the atomic gas is compressed almost along the mean magnetic field ($θ\llθ_\mathrm{cr}$), super-Alfvénic anisotropic turbulence is maintained by the accretion of the highly inhomogeneous upstream atomic gas. As a result, a greatly extended turbulence-dominated post-shock layer is generated. Around $θ\sim θ_\mathrm{cr}$, the shock-amplified magnetic field weakens the post-shock turbulence, leading to a dense post-shock layer. For $θ\gg θ_\mathrm{cr}$, the strong magnetic pressure suppresses the formation of cold dense clouds. Efficient molecular cloud formation is expected if $θ$ is less than a few times $θ_\mathrm{cr}$. Developing an analytic model and performing a parameter survey, we obtain an analytic formula for the critical angle as a function of the mean density, collision speed, and field strength of the upstream atomic gas. The critical angle is found to be less than $\sim 15^\circ$ as long as the field strength is larger than $1~μ$G, indicating that the probability of occurrence of compression with $θ<θ_\mathrm{cr}$ is limited if shock waves come from various directions.
△ Less
Submitted 15 April, 2019; v1 submitted 11 June, 2018;
originally announced June 2018.
-
Magnetic fields towards Ophiuchus-B derived from SCUBA-2 polarization measurements
Authors:
Archana Soam,
Kate Pattle,
Derek Ward-Thompson,
Chang Won Lee,
Sarah Sadavoy,
Patrick M. Koch,
Gwanjeong Kim,
Jungmi Kwon,
Woojin Kwon,
Doris Arzoumanian,
David Berry,
Thiem Hoang,
Motohide Tamura,
Sang-Sung Lee,
Tie Liu,
Kee-Tae Kim,
Doug Johnstone,
Fumitaka Nakamura,
A-Ran Lyo,
Takashi Onaka,
Jongsoo Kim,
Ray S. Furuya,
Tetsuo Hasegawa,
Shih-Ping Lai,
Pierre Bastien
, et al. (99 additional authors not shown)
Abstract:
We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fi…
▽ More
We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $μ$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $λ$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.
△ Less
Submitted 7 June, 2018; v1 submitted 16 May, 2018;
originally announced May 2018.
-
A First Look at BISTRO Observations of The $ρ$ Oph-A core
Authors:
Jungmi Kwon,
Yasuo Doi,
Motohide Tamura,
Masafumi Matsumura,
Kate Pattle,
David Berry,
Sarah Sadavoy,
Brenda C. Matthews,
Derek Ward-Thompson,
Tetsuo Hasegawa,
Ray S. Furuya,
Andy Pon,
James Di Francesco,
Doris Arzoumanian,
Saeko S. Hayashi,
Koji S. Kawabata,
Takashi Onaka,
Minho Choi,
Miju Kang,
Thiem Hoang,
Chang Won Lee,
Sang-Sung Lee,
Hong-Li Liu,
Tie Liu,
Shu-Ichiro Inutsuka
, et al. (97 additional authors not shown)
Abstract:
We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We iden…
▽ More
We present 850 $μ$m imaging polarimetry data of the $ρ$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $ρ$ Oph A core; some of them can be part of a coherent magnetic field structure in the $ρ$ Oph region. The results are consistent with previous observations of the brightest regions of $ρ$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.
△ Less
Submitted 24 April, 2018;
originally announced April 2018.
-
Incidence structures near configurations of type $(n_3)$
Authors:
Peter Dukes,
Kaoruko Iwasaki
Abstract:
An $(n_3)$ configuration is an incidence structure equivalent to a linear hypergraph on $n$ vertices which is both 3-regular and 3-uniform. We investigate a variant in which one constraint, say 3-regularity, is present, and we allow exactly one line to have size four, exactly one line to have size two, and all other lines to have size three. In particular, we study planar (Euclidean or projective)…
▽ More
An $(n_3)$ configuration is an incidence structure equivalent to a linear hypergraph on $n$ vertices which is both 3-regular and 3-uniform. We investigate a variant in which one constraint, say 3-regularity, is present, and we allow exactly one line to have size four, exactly one line to have size two, and all other lines to have size three. In particular, we study planar (Euclidean or projective) representations, settling the existence question and adapting Steinitz' theorem for this setting.
△ Less
Submitted 24 April, 2018; v1 submitted 31 January, 2018;
originally announced January 2018.
-
A Three-dimensional Simulation of a Magnetized Accretion Disk: Fast Funnel Accretion onto a Weakly Magnetized Star
Authors:
Shinsuke Takasao,
Kengo Tomida,
Kazunari Iwasaki,
Takeru K. Suzuki
Abstract:
We present the results of a global, three-dimensional magnetohydrodynamics simulation of an accretion disk with a rotating, weakly magnetized central star. The disk is threaded by a weak, large-scale poloidal magnetic field, and the central star has no strong stellar magnetosphere initially. Our simulation investigates the structure of the accretion flows from a turbulent accretion disk onto the s…
▽ More
We present the results of a global, three-dimensional magnetohydrodynamics simulation of an accretion disk with a rotating, weakly magnetized central star. The disk is threaded by a weak, large-scale poloidal magnetic field, and the central star has no strong stellar magnetosphere initially. Our simulation investigates the structure of the accretion flows from a turbulent accretion disk onto the star. The simulation reveals that fast accretion onto the star at high latitudes occurs even without a stellar magnetosphere. We find that the failed disk wind becomes the fast, high-latitude accretion as a result of angular momentum exchange mediated by magnetic fields well above the disk, where the Lorentz force that decelerates the rotational motion of gas can be comparable to the centrifugal force. Unlike the classical magnetospheric accretion scenario, fast accretion streams are not guided by magnetic fields of the stellar magnetosphere. Nevertheless, the accretion velocity reaches the free-fall velocity at the stellar surface due to the efficient angular momentum loss at a distant place from the star. This study provides a possible explanation why Herbig Ae/Be stars whose magnetic fields are generally not strong enough to form magnetospheres also show indications of fast accretion. A magnetically driven jet is not formed from the disk in our model. The differential rotation cannot generate sufficiently strong magnetic fields for the jet acceleration because the Parker instability interrupts the field amplification.
△ Less
Submitted 6 February, 2019; v1 submitted 22 January, 2018;
originally announced January 2018.
-
Contiguous Relations, Laplace's Methods and Continued Fractions for 3F2(1)
Authors:
Akihito Ebisu,
Katsunori Iwasaki
Abstract:
Using contiguous relations we construct an infinite number of continued fraction expansions for ratios of generalized hypergeometric series 3F2(1). We establish exact error term estimates for their approximants and prove their rapid convergences. To do so we develop a discrete version of Laplace's method for hypergeometric series in addition to the use of ordinary (continuous) Laplace's method for…
▽ More
Using contiguous relations we construct an infinite number of continued fraction expansions for ratios of generalized hypergeometric series 3F2(1). We establish exact error term estimates for their approximants and prove their rapid convergences. To do so we develop a discrete version of Laplace's method for hypergeometric series in addition to the use of ordinary (continuous) Laplace's method for Euler's hypergeometric integrals.
△ Less
Submitted 1 September, 2017;
originally announced September 2017.
-
The Formation of Massive Molecular Filaments and Massive Stars Triggered by a MHD Shock Wave
Authors:
Tsuyoshi Inoue,
Patrick Hennebelle,
Yasuo Fukui,
Tomoaki Matsumoto,
Kazunari Iwasaki,
Shu-ichiro Inutsuka
Abstract:
Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adapt…
▽ More
Recent observations suggest that intensive molecular cloud collision can trigger massive star/cluster formation. The most important physical process caused by the collision is a shock compression. In this paper, the influence of a shock wave on the evolution of a molecular cloud is studied numerically by using isothermal magnetohydrodynamics (MHD) simulations with the effect of self-gravity. Adaptive-mesh-refinement and sink particle techniques are used to follow long-time evolution of the shocked cloud. We find that the shock compression of turbulent inhomogeneous molecular cloud creates massive filaments, which lie perpendicularly to the background magnetic field as we have pointed out in a previous paper. The massive filament shows global collapse along the filament, which feeds a sink particle located at the collapse center. We observe high accretion rate dot{M}_acc > 10^{-4} M_sun/yr that is high enough to allow the formation of even O-type stars. The most massive sink particle achieves M>50 M_sun in a few times 10^5 yr after the onset of the filament collapse.
△ Less
Submitted 3 August, 2017; v1 submitted 7 July, 2017;
originally announced July 2017.
-
The impact of the Hall effect during cloud core collapse:implications for circumstellar disk evolution
Authors:
Yusuke Tsukamoto,
Satoshi Okuzumi,
Kazunari Iwasaki,
Masahiro N. Machida,
Shu-ichiro Inutsuka
Abstract:
We perform three-dimensional radiation non-ideal magnetohydrodynamics simulations and investigate the impact of the Hall effect on the angular momentum evolution in the collapsing cloud cores in which the magnetic field $\magB$ and angular momentum $\angJ$ are misaligned with each other. We find that the Hall effect notably changes the magnetic torques in the pseudo-disk, and strengthens and weake…
▽ More
We perform three-dimensional radiation non-ideal magnetohydrodynamics simulations and investigate the impact of the Hall effect on the angular momentum evolution in the collapsing cloud cores in which the magnetic field $\magB$ and angular momentum $\angJ$ are misaligned with each other. We find that the Hall effect notably changes the magnetic torques in the pseudo-disk, and strengthens and weakens the magnetic braking in cores with an acute and obtuse relative angles between $\magB$ and $\angJ$, respectively. This suggests that the bimodal evolution of the disk size may occur in early disk evolutionary phase even if $\magB$ and $\angJ$ are randomly distributed. We show that a counter-rotating envelope form in the upper envelope of the pseudo-disk in cloud cores with obtuse relative angles. We also find that a counter-rotating region forms at the midplane of the pseudo-disk in cloud cores with acute relative angles. The former and latter types of counter-rotating envelopes may be associated with the YSOs with a large ($r\sim100$ AU) and small ($r\lesssim10$ AU) disks, respectively.
△ Less
Submitted 4 October, 2017; v1 submitted 14 June, 2017;
originally announced June 2017.
-
First results from BISTRO -- a SCUBA-2 polarimeter survey of the Gould Belt
Authors:
Derek Ward-Thompson,
Kate Pattle,
Pierre Bastien,
Ray S. Furuya,
Woojin Kwon,
Shih-Ping Lai,
Keping Qiu,
David Berry,
Minho Choi,
Simon Coudé,
James Di Francesco,
Thiem Hoang,
Erica Franzmann,
Per Friberg,
Sarah F. Graves,
Jane S. Greaves,
Martin Houde,
Doug Johnstone,
Jason M. Kirk,
Patrick M. Koch,
Jungmi Kwon,
Chang Won Lee,
Di Li,
Brenda C. Matthews,
Joseph C. Mottram
, et al. (89 additional authors not shown)
Abstract:
We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will a…
▽ More
We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.
△ Less
Submitted 27 April, 2017;
originally announced April 2017.
-
Detection of submillimeter-wave [C I] emission in gaseous debris disks of 49 Ceti and Beta Pictoris
Authors:
Aya E. Higuchi,
Aki Sato,
Takashi Tsukagoshi,
Nami Sakai,
Kazunari Iwasaki,
Munetake Momose,
Hiroshi Kobayashi,
Daisuke Ishihara,
Sakae Watanabe,
Hidehiro Kaneda,
Satoshi Yamamoto
Abstract:
We have detected [C I] 3P1-3P0 emissions in the gaseous debris disks of 49 Ceti and Beta Pictoris with the 10 m telescope of the Atacama Submillimeter Telescope Experiment, which is the first detection of such emissions. The line profiles of [C I] are found to resemble those of CO(J=3-2) observed with the same telescope and the Atacama Large Millimeter/submillimeter Array. This result suggests tha…
▽ More
We have detected [C I] 3P1-3P0 emissions in the gaseous debris disks of 49 Ceti and Beta Pictoris with the 10 m telescope of the Atacama Submillimeter Telescope Experiment, which is the first detection of such emissions. The line profiles of [C I] are found to resemble those of CO(J=3-2) observed with the same telescope and the Atacama Large Millimeter/submillimeter Array. This result suggests that atomic carbon (C) coexists with CO in the debris disks, and is likely formed by the photodissociation of CO. Assuming an optically thin [C I] emission with the excitation temperature ranging from 30 to 100 K, the column density of C is evaluated to be (2.2+-0.2)x10^17 and (2.5+-0.7)x10^16 cm^-2 for 49 Ceti and Beta Pictoris, respectively. The C/CO column density ratio is thus derived to be 54+-19 and 69+-42 for 49 Ceti and Beta Pictoris, respectively. These ratios are higher than those of molecular clouds and diffuse clouds by an order of magnitude. The unusually high ratios of C to CO are likely attributed to a lack of H2 molecules needed to reproduce CO molecules efficiently from C. This result implies a small number of H2 molecules in the gas disk; i.e., there is an appreciable contribution of secondary gas from dust grains.
△ Less
Submitted 20 March, 2017;
originally announced March 2017.
-
Three-Term Relations for 3F2(1)
Authors:
Akihito Ebisu,
Katsunori Iwasaki
Abstract:
For the hypergeometric function of unit argument 3F2(1) we prove the existence and uniqueness of three-term relations with arbitrary integer shifts. We show that not only the original 3F2(1) function but also other five functions related to it satisfy one and the same three-term relation. This fact is referred to as simultaneousness. The uniqueness and simultaneousness provide three-term relations…
▽ More
For the hypergeometric function of unit argument 3F2(1) we prove the existence and uniqueness of three-term relations with arbitrary integer shifts. We show that not only the original 3F2(1) function but also other five functions related to it satisfy one and the same three-term relation. This fact is referred to as simultaneousness. The uniqueness and simultaneousness provide three-term relations with a group symmetry of order 72.
△ Less
Submitted 3 September, 2017; v1 submitted 2 April, 2016;
originally announced April 2016.
-
Minimizing Dispersive Errors in Smoothed Particle Magnetohydrodynamics for Strongly Magnetized Medium
Authors:
Kazunari Iwasaki
Abstract:
In this study, we investigate the dispersive properties of smoothed particle magnetohydrodynamics (SPM) in a strongly magnetized medium by using linear analysis. In modern SPM, a correction term proportional to the divergence of the magnetic fields is subtracted from the equation of motion to avoid a numerical instability arising in a strongly magnetized medium. From the linear analysis, it is fou…
▽ More
In this study, we investigate the dispersive properties of smoothed particle magnetohydrodynamics (SPM) in a strongly magnetized medium by using linear analysis. In modern SPM, a correction term proportional to the divergence of the magnetic fields is subtracted from the equation of motion to avoid a numerical instability arising in a strongly magnetized medium. From the linear analysis, it is found that SPM with the correction term suffer from significant dispersive errors, especially for slow waves propagating along magnetic fields. The phase velocity for all wave numbers is significantly larger than the exact solution and has a peak at a finite wavenumber. These excessively large dispersive errors occur because magnetic fields contribute an unphysical repulsive force along magnetic fields. The dispersive errors cannot be reduced, even with a larger smoothing length and smoother kernel functions such as the Gaussian or quintic spline kernels. We perform the linear analysis for this problem and find that the dispersive errors can be removed completely while keeping SPM stable if the correction term is reduced by half. These findings are confirmed by several simple numerical experiments.
△ Less
Submitted 14 September, 2015;
originally announced September 2015.
-
Bimodality of circumstellar disk evolution induced by Hall current
Authors:
Y. Tsukamoto,
K. Iwasaki,
S. Okuzumi,
M. N. Machida,
S. Inutsuka
Abstract:
The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of the circumstellar disk using three- dimensional simulations. In our simulations, all non-ideal effects as well as the radiat…
▽ More
The formation process of circumstellar disks is still controversial because of the interplay of complex physical processes that occurs during the gravitational collapse of prestellar cores. In this study, we investigate the effect of the Hall current term on the formation of the circumstellar disk using three- dimensional simulations. In our simulations, all non-ideal effects as well as the radiation transfer are considered. The size of the disk is significantly affected by a simple difference in the inherent properties of the prestellar core, namely whether the rotation vector and the magnetic field are parallel or anti- parallel. In the former case, only a very small disk (< 1 AU) is formed. On the other hand, in the latter case, a massive and large (> 20 AU) disk is formed in the early phase of protostar formation. Since the parallel and anti-parallel properties do not readily change, we expect that the parallel and anti-parallel properties are also important in the subsequent disk evolution and the difference between the two cases is maintained or enhanced. This result suggests that the disk size distribution of the Class 0 young stellar objects is bimodal. Thus, the disk evolution can be categorized into two cases and we may call the parallel and anti-parallel systems as Ortho-disk and Para-disk, respectively. We also show that the anti-rotating envelopes against the disk-rotation appear with a size of 200 AU. We predict that the anti-rotating envelope will be found in the future observations.
△ Less
Submitted 18 August, 2015; v1 submitted 24 June, 2015;
originally announced June 2015.
-
The Formation and Destruction of Molecular Clouds and Galactic Star Formation
Authors:
Shu-ichiro Inutsuka,
Tsuyoshi Inoue,
Kazunari Iwasaki,
Takashi Hosokawa
Abstract:
We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubble…
▽ More
We describe an overall picture of galactic-scale star formation. Recent high-resolution magneto-hydrodynamical simulations of two-fluid dynamics with cooling/heating and thermal conduction have shown that the formation of molecular clouds requires multiple episodes of supersonic compression. This finding enables us to create a scenario in which molecular clouds form in interacting shells or bubbles on a galactic scale. First we estimate the ensemble-averaged growth rate of molecular clouds over a timescale larger than a million years. Next we perform radiation hydrodynamics simulations to evaluate the destruction rate of magnetized molecular clouds by the stellar FUV radiation. We also investigate the resultant star formation efficiency within a cloud which amounts to a low value (a few percent) if we adopt the power-law exponent -2.5 for the mass distribution of stars in the cloud. We finally describe the time evolution of the mass function of molecular clouds over a long timescale (>1Myr) and discuss the steady state exponent of the power-law slope in various environments.
△ Less
Submitted 18 May, 2015;
originally announced May 2015.
-
On an Orbifold Hamiltonian Structure for the First Painleve Equation
Authors:
Katsunori Iwasaki,
Shu Okada
Abstract:
For the first Painleve equation we establish an orbifold polynomial Hamiltonian structure on the fibration of Okamoto's spaces and show that this geometric structure uniquely recovers the original Painleve equation, thereby solving a problem posed by K. Takano.
For the first Painleve equation we establish an orbifold polynomial Hamiltonian structure on the fibration of Okamoto's spaces and show that this geometric structure uniquely recovers the original Painleve equation, thereby solving a problem posed by K. Takano.
△ Less
Submitted 15 April, 2015;
originally announced April 2015.