In response to the need for the Astro2020 Decadal Survey to explicitly engage early career astronomers, the National Academies of Sciences, Engineering, and Medicine hosted the Early Career Astronomer and Astrophysicist Focus Session (ECFS) on October 8-9, 2018 under the auspices of Committee of Astronomy and Astrophysics. The meeting was attended by fifty six pre-tenure faculty, research scientists, postdoctoral scholars, and senior graduate students, as well as eight former decadal survey committee members, who acted as facilitators. The event was designed to educate early career astronomers about the decadal survey process, to solicit their feedback on the role that early career astronomers should play in Astro2020, and to provide a forum for the discussion of a wide range of topics regarding the astrophysics career path. This white paper presents highlights and themes that emerged during two days of discussion. In Section 1, we discuss concerns that emerged regarding the coming decade and the astrophysics career path, as well as specific recommendations from participants regarding how to address them. We have organized these concerns and suggestions into five broad themes. These include (sequentially): (1) adequately training astronomers in the statistical and computational techniques necessary in an era of "big data", (2) responses to the growth of collaborations and telescopes, (3) concerns about the adequacy of graduate and postdoctoral training, (4) the need for improvements in equity and inclusion in astronomy, and (5) smoothing and facilitating transitions between early career stages. Section 2 is focused on ideas regarding the decadal survey itself, including: incorporating early career voices, ensuring diverse input from a variety of stakeholders, and successfully and broadly disseminating the results of the survey.

A suite of discoveries in the last two decades demonstrate that we are now at a point where incorporating magnetic behavior is key for advancing our ability to characterize substellar and planetary systems. The next decade heralds the exciting maturation of the now-burgeoning field of brown dwarf magnetism, and investing now in brown dwarf magnetism will provide a key platform for exploring exoplanetary magnetism and habitability beyond the solar system. We anticipate significant discoveries including: the nature of substellar and planetary magnetic dynamos, the characterization of exo-aurora physics and brown dwarf magnetospheric environments, and the role of satellites in manifestations of substellar magnetic activity. These efforts will require significant new observational capabilities at radio and near infrared wavelengths, dedicated long-term monitoring programs, and committed support for the theoretical modeling efforts underpinning the physical processes of the magnetic phenomena

We have used NSF's Karl G. Jansky Very Large Array (VLA) to observe a sample of five known radio-emitting late L and T dwarfs ranging in age from ~0.2 - 3.4Gyr. We observed each target for seven hours, extending to higher frequencies than previously attempted and establishing proportionally higher limits on maximum surface magnetic field strengths. Detections of circularly polarized pulses at 8 - 12~GHz yield measurements of 3.2 - 4.1 kG localized magnetic fields on four of our targets, including the archetypal cloud variable and likely planetary-mass object T2.5 dwarf SIMP J01365663+0933473. We additionally detect a pulse at 15 - 16.5 GHz for the T6.5 dwarf 2MASS 10475385+2124234, corresponding to a localized 5.6 kG field strength. For the same object, we tentatively detect a 16.5 - 18 GHz pulse, corresponding to a localized 6.2 kG field strength. We measure rotation periods between ~1.47 - 2.28 hr for 2MASS J10430758+2225236, 2MASS J12373919+6526148, and SDSS J04234858-0414035, supporting (i) an emerging consensus that rapid rotation may be important for producing strong dipole fields in convective dynamos and/or (ii) rapid rotation is a key ingredient for driving the current systems powering auroral radio emission. We observe evidence of variable structure in the frequency-dependent time series of our targets on timescales shorter than a rotation period, suggesting a higher degree of variability in the current systems near the surfaces of brown dwarfs. Finally, we find that age, mass, and temperature together cannot account for the strong magnetic fields produced by our targets.

Stellar coronal activity has been shown to persist into the low-mass star regime, down to late M-dwarf spectral types. However, there is now an accumulation of evidence suggesting that at the end of the main sequence there is a transition in the nature of the magnetic activity from chromospheric and coronal to planet-like and auroral, from local impulsive heating via flares and MHD wave dissipation to energy dissipation from strong large-scale magnetospheric current systems. We examine this transition and the prevalence of auroral activity in brown dwarfs through a compilation of multi-wavelength surveys of magnetic activity, including radio, X-ray, and optical. We compile the results of those surveys and place their conclusions in the context of auroral emission as the consequence of large-scale magnetospheric current systems that accelerate energetic electron beams and drive the particles to impact the cool atmospheric gas. We explore the different manifestation of auroral phenomena in brown dwarf atmospheres, like H$\alpha$, and define their distinguishing characteristics. We conclude that large amplitude photometric variability in the near infrared is most likely a consequence of clouds in brown dwarf atmospheres, but that auroral activity may be responsible for long-lived stable surface features. We report a connection between auroral H$\alpha$ emission and quiescent radio emission in ECMI pulsing brown dwarfs, suggesting a potential underlying physical connection between the quiescent and auroral emissions. We also discuss the electrodynamic engines powering brown dwarf aurorae and the possible role of satellites around these systems to both power the aurorae and seed the magnetosphere with plasma.

Recently, studies of brown dwarfs have demonstrated that they possess strong magnetic fields and have the potential to produce radio and optical auroral emissions powered by magnetospheric currents. This emission provides the only window on magnetic fields in the coolest brown dwarfs and identifying additional benchmark objects is key to constraining dynamo theory in this regime. To this end, we conducted a new red optical (6300 - 9700 Angstrom) survey with the Keck telescopes looking for H-alpha emission from a sample of late L dwarfs and T dwarfs. Our survey gathered optical spectra for 29 targets, 18 of which did not have previous optical spectra in the literature, greatly expanding the number of moderate resolution (R~2000) spectra available at these spectral types. Combining our sample with previous surveys, we confirm an H-alpha detection rate of 9.2 (+3.5/-2.1) % for L and T dwarfs in the optical spectral range of L4 - T8. This detection rate is consistent with the recently measured detection rate for auroral radio emission from Kao et al. (2016), suggesting that geometrical selection effects due to the beaming of the radio emission are small or absent. We also provide the first detection of H-alpha emission from 2MASS 0036+1821, previously notable as the only electron cyclotron maser radio source without a confirmed detection of H-alpha emission. Finally, we also establish optical standards for spectral types T3 and T4, filling in the previous gap between T2 and T5.

We have observed 6 late-L and T dwarfs with the Karl G. Jansky Very Large Array (VLA) to investigate the presence of highly circularly polarized radio emission, associated with large-scale auroral currents. Previous surveys encompassing ~60 L6 or later targets in this spectral range have yielded only one detection. Our sample includes the previously detected T6.5 dwarf 2MASS 10475385+2124234 as well as 5 new targets selected for the presence of H-alpha emission or optical/infrared photometric variability, which are possible manifestations of auroral activity. We detect 2MASS 10475385+2124234, as well as 4 of the 5 targets in our biased sample, including the strong IR variable SIMP J01365662+0933473 and bright H-alpha emitter 2MASS 12373919+6526148, reinforcing the possibility that activity at these disparate wavelengths is related. The radio emission frequency corresponds to a precise determination of the lower-bound magnetic field strength near the surface of each dwarf and this new sample provides robust constraints on dynamo theory in the low mass brown dwarf regime. Magnetic fields >2.5 kG are confirmed for 5/6 targets. Our results provide tentative evidence that the dynamo operating in this mass regime may be inconsistent with predicted values from a recently proposed model. Further observations at higher radio frequencies are essential for verifying this assertion.

Aurorae are detected from all the magnetized planets in our Solar System, including Earth. They are powered by magnetospheric current systems that lead to the precipitation of energetic electrons into the high-latitude regions of the upper atmosphere. In the case of the gas-giant planets, these aurorae include highly polarized radio emission at kilohertz and megahertz frequencies produced by the precipitating electrons, as well as continuum and line emission in the infrared, optical, ultraviolet and X-ray parts of the spectrum, associated with the collisional excitation and heating of the hydrogen-dominated atmosphere. Here we report simultaneous radio and optical spectroscopic observations of an object at the end of the stellar main sequence, located right at the boundary between stars and brown dwarfs, from which we have detected radio and optical auroral emissions both powered by magnetospheric currents. Whereas the magnetic activity of stars like our Sun is powered by processes that occur in their lower atmospheres, these aurorae are powered by processes originating much further out in the magnetosphere of the dwarf star that couple energy into the lower atmosphere. The dissipated power is at least four orders of magnitude larger than what is produced in the Jovian magnetosphere, revealing aurorae to be a potentially ubiquitous signature of large-scale magnetospheres that can scale to luminosities far greater than those observed in our Solar System. These magnetospheric current systems may also play a part in powering some of the weather phenomena reported on brown dwarfs.

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1-500 M_Jup, with orbital semi-major axes typically between 1-75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51 +/- 10% for companions with masses between 1-13 M_Jup and orbital semi-major axes between 1-20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

We have performed an analysis of over 34,000 MgII doublets at 0.36 < z < 2.29 in Sloan Digital Sky Survey (SDSS) Data-Release 7 quasar spectra; the catalog, advanced data products, and tools for analysis are publicly available. The catalog was divided into 14 small redshift bins with roughly 2500 doublets in each, and from Monte-Carlo simulations, we estimate 50% completeness at rest equivalent width W_r ~ 0.8 Angstrom. The equivalent-width frequency distribution is described well by an exponential model at all redshifts, and the distribution becomes flatter with increasing redshift, i.e., there are more strong systems relative to weak ones. Direct comparison with previous SDSS MgII surveys reveal that we recover at least 70% of the doublets in these other catalogs, in addition to detecting thousands of new systems. We discuss how these surveys come by their different results, which qualitatively agree but, due to the very small uncertainties, differ by a statistically significant amount. The estimated physical cross-section of MgII-absorbing galaxy halos increased three-fold, approximately, from z = 0.4 --> 2.3, while the W_r >= 1 Angstrom absorber line density grew, dN_MgII/dX, by roughly 45%. Finally, we explore the different evolution of various absorber populations - damped Lyman-alpha absorbers, Lyman-limit systems, strong CIV absorbers, and strong and weaker MgII systems - across cosmic time (0 < z < 6).

We present the first secondary eclipse and phase curve observations for the highly eccentric hot Jupiter HAT-P-2b in the 3.6, 4.5, 5.8, and 8.0 \mu m bands of the Spitzer Space Telescope. The 3.6 and 4.5 \mu m data sets span an entire orbital period of HAT-P-2b, making them the longest continuous phase curve observations obtained to date and the first full-orbit observations of a planet with an eccentricity exceeding 0.2. We present an improved non-parametric method for removing the intrapixel sensitivity variations in Spitzer data at 3.6 and 4.5 \mu m that robustly maps position-dependent flux variations. We find that the peak in planetary flux occurs at 4.39+/-0.28, 5.84+/-0.39, and 4.68+/-0.37 hours after periapse passage with corresponding maxima in the planet/star flux ratio of 0.1138%+/-0.0089%, 0.1162%+/-0.0080%, and 0.1888%+/-0.0072% in the 3.6, 4.5, and 8.0 \mu m bands respectively. We compare our measured secondary eclipse depths to the predictions from a one-dimensional radiative transfer model, which suggests the possible presence of a transient day side inversion in HAT-P-2b's atmosphere near periapse. We also derive improved estimates for the system parameters, including its mass, radius, and orbital ephemeris. Our simultaneous fit to the transit, secondary eclipse, and radial velocity data allows us to determine the eccentricity and argument of periapse of HAT-P-2b's orbit with a greater precision than has been achieved for any other eccentric extrasolar planet. We also find evidence for a long-term linear trend in the radial velocity data. This trend suggests the presence of another substellar companion in the HAT-P-2 system, which could have caused HAT-P-2b to migrate inward to its present-day orbit via the Kozai mechanism.

We have vastly increased the CIV statistics at intermediate redshift by surveying the thousands of quasars in the Sloan Digital Sky Survey Data-Release 7. We visually verified over 16,000 CIV systems with 1.46 < z < 4.55---a sample size that renders Poisson error negligible. Detailed Monte Carlo simulations show we are approximately 50% complete down to rest equivalent widths W_r ~ 0.6 \AA. We analyzed the sample as a whole and in ten small redshift bins with approximately 1500 doublets each. The equivalent width frequency distributions f(W_r) were well modeled by an exponential, with little evolution in shape. In contrast with previous studies that modeled the frequency distribution as a single power law, the fitted exponential gives a finite mass density for the CIV ions. The co-moving line density dN_CIV/dX evolved smoothly with redshift, increasing by a factor of 2.37+/-0.09 from z = 4.55 to 1.96, then plateauing at dN_CIV/dX ~ 0.34 for z = 1.96 to 1.46. Comparing our SDSS sample with z < 1 (ultraviolet) and z > 5 (infrared) surveys, we see an approximately 10-fold increase in dN_CIV/dX over z ~ 6 --> 0, for W_r >= 0.6 \AA. This suggests a monotonic and significant increase in the enrichment of gas outside galaxies over the 12 Gyr lifetime of the universe.

In typical astrophysical environments, the abundance of heavy elements ranges from 0.001 to 2 times the solar concentration. Lower abundances have been seen in select stars in the Milky Way's halo and in two quasar absorption systems at redshift z=3. These are widely interpreted as relics from the early universe, when all gas possessed a primordial chemistry. Before now there have been no direct abundance measurements from the first Gyr after the Big Bang, when the earliest stars began synthesizing elements. Here we report observations of hydrogen and heavy element absorption in a quasar spectrum at z=7.04, when the universe was just 772 Myr old (5.6% its present age). We detect a large column of neutral hydrogen but no corresponding heavy elements, limiting the chemical abundance to less than 1/10,000 the solar level if the gas is in a gravitationally bound protogalaxy, or less than 1/1,000 solar if it is diffuse and unbound. If the absorption is truly intergalactic, it would imply that the universe was neither ionized by starlight nor chemically enriched in this neighborhood at z~7. If it is gravitationally bound, the inferred abundance is too low to promote efficient cooling, and the system would be a viable site to form the predicted but as-yet unobserved massive population III stars in the early universe.