(a) shows the equilibrium height of the load (assuming M _⋆ = 10 × solar mass, other solar masses yield similar results). (b)–(d) show the nondimensional load size Q as a function of insolation level (measured by the substellar point temperature T _s0) and planetary mass M _p , for M _⋆ = 1, 3, and 10 M _Sun, respectively. Red contours denote positive Q ^L and black contours denote negative Q ^L . The black shading in (b)–(d) marks the parameter combinations where the timescale for a specific planet to get tidally locked τ _lock is greater than 10% of the lifespan of the star τ _star. The colorful shading overlaid shows atmospheric escape rate obtained from the analytical model in Kang et al. (2021). Red color corresponds to escape rate beyond 10¹² kg s⁻¹ and blue color corresponds to escape rate below 10² kg s⁻¹. With an escape rate of 10⁷ kg s⁻¹, it will take 0.2 Gyr to completely vaporize a planet of 0.01 M _e . In panels (b)–(d), areas encompassed by the bold 0-contour are likely to undergo TPW.