Source author record

Pin-Gao Gu

Pin-Gao Gu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

13works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

13 published item(s)

preprint2022arXiv

Revealing the drag instability in one-fluid nonideal MHD simulations of a 1D isothermal C-shock

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.

preprint2020arXiv

Formation of multiple-planet systems in resonant chains around M dwarfs

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.

preprint2020arXiv

The Drag Instability in a 1D Isothermal C-Shock

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.

preprint2019arXiv

Modeling the Thermal Bulge of A Hot Jupiter with the Two-Stream Approximation

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.

preprint2016arXiv

Gas Gaps in the Protoplanetary Disk around the Young Protostar HL Tau

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.

preprint2015arXiv

The Effect of Vertical Temperature Gradient on the Propagation of Three-dimensional Waves in a Protoplanetary Disk

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.

preprint2014arXiv

3D Structures of equatorial waves and the resulting superrotation in the atmosphere of a tidally locked hot Jupiter

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.

preprint2014arXiv

Formation of Isothermal Disks around Protoplanets. I. Introductory Three-Dimensional Global Simulations for Sub-Neptune-Mass Protoplanets

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.

preprint2014arXiv

Surface Geometry of Protoplanetary Disks Inferred From Near-Infrared Imaging Polarimetry

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.

preprint2012arXiv

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

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.

preprint2012arXiv

On the secular behavior of dust particles in an eccentric protoplanetary disk with an embedded massive gas giant planet

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.

preprint2009arXiv

Thermal Response of A Solar-like Atmosphere to An Electron Beam from A Hot Jupiter: A Numerical Experiment

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.

preprint2009arXiv

Tidal and Magnetic Interactions between a Hot Jupiter and its Host Star in the Magnetospheric Cavity of a Protoplanetary Disk

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.