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

Soko Matsumura

Soko Matsumura contributes to research discovery and scholarly infrastructure.

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

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
15works
0followers
4topics
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

15 published item(s)

preprint2020arXiv

A new and simple prescription for planet orbital migration and eccentricity damping by planet-disc interactions based on dynamical friction

During planet formation gravitational interaction between a planetary embryo and the protoplanetary gas disc causes orbital migration of the planetary embryo, which plays an important role in shaping the final planetary system. While migration sometimes occurs in the supersonic regime, wherein the relative velocity between the planetary embryo and the gas is higher than the sound speed, migration prescriptions proposed thus far describing the planet-disc interaction force and the timescales of orbital change in the supersonic regime are inconsistent with one another. Here we discuss the details of existing prescriptions in the literature and derive a new simple and intuitive formulation for planet-disc interactions based on dynamical friction that can be applied in both supersonic and subsonic cases. While the existing prescriptions assume particular disc models, ours include the explicit dependence on the disc parameters; hence it can be applied to discs with any radial surface density and temperature dependence (except for the local variations with radial scales less than the disc scale height). Our prescription will reduce the uncertainty originating from different literature formulations of planet migration and will be an important tool to study planet accretion processes, especially when studying the formation of close-in low-mass planets that are commonly found in exoplanetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Effects of Dynamical Evolution of Giant Planets on the Delivery of Atmophile Elements During Terrestrial Planet Formation

Recent observations started revealing the compositions of protostellar discs and planets beyond the Solar System. In this paper, we explore how the compositions of terrestrial planets are affected by dynamical evolution of giant planets. We estimate the initial compositions of building blocks of these rocky planets by using a simple condensation model, and numerically study the compositions of planets formed in a few different formation models of the Solar System. We find that the abundances of refractory and moderately volatile elements are nearly independent of formation models, and that all the models could reproduce the abundances of these elements of the Earth. The abundances of atmophile elements, on the other hand, depend on the scattering rate of icy planetesimals into the inner disc as well as the mixing rate of the inner planetesimal disc. For the classical formation model, neither of these mechanisms are efficient and the accretion of atmophile elements during the final assembly of terrestrial planets appears to be difficult. For the Grand Tack model, both of these mechanisms are efficient, which leads to a relatively uniform accretion of atmophile elements in the inner disc. It is also possible to have a "hybrid" scenario where the mixing is not very efficient but the scattering is efficient. The abundances of atmophile elements in this case increases with orbital radii. Such a scenario may occur in some of the extrasolar planetary systems which are not accompanied by giant planets or those without strong perturbations from giants. We also confirm that the Grand Tack scenario leads to the distribution of asteroid analogues where rocky planetesimals tend to exist interior to icy ones, and show that their overall compositions are consistent with S-type and C-type chondrites, respectively.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Deep Thermal Infrared Imaging of HR 8799 bcde: New Atmospheric Constraints and Limits on a Fifth Planet

We present new $L^\prime$ (3.8 $μm$) and Br-$α$ (4.05 $μm$) data and reprocessed archival $L^\prime$ data for the young, planet-hosting star HR 8799 obtained with Keck/NIRC2, VLT/NaCo and Subaru/IRCS. We detect all four HR 8799 planets in each dataset at a moderate to high signal-to-noise (SNR $\gtrsim$ 6-15). We fail to identify a fifth planet, &#34;HR 8799 f&#34;, at $r$ $<$ 15 $AU$ at a 5-$σ$ confidence level: one suggestive, marginally significant residual at 0.2&#34; is most likely a PSF artifact. Assuming companion ages of 30 $Myr$ and the Baraffe (Spiegel \& Burrows) planet cooling models, we rule out an HR 8799 f with mass of 5 $M_{J}$ (7 $M_{J}$), 7 $M_{J}$ (10 $M_{J}$), and 12 $M_{J}$ (13 $M_{J}$) at $r_{proj}$ $\sim$ 12 $AU$, 9 $AU$, and 5 $AU$, respectively. All four HR 8799 planets have red early T dwarf-like $L^\prime$ - [4.05] colors, suggesting that their SEDs peak in between the $L^\prime$ and $M^\prime$ broadband filters. We find no statistically significant difference in HR 8799 cde&#39;s colors. Atmosphere models assuming thick, patchy clouds appear to better match HR 8799 bcde&#39;s photometry than models assuming a uniform cloud layer. While non-equilibrium carbon chemistry is required to explain HR 8799 bc&#39;s photometry/spectra, evidence for it from HR 8799 de&#39;s photometry is weaker. Future, deep IR spectroscopy/spectrophotometry with the Gemini Planet Imager, SCExAO/CHARIS, and other facilities may clarify whether the planets are chemically similar or heterogeneous.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Brazil-nut effect and its application to asteroids

Out of the handful of asteroids that have been imaged, some have distributions of blocks that are not easily explained. In this paper, we investigate the possibility that seismic shaking leads to the size sorting of particles in asteroids. In particular, we focus on the so-called Brazil Nut Effect (BNE) that separates large particles from small ones under vibrations. We study the BNE over a wide range of parameters by using the N-body code PKDGRAV, and find that the effect is largely insensitive to the coefficients of restitution, but sensitive to friction constants and oscillation speeds. Agreeing with the previous results, we find that convection drives the BNE, where the intruder rises to the top of the particle bed. For the wide-cylinder case, we also observe a &#34;whale&#34; effect, where the intruder follows the convective current and does not stay at the surface. We show that the non-dimensional critical conditions for the BNE agree well with previous studies. We also show that the BNE is scalable for low-gravity environments and that the rise speed of an intruder is proportional to the square root of the gravitational acceleration. Finally, we apply the critical conditions to observed asteroids, and find that the critical oscillation speeds are comparable to the seismic oscillation speeds that are expected from non-destructive impacts.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Combined VLT and Gemini Study of the Atmosphere of the Directly-Imaged Planet, beta Pictoris b

We analyze new/archival VLT/NaCo and Gemini/NICI high-contrast imaging of the young, self-luminous planet $β$ Pictoris b in seven near-to-mid IR photometric filters, using advanced image processing methods to achieve high signal-to-noise, high precision measurements. While $β$ Pic b&#39;s near-IR colors mimick that of a standard, cloudy early-to-mid L dwarf, it is overluminous in the mid-infrared compared to the field L/T dwarf sequence. Few substellar/planet-mass objects -- i.e. $κ$ And b and 1RXJ 1609B -- match $β$ Pic b&#39;s $JHK_{s}L^\prime$ photometry, and its 3.1 $μm$ and 5 $μm$ photometry are particularly difficult to reproduce. Atmosphere models adopting cloud prescriptions and large ($\sim$ 60 $μm$) dust grains fail to reproduce the $β$ Pic b spectrum. However, models incorporating thick clouds similar to those found for HR 8799 bcde but also with small (a few microns) modal particle sizes yield fits consistent with the data within uncertainties. Assuming solar abundance models, thick clouds, and small dust particles ($<a>$ = 4 $μm$) we derive atmosphere parameters of log(g) = 3.8 $\pm$ 0.2 and $T_{eff}$ = 1575--1650 $K$, an inferred mass of 7$^{+4}_{-3}$ $M_{J}$, and a luminosity of log(L/L$_{\odot}$) $\sim$ -3.80 $\pm$ 0.02. The best-estimated planet radius, $\approx$ 1.65 $\pm$ 0.06 $R_{J}$, is near the upper end of allowable planet radii for hot-start models given the host star&#39;s age and likely reflects challenges with constructing accurate atmospheric models. Alternatively, these radii are comfortably consistent with hot-start model predictions if $β$ Pic b is younger than $\approx$ 7 Myr, consistent with a late formation, well after its host star&#39;s birth $\sim$ 12$^{+8}_{-4}$ Myr ago.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Effects of Dynamical Evolution of Giant Planets on Survival of Terrestrial Planets

The orbital distributions of currently observed extrasolar giant planets allow marginally stable orbits for hypothetical, terrestrial planets. In this paper, we propose that many of these systems may not have additional planets on these &#34;stable&#34; orbits, since past dynamical instability among giant planets could have removed them. We numerically investigate the effects of early evolution of multiple giant planets on the orbital stability of the inner, sub-Neptune-like planets which are modeled as test particles, and determine their dynamically unstable region. Previous studies have shown that the majority of such test particles are ejected out of the system as a result of close encounters with giant planets. Here, we show that secular perturbations from giant planets can remove test particles at least down to 10 times smaller than their minimum pericenter distance. Our results indicate that, unless the dynamical instability among giant planets is either absent or quiet like planet-planet collisions, most test particles down to ~0.1 AU within the orbits of giant planets at a few AU may be gone. In fact, out of ~30 % of survived test particles, about three quarters belong to the planet-planet collision cases. We find a good agreement between our numerical results and the secular theory, and present a semi-analytical formula which estimates the dynamically unstable region of the test particles just from the evolution of giant planets. Finally, our numerical results agree well with the observations, and also predict the existence of hot rocky planets in eccentric giant planet systems.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Secular Orbital Evolution of Compact Planet Systems

Recent observations have shown that at least some close-in exoplanets maintain eccentric orbits despite tidal circularization timescales that are typically shorter than stellar ages. We explore gravitational interactions with a distant planetary companion as a possible cause of these non-zero eccentricities. For simplicity, we focus on the evolution of a planar two-planet system subject to slow eccentricity damping and provide an intuitive interpretation of the resulting long-term orbital evolution. We show that dissipation shifts the two normal eigenmode frequencies and eccentricity ratios of the standard secular theory slightly, and that each mode decays at its own rate. Tidal damping of the eccentricities drives orbits to transition between periods of pericenter circulation and libration, and the planetary system settles into a locked state where the pericenters are nearly aligned or anti-aligned. Once in the locked state, the eccentricities of the two orbits decrease very slowly due to tides rather than at the much more rapid single-planet rate, and thus eccentric orbits, even for close-in planets, can often survive much longer than the age of the system. Assuming that an observed close-in planet on an elliptical orbit is apsidally-locked to a more distant, and perhaps unseen companion, we provide a constraint on the mass, semi-major axis, and eccentricity of the companion. We find the observed two-planet system HAT-P-13 might be in just such an apsidally-locked state, with parameters that obey our constraint well. We also survey close-in single planets, and found that none provide compelling evidence for unseen companions. Instead, we suspect that (1) orbits are circular, (2) tidal damping rates are slower than our assumption, or (3) a recent event has excited these eccentricities. Our method should prove useful for interpreting the results of current and future planet searches.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct Detection and Orbit Analysis of the Exoplanets HR 8799 bcd from Archival 2005 Keck/NIRC2 Data

We present previously unpublished July 2005 $H$-band coronagraphic data of the young, planet-hosting star HR 8799 from the newly-released Keck/NIRC2 archive. Despite poor observing conditions, we detect three of the planetary companions (HR 8799 bcd), two of them (HR 8799 bc) without advanced image processing. Comparing these data with previously published 1998-2011 astrometry and that from re-reduced October 2010 Keck data constrains the orbits of the planets. Analyzing the planets&#39; astrometry separately, HR 8799 d&#39;s orbit is likely inclined at least 25$^\circ$ from face-on and the others may be on in inclined orbits. For semimajor axis ratios consistent with a 4:2:1 mean-motion resonance, our analysis yields precise values for HR 8799 bcd&#39;s orbital parameters and strictly constrains the planets&#39; eccentricities to be less than 0.18--0.3. However, we find no acceptable orbital solutions with this resonance that place the planets in face-on orbits; HR 8799 d shows the largest deviation from such orbits. Moreover, few orbits make HR 8799 d coplanar with b and c, whereas dynamical stability analyses used to constrain the planets&#39; masses typically assume coplanar and/or face-on orbits. This paper illustrates the significant science gain enabled with the release of the NIRC2 archive.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Direct Imaging Confirmation and Characterization of a Dust-Enshrouded Candidate Exoplanet Orbiting Fomalhaut

We present Subaru/IRCS J band data for Fomalhaut and a (re)reduction of archival 2004--2006 HST/ACS data first presented by Kalas et al. (2008). We confirm the existence of a candidate exoplanet, Fomalhaut b, in both the 2004 and 2006 F606W data sets at a high signal-to-noise. Additionally, we confirm the detection at F814W and present a new detection in F435W. Fomalhaut b&#39;s space motion may be consistent with it being in an apsidally-aligned, non debris ring-crossing orbit, although new astrometry is required for firmer conclusions. We cannot confirm that Fomalhaut b exhibits 0.7-0.8 mag variability cited as evidence for planet accretion or a semi-transient dust cloud. The new, combined optical SED and IR upper limits confirm that emission identifying Fomalhaut b originates from starlight scattered by small dust, but this dust is most likely associated with a massive body. The Subaru and IRAC/4.5 micron upper limits imply M < 2 Mj, still consistent with the range of Fomalhaut b masses needed to sculpt the disk. Fomalhaut b is very plausibly &#34;a planet identified from direct imaging&#34; even if current images of it do not, strictly speaking, show thermal emission from a directly imaged planet.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A 5 Micron Image of beta Pictoris b at a Sub-Jupiter Projected Separation: Evidence for a Misalignment Between the Planet and the Inner, Warped Disk

We present and analyze a new M&#39; detection of the young exoplanet beta Pictoris b from 2008 VLT/NaCo data at a separation of ~ 4 AU and a high signal-to-noise rereduction of L&#39; data taken in December 2009. Based on our orbital analysis, the planet&#39;s orbit is viewed almost perfectly edge-on (i ~ 89 degrees) and has a Saturn-like semimajor axis of 9.50 (+3.93, -1.7) AU. Intriguingly, the planet&#39;s orbit is aligned with the major axis of the outer disk (Omega ~ 31 degrees) but probably misaligned with the warp/inclined disk at 80 AU often cited as a signpost for the planet&#39;s existence. Our results motivate new studies to clarify how $β$ Pic b sculpts debris disk structures and whether a second planet is required to explain the warp/inclined disk.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Combined Subaru/VLT/MMT 1--5 Micron Study of Planets Orbiting HR 8799: Implications for Atmospheric Properties, Masses, and Formation

We present new 1--1.25 micron (z and J band) Subaru/IRCS and 2 micron (K band) VLT/NaCo data for HR 8799 and a rereduction of the 3--5 micron MMT/Clio data first presented by Hinz et al. (2010). Our VLT/NaCo data yields a detection of a fourth planet at a projected separation of ~ 15 AU -- &#34;HR 8799e&#34;. We also report new, albeit weak detections of HR 8799b at 1.03 microns and 3.3 microns. Empirical comparisons to field brown dwarfs show that at least HR 8799b and HR8799c, and possibly HR 8799d, have near-to-mid IR colors/magnitudes significantly discrepant from the L/T dwarf sequence. Standard cloud deck atmosphere models appropriate for brown dwarfs provide only (marginally) statistically meaningful fits to HR 8799b and c for unphysically small radii. Models with thicker cloud layers not present in brown dwarfs reproduce the planets&#39; SEDs far more accurately and without the need for rescaling the planets&#39; radii. Our preliminary modeling suggests that HR 8799b has log(g) = 4--4.5, Teff = 900K, while HR 8799c, d, and (by inference) e have log(g) = 4--4.5, Teff = 1000--1200K. Combining results from planet evolution models and new dynamical stability limits implies that the masses of HR 8799b, c, d, and e are 6--7 Mj, 7--10 Mj, 7--10 Mj and 7--10 Mj. &#39;Patchy&#34; cloud prescriptions may provide even better fits to the data and may lower the estimated surface gravities and masses. Finally, contrary to some recent claims, forming the HR 8799 planets by core accretion is still plausible, although such systems are likely rare.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Tidal Evolution of Close-in Planets

Recent discoveries of several transiting planets with clearly non-zero eccentricities and some large inclinations started changing the simple picture of close-in planets having circular and well-aligned orbits. Two major scenarios to form such planets are planet migration in a disk, and planet--planet interactions combined with tidal dissipation. The former scenario can naturally produce a circular and low-obliquity orbit, while the latter implicitly assumes an initially highly eccentric and possibly high-obliquity orbit, which are then circularized and aligned via tidal dissipation. We investigate the tidal evolution of transiting planets on eccentric orbits. We show that the current and future orbital evolution of these systems is likely dominated by tidal dissipation, and not by a more distant companion. Although most of these close-in planets experience orbital decay all the way to the Roche limit, there are two characteristic evolution paths for them, depending on the relative efficiency of tidal dissipation inside the star and the planet. We point out that the current observations may be consistent with one of them. Our results suggest that at least some of the close-in planets with non-zero orbital eccentricity may have been formed by tidally circularizing an initially eccentric orbit. We also find that even when the stellar spin-orbit misalignment is observed to be small at present, some systems could have had a highly misaligned orbit in the past. Finally, we also re-examine the recent claim by Levrard et. al., who found that all orbital and spin parameters evolve on a similar timescale to orbital decay.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Unstable Planetary Systems Emerging Out Of Gas Disks

The discovery of over 400 extrasolar planets allows us to statistically test our understanding of formation and dynamics of planetary systems via numerical simulations. Traditional N-body simulations of multiple-planet systems without gas disks have successfully reproduced the eccentricity (e) distribution of the observed systems, by assuming that the planetary systems are relatively closely packed when the gas disk dissipates, so that they become dynamically unstable within the stellar lifetime. However, such studies cannot explain the small semi-major axes (a) of extrasolar planetary systems, if planets are formed, as the standard planet formation theory suggests, beyond the ice line. In this paper, we numerically study the evolution of three-planet systems in dissipating gas disks, and constrain the initial conditions that reproduce the observed semi-major axis and eccentricity distributions simultaneously. We adopt the initial conditions that are motivated by the standard planet formation theory, and self-consistently simulate the disk evolution, and planet migration by using a hybrid N-body and 1D gas disk code. We also take account of eccentricity damping, and investigate the effect of saturation of corotation resonances on the evolution of planetary systems. We find that the semi-major axis distribution is largely determined in a gas disk, while the eccentricity distribution is determined after the disk dissipation. We also find that there may be an optimum disk mass which leads to the observed a-e distribution. Our simulations generate a larger fraction of planetary systems trapped in mean-motion resonances (MMRs) than the observations, indicating that the disk&#39;s perturbation to the planetary orbits may be important to explain the observed rate of MMRs. We also find much lower occurrence of planets on retrograde orbits than the current observations of close-in planets suggest.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

The Growth & Migration of Jovian Planets in Evolving Protostellar Disks with Dead Zones

The growth of Jovian mass planets during migration in their protoplanetary disks is one of the most important problems that needs to be solved in light of observations of the exosolar planets. Studies of the migration of planets in standard gas disk models routinely show that migration is too fast to form Jovian planets, and that such migrating planetary cores generally plunge into the central stars in less than a Myr. In previous work, we have shown that a poorly ionized, less viscous region in a protoplanetary disk called a dead zone slows down the migration of fixed-mass planets. In this paper, we extend our numerical calculations to include dead zone evolution along with the disk, as well as planet formation via accretion of rocky and gaseous materials. Using our symplectic-integrator-gas dynamics code, we find that dead zones, even in evolving disks wherein migrating planets grow by accretion, still play a fundamental role in saving planetary systems. We demonstrate that Jovian planets form within 2.5 Myr for disks that are ten times more massive than a minimum mass solar nebula (MMSN) with an opacity reduction and without slowing down migration artificially. Our simulations indicate that protoplanetary disks with an initial mass comparable to the MMSN only produce Neptunian mass planets. We also find that planet migration does not help core accretion as much in the oligarchic planetesimal accretion scenario as it was expected in the runaway accretion scenario. Therefore we expect that an opacity reduction (or some other mechanisms) is needed to solve the formation timescale problem even for migrating protoplanets, as long as we consider the oligarchic growth. We also point out a possible role of a dead zone in explaining long-lived, strongly accreting gas disks.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

Saving Planetary Systems: Dead Zones & Planetary Migration

The tidal interaction between a disk and a planet leads to the planet&#39;s migration. A long-standing question regarding this mechanism is how to stop the migration before planets plunge into their central stars. In this paper, we propose a new, simple mechanism to significantly slow down planet migration, and test the possibility by using a hybrid numerical integrator to simulate the disk-planet interaction. The key component of the scenario is the role of low viscosity regions in protostellar disks known as dead zones, which affect planetary migration in two ways. First of all, it allows a smaller-mass planet to open a gap, and hence switch the faster type I migration to the slower type II migration. Secondly, a low viscosity slows down type II migration itself, because type II migration is directly proportional to the viscosity. We present numerical simulations of planetary migration by using a hybrid symplectic integrator-gas dynamics code. Assuming that the disk viscosity parameter inside the dead zone is (alpha=1e-4-1e-5), we find that, when a low-mass planet (e.g. 1-10 Earth masses) migrates from outside the dead zone, its migration is stopped due to the mass accumulation inside the dead zone. When a low-mass planet migrates from inside the dead zone, it opens a gap and slows down its migration. A massive planet like Jupiter, on the other hand, opens a gap and slows down inside the dead zone, independent of its initial orbital radius. The final orbital radius of a Jupiter mass planet depends on the dead zone&#39;s viscosity. For the range of alpha&#39;s noted above, this can vary anywhere from 7 AU, to an orbital radius of 0.1 AU that is characteristic of the hot Jupiters.

0 citations0 reviewsNo code linkedOpen access