The James Webb Space Telescope (JWST) presents the opportunity to transform our understanding of planets and the origins of life by revealing the atmospheric compositions, structures, and dynamics of transiting exoplanets in unprecedented detail. However, the high-precision, time-series observations required for such investigations have unique technical challenges, and prior experience with other facilities indicates that there will be a steep learning curve when JWST becomes operational. In this paper we describe the science objectives and detailed plans of the Transiting Exoplanet Community Early Release Science (ERS) Program, which is a recently approved program for JWST observations early in Cycle 1. The goal of this project, for which the obtained data will have no exclusive access period, is to accelerate the acquisition and diffusion of technical expertise for transiting exoplanet observations with JWST, while also providing a compelling set of representative datasets that will enable immediate scientific breakthroughs. The Transiting Exoplanet Community ERS Program will exercise the time-series modes of all four JWST instruments that have been identified as the consensus highest priorities, observe the full suite of transiting planet characterization geometries (transits, eclipses, and phase curves), and target planets with host stars that span an illustrative range of brightnesses. The observations in this program were defined through an inclusive and transparent process that had participation from JWST instrument experts and international leaders in transiting exoplanet studies. Community engagement in the project will be centered on a two-phase Data Challenge that culminates with the delivery of planetary spectra, time-series instrument performance reports, and open-source data analysis toolkits in time to inform the agenda for Cycle 2 of the JWST mission.

Planet-planet scattering is the leading mechanism to explain the broad eccentricity distribution of observed giant exoplanets. Here we study the orbital stability of primordial giant planet moons in this scenario. We use N-body simulations including realistic oblateness and evolving spin evolution for the giant planets. We find that the vast majority (80-90% across all our simulations) of orbital parameter space for moons is destabilized. There is a strong radial dependence, as moons past 0.1 Hill radii are systematically removed. Closer-in moons on Galilean-moon-like orbits (< 0.04 Hill radii) have a good ( 20-40%) chance of survival. Destabilized moons may undergo a collision with the star or a planet, be ejected from the system, be captured by another planet, be ejected but still orbiting its free-floating host planet, or survive on heliocentric orbits as "planets." The survival rate of moons increases with the host planet mass but is independent of the planet's final (post-scattering) orbits. Based on our simulations we predict the existence of an abundant galactic population of free- floating (former) moons.

We use Hubble Space Telescope imaging to study the structural properties of ten of the most massive ($M \geq 10^{11.25}$ Msun) quiescent galaxies (QGs) in the UKIDSS UDS at $2.5<z<3.0$. The low spatial density of these galaxies required targeted WFC3 $H_{160}$ imaging, as such systems are rare in existing surveys like CANDELS. We fit Sersic models to the 2D light profiles and find that the median half-light radius is $R_e \sim 3$ kpc, a factor of $\sim 3$ smaller than QGs with similar masses at $z \sim 0$. Complementing our sample with similarly massive QGs at lower redshifts, we find that the median size evolves as $R_e \propto H(z)^{-0.85 \pm 0.12}$ (or alternatively, $R_e \propto (1+z)^{-0.90 \pm 0.12}$). This rate of evolution is slower than that for lower mass QGs. When compared to low redshift QGs, the axis ratio distribution for our high redshift massive QG sample is most consistent with those in which spheroids are dominant. These observations point to earlier size growth among massive QGs that also resulted in spheroidal systems. Finally, we measured residual-corrected surface brightness profiles for our sample. These show that the Sersic parameterization is generally representative out to several effective radii and does not miss excess low surface brightness light. The sizes inferred from the light profiles therefore confirm the compactness of these most massive high redshift QGs.

This work shows the dynamical instability that can happen to close-in satellites when planet oblateness is not accounted for in non-coplanar multiplanet systems. Simulations include two secularly interacting Jupiter-mass planets mutually inclined by 10 degrees, with the host planet either oblate or spherical. With a spherical host planet, moons within a critical planetocentric distance experience high inclinations and in some cases high eccentricities, while more distant moons orbit stably with low inclinations and eccentricities, as expected. These counter-intuitive dynamical phenomena disappear with an oblate host planet, in which case the moons' Laplace plane transitions from the host planet's equatorial plane to the host planet's precessing orbital plane as their semi-major axes increase, and all moons are dynamically stable with very mild changes in orbits. Direct perturbation from the perturbing planet has been investigated and ruled out as an explanation for the behavior of the innermost satellites, therefore leaving the central star's perturbation as the cause. Instability occurs while the nodal precession of the satellite and the central star (as seen from the host planet's frame) approaches the 1:1 secular resonance. In non-coplanar systems, around a non-oblate planet, the nodal precession of the moon becomes slow and comparable to that of the planet, giving rise to resonant configurations. The above effect needs to be taken into account in setting up numerical simulations.