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.

Global 3D simulations of solar giant-cell convection have provided significant insight into the processes which yield the Sun's observed differential rotation and cyclic dynamo action. However, as we move to higher resolution simulations a variety of codes have encountered what has been termed the convection conundrum. As these simulations increase in resolution and hence the level of turbulence achieved, they tend to produce weak or even anti-solar differential rotation patterns associated with a weak rotational influence (high Rossby number) due to large convective velocities. One potential culprit for this convection conundrum is the upper boundary condition applied in most simulations which is generally impenetrable. Here we present an alternative stochastic plume boundary condition which imposes small-scale convective plumes designed to mimic near-surface convective downflows, thus allowing convection to carry the majority of the outward solar energy flux up to and through our simulated upper boundary. The use of a plume boundary condition leads to significant changes in the convective driving realized in the simulated domain and thus to the convective energy transport, the dominant scale of the convective enthalpy flux, and the relative strength of the strongest downflows, the downflow network, and the convective upflows. These changes are present even far from the upper boundary layer. Additionally, we demonstrate that in spite of significant changes, giant cell morphology in the convective patterns is still achieved with self-organization of the imposed boundary plumes into downflow lanes, cellular patterns, and even rotationally-aligned banana cells in equatorial regions. This plume boundary presents an alternative pathway for 3D global convection simulations where driving is non-local and may provide a new approach towards addressing the convection conundrum.

The Astropy Project (http://astropy.org) is, in its own words, "a community effort to develop a single core package for Astronomy in Python and foster interoperability between Python astronomy packages." For five years this project has been managed, written, and operated as a grassroots, self-organized, almost entirely volunteer effort while the software is used by the majority of the astronomical community. Despite this, the project has always been and remains to this day effectively unfunded. Further, contributors receive little or no formal recognition for creating and supporting what is now critical software. This paper explores the problem in detail, outlines possible solutions to correct this, and presents a few suggestions on how to address the sustainability of general purpose astronomical software.

We present multidimensional modeling of convection and oscillations in main-sequence stars somewhat more massive than the Sun, using three separate approaches: 1) Using the 3-D planar StellarBox radiation hydrodynamics code to model the envelope convection zone and part of the radiative zone. Our goals are to examine the interaction of stellar pulsations with turbulent convection in the envelope, excitation of acoustic modes, and the role of convective overshooting; 2) Applying the spherical 3-D MHD ASH (Anelastic Spherical Harmonics) code to simulate the core convection and radiative zone. Our goal is to determine whether core convection can excite low-frequency gravity modes, and thereby explain the presence of low frequencies for some hybrid gamma Doradus/delta Scuti variables for which the envelope convection zone is too shallow for the convective blocking mechanism to drive gravity modes; 3) Applying the ROTORC 2-D stellar evolution and dynamics code to calculate evolution with a variety of initial rotation rates and extents of core convective overshooting. The nonradial adiabatic pulsation frequencies of these nonspherical models are calculated using the 2-D pulsation code NRO. We present new insights into pulsations for stars of one to two solar masses gained by multidimensional modeling.

The mixing of materials due to the Richtmyer-Meshkov instability and the ensuing turbulent behavior is of intense interest in a variety of physical systems including inertial confinement fusion, combustion, and the final stages of stellar evolution. Extensive numerical and laboratory studies of shock-driven mixing have demonstrated the rich behavior associated with the onset of turbulence due to the shocks. Here we report on progress in understanding shock-driven mixing at interfaces between fluids of differing densities through 3D numerical simulations using the RAGE code in the implicit large eddy simulation context. We consider a shock tube configuration with a band of high density gas (SF$_6$) embedded in low density gas (air). Shocks with a Mach number of 1.26 are passed through SF$_6$ bands, resulting in transition to turbulence driven by the Richtmyer-Meshkov instability. The system is followed as a rarefaction wave and a reflected secondary shock from the back wall pass through the SF$_6$ band. We apply a variety of initial perturbations to the interfaces between the two fluids in which the physical standard deviation, wave number range, and the spectral slope of the perturbations are held constant, but the number of modes initially present is varied. By thus decreasing the density of initial spectral modes of the interface, we find that we can achieve as much as 25\% less total mixing at late times. This has potential direct implications for the treatment of initial conditions applied to material interfaces in both 3D and reduced dimensionality simulation models.

Our Sun exhibits strong convective dynamo action which results in magnetic flux bundles emerging through the stellar surface as magnetic spots. Global-scale dynamo action is believed to generate large-scale magnetic structures in the deep solar interior through the interplay of convection, rotation, and shear. Portions of these large-scale magnetic structures are then believed to rise through the convective layer, forming magnetic loops which then pierce the photosphere as sunspot pairs. Previous global simulations of 3D MHD convection in rotating spherical shells have demonstrated mechanisms whereby large-scale magnetic wreaths can be generated in the bulk of the convection zone. Our recent simulations have achieved sufficiently high levels of turbulence to permit portions of these wreaths to become magnetically buoyant and rise through the simulated convective layer through a combination of magnetic buoyancy and advection by convective giant cells. These buoyant magnetic loops are created in the bulk of the convective layer as strong Lorentz force feedback in the cores of the magnetic wreaths dampen small-scale convective motions, permitting the amplification of local magnetic energies to over 100 times the local kinetic energy. While the magnetic wreaths are largely generated the shearing of axisymmetric poloidal magnetic fields by axisymmetric rotational shear (the $\Omega$-effect), the loops are amplified to their peak field strengths before beginning to rise by non-axisymmetric processes. This further extends and enhances a new paradigm for the generation of emergent magnetic flux bundles, which we term turbulence-enabled magnetic buoyancy.

Our global 3D simulations of convection and dynamo action in a Sun-like star reveal that persistent wreaths of strong magnetism can be built within the bulk of the convention zone. Here we examine the characteristics of buoyant magnetic structures that are self-consistently created by dynamo action and turbulent convective motions in a simulation with solar stratification but rotating at three times the current solar rate. These buoyant loops originate within sections of the magnetic wreaths in which turbulent flows amplify the fields to much larger values than is possible through laminar processes. These amplified portions can rise through the convective layer by a combination of magnetic buoyancy and advection by convective giant cells, forming buoyant loops. We measure statistical trends in the polarity, twist, and tilt of these loops. Loops are shown to preferentially arise in longitudinal patches somewhat reminiscent of active longitudes in the Sun, although broader in extent. We show that the strength of the axisymmetric toroidal field is not a good predictor of the production rate for buoyant loops or the amount of magnetic flux in the loops that are produced.

Solar-type stars exhibit a rich variety of magnetic activity. Seeking to explore the convective origins of this activity, we have carried out a series of global 3D magnetohydrodynamic (MHD) simulations with the anelastic spherical harmonic (ASH) code. Here we report on the dynamo mechanisms achieved as the effects of artificial diffusion are systematically decreased. The simulations are carried out at a nominal rotation rate of three times the solar value (3$\Omega_\odot$), but similar dynamics may also apply to the Sun. Our previous simulations demonstrated that convective dynamos can build persistent toroidal flux structures (magnetic wreaths) in the midst of a turbulent convection zone and that high rotation rates promote the cyclic reversal of these wreaths. Here we demonstrate that magnetic cycles can also be achieved by reducing the diffusion, thus increasing the Reynolds and magnetic Reynolds numbers. In these more turbulent models, diffusive processes no longer play a significant role in the key dynamical balances that establish and maintain the differential rotation and magnetic wreaths. Magnetic reversals are attributed to an imbalance in the poloidal magnetic induction by convective motions that is stabilized at higher diffusion levels. Additionally, the enhanced levels of turbulence lead to greater intermittency in the toroidal magnetic wreaths, promoting the generation of buoyant magnetic loops that rise from the deep interior to the upper regions of our simulated domain. The implications of such turbulence-induced magnetic buoyancy for solar and stellar flux emergence are also discussed.

The current dynamo paradigm for the Sun and sun-like stars places the generation site for strong toroidal magnetic structures deep in the solar interior. Sunspots and star-spots on sun-like stars are believed to arise when sections of these magnetic structures become buoyantly unstable and rise from the deep interior to the photosphere. Here we present the first 3-D global magnetohydrodynamic (MHD) simulation in which turbulent convection, stratification, and rotation combine to yield a dynamo that self-consistently generates buoyant magnetic loops. We simulate stellar convection and dynamo action in a spherical shell with solar stratification, but rotating three times faster than the current solar rate. Strong wreaths of toroidal magnetic field are realized by dynamo action in the convection zone. By turning to a dynamic Smagorinsky model for subgrid-scale turbulence, we here attain considerably reduced diffusion in our simulation. This permits the regions of strongest magnetic field in these wreaths to rise toward the top of the convection zone via a combination of magnetic buoyancy instabilities and advection by convective giant cells. Such a global simulation yielding buoyant loops represents a significant step forward in combining numerical models of dynamo action and flux emergence.

Observations of sun-like stars rotating faster than our current sun tend to exhibit increased magnetic activity as well as magnetic cycles spanning multiple years. Using global simulations in spherical shells to study the coupling of large-scale convection, rotation, and magnetism in a younger sun, we have probed effects of rotation on stellar dynamos and the nature of magnetic cycles. Major 3-D MHD simulations carried out at three times the current solar rotation rate reveal hydromagnetic dynamo action that yields wreaths of strong toroidal magnetic field at low latitudes, often with opposite polarity in the two hemispheres. Our recent simulations have explored behavior in systems with considerably lower diffusivities, achieved with sub-grid scale models including a dynamic Smagorinsky treatment of unresolved turbulence. The lower diffusion promotes the generation of magnetic wreaths that undergo prominent temporal variations in field strength, exhibiting global magnetic cycles that involve polarity reversals. In our least diffusive simulation, we find that magnetic buoyancy coupled with advection by convective giant cells can lead to the rise of coherent loops of magnetic field toward the top of the simulated domain.

Stellar dynamos are driven by complex couplings between rotation and turbulent convection, which drive global-scale flows and build and rebuild stellar magnetic fields. When stars like our sun are young, they rotate much more rapidly than the current solar rate. Observations generally indicate that more rapid rotation is correlated with stronger magnetic activity and perhaps more effective dynamo action. Here we examine the effects of more rapid rotation on dynamo action in a star like our sun. We find that vigorous dynamo action is realized, with magnetic field generated throughout the bulk of the convection zone. These simulations do not possess a penetrative tachocline of shear where global-scale fields are thought to be organized in our sun, but despite this we find strikingly ordered fields, much like sea-snakes of toroidal field, which are organized on global scales. We believe this to be a novel finding.