Pa Chia Thao
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Previous Organizations / Institutes
The organizations / institutes from the scholar's other articles in the past ten years
Research Subjects
Research Interests
Pa Chia Thao has a total of 319 co-authors. Here are the top co-authors by number of shared publications:
- Andrew W. Mann
The University of North Carolina at Chapel Hill
32 shared publications - Elisabeth R. Newton
Dartmouth College
23 shared publications - Adam L. Kraus
University of Texas at Austin
22 shared publications - Andrew Vanderburg
Center for Astrophysics, Harvard & Smithsonian, Cambridge, MA 02138, USA
22 shared publications - Benjamin M. Tofflemire
Department of Astronomy, The University of Texas at Austin, Austin, TX 78712, USA
18 shared publications - Mackenna L. Wood
Department of Physics and Astronomy, The University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
17 shared publications
Latest Publications
- Open Access
- 5 January 2026
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 171 (2) , 63
Searching for Transit Timing Variations in Young Transiting Systems
The discovery of young (4 σ ), although differences in methodology make a totally fair comparison challenging. We show that spots have a weak-to-negligible impact on our results, and similarly cannot explain the higher TTV fraction. Longer-term monitoring will be required to validate these TTVs as planetary in nature and confirm the high TTV rate. While the candidate TTV signals detected here are sparsely sampled, our work provides a clear priority list for additional ground-based observations, and for multiplanet TTVs, to measure the masses and eccentricities of these planets.
- 5 January 2026
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 171 (2) , 63
Open AccessSearching for Transit Timing Variations in Young Transiting Systems
The discovery of young (4 σ ), although differences in methodology make a totally fair comparison challenging. We show that spots have a weak-to-negligible impact on our results, and similarly cannot explain the higher TTV fraction. Longer-term monitoring will be required to validate these TTVs as planetary in nature and confirm the high TTV rate. While the candidate TTV signals detected here are sparsely sampled, our work provides a clear priority list for additional ground-based observations, and for multiplanet TTVs, to measure the masses and eccentricities of these planets.
- Open Access
Searching for Transit Timing Variations in young transiting systems
The discovery of young ( 4$σ$), although differences in methodology make a totally fair comparison challenging. We show that spots have a weak-to-negligible impact on our results, and similarly cannot explain the higher TTV fraction. Longer-term monitoring will be required to validate these TTVs as planetary in nature and confirm the high TTV rate. While the candidate TTV signals detected here are sparsely sampled, our work provides a clear priority list for additional ground-based observations, and for multi-planet TTVs, attempt to measure the masses and eccentricities of these planets.
Open AccessSearching for Transit Timing Variations in young transiting systems
The discovery of young ( 4$σ$), although differences in methodology make a totally fair comparison challenging. We show that spots have a weak-to-negligible impact on our results, and similarly cannot explain the higher TTV fraction. Longer-term monitoring will be required to validate these TTVs as planetary in nature and confirm the high TTV rate. While the candidate TTV signals detected here are sparsely sampled, our work provides a clear priority list for additional ground-based observations, and for multi-planet TTVs, attempt to measure the masses and eccentricities of these planets.
- Open Access
- 1 June 2025
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 986 (1)
Statistical Trends in JWST Transiting Exoplanet Atmospheres
Our brains are hardwired for pattern recognition as correlations are useful for predicting and understanding nature. As more exoplanet atmospheres are being characterized with JWST, we are starting to unveil their properties on a population level. Here we present a framework for comparing exoplanet transmission spectroscopy from 3 to 5 μm with four bands: L (2.9–3.7 μm), SO2 (3.95–4.1 μm), CO2 (4.25–4.4 μm), and CO (4.5–4.9 μm). Together, the four bands cover the major carbon-, oxygen-, nitrogen-, and sulfur-bearing molecules including H2O, CH4, NH3, H2S, SO2, CO2, and CO. Among the eight high-precision gas giant exoplanet planet spectra we collect, we find strong correlations between the SO2 – L index and planet mass (r = −0.41 ± 0.09) and temperature (r = −0.64 ± 0.08), indicating SO2 preferably exists (SO2 – L > −0.5) among low-mass (∼ M J) and cooler (∼ 2 – L and CO – L indices. Under equilibrium chemistry and isothermal thermal structure assumptions, we find that the planet sample favors supersolar metallicity and a low C/O ratio (2 – L versus CO2 – L diagram, similar to the color–magnitude diagram for stars and brown dwarfs. All reported trends here will be testable and will be further quantified with existing and future JWST observations within the next few years.
- 1 June 2025
- Published by American Astronomical Society in The Astrophysical Journal
- Vol. 986 (1)
Open AccessStatistical Trends in JWST Transiting Exoplanet Atmospheres
Our brains are hardwired for pattern recognition as correlations are useful for predicting and understanding nature. As more exoplanet atmospheres are being characterized with JWST, we are starting to unveil their properties on a population level. Here we present a framework for comparing exoplanet transmission spectroscopy from 3 to 5 μm with four bands: L (2.9–3.7 μm), SO2 (3.95–4.1 μm), CO2 (4.25–4.4 μm), and CO (4.5–4.9 μm). Together, the four bands cover the major carbon-, oxygen-, nitrogen-, and sulfur-bearing molecules including H2O, CH4, NH3, H2S, SO2, CO2, and CO. Among the eight high-precision gas giant exoplanet planet spectra we collect, we find strong correlations between the SO2 – L index and planet mass (r = −0.41 ± 0.09) and temperature (r = −0.64 ± 0.08), indicating SO2 preferably exists (SO2 – L > −0.5) among low-mass (∼ M J) and cooler (∼ 2 – L and CO – L indices. Under equilibrium chemistry and isothermal thermal structure assumptions, we find that the planet sample favors supersolar metallicity and a low C/O ratio (2 – L versus CO2 – L diagram, similar to the color–magnitude diagram for stars and brown dwarfs. All reported trends here will be testable and will be further quantified with existing and future JWST observations within the next few years.
- Open Access
- 1 December 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (6) , 297
The Featherweight Giant: Unraveling the Atmosphere of a 17 Myr Planet with JWST
The characterization of young planets (μm transmission spectrum of the 17 Myr, Jupiter-size (R ∼10R ⊕) planet, HIP 67522b, observed with JWST NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous g-band transit with the Southern Astrophysical Research Telescope. The spectrum exhibits absorption features 30%–50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522b’s mass is M ⊕ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of 13.8 ± 1.0 M ⊕. This challenges the previous classification of HIP 67522b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of −3, HIP 67522 b is one of the lowest-density planets known. We find strong absorption from H2O and CO2 (≥7σ), a modest detection of CO (3.5σ), and weak detections of H2S and SO2 (≃2σ). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3 × 10 solar due to the amplitude of the SO2 feature. These results point to the formation of HIP 67522b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01–0.03 M ⊕ Myr−1), sufficient for envelope destruction within a gigayear. This highlights the dramatic evolution occurring within the first 100 Myr of its existence.
- 1 December 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (6) , 297
Open AccessThe Featherweight Giant: Unraveling the Atmosphere of a 17 Myr Planet with JWST
The characterization of young planets (μm transmission spectrum of the 17 Myr, Jupiter-size (R ∼10R ⊕) planet, HIP 67522b, observed with JWST NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous g-band transit with the Southern Astrophysical Research Telescope. The spectrum exhibits absorption features 30%–50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522b’s mass is M ⊕ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of 13.8 ± 1.0 M ⊕. This challenges the previous classification of HIP 67522b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of −3, HIP 67522 b is one of the lowest-density planets known. We find strong absorption from H2O and CO2 (≥7σ), a modest detection of CO (3.5σ), and weak detections of H2S and SO2 (≃2σ). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3 × 10 solar due to the amplitude of the SO2 feature. These results point to the formation of HIP 67522b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01–0.03 M ⊕ Myr−1), sufficient for envelope destruction within a gigayear. This highlights the dramatic evolution occurring within the first 100 Myr of its existence.
- Open Access
The Featherweight Giant: Unraveling the Atmosphere of a 17 Myr Planet with JWST
The characterization of young planets (< 300 Myr) is pivotal for understanding planet formation and evolution. We present the 3-5$μ$m transmission spectrum of the 17 Myr, Jupiter-size ($R$ $\sim$10$R_{\oplus}$) planet, HIP 67522 b, observed with JWST/NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous $g$-band transit with SOAR. The spectrum exhibits absorption features 30-50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522 b's mass is $<$20 $M_{\oplus}$ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of $13.8\pm1.0M_{\oplus}$. This challenges the previous classification of HIP 67522 b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of $<$0.10g/cm$^{3}$, HIP 67522 b is one of the lowest density planets known. We find strong absorption from H$_{2}$O and CO$_{2}$ ($\ge7σ$), a modest detection of CO (3.5$σ$), and weak detections of H$_2$S and SO$_2$ ($\simeq2σ$). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3$\times$10 Solar due to the amplitude of the SO$_2$ feature. These results point to the formation of HIP 67522 b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01-0.03 M$_{\oplus}$ Myr$^{-1}$), sufficient for envelope destruction within a Gyr. This highlights the dramatic evolution occurring within the first 100 Myr of its existence.Open AccessThe Featherweight Giant: Unraveling the Atmosphere of a 17 Myr Planet with JWST
The characterization of young planets (< 300 Myr) is pivotal for understanding planet formation and evolution. We present the 3-5$μ$m transmission spectrum of the 17 Myr, Jupiter-size ($R$ $\sim$10$R_{\oplus}$) planet, HIP 67522 b, observed with JWST/NIRSpec/G395H. To check for spot contamination, we obtain a simultaneous $g$-band transit with SOAR. The spectrum exhibits absorption features 30-50% deeper than the overall depth, far larger than expected from an equivalent mature planet, and suggests that HIP 67522 b's mass is $<$20 $M_{\oplus}$ irrespective of cloud cover and stellar contamination. A Bayesian retrieval analysis returns a mass constraint of $13.8\pm1.0M_{\oplus}$. This challenges the previous classification of HIP 67522 b as a hot Jupiter and instead, positions it as a precursor to the more common sub-Neptunes. With a density of $<$0.10g/cm$^{3}$, HIP 67522 b is one of the lowest density planets known. We find strong absorption from H$_{2}$O and CO$_{2}$ ($\ge7σ$), a modest detection of CO (3.5$σ$), and weak detections of H$_2$S and SO$_2$ ($\simeq2σ$). Comparisons with radiative-convective equilibrium models suggest supersolar atmospheric metallicities and solar-to-subsolar C/O ratios, with photochemistry further constraining the inferred atmospheric metallicity to 3$\times$10 Solar due to the amplitude of the SO$_2$ feature. These results point to the formation of HIP 67522 b beyond the water snowline, where its envelope was polluted by icy pebbles and planetesimals. The planet is likely experiencing substantial mass loss (0.01-0.03 M$_{\oplus}$ Myr$^{-1}$), sufficient for envelope destruction within a Gyr. This highlights the dramatic evolution occurring within the first 100 Myr of its existence. - Open Access
- 1 September 2024
- Published by American Astronomical Society in The Astrophysical Journal Letters
- Vol. 973 (1) , L30
TESS Investigation—Demographics of Young Exoplanets (TI-DYE). II. A Second Giant Planet in the 17 Myr System HIP 67522
The youngest (R ⊕) planet on a tight orbit. In their discovery paper, Rizzuto et al. reported a tentative single-transit detection of an additional planet in the system using TESS. Here, we report the discovery of HIP 67522c, a 7.9 R ⊕ planet that matches with that single-transit event. We confirm the signal with ground-based multiwavelength photometry from Sinistro and MuSCAT4. At a period of 14.33 days, planet c is close to a 2:1 mean-motion resonance with b (6.96 days or 2.06:1). The light curve shows distortions during many of the transits, which are consistent with spot-crossing events and/or flares. Fewer stellar activity events are seen in the transits of planet b, suggesting that planet c is crossing a more active latitude. Such distortions, combined with systematics in the TESS light-curve extraction, likely explain why planet c was previously missed.
- 1 September 2024
- Published by American Astronomical Society in The Astrophysical Journal Letters
- Vol. 973 (1) , L30
Open AccessTESS Investigation—Demographics of Young Exoplanets (TI-DYE). II. A Second Giant Planet in the 17 Myr System HIP 67522
The youngest (R ⊕) planet on a tight orbit. In their discovery paper, Rizzuto et al. reported a tentative single-transit detection of an additional planet in the system using TESS. Here, we report the discovery of HIP 67522c, a 7.9 R ⊕ planet that matches with that single-transit event. We confirm the signal with ground-based multiwavelength photometry from Sinistro and MuSCAT4. At a period of 14.33 days, planet c is close to a 2:1 mean-motion resonance with b (6.96 days or 2.06:1). The light curve shows distortions during many of the transits, which are consistent with spot-crossing events and/or flares. Fewer stellar activity events are seen in the transits of planet b, suggesting that planet c is crossing a more active latitude. Such distortions, combined with systematics in the TESS light-curve extraction, likely explain why planet c was previously missed.
- Open Access
- 23 August 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (3) , 129
Revising Properties of Planet–Host Binary Systems. IV. The Radius Distribution of Small Planets in Binary Star Systems Is Dependent on Stellar Separation*
Small planets (R p ≤ 4 R ⊕) are divided into rocky super-Earths and gaseous sub-Neptunes separated by a radius gap, but the mechanisms that produce these distinct planet populations remain unclear. Binary stars are the only main-sequence systems with an observable record of the protoplanetary disk lifetime and mass reservoir, and the demographics of planets in binaries may provide insights into planet formation and evolution. To investigate the radius distribution of planets in binary star systems, we observed 207 binary systems hosting 283 confirmed and candidate transiting planets detected by the Kepler mission, then recharacterized the planets while accounting for the observational biases introduced by the secondary star. We found that the population of planets in close binaries (ρ ≤ 100 au) is significantly different from the planet population in wider binaries (ρ > 300 au) or single stars. In contrast to planets around single stars, planets in close binaries appear to have a unimodal radius distribution with a peak near the expected super-Earth peak of R p ∼ 1.3 R ⊕ and a suppressed population of sub-Neptunes. We conclude that we are observing the direct impact of a reduced disk lifetime, smaller mass reservoir, and possible altered distribution of solids reducing the sub-Neptune formation efficiency. Our results demonstrate the power of binary stars as a laboratory for exploring planet formation and as a controlled experiment of the impact of varied initial conditions on mature planet populations.
- 23 August 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (3) , 129
Open AccessRevising Properties of Planet–Host Binary Systems. IV. The Radius Distribution of Small Planets in Binary Star Systems Is Dependent on Stellar Separation*
Small planets (R p ≤ 4 R ⊕) are divided into rocky super-Earths and gaseous sub-Neptunes separated by a radius gap, but the mechanisms that produce these distinct planet populations remain unclear. Binary stars are the only main-sequence systems with an observable record of the protoplanetary disk lifetime and mass reservoir, and the demographics of planets in binaries may provide insights into planet formation and evolution. To investigate the radius distribution of planets in binary star systems, we observed 207 binary systems hosting 283 confirmed and candidate transiting planets detected by the Kepler mission, then recharacterized the planets while accounting for the observational biases introduced by the secondary star. We found that the population of planets in close binaries (ρ ≤ 100 au) is significantly different from the planet population in wider binaries (ρ > 300 au) or single stars. In contrast to planets around single stars, planets in close binaries appear to have a unimodal radius distribution with a peak near the expected super-Earth peak of R p ∼ 1.3 R ⊕ and a suppressed population of sub-Neptunes. We conclude that we are observing the direct impact of a reduced disk lifetime, smaller mass reservoir, and possible altered distribution of solids reducing the sub-Neptune formation efficiency. Our results demonstrate the power of binary stars as a laboratory for exploring planet formation and as a controlled experiment of the impact of varied initial conditions on mature planet populations.
- Open Access
- 28 June 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (1) , 41
TESS Hunt for Young and Maturing Exoplanets (THYME). X. A Two-planet System in the 210 Myr MELANGE-5 Association
Young (R ⊕ and 2.88 ± 0.10 R ⊕ and orbit their host star every 4.18 and 17.95 days, respectively. With their bright (K = 9.1 mag), small (R * = 0.44 R ⊙), and cool (T eff = 3326 K) host star, these planets represent excellent candidates for atmospheric characterization with JWST.
- 28 June 2024
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 168 (1) , 41
Open AccessTESS Hunt for Young and Maturing Exoplanets (THYME). X. A Two-planet System in the 210 Myr MELANGE-5 Association
Young (R ⊕ and 2.88 ± 0.10 R ⊕ and orbit their host star every 4.18 and 17.95 days, respectively. With their bright (K = 9.1 mag), small (R * = 0.44 R ⊙), and cool (T eff = 3326 K) host star, these planets represent excellent candidates for atmospheric characterization with JWST.
- Open Access
- 15 November 2023
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 166 (6) , 247
A Lithium Depletion Age for the Carina Association
The dispersed remnants of stellar nurseries, stellar associations, provide unparalleled samples of coeval stars critical for studies of stellar and planetary formation and evolution. The Carina Stellar Association is one of the closest stellar associations to Earth, and yet measurements of its age have varied from 13 to 45 Myr. We aim to update the age of Carina using the lithium depletion boundary (LDB) method. We obtain new measurements of the Li 6708 Å absorption feature in likely members using optical spectra from the Goodman High Throughput Spectrograph on SOAR and NRES on LCO. We detect the depletion boundary at M K ≃ 6.8 (M5). This age is consistent within uncertainties across six different models, including those that account for magnetic fields and spots. We also estimate the age through analysis of the group’s overall variability, and by comparing the association members’ color–magnitude diagram to stellar evolutionary models using a Gaussian Mixture Model, recovering ages consistent with the LDB. Combining these age measures we obtain an age for the Carina association of 41 − 5 + 3 Myr. The resulting age agrees with the older end of previous age measurements and is consistent with the lithium depletion age for the neighboring Tucana-Horologium moving group.
- 15 November 2023
- Published by American Astronomical Society in The Astronomical Journal
- Vol. 166 (6) , 247
Open AccessA Lithium Depletion Age for the Carina Association
The dispersed remnants of stellar nurseries, stellar associations, provide unparalleled samples of coeval stars critical for studies of stellar and planetary formation and evolution. The Carina Stellar Association is one of the closest stellar associations to Earth, and yet measurements of its age have varied from 13 to 45 Myr. We aim to update the age of Carina using the lithium depletion boundary (LDB) method. We obtain new measurements of the Li 6708 Å absorption feature in likely members using optical spectra from the Goodman High Throughput Spectrograph on SOAR and NRES on LCO. We detect the depletion boundary at M K ≃ 6.8 (M5). This age is consistent within uncertainties across six different models, including those that account for magnetic fields and spots. We also estimate the age through analysis of the group’s overall variability, and by comparing the association members’ color–magnitude diagram to stellar evolutionary models using a Gaussian Mixture Model, recovering ages consistent with the LDB. Combining these age measures we obtain an age for the Carina association of 41 − 5 + 3 Myr. The resulting age agrees with the older end of previous age measurements and is consistent with the lithium depletion age for the neighboring Tucana-Horologium moving group.
- Open Access
A Lithium Depletion Age for the Carina Association
The dispersed remnants of stellar nurseries, stellar associations provide unparalleled samples of coeval stars critical for studies of stellar and planetary formation and evolution. The Carina Stellar Association is one of the closest stellar associations to Earth, and yet measurements of its age have varied from 13 to 45 Myr. We aim to update the age of Carina using the Lithium Depletion Boundary method. We obtain new measurements of the Li 6708 Angstrom, absorption feature in likely members using optical spectra from the Goodman HTS on SOAR and NRES on LCO. We detect the depletion boundary at M_K ~= 6.8 (M5), which corresponds to an age of 41(+3,-5) Myr. The age is consistent within uncertainties across six different models, including those that account for magnetic fields and spots. We also estimate the age through analysis of the group's overall variability, and by comparing the association members' CMD to stellar evolutionary models using a Gaussian Mixture Model, recovering ages consistent with the LDB. The resulting age agrees with the older end of previous age measurements and is consistent with the lithium depletion age for the neighboring Tucana-Horologium Moving Group.
Open AccessA Lithium Depletion Age for the Carina Association
The dispersed remnants of stellar nurseries, stellar associations provide unparalleled samples of coeval stars critical for studies of stellar and planetary formation and evolution. The Carina Stellar Association is one of the closest stellar associations to Earth, and yet measurements of its age have varied from 13 to 45 Myr. We aim to update the age of Carina using the Lithium Depletion Boundary method. We obtain new measurements of the Li 6708 Angstrom, absorption feature in likely members using optical spectra from the Goodman HTS on SOAR and NRES on LCO. We detect the depletion boundary at M_K ~= 6.8 (M5), which corresponds to an age of 41(+3,-5) Myr. The age is consistent within uncertainties across six different models, including those that account for magnetic fields and spots. We also estimate the age through analysis of the group's overall variability, and by comparing the association members' CMD to stellar evolutionary models using a Gaussian Mixture Model, recovering ages consistent with the LDB. The resulting age agrees with the older end of previous age measurements and is consistent with the lithium depletion age for the neighboring Tucana-Horologium Moving Group.