BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Satoshi Okuzumi

Satoshi Okuzumi contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.EPastro-ph.SRastro-ph.GA

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
3topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 published item(s)

preprint2026arXiv

A Dichotomy of the Mass-Metallicity Relation of Exoplanetary Atmospheres Demarcated by their Birthplace

Atmospheric observations by JWST raise a growing evidence that atmospheric metallicity exhibits an anti-correlation with masses of giant exoplanets. While such a trend was anticipated by planetesimal-based planet formation models, it remains unclear what kind of atmospheric metallicity trends emerge from pebble-based planet formation. Moreover, while recent studies of solar system Jupiter suggest that uppermost observable atmosphere may not represent the bulk envelope composition, it remains uncertain how the envelope inhomogeneity influences the atmospheric metallicity trend. In this study, we develop disk evolution and planet formation models to investigate the possible atmospheric metallicity trends of giant exoplanets formed via pebble accretion and how they depend on the metallicity inhomogeneity within the envelope. We find that pebble-based planet formation produces two distinct mass-metallicity relations depending on planetary birthplace. Planets formed beyond the H2O snowline exhibit a mass-metallicity anti-correlation similar to that predicted by planetesimal-based models if their atmospheres are fully convective. This anti-correlation disappears if the convective mixing is inefficient. In contrast, planets formed inside the H2O snowline show a shallower mass-metallicity anti-correlation, regardless of the efficiency of atmospheric mixing. Many gas giants observed by JWST observations lie around the mass-metallicity relation predicted for formation at close-in orbits, although some planets with sub-stellar atmospheric metallicity appear to require unmixed envelopes and formation beyond the H2O snowline. We also examine the relationship between bulk and atmospheric metallicity and find a clear correlation that closely follows atmospheric metallicity that is comparable to bulk metallicity.

0 citations0 reviewsNo code linkedOpen access
preprint2026arXiv

Bridging the gap: consistent modeling of protoplanetary disk heating and gap formation by planet-induced spiral shocks

A giant planet embedded in a protoplanetary disk excites spiral density waves, which steepen into shocks as they propagate away from the planet. These shocks lead to secular disk heating and gap opening, both of which can have important implications for the evolution of solids near the planet. To date, these two effects have largely been modeled independently. In this study, we present a self-consistent model that unifies these processes by linking shock heating and angular momentum deposition through the entropy jumps across the spiral shocks. We show that this model accurately reproduces the temperature and surface density profiles around the planet's orbit, as obtained from two-dimensional hydrodynamic simulations with standard $α$ viscosity and $β$ thermal relaxation prescriptions. Furthermore, by incorporating an empirically derived scaling law for the radial distribution of the entropy jump, we construct a fully analytic model that self-consistently predicts the temperature and surface density structures of disks hosting a giant planet. This work represents a first step toward understanding how a giant planet forming in the inner disk region influences the distribution and composition of second-generation planets and planetesimals in its vicinity.

0 citations0 reviewsNo code linkedOpen access
preprint2023arXiv

Structured Distributions of Gas and Solids in Protoplanetary Disks

Recent spatially-resolved observations of protoplanetary disks revealed a plethora of substructures, including concentric rings and gaps, inner cavities, misalignments, spiral arms, and azimuthal asymmetries. This is the major breakthrough in studies of protoplanetary disks since Protostars and Planets VI and is reshaping the field of planet formation. However, while the capability of imaging substructures in protoplanetary disks has been steadily improving, the origin of many substructures are still largely debated. The structured distributions of gas and solids in protoplanetary disks likely reflect the outcome of physical processes at work, including the formation of planets. Yet, the diverse properties among the observed protoplanetary disk population, for example, the number and radial location of rings and gaps in the dust distribution, suggest that the controlling process may differ between disks and/or the outcome may be sensitive to stellar or disk properties. In this review, we (1) summarize the existing observations of protoplanetary disk substructures collected from the literature; (2) provide a comprehensive theoretical review of various processes proposed to explain observed protoplanetary disk substructures; (3) compare current theoretical predictions with existing observations and highlight future research directions to distinguish between different origins; and (4) discuss implications of state-of-the-art protoplanetary disk observations to protoplanetary disk and planet formation theory.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

A global two-layer radiative transfer model for axisymmetric, shadowed protoplanetary disks

Understanding the thermal structure of protoplanetary disks is crucial for modeling planet formation and interpreting disk observations. We present a new two-layer radiative transfer model for computing the thermal structure of axisymmetric irradiated disks. Unlike the standard two-layer model, our model accounts for the radial as well as vertical transfer of the starlight reprocessed at the disk surface. The model thus allows us to compute the temperature below "shadowed" surfaces receiving no direct starlight. Thanks to the assumed axisymmetry, the reprocessed starlight flux is given in one-dimensional integral form that can be computed at a low cost. Furthermore, our model evolves the midplane temperature using a time-dependent energy equation and can therefore treat thermal instabilities. We apply our global two-layer model to disks with a planetary induced gap and confirm that the model reproduces the disks' temperature profiles obtained from more computationally expensive Monte Carlo radiative transfer calculations to an accuracy of less than 20%. We also apply the model to study the long-term behavior of the thermal wave instability in irradiated disks. Being simple and computationally efficient, the global two-layer model will be suitable for studying the interplay between disks' thermal evolution and dust evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

ALMA High-resolution Multiband Analysis for the Protoplanetary Disk around TW Hya

We present a high-resolution (2.5 au) multiband analysis of the protoplanetary disk around TW Hya using ALMA long baseline data at Bands 3, 4, 6, and 7. We aim to reconstruct a high-sensitivity millimeter continuum image and revisit the spectral index distribution. The imaging is performed by combining new ALMA data at Bands 4 and 6 with available archive data. Two methods are employed to reconstruct the images; multi-frequency synthesis (MFS) and the fiducial image-oriented method, where each band is imaged separately and the frequency dependence is fitted pixel by pixel. We find that the MFS imaging with the second order of Taylor expansion can reproduce the frequency dependence of the continuum emission between Bands 3 and 7 in a manner consistent with previous studies and is a reasonable method to reconstruct the spectral index map. The image-oriented method provides a spectral index map consistent with the MFS imaging, but with a two times lower resolution. Mock observations of an intensity model were conducted to validate the images from the two methods. We find that the MFS imaging provides a high-resolution spectral index distribution with an uncertainty of $<10$~\%. Using the submillimeter spectrum reproduced from our MFS images, we directly calculated the optical depth, power-law index of the dust opacity coefficient ($β$), and dust temperature. The derived parameters are consistent with previous works, and the enhancement of $β$ within the intensity gaps is also confirmed, supporting a deficit of millimeter-sized grains within the gaps.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Color Dependence of the Transit Detectability for Young Active M-dwarfs

We investigate the planetary transit detectability in the presence of stellar rotational activity from light curves for young M-dwarfs and estimate improvements of the detection at near-infrared (NIR) wavelengths. Making maps of the transit signal detection efficiency over the orbital period and planetary radius with light curves of members of four clusters, Hyades, Praesepe, Pleiades, and Upper Scorpius observed by the K2 mission, we evaluate the detectability for the rotation period and modulation semi-amplitude. We find that the detection efficiency remarkably decreases to about 20% for rapidly rotators with P_{rot} <= 1 d and the lack of planets in Pleiades is likely due to the high fraction of rapidly rotating M-dwarfs. We also evaluate the improvements of the planet detection with NIR photometry via tests using mock light curves assuming that the signal amplitude of stellar rotation decreases at NIR wavelengths. Our results suggest that NIR photometric monitoring would double relative detection efficiency for transiting planetary candidates with P_{rot} <= 1 d and find planets around M-dwarfs with approximately 100 Myr missing in the past transit surveys from the space.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Impact of dust size distribution including large dust grains on magnetic resistivity: an analytical approach

This paper investigates the impact of dust size distribution on magnetic resistivity. In particular, we focus on its impact when the maximum dust size significantly increases from sub-micron. The first half of the paper describes our calculation method for magnetic resistivity based on the model of \citet{1987ApJ...320..803D} and shows that the method reproduces the results of a more realistic chemical reaction network calculations reasonably well. Then, we describe the results of the resistivity calculations for dust distributions with large maximum dust grains. Our results show that resistivity tends to decrease with dust growth, which is particularly true when the dust size power exponent $q$ is $q=2.5$. On the other hand, the decrease is less pronounced when the dust size power exponent $q$ is $q=3.5$, i.e., when the small dust is also responsible for the dust cross-section. Our results suggest that detailed dust coagulation and fragmentation processes play a vital role in the magnetic resistivities in protostar formation.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Planetesimal formation around the snow line. II. Dust or pebbles?

Around the snow line, icy pebbles and silicate dust may locally pile-up and form icy and rocky planetesimals via streaming instability and/or gravitational instability. We perform 1D diffusion-advection simulations that include the back-reaction to radial drift and diffusion of icy pebbles and silicate dust, ice sublimation, release of silicate dust, and their recycling through recondensation and sticking onto pebbles outside the snow line. We use a realistic description of the scale height of silicate dust obtained from Ida et al. and that of pebbles including the effects of a Kelvin-Helmholtz instability. We study the dependence of solid pile-up on distinct effective viscous parameters for turbulent diffusions in the radial and vertical directions ($α_{\rm Dr}$ and $α_{\rm Dz}$) and for the gas accretion to the star ($α_{\rm acc}$) as well as that on the pebble-to-gas mass flux ($F_{\rm p/g}$). We derive the sublimation width of drifting icy pebbles which is a critical parameter to characterize the pile-up of silicate dust and pebbles around the snow line. We identify a parameter space (in the $F_{\rm p/g}-α_{\rm acc}-α_{\rm Dz}(=α_{\rm Dr})$ space) where pebbles no longer drift inward to reach the snow line due to the back-reaction that slows down radial velocity of pebbles. We show that the pile-up of solids around the snow line occurs in a broader range of parameters for $α_{\rm acc}=10^{-3}$ than for $α_{\rm acc}=10^{-2}$. Above a critical $F_{\rm p/g}$ value, the runaway pile-up of silicate dust inside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \ll 1$, while that of pebbles outside the snow line is favored for $α_{\rm Dr}/α_{\rm acc} \sim 1$. Our results imply that a distinct evolutionary path could produce a diversity of outcomes in terms of planetesimal formation around the snow line.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Steady-state accretion in magnetized protoplanetary disks

[abridged] We present a 1+1D global magnetically-driven disk accretion model that captures the essence of the MRI-driven accretion, without resorting to 3D global non-ideal MHD simulations. The gas dynamics is assumed to be solely controlled by the MRI and hydrodynamic instabilities. For given stellar and disk parameters, the Shakura-Sunyaev viscosity parameter $α$ is determined self-consistently under the framework of viscously-driven accretion from detailed considerations of the MRI with non-ideal MHD effects (Ohmic resistivity and ambipolar diffusion), accounting for disk heating by stellar irradiation, non-thermal sources of ionization, and dust effects on the ionization chemistry. Additionally, the magnetic field strength is constrained and adopted to maximize the MRI activity. We demonstrate the use of our framework by investigating steady-state MRI-driven accretion in a fiducial protoplanetary disk model around a solar-type star. We find that the equilibrium solution displays no pressure maximum at the dead zone outer edge, except if a sufficient amount of dust particles have accumulated there before the disk reaches a steady-state accretion regime. Furthermore, the steady-state accretion solution describes a disk that displays a spatially extended long-lived inner disk gas reservoir (the dead zone) accreting a few $10^{-9}\, M_{\odot}.\rm{yr}^{-1}$. By conducting a detailed parameter study, we find that the extend to which the MRI can drive efficient accretion is primarily determined by the total disk gas mass, the representative grain size, the vertically-integrated dust-to-gas mass ratio, and the stellar X-ray luminosity. A self-consistent time-dependent coupling between gas, dust, stellar evolution models and our general framework on million-year timescales is required to fully understand the formation of dead zones and their potential to trap dust particles.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Clouds of Fluffy Aggregates: How They Form in Exoplanetary Atmospheres and Influence Transmission Spectra

Transmission spectrum surveys have suggested the ubiquity of high-altitude clouds in exoplanetary atmospheres. Theoretical studies have investigated the formation processes of the high-altitude clouds; however, cloud particles have been commonly approximated as compact spheres, which is not always true for solid mineral particles that likely constitute exoplanetary clouds. Here, we investigate how the porosity of cloud particles evolve in exoplanetary atmospheres and influence the cloud vertical profiles. We first construct a porosity evolution model that takes into account the fractal aggregation and the compression of cloud particle aggregates. Using a cloud microphysical model coupled with the porosity model, we demonstrate that the particle internal density can significantly decrease during the cloud formation. As a result, fluffy-aggregate clouds ascend to altitude much higher than that for compact-sphere clouds assumed so far. We also examine how the fluffy-aggregate clouds affect transmission spectra. We find that the clouds largely obscure the molecular features and produce a spectral slope originated by the scattering properties of aggregates. Finally, we compare the synthetic spectra with the observations of GJ1214 b and find that its flat spectrum could be explained if the atmospheric metallicity is sufficiently high ($\ge100\times$ solar) and the monomer size is sufficiently small ($r_{\rm mon}<1~{\rm μm}$). The high-metallicity atmosphere may offer the clues to explore the gas accretion processes onto past GJ1214b.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planetesimal formation around the snow line: I. Monte Carlo simulations of silicate dust pile-up in a turbulent disk

Context: The formation of rocky planetesimals is a long-standing problem in planet formation theory. One of the possibilities is that it results from gravitational instability as a result of pile-up of small silicate dust particles released from sublimating icy pebbles that pass the snow line. Aims: We want to understand and quantify the role of the water snow line for the formation of rock-rich and ice-rich planetesimals. In this paper, we focus on the formation of rock-rich planetesimals. A companion paper examines the combined formation of both rock-rich and ice-rich planetesimals. Methods: We develop a new Monte Carlo code to calculate the radial evolution of silicate particles in a turbulent accretion disk, accounting for the back-reaction (i.e., inertia) of the particles on their radial drift velocity and diffusion. Results depend in particular on the particle injection width (determined from the radial sublimation width of icy pebbles), the pebble scale height and the pebble mass flux through the disk. The scale height evolution of the silicate particles, which is the most important factor for the runaway pile-up, is automatically calculated in this Lagrange method. Results: From the numerical results, we derive semi-analytical relations for the scale height of the silicate dust particles and the particles-to-gas density ratio at the midplane, as functions of a pebble-to-gas mass flux ratio and the $α$ parameters for disk gas accretion and vertical/radial diffusion. We find that the runaway pile-up of the silicate particles (formation of rocky planetesimals) occurs if the pebble-to-gas mass flux ratio is $> [(α_{Dz}/α_{acc})/3 \times 10^{-2}]^{1/2}$ where $α_{Dz}$ and $α_{acc}$ are the $α$ parameters for vertical turbulent diffusion and disk gas accretion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Galilean Satellites Formed Slowly from Pebbles

It is generally accepted that the four major (Galilean) satellites formed out of the gas disk that accompanied Jupiter&#39;s formation. However, understanding the specifics of the formation process is challenging as both small particles (pebbles) as well as the satellites are subject to fast migration processes. Here, we hypothesize a new scenario for the origin of the Galilean system, based on the capture of several planetesimal seeds and subsequent slow accretion of pebbles. To halt migration, we invoke an inner disk truncation radius, and other parameters are tuned for the model to match physical, dynamical, compositional, and structural constraints. In our scenario it is natural that Ganymede&#39;s mass is determined by pebble isolation. Our slow-pebble-accretion scenario then reproduces the following characteristics: (1) the mass of all the Galilean satellites; (2) the orbits of Io, Europa, and Ganymede captured in mutual 2:1 mean motion resonances; (3) the ice mass fractions of all the Galilean satellites; (4) the unique ice-rock partially differentiated Callisto and the complete differentiation of the other satellites. Our scenario is unique to simultaneously reproduce these disparate properties.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization

Dust coagulation in a protoplanetary disk is the first step of planetesimal formation. However, a pathway from dust aggregates to planetesimals remains unclear. Both numerical simulations and laboratory experiments have suggested the importance of dust structure in planetesimal formation, but it is not well constrained by observations. We study how dust structure and porosity alters polarimetric images at millimeter wavelength by performing 3D radiative transfer simulations. Aggregates with different porosity and fractal dimension are considered. As a result, we find that dust aggregates with lower porosity and/or higher fractal dimension are favorable to explain observed millimeter-wave scattering polarization of disks. Aggregates with extremely high porosity fail to explain the observations. In addition, we also show that particles with moderate porosity show weak wavelength dependence of scattering polarization, indicating that multi-wavelength polarimetry is useful to constrain dust porosity. Finally, we discuss implications for dust evolution and planetesimal formation in disks.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Opacity of fluffy dust aggregates

Context. Dust grains coagulate to form dust aggregates in protoplanetary disks. Their porosity can be extremely high in the disks. Although disk emission may come from fluffy dust aggregates, the emission has been modeled with compact grains. Aims. We aim to reveal the mass opacity of fluffy aggregates from infrared to millimeter wavelengths with the filling factor ranging from 1 down to $10^{-4}$. Methods. We use Mie calculations with an effective medium theory. The monomers are assumed to be 0.1 ${\rm μm}$ sized grains, which is much shorter than the wavelengths that we focus on. Results. We find that the absorption mass opacity of fluffy aggregates are characterized by the product $a\times f$, where $a$ is the dust radius and $f$ is the filling factor, except for the interference structure. The scattering mass opacity is also characterized by $af$ at short wavelengths while it is higher in more fluffy aggregates at long wavelengths. We also derive the analytic formula of the mass opacity and find that it reproduces the Mie calculations. We also calculate the expected difference of the emission between compact and fluffy aggregates in protoplanetary disks with a simple dust growth and drift model. We find that compact grains and fluffy aggregates can be distinguished by the radial distribution of the opacity index $β$. The previous observation of the radial distribution of $β$ is consistent with the fluffy case, but more observations are required to distinguish between fluffy or compact. In addition, we find that the scattered light would be another way to distinguish between compact grains and fluffy aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Electric charging of dust aggregates and its effect on dust coagulation in protoplanetary disks

Mutual sticking of dust aggregates is the first step toward planetesimal formation in protoplanetary disks. In spite that the electric charging of dust particles is well recognized in some contexts, it has been largely ignored in the current modeling of dust coagulation. In this study, we present a general analysis of the dust charge state in protoplanetary disks, and then demonstrate how the electric charging could dramatically change the currently accepted scenario of dust coagulation. First, we describe a new semianalytical method to calculate the dust charge state and gas ionization state self-consistently. This method is far more efficient than previous numerical methods, and provides a general and clear description of the charge state of gas-dust mixture. Second, we apply this analysis to early evolutionary stages where the dust has been thought to grow into fractal ($D \sim 2$) aggregates with a quasi-monodisperse (i.e., narrow) size distribution. We find that, for a wide range of model parameters, the fractal growth is strongly inhibited by the electric repulsion between colliding aggregates and eventually &#34;freezes out&#34; on its way to the subsequent growth stage involving collisional compression. Strong disk turbulence would help the aggregates to overcome this growth barrier, but then it would cause catastrophic collisional fragmentation in later growth stages. These facts suggest that the combination of electric repulsion and collisional fragmentation would impose a serious limitation on dust growth in protoplanetary disks. We propose a possible scenario of dust evolution after the freeze-out. Finally, we point out that the fractal growth of dust aggregates tends to maintain a low ionization degree and, as a result, a large magnetorotationally stable region in the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates

Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. We point out an important implication of this discovery for dust growth in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access