-
Probing atmospheric escape through metastable He I triplet lines in 15 exoplanets observed with SPIRou
Authors:
A. Masson,
S. Vinatier,
B. Bézard,
M. López-Puertas,
M. Lampón,
F. Debras,
A. Carmona,
B. Klein,
E. Artigau,
W. Dethier,
S. Pelletier,
T. Hood,
R. Allart,
V. Bourrier,
C. Cadieux,
B. Charnay,
N. B. Cowan,
N. J. Cook,
X. Delfosse,
J. -F. Donati,
P. -G. Gu,
G. Hébrard,
E. Martioli,
C. Moutou,
O. Venot
, et al. (1 additional authors not shown)
Abstract:
For several years, the metastable helium triplet line has been successfully used as a tracer to probe atmospheric escape in transiting exoplanets. This absorption in the near-infrared (1083.3 nm) can be observed from the ground using high-resolution spectroscopy, providing new constraints on the mass-loss rate and the temperature characterizing the upper atmosphere of close-in exoplanets.
The ai…
▽ More
For several years, the metastable helium triplet line has been successfully used as a tracer to probe atmospheric escape in transiting exoplanets. This absorption in the near-infrared (1083.3 nm) can be observed from the ground using high-resolution spectroscopy, providing new constraints on the mass-loss rate and the temperature characterizing the upper atmosphere of close-in exoplanets.
The aim of this work is to search for the He triplet signature in 15 transiting exoplanets -- ranging from super-Earths to ultrahot Jupiters -- observed with SPIRou, a high-resolution (R~70 000) near-infrared spectropolarimeter at the CFHT, in order to bring new constraints or to improve existing ones regarding atmospheric escape through a homogeneous study.
We developed a full data processing and analysis pipeline to correct for the residual telluric and stellar contributions. We then used two different 1D models based on the Parker-wind equations and nonlocal thermodynamic equilibrium (NLTE) radiative transfer to interpret the observational results.
We confirm published He triplet detections for HAT-P-11 b, HD 189733 b, and WASP-69 b. We tentatively detect the signature of escaping He in HD 209458 b, GJ 3470 b, and WASP-76 b. We report new constraints on the mass-loss rate and temperature for our three detections and set upper limits for the tentative and nondetections. We notably report improved constraints on the mass-loss rate and temperature of the escaping gas for TOI-1807 b, and report a nondetection for the debated atmospheric escape in GJ 1214 b. We also conducted the first search for the He signature in GJ 486 b since its discovery and report a nondetection of the He triplet. Finally, we studied the impact of important model assumptions on our retrieved parameters, notably the limitations of 1D models and the influence of the H/He ratio on the derived constraints.
△ Less
Submitted 13 June, 2024;
originally announced June 2024.
-
Semi-analytic calculations for extended mid-Infrared emission associated with FU Ori-type objects
Authors:
Michihiro Takami,
Pin-Gao Gu,
Gilles Otten,
Christian Delacroix,
Sheng-Yuan Liu,
Shiang-Yu Wang,
Jennifer Karr
Abstract:
(Aims) Near-infrared imaging polarimetry at high-angular resolutions has revealed intriguing distribution of circumstellar dust towards FU Ori-type objects (FUors). These dust grains are probably associated with either an accretion disk or an infalling envelope. Follow-up observations in the mid-infrared would lead us to a better understanding of the hierarchy of the mass accretion processes onto…
▽ More
(Aims) Near-infrared imaging polarimetry at high-angular resolutions has revealed intriguing distribution of circumstellar dust towards FU Ori-type objects (FUors). These dust grains are probably associated with either an accretion disk or an infalling envelope. Follow-up observations in the mid-infrared would lead us to a better understanding of the hierarchy of the mass accretion processes onto FUors (i.e., envelope and disk accretion), which hold keys for understanding the mechanism of their accretion outbursts and the growth of low-mass young stellar objects in general. (Method) We have developed a semi-analytic method to estimate the mid-infrared intensity distributions using the observed polarized intensity (PI) distributions in H-band (lambda=1.65 micron). We have derived intensity distributions for two FUors, FU Ori and V1735 Cyg, at three wavelengths (lambda=3.5/4.8/12 micron) for various cases, i.e. with a star or a flat compact self-luminous disk as an illuminating source; an optically thick disk or an optically thin envelope for circumstellar dust grains; and three different dust models. (Results) We have been able to obtain self-consistent results for many cases and regions, in particular when the viewing angle of the disk/envelope is zero (face-on). Our calculations suggest that the mid-infrared extended emission at the above wavelengths is dominated by the single scattering process. The contribution of thermal emission is negligible unless we add an additional heating mechanism such as adiabatic heating in spiral structures and/or fragments. The uncertain nature of the central illuminating source, the distribution of circumstellar dust grains and the optical properties of dust grains yield uncertainties in the intensity levels on orders of magnitude, e.g., 20-800, for the disk/envelope aspect ratios of ~0.2 and lambda=3-13 micron.
△ Less
Submitted 3 October, 2023;
originally announced October 2023.
-
Deuterium Escape on Photoevaporating Sub-Neptunes
Authors:
Pin-Gao Gu,
Howard Chen
Abstract:
We investigate the evolution of the deuterium-to-hydrogen (D/H) mass ratio driven by EUV photoevaporation of hydrogen-rich atmospheres of close-in sub-Neptunes around solar-type stars. For the first time, the diffusion-limited approach in conjunction with energy-limited photoevaporation is considered in evaluating deuterium escape from evolving exoplanet H/He envelopes. We find that the planets wi…
▽ More
We investigate the evolution of the deuterium-to-hydrogen (D/H) mass ratio driven by EUV photoevaporation of hydrogen-rich atmospheres of close-in sub-Neptunes around solar-type stars. For the first time, the diffusion-limited approach in conjunction with energy-limited photoevaporation is considered in evaluating deuterium escape from evolving exoplanet H/He envelopes. We find that the planets with smaller initial gas envelopes and thus smaller sizes can lead to weaker atmospheric escape, which facilitates hydrogen-deuterium fractionation. Specifically, in our grid of simulations with a low envelope mass fraction less than 0.005, a low-mass sub-Neptune (4-$5M_\oplus$) at about 0.25-0.4 au or a high-mass sub-Neptune (10-$15M_\oplus$) at about 0.1-0.25 au can increase the D/H values by greater than 20% over 7.5 Gyr. Akin to the helium-enhanced envelopes of sub-Neptunes due to photoevaporating escape, the planets along the upper boundary of the radius valley are the best targets to detect high D/H ratios. The ratio can rise by a factor of $\lesssim$ 1.65 within 7.5 Gyrs in our grid of evolutionary calculations. The D/H ratio is expected to be higher in thinner envelopes as long as the planets do not become bare rocky cores.
△ Less
Submitted 18 August, 2023; v1 submitted 9 August, 2023;
originally announced August 2023.
-
Homogeneous search for helium in the atmosphere of 11 gas giant exoplanets with SPIRou
Authors:
R. Allart,
P. -B. Lemée-Joliecoeur,
A. Y. Jaziri,
D. Lafrenière,
E. Artigau,
N. Cook,
A. Darveau-Bernier,
L. Dang,
C. Cadieux,
A. Boucher,
V. Bourrier,
E. K. Deibert,
S. Pelletier,
M. Radica,
B. Benneke,
A. Carmona,
R. Cloutier,
N. B. Cowan,
X. Delfosse,
J. -F. Donati,
R. Doyon,
P. Figueira,
T. Forveille,
P. Fouqué,
E. Gaidos
, et al. (9 additional authors not shown)
Abstract:
The metastable helium triplet in the near-infrared (10833Å) is among the most important probes of exoplanet atmospheres. It can trace their extended outer layers and constrain mass-loss. We use the near-infrared high-resolution spectropolarimeter SPIRou on the CFHT to search for the spectrally resolved helium triplet in the atmospheres of eleven exoplanets, ranging from warm mini-Neptunes to hot J…
▽ More
The metastable helium triplet in the near-infrared (10833Å) is among the most important probes of exoplanet atmospheres. It can trace their extended outer layers and constrain mass-loss. We use the near-infrared high-resolution spectropolarimeter SPIRou on the CFHT to search for the spectrally resolved helium triplet in the atmospheres of eleven exoplanets, ranging from warm mini-Neptunes to hot Jupiters and orbiting G, K, and M dwarfs. Observations were obtained as part of the SPIRou Legacy Survey and complementary open-time programs. We apply a homogeneous data reduction to all datasets and set constraints on the presence of metastable helium, despite the presence of systematics in the data. We confirm published detections for HAT-P-11b, HD189733b, and WASP-69b and set upper limits for the other planets. We apply the p-winds open source code to set upper limits on the mass-loss rate for the non-detections and to constrain the thermosphere temperature, mass-loss rate, line-of-sight velocity, and the altitude of the thermosphere for the detections. We confirm that the presence of metastable helium correlates with the stellar mass and the XUV flux received by the planets. We investigated the correlation between the mass-loss rate and the presence of metastable helium, but it remains difficult to draw definitive conclusions. Finally, some of our results are in contradiction with previous results in the literature, therefore we stress the importance of repeatable, homogeneous, and larger-scale analyses of the helium triplet to obtain robust statistics, study temporal variability, and better understand how the helium triplet can be used to explore the evolution of exoplanets.
△ Less
Submitted 10 July, 2023;
originally announced July 2023.
-
Kinematical Constraint on Eccentricity in the Protoplanetary Disk MWC 758 with ALMA
Authors:
I-Hsuan Genevieve Kuo,
Hsi-Wei Yen,
Pin-Gao Gu,
Tze-En Chang
Abstract:
We analyzed the archival data of the $^{13}\mathrm{CO}$ and $\mathrm{C}^{18}\mathrm{O}$ $J=3-2$ emission lines in the protoplanetary disk around MWC 758 obtained with the Atacama Large Millimeter/submillimeter Array to discuss possible planet-disk interaction and non-Keplerian motion in the disk. We performed fitting of a Keplerian disk model to the observational data and measured the velocity dev…
▽ More
We analyzed the archival data of the $^{13}\mathrm{CO}$ and $\mathrm{C}^{18}\mathrm{O}$ $J=3-2$ emission lines in the protoplanetary disk around MWC 758 obtained with the Atacama Large Millimeter/submillimeter Array to discuss possible planet-disk interaction and non-Keplerian motion in the disk. We performed fitting of a Keplerian disk model to the observational data and measured the velocity deviations from the Keplerian rotation. We found significant velocity deviations around the inner cavity in the MWC 758 disk. We examined several possibilities that may cause the velocity deviations, such as pressure gradient, height of the emitting layer, infall motion, inner warp, and eccentricity in the disk. We found that the combination of an eccentric orbital motion with eccentricity of $0.1\pm0.04$ at the radius of the inner cavity and an infalling flow best explains the observed velocity deviations. Our kinematically constrained eccentricity of the gas orbital motion close to the inner cavity is consistent with the eccentricity of the dust ring around the inner cavity measured in the submillimeter continuum emission. Our results hint at strong dust-gas coupling around the inner cavity and presence of a gas giant planet inside the inner cavity in the MWC 758 disk.
△ Less
Submitted 26 September, 2022;
originally announced September 2022.
-
Revealing the drag instability in one-fluid nonideal MHD simulations of a 1D isothermal C-shock
Authors:
Pin-Gao Gu,
Che-Yu Chen,
Emma Shen,
Chien-Chang Yen,
Min-Kai Lin
Abstract:
C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical enviro…
▽ More
C-type shocks are believed to be ubiquitous in turbulent molecular clouds thanks to ambipolar diffusion. We investigate whether the drag instability in 1D isothermal C-shocks, inferred from the local linear theory of Gu & Chen, can appear in non-ideal magnetohydrodynamic simulations. Two C-shock models (with narrow and broad steady-state shock widths) are considered to represent the typical environment of star-forming clouds. The ionization-recombination equilibrium is adopted for the one-fluid approach. In the 1D simulation, the inflow gas is continuously perturbed by a sinusoidal density fluctuation with a constant frequency. The perturbations clearly grow after entering the C-shock region until they start being damped at the transition to the postshock region. We show that the profiles of a predominant Fourier mode extracted locally from the simulated growing perturbation match those of the growing mode derived from the linear analysis. Moreover, the local growth rate and wave frequency derived from the predominant mode generally agree with those from the linear theory. Therefore, we confirm the presence of the drag instability in simulated 1D isothermal C-shocks. We also explore the nonlinear behavior of the instability by imposing larger-amplitude perturbations to the simulation. We find that the drag instability is subject to wave steepening, leading to saturated perturbation growth. Issues concerning local analysis, nonlinear effects, one-fluid approach, and astrophysical applications are discussed.
△ Less
Submitted 19 August, 2022; v1 submitted 9 July, 2022;
originally announced July 2022.
-
TOI-1759 b: a transiting sub-Neptune around a low mass star characterized with SPIRou and TESS
Authors:
Eder Martioli,
Guillaume Hébrard,
Pascal Fouqué,
Étienne Artigau,
Jean-François Donati,
Charles Cadieux,
Stefano Bellotti,
Alain Lecavelier des Etangs,
Réne Doyon,
J. -D. do Nascimento Jr.,
L. Arnold,
A. Carmona,
N. J. Cook,
P. Cortes-Zuleta,
L. de Almeida,
X. Delfosse,
C. P. Folsom,
P. -C. König,
C. Moutou,
M. Ould-Elhkim,
P. Petit,
K. G. Stassun,
A. A. Vidotto,
T. Vandal,
B. Benneke
, et al. (35 additional authors not shown)
Abstract:
We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial…
▽ More
We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial velocities were extracted from the SPIRou spectra using both the CCF and the LBL methods, optimizing the velocity measurements in the near infrared domain. We analyzed the broadband SED of the star and the high-resolution SPIRou spectra to constrain the stellar parameters and thus improve the accuracy of the derived planet parameters. A LSD analysis of the SPIRou Stokes $V$ polarized spectra detects Zeeman signatures in TOI-1759. We model the rotational modulation of the magnetic stellar activity using a GP regression with a quasi-periodic covariance function, and find a rotation period of $35.65^{+0.17}_{-0.15}$ d. We reconstruct the large-scale surface magnetic field of the star using ZDI, which gives a predominantly poloidal field with a mean strength of $18\pm4$ G. Finally, we perform a joint Bayesian MCMC analysis of the TESS photometry and SPIRou RVs to optimally constrain the system parameters. At $0.1176\pm0.0013$ au from the star, the planet receives $6.4$ times the bolometric flux incident on Earth, and its equilibrium temperature is estimated at $433\pm14$ K. TOI-1759 b is a likely gas-dominated sub-Neptune with an expected high rate of photoevaporation. Therefore, it is an interesting target to search for neutral hydrogen escape, which may provide important constraints on the planetary formation mechanisms responsible for the observed sub-Neptune radius desert.
△ Less
Submitted 2 February, 2022;
originally announced February 2022.
-
Size evolution of close-in super-Earths through giant impacts and photoevaporation
Authors:
Yuji Matsumoto,
Eiichiro Kokubo,
Pin-Gao Gu,
Kenji Kurosaki
Abstract:
The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the siz…
▽ More
The Kepler transit survey with follow-up spectroscopic observations has discovered numerous super-Earth sized planets and revealed intriguing features of their sizes, orbital periods, and their relations between adjacent planets. For the first time, we investigate the size evolution of planets via both giant impacts and photoevaporation to compare with these observed features. We calculate the size of a protoplanet, which is the sum of its core and envelope sizes, by analytical models. $N$-body simulations are performed to evolve planet sizes during the giant impact phase with envelope stripping via impact shocks. We consider the initial radial profile of the core mass and the initial envelope mass fractions as parameters. Inner planets can lose their whole envelopes via giant impacts, while outer planets can keep their initial envelopes since they do not experience giant impacts. Photoevaporation is simulated to evolve planet sizes afterward. Our results suggest that the period-radius distribution of the observed planets would be reproduced if we perform simulations in which the initial radial profile of the core mass follows a wide range of power-law distributions and the initial envelope mass fractions are $\sim0.1$. Moreover, our model shows that the adjacent planetary pairs have similar sizes and regular spacings, with slight differences from detailed observational results such as the radius gap.
△ Less
Submitted 28 September, 2021;
originally announced September 2021.
-
The Drag Instability in a 2D Isothermal C-shock
Authors:
Pin-Gao Gu
Abstract:
We extend the linear analysis of the drag instability in a 1D perpendicular isothermal C-shock by Gu & Chen to 2D perpendicular and oblique C-shocks in the typical environment of star-forming clouds. Simplified dispersion relations are derived for the unstable modes. We find that the mode property of the drag instability generally depends on the ratio of the transverse (normal to the shock flow) t…
▽ More
We extend the linear analysis of the drag instability in a 1D perpendicular isothermal C-shock by Gu & Chen to 2D perpendicular and oblique C-shocks in the typical environment of star-forming clouds. Simplified dispersion relations are derived for the unstable modes. We find that the mode property of the drag instability generally depends on the ratio of the transverse (normal to the shock flow) to longitudinal (along the shock flow) wavenumber. For the transversely large-scale mode, the growth rate and wave frequency of the drag instability in a 2D shock resemble those in a 1D shock. For the transversely small-scale mode, the drag instability is characterized by an unstable mode coupled with an acoustic mode primarily along the transverse direction. When the shock is perpendicular or less oblique, there exists a slowly propagating mode, which can potentially grow into a nonlinear regime and contribute to the maximum growth of the instability. In contrast, when the shock is more oblique, this slowly propagating unstable mode disappears, and the maximum growth of the drag instability is likely contributed from the transversely large-scale mode (i.e., almost 1D mode). In all cases that we consider, the magnitude of the density perturbations is significantly larger than that of the velocity and magnetic field perturbations, implying that the density enhancement governs the dynamics in the linear regime of the instability. A few issues in the linear analysis, as well as the possible astrophysical implications, are also briefly discussed.
△ Less
Submitted 18 June, 2021; v1 submitted 21 April, 2021;
originally announced April 2021.
-
Formation of multiple-planet systems in resonant chains around M dwarfs
Authors:
Yu-Chia Lin,
Yuji Matsumoto,
Pin-Gao Gu
Abstract:
Recent observations have revealed the existence of multiple-planet systems composed of Earth-mass planets around late M dwarfs. Most of their orbits are close to commensurabilities, which suggests that planets were commonly trapped in resonant chains in their formation around low-mass stars. We investigate the formation of multiple-planet systems in resonant chains around low-mass stars. A time-ev…
▽ More
Recent observations have revealed the existence of multiple-planet systems composed of Earth-mass planets around late M dwarfs. Most of their orbits are close to commensurabilities, which suggests that planets were commonly trapped in resonant chains in their formation around low-mass stars. We investigate the formation of multiple-planet systems in resonant chains around low-mass stars. A time-evolution model of the multiple-planet formation via pebble accretion in the early phase of the disk evolution is constructed based on the formation model for the TRAPPIST-1 system by Ormel et al. (2017). Our simulations show that knowing the protoplanet appearance timescale is important for determining the number of planets and their trapped resonances: as the protoplanet appearance timescale increases, fewer planets are formed, which are trapped in more widely separated resonances. We find that there is a range of the protoplanet appearance timescale for forming the stable multiple-planet systems in resonant chains. This range depends on the stellar mass and disk size. We suggest that the protoplanet appearance timescale is a key parameter for studying the formation of multiple-planet systems with planets in resonant chains around low-mass stars. The composition of the planets in our model is also discussed.
△ Less
Submitted 7 December, 2020;
originally announced December 2020.
-
Kinematical signs of dust trapping and feedback in a local pressure bump in the protoplanetary disk around HD 142527 revealed with ALMA
Authors:
Hsi-Wei Yen,
Pin-Gao Gu
Abstract:
We analyzed the archival data of the continuum emission at six wavelengths from 3 to 0.4 mm and 13CO and C18O (1-0, 2-1, and 3-2) lines in the protoplanetary disk around HD 142527 obtained with the Atacama Large Millimeter/submillimeter Array. We performed fitting to the spectral energy distributions obtained at the six wavelengths with the gray-body slab models to estimate the distributions of th…
▽ More
We analyzed the archival data of the continuum emission at six wavelengths from 3 to 0.4 mm and 13CO and C18O (1-0, 2-1, and 3-2) lines in the protoplanetary disk around HD 142527 obtained with the Atacama Large Millimeter/submillimeter Array. We performed fitting to the spectral energy distributions obtained at the six wavelengths with the gray-body slab models to estimate the distributions of the dust surface density and spectral index of dust absorption coefficient beta. We also estimated the distribution of the gas column density by fitting the C18O spectra and measured the disk rotation by fitting the Keplerian disk models to the C18O data. We found super- and sub-Keplerian rotation inside and outside the dust ring in the northwest in the HD 142527 disk, suggestive of the presence of a local pressure bump. In comparison with our estimated dust and gas distributions, the location of the pressure bump is coincident with the region showing a three times higher dust density and a three times lower gas-to-dust mass ratio than the mean values in the disk, suggesting dust trapping in the pressure bump. Nevertheless, there is no correlation between our derived beta distribution and the location of the pressure bump. In addition, we found that the width of the dust ring is comparable or larger than the width of the pressure bump, which could suggest that dust feedback is significant in the pressure bump.
△ Less
Submitted 26 October, 2020;
originally announced October 2020.
-
Ejection of close-in super-Earths around low-mass stars in the giant impact stage
Authors:
Yuji Matsumoto,
Pin-Gao Gu,
Eiichiro Kokubo,
Shoichi Oshino,
Masashi Omiya
Abstract:
Earth-sized planets were observed in close-in orbits around M dwarfs. While more and more planets are expected to be uncovered around M dwarfs, theories of their formation and dynamical evolution are still in their infancy. We investigate the giant impact growth of protoplanets, which includes strong scattering around low-mass stars. The aim is to clarify whether strong scattering around low-mass…
▽ More
Earth-sized planets were observed in close-in orbits around M dwarfs. While more and more planets are expected to be uncovered around M dwarfs, theories of their formation and dynamical evolution are still in their infancy. We investigate the giant impact growth of protoplanets, which includes strong scattering around low-mass stars. The aim is to clarify whether strong scattering around low-mass stars affects the orbital and mass distributions of the planets. We perform $N$-body simulation of protoplanets by systematically surveying the parameter space of the stellar mass and surface density of protoplanets. We find that protoplanets are often ejected after twice or three times close-scattering around late M dwarfs. The ejection sets the upper limit of the largest planet mass. Adopting the surface density scaling linearly with the stellar mass, we find that as the stellar mass decreases less massive planets are formed in orbits with higher eccentricities and inclinations. Under this scaling, we also find that a few close-in protoplanets are generally ejected. The ejection of protoplanets plays an important role in the mass distribution of super-Earths around late M dwarfs. The mass relation of observed close-in super-Earths and their central star mass is well reproduced by ejection.
△ Less
Submitted 28 July, 2020;
originally announced July 2020.
-
The Drag Instability in a 1D Isothermal C-Shock
Authors:
Pin-Gao Gu,
Che-Yu Chen
Abstract:
We investigate whether the drag instability, proposed by Gu et al., occurs in a one-dimensional (1D) C-shock. The 1D background model proposed by Chen and Ostriker for a steady isothermal C-shock is adopted, and a 1D isothermal linear analysis is performed. We confirm the postulation of Gu et al. that the drift velocity between the ions and the neutrals is sufficiently high within a C-shock to all…
▽ More
We investigate whether the drag instability, proposed by Gu et al., occurs in a one-dimensional (1D) C-shock. The 1D background model proposed by Chen and Ostriker for a steady isothermal C-shock is adopted, and a 1D isothermal linear analysis is performed. We confirm the postulation of Gu et al. that the drift velocity between the ions and the neutrals is sufficiently high within a C-shock to allow for the drag instability. We also study the underlying physics of the decaying modes in the shock and post-shock regions. The drag instability is an overstability phenomenon associated with an exponentially growing mode of a propagating wave. We find that the growing wave mode can only propagate downstream within the shock and subsequently decay in the post-shock region. The maximum total growth (MTG) for such an unstable wave before it is damped is estimated in typical environments of star-forming clouds, which is approximately 10-30 times larger than the initial perturbation at the modest shock velocities and can be significantly enhanced several hundred times for a stronger C-shock with a larger width.
△ Less
Submitted 10 July, 2020;
originally announced July 2020.
-
Effect of Dust Rotational Disruption by Radiative Torques and Implications for F-corona decrease revealed by the Parker Solar Probe
Authors:
Thiem Hoang,
Alex Lazarian,
Hyeseung Lee,
Kyungsuk Cho,
Pin-Gao Gu,
Chi-Hang Ng
Abstract:
The first-year results from the Parker Solar Probe (PSP) reveal a gradual decrease of F-coronal dust from distances of $D=0.166-0.336$ AU (or the inner elongations of $\sim 9.22- 18.69~R_{\odot}$) to the Sun (Howard et al. 2019). Such a F-corona decrease cannot be explained by the dust sublimation scenario of the popular silicate composition that implies a dust-free-zone of boundary at heliocentri…
▽ More
The first-year results from the Parker Solar Probe (PSP) reveal a gradual decrease of F-coronal dust from distances of $D=0.166-0.336$ AU (or the inner elongations of $\sim 9.22- 18.69~R_{\odot}$) to the Sun (Howard et al. 2019). Such a F-corona decrease cannot be explained by the dust sublimation scenario of the popular silicate composition that implies a dust-free-zone of boundary at heliocentric radius $R\lesssim 4-5R_{\odot}$, but may be explained by appealing to various dust compositions with different sublimation fronts. In this paper, we present an additional explanation for the F-corona decrease using our newly introduced mechanism of dust destruction so-called Radiative Torque Disruption (RATD) mechanism. We demonstrate that RATD rapidly breaks large grains into nanoparticles so that they can be efficiently destroyed by nonthermal sputtering induced by bombardment of energetic protons from slow solar winds, which extends the dust-free-zone established by thermal sublimation to $R_{dfz}\sim 8R_{\odot}$. Beyond this extended dust-free-zone, we find that the dust mass decreases gradually from $R\sim 42R_{\odot}$ toward the Sun due to partial removal of nanodust by nonthermal sputtering. The joint effect of RATD and nonthermal sputtering can successfully reproduce the gradual decrease of the F-corona between $19-9R_{\odot}$ observed by the PSP. Finally, the RATD mechanism can efficiently produce nanoparticles usually observed in the inner solar system.
△ Less
Submitted 6 July, 2021; v1 submitted 13 April, 2020;
originally announced April 2020.
-
Modeling the Thermal Bulge of A Hot Jupiter with the Two-Stream Approximation
Authors:
Pin-Gao Gu,
Da-Kai Peng,
Chien-Chang Yen
Abstract:
We revisit the problem of thermal bulge of asynchronous hot Jupiters, using HD 209458 b as a fiducial study. We improve upon previous works by using a double-gray atmosphere model and interior structure from MESA as the background state, and then solve for the thermal bulge in response to the semidiurnal component of stellar insolation. The atmosphere model is based on the radiative transfer with…
▽ More
We revisit the problem of thermal bulge of asynchronous hot Jupiters, using HD 209458 b as a fiducial study. We improve upon previous works by using a double-gray atmosphere model and interior structure from MESA as the background state, and then solve for the thermal bulge in response to the semidiurnal component of stellar insolation. The atmosphere model is based on the radiative transfer with Eddington's two-stream approximation. Two opacity cases are considered: the first introduces a greenhouse effect and the second exhibits a strong temperature inversion. We find that for the predominant thermal bulges excited by g-modes of lower orders, our results are qualitatively similar to the adiabatic results from Arras and Socrates (2010). It arises because the perturbed heating due to self-absorption of thermal emissions can be significant (i.e., greenhouse effect) against Newtonian damping, thereby leading to almost undamped thermal bulges. We also find that the contribution to the thermal bulge from the evanescent waves in the convective zone is not negligible, implying that the thermal bulge is not merely confined in the atmosphere and radiative envelope. Assuming the torque balance between the thermal and gravitational bulges, we estimate the tidal quality factor of the planet for gravitational tides to match the observed radius. The limitations of our model are also briefly discussed.
△ Less
Submitted 28 November, 2019;
originally announced November 2019.
-
HL Tau disk in HCO+ (3-2) and (1-0) with ALMA: gas density, temperature, gap, and one-arm spiral
Authors:
Hsi-Wei Yen,
Pin-Gao Gu,
Naomi Hirano,
Patrick M. Koch,
Chin-Fei Lee,
Hauyu Baobab Liu,
Shigehisa Takakuwa
Abstract:
We present our observational results of the 1.1 mm continuum and the HCO+ (3-2) line in HL Tau at angular resolutions of 0.1" obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+ (3-2) and (1-0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+ (3-2) emission, and the stellar mass is estimated to be 2.1+/-0.2 Msun w…
▽ More
We present our observational results of the 1.1 mm continuum and the HCO+ (3-2) line in HL Tau at angular resolutions of 0.1" obtained with ALMA and our data analysis of the 2.9 mm and 1.1 mm continuum and the HCO+ (3-2) and (1-0) lines of the HL Tau disk. The Keplerian rotation of the HL Tau disk is well resolved in the HCO+ (3-2) emission, and the stellar mass is estimated to be 2.1+/-0.2 Msun with a disk inclination angle of 47 deg. The radial profiles of the HCO+ column density and excitation temperature are measured with the LTE analysis of the two transitions of the HCO+ emission. An HCO+ gas gap at a radius of 30 au, where the column density drops by a factor of 4-8, is found in the HCO+ column density profile, coincident with the dust gap traced by the continuum emission. No other clear HCO+ gas gap is seen. This HCO+ gas gap can be opened by a planet with mass of 0.5-0.8 Mj, which is comparable to the planet mass adopted in numerical simulations to form the dust gap at the same radius in the HL Tau disk. In addition to the disk component, a one-arm spiral with a length of ~3" (520 au) stretching out from the inner disk is observed in the HCO+ (3-2) emission. The observed velocity structures along the spiral suggest an infalling and rotating gas stream toward the inner disk.
△ Less
Submitted 13 June, 2019;
originally announced June 2019.
-
Signs of outflow feedback from a nearby young stellar object on the protostellar envelope around HL Tau
Authors:
Hsi-Wei Yen,
Shigehisa Takakuwa,
Pin-Gao Gu,
Naomi Hirano,
Chin-Fei Lee,
Hauyu Baobab Liu,
Sheng-Yuan Liu,
Chun-Ju Wu
Abstract:
HL Tau is a Class I-II protostar embedded in an infalling and rotating envelope and possibly associated with a planet forming disk, and it is co-located in a 0.1 pc molecular cloud with two nearby young stellar objects. Our ALMA observations revealed two arc-like structures on a 1000 au scale connected to the disk, and their kinematics could not be explained with any conventional model of infallin…
▽ More
HL Tau is a Class I-II protostar embedded in an infalling and rotating envelope and possibly associated with a planet forming disk, and it is co-located in a 0.1 pc molecular cloud with two nearby young stellar objects. Our ALMA observations revealed two arc-like structures on a 1000 au scale connected to the disk, and their kinematics could not be explained with any conventional model of infalling and rotational motions. In this work, we investigate the nature of these arc-like structures connected to the HL Tau disk. We conducted new observations in the 13CO and C18O (3-2; 2-1) lines with JCMT and IRAM 30m, and obtained the ACA data with the 7-m array. With the single-dish, ACA, and ALMA data, we analyzed the gas motions on both 0.1 pc and 1000 au scales in the HL Tau region. We constructed new kinematical models of an infalling and rotating envelope with the consideration of relative motion between HL Tau and the envelope. By including the relative motion between HL Tau and its protostellar envelope, our kinematical model can explain the observed velocity features in the arc-like structures. The morphologies of the arc-like structures can also be explained with an asymmetric initial density distribution in our model envelope. In addition, our single-dish results support that HL Tau is located at the edge of a large-scale (0.1 pc) expanding shell driven by the wind or outflow from XZ Tau, as suggested in the literature. The estimated expanding velocity of the shell is comparable to the relative velocity between HL Tau and its envelope in our kinematical model. These results hints that the large-scale expanding motion likely impacts the protostellar envelope around HL Tau and affects its gas kinematics. We found that the mass infalling rate from the envelope onto the HL Tau disk can be decreased by a factor of two due to this impact by the large-scale expanding shell.
△ Less
Submitted 10 January, 2019;
originally announced January 2019.
-
TeV Scale Neutrino Mass Generation, Minimal Inelastic Dark Matter, and High Scale Leptogenesis
Authors:
Pei-Hong Gu,
Hong-Jian He
Abstract:
The seesaw and leptogenesis commonly depend on the masses of same particles, and thus are both realized at the same scale. In this work, we demonstrate a new possibility to realize a TeV-scale neutrino seesaw and a natural high-scale leptogenesis. We extend the standard model by two gauge-singlet scalars, a vector-like iso-doublet fermion and one iso-triplet Higgs scalar. Our model respects a soft…
▽ More
The seesaw and leptogenesis commonly depend on the masses of same particles, and thus are both realized at the same scale. In this work, we demonstrate a new possibility to realize a TeV-scale neutrino seesaw and a natural high-scale leptogenesis. We extend the standard model by two gauge-singlet scalars, a vector-like iso-doublet fermion and one iso-triplet Higgs scalar. Our model respects a softly broken lepton number and an exactly conserved $Z_2$ discrete symmetry. It can achieve three things altogether: (i) realizing a testable type-II seesaw at TeV scale with two nonzero neutrino mass-eigenvalues, (ii) providing a minimal inelastic dark matter from the new fermion doublets, and (iii) accommodating a thermal or nonthermal leptogenesis through the singlet scalar decays. We further analyze the current experimental constraints on our model and discuss the implications for the dark matter direct detections and the LHC searches.
△ Less
Submitted 4 January, 2019; v1 submitted 28 August, 2018;
originally announced August 2018.
-
SMA and ALMA Studies of Disk- and Planet Formation around Low-mass Protostars
Authors:
Shigehisa Takakuwa,
Hsi-Wei Yen,
Ti-Lin Chou,
Nagayoshi Ohashi,
Yusuke Aso,
Patrick M. Koch,
Ruben Krasnopolsky,
Paul T. P. Ho,
Hauyu Baobab Liu,
Naomi Hirano,
Pin-Gao Gu,
Chin-Fei Lee,
Evaria Puspitaningrum,
Yuri Aikawa,
Masahiro N. Machida,
Kazuya Saigo,
Masao Saito,
Kengo Tomida,
Kohji Tomisaka
Abstract:
We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the…
▽ More
We report our current SMA and ALMA studies of disk and planet formation around protostars. We have revealed that $r \gtrsim$100 AU scale disks in Keplerian rotation are ubiquitous around Class I sources. These Class I Keplerian disks are often embedded in rotating and infalling protostellar envelopes. The infalling speeds of the protostellar envelopes are typically $\sim$ 3-times smaller than the free-fall velocities, and the rotational profiles follow the $r^{-1}$ profile, that is, rotation with the conserved specific angular momentum. Our latest high-resolution ($\sim$0$\farcs$5) ALMA studies, as well as the other studies in the literature, have unveiled that $r \sim$100-AU scale Keplerian disks are also present in several Class 0 protostars, while in the other Class 0 sources the inferred upper limits of the Keplerian disks are very small ($r \lessim$20 AU). Our recent data analyses of the ALMA long baseline data of the Class I-II source HL Tau have revealed gaps in molecular gas as well as in dust in the surrounding disk, suggesting the presence of sub-Jovian planets in the disk. These results imply that disk and planet formation should be completed in the protostellar stage.
△ Less
Submitted 1 November, 2017;
originally announced November 2017.
-
Planet Formation in AB Aurigae: Imaging of the inner gaseous Spirals observed inside the Dust Cavity
Authors:
Ya-Wen Tang,
Stephane Guilloteau,
Anne Dutrey,
Takayuki Muto,
Bo-Ting Shen,
Pin-Gao Gu,
Shu-ichiro Inutsuka,
Munetake Momose,
Vincent Pietu,
Misato Fukagawa,
Edwige Chapillon,
Paul T. P. Ho,
Emmanuel di Folco,
Stuartt Corder,
Nagayoshi Ohashi,
Jun Hashimoto
Abstract:
We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within t…
▽ More
We report the results of ALMA observations of a protoplanetary disk surrounding the Herbig Ae star AB Aurigae. We obtained high-resolution (0.1"; 14 au) images in $^{12}$CO (J=2-1) emission and in dust continuum at the wavelength of 1.3 mm. The continuum emission is detected at the center and at the ring with a radius of $\sim$ 120 au. The CO emission is dominated by two prominent spirals within the dust ring. These spirals are trailing and appear to be about 4 times brighter than their surrounding medium. Their kinematics is consistent with Keplerian rotation at an inclination of 23 degree. The apparent two-arm-spiral pattern is best explained by tidal disturbances created by an unseen companion located at 60--80 au, with dust confined in the pressure bumps created outside this companion orbit. An additional companion at r of 30 au, coinciding with the peak CO brightness and a large pitch angle of the spiral, would help to explain the overall emptiness of the cavity. Alternative mechanisms to excite the spirals are discussed. The origin of the large pitch angle detected here remain puzzling.
△ Less
Submitted 9 April, 2017;
originally announced April 2017.
-
Gas Gaps in the Protoplanetary Disk around the Young Protostar HL Tau
Authors:
Hsi-Wei Yen,
Hauyu Baobab Liu,
Pin-Gao Gu,
Naomi Hirano,
Chin-Fei Lee,
Evaria Puspitaningrum,
Shigehisa Takakuwa
Abstract:
We have analyzed the HCO+ (1-0) data of the Class I-II protostar, HL Tau, obtained from the Atacama Large Millimeter/Submillimeter Array long baseline campaign. We generated the HCO+ image cube at an angular resolution of ~0.07 (~10 AU), and performed azimuthal averaging on the image cube to enhance the signal-to-noise ratio and measure the radial profile of the HCO+ integrated intensity. Two gaps…
▽ More
We have analyzed the HCO+ (1-0) data of the Class I-II protostar, HL Tau, obtained from the Atacama Large Millimeter/Submillimeter Array long baseline campaign. We generated the HCO+ image cube at an angular resolution of ~0.07 (~10 AU), and performed azimuthal averaging on the image cube to enhance the signal-to-noise ratio and measure the radial profile of the HCO+ integrated intensity. Two gaps at radii of ~28 AU and ~69 AU and a central cavity are identified in the radial intensity profile. The inner HCO+ gap is coincident with the millimeter continuum gap at a radius of 32 AU. The outer HCO+ gap is located at the millimeter continuum bright ring at a radius of 69 AU and overlaps with the two millimeter continuum gaps at radii of 64 AU and 74 AU. On the contrary, the presence of the central cavity is likely due to the high optical depth of the 3 mm continuum emission and not the depletion of the HCO+ gas. We derived the HCO+ column density profile from its intensity profile. From the column density profile, the full-width-half-maximum widths of the inner and outer HCO+ gaps are both estimated to be ~14 AU, and their depths are estimated to be ~2.4 and ~5.0. These results are consistent with the expectation from the gaps opened by forming (sub-)Jovian mass planets, while placing tight constraints on the theoretical models solely incorporating the variation of dust properties and grain sizes.
△ Less
Submitted 4 March, 2016;
originally announced March 2016.
-
The Effect of Vertical Temperature Gradient on the Propagation of Three-dimensional Waves in a Protoplanetary Disk
Authors:
Wing-Kit Lee,
Pin-Gao Gu
Abstract:
Excitation and propagation of waves in a thermally stratified disk with an arbitrary vertical temperature profile are studied. Previous analytical studies of three-dimensional waves had been focused on either isothermal or polytropic vertical disk structures. However, at the location in a protoplanetary disk where the dominant heating source is stellar irradiation, the temperature gradient may bec…
▽ More
Excitation and propagation of waves in a thermally stratified disk with an arbitrary vertical temperature profile are studied. Previous analytical studies of three-dimensional waves had been focused on either isothermal or polytropic vertical disk structures. However, at the location in a protoplanetary disk where the dominant heating source is stellar irradiation, the temperature gradient may become positive in the vertical direction. We extend the analysis to study the effects of the vertical temperature structure on the waves that are excited at the Lindblad resonances. For a hotter disk atmosphere, the $f$-mode contributes less to the torque and remains confined near the midplane as it propagates away from the resonances. On the other hand, the excitation of the $g$-modes is stronger. As they propagate, they channel to the top of disk atmosphere and their group velocities decrease. The differences compared to previous studies may have implications in understanding the wave dynamics in a realistic disk structure.
△ Less
Submitted 27 August, 2015;
originally announced August 2015.
-
Surface Geometry of Protoplanetary Disks Inferred From Near-Infrared Imaging Polarimetry
Authors:
Michihiro Takami,
Yasuhiro Hasegawa,
Takayuki Muto,
Pin-Gao Gu,
Ruobing Dong,
Jennifer L. Karr,
Jun Hashimoto,
Nobuyuki Kusakabe,
Edwige Chapillon,
Ya-Wen Tang,
Youchi Itoh,
Joseph Carson,
Katherine B. Follette,
Satoshi Mayama,
Michael Sitko,
Markus Janson,
Carol A. Grady,
Tomoyuki Kudo,
Eiji Akiyama,
Jungmi Kwon,
Yasuhiro Takahashi,
Takuya Suenaga,
Lyu Abe,
Wolfgang Brandner,
Timothy D. Brandt
, et al. (36 additional authors not shown)
Abstract:
We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geom…
▽ More
We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.
△ Less
Submitted 18 September, 2014; v1 submitted 4 September, 2014;
originally announced September 2014.
-
Formation of Isothermal Disks around Protoplanets. I. Introductory Three-Dimensional Global Simulations for Sub-Neptune-Mass Protoplanets
Authors:
Hsiang-Hsu Wang,
Defu Bu,
Hsien Shang,
Pin-Gao Gu
Abstract:
The regular satellites found around Neptune ($\approx 17~M_{\Earth}$) and Uranus ($\approx 14.5~M_{\Earth}$) suggest that past gaseous circumplanetary disks may have co-existed with solids around sub-Neptune-mass protoplanets ($< 17~M_{\Earth}$). These disks have been shown to be cool, optically thin, quiescent, with low surface density and low viscosity. Numerical studies of the formation are dif…
▽ More
The regular satellites found around Neptune ($\approx 17~M_{\Earth}$) and Uranus ($\approx 14.5~M_{\Earth}$) suggest that past gaseous circumplanetary disks may have co-existed with solids around sub-Neptune-mass protoplanets ($< 17~M_{\Earth}$). These disks have been shown to be cool, optically thin, quiescent, with low surface density and low viscosity. Numerical studies of the formation are difficult and technically challenging. As an introductory attempt, three-dimensional global simulations are performed to explore the formation of circumplanetary disks around sub-Neptune-mass protoplanets embedded within an isothermal protoplanetary disk at the inviscid limit of the fluid in the absence of self-gravity. Under such conditions, a sub-Neptune-mass protoplanet can reasonably have a rotationally supported circumplanetary disk. The size of the circumplanetary disk is found to be roughly one-tenth of the corresponding Hill radius, which is consistent with the orbital radii of irregular satellites found for Uranus. The protoplanetary gas accretes onto the circumplanetary disk vertically from high altitude and returns to the protoplanetary disk again near the midplane. This implies an open system in which the circumplanetary disk constantly exchanges angular momentum and material with its surrounding prenatal protoplanetary gas.
△ Less
Submitted 3 June, 2014;
originally announced June 2014.
-
3D Structures of equatorial waves and the resulting superrotation in the atmosphere of a tidally locked hot Jupiter
Authors:
Shang-Min Tsai,
Ian Dobbs-Dixon,
Pin-Gao Gu
Abstract:
Three-dimensional equatorial trapped waves excited by stellar isolation and the resulting equatorial superrotating jet in a vertical stratified atmosphere of a tidally-locked hot Jupiter are investigated. Taking the hot Jupiter HD 189733b as a fiducial example, we analytically solve a set of linear equations subject to stationary stellar heating with a uniform zonal-mean flow included. We also ext…
▽ More
Three-dimensional equatorial trapped waves excited by stellar isolation and the resulting equatorial superrotating jet in a vertical stratified atmosphere of a tidally-locked hot Jupiter are investigated. Taking the hot Jupiter HD 189733b as a fiducial example, we analytically solve a set of linear equations subject to stationary stellar heating with a uniform zonal-mean flow included. We also extract wave information in the final equilibrium state of the atmosphere from the radiative hydrodynamical simulation for HD 189733b by Dobbs-Dixon & Agol (2013). We find that the analytic wave solutions are able to qualitatively explain the three-dimensional simulation results. Studying the vertical structure of waves allows us to explore new wave features such as the westward tilt of wavefronts related to the Rossby-wave resonance as well as double gyres of dispersive Rossby waves. We also make an attempt to apply our linear wave analysis to explain some numerical features associated with the equatorial jet development seen in the GCM by Showman & Polvani (2011). During the spin-up phase of the equatorial jet, the acceleration of the jet as a result of the divergence of the wave momentum flux can be in principle boosted by the Rossby-wave resonance as the superrotating jet speed is close to the phase speed of the Rossby waves. However, we also find that as the jet speed increases, the Rossby-wave structure shifts eastward, while the Kelvin-wave structure remains approximately stationary, leading to the decline of the acceleration rate. Our analytic model of jet evolution implies that there exists only one stable equilibrium state of the atmosphere, possibly implying that the final state of the atmosphere is independent of initial conditions in the linear regime. Limitations and future improvements are also discussed.
△ Less
Submitted 20 August, 2014; v1 submitted 30 April, 2014;
originally announced May 2014.
-
On the secular behavior of dust particles in an eccentric protoplanetary disk with an embedded massive gas giant planet
Authors:
He-Feng Hsieh,
Pin-Gao Gu
Abstract:
We investigate the dust velocity and spatial distribution in an eccentric protoplanetary disk under the secular gravitational perturbation of an embedded planet of about 5 Jupiter masses. We first employ the FARGO code to obtain the two-dimensional density and velocity profiles of the eccentric gas disk exterior to the gap opened up by the embedded planet in the quasi-steady state. We then apply t…
▽ More
We investigate the dust velocity and spatial distribution in an eccentric protoplanetary disk under the secular gravitational perturbation of an embedded planet of about 5 Jupiter masses. We first employ the FARGO code to obtain the two-dimensional density and velocity profiles of the eccentric gas disk exterior to the gap opened up by the embedded planet in the quasi-steady state. We then apply the secular perturbation theory and incorporate the gas drag to estimate the dust velocity and density on the secular timescale. The dust-to-gas ratio of the unperturbed disk is simply assumed to be 0.01. In our fiducial disk model with the planet at 5 AU, we find that 0.01 cm- to 1 m-sized dust particles are well coupled to the gas. Consequently, the particles behave similarly to the gas and exhibit asymmetric dynamics as a result of eccentric orbits. The dust surface density is enhanced around the apocenter of the disk. However, for the case of a low-density gaseous disk (termed "transition disk" henceforth in this work) harboring the planet at 100 AU, the azimuthal distributions of dust of various sizes can deviate significantly. Overall, the asymmetric structure exhibits a phase correlation between the gas velocity fields and dust density distribution. Therefore, our study potentially provides a reality check as to whether an asymmetric disk gap detected at sub-millimeter and centimeter wavelengths is a signpost of a massive gas giant planet.
△ Less
Submitted 5 October, 2012;
originally announced October 2012.
-
Coupled evolutions of the stellar obliquity, orbital distance, and planet's radius due to the Ohmic dissipation induced in a diamagnetic hot Jupiter around a magnetic T Tauri star
Authors:
Yu-Ling Chang,
Peter H. Bodenheimer,
Pin-Gao Gu
Abstract:
We revisit the calculation of the Ohmic dissipation in a hot Jupiter presented in Laine et al. (2008) by considering more realistic interior structures, stellar obliquity, and the resulting orbital evolution. In this simplified approach, the young hot Jupiter of one Jupiter mass is modelled as a diamagnetic sphere with a finite resistivity, orbiting across tilted stellar magnetic dipole fields in…
▽ More
We revisit the calculation of the Ohmic dissipation in a hot Jupiter presented in Laine et al. (2008) by considering more realistic interior structures, stellar obliquity, and the resulting orbital evolution. In this simplified approach, the young hot Jupiter of one Jupiter mass is modelled as a diamagnetic sphere with a finite resistivity, orbiting across tilted stellar magnetic dipole fields in vacuum. Since the induced Ohmic dissipation occurs mostly near the planet's surface, we find that the dissipation is unable to significantly expand the young hot Jupiter. Nevertheless, the planet inside a small co-rotation orbital radius can undergo orbital decay by the dissipation torque and finally overfill its Roche lobe during the T Tauri star phase. The stellar obliquity can evolve significantly if the magnetic dipole is parallel/anti-parallel to the stellar spin. Our results are validated by the general torque-dissipation relation in the presence of the stellar obliquity. We also run the fiducial model in Laine et al. (2008) and find that the planet's radius is sustained at a nearly constant value by the Ohmic heating, rather than being thermally expanded to the Roche radius as suggested by the authors.
△ Less
Submitted 7 August, 2012; v1 submitted 4 August, 2012;
originally announced August 2012.
-
Natural Inflation and Flavor Mixing from Peccei-Quinn Symmetry Breaking
Authors:
Pei-Hong Gu,
Manfred Lindner
Abstract:
We propose a left-right symmetric model to simultaneously give natural inflation and flavor mixing from a Peccei-Quinn symmetry breaking at the Planck scale. Our model can be embedded in SO(10) grand unification theories.
We propose a left-right symmetric model to simultaneously give natural inflation and flavor mixing from a Peccei-Quinn symmetry breaking at the Planck scale. Our model can be embedded in SO(10) grand unification theories.
△ Less
Submitted 22 November, 2010;
originally announced November 2010.
-
Tidal and Magnetic Interactions between a Hot Jupiter and its Host Star in the Magnetospheric Cavity of a Protoplanetary Disk
Authors:
Shih-Hsin Chang,
Pin-Gao Gu,
Peter Bodenheimer
Abstract:
We present a simplified model to study the orbital evolution of a young hot Jupiter inside the magnetospheric cavity of a proto-planetary disk. The model takes into account the disk locking of stellar spin as well as the tidal and magnetic interactions between the star and the planet. We focus on the orbital evolution starting from the orbit in the 2:1 resonance with the inner edge of the disk,…
▽ More
We present a simplified model to study the orbital evolution of a young hot Jupiter inside the magnetospheric cavity of a proto-planetary disk. The model takes into account the disk locking of stellar spin as well as the tidal and magnetic interactions between the star and the planet. We focus on the orbital evolution starting from the orbit in the 2:1 resonance with the inner edge of the disk, followed by the inward and then outward orbital migration driven by the tidal and magnetic torques as well as the Roche-lobe overflow of the tidally inflated planet. The goal in this paper is to study how the orbital evolution inside the magnetospheric cavity depends on the cavity size, planet mass, and orbital eccentricity. In the present work, we only target the mass range from 0.7 to 2 Jupiter masses. In the case of the large cavity corresponding to the rotational period ~ 7 days, the planet of mass >1 Jupiter mass with moderate initial eccentricities (>~ 0.3) can move to the region < 0.03 AU from its central star in 10^7 years, while the planet of mass <1 Jupiter mass cannot. We estimate the critical eccentricity beyond which the planet of a given mass will overflow its Roche radius and finally lose all of its gas onto the star due to runaway mass loss. In the case of the small cavity corresponding to the rotational period ~ 3 days, all of the simulated planets lose all of their gas even in circular orbits. Our results for the orbital evolution of young hot Jupiters may have the potential to explain the absence of low-mass giant planets inside ~ 0.03 AU from their dwarf stars revealed by transit surveys.
△ Less
Submitted 23 November, 2009;
originally announced November 2009.
-
Thermal Response of A Solar-like Atmosphere to An Electron Beam from A Hot Jupiter: A Numerical Experiment
Authors:
Pin-Gao Gu,
Takeru K. Suzuki
Abstract:
We investigate the thermal response of the atmosphere of a solar-type star to an electron beam injected from a hot Jupiter by performing a 1-dimensional magnetohydrodynamic numerical experiment with non-linear wave dissipation, radiative cooling, and thermal conduction. In our experiment, the stellar atmosphere is non-rotating and is modelled as a 1-D open flux tube expanding super-radially from…
▽ More
We investigate the thermal response of the atmosphere of a solar-type star to an electron beam injected from a hot Jupiter by performing a 1-dimensional magnetohydrodynamic numerical experiment with non-linear wave dissipation, radiative cooling, and thermal conduction. In our experiment, the stellar atmosphere is non-rotating and is modelled as a 1-D open flux tube expanding super-radially from the stellar photosphere to the planet. An electron beam is assumed to be generated from the reconnection site of the planet's magnetosphere. The effects of the electron beam are then implemented in our simulation as dissipation of the beam momentum and energy at the base of the corona where the Coulomb collisions become effective. When the sufficient energy is supplied by the electron beam, a warm region forms in the chromosphere. This warm region greatly enhances the radiative fluxes corresponding to the temperature of the chromosphere and transition region. The warm region can also intermittently contribute to the radiative flux associated with the coronal temperature due to the thermal instability. However, owing to the small area of the heating spot, the total luminosity of the beam-induced chromospheric radiation is several orders of magnitude smaller than the observed Ca II emissions from HD 179949.
△ Less
Submitted 25 September, 2009;
originally announced September 2009.
-
Double Type-II Seesaw, Baryon Asymmetry and Dark Matter for Cosmic e^\pm Excesses
Authors:
Pei-Hong Gu,
Hong-Jian He,
Utpal Sarkar,
Xinmin Zhang
Abstract:
We construct a new realization of type-II seesaw for neutrino masses and baryon asymmetry by extending the standard model with one light and two heavy singlet scalars besides one Higgs triplet. The heavy singlets pick up small vacuum expectation values to give a suppressed trilinear coupling between the triplet and doublet Higgs bosons after the light singlet drives the spontaneous breaking of l…
▽ More
We construct a new realization of type-II seesaw for neutrino masses and baryon asymmetry by extending the standard model with one light and two heavy singlet scalars besides one Higgs triplet. The heavy singlets pick up small vacuum expectation values to give a suppressed trilinear coupling between the triplet and doublet Higgs bosons after the light singlet drives the spontaneous breaking of lepton number. The Higgs triplet can thus remain light and be accessible at the LHC. The lepton number conserving decays of the heavy singlets can generate a lepton asymmetry stored in the Higgs triplet to account for the matter-antimatter asymmetry in the Universe. We further introduce stable gauge bosons from a hidden sector, which obtain masses and annihilate into the Higgs triplet after spontaneous breaking of the associated non-Abelian gauge symmetry. With Breit-Wigner enhancement, the stable gauge bosons can simultaneously explain the relic density of dark matter and the cosmic positron/electron excesses.
△ Less
Submitted 27 August, 2009; v1 submitted 2 June, 2009;
originally announced June 2009.
-
Diurnal Thermal Tides in a Non-synchronized Hot Jupiter
Authors:
Pin-Gao Gu,
Gordon I. Ogilvie
Abstract:
We perform a linear analysis to investigate the dynamical response of a non-synchronized hot Jupiter to stellar irradiation. In this work, we consider the diurnal Fourier harmonic of the stellar irradiation acting at the top of a radiative layer of a hot Jupiter with no clouds and winds. In the absence of the Coriolis force, the diurnal thermal forcing can excite internal waves propagating into…
▽ More
We perform a linear analysis to investigate the dynamical response of a non-synchronized hot Jupiter to stellar irradiation. In this work, we consider the diurnal Fourier harmonic of the stellar irradiation acting at the top of a radiative layer of a hot Jupiter with no clouds and winds. In the absence of the Coriolis force, the diurnal thermal forcing can excite internal waves propagating into the planet's interior when the thermal forcing period is longer than the sound crossing time of the planet's surface. When the Coriolis effect is taken into consideration, the latitude-dependent stellar heating can excite weak internal waves (g modes) and/or strong baroclinic Rossby waves (buoyant r modes) depending on the asynchrony of the planet. When the planet spins faster than its orbital motion (i.e. retrograde thermal forcing), these waves carry negative angular momentum and are damped by radiative loss as they propagate downwards from the upper layer of the radiative zone. As a result, angular momentum is transferred from the lower layer of the radiative zone to the upper layer and generates a vertical shear. We estimate the resulting internal torques for different rotation periods based on the parameters of HD 209458b.
△ Less
Submitted 21 January, 2009;
originally announced January 2009.
-
Neutrino Mass and Baryon Asymmetry from Dirac Seesaw
Authors:
Pei-Hong Gu,
Hong-Jian He
Abstract:
We extend the visible content of the standard model (SM) with a hidden sector composed of three right-handed singlet neutrinos and one singlet Higgs. These extra singlets are charged under a new U(1)_X gauge symmetry while the SM particles are not. Two heavy scalar doublets are introduced to play the role of the messengers between the visible and hidden sectors. The neutrinos naturally acquire t…
▽ More
We extend the visible content of the standard model (SM) with a hidden sector composed of three right-handed singlet neutrinos and one singlet Higgs. These extra singlets are charged under a new U(1)_X gauge symmetry while the SM particles are not. Two heavy scalar doublets are introduced to play the role of the messengers between the visible and hidden sectors. The neutrinos naturally acquire tiny Dirac masses because the ratio of weak scale over the heavy messenger masses is highly suppressed. Furthermore, the heavy messengers simultaneously generate baryon asymmetry of the universe through their out-of-equilibrium CP-violating decays.
△ Less
Submitted 14 January, 2007; v1 submitted 22 October, 2006;
originally announced October 2006.
-
Possible Effects of Dark Energy on the Detection of Dark Matter Particles
Authors:
Peihong Gu,
Xiao-Jun Bi,
Zhi-Hai Lin,
Xinmin Zhang
Abstract:
We study in this paper the possible influence of the dark energy on the detection of the dark matter particles. In models of dark energy described by a dynamical scalar field such as the Quintessence, its interaction with the dark matter will cause the dark matter particles such as the neutralino vary as a function of space and time. Given a specific model of the Quintessence and its interaction…
▽ More
We study in this paper the possible influence of the dark energy on the detection of the dark matter particles. In models of dark energy described by a dynamical scalar field such as the Quintessence, its interaction with the dark matter will cause the dark matter particles such as the neutralino vary as a function of space and time. Given a specific model of the Quintessence and its interaction in this paper we calculate numerically the corrections to the neutralino masses and the induced spectrum of the neutrinos from the annihilation of the neutralinos pairs in the core of the Sun. This study gives rise to a possibility of probing for dark energy in the experiments of detecting the dark matter particles.
△ Less
Submitted 14 March, 2005;
originally announced March 2005.
-
Hot Jupiters and Hot Spots: The Short- and Long-term Chromospheric Activity on Stars with Giant Planets
Authors:
E. Shkolnik,
G. A. H. Walker,
D. A. Bohlender,
P. -G. Gu,
M. Kürster
Abstract:
We monitored the chromospheric activity in the Ca II H & K lines of 13 solar-type stars (including the Sun); 8 of them over three years at the CFHT and 5 in a single run at the VLT. Ten of the 13 targets have close planetary companions. All of the stars observed at the CFHT show long-term (months to years) changes in H & K intensity levels. Four stars display short-term (days) cyclical activity.…
▽ More
We monitored the chromospheric activity in the Ca II H & K lines of 13 solar-type stars (including the Sun); 8 of them over three years at the CFHT and 5 in a single run at the VLT. Ten of the 13 targets have close planetary companions. All of the stars observed at the CFHT show long-term (months to years) changes in H & K intensity levels. Four stars display short-term (days) cyclical activity. For two, HD 73256 and kappa^1 Ceti, the activity is likely associated with an active region rotating with the star, however, the flaring in excess of the rotational modulation may be associated with a hot jupiter. A planetary companion remains a possibility for kappa^1 Ceti. For the other two, HD 179949 and upsilon And, the cyclic variation is synchronized to the hot jupiter's orbit. For both stars this synchronicity with the orbit is clearly seen in two out of three epochs. The effect is only marginal in the third epoch at which the seasonal level of chromospheric activity had changed for both stars. Short-term chromospheric activity appears weakly dependent on the mean K-line reversal intensities for the sample of 13 stars. Also, a suggestive correlation exists between this activity and the M_p sin(i) of the star's hot jupiter. Because of their small separation (<= 0.1 AU), many of the hot jupiters lie within the Alfvén radius of their host stars which allows a direct magnetic interaction with the stellar surface. We discuss the conditions under which a planet's magnetic field might induce activity on the stellar surface and why no such effect was seen for the prime candidate, tau Boo. This work opens up the possibility of characterizing planet-star interactions, with implications for extrasolar planet magnetic fields and the energy contribution to stellar atmospheres.
△ Less
Submitted 23 November, 2004;
originally announced November 2004.
-
The Internal Structural Adjustment due to Tidal Heating of Short-Period Inflated Giant Planets
Authors:
Pin-Gao Gu,
Peter H. Bodenheimer,
Douglas N. C. Lin
Abstract:
Several short-period Jupiter-mass planets have been discovered around nearby solar-type stars. During the circularization of their orbits, the dissipation of tidal disturbance by their host stars heats the interior and inflates the sizes of these planets. Based on a series of internal structure calculations for giant planets, we examine the physical processes which determine their luminosity-rad…
▽ More
Several short-period Jupiter-mass planets have been discovered around nearby solar-type stars. During the circularization of their orbits, the dissipation of tidal disturbance by their host stars heats the interior and inflates the sizes of these planets. Based on a series of internal structure calculations for giant planets, we examine the physical processes which determine their luminosity-radius relation. In order for young or intensely heated gas giant planets to attain quasi-hydrostatic equilibria, with sizes comparable to or larger than two Jupiter radii, their interiors must have sufficiently high temperature and low density such that degeneracy effects are relatively weak compared to those in a mature and compact Jupiter. Consequently, the polytropic index monotonically increases whereas the central temperature increases and then decreases with the planets' size. These effects, along with a temperature-sensitive opacity for the radiative surface layers of giant planets, cause the power index of the luminosity's dependence on radius to decrease with increasing radius. For planets larger than twice Jupiter's radius, this index is sufficiently small that they become unstable to tidal inflation. We make comparisons between cases of uniform heating and cases in which the heating is concentrated in various locations within the giant planet. Based on these results we suggest that accurate measurement of the sizes of close-in young Jupiters can be used to probe their internal structure under the influence of tidal heating.
△ Less
Submitted 2 March, 2004;
originally announced March 2004.
-
Thermal leptogenesis in a model with mass varying neutrinos
Authors:
Xiao-Jun Bi,
Peihong Gu,
Xiulian Wang,
Xinmin Zhang
Abstract:
In this paper we consider the possibility of neutrino mass varying during the evolution of the Universe and study its implications on leptogenesis. Specifically, we take the minimal seesaw model of neutrino masses and introduce a coupling between the right-handed neutrinos and the dark energy scalar field, the Quintessence. In our model, the right-handed neutrino masses change as the Quintessenc…
▽ More
In this paper we consider the possibility of neutrino mass varying during the evolution of the Universe and study its implications on leptogenesis. Specifically, we take the minimal seesaw model of neutrino masses and introduce a coupling between the right-handed neutrinos and the dark energy scalar field, the Quintessence. In our model, the right-handed neutrino masses change as the Quintessence scalar evolves. We then examine in detail the parameter space of this model allowed by the observed baryon number asymmetry. Our results show that it is possible to lower the reheating temperature in this scenario in comparison with the case that the neutrino masses are unchanged, which helps solve the gravitino problem. Furthermore, a degenerate neutrino mass patten with $m_i$ larger than the upper limit given in the minimal leptogenesis scenario is permitted.
△ Less
Submitted 22 April, 2004; v1 submitted 3 November, 2003;
originally announced November 2003.
-
Dark Energy and Neutrino Mass Limits from Baryogenesis
Authors:
Peihong Gu,
Xiulian Wang,
Xinmin Zhang
Abstract:
In this brief report we consider couplings of the dark energy scalar, such as Quintessence to the neutrinos and discuss its implications in studies on the neutrino mass limits from Baryogenesis. During the evolution of the dark energy scalar, the neutrino masses vary, consequently the bounds on the neutrino masses we have here differ from those obtained before.
In this brief report we consider couplings of the dark energy scalar, such as Quintessence to the neutrinos and discuss its implications in studies on the neutrino mass limits from Baryogenesis. During the evolution of the dark energy scalar, the neutrino masses vary, consequently the bounds on the neutrino masses we have here differ from those obtained before.
△ Less
Submitted 2 August, 2003; v1 submitted 10 July, 2003;
originally announced July 2003.
-
The Effect of Tidal Inflation Instability on the Mass and Dynamical Evolution of Extrasolar Planets with Ultra-Short Periods
Authors:
Pin-Gao Gu,
Doug Lin,
Peter Bodenheimer
Abstract:
We investigate the possibility of substantial inflation of short-period Jupiter-mass planets, as a result of their internal tidal dissipation associated with the synchronization and circularization of their orbits. We employ the simplest prescription based on an equilibrium model with a constant lag angle for all components of the tide. We show that for young Jupiter-mass planets, with a period…
▽ More
We investigate the possibility of substantial inflation of short-period Jupiter-mass planets, as a result of their internal tidal dissipation associated with the synchronization and circularization of their orbits. We employ the simplest prescription based on an equilibrium model with a constant lag angle for all components of the tide. We show that for young Jupiter-mass planets, with a period less than 3 days, an initial radius about 2 Jupiter radii, and an orbital eccentricity greater than 0.2, the energy dissipated during the circularization of their orbits is sufficiently intense and protracted to inflate their sizes up to their Roche radii.
△ Less
Submitted 16 March, 2003;
originally announced March 2003.
-
Vorticity Budget of Weak Thermal Convection in Keplerian disks
Authors:
Pin-Gao Gu
Abstract:
By employing the equations of mean-square vorticity (enstrophy) fluctuations in strong shear flows, we demonstrate that unlike energy production of turbulent vorticity in nonrotating shear flows, the turbulent vorticity of weak convection in Keplerian disks cannot gain energy from vortex stretching/tilting by background shear unless the asscoiated Reynolds stresses are negative. This is because…
▽ More
By employing the equations of mean-square vorticity (enstrophy) fluctuations in strong shear flows, we demonstrate that unlike energy production of turbulent vorticity in nonrotating shear flows, the turbulent vorticity of weak convection in Keplerian disks cannot gain energy from vortex stretching/tilting by background shear unless the asscoiated Reynolds stresses are negative. This is because the epicyclic motion is an energy sink of the radial component of mean-square turbulent vorticity in Keplerian disks when Reynolds stresses are positive. Consequently, weak convection cannot be self-sustained in Keplerian flows. This agrees with the results implied from the equations of mean-square velocity fluctuations in strong shear flows. Our analysis also sheds light on the explanation of the simulation result in which positive kinetic helicity is produced by the Balbus-Hawley instability in a vertically stratified Keplerian disk. We also comment on the possibility of outward angular momentum transport by strong convection based on azimuthal pressure perturbations and directions of energy cascade.
△ Less
Submitted 20 April, 2000; v1 submitted 20 January, 2000;
originally announced January 2000.
-
Thermal Equilibrium Curves and Turbulent Mixing in Keplerian Accretion Disks
Authors:
Pin-Gao Gu,
Ethan T. Vishniac,
John K. Cannizzo
Abstract:
We consider vertical heat transport in Keplerian accretion disks, including the effects of radiation, convection, and turbulent mixing driven by the Balbus-Hawley instability, in astronomical systems ranging from dwarf novae (DNe), and soft X-ray transients (SXTs), to active galactic nuclei (AGN). We propose a modified, anisotropic form of mixing-length theory, which includes radiative and turbu…
▽ More
We consider vertical heat transport in Keplerian accretion disks, including the effects of radiation, convection, and turbulent mixing driven by the Balbus-Hawley instability, in astronomical systems ranging from dwarf novae (DNe), and soft X-ray transients (SXTs), to active galactic nuclei (AGN). We propose a modified, anisotropic form of mixing-length theory, which includes radiative and turbulent damping. We also include turbulent heat transport, which acts everywhere within disks, regardless of whether or not they are stably stratified, and can move entropy in either direction. We have generated a series of vertical structure models and thermal equilibrium curves using the scaling law for the viscosity parameter $α$ suggested by the exponential decay of the X-ray luminosity in SXTs. We have also included equilibrium curves for DNe using an $α$ which is constant down to a small magnetic Reynolds number ($\sim 10^4$). Our models indicate that weak convection is usually eliminated by turbulent radial mixing. The substitution of turbulent heat transport for convection is more important on the unstable branches of thermal equilibrium S-curves when $α$ is larger. The low temperature turnover points $Σ_{max}$ on the equilibrium S-curves are significantly reduced by turbulent mixing in DNe and SXT disks. However, in AGN disks the standard mixing-length theory for convection is still a useful approximation when we use the scaling law for $α$, since these disks are very thin at the relevant radii. In accordance with previous work, we find that constant $α$ models give almost vertical S-curves in the $Σ-T$ plane and consequently imply very slow, possibly oscillating, cooling waves.
△ Less
Submitted 16 August, 1999;
originally announced August 1999.