Support for early career astronomers who are just beginning to explore astronomy research is imperative to increase retention of diverse practitioners in the field. Since 2010, Astrobites has played an instrumental role in engaging members of the community -- particularly undergraduate and graduate students -- in research. In this white paper, the Astrobites collaboration outlines our multi-faceted online education platform that both eases the transition into astronomy research and promotes inclusive professional development opportunities. We additionally offer recommendations for how the astronomy community can reduce barriers to entry to astronomy research in the coming decade.

A primary goal of many undergraduate- and graduate-level courses in the physical sciences is to prepare students to engage in scientific research, or to prepare students for careers that leverage skillsets similar to those used by research scientists. Even for students who may not intend to pursue a career with these characteristics, exposure to the context of applications in modern research can be a valuable tool for teaching and learning. However, a persistent barrier to student participation in research is familiarity with the technical language, format, and context that academic researchers use to communicate research methods and findings with each other: the literature of the field. Astrobites, an online web resource authored by graduate students, has published brief and accessible summaries of more than 1,300 articles from the astrophysical literature since its founding in 2010. This article presents three methods for introducing students at all levels within the formal higher education setting to approaches and results from modern research. For each method, we provide a sample lesson plan that integrates content and principles from Astrobites, including step-by-step instructions for instructors, suggestions for adapting the lesson to different class levels across the undergraduate and graduate spectrum, sample student handouts, and a grading rubric.

We investigate the spectral properties of the brightest gamma-ray flares of blazars detected by the Fermi Large Area Telescope. We search for the presence of spectral breaks and measure the spectral curvature on typical time scales of a few days. We identify significant spectral breaks in fewer than half of the analyzed flares, but their parameters do not show any discernible regularities, and in particular there is no indication for gamma-ray absorption at any fixed source-frame photon energy. More interestingly, we find that the studied blazars are characterized by significant spectral variability. Gamma-ray flares of short duration are often characterized by strong spectral curvature, with the spectral peak located above 100 MeV. Since these spectral variations are observed despite excellent photon statistics, they must reflect temporal fluctuations in the energy distributions of the emitting particles. We suggest that highly regular gamma-ray spectra of blazars integrated over long time scales emerge from a superposition of many short-lived irregular components with relatively narrow spectra. This would imply that the emitting particles are accelerated in strongly turbulent environments.

Hot relativistic jets, passing through a background medium with a pressure gradient p \propto r^{-\eta} where 2 < \eta <= 8/3, develop a shocked boundary layer containing a significant fraction of the jet power. In previous work, we developed a self-similar description of the boundary layer assuming isentropic flow, but we found that such models respect global energy conservation only for the special case \eta = 8/3. Here we demonstrate that models with \eta < 8/3 can be made self-consistent if we relax the assumption of constant specific entropy. Instead, the entropy must increase with increasing r along the boundary layer, presumably due to multiple shocks driven into the flow as it gradually collimates. The increase in specific entropy slows the acceleration rate of the flow and provides a source of internal energy that could be channeled into radiation. We suggest that this process may be important for determining the radiative characteristics of tidal disruption events and gamma-ray bursts from collapsars.

Because undergraduate participation in research is a longstanding and increasingly important aspect of the career path for future scientists, students can benefit from additional resources to introduce them to the culture and process of research. We suggest the adoption of the web resource Astrobites as a classroom tool to increase the preparation of undergraduate physics and astronomy students for careers in research. We describe the content and development of the website, discuss previous university courses that have made use of Astrobites, and suggest additional strategies for using Astrobites in the classroom.

We study the collimation of relativistic magnetohydrodynamic jets by the pressure of an ambient medium, in the limit where the jet interior loses causal contact with its surroundings. This follows up a hydrodynamic study in a previous paper, adding the effects of a toroidal magnetic field threading the jet. As the ultrarelativistic jet encounters an ambient medium with a pressure profile with a radial scaling of p ~ r^-eta where 2<eta<4, it loses causal contact with its surroundings and forms a boundary layer with a large pressure gradient. By constructing self-similar solutions to the fluid equations within this boundary layer, we examine the structure of this layer as a function of the external pressure profile. We show that the boundary layer always becomes magnetically dominated far from the source, and that in the magnetic limit, physical self-similar solutions are admitted in which the total pressure within the layer decreases linearly with distance from the contact discontinuity inward. These solutions suggest a `hollow cone' behavior of the jet, with the boundary layer thickness prescribed by the value of eta. In contrast to the hydrodynamical case, however, the boundary layer contains an asymptotically vanishing fraction of the jet energy flux.

We study the collimation of relativistic hydrodynamic jets by the pressure of an ambient medium in the limit where the jet interior has lost causal contact with its surroundings. For a jet with an ultrarelativistic equation of state and external pressure that decreases as a power of spherical radius, p \propto r^(-eta), the jet interior will lose causal contact when eta > 2. However, the outer layers of the jet gradually collimate toward the jet axis as long as eta < 4, leading to the formation of a shocked boundary layer. Assuming that pressure-matching across the shock front determines the shape of the shock, we study the resulting structure of the jet in two ways: first by assuming that the pressure remains constant across the shocked boundary layer and looking for solutions to the shock jump equations, and then by constructing self-similar boundary-layer solutions that allow for a pressure gradient across the shocked layer. We demonstrate that the constant-pressure solutions can be characterized by four initial parameters that determine the jet shape and whether the shock closes to the axis. We show that self-similar solutions for the boundary layer can be constructed that exhibit a monotonic decrease in pressure across the boundary layer from the contact discontinuity to the shock front, and that the addition of this pressure gradient in our initial model generally causes the shock front to move outwards, creating a thinner boundary layer and decreasing the tendency of the shock to close. We discuss trends based on the value of the pressure power-law index eta.