Source author record

Kento Masuda

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

24works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

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

24 published item(s)

preprint2026arXiv

A young progenitor for the most common planetary systems in the Galaxy

The Galaxy's most common known planetary systems have several Earth-to-Neptune-size planets in compact orbits. At small orbital separations, larger planets are less common than their smaller counterparts by an order of magnitude. The young star V1298 Tau hosts one such compact planetary system, albeit with four planets that are uncommonly large (5 to 10 Earth radii). The planets form a chain of near-resonances that result in transit-timing variations of several hours. Here we present a multi-year campaign to characterize this system with transit-timing variations, a method insensitive to the intense magnetic activity of the star. Through targeted observations, we first resolved the previously unknown orbital period of the outermost planet. The full 9-year baseline from these and archival data then enabled robust determination of the masses and orbital parameters for all four planets. We find the planets have low, sub-Neptune masses and nearly circular orbits, implying a dynamically tranquil history. Their low masses and large radii indicate that the inner planets underwent a period of rapid cooling immediately after dispersal of the protoplanetary disk. Still, they are much less dense than mature planets of comparable size. We predict the planets will contract to 1.5-4.0 Earth radii and join the population of super-Earths and sub-Neptunes that nature produces in abundance.

preprint2026arXiv

PyIRD: A Python-Based Data Reduction Pipeline for Subaru/IRD and REACH

PyIRD is a Python-based pipeline for reducing spectroscopic data obtained with IRD (InfraRed Doppler; Kotani et al. (2018)) and REACH (Rigorous Exoplanetary Atmosphere Characterization with High dispersion coronagraphy; Kotani et al. (2020)) on the Subaru Telescope. It is designed to process raw images into one-dimensional spectra in a semi-automatic manner. Unlike traditional methods, it does not rely on IRAF (Tody, 1986; 1993), a software used for astronomical data reduction. This approach simplifies the workflow while maintaining efficiency and accuracy. Additionally, the pipeline includes an updated method for removing readout noise patterns from raw images, enabling efficient extraction of spectra even for faint targets such as brown dwarfs. The code is open source and available at https://github.com/prvjapan/pyird .

preprint2026arXiv

TOI-4495: A Pair of Aligned, Near-Resonant Sub-Neptunes that Likely Experienced Overstable Migration

We report the discovery of a sub-Neptune and a Neptune-like planet ($R_b = 2.48^{+0.14}_{-0.10}\,R_\oplus$, $R_c = 4.03^{+0.23}_{-0.15}\,R_\oplus$) orbiting the F-type star TOI-4495. The planets have orbital periods of 2.567 days and 5.185 days, lying close to a 2:1 mean-motion resonance (MMR). Our photodynamical analysis of the TESS light curves constrains the planetary masses to $M_b = 7.7 \pm 1.4\,M_\oplus$ and $M_c = 23.2 \pm 4.7\,M_\oplus$. The measured masses and radii indicate the presence of volatile-rich gaseous envelopes on both planets. The Rossiter-McLaughlin effect and the Doppler shadow of TOI-4495 c reveal a well-aligned orbit with a projected stellar obliquity of $λ= -2.3^{+8.3}_{-7.8}\,\mathrm{deg}$. Combined with the low mutual inclination constrained by the photodynamical analysis ($ΔI < 8.7\,\mathrm{deg}$), the planetary orbits are likely coplanar and aligned with the host star's spin axis. We show that the planets are near, but not in, the 2:1 MMR, with a circulating resonant angle. We also find substantial free eccentricity for the inner planet, TOI-4495 b ($e_b = 0.078^{+0.020}_{-0.013}$). Given the observed proximity to the 2:1 resonance and the more massive outer planet, TOI-4495 b and c are particularly susceptible to resonant overstability, which can convert resonantly excited eccentricity into free eccentricity. However, additional mechanisms (e.g., planetesimal scattering) may be required to further excite the eccentricity by $\sim 4\%$. To prevent tidal damping from reducing the eccentricity below the observed level over the star's lifetime (1.9 Gyr), the reduced tidal quality factor of TOI-4495 b must be $Q' \gtrsim 10^5$, consistent with the presence of a thick envelope on the planet.

preprint2022arXiv

Dippers from TESS Full-frame Images. II. Spectroscopic Characterization of Four Young Dippers

Photometric monitoring by the Transiting Exoplanet Survey Satellite (TESS) has discovered not only periodic signals by transiting exoplanets but also episodic or quasi-periodic dimming around young stellar objects. The mechanisms of the dimming of these objects, so-called "dippers", are thought to be related to the property of the accretion or the structure of protoplanetary disks especially in regions close to the host star. Recently, we have created the catalog of dippers from the one year of TESS Full Frame Image (FFI) data. In this paper, we report spectral features of four newly found dippers in that catalog and show that they potentially shed light on the dimming mechanisms. We found that all of the targets exhibit the H alpha emission line, which is an indicator of an ongoing accretion. Based on its line profiles and/or their variability, we characterized the properties of the disks of each source, which can support the dimming mechanisms due to a dusty disk wind or an accretion warp. Also, we found an interesting dipper (TIC 317873721), "close-in binary dipper," showing the complex variability of the line profile and the large radial velocity variation. Since the dimming intervals are similar to the orbital period of the binary, we suggest that the dips are caused by dust in the warp of accretion from a circumbinary disk onto stars. Such a close-in ($<$ 0.1 au) binary dipper is rarely reported so far, further investigation will reveal the new aspect of the disk evolution and planetary formation.

preprint2022arXiv

On the Evolution of Rotational Modulation Amplitude in Solar-mass Main-sequence Stars

We investigate the relation between rotation periods $P_\mathrm{rot}$ and photometric modulation amplitudes $R_\mathrm{per}$ for $\approx 4,000$ Sun-like main-sequence stars observed by Kepler, using $P_\mathrm{rot}$ and $R_\mathrm{per}$ from McQuillan et al. (2014), effective temperature $T_\mathrm{eff}$ from LAMOST DR6, and parallax data from Gaia EDR3. As has been suggested in previous works, we find that $P_\mathrm{rot}$ scaled by the convective turnover time $τ_\mathrm{c}$, or the Rossby number $\mathrm{Ro}=P_\mathrm{rot}/τ_\mathrm{c}$, serves as a good predictor of $R_\mathrm{per}$: $R_\mathrm{per}$ plateaus around $1\%$ in relative flux for $0.2 \lesssim \mathrm{Ro}/\mathrm{Ro}_\odot \lesssim 0.4$, and decays steeply with increasing $\mathrm{Ro}$ for $0.4 \lesssim \mathrm{Ro}/\mathrm{Ro}_\odot \lesssim 0.8$, where $\mathrm{Ro}_\odot$ denotes $\mathrm{Ro}$ of the Sun. In the latter regime we find $\mathrm{d}\ln R_\mathrm{per}/\mathrm{d}\ln\mathrm{Ro} \sim -4.5$ to $-2.5$, although the value is sensitive to detection bias against weak modulation and may depend on other parameters including $T_\mathrm{eff}$ and surface metallicity. The existing X-ray and Ca II H&K flux data also show transitions at $\mathrm{Ro}/\mathrm{Ro}_\odot\sim 0.4$, suggesting that all these transitions share the same physical origin. We also find that the rapid decrease of $R_\mathrm{per}$ with increasing $\mathrm{Ro}$ causes rotational modulation of fainter Kepler stars with $\mathrm{Ro}/\mathrm{Ro}_\odot \gtrsim 0.6$ to be buried under the photometric noise. This effect sets the longest $P_\mathrm{rot}$ detected in the McQuillan et al. (2014) sample as a function of $T_\mathrm{eff}$, and obscures the signature of stalled spin down that has been proposed to set in around $\mathrm{Ro}/\mathrm{Ro}_\odot \sim 1$.

preprint2022arXiv

TOI-1136 is a Young, Coplanar, Aligned Planetary System in a Pristine Resonant Chain

Convergent disk migration has long been suspected to be responsible for forming planetary systems with a chain of mean-motion resonances (MMR). Dynamical evolution over time could disrupt the delicate resonant configuration. We present TOI-1136, a 700-Myr-old G star hosting at least 6 transiting planets between $\sim$2 and 5 $R_\oplus$. The orbital period ratios deviate from exact commensurability by only $10^{-4}$, smaller than the $\sim$\,$10^{-2}$ deviations seen in typical Kepler near-resonant systems. A transit-timing analysis measured the masses of the planets (3-8$M_\oplus$) and demonstrated that the planets in TOI-1136 are in true resonances with librating resonant angles. Based on a Rossiter-McLaughlin measurement of planet d, the star's rotation appears to be aligned with the planetary orbital planes. The well-aligned planetary system and the lack of detected binary companion together suggest that TOI-1136's resonant chain formed in an isolated, quiescent disk with no stellar fly-by, disk warp, or significant axial asymmetry. With period ratios near 3:2, 2:1, 3:2, 7:5, and 3:2, TOI-1136 is the first known resonant chain involving a second-order MMR (7:5) between two first-order MMR. The formation of the delicate 7:5 resonance places strong constraints on the system's migration history. Short-scale (starting from $\sim$0.1 AU) Type-I migration with an inner disk edge is most consistent with the formation of TOI-1136. A low disk surface density ($Σ_{\rm 1AU}\lesssim10^3$g~cm$^{-2}$; lower than the minimum-mass solar nebula) and the resultant slower migration rate likely facilitated the formation of the 7:5 second-order MMR. TOI-1136's deep resonance suggests that it has not undergone much resonant repulsion during its 700-Myr lifetime. One can rule out rapid tidal dissipation within a rocky planet b or obliquity tides within the largest planets d and f.

preprint2021arXiv

Auto-Differentiable Spectrum Model for High-Dispersion Characterization of Exoplanets and Brown Dwarfs

We present an auto-differentiable spectral modeling of exoplanets and brown dwarfs. This model enables a fully Bayesian inference of the high--dispersion data to fit the ab initio line-by-line spectral computation to the observed spectrum by combining it with the Hamiltonian Monte Carlo in recent probabilistic programming languages. An open source code, exojax, developed in this study, was written in Python using the GPU/TPU compatible package for automatic differentiation and accelerated linear algebra, JAX (Bradbury et al. 2018). We validated the model by comparing it with existing opacity calculators and a radiative transfer code and found reasonable agreements of the output. As a demonstration, we analyzed the high-dispersion spectrum of a nearby brown dwarf, Luhman 16 A and found that a model including water, carbon monoxide, and $\mathrm{H_2/He}$ collision induced absorption was well fitted to the observed spectrum ($R=10^5$ and 2.28-2.30 $μ$m). As a result, we found that $T_0=1295_{-32}^{+35}$ K at 1 bar and C/O $=0.62 \pm 0.03$, which is slightly higher than the solar value. This work demonstrates the potential of full Bayesian analysis of brown dwarfs and exoplanets as observed by high-dispersion spectrographs and also directly-imaged exoplanets as observed by high-dispersion coronagraphy.

preprint2020arXiv

Bayesian Dynamic Mapping of an Exo-Earth from Photometric Variability

Photometric variability of a directly imaged exo-Earth conveys spatial information on its surface and can be used to retrieve a two-dimensional geography and axial tilt of the planet (spin-orbit tomography). In this study, we relax the assumption of the static geography and present a computationally tractable framework for dynamic spin-orbit tomography applicable to the time-varying geography. First, a Bayesian framework of static spin-orbit tomography is revisited using analytic expressions of the Bayesian inverse problem with a Gaussian prior. We then extend this analytic framework to a time-varying one through a Gaussian process in time domain, and present analytic expressions that enable efficient sampling from a full joint posterior distribution of geography, axial tilt, spin rotation period, and hyperparameters in the Gaussian-process priors. Consequently, it only takes 0.3 s for a laptop computer to sample one posterior dynamic map conditioned on the other parameters with 3,072 pixels and 1,024 time grids, for a total of $\sim 3 \times 10^6$ parameters. We applied our dynamic mapping method on a toy model and found that the time-varying geography was accurately retrieved along with the axial-tilt and spin rotation period. In addition, we demonstrated the use of dynamic spin-orbit tomography with a real multi-color light curve of the Earth as observed by the Deep Space Climate Observatory. We found that the resultant snapshots from the dominant component of a principle component analysis roughly captured the large-scale, seasonal variations of the clear-sky and cloudy areas on the Earth.

preprint2020arXiv

Hot Stars With Kepler Planets Have High Obliquities

It has been known for a decade that hot stars with hot Jupiters tend to have high obliquities. Less is known about the degree of spin-orbit alignment for hot stars with other kinds of planets. Here, we re-assess the obliquities of hot Kepler stars with transiting planets smaller than Neptune, based on spectroscopic measurements of their projected rotation velocities (vsini). The basis of the method is that a lower obliquity -- all other things being equal -- causes sini to be closer to unity and increases the value of vsini. We sought evidence for this effect using a sample of 150 Kepler stars with effective temperatures between 5950 and 6550K and a control sample of 101 stars with matching spectroscopic properties and random orientations. The planet hosts have systematically higher values of vsini than the control stars, but not by enough to be compatible with perfect spin-orbit alignment. The mean value of sini is 0.856 +/- 0.036, which is 4-sigma away from unity (perfect alignment), and 2-sigma away from pi/4 (random orientations). There is also evidence that the hottest stars have a broader obliquity distribution: when modeled separately, the stars cooler than 6250K have <sini> = 0.928 +/- 0.042, while the hotter stars are consistent with random orientations. This is similar to the pattern previously noted for stars with hot Jupiters. Based on these results, obliquity excitation for early-G and late-F stars appears to be a general outcome of star and planet formation, rather than being exclusively linked to hot Jupiter formation.

preprint2020arXiv

Mutual Orbital Inclinations Between Cold Jupiters and Inner Super-Earths

Previous analyses of Doppler and Kepler data have found that Sun-like stars hosting "cold Jupiters" (giant planets with $a\gtrsim 1\,\mathrm{AU}$) almost always host "inner super-Earths" (1-$4\,R_\oplus$, $a\lesssim1\,\mathrm{AU}$). Here, we attempt to determine the degree of alignment between the orbital planes of the cold Jupiters and the inner super-Earths. The key observational input is the fraction of Kepler stars with transiting super-Earths that also have transiting cold Jupiters. This fraction depends on both the probability for cold Jupiters to occur in such systems, and on the mutual orbital inclinations. Since the probability of occurrence has already been measured in Doppler surveys, we can use the data to constrain the dispersion of the mutual inclination distribution. We find $σ=11.8^{+12.7}_{-5.5}\,\mathrm{deg}$ (68% confidence) and $σ>3.5\,\mathrm{deg}$ (95% confidence), where $σ$ is the scale parameter of the Rayleigh distribution. This suggests that planetary orbits in systems with cold Jupiters tend to be coplanar - although not quite as coplanar as those in the Solar System, which have a mean inclination from the invariable plane of $1.8\,\mathrm{deg}$. We also find evidence that cold Jupiters have lower mutual inclinations relative to inner systems with higher transit multiplicity. This suggests a link between the dynamical excitation in the inner and outer systems. For example, perturbations from misaligned cold Jupiters may dynamically heat or destabilize systems of inner super-Earths.

preprint2020arXiv

On the Inference of a Star's Inclination Angle from its Rotation Velocity and Projected Rotation Velocity

It is possible to learn about the orientation of a star's rotation axis by combining measurements of the star's rotation velocity ($v$) and its projection onto our line of sight ($v\sin i$). This idea has found many applications, including the investigation of the obliquities of stars with transiting planets. Here, we present a method for the probabilistic inference of the inclination of the star's rotation axis based on independent data sets that constrain $v$ and $v\sin i$. We also correct several errors and misconceptions that appear in the literature.

preprint2019arXiv

It takes two planets in resonance to tango around K2-146

K2-146 is a cool, 0.358 M_sun dwarf that was found to host a mini-Neptune with a 2.67-days period. The planet exhibited strong transit timing variations (TTVs) of greater than 30 minutes, indicative of the presence of a further object in the system. Here we report the discovery of the previously undetected outer planet, K2-146 c, in the system using additional photometric data. K2-146 c was found to have a grazing transit geometry and a 3.97-day period. The outer planet was only significantly detected in the latter K2 campaigns presumably because of precession of its orbital plane. The TTVs of K2-146 b and c were measured using observations spanning a baseline of almost 1200 days. We found strong anti-correlation in the TTVs, suggesting the two planets are gravitationally interacting. Our TTV and transit model analyses revealed that K2-146 b has a radius of 2.25 $\pm$ 0.10 \R_earth and a mass of 5.6 $\pm$ 0.7 M_earth, whereas K2-146 c has a radius of $2.59_{-0.39}^{+1.81}$ R_earth and a mass of 7.1 $\pm$ 0.9 M_earth. The inner and outer planets likely have moderate eccentricities of $e = 0.14 \pm 0.07$ and $0.16 \pm 0.07$, respectively. Long-term numerical integrations of the two-planet orbital solution show that it can be dynamically stable for at least 2 Myr. The evaluation of the resonance angles of the planet pair indicates that K2-146 b and c are likely trapped in a 3:2 mean motion resonance. The orbital architecture of the system points to a possible convergent migration origin.

preprint2019arXiv

The Featureless Transmission Spectra of Two Super-Puff Planets

The Kepler mission revealed a class of planets known as ''super-puffs,'' with masses only a few times larger than Earth's but radii larger than Neptune, giving them very low mean densities. All three of the known planets orbiting the young solar-type star Kepler 51 are super-puffs. The Kepler 51 system thereby provides an opportunity for a comparative study of the structures and atmospheres of this mysterious class of planets, which may provide clues about their formation and evolution. We observed two transits each of Kepler 51b and 51d with the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. Combining new WFC3 transit times with re-analyzed Kepler data and updated stellar parameters, we confirmed that all three planets have densities lower than 0.1 g/cm$^{3}$. We measured the WFC3 transmission spectra to be featureless between 1.15 and 1.63 $μ$m, ruling out any variations greater than 0.6 scale heights (assuming a H/He dominated atmosphere), thus showing no significant water absorption features. We interpreted the flat spectra as the result of a high-altitude aerosol layer (pressure $<$3 mbar) on each planet. Adding this new result to the collection of flat spectra that have been observed for other sub-Neptune planets, we find support for one of the two hypotheses introduced by Crossfield and Kreidberg (2017), that planets with cooler equilibrium temperatures have more high-altitude aerosols. We strongly disfavor their other hypothesis that the H/He mass fraction drives the appearance of large amplitude transmission features.

preprint2016arXiv

Transiting Planet Candidates Beyond the Snow Line Detected by Visual Inspection of 7557 Kepler Objects of Interest

We visually inspected the light curves of 7557 Kepler Objects of Interest (KOIs) to search for single transit events (STEs) possibly due to long-period giant planets. We identified 28 STEs in 24 KOIs, among which 14 events are newly reported in this paper. We estimate the radius and orbital period of the objects causing STEs by fitting the STE light curves simultaneously with the transits of the other planets in the system or with the prior information on the host star density. As a result, we found that STEs in seven of those systems are consistent with Neptune- to Jupiter-sized objects of orbital periods ranging from a few to $\sim$ $20\,\mathrm{yr}$. We also estimate that $\gtrsim20\%$ of the compact multi-transiting systems host cool giant planets with periods $\gtrsim 3\,\mathrm{yr}$ on the basis of their occurrence in the KOIs with multiple candidates, assuming the small mutual inclination between inner and outer planetary orbits.

preprint2016arXiv

Transiting planets as a precision clock to constrain the time variation of the gravitational constant

Analysis of transit times in exoplanetary systems accurately provides an instantaneous orbital period, $P(t)$, of their member planets. A long-term monitoring of those transiting planetary systems puts limits on the variability of $P(t)$, which are translated into the constraints on the time variation of the gravitational constant $G$. We apply this analysis to $10$ transiting systems observed by the Kepler spacecraft, and find that $ΔG/G\lesssim 5\times10^{-6}$ for 2009-2013, or $\dot{G}/G \lesssim 10^{-6}\,\mathrm{yr}^{-1}$ if $\dot{G}$ is constant. While the derived limit is weaker than those from other analyses, it is complementary to them and can be improved by analyzing numerous transiting systems that are continuously monitored.

preprint2015arXiv

Absolute Dimensions of a Flat Hierarchical Triple System KIC 6543674 from the Kepler Photometry

Many of the Kepler close binaries are suggested to constitute hierarchical triple systems through their eclipse timing variations (ETVs). Eclipses by the third body in those systems, if observed, provide precise constraints on its physical and orbital properties, which are otherwise difficult to obtain. In this Letter, we analyze such a "tertiary event" observed only once in the KIC 6543674 system. The system consists of a short-period ($2.4\,\mathrm{days}$) inner eclipsing binary and a third body on a wide ($1100\,\mathrm{days}$) and eccentric ($e\simeq0.6$) orbit. Analysis of three tertiary eclipses around a single inferior conjunction of the third body yields the mutual inclination between the inner and outer binary planes to be $3.3^\circ\pm0.6^\circ$, indicating an extremely flat geometry. Furthermore, combining the timings and shapes of the tertiary eclipses with the phase curve and ETVs of the inner binary, we determine the mass and radius ratios of all three bodies in the system using the Kepler photometry alone. With the primary mass and temperature from the Kepler Input Catalog, the absolute masses, radii, and effective temperatures of the three stars are obtained as follows: $M_\mathrm{A}=1.2\pm0.3\,M_\odot$, $R_\mathrm{A}=1.8\pm0.1\,R_\odot$, $M_\mathrm{B}=1.1_{-0.2}^{+0.3}\,M_\odot$, $R_\mathrm{B}=1.4\pm0.1\,R_\odot$, $M_\mathrm{C}=0.50_{-0.08}^{+0.07}\,M_\odot$, $R_\mathrm{C}=0.50\pm0.04\,R_\odot$, $T_\mathrm{A} \simeq T_\mathrm{B}\simeq 6100\,\mathrm{K}$, and $T_\mathrm{C}<5000\,\mathrm{K}$. Implication for the formation scenario of close binaries is briefly discussed.

preprint2015arXiv

Precise Radial Velocity Measurements for Kepler Giants Hosting Planetary Candidates: Kepler-91 and KOI-1894

We present results of radial-velocity follow-up observations for the two Kepler evolved stars Kepler-91 (KOI-2133) and KOI-1894, which had been announced as candidates to host transiting giant planets, with the Subaru 8.2m telescope and the High Dispersion Spectrograph (HDS). By global modeling of the high-precision radial-velocity data taken with Subaru/HDS and photometric ones taken by Kepler mission taking account of orbital brightness modulations (ellipsoidal variations, reflected/emitted light, etc.) of the host stars, we independently confirmed that Kepler-91 hosts a transiting planet with a mass of 0.66 M_Jup (Kepler-91b), and newly detected an offset of ~20 m s$^{-1}$ between the radial velocities taken at ~1-yr interval, suggesting the existence of additional companion in the system. As for KOI-1894, we detected possible phased variations in the radial velocities and light curves with 2--3 sigma confidence level which could be explained as a reflex motion and ellipsoidal variation of the star caused by the transiting sub-saturn-mass (~0.18 M_Jup) planet.

preprint2015arXiv

Revisiting a gravity-darkened and precessing planetary system PTFO 8-8695: spin-orbit non-synchronous case

We reanalyse the time-variable lightcurves of the transiting planetary system PTFO 8-8695, in which a planet of 3 to 4 Jupiter mass orbits around a rapidly rotating pre-main-sequence star. Both the planetary orbital period of 0.448 days and the stellar spin period less than 0.671 days are unusually short, which makes PTFO 8-8695 an ideal system to check the model of gravity darkening and nodal precession. While the previous analysis of PTFO 8-8695 assumed that the stellar spin and planetary orbital periods are the same, we extend the analysis by discarding the spin-orbit synchronous condition, and find three different classes of solutions roughly corresponding to the nodal precession periods of 199$\pm$16, 475$\pm$21, and 827$\pm$53 days that reproduce the transit lightcurves observed in 2009 and 2010. We compare the predicted lightcurves of the three solutions against the photometry data of a few percent accuracy obtained at Koyama Astronomical Observatory in 2014 and 2015, and find that the solution with the precession period of 199$\pm$16 days is preferred even though preliminary. Future prospect and implications to other transiting systems are briefly discussed.

preprint2015arXiv

Spin-Orbit Angles of Kepler-13Ab and HAT-P-7b from Gravity-Darkened Transit Light Curves

Analysis of the transit light curve deformed by the stellar gravity darkening allows us to photometrically measure both components of the spin-orbit angle $ψ$, its sky projection $λ$ and inclination of the stellar spin axis $i_\star$. In this paper, we apply the method to two transiting hot Jupiter systems monitored with the Kepler spacecraft, Kepler-13A and HAT-P-7. For Kepler-13A, we find $i_\star=81^\circ\pm5^\circ$ and $ψ=60^\circ\pm2^\circ$ adopting the spectroscopic constraint $λ=58.6^\circ\pm2.0^\circ$ by Johnson et al. (2014). In our solution, the discrepancy between the above $λ$ and that previously reported by Barnes et al. (2011) is solved by fitting both of the two parameters in the quadratic limb-darkening law. We also report the temporal variation in the orbital inclination of Kepler-13Ab, $\mathrm{d} |\cos i_{\rm orb}|/\mathrm{d}t=(-7.0\pm0.4)\times10^{-6}\,\mathrm{day}^{-1}$, providing further evidence for the spin-orbit precession in this system. By fitting the precession model to the time series of $i_{\rm orb}$, $λ$, and $i_\star$ obtained with the gravity-darkened model, we constrain the stellar quadrupole moment $J_2=(6.1\pm0.3)\times10^{-5}$ for our new solution, which is several times smaller than $J_2=(1.66\pm0.08)\times10^{-4}$ obtained for the previous one. We show that the difference can be observable in the future evolution of $λ$, thus providing a possibility to test our solution with follow-up observations. The second target, HAT-P-7, is the first F-dwarf star analyzed with the gravity-darkening method. Our analysis points to a nearly pole-on configuration with $ψ=101^\circ\pm2^\circ$ or $87^\circ\pm2^\circ$ and the gravity-darkening exponent $β$ consistent with $0.25$. Such an observational constraint on $β$ can be useful for testing the theory of gravity darkening.

preprint2014arXiv

Global Analysis of KOI-977: Spectroscopy, Asteroseismology, and Phase-curve Analysis

We present a global analysis of KOI-977, one of the planet host candidates detected by {\it Kepler}. Kepler Input Catalog (KIC) reports that KOI-977 is a red giant, for which few close-in planets have been discovered. Our global analysis involves spectroscopic and asteroseismic determinations of stellar parameters (e.g., mass and radius) and radial velocity (RV) measurements. Our analyses reveal that KOI-977 is indeed a red giant in the red clump, but its estimated radius ($\gtrsim 20R_\odot=0.093$ AU) is much larger than KOI-977.01's orbital distance ($\sim 0.027$ AU) estimated from its period ($P_\mathrm{orb}\sim 1.35$ days) and host star's mass. RV measurements show a small variation, which also contradicts the amplitude of ellipsoidal variations seen in the light-curve folded with KOI-977.01's period. Therefore, we conclude that KOI-977.01 is a false positive, meaning that the red giant, for which we measured the radius and RVs, is different from the object that produces the transit-like signal (i.e., an eclipsing binary). On the basis of this assumption, we also perform a light-curve analysis including the modeling of transits/eclipses and phase-curve variations, adopting various values for the dilution factor $D$, which is defined as the flux ratio between the red giant and eclipsing binary. Fitting the whole folded light-curve as well as individual transits in the short cadence data simultaneously, we find that the estimated mass and radius ratios of the eclipsing binary are consistent with those of a solar-type star and a late-type star (e.g., an M dwarf) for $D\gtrsim 20$.

preprint2014arXiv

Tidal evolution of the spin-orbit angle in exoplanetary systems

The angle between the stellar spin and the planetary orbit axes (spin-orbit angle) is supposed to carry valuable information on the initial condition of the planet formation and the subsequent migration history. Indeed current observations of the Rossiter- McLaughlin effect have revealed a wide range of spin-orbit misalignments for transiting exoplanets. We examine in detail the tidal evolution of a simple system comprising a Sun-like star and a hot Jupiter adopting the equilibrium tide and the inertial wave dissipation effects simultaneously. We find that the combined tidal model works as a very efficient realignment mechanism; it predicts three distinct states of the spin-orbit angle (i.e., parallel, polar, and anti-parallel orbits) for a while, but the latter two states eventually approach the parallel spin-orbit configuration. The intermediate spin-orbit angles as measured in recent observations are difficult to be achieved. Therefore the current model cannot reproduce the observed broad distribution of the spin-orbit angles, at least in its simple form. This indicates that the observed diversity of the spin-orbit angles may emerge from more complicated interactions with outer planets and/or may be the consequence of the primordial misalignment between the proto-planetary disk and the stellar spin, which requires future detailed studies.

preprint2014arXiv

Very Low-Density Planets around Kepler-51 Revealed with Transit Timing Variations and an Anomaly Similar to a Planet-Planet Eclipse Event

We present an analysis of the transit timing variations (TTVs) in the multi-transiting planetary system around Kepler-51 (KOI-620). This system consists of two confirmed transiting planets, Kepler-51b ($P_{\rm b} = 45.2\,\mathrm{days}$) and Kepler-51c ($P_{\rm c} = 85.3\,\mathrm{days}$), and one transiting planet candidate KOI-620.02 ($P_{02} = 130.2\,\mathrm{days}$), which lie close to a $1:2:3$ resonance chain. Our analysis shows that their TTVs are consistently explained by the three-planet model, and constrains their masses as $M_{\rm b} = 2.1_{-0.8}^{+1.5} M_{\oplus}$ (Kepler-51b), $M_{\rm c} = 4.0 \pm 0.4 M_{\oplus}$ (Kepler-51c), and $M_{02} = 7.6 \pm 1.1 M_{\oplus}$ (KOI-620.02), thus confirming KOI-620.02 as a planet in this system. The masses inferred from the TTVs are rather small compared to the planetary radii based on the stellar density and planet-to-star radius ratios determined from the transit light curves. Combining these estimates, we find that all three planets in this system have densities among the lowest determined, $ρ_p \lesssim 0.05\,{\rm g\,cm^{-3}}$. With this feature, the Kepler-51 system serves as another example of low-density compact multi-transiting planetary systems. We also identify a curious feature in the archived Kepler light curve during the double transit of Kepler-51b and KOI-620.02, which could be explained by their overlapping on the stellar disk (a planet-planet eclipse). If this is really the case, the sky-plane inclination of KOI-620.02's orbit relative to that of Kepler-51b is given by $ΔΩ= -25.3_{-6.8}^{+6.2}\mathrm{deg}$, implying significant misalignment of their orbital planes. This interpretation, however, seems unlikely because such an event that is consistent with all of the observations is found to be exceedingly rare.

preprint2013arXiv

Characterization of the KOI-94 System with Transit Timing Variation Analysis: Implication for the Planet-Planet Eclipse

The KOI-94 system is a closely-packed, multi-transiting planetary system discovered by the Kepler space telescope. It is known as the first system that exhibited a rare event called a "planet-planet eclipse (PPE)," in which two planets partially overlap with each other in their double-transit phase. In this paper, we constrain the parameters of the KOI-94 system with an analysis of the transit timing variations (TTVs). Such constraints are independent of the radial velocity (RV) analysis recently performed by Weiss and coworkers, and valuable in examining the reliability of the parameter estimate using TTVs. We numerically fit the observed TTVs of KOI-94c, KOI-94d, and KOI-94e for their masses, eccentricities, and longitudes of periastrons, and obtain the best-fit parameters including $m_{\rm c} = 9.4_{-2.1}^{+2.4} M_{\oplus}$, $m_{\rm d} = 52.1_{-7.1}^{+6.9} M_{\oplus}$, $m_{\rm e} = 13.0_{-2.1}^{+2.5} M_{\oplus}$, and $e \lesssim 0.1$ for all the three planets. While these values are mostly in agreement with the RV result, the mass of KOI-94d estimated from the TTV is significantly smaller than the RV value $m_{\rm d} = 106 \pm 11 M_{\oplus}$. In addition, we find that the TTV of the outermost planet KOI-94e is not well reproduced in the current modeling. We also present analytic modeling of the PPE and derive a simple formula to reconstruct the mutual inclination of the two planets from the observed height, central time, and duration of the brightening caused by the PPE. Based on this model, the implication of the results of TTV analysis for the time of the next PPE is discussed.

preprint2012arXiv

Planet-Planet Eclipse and the Rossiter-McLaughlin Effect of a Multiple Transiting System: Joint Analysis of the Subaru Spectroscopy and the Kepler Photometry

We report a joint analysis of the Rossiter-McLaughlin (RM) effect with Subaru and the Kepler photometry for Kepler Object of Interest (KOI) 94 system. The system comprises four transiting planet candidates with orbital periods of 22.3 (KOI-94.01), 10.4 (KOI-94.02), 54.3 (KOI-94.03), and 3.7 (KOI-94.04) days from the Kepler photometry. We performed the radial velocity (RV) measurement of the system with the Subaru 8.2 m telescope on August 10, 2012 (UT), covering a complete transit of KOI-94.01 for $\sim 6.7$ hours. The resulting RV variation due to the RM effect spectroscopically confirms that KOI-94.01 is indeed the transiting planet, and implies that its orbital axis is well aligned with the stellar spin axis; the projected spin-orbit angle $λ$ is estimated as $-6_{-11}^{+13}$ deg. This is the first measurement of the RM effect for a multiple transiting system. Remarkably, the archived Kepler lightcurve around BJD=2455211.5 (date in UT January 14/15, 2010) indicates a "ouble transit" event of KOI-94.01 and KOI-94.03, in which the two planets transit the stellar disk simultaneously. Moreover, the two planets partially overlap each other, and exhibit a "planet-planet eclipse" around the transit center. This provides a rare opportunity to put tight constraints on the configuration of the two transiting planets by joint analysis with our Subaru RM measurement. Indeed, we find that the projected mutual inclination of KOI-94.01 and KOI-94.03 is estimated to be $δ= -1.15 \pm 0.55$ deg. Implications for the migration model of multiple planet systems are also discussed.