BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Tomoyuki Hanawa

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.EPastro-phastro-ph.HEastro-ph.IM

Catalog footprint

What is connected

20works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

20 published item(s)

preprint2025arXiv

FAUST XXX: Dust enhancement in the young binary L1551 IRS 5

Young binary stars with discs provide unique laboratories to study the earliest stages of planet formation in star-forming environments. The detection of substructure in discs around Class I protostars challenges current models of disc evolution, suggesting that planets may form earlier than previously expected ($<1$ Myr). In the context of the ALMA Large Program FAUST, we present observations of the circumbinary disc (CBD) around the young binary system L1551 IRS 5. The CBD exhibits two prominent over-densities in the continuum emission at the edge of the cavity, with the Northern over-density being about 20% brighter than the Southern one. By analysing the disc morphology and kinematics of L1551 IRS 5, we delineate dynamical constraints on the binary's orbital parameters. Additionally, we present 3D hydrodynamical models of the CBD to predict both the dust and the gas surface densities. Then, we compare the resulting synthetic observations with ALMA observations of the continuum emission at 1.3 mm and the C$^{18}$O line emission. Our analysis suggests that the density enhancements observed with ALMA in L1551 IRS 5 can be caused by interactions between the binary stars and the CBD, leading to dust concentration within the disc. We conclude that the observed over-density corresponds to a location where could potentially grow under favourable conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Cloudlet Capture Model for Asymmetric Molecular Emission Lines Observed in TMC-1A with ALMA

TMC-1A is a protostellar source harboring a young protostar, IRAS 04365+2353, and shows a highly asymmetric features of a few 100 au scale in the molecular emission lines. Blue-shifted emission is much stronger in the CS ($J=5$-4) line than red-shifted one. The asymmetry can be explained if the gas accretion is episodic and takes the form of cloudlet capture, given the cloudlet approached toward us. The gravity of the protostar transforms the cloudlet into a stream and changes its velocity along the flow. The emission from the cloudlet should be blue-shifted before the periastron, while it should be red-shifted after the periastron. If a major part of cloudlet has not reached the periastron, the former should be dominant. We perform hydrodynamical simulations to examine the validity of the scenario. Our numerical simulations can reproduce the observed asymmetry if the orbit of the cloudlet is inclined to the disk plane. The inclination can explain the slow infall velocity observed in the C$^{18}$O ($J$=2-1) line emission. Such episodic accretion may occur in various protostellar cores since actual clouds could have inhomogeneous density distribution. We also discuss the implication of the cloudlet capture on observations of related objects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

FAUST III. Misaligned rotations of the envelope, outflow, and disks in the multiple protostellar system of VLA 1623$-$2417

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Formation of dust clumps with sub-Jupiter mass and cold shadowed region in gravitationally unstable disk around Class 0/I protostar in L1527 IRS

We have investigated the protostellar disk around a Class 0/I protostar, L1527 IRS, using multi-wavelength observations of the dust continuum emission at $λ=0.87$, 2.1, 3.3, and 6.8 mm obtained by the Atacama Large Millimeter/submillimeter Array (ALMA) and the Jansky Very Large Array (VLA). Our observations achieved a spatial resolution of $3-13$ au and revealed an edge-on disk structure with a size of $\sim80-100$ au. The emission at 0.87 and 2.1 mm is found to be optically thick within a projected disk radius of $ r_{\rm proj}\lesssim50$ au. The emission at 3.3 and 6.8 mm shows that the power-law index of the dust opacity ($β$) is $β\sim1.7$ around $ r_{\rm proj}\sim 50$ au, suggesting that grain growth has not yet begun. The dust temperature ($T_{\rm dust}$) shows a steep decrease with $T_{\rm dust}\propto r_{\rm proj}^{-2}$ outside of the VLA clumps previously identified at $r_{\rm proj}\sim20$ au. Furthermore, the disk is gravitationally unstable at $r_{\rm proj}\sim20$ au, as indicated by a Toomre {\it Q} parameter value of $Q\lesssim1.0$. These results suggest that the VLA clumps are formed via gravitational instability, which creates a shadow on the outside of the substructure, resulting in the sudden drop in temperature. The derived dust masses for the VLA clumps are $\gtrsim0.1$ $M_{\rm J}$. Thus, we suggest that Class 0/I disks can be massive enough to be gravitationally unstable, which might be the origin of gas-giant planets in a 20 au radius. Furthermore, the protostellar disks can be cold due to shadowing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

An Extension of the Athena++ Framework for Fully Conservative Self-Gravitating Hydrodynamics

Numerical simulations of self-gravitating flows evolve a momentum equation and an energy equation that account for accelerations and gravitational energy releases due to a time-dependent gravitational potential. In this work, we implement a fully conservative numerical algorithm for self-gravitating flows, using source terms, in the astrophysical magnetohydrodynamics framework Athena++. We demonstrate that properly evaluated source terms are conservative when they are equivalent to the divergence of a corresponding "gravity flux" (i.e., a gravitational stress tensor or a gravitational energy flux). We provide test problems that demonstrate several advantages of the source-term-based algorithm, including second order convergence and round-off error total momentum and total energy conservation. The fully conservative scheme suppresses anomalous accelerations that arise when applying a common numerical discretization of the gravitational stress tensor that does not guarantee curl-free gravity.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Substructure Formation in a Protostellar Disk of L1527 IRS

We analyze multi-frequency, high-resolution continuum data obtained by ALMA and JVLA to study detailed structure of the dust distribution in the infant disk of a Class~0/I source, L1527 IRS. We find three clumps aligning in the north-south direction in the $7 {\rm \,mm}$ radio continuum image. The three clumps remain even after subtracting free-free contamination, which is estimated from the $1.3{\rm \,cm}$ continuum observations. The northern and southern clumps are located at a distance of $\sim 15{\rm \,au}$ from the central clump and are likely optically thick at $7{\rm \,mm}$ wavelength. The clumps have similar integrated intensities. The symmetric physical properties could be realized when a dust ring or spiral arms around the central protostar is projected to the plane of the sky. We demonstrates for the first time that such substructure may form even in the disk-forming stage, where the surrounding materials actively accrete toward a disk-protostar system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Approximate Riemann Solvers for the Cosmic Ray Magnetohydrodynamical Equations

We analyze the cosmic-ray magnetohydrodynamic (CR MHD) equations to improve the numerical simulations. We propose to solve them in the fully conservation form, which is equivalent to the conventional CR MHD equations. In the fully conservation form, the CR energy equation is replaced with the CR "number" conservation, where the CR number density is defined as the three fourths power of the CR energy density. The former contains an extra source term, while latter does not. An approximate Riemann solver is derived from the CR MHD equations in the fully conservation form. Based on the analysis, we propose a numerical scheme of which solutions satisfy the Rankine-Hugoniot relation at any shock. We demonstrate that it reproduces the Riemann solution derived by Pfrommer et al. (2006) for a 1D CR hydrodynamic shock tube problem. We compare the solution with those obtained by solving the CR energy equation. The latter solutions deviate from the Riemann solution seriously, when the CR pressure dominates over the gas pressure in the post-shocked gas. The former solutions converge to the Riemann solution and are of the second order accuracy in space and time. Our numerical examples include an expansion of high pressure sphere in an magnetized medium. Fast and slow shocks are sharply resolved in the example. We also discuss possible extension of the CR MHD equations to evaluate the average CR energy.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Formation of the Unequal-Mass Binary Protostars in L1551 NE by Rotationally-Driven Fragmentation

We present observations at 7 mm that fully resolve the two circumstellar disks, and a reanalyses of archival observations at 3.5 cm that resolve along their major axes the two ionized jets, of the class I binary protostellar system L1551 NE. We show that the two circumstellar disks are better fit by a shallow inner and steep outer power-law than a truncated power-law. The two disks have very different transition radii between their inner and outer regions of $\sim$18.6 AU and $\sim$8.9 AU respectively. Assuming that they are intrinsically circular and geometrically thin, we find that the two circumstellar disks are parallel with each other and orthogonal in projection to their respective ionized jets. Furthermore, the two disks are closely aligned if not parallel with their circumbinary disk. Over an interval of $\sim$10 yr, source B (possessing the circumsecondary disk) has moved northwards with respect to and likely away from source A, indicating an orbital motion in the same direction as the rotational motion of their circumbinary disk. All the aforementioned elements therefore share the same axis for their angular momentum, indicating that L1551 NE is a product of rotationally-driven fragmentation of its parental core. Assuming a circular orbit, the relative disk sizes are compatible with theoretical predictions for tidal truncation by a binary system having a mass ratio of $\sim$0.2, in agreement with the reported relative separations of the two protostars from the center of their circumbinary disk. The transition radii of both disks, however, are a factor of $\gtrsim$1.5 smaller than their predicted tidally-truncated radii.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Rotationally-Driven Fragmentation for the Formation of the Binary Protostellar System L1551 IRS 5

Either bulk rotation or local turbulence is widely invoked to drive fragmentation in collapsing cores so as to produce multiple star systems. Even when the two mechanisms predict different manners in which the stellar spins and orbits are aligned, subsequent internal or external interactions can drive multiple systems towards or away from alignment thus masking their formation process. Here, we demonstrate that the geometrical and dynamical relationship between the binary system and its surrounding bulk envelope provide the crucial distinction between fragmentation models. We find that the circumstellar disks of the binary protostellar system L1551 IRS 5 are closely parallel not just with each other but also with their surrounding flattened envelope. Measurements of the relative proper motion of the binary components spanning nearly 30 yr indicate an orbital motion in the same sense as the envelope rotation. Eliminating orbital solutions whereby the circumstellar disks would be tidally truncated to sizes smaller than are observed, the remaining solutions favor a circular or low-eccentricity orbit tilted by up to $\sim$25$^\circ$ from the circumstellar disks. Turbulence-driven fragmentation can generate local angular momentum to produce a coplanar binary system, but which bears no particular relationship with its surrounding envelope. Instead, the observed properties conform with predictions for rotationally-driven fragmentation. If the fragments were produced at different heights or on opposite sides of the midplane in the flattened central region of a rotating core, the resulting protostars would then exhibit circumstellar disks parallel with the surrounding envelope but tilted from the orbital plane as is observed.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Submillimeter polarization observation of the protoplanetary disk around HD 142527

We present the polarization observations toward the circumstellar disk around HD 142527 by using Atacama Large Millimeter/submillimeter Array (ALMA) at the frequency of 343 GHz. The beam size is $0.51 " \times 0.44 "$, which corresponds to the spatial resolution of $\sim$ 71 $\times$ 62 AU. The polarized intensity displays a ring-like structure with a peak located on the east side with a polarization fraction of $P= 3.26 \pm 0.02$ %, which is different from the peak of the continuum emission from the northeast region. The polarized intensity is significantly weaker at the peak of the continuum where $P= 0.220 \pm 0.010$ %. The polarization vectors are in the radial direction in the main ring of the polarized intensity, while there are two regions outside at the northwest and northeast areas where the vectors are in the azimuthal direction. If the polarization vectors represent the magnetic field morphology, the polarization vectors indicate the toroidal magnetic field configuration on the main ring and the poloidal fields outside. On the other hand, the flip of the polarization vectors is predicted by the self-scattering of thermal dust emission due to the change of the direction of thermal radiation flux. Therefore, we conclude that self-scattering of thermal dust emission plays a major role in producing polarization at millimeter wavelengths in this protoplanetary disk. Also, this puts a constraint on the maximum grain size to be approximately 150 ${\rm μm}$ if we assume compact spherical dust grains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Millimeter-wave polarization of protoplanetary disks due to dust scattering

We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~μm}$ and the observing wavelength to be 870 ${\rm μm}$. We find that the polarization degree is as high as 2.5 % with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarization degree is the highest if dust grains have a maximum size of $a_{\rm max}\simλ/2π$, where $λ$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Significant Gas-to-Dust Ratio Asymmetry and Variation in the Disk of HD 142527 and the Indication of Gas Depletion

We investigate the dust and gas distribution in the disk around HD 142527 based on ALMA observations of dust continuum, 13CO(3-2), and C18O(3-2) emission. The disk shows strong azimuthal asymmetry in the dust continuum emission, while gas emission is more symmetric. In this paper, we investigate how gas and dust are distributed in the dust-bright northern part of the disk and in the dust-faint southern part. We construct two axisymmetric disk models. One reproduces the radial profiles of the continuum and the velocity moments 0 and 1 of CO lines in the north and the other reproduces those in the south. We have found that the dust is concentrated in a narrow ring having ~50AU width (in FWHM; w_d=30AU in our parameter definition) located at ~170-200AU from the central star. The dust particles are strongly concentrated in the north. We have found that the dust surface density contrast between the north and south amounts to ~70. Compared to the dust, the gas distribution is more extended in the radial direction. We find that the gas component extends at least from ~100AU to ~250AU from the central star, and there should also be tenuous gas remaining inside and outside of these radii. The azimuthal asymmetry of gas distribution is much smaller than dust. The gas surface density differs only by a factor of ~3-10 between the north and south. Hence, gas-to-dust ratio strongly depends on the location of the disk: ~30 at the location of the peak of dust distribution in the south and ~3 at the location of the peak of dust distribution in the north. Despite large uncertainties, the overall gas-to-dust ratio is inferred to be ~10-30, indicating that the gas depletion may have already been under way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Structure and Stability of Filamentary Clouds Supported by Lateral Magnetic Field

We have constructed two types of analytical models for an isothermal filamentary cloud supported mainly by magnetic tension. The first one describes an isolated cloud while the second considers filamentary clouds spaced periodically. Both the models assume that the filamentary clouds are highly flattened. The former is proved to be the asymptotic limit of the latter in which each filamentary cloud is much thinner than the distance to the neighboring filaments. We show that these models reproduce main features of the 2D equilibrium model of Tomisaka (2014) for filamentary cloud threaded by perpendicular magnetic field. It is also shown that the critical mass to flux ratio is $ M /Φ= (2 π\sqrt{G}) ^{-1} $, where $ M $, $ Φ$ and $ G $ denote the cloud mass, the total magnetic flux of the cloud, and the gravitational constant, respectively. This upper bound coincides with that for an axisymmetric cloud supported by poloidal magnetic fields. We applied the variational principle for studying the Jeans instability of the first model. Our model cloud is unstable against fragmentation as well as the filamentary clouds threaded by longitudinal magnetic field. The fastest growing mode has a wavelength several times longer than the cloud diameter. The second model describes quasi-static evolution of filamentary molecular cloud by ambipolar diffusion.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Angular Momentum Exchange by Gravitational Torques and Infall in the Circumbinary Disk of the Protostellar System L1551 NE

We report the ALMA observation of the Class I binary protostellar system L1551 NE in the 0.9-mm continuum, C18O (3-2), and 13CO (3-2) lines at a ~1.6 times higher resolution and a ~6 times higher sensitivity than those of our previous SMA observations, which revealed a r ~300 AU-scale circumbinary disk in Keplerian rotation. The 0.9-mm continuum shows two opposing U-shaped brightenings in the circumbinary disk, and exhibits a depression between the circumbinary disk and the circumstellar disk of the primary protostar. The molecular lines trace non-axisymmetric deviations from Keplerian rotation in the circumbinary disk at higher velocities relative to the systemic velocity, where our previous SMA observations could not detect the lines. In addition, we detect inward motion along the minor axis of the circumbinary disk. To explain the newly-observed features, we performed a numerical simulation of gas orbits in a Roche potential tailored to the inferred properties of L1551 NE. The observed U-shaped dust features coincide with locations where gravitational torques from the central binary system are predicted to impart angular momentum to the circumbinary disk, producing shocks and hence density enhancements seen as a pair of spiral arms. The observed inward gas motion coincides with locations where angular momentum is predicted to be lowered by the gravitational torques. The good agreement between our observation and model indicates that gravitational torques from the binary stars constitute the primary driver for exchanging angular momentum so as to permit infall through the circumbinary disk of L1551 NE.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Kinetic Scheme for Solving M1 Model of Radiative Transfer

We show a numerical scheme to solve the moment equations of the radiative transfer, i.e., M1 model which follows the evolution of the energy density, $ E $, and the energy flux, $ \mbox{\boldmath$F$} $. In our scheme we reconstruct the intensity from $ E $ and $ \mbox{\boldmath$F$} $ so that it is consistent with the closure relation, relation, $ χ= (3 + 4 f ^2)/(5 + 2 \sqrt{4 - 3 f ^2}) $. Here the symbols, $ χ$, $ f = |\mbox{\boldmath$F$}|/(cE) $, and $ c $, denote the Eddington factor, the reduced flux, and the speed of light, respectively. We evaluate the numerical flux across the cell surface from the kinetically reconstructed intensity. It is an explicit function of $ E $ and $ \mbox{\boldmath$F$} $ in the neighboring cells across the surface considered. We include absorption and reemission within a numerical cell in the evaluation of the numerical flux. The numerical flux approaches to the diffusion approximation when the numerical cell itself is optically thick. Our numerical flux gives a stable solution even when some regions computed are very optically thick. We show the advantages of the numerical flux with examples. They include flash of beamed photons and irradiated protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Local Enhancement of Surface Density in the Protoplanetary Ring Surrounding HD 142527

We report ALMA observations of dust continuum, 13CO J=3--2, and C18O J=3--2 line emission toward a gapped protoplanetary disk around HD 142527. The outer horseshoe-shaped disk shows the strong azimuthal asymmetry in dust continuum with the contrast of about 30 at 336 GHz between the northern peak and the southwestern minimum. In addition, the maximum brightness temperature of 24 K at its northern area is exceptionally high at 160 AU from a star. To evaluate the surface density in this region, the grain temperature needs to be constrained and was estimated from the optically thick 13CO J=3--2 emission. The lower limit of the peak surface density was then calculated to be 28 g cm-2 by assuming a canonical gas-to-dust mass ratio of 100. This finding implies that the region is locally too massive to withstand self-gravity since Toomre's Q <~1--2, and thus, it may collapse into a gaseous protoplanet. Another possibility is that the gas mass is low enough to be gravitationally stable and only dust grains are accumulated. In this case, lower gas-to-dust ratio by at least 1 order of magnitude is required, implying possible formation of a rocky planetary core.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Direct Imaging of Bridged Twin Protoplanetary Disks in a Young Multiple Star

Studies of the structure and evolution of protoplanetary disks are important for understanding star and planet formation. Here, we present the direct image of an interacting binary protoplanetary system. Both circumprimary and circumsecondary disks are resolved in the near-infrared. There is a bridge of infrared emission connecting the two disks and a long spiral arm extending from the circumprimary disk. Numerical simulations show that the bridge corresponds to gas flow and a shock wave caused by the collision of gas rotating around the primary and secondary stars. Fresh material streams along the spiral arm, consistent with the theoretical scenarios where gas is replenished from a circummultiple reservoir.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Protostellar collapse of magneto-turbulent cloud cores: shape during collapse and outflow formation

We investigate protostellar collapse of molecular cloud cores by numerical simulations, taking into account turbulence and magnetic fields. By using the adaptive mesh refinement technique, the collapse is followed over a wide dynamic range from the scale of a turbulent cloud core to that of the first core. The cloud core is lumpy in the low density region owing to the turbulence, while it has a smooth density distribution in the dense region produced by the collapse. The shape of the dense region depends mainly on the mass of the cloud core; a massive cloud core tends to be prolate while a less massive cloud core tends to be oblate. In both cases, anisotropy of the dense region increases during the isothermal collapse. The minor axis of the dense region is always oriented parallel to the local magnetic field. All the models eventually yield spherical first cores supported mainly by the thermal pressure. Most of turbulent cloud cores exhibit protostellar outflows around the first cores. These outflows are classified into two types, bipolar and spiral flows, according to the morphology of the associated magnetic field. Bipolar flow often appears in the less massive cloud core. The rotation axis of the first core is oriented parallel to the local magnetic field for bipolar flow, while the orientation of the rotation axis from the global magnetic field depends on the magnetic field strength. In spiral flow, the rotation axis is not aligned with the local magnetic field.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Gas Accretion from a Circumbinary Disk

A new computational scheme is developed to study gas accretion from a circumbinary disk. The scheme decomposes the gas velocity into two components one of which denotes the Keplerian rotation and the other of which does the deviation from it. This scheme enables us to solve the centrifugal balance of a gas disk against gravity with better accuracy, since the former inertia force cancels the gravity. It is applied to circumbinary disk rotating around binary of which primary and secondary has mass ratio, 1.4:0.95. The gravity is reduced artificially softened only in small circular regions around the primary and secondary. The radii are 7% of the binary separation and much smaller than those in the previous grid based simulations. 7 Models are constructed to study dependence on the gas temperature and the initial inner radius of the disk. The gas accretion shows both fast and slow time variations while the binary is assumed to have a circular orbit. The time variation is due to oscillation of spiral arms in the circumbinary disk. The masses of primary and secondary disks increase while oscillating appreciably. The mass accretion rate tends to be higher for the primary disk although the secondary disk has a higher accretion rate in certain periods. The primary disk is perturbed intensely by the impact of gas flow so that the outer part is removed. The secondary disk is quiet in most of time on the contrary. Both the primary and secondary disks have traveling spiral waves which transfer angular momentum within them.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Stability of Dynamically Collapsing Gas Sphere

We discuss stability of dynamically collapsing gas spheres. We use a similarity solution for a dynamically collapsing sphere as the unperturbed state. In the similarity solution the gas pressure is approximated by a polytrope of $ P = K ρ^γ$. We examine three types of perturbations: bar ($ \ell = 2$) mode, spin-up mode, and Ori-Piran mode. When $ γ< 1.097 $, it is unstable against bar-mode. It is unstable against spin-up mode for any $ γ$. When $ γ< 0.961 $, the similarity solution is unstable against Ori-Piran mode. The unstable mode grows in proportion to $ | t - t_0 | ^{-σ} $ while the central density increases in proportion to $ ρ_c \propto (t - t_0) ^{-2} $ in the similarity solution. The growth rate, $ σ$ is obtained numerically as a function of $ γ$ for bar mode and Ori-Piran mode. The growth rate of the bar mode is larger for a smaller $ γ$. The spin-up mode has the growth rate of $ σ= 1/3 $ for any $ γ$.

0 citations0 reviewsNo code linkedOpen access