Showing 1–50 of 189 results for author: Suto, Y

Cosmic evolution of black hole-spin and galaxy orientations: clues from the NewHorizon and Galactica simulations

Authors: Sebastien Peirani, Yasushi Suto, Ricarda S. Beckmann, Marta Volonteri, Yen-Ting Lin, Yohan Dubois, Sukyoung K. Yi, Christophe Pichon, Katarina Kraljic, Minjung Park, Julien Devriendt, San Han, Wei-Huai Chen

Abstract: (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the a… ▽ More (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses. △ Less

Submitted 25 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A (12/03/2024)

Journal ref: A&A 686, A233 (2024)

Formation of unipolar outflow and $\textit{protostellar rocket effect}$ in magnetized turbulent molecular cloud cores

Authors: Daisuke Takaishi, Yusuke Tsukamoto, Miyu Kido, Shigehisa Takakuwa, Yoshiaki Misugi, Yuki Kudoh, Yasushi Suto

Abstract: Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular m… ▽ More Observed protostellar outflows exhibit a variety of asymmetrical features, including remarkable unipolar outflows and bending outflows. Revealing the formation and early evolution of such asymmetrical protostellar outflows, especially the unipolar outflows, is essential for a better understanding of the star and planet formation because they can dramatically change the mass accretion and angular momentum transport to the protostars and protoplanetary disks. Here, we perform the three-dimensional non-ideal magnetohydrodynamics simulations to investigate the formation and early evolution of the asymmetrical protostellar outflows in magnetized turbulent isolated molecular cloud cores. We find, for the first time to our knowledge, that the unipolar outflow forms even in the single low-mass protostellar system. The results show that the unipolar outflow is driven in the weakly magnetized cloud cores with the dimensionless mass-to-flux ratios of $μ=8$ and $16$. Furthermore, we find the $\textit{protostellar rocket effect}$ of the unipolar outflow, which is similar to the launch and propulsion of a rocket. The unipolar outflow ejects the protostellar system from the central dense region to the outer region of the parent cloud core, and the ram pressure caused by its ejection suppresses the driving of additional new outflows. In contrast, the bending bipolar outflow is driven in the moderately magnetized cloud core with $μ=4$. The ratio of the magnetic to turbulent energies of a parent cloud core may play a key role in the formation of asymmetrical protostellar outflows. △ Less

Submitted 4 January, 2024; originally announced January 2024.

Comments: 24 pages, 6 figures, accepted for publication in ApJ

Constraining the binarity of black hole candidates: a proof-of-concept study of Gaia BH1 and Gaia BH2

Authors: Toshinori Hayashi, Yasushi Suto, Alessandro A. Trani

Abstract: Nearly a hundred of binary black holes (BBHs) have been discovered with gravitational-wave signals emitted at their merging events. Thus, it is quite natural to expect that significantly more abundant BBHs with wider separations remain undetected in the universe, or even in our Galaxy. We consider a possibility that star-BH binary candidates may indeed host an inner BBH, instead of a single BH. We… ▽ More Nearly a hundred of binary black holes (BBHs) have been discovered with gravitational-wave signals emitted at their merging events. Thus, it is quite natural to expect that significantly more abundant BBHs with wider separations remain undetected in the universe, or even in our Galaxy. We consider a possibility that star-BH binary candidates may indeed host an inner BBH, instead of a single BH. We present a detailed feasibility study of constraining the binarity of the currently available two targets, Gaia BH1 and Gaia BH2. Specifically, we examine three types of radial velocity (RV) modulations of a tertiary star in star-BBH triple systems; short-term RV modulations induced by the inner BBH, long-term RV modulations induced by the nodal precession, and long-term RV modulations induced by the von Zeipel-Kozai-Lidov oscillations. Direct three-body simulations combined with approximate analytic models reveal that Gaia BH1 system may exhibit observable signatures of the hidden inner BBH if it exists at all. The methodology that we examine here is quite generic, and is expected to be readily applicable to future star-BH binary candidates in a straightforward manner. △ Less

Submitted 29 August, 2023; v1 submitted 4 July, 2023; originally announced July 2023.

Comments: 21 pages, 11 figures, 1 table, ApJ, in press

Modeling photometric variations due to a global inhomogeneity on an obliquely rotating star: application to lightcurves of white dwarfs

Authors: Yasushi Suto, Shin Sasaki, Masataka Aizawa, Kotaro Fujisawa, Kazumi Kashiyama

Abstract: We develop a general framework to compute photometric variations induced by the oblique rotation of a star with an axisymmetric inhomogeneous surface. We apply the framework to compute lightcurves of white dwarfs adopting two simple models of their surface inhomogeneity. Depending on the surface model and the location of the observer, the resulting lightcurve exhibits a departure from a purely sin… ▽ More We develop a general framework to compute photometric variations induced by the oblique rotation of a star with an axisymmetric inhomogeneous surface. We apply the framework to compute lightcurves of white dwarfs adopting two simple models of their surface inhomogeneity. Depending on the surface model and the location of the observer, the resulting lightcurve exhibits a departure from a purely sinusoidal curve that are observed for a fraction of white dwarfs. As a specific example, we fit our model to the observed phase-folded lightcurve of a fast-spinning white dwarf ZTF J190132.9+145808.7 (with the rotation period of 419s). We find that the size and obliquity angle of the spot responsible for the photometric variation are $\dts \approx 60^\circ$ and $\thetaS \approx 60^\circ$ or $90^\circ$, respectively, implying an interesting constraint on the surface distribution of the magnetic field on white dwarfs. △ Less

Submitted 7 November, 2022; v1 submitted 11 October, 2022; originally announced October 2022.

Comments: 25 pages, 9 figures, accepted for publication in PASJ

Lagrange vs. Lyapunov stability of hierarchical triple systems: dependence on the mutual inclination between inner and outer orbits

Authors: Toshinori Hayashi, Alessandro A. Trani, Yasushi Suto

Abstract: While there have been many studies examining the stability of hierarchical triple systems, the meaning of ``stability'' is somewhat vague and has been interpreted differently in previous literatures. The present paper focuses on ``Lagrange stability'', which roughly refers to the stability against the escape of a body from the system, or ``disruption'' of the triple system, in contrast to ``Lyapun… ▽ More While there have been many studies examining the stability of hierarchical triple systems, the meaning of ``stability'' is somewhat vague and has been interpreted differently in previous literatures. The present paper focuses on ``Lagrange stability'', which roughly refers to the stability against the escape of a body from the system, or ``disruption'' of the triple system, in contrast to ``Lyapunov-like stability'' that is related to the chaotic nature of the system dynamics. We compute the evolution of triple systems using direct $N$-body simulations up to $10^7 P_\mathrm{out}$, which is significantly longer than previous studies (with $P_\mathrm{out}$ being the initial orbital period of the outer body). We obtain the resulting disruption timescale $T_\mathrm{d}$ as a function of the triple orbital parameters with particular attention to the dependence on the mutual inclination between the inner and outer orbits, $i_\mathrm{mut}$. By doing so, we have clarified explicitly the difference between Lagrange and Lyapunov stabilities in astronomical triples. Furthermore, we find that the von Zeipel-Kozai-Lidov oscillations significantly destabilize inclined triples (roughly with $60^\circ < i_\mathrm{mut} < 150^\circ$) relative to those with $i_\mathrm{mut}=0^\circ$. On the other hand, retrograde triples with $i_\mathrm{mut}>160^\circ$ become strongly stabilized with much longer disruption timescales. We show the sensitivity of the normalized disruption timescale $T_\mathrm{d}/P_\mathrm{out}$ to the orbital parameters of triple system. The resulting $T_\mathrm{d}/P_\mathrm{out}$ distribution is practically more useful in a broad range of astronomical applications than the stability criterion based on the Lyapunov divergence. △ Less

Submitted 16 December, 2022; v1 submitted 18 September, 2022; originally announced September 2022.

Comments: 13 pages, 3 figures, ApJ, in press

Dynamical disruption timescales and chaotic behavior of hierarchical triple systems

Authors: Toshinori Hayashi, Alessandro A. Trani, Yasushi Suto

Abstract: We examine the stability of hierarchical triple systems using direct $N$-body simulations without adopting a secular perturbation approximation. We estimate their disruption timescales in addition to the mere stable/unstable criterion, with particular attention to the mutual inclination between the inner and outer orbits. First, we improve the fit to the dynamical stability criterion by \citet{Mar… ▽ More We examine the stability of hierarchical triple systems using direct $N$-body simulations without adopting a secular perturbation approximation. We estimate their disruption timescales in addition to the mere stable/unstable criterion, with particular attention to the mutual inclination between the inner and outer orbits. First, we improve the fit to the dynamical stability criterion by \citet{Mardling1999,Mardling2001} widely adopted in the previous literature. Especially, we find that that the stability boundary is very sensitive to the mutual inclination; coplanar retrograde triples and orthogonal triples are much more stable and unstable, respectively, than coplanar prograde triples. Next, we estimate the disruption timescales of triples satisfying the stability condition up to $10^9$ times the inner orbital period. The timescales follow the scaling predicted by \citet{Mushkin2020}, especially at high $e_\mathrm{out}$ where their random walk model is most valid. We obtain an improved empirical fit to the disruption timescales, which indicates that the coplanar retrograde triples are significantly more stable than the previous prediction. We furthermore find that the dependence on the mutual inclination can be explained by the energy transfer model based on a parabolic encounter approximation. We also show that the disruption timescales of triples are highly sensitive to the tiny change of the initial parameters, reflecting the genuine chaotic nature of the dynamics of those systems. △ Less

Submitted 5 September, 2022; v1 submitted 26 July, 2022; originally announced July 2022.

Comments: 33 pages, 19 figures, 1 table, ApJ, in press. Comments welcome

Meta-analysis of photometric and asteroseismic measurements of stellar rotation periods: the Lomb-Scargle periodogram, autocorrelation function, wavelet and rotational splitting analysis for 92 Kepler asteroseismic targets

Authors: Yuting Lu, Othman Benomar, Shoya Kamiaka, Yasushi Suto

Abstract: We perform photometric (the Lomb-Scargle periodogram, autocorrelation, and wavelet) and asteroseismic analyses of 92 Kepler solar-like main-sequence stars to understand the reliability of the measured stellar rotation periods. We focus on the 70 stars without reported stellar companions, and classify them into four groups according to the quarter-to-quarter variance of the Lomb-Scargle period and… ▽ More We perform photometric (the Lomb-Scargle periodogram, autocorrelation, and wavelet) and asteroseismic analyses of 92 Kepler solar-like main-sequence stars to understand the reliability of the measured stellar rotation periods. We focus on the 70 stars without reported stellar companions, and classify them into four groups according to the quarter-to-quarter variance of the Lomb-Scargle period and the precision of the asteroseismic period. We present detailed individual comparison among photometric and asteroseismic constraints for these stars. We find that most of our targets exhibit significant quarter-to-quarter variances in the photometric periods, suggesting that the photometrically estimated period should be regarded as a simplified characterization of the true stellar rotation period, especially under the presence of the latitudinal differential rotation. On the other hand, there are a fraction of stars with a relatively small quarter-to-quarter variance in the photometric periods, most of which have consistent values for asteroseismically and photometrically estimated rotation periods. We also identify over ten stars whose photometric and asteroseismic periods significantly disagree, which would be potentially interesting targets for further individual investigations. △ Less

Submitted 25 October, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Analytic model for photometric variation due to starspots on a differentially rotating star

Authors: Yasushi Suto, Shin Sasaki, Yuta Nakagawa, Othman Benomar

Abstract: We present an analytic model of the lightcurve variation for stars with non-evolving starspots on a differentially rotating surface. The Fourier coefficients of the harmonics of the rotation period are expressed in terms of the latitude of the spot, $\ell_{s}$, and the observer's line-of-sight direction, $\ell_{o}$, including the limb darkening effect. We generate different realizations of multi-s… ▽ More We present an analytic model of the lightcurve variation for stars with non-evolving starspots on a differentially rotating surface. The Fourier coefficients of the harmonics of the rotation period are expressed in terms of the latitude of the spot, $\ell_{s}$, and the observer's line-of-sight direction, $\ell_{o}$, including the limb darkening effect. We generate different realizations of multi-spots according to the model, and perform mock observations of the resulting lightcurve modulations. We discuss to what extent one can recover the properties of the spots and the parameters for the differential rotation law from the periodogram analysis. Although our analytical model neglects the evolution of spots on the stellar surface (dynamical motion, creation and annihilation), it provides a basic framework to interpret the photometric variation of stars, in particular from the existing Kepler data and the future space-born mission. It is also applicable to photometric modulations induced by rotation of various astronomical objects. △ Less

Submitted 10 May, 2022; originally announced May 2022.

Comments: 31 pages, 15 figures, accepted for publication in PASJ

Architecture of planetary systems predicted from protoplanetary disks observed with ALMA II: evolution outcomes and dynamical stability

Authors: Shijie Wang, Kazuhiro D. Kanagawa, Yasushi Suto

Abstract: Recent ALMA observations on disk substructures suggest the presence of embedded protoplanets in a large number disks. The primordial configurations of these planetary systems can be deduced from the morphology of the disk substructure and serve as initial conditions for numerical investigation of their future evolution. Starting from the initial configurations of 12 multi-planetary systems deduced… ▽ More Recent ALMA observations on disk substructures suggest the presence of embedded protoplanets in a large number disks. The primordial configurations of these planetary systems can be deduced from the morphology of the disk substructure and serve as initial conditions for numerical investigation of their future evolution. Starting from the initial configurations of 12 multi-planetary systems deduced from ALMA disks, we carried out two-stage N-body simulation to investigate the evolution of the planetary systems at the disk stage as well as the long term orbital stability after the disk dispersal. At the disk stage, our simulation includes both the orbital migration and pebble/gas accretion effects. We found a variety of planetary systems are produced and can be categorised into distant giant planets, Jupiter-like planets, Neptune-like planets and distant small planets. We found the disk stage evolution as well as the final configurations are sensitive to both the initial mass assignments and viscosity. After the disk stage, we implement only mutual gravity between star and planets and introduce stochastic perturbative forces. All systems are integrated for up to 10 Gyr to test their orbital stability. Most planetary systems are found to be stable for at least 10 Gyr with perturbative force in a reasonable range. Our result implies that a strong perturbation source such as stellar flybys is required to drive the planetary system unstable. We discuss the implications of our results on both the disk and planet observation, which may be confirmed by the next generation telescopes such as JWST and ngVLA. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 41 pages, 16 figures, accepted for publication in ApJ

Architecture of planetary systems predicted from protoplanetary disks observed with ALMA I: mass of the possible planets embedded in the dust gap

Authors: Shijie Wang, Kazuhiro D. Kanagawa, Yasushi Suto

Abstract: Recent ALMA observations have identified a variety of dust gaps in protoplanetary disks, which are commonly interpreted to be generated by unobserved planets. Predicting mass of such embedded planets is of fundamental importance in comparing those disk architectures with the observed diversity of exoplanets. The prediction, however, depends on the assumption that whether the same gap structure exi… ▽ More Recent ALMA observations have identified a variety of dust gaps in protoplanetary disks, which are commonly interpreted to be generated by unobserved planets. Predicting mass of such embedded planets is of fundamental importance in comparing those disk architectures with the observed diversity of exoplanets. The prediction, however, depends on the assumption that whether the same gap structure exists in the dust component alone or in the gas component as well. We assume a planet can only open a gap in the gas component when its mass exceeds the pebble isolation mass by considering the core accretion scenario. We then propose two criteria to distinguish if a gap is opened in the dust disk alone or the gas gap as well when observation data on the gas profile is not available. We apply the criteria to 35 disk systems with a total of 55 gaps compiled from previous studies, and classify each gap into four different groups. The classification of the observed gaps allows us to predict the mass of embedded planets in a consistent manner with the pebble isolation mass. We find that outer gaps are mostly dust alone, while inner gaps are more likely to be associated with a gas gap as well. The distribution of such embedded planets is very different from the architecture of the observed planetary systems, suggesting that the significant inward migration is required in their evolution. △ Less

Submitted 31 October, 2021; v1 submitted 10 October, 2021; originally announced October 2021.

Comments: 25 pages, 8 figures, accepted for publication in ApJ; sec 2 revised, figures and tables updated

Disentangling the stellar inclination of transiting planetary systems: fully analytic approach to the Rossiter-McLaughlin effect incorporating the stellar differential rotation

Authors: Shin Sasaki, Yasushi Suto

Abstract: The Rossiter-McLaughlin (RM) effect has been widely used to estimate the sky-projected spin-orbit angle, $λ$, of transiting planetary systems. Most of the previous analysis assume that the host stars are rigid rotators in which the amplitude of the RM velocity anomaly is proportional to $v_\star \sin i_\star$. When their latitudinal differential rotation is taken into account, one can break the de… ▽ More The Rossiter-McLaughlin (RM) effect has been widely used to estimate the sky-projected spin-orbit angle, $λ$, of transiting planetary systems. Most of the previous analysis assume that the host stars are rigid rotators in which the amplitude of the RM velocity anomaly is proportional to $v_\star \sin i_\star$. When their latitudinal differential rotation is taken into account, one can break the degeneracy, and determine separately the equatorial rotation velocity $v_\star$ and the inclination $i_{\star}$ of the host star. We derive a fully analytic approximate formula for the RM effect adopting a parameterized model for the stellar differential rotation. For those stars that exhibit the differential rotation similar to that of the Sun, the corresponding RM velocity modulation amounts to several m/s. We conclude that the latitudinal differential rotation offers a method to estimate $i_\star$, and thus the full spin-orbit angle $ψ$, from the RM data analysis alone. △ Less

Submitted 6 October, 2021; originally announced October 2021.

Comments: 20 pages, 11 figures, accepted for publication in PASJ

Third Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Sumiko Harasawa, Yuichi Harikane, Bau-Ching Hsieh, Hiroyuki Ikeda, Kei Ito, Ikuru Iwata, Tadayuki Kodama, Michitaro Koike, Mitsuru Kokubo, Yutaka Komiyama, Xiangchong Li, Yongming Liang, Yen-Ting Lin, Robert H. Lupton, Nate B Lust, Lauren A. MacArthur, Ken Mawatari , et al. (42 additional authors not shown)

Abstract: The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of ob… ▽ More The paper presents the third data release of Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP), a wide-field multi-band imaging survey with the Subaru 8.2m telescope. HSC-SSP has three survey layers (Wide, Deep, and UltraDeep) with different area coverages and depths, designed to address a wide array of astrophysical questions. This third release from HSC-SSP includes data from 278 nights of observing time and covers about 670 square degrees in all five broad-band filters at the full depth ($\sim26$~mag at $5σ$) in the Wide layer. If we include partially observed area, the release covers 1,470 square degrees. The Deep and UltraDeep layers have $\sim80\%$ of the originally planned integration times, and are considered done, as we have slightly changed the observing strategy in order to compensate for various time losses. There are a number of updates in the image processing pipeline. Of particular importance is the change in the sky subtraction algorithm; we subtract the sky on small scales before the detection and measurement stages, which has significantly reduced false detections. Thanks to this and other updates, the overall quality of the processed data has improved since the previous release. However, there are limitations in the data (for example, the pipeline is not optimized for crowded fields), and we encourage the user to check the quality assurance plots as well as a list of known issues before exploiting the data. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 30 August, 2021; originally announced August 2021.

Comments: 25 pages, 19 figures, submitted to PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

A new formation scenario of a counter-rotating circumstellar disk: spiral-arm accretion from a circumbinary disk in a triple protostar system

Authors: Daisuke Takaishi, Yusuke Tsukamoto, Yasushi Suto

Abstract: We present the evolution of rotational directions of circumstellar disks in a triple protostar system simulated from a turbulent molecular cloud core with no magnetic field. We find a new formation pathway of a counter-rotating circumstellar disk in such triple systems. The tertiary protostar forms via the circumbinary disk fragmentation and the initial rotational directions of all the three circu… ▽ More We present the evolution of rotational directions of circumstellar disks in a triple protostar system simulated from a turbulent molecular cloud core with no magnetic field. We find a new formation pathway of a counter-rotating circumstellar disk in such triple systems. The tertiary protostar forms via the circumbinary disk fragmentation and the initial rotational directions of all the three circumstellar disks are almost parallel to that of the orbital motion of the binary system. Their mutual gravito-hydrodynamical interaction for the subsequent $\sim10^4\thinspace\rm{yr}$ greatly disturbs the orbit of the tertiary, and the rotational directions of the tertiary disk and star are reversed due to the spiral-arm accretion of the circumbinary disk. The counter-rotation of the tertiary circumstellar disk continues to the end of the simulation ($\sim6.4\times10^4\thinspace\rm{yr}$ after its formation), implying that the counter-rotating disk is long-lived. This new formation pathway during the disk evolution in Class 0/I Young Stellar Objects possibly explains the counter-rotating disks recently discovered by ALMA. △ Less

Submitted 3 August, 2021; originally announced August 2021.

Comments: Accepted for publication in PASJ

Unveiling the architecture of a pulsar - binary black-hole triple system with pulsar arrival time analysis

Authors: Toshinori Hayashi, Yasushi Suto

Abstract: A large number of binary black holes (BBHs) with longer orbital periods are supposed to exist as progenitors of BBH mergers recently discovered with gravitational wave (GW) detectors. In our previous papers, we proposed to search for such BBHs in triple systems through the radial-velocity modulation of the tertiary orbiting star. If the tertiary is a pulsar, high precision and cadence observations… ▽ More A large number of binary black holes (BBHs) with longer orbital periods are supposed to exist as progenitors of BBH mergers recently discovered with gravitational wave (GW) detectors. In our previous papers, we proposed to search for such BBHs in triple systems through the radial-velocity modulation of the tertiary orbiting star. If the tertiary is a pulsar, high precision and cadence observations of its arrival time enable an unambiguous characterization of the pulsar -- BBH triples located at several kpc, which are inaccessible with the radial velocity of stars. The present paper shows that such inner BBHs can be identified through the short-term Rømer delay modulation, on the order of $10$ msec for our fiducial case, a triple consisting of $20~M_\odot$ BBH and $1.4~M_\odot$ pulsar with $P_\mathrm{in}=10$ days and $P_\mathrm{out}=100$ days. If the relativistic time delays are measured as well, one can determine basically all the orbital parameters of the triple. For instance, this method is applicable to inner BBHs of down to $\sim 1$ hr orbital periods if the orbital period of the tertiary pulsar is around several days. Inner BBHs with $\lesssim 1$ hr orbital period emit the GW detectable by future space-based GW missions including LISA, DECIGO, and BBO, and very short inner BBHs with sub-second orbital period can be even probed by the existing ground-based GW detectors. Therefore, our proposed methodology provides a complementary technique to search for inner BBHs in triples, if exist at all, in the near future. △ Less

Submitted 28 November, 2020; originally announced November 2020.

Comments: 24 pages, 8 figures, 3 tables, ApJ, in press

Search for Alignment of Disk Orientations in Nearby Star-Forming Regions: Lupus, Taurus, Upper Scorpius, $ρ$ Ophiuchi, and Orion

Authors: Masataka Aizawa, Yasushi Suto, Yoko Oya, Shiro Ikeda, Takeshi Nakazato

Abstract: Spatial correlations among proto-planetary disk orientations carry unique information on physics of multiple star formation processes. We select five nearby star-forming regions that comprise a number of proto-planetary disks with spatially-resolved images with ALMA and HST, and search for the mutual alignment of the disk axes. Specifically, we apply the Kuiper test to examine the statistical unif… ▽ More Spatial correlations among proto-planetary disk orientations carry unique information on physics of multiple star formation processes. We select five nearby star-forming regions that comprise a number of proto-planetary disks with spatially-resolved images with ALMA and HST, and search for the mutual alignment of the disk axes. Specifically, we apply the Kuiper test to examine the statistical uniformity of the position angle (PA: the angle of the major axis of the projected disk ellipse measured counter-clockwise from the north) distribution. The disks located in the star-forming regions, except the Lupus clouds, do not show any signature of the alignment, supporting the random orientation. Rotational axes of 16 disks with spectroscopic measurement of PA in the Lupus III cloud, a sub-region of the Lupus field, however, exhibit a weak and possible departure from the random distribution at a $2σ$ level, and the inclination angles of the 16 disks are not uniform as well. Furthermore, the mean direction of the disk PAs in the Lupus III cloud is parallel to the direction of its filament structure, and approximately perpendicular to the magnetic field direction. We also confirm the robustness of the estimated PAs in the Lupus clouds by comparing the different observations and estimators based on three different methods including sparse modeling. The absence of the significant alignment of the disk orientation is consistent with the turbulent origin of the disk angular momentum. Further observations are required to confirm/falsify the possible disk alignment in the Lupus III cloud. △ Less

Submitted 7 July, 2020; originally announced July 2020.

Comments: 26 pages, 8 figures, accepted for publication in ApJ

Beyond a pale blue dot : how to search for possible bio-signatures on earth-like planets

Authors: Yasushi Suto

Abstract: The Earth viewed from outside the Solar system would be identified merely like a pale blue dot, as coined by Carl Sagan. In order to detect possible signatures of the presence of life on a second earth among several terrestrial planets discovered in a habit-able zone, one has to develop and establish a methodology to characterize the planet as something beyond a mere pale blue dot. We pay particul… ▽ More The Earth viewed from outside the Solar system would be identified merely like a pale blue dot, as coined by Carl Sagan. In order to detect possible signatures of the presence of life on a second earth among several terrestrial planets discovered in a habit-able zone, one has to develop and establish a methodology to characterize the planet as something beyond a mere pale blue dot. We pay particular attention to the periodic change of the color of the dot according to the rotation of the planet. Because of the large-scale inhomogeneous distribution of the planetary surface, the reflected light of the dot comprises different color components corresponding to land, ocean, ice, and cloud that cover the surface of the planet. If we decompose the color of the dot into several principle components, in turn, one can identify the presence of the different surface components. Furthermore, the vegetation on the earth is known to share a remarkable reflection signature; the reflection becomes significantly enhanced at wave-lengths longer than 760nm, which is known as a red-edge of the vegetation. If one can identify the corresponding color signature in a pale blue dot, it can be used as a unique probe of the presence of life. I will describe the feasibility of the methodology for future space missions, and consider the direction towards astrobiology from an astrophysicist's point of view. △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: 11 pages, 5 figures, published in Yamagishi A., Kakegawa T., Usui T. (eds) Astrobiology. Springer, Singapore (2019)

Obliquity of an Earth-like planet from frequency modulation of its direct imaged lightcurve: mock analysis from general circulation model simulation

Authors: Yuta Nakagawa, Takanori Kodama, Masaki Ishiwatari, Hajime Kawahara, Yasushi Suto, Yoshiyuki O. Takahashi, George L. Hashimoto, Kiyoshi Kuramoto, Kensuke Nakajima, Shin-ichi Takehiro, Yoshi-Yuki Hayashi

Abstract: Direct-imaging techniques of exoplanets have made significant progress recently, and will eventually enable to monitor photometric and spectroscopic signals of earth-like habitable planets in the future. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produc… ▽ More Direct-imaging techniques of exoplanets have made significant progress recently, and will eventually enable to monitor photometric and spectroscopic signals of earth-like habitable planets in the future. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produce different cloud patterns by performing a series of GCM (General Circulation Model) simulation runs using a set of parameters relevant for our Earth. Then we use the simulated photometric lightcurves to compute their frequency modulation due to the planetary spin-orbit coupling over an entire orbital period, and attempt to see to what extent one can estimate the obliquity of an Earth-twin. We find that it is possible to estimate the obliquity of an Earth-twin within the uncertainty of several degrees with a dedicated 4 m space telescope at 10 pc away from the system if the stellar flux is completely blocked. While our conclusion is based on several idealized assumptions, a frequency modulation of a directly-imaged earth-like planet offers a unique methodology to determine its obliquity. △ Less

Submitted 19 June, 2020; originally announced June 2020.

Comments: 29 pages, 18 figures, accepted for publication in ApJ

Radial-velocity variation of a tertiary star orbiting a binary black hole in coplanar and non-coplanar triples: short- and long-term anomalous behavior

Authors: Toshinori Hayashi, Yasushi Suto

Abstract: A number of ongoing surveys are likely to discover star-black hole binaries in our Galaxy in the near future. A fraction of them may be triple systems comprising an inner binary, instead of a single black hole, which might be progenitors of binary black holes (BBHs) routinely discovered now from the gravitational wave. We extend our previous proposal to locate inner BBHs from the short-term radial… ▽ More A number of ongoing surveys are likely to discover star-black hole binaries in our Galaxy in the near future. A fraction of them may be triple systems comprising an inner binary, instead of a single black hole, which might be progenitors of binary black holes (BBHs) routinely discovered now from the gravitational wave. We extend our previous proposal to locate inner BBHs from the short-term radial-velocity (RV) variation of a tertiary star in coplanar triples, and we consider noncoplanar triples and their long-term RV variations as well. Specifically, we assume coplanar and noncoplanar triples with an inner BBH of the total mass $20~M_\odot$, whose outer and inner orbital periods are 80 days and 10 days, respectively. We perform a series of N-body simulations and compare the results with analytic approximate solutions based on quadrupole perturbation theory. For coplanar triples, the pericenter shift of the outer star can be used to detect the hidden inner BBH. For noncoplanar triples, the total RV semi-amplitude of the outer star is modulated periodically on the order of $100$km/s due to its precession over roughly the Kozai-Lidov oscillation timescale. Such long-term modulations would be detectable within a decade, independent of the short-term RV variations on the order of of $100$ m/s at roughly twice the orbital frequency of the inner binary. Thus the RV monitoring of future star-black hole binary candidates offers a promising method for searching for their inner hidden BBHs in optical bands. △ Less

Submitted 20 April, 2022; v1 submitted 29 May, 2020; originally announced June 2020.

Comments: 30 pages, 14 figures, 3 tables; this is the updated version of arXiv:2006.02210; Fig.8, Fig.9, Fig.10, and Fig.11 are revised according to the erratum, ApJ, 2022, 929, 109

Journal ref: ApJ, 2020, 897, 29; Erratum, ApJ, 2022, 929, 109

Shapes and alignments of dark matter haloes and their brightest cluster galaxies in 39 strong lensing clusters

Authors: Taizo Okabe, Masamune Oguri, Sébastien Peirani, Yasushi Suto, Yohan Dubois, Christophe Pichon, Tetsu Kitayama, Shin Sasaki, Takahiro Nishimichi

Abstract: We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light… ▽ More We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter. △ Less

Submitted 23 May, 2020; originally announced May 2020.

Comments: 18 pages, 16 figures, 6 Tables, accepted by MNRAS

Report number: YITP-20-76

Architecture of three-planet systems predicted from the observed protoplanetary disk of HL Tau

Authors: Shijie Wang, Kazuhiro D. Kanagawa, Toshinori Hayashi, Yasushi Suto

Abstract: A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolut… ▽ More A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolution of planets within the protoplanetary disk and after the disk dispersal. Before the disk dispersal, our N-body simulations include both migration and mass-growth of the planet coupled with evolution of the disk. By varying the disk parameters, we produce a variety of widely-separated planetary systems consisting of three super-Jupiters at the end of disk dispersal. We found the outer planet is more massive than the inner one, and the migration of the innermost planet is inefficient due to the accretion of outer planet(s). We also showed how the final configuration and the final planetary mass depend on disk parameters. The migration is found to be convergent and no planet-pair has a period ratio less than 2. After the disk dispersal, we switch to pure gravitational N-body simulations and integrate the orbits up to 10 Gyr. Most simulated systems remain stable for at least 10 Gyr. We discuss implications of our result in terms of the observed widely-separated planetary systems HR 8799 and PDS 70. △ Less

Submitted 19 February, 2020; originally announced February 2020.

Comments: 16 pages, 9 figures. Accepted for publication in the ApJ

Star-disk alignment in the protoplanetary disks: SPH simulation of the collapse of turbulent molecular cloud cores

Authors: Daisuke Takaishi, Yusuke Tsukamoto, Yasushi Suto

Abstract: We perform a series of three-dimensional smoothed particle hydrodynamics (SPH) simulations to study the evolution of the angle between the protostellar spin and the protoplanetary disk rotation axes (the star-disk angle $ψ_{\rm sd}$) in turbulent molecular cloud cores. While $ψ_{\rm sd}$ at the protostar formation epoch exhibits broad distribution up to $\sim 130^{\circ}$, $ψ_{\rm sd}$ decreases (… ▽ More We perform a series of three-dimensional smoothed particle hydrodynamics (SPH) simulations to study the evolution of the angle between the protostellar spin and the protoplanetary disk rotation axes (the star-disk angle $ψ_{\rm sd}$) in turbulent molecular cloud cores. While $ψ_{\rm sd}$ at the protostar formation epoch exhibits broad distribution up to $\sim 130^{\circ}$, $ψ_{\rm sd}$ decreases ($\lesssim 20^{\circ}$) in a timescale of $\sim 10^{4}$ yr. This timescale of the star-disk alignment, $t_{\rm alignment}$, corresponds basically to the mass doubling time of the central protostar, in which the protostar forgets its initial spin direction due to the mass accretion from the disk. Values of $ψ_{\rm sd}$ both at $t=10^2$ yr and $t=10^5$ yr after the protostar formation are independent of the ratios of thermal and turbulent energies to gravitational energy of the initial cloud cores: $α=E_{\rm thermal}/|E_{\rm gravity}|$ and $γ_{\rm turb}=E_{\rm turbulence}/|E_{\rm gravity}|$. We also find that a warped disk is possibly formed by the turbulent accretion flow from the circumstellar envelope. △ Less

Submitted 15 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS (16 pages, 16 figures, 2 tables)

Cosmological evolution of orientations of cluster-sized dark matter haloes and their central galaxies in the Horizon-AGN simulation

Authors: Taizo Okabe, Takahiro Nishimichi, Masamune Oguri, Sébastien Peirani, Tetsu Kitayama, Shin Sasaki, Yasushi Suto, Christophe Pichon, Yohan Dubois

Abstract: It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We t… ▽ More It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: 14 pages, 11 figures

Second Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Yusra AlSayyad, Makoto Ando, Robert Armstrong, James Bosch, Eiichi Egami, Hisanori Furusawa, Junko Furusawa, Andy Goulding, Yuichi Harikane, Chiaki Hikage, Paul T. P. Ho, Bau-Ching Hsieh, Song Huang, Hiroyuki Ikeda, Masatoshi Imanishi, Kei Ito, Ikuru Iwata, Anton T. Jaelani, Ryota Kakuma, Kojiro Kawana, Satoshi Kikuta, Umi Kobayashi, Michitaro Koike, Yutaka Komiyama , et al. (40 additional authors not shown)

Abstract: This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially ob… ▽ More This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 22 August, 2019; v1 submitted 29 May, 2019; originally announced May 2019.

Comments: 26 pages, 23 figures, accepted for publication in PASJ. Data available at https://hsc-release.mtk.nao.ac.jp/

A strategy to search for an inner binary black hole from the motion of the tertiary star

Authors: Toshinori Hayashi, Shijie Wang, Yasushi Suto

Abstract: There are several on-going projects to search for stars orbiting around an invisible companion. A fraction of such candidates may be a triple, instead of a binary, consisting of an inner binary black hole (BBH) and an outer orbiting star. In this paper, we propose a methodology to search for a signature of such an inner BBH, possibly a progenitor of gravitational-wave sources discovered by {\it LI… ▽ More There are several on-going projects to search for stars orbiting around an invisible companion. A fraction of such candidates may be a triple, instead of a binary, consisting of an inner binary black hole (BBH) and an outer orbiting star. In this paper, we propose a methodology to search for a signature of such an inner BBH, possibly a progenitor of gravitational-wave sources discovered by {\it LIGO}, from the precise radial velocity (RV) follow-up of the outer star. We first describe a methodology using an existing approximate RV formula for coplanar circular triples. We apply this method and constrain the parameters of a possible inner binary objects in 2M05215658+4359220, which consists of a red giant and an unseen companion. Next we consider co-planar but non-circular triples. We compute numerically the RV variation of a tertiary star orbiting around an inner BBH, generate mock RV curves, and examine the feasibility of the BBH detection for our fiducial models. We conclude that the short-cadence RV monitoring of a star-BH binary provides an interesting and realistic method to constrain and/or search for possible inner BBHs. Indeed a recent discovery of a star--BH binary system LB-1 may imply that there are a large number of such unknown objects in our Galaxy, which are ideal targets for the methodology proposed here. △ Less

Submitted 19 January, 2020; v1 submitted 16 May, 2019; originally announced May 2019.

Comments: 20 pages, 10 figures, 5 tables Substantially revised and now accepted for publication in ApJ

Asteroseismic determination of the stellar rotation period of the Kepler transiting planetary systems and its implications for the spin-orbit architecture

Authors: Yasushi Suto, Shoya Kamiaka, Othman Benomar

Abstract: We measure the rotation periods of 19 stars in the {\it Kepler} transiting planetary systems, $P_{\rm rot, astero}$ from asteroseismology and $P_{\rm rot, phot}$ from photometric variation of their lightcurve. Two stars exhibit two clear peaks in the Lomb-Scargle periodogram, neither of which agrees with the seismic rotation period. Other four systems do not show any clear peak, whose stellar rota… ▽ More We measure the rotation periods of 19 stars in the {\it Kepler} transiting planetary systems, $P_{\rm rot, astero}$ from asteroseismology and $P_{\rm rot, phot}$ from photometric variation of their lightcurve. Two stars exhibit two clear peaks in the Lomb-Scargle periodogram, neither of which agrees with the seismic rotation period. Other four systems do not show any clear peak, whose stellar rotation period is impossible to estimate reliably from the photometric variation; their stellar equators may be significantly inclined with respect to the planetary orbital plane. For the remaining 13 systems, $P_{\rm rot, astero}$ and $P_{\rm rot, phot}$ agree within 30\%. Interestingly, three out of the 13 systems are in the spin-orbit resonant state in which $P_{\rm orb, b}/P_{\rm rot, astero} \approx 1$ with $P_{\rm orb, b}$ being the orbital period of the inner-most planet of each system. The corresponding chance probability is ($0.2$-$4.7$) \% based on the photometric rotation period data for 464 {\it Kepler} transiting planetary systems. While further analysis of stars with reliable rotation periods is required to examine the statistical significance, the spin-orbit resonance between the star and planets, if confirmed, have important implications for the star-planet tidal interaction, in addition to the origin of the spin-orbit (mis-)alignment of transiting planetary systems. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 25 pages, 14 figures, 3 tables, accepted for publication in AJ

The misaligned orbit of the Earth-sized planet Kepler-408b

Authors: Shoya Kamiaka, Othman Benomar, Yasushi Suto, Fei Dai, Kento Masuda, Joshua N. Winn

Abstract: Kepler-408 is one of the 33 planet-hosting {\it Kepler} stars for which asteroseismology has been used to investigate the orientation of the stellar rotation axis relative to the planetary orbital plane. The transiting "hot Earth," Kepler-408b, has an orbital period of 2.5 days and a radius of $0.86$~$R_\oplus$, making it much smaller than the planets for which spin-orbit alignment has been studie… ▽ More Kepler-408 is one of the 33 planet-hosting {\it Kepler} stars for which asteroseismology has been used to investigate the orientation of the stellar rotation axis relative to the planetary orbital plane. The transiting "hot Earth," Kepler-408b, has an orbital period of 2.5 days and a radius of $0.86$~$R_\oplus$, making it much smaller than the planets for which spin-orbit alignment has been studied using the Rossiter-McLaughlin effect. Because conflicting asteroseismic results have been reported in the literature, we undertake a thorough re-appraisal of this system and perform numerous checks for consistency and robustness. We find that the conflicting results are due to the different models for the low-frequency noise in the power spectrum. A careful treatment of the background noise resolves these conflicts, and shows that the stellar inclination is $\is=42^{+5}_{-4}$ degrees. Kepler-408b is, by far, the smallest planet known to have a significantly misaligned orbit. △ Less

Submitted 9 February, 2019; v1 submitted 6 February, 2019; originally announced February 2019.

Comments: 20 pages, 12 figures, accepted for publication in AJ

Cool core disturbed: Observational evidence for coexistence of sub-sonic sloshing gas and stripped shock-heated gas around the core of RX J1347.5-1145

Authors: Shutaro Ueda, Tetsu Kitayama, Masamune Oguri, Eiichiro Komatsu, Takuya Akahori, Daisuke Iono, Takumi Izumi, Ryohei Kawabe, Kotaro Kohno, Hiroshi Matsuo, Naomi Ota, Yasushi Suto, Shigehisa Takakuwa, Motokazu Takizawa, Takahiro Tsutsumi, Kohji Yoshikawa

Abstract: RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern… ▽ More RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5σ. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL). △ Less

Submitted 28 August, 2018; originally announced August 2018.

Comments: 18 pages, 10 figures, accepted for publication in ApJ

Reliability of stellar inclination estimated from asteroseismology: analytical criteria, mock simulations and Kepler data analysis

Authors: Shoya Kamiaka, Othman Benomar, Yasushi Suto

Abstract: Advances in asteroseismology of solar-like stars, now provide a unique method to estimate the stellar inclination $i_{\star}$. This enables to evaluate the spin-orbit angle of transiting planetary systems, in a complementary fashion to the Rossiter-McLaughlin effect, a well-established method to estimate the projected spin-orbit angle $λ$. Although the asteroseismic method has been broadly applied… ▽ More Advances in asteroseismology of solar-like stars, now provide a unique method to estimate the stellar inclination $i_{\star}$. This enables to evaluate the spin-orbit angle of transiting planetary systems, in a complementary fashion to the Rossiter-McLaughlin effect, a well-established method to estimate the projected spin-orbit angle $λ$. Although the asteroseismic method has been broadly applied to the Kepler data, its reliability has yet to be assessed intensively. In this work, we evaluate the accuracy of $i_{\star}$ from asteroseismology of solar-like stars using 3000 simulated power spectra. We find that the low signal-to-noise ratio of the power spectra induces a systematic under-estimate (over-estimate) bias for stars with high (low) inclinations. We derive analytical criteria for the reliable asteroseismic estimate, which indicates that reliable measurements are possible in the range of $20^\circ \lesssim i_{\star} \lesssim 80^\circ$ only for stars with high signal-to-noise ratio. We also analyse and measure the stellar inclination of 94 Kepler main-sequence solar-like stars, among which 33 are planetary hosts. According to our reliability criteria, a third of them (9 with planets, 22 without) have accurate stellar inclination. Comparison of our asteroseismic estimate of $v\sin{i_{\star}}$ against spectroscopic measurements indicates that the latter suffers from a large uncertainty possibly due to the modeling of macro-turbulence, especially for stars with projected rotation speed $v\sin{i_{\star}} \lesssim 5$ km/s. This reinforces earlier claims, and the stellar inclination estimated from the combination of measurements from spectroscopy and photometric variation for slowly rotating stars needs to be interpreted with caution. △ Less

Submitted 18 May, 2018; originally announced May 2018.

Comments: 17 pages, 10 figures, 3 tables, accepted for publication in MNRAS

Projected alignment of non-sphericities of stellar, gas, and dark matter distributions in galaxy clusters: analysis of the Horizon-AGN simulation

Authors: Taizo Okabe, Takahiro Nishimichi, Masamune Oguri, Sébastien Peirani, Tetsu Kitayama, Shin Sasaki, Yasushi Suto

Abstract: While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various component… ▽ More While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters. △ Less

Submitted 28 June, 2018; v1 submitted 24 April, 2018; originally announced April 2018.

Comments: 21 pages, 12 figures, 2 tables, accepted MNRAS

Systematic Search for Rings around Kepler Planet Candidates: Constraints on Ring Size and Occurrence Rate

Authors: Masataka Aizawa, Kento Masuda, Hajime Kawahara, Yasushi Suto

Abstract: We perform a systematic search for rings around 168 Kepler planet candidates with sufficient signal-to-noise ratios that are selected from all the short-cadence data. We fit ringed and ringless models to their lightcurves, and compare the fitting results to search for the signatures of planetary rings. First, we identify 29 tentative systems, for which the ringed models exhibit statistically signi… ▽ More We perform a systematic search for rings around 168 Kepler planet candidates with sufficient signal-to-noise ratios that are selected from all the short-cadence data. We fit ringed and ringless models to their lightcurves, and compare the fitting results to search for the signatures of planetary rings. First, we identify 29 tentative systems, for which the ringed models exhibit statistically significant improvement over the ringless models. The lightcurves of those systems are individually examined, but we are not able to identify any candidate that indicates evidence for rings. In turn, we find out several mechanisms of false-positives that would produce ring-like signals, and the null detection enables us to place upper limits on the size of rings. Furthermore, assuming the tidal alignment between axes of the planetary rings and orbits, we conclude that the occurrence rate of rings larger than twice the planetary radius is less than 15 percent. Even though the majority of our targets are short-period planets, our null detection provides statistical and quantitative constraints on largely uncertain theoretical models of origin, formation, and evolution of planetary rings. △ Less

Submitted 24 March, 2018; originally announced March 2018.

Comments: 36 pages, 13 figures, 3 tables, accepted for publication in the Astronomical Journal

The Hyper Suprime-Cam SSP Survey: Overview and Survey Design

Authors: H. Aihara, N. Arimoto, R. Armstrong, S. Arnouts, N. A. Bahcall, S. Bickerton, J. Bosch, K. Bundy, P. L. Capak, J. H. H. Chan, M. Chiba, J. Coupon, E. Egami, M. Enoki, F. Finet, H. Fujimori, S. Fujimoto, H. Furusawa, J. Furusawa, T. Goto, A. Goulding, J. P. Greco, J. E. Greene, J. E. Gunn, T. Hamana , et al. (118 additional authors not shown)

Abstract: Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), w… ▽ More Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,σ$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey. △ Less

Submitted 15 March, 2018; v1 submitted 19 April, 2017; originally announced April 2017.

Comments: 14 pages, 7 figures, 5 tables. Corrected for a typo in the coordinates of HSC-Wide spring equatorial field in Table 5

Journal ref: Publications of the Astronomical Society of Japan, Volume 70, Issue SP1, S4 (2018)

First Data Release of the Hyper Suprime-Cam Subaru Strategic Program

Authors: Hiroaki Aihara, Robert Armstrong, Steven Bickerton, James Bosch, Jean Coupon, Hisanori Furusawa, Yusuke Hayashi, Hiroyuki Ikeda, Yukiko Kamata, Hiroshi Karoji, Satoshi Kawanomoto, Michitaro Koike, Yutaka Komiyama, Robert H. Lupton, Sogo Mineo, Hironao Miyatake, Satoshi Miyazaki, Tomoki Morokuma, Yoshiyuki Obuchi, Yukie Oishi, Yuki Okura, Paul A. Price, Tadafumi Takata, Manobu M. Tanaka, Masayuki Tanaka , et al. (83 additional authors not shown)

Abstract: The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This… ▽ More The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/. △ Less

Submitted 28 July, 2017; v1 submitted 27 February, 2017; originally announced February 2017.

Comments: 34 pages, 20 figures, 7 tables, moderate revision, accepted for publication in PASJ

Towards Detection of Exoplanetary Rings Via Transit Photometry: Methodology and a Possible Candidate

Authors: Masataka Aizawa, Sho Uehara, Kento Masuda, Hajime Kawahara, Yasushi Suto

Abstract: Detection of a planetary ring of exoplanets remains as one of the most attractive but challenging goals in the field. We present a methodology of a systematic search for exoplanetary rings via transit photometry of long-period planets. The methodology relies on a precise integration scheme we develop to compute a transit light curve of a ringed planet. We apply the methodology to 89 long-period pl… ▽ More Detection of a planetary ring of exoplanets remains as one of the most attractive but challenging goals in the field. We present a methodology of a systematic search for exoplanetary rings via transit photometry of long-period planets. The methodology relies on a precise integration scheme we develop to compute a transit light curve of a ringed planet. We apply the methodology to 89 long-period planet candidates from the Kepler data so as to estimate, and/or set upper limits on, the parameters of possible rings. While a majority of our samples do not have a sufficiently good signal-to-noise ratio for meaningful constraints on ring parameters, we find that six systems with a higher signal-to-noise ratio are inconsistent with the presence of a ring larger than 1.5 times the planetary radius assuming a grazing orbit and a tilted ring. Furthermore, we identify five preliminary candidate systems whose light curves exhibit ring-like features. After removing four false positives due to the contamination from nearby stars, we identify KIC 10403228 as a reasonable candidate for a ringed planet. A systematic parameter fit of its light curve with a ringed planet model indicates two possible solutions corresponding to a Saturn-like planet with a tilted ring. There also remain other two possible scenarios accounting for the data; a circumstellar disk and a hierarchical triple. Due to large uncertain factors, we cannot choose one specific model among the three. △ Less

Submitted 27 February, 2017; originally announced February 2017.

Comments: 45 pages, 20 figures, accepted for publication in AJ

Possible Outcomes of Coplanar High-eccentricity Migration: Hot Jupiters, Close-in Super-Earths, and Counter-orbiting Planets

Authors: Yuxin Xue, Kento Masuda, Yasushi Suto

Abstract: We investigate the formation of close-in planets in near-coplanar eccentric hierarchical triple systems via the secular interaction between an inner planet and an outer perturber (Coplanar High-eccentricity Migration, CHEM). We generalize the previous work on the analytical condition for successful CHEM for point masses interacting only through gravity by taking into account the finite mass effect… ▽ More We investigate the formation of close-in planets in near-coplanar eccentric hierarchical triple systems via the secular interaction between an inner planet and an outer perturber (Coplanar High-eccentricity Migration, CHEM). We generalize the previous work on the analytical condition for successful CHEM for point masses interacting only through gravity by taking into account the finite mass effect of the inner planet. We find that efficient CHEM requires that the systems should have m_1<<m_0 and m_1<<m_2. In addition to the gravity for point masses, we examine the importance of the short-range forces, and provide an analytical estimate of the migration time scale. We perform a series of numerical simulations in CHEM for systems consisting of a sun-like central star, giant gas inner planet and planetary outer perturber, including the short-range forces and stellar and planetary dissipative tides. We find that most of such systems end up with a tidal disruption, a small fraction of the systems produce prograde hot Jupiters (HJs), but no retrograde one. In addition, we extend CHEM to super-Earth mass range, and show that the formation of close-in super-Earths in prograde orbits is also possible. Finally, we carry out CHEM simulation for the observed hierarchical triple and counter-orbiting HJ systems. We find that CHEM can explain a part of the former systems, but it is generally very difficult to reproduce counter-orbiting HJ systems. △ Less

Submitted 19 December, 2016; originally announced December 2016.

Comments: Accepted for publication in ApJ

Projected Axis Ratios of Galaxy Clusters in the Horizon-AGN Simulation: Impact of Baryon Physics and Comparison with Observations

Authors: Daichi Suto, Sebastien Peirani, Yohan Dubois, Tetsu Kitayama, Takahiro Nishimichi, Shin Sasaki, Yasushi Suto

Abstract: We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5E13 Msun from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the AGN feedback. We find that the baryonic physics a… ▽ More We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5E13 Msun from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the AGN feedback. We find that the baryonic physics around the central region of galaxy clusters significantly affects the non-sphericity of dark matter distribution even beyond the central region, approximately up to the half of the virial radius. Therefore it is very difficult to predict the the probability density function (PDF) of the projected axis ratio of XSB from dark-matter only N-body simulations as attempted in previous studies. Indeed we find that the PDF derived from our simulated clusters exhibits much better agreement with that from the observed X-ray clusters. This indicates that our present methodology to estimate the non-sphericity directly from the Horizon simulation is useful and promising. Further improvements in both numerical modeling and observational data will establish the non-sphericity of clusters as a cosmological test complementary to more conventional statistics based on spherically averaged quantities. △ Less

Submitted 16 November, 2016; originally announced November 2016.

Comments: 23 pages, 31 figures, accepted for publication in PASJ

Evolution and Statistics of Non-Sphericity of Dark Matter Halos from Cosmological N-Body Simulation

Authors: Daichi Suto, Tetsu Kitayama, Takahiro Nishimichi, Shin Sasaki, Yasushi Suto

Abstract: We revisit the non-sphericity of cluster-mass scale halos from cosmological N-body simulation on the basis of triaxial modelling. In order to understand the difference between the simulation results and the conventional ellipsoidal collapse model (EC), we first consider the evolution of individual simulated halos. The major difference between EC and the simulation becomes appreciable after the tur… ▽ More We revisit the non-sphericity of cluster-mass scale halos from cosmological N-body simulation on the basis of triaxial modelling. In order to understand the difference between the simulation results and the conventional ellipsoidal collapse model (EC), we first consider the evolution of individual simulated halos. The major difference between EC and the simulation becomes appreciable after the turn-around epoch. Moreover, it is sensitive to the individual evolution history of each halo. Despite such strong dependence on individual halos, the resulting nonsphericity of halos exhibits weak but robust mass dependence in a statistical fashion; massive halos are more spherical up to the turn-around, but gradually become less spherical by z = 0. This is clearly inconsistent with the EC prediction; massive halos are usually more spherical. In addition, at z=0, inner regions of the halos are less spherical than outer regions, i.e., the density distribution inside the halos is highly inhomogeneous and therefore not self-similar. Since most of previous fitting formulae for the PDF of axis ratio of triaxial ellipsoids have been constructed under the self-similarity assumption, they are not accurate. Indeed, we compute the PDF of projected axis ratio a1/a2 directly from the simulation data without the self-similarity assumption, and find that it is very sensitive to the assumption. The latter needs to be carefully taken into account in direct comparison with observations, and therefore we provide an empirical fitting formula for the PDF of a1/a2. Our preliminary analysis suggests that the derived PDF of a1/a2 roughly agrees with the current weak-lensing observations. More importantly, the present results will be useful in future exploration of the non-sphericity of clusters in X-ray and optical observations. △ Less

Submitted 23 August, 2016; originally announced August 2016.

Comments: 35 pages, 12 figures, 2 tables, accepted for publication in PASJ

Prime Focus Spectrograph (PFS) for the Subaru Telescope: Overview, recent progress, and future perspectives

Authors: Naoyuki Tamura, Naruhisa Takato, Atsushi Shimono, Yuki Moritani, Kiyoto Yabe, Yuki Ishizuka, Akitoshi Ueda, Yukiko Kamata, Hrand Aghazarian, Stephane Arnouts, Gabriel Barban, Robert H. Barkhouser, Renato C. Borges, David F. Braun, Michael A. Carr, Pierre-Yves Chabaud, Yin-Chang Chang, Hsin-Yo Chen, Masashi Chiba, Richard C. Y. Chou, You-Hua Chu, Judith G. Cohen, Rodrigo P. de Almeida, Antonio C. de Oliveira, Ligia S. de Oliveira , et al. (75 additional authors not shown)

Abstract: PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simult… ▽ More PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward. △ Less

Submitted 3 August, 2016; originally announced August 2016.

Comments: 17 pages, 10 figures. Proceeding of SPIE Astronomical Telescopes and Instrumentation 2016

The Sunyaev-Zel'dovich Effect at Five Arc-seconds: RXJ1347.5-1145 Imaged by ALMA

Authors: T. Kitayama, S. Ueda, S. Takakuwa, T. Tsutsumi, E. Komatsu, T. Akahori, D. Iono, T. Izumi, R. Kawabe, K. Kohno, H. Matsuo, N. Ota, Y. Suto, M. Takizawa, K. Yoshikawa

Abstract: We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for ima… ▽ More We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for imaging the SZE, while retaining extended signals out to 40 arc-seconds. The 1-sigma statistical sensitivity of the image is 0.017 mJy/beam or 0.12 mK_CMB at the 5 arc-second full width at half maximum. The SZE image shows a good agreement with an electron pressure map reconstructed independently from the X-ray data and offers a new probe of the small-scale structure of the intracluster medium. Our results demonstrate that ALMA is a powerful instrument for imaging the SZE in compact galaxy clusters with unprecedented angular resolution and sensitivity. As the first report on the detection of the SZE by ALMA, we present detailed analysis procedures including corrections for the missing flux, to provide guiding methods for analyzing and interpreting future SZE images by ALMA. △ Less

Submitted 29 July, 2016; originally announced July 2016.

Comments: 20 pages, 13 figures. Accepted for publication in PASJ

Difficulty in Formation of Counter-orbiting Hot Jupiters from Near-coplanar Hierarchical Triple Systems: A Sub-stellar Perturber

Authors: Yuxin Xue, Yasushi Suto

Abstract: Among a hundred transiting planets with a measured projected spin-orbit angle λ, several systems are suggested to be counter-orbiting. While they may be due to the projection effect, the mechanism to produce a counter-orbiting planet is not established. A promising scenario for the counter-orbiting planets is the extreme eccentricity evolution in near-coplanar hierarchical triple systems with ecce… ▽ More Among a hundred transiting planets with a measured projected spin-orbit angle λ, several systems are suggested to be counter-orbiting. While they may be due to the projection effect, the mechanism to produce a counter-orbiting planet is not established. A promising scenario for the counter-orbiting planets is the extreme eccentricity evolution in near-coplanar hierarchical triple systems with eccentric inner and outer or- bits. We examine this scenario in detail by performing a series of systematic numerical simulations, and consider the possibility of forming hot Jupiters, especially counter- orbiting one under this mechanism with a distant sub-stellar perturber. We incorporate quadrupole and octupole secular gravitational interaction between the two orbits, and also short-range forces (correction for general relativity, star and inner planetary tide and rotational distortion) simultaneously. We find that most of systems are tidally disrupted and that a small fraction of survived planets turns out to be prograde. The formation of counter-orbiting hot Jupiters in this scenario is possible only in a very restricted parameter region, and thus very unlikely in practice. △ Less

Submitted 14 February, 2016; originally announced February 2016.

Comments: accepted for publication in ApJ

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

Authors: Kento Masuda, Yasushi Suto

Abstract: 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… ▽ More 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. △ Less

Submitted 8 February, 2016; originally announced February 2016.

Comments: 9 pages, 4 figures, accepted for publication in PASJ

Image Stacking Analysis of SDSS Galaxies with AKARI Far-Infrared Surveyor Maps at 65μm, 90μm, and 140μm

Authors: Taizo Okabe, Toshiya Kashiwagi, Yasushi Suto, Shuji Matsuura, Yasuo Doi, Satoshi Takita, Takafumi Ootsubo

Abstract: We perform image stacking analysis of Sloan Digital Sky Survey (SDSS) photometric galaxies over the AKARI Far-Infrared Surveyor (FIS) maps at 65μm, 90μm, and 140μm. The resulting image profiles are decomposed into the central galaxy component (single term) and the nearby galaxy component (clustering term), as a function of the r-band magnitude, m_r of the central galaxy. We find that the mean far-… ▽ More We perform image stacking analysis of Sloan Digital Sky Survey (SDSS) photometric galaxies over the AKARI Far-Infrared Surveyor (FIS) maps at 65μm, 90μm, and 140μm. The resulting image profiles are decomposed into the central galaxy component (single term) and the nearby galaxy component (clustering term), as a function of the r-band magnitude, m_r of the central galaxy. We find that the mean far-infrared (FIR) flux of a galaxy with magnitude m_r is well fitted with f^s_{90μm}=13*10^{0.306(18-m_r)}[mJy]. The FIR amplitude of the clustering term is consistent with that expected from the angular-correlation function of the SDSS galaxies, but galaxy morphology dependence needs to be taken into account for a more quantitative conclusion. We also fit the spectral energy distribution of stacked galaxies at 65μm, 90μm, and 140μm, and derive a mean dust temperature of ~30K. This is consistent with the typical dust temperature of galaxies that are FIR luminous and individually detected. △ Less

Submitted 9 December, 2015; v1 submitted 9 December, 2015; originally announced December 2015.

Comments: 19 pages, 17 figures, 3tables, accepted for publication in PASJ

Confrontation of Top-Hat Spherical Collapse against Dark Halos from Cosmological N-Body Simulations

Authors: Daichi Suto, Tetsu Kitayama, Ken Osato, Shin Sasaki, Yasushi Suto

Abstract: The top-hat spherical collapse model (TSC) is one of the most fundamental analytical frameworks to describe the non-linear growth of cosmic structure. TSC has motivated, and been widely applied in, various researches even in the current era of precision cosmology. While numerous studies exist to examine its validity against numerical simulations in a statistical fashion, there are few analyses to… ▽ More The top-hat spherical collapse model (TSC) is one of the most fundamental analytical frameworks to describe the non-linear growth of cosmic structure. TSC has motivated, and been widely applied in, various researches even in the current era of precision cosmology. While numerous studies exist to examine its validity against numerical simulations in a statistical fashion, there are few analyses to compare the TSC dynamics in an individual object-wise basis, which is what we attempt in the present paper. We extract 100 halos at z = 0 from a cosmological N-body simulation according to the conventional TSC criterion for the spherical over-density. Then we trace back their spherical counter-parts at earlier epochs. Just prior to the turn-around epoch of the halos, their dynamics is well approximated by TSC, but their turn-around epochs are systematically delayed and the virial radii are larger by ~ 20 percent on average relative to the TSC predictions. We find that this systematic deviation is mainly ascribed to the non-uniformity/inhomogeneity of dark matter density profiles and the non-zero velocity dispersions, both of which are neglected in TSC. In particular, the inside-out-collapse and shell-crossing of dark matter halos play an important role in generating the significant velocity dispersion. The implications of the present result are briefly discussed. △ Less

Submitted 21 November, 2015; originally announced November 2015.

Comments: 17 pages, 40 figures, Accepted for publication in PASJ

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

Authors: Shoya Kamiaka, Kento Masuda, Yuxin Xue, Yasushi Suto, Tsubasa Nishioka, Risa Murakami, Koichiro Inayama, Madoka Saitoh, Michisuke Tanaka, Atsunori Yonehara

Abstract: 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 pre… ▽ More 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. △ Less

Submitted 16 June, 2015; originally announced June 2015.

Comments: 9 pages, 6 Figures, 2 Tables, accepted for publication in PASJ

Constraining spatial extent and temperature of dust around galaxies from far-infrared image stacking analysis

Authors: Toshiya Kashiwagi, Yasushi Suto

Abstract: We propose a novel method to constrain the spatial extent of dust around galaxies through the measurement of dust temperature. Our method combines the dust emission of galaxies from far-infrared (FIR) image stacking analysis and the quasar reddening due to the dust absorption around galaxies. As a specific application of our method, we use the stacked FIR emission profiles of SDSS photometric gala… ▽ More We propose a novel method to constrain the spatial extent of dust around galaxies through the measurement of dust temperature. Our method combines the dust emission of galaxies from far-infrared (FIR) image stacking analysis and the quasar reddening due to the dust absorption around galaxies. As a specific application of our method, we use the stacked FIR emission profiles of SDSS photometric galaxies over the IRAS 100$μ$m map, and the recent measurement of the SDSS galaxy-quasar cross-correlation. If we adopt a single-temperature dust model, the resulting temperature is around 18K, which is consistent with a typical dust temperature for a central part of galaxies. If we assume an additional dust component with much lower temperature, the current data imply the temperature of the galactic dust needs to be higher, 20K to 30K. Since the model of the density and temperature distribution of dust adopted in the current paper is very simple, we cannot draw any strong conclusion at this point. Nevertheless our novel method with the elaborated theoretical model and multi-band measurement of dust will offer an interesting constraint on the statistical nature of galactic dust. △ Less

Submitted 2 June, 2015; originally announced June 2015.

Comments: 8 pages, 4 figures, accepted for publication in MNRAS

Vlasov versus N-body: the Hénon sphere

Authors: S. Colombi, T. Sousbie, S. Peirani, G. Plum, Y. Suto

Abstract: We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized Hénon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual an… ▽ More We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized Hénon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles. △ Less

Submitted 27 April, 2015; originally announced April 2015.

Comments: 19 pages, 11 figures, MNRAS, in press

Modeling the anomaly of surface number densities of galaxies on the Galactic extinction map due to their FIR emission contamination

Authors: Toshiya Kashiwagi, Yasushi Suto, Atsushi Taruya, Issha Kayo, Takahiro Nishimichi, Kazuhiro Yahata

Abstract: The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of th… ▽ More The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of the $r$-band extinction, $A_{r,\rm SFD}$. We confirm that the surface densities of those galaxies positively correlate with $A_{r,\rm SFD}$ for $A_{r,\rm SFD}<0.1$, as first discovered by Yahata et al. (2007) for SDSS DR4 galaxies. Next we construct an analytic model to compute the surface density of galaxies taking account of the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of $100μ{\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (νL)_{100μ{\rm m}}/(νL)_r$. Then we search for the mean and r.m.s values of $y$ that fit the observed anomaly using the analytic model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies (Kashiwagi, Yahata, & Suto 2013). Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. While the corresponding systematic error in the Galactic extinction map is 0.1 to 1mmag, it is directly correlated with galaxy clustering, and thus needs to be carefully examined in precision cosmology. △ Less

Submitted 1 December, 2014; originally announced December 2014.

Comments: 20 pages, 12 figures, accepted for publication in ApJ

Determination of Three-dimensional Spin-orbit Angle with Joint Analysis of Asteroseismology, Transit Lightcurve, and the Rossiter-McLaughlin Effect: Cases of HAT-P-7 and Kepler-25

Authors: O. Benomar, K. Masuda, H. Shibahashi, Y. Suto

Abstract: We develop a detailed methodology of determining three-dimensionally the angle between the stellar spin and the planetary orbit axis vectors, $ψ$, for transiting planetary systems. The determination of $ψ$ requires the independent estimates of the inclination angles of the stellar spin axis and of the planetary orbital axis with respect to the line-of-sight, $i_\star$ and $i_{\rm orb}$, and the pr… ▽ More We develop a detailed methodology of determining three-dimensionally the angle between the stellar spin and the planetary orbit axis vectors, $ψ$, for transiting planetary systems. The determination of $ψ$ requires the independent estimates of the inclination angles of the stellar spin axis and of the planetary orbital axis with respect to the line-of-sight, $i_\star$ and $i_{\rm orb}$, and the projection of the spin--orbit angle onto the plane of the sky, $λ$. These are mainly derived from asteroseismology, transit lightcurve and the Rossiter-McLaughlin effect, respectively. The detailed joint analysis of those three datasets enables an accurate and precise determination of the numerous parameters characterizing the planetary system, in addition to $ψ$. We demonstrate the power of the joint analysis for the two specific systems, HAT-P-7 and Kepler-25. HAT-P-7b is the first exoplanet suspected to be a retrograde (or polar) planet because of the significant misalignment $λ\approx 180^\circ$. Our joint analysis indicates $i_\star \approx {30D}$ and $ψ\approx 120^\circ$, suggesting that the planetary orbit is closer to polar rather than retrograde. Kepler-25 is one of the few multi-transiting planetary systems with measured $λ$, and hosts two short-period transiting planets and one outer non-transiting planet. The projected spin--orbit angle of the larger transiting planet, Kepler-25c, has been measured to be $λ\approx 0^\circ$, implying that the system is well-aligned. With the help of the tight constraint from asteroseismology, however, we obtain $i_\star={65.4}^{+{10.6}}_{-{6.4}}$ and $ψ={26.9}^{+{7.0}}_{-{9.2}}$, and thus find that the system is actually mildly misaligned. △ Less

Submitted 28 July, 2014; originally announced July 2014.

Comments: 20 pages, 13 Figures, 5 Tables

Tidal evolution of the spin-orbit angle in exoplanetary systems

Authors: Yuxin Xue, Yasushi Suto, Atsushi Taruya, Teruyuki Hirano, Yuka Fujii, Kento Masuda

Abstract: 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… ▽ More 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. △ Less

Submitted 23 January, 2014; originally announced January 2014.

Comments: accepted to ApJ

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

Authors: Kento Masuda, Teruyuki Hirano, Atsushi Taruya, Makiko Nagasawa, Yasushi Suto

Abstract: 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 tim… ▽ More 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. △ Less

Submitted 14 November, 2013; v1 submitted 21 October, 2013; originally announced October 2013.

Comments: 18 pages, 17 figures, accepted for publication in ApJ

Report number: YITP-13-90

Journal ref: Astrophys.J.778:185-200,2013

Validity of Hydrostatic Equilibrium in Galaxy Clusters from Cosmological Hydrodynamical Simulations

Authors: Daichi Suto, Hajime Kawahara, Tetsu Kitayama, Shin Sasaki, Yasushi Suto, Renyue Cen

Abstract: We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of the gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the H… ▽ More We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of the gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the HSE assumption (the HSE mass) deviates from the true mass by up to ~ 30 %. This level of departure from HSE is consistent with the previous claims, but our physical interpretation is rather different. We demonstrate that the inertial term in the Euler equations makes a negligible contribution to the total mass, and the overall gravity of the cluster is balanced by the thermal gas pressure gradient and the gas acceleration term. Indeed the deviation from the HSE mass is well explained by the acceleration term at almost all radii. We also clarify the confusion of previous work due to the inappropriate application of the Jeans equations in considering the validity of HSE from the gas dynamics extracted from cosmological hydrodynamical simulations. △ Less

Submitted 20 February, 2013; originally announced February 2013.

Comments: 14 pages, 8 figures, to appear in ApJ