We present numerical simulations, using two complementary setups, of rotating Boussinesq thermal convection in a three-dimensional Cartesian geometry with misaligned gravity and rotation vectors. This model represents a small region at a non-polar latitude in the convection zone of a star or planet. We investigate the effects of rotation on the bulk properties of convection at different latitudes, focusing on determining the relation between the heat flux and temperature gradient. We show that our results may be interpreted using rotating mixing length theory (RMLT). The simplest version of RMLT (due to Stevenson) considers the single mode that transports the most heat. This works reasonably well in explaining our results, but there is a systematic departure from these predictions (up to approximately $30\%$ in the temperature gradient) at mid-latitudes. We develop a more detailed treatment of RMLT that includes the transport afforded by multiple modes, and we show that this accounts for most of the systematic differences. We also show that convectively-generated zonal flows and meridional circulations are produced in our simulations, and that their properties depend strongly on the dimensions of the box. These flows also affect the heat transport, contributing to departures from RMLT at some latitudes. However, we find the theoretical predictions of the multi-mode theory for the mid-layer temperature gradient, the root-mean-square (RMS) vertical velocity, the RMS temperature fluctuation, and the spatial spectrum of the heat transport at different latitudes, are all in reasonably good agreement with our numerical results when zonal flows are small.

It is now recognized that energetic stellar photon and particle radiation evaporates and erodes planetary atmospheres and controls upper atmospheric chemistry. Key exoplanet host stars will be too faint at X-ray wavelengths for accurate characterization using existing generation and future slated X-ray telescopes. Observation of stellar coronal mass ejections and winds are also beyond current instrumentation. In line with theCommittee on an Exoplanet Science Strategy recognition that holistic observational approaches are needed, we point out here that a full understanding of exoplanet atmospheres, their evolution and determination of habitability requires a powerful high-resolution X-ray imaging and spectroscopic observatory. This is the only capability that can: (1) characterize by proxy the crucial, difficult to observe, EUV stellar flux, its history and its variations for planet hosting stars; (2) observe the stellar wind; (3) detect the subtle Doppler signatures of coronal mass ejections.

The lowest-mass stars, brown dwarfs and giant exoplanets span a minimum in the mass-radius relationship that probes the fundamental physics of extreme states of matter, magnetism, and fusion. This White Paper outlines scientific opportunities and the necessary resources for modeling and measuring the mass-radius relationship in this regime.

Some low-mass stars appear to have larger radii than predicted by standard 1D structure models; prior work has suggested that inefficient convective heat transport, due to rotation and/or magnetism, may ultimately be responsible. We examine this issue using 1D stellar models constructed using Modules for Experiments in Stellar Astrophysics (MESA). First, we consider standard models that do not explicitly include rotational/magnetic effects, with convective inhibition modeled by decreasing a depth-independent mixing length theory (MLT) parameter $\alpha_{\text{MLT}}$ (following Cox et al. 1981; Chabrier et al. 2007). We provide formulae linking changes in $\alpha_{\text{MLT}}$ to changes in the interior specific entropy, and hence to the stellar radius. Next, we modify the MLT formulation in MESA to mimic explicitly the influence of rotation and magnetism, using formulations suggested by Stevenson (1979) and MacDonald & Mullan (2014) respectively. We find rapid rotation in these models has a negligible impact on stellar structure, primarily because a star's adiabat, and hence its radius, is predominantly affected by layers near the surface; convection is rapid and largely uninfluenced by rotation there. Magnetic fields, if they influenced convective transport in the manner described by MacDonald & Mullan (2014), could lead to more noticeable radius inflation. Finally, we show that these non-standard effects on stellar structure can be fabricated using a depth-dependent $\alpha_{\text{MLT}}$: a non-magnetic, non-rotating model can be produced that is virtually indistinguishable from one that explicitly parameterizes rotation and/or magnetism using the two formulations above. We provide formulae linking the radially-variable $\alpha_{\text{MLT}}$ to these putative MLT reformulations.

Convection in astrophysical systems must be maintained against dissipation. Although the effects of dissipation are often assumed to be negligible, theory suggests that in strongly stratified convecting fluids, the dissipative heating rate can exceed the luminosity carried by convection. Here we explore this possibility using a series of numerical simulations. We consider two-dimensional numerical models of hydrodynamic convection in a Cartesian layer under the anelastic approximation and demonstrate that the dissipative heating rate can indeed exceed the imposed luminosity. We establish a theoretical expression for the ratio of the dissipative heating rate to the luminosity emerging at the upper boundary, in terms only of the depth of the layer and the thermal scale height. In particular, we show that this ratio is independent of the diffusivities and confirm this with a series of numerical simulations. Our results suggest that dissipative heating may significantly alter the internal dynamics of stars and planets.

We present highlights from a large set of simulations of a hot Jupiter atmosphere, nominally based on HD 209458b, aimed at exploring both the evolution of the deep atmosphere, and the acceleration of the zonal flow or jet. We find the occurrence of a super-rotating equatorial jet is robust to changes in various parameters, and over long timescales, even in the absence of strong inner or bottom boundary drag. This jet is diminished in one simulation only, where we strongly force the deep atmosphere equator-to-pole temperature gradient over long timescales. Finally, although the eddy momentum fluxes in our atmosphere show similarities with the proposed mechanism for accelerating jets on tidally-locked planets, the picture appears more complex. We present tentative evidence for a jet driven by a combination of eddy momentum transport and mean flow.

Evidence of surface magnetism is now observed on an increasing number of cool stars. The detailed manner by which dynamo-generated magnetic fields giving rise to starspots traverse the convection zone still remains unclear. Some insight into this flux emergence mechanism has been gained by assuming bundles of magnetic field can be represented by idealized thin flux tubes (TFTs). Weber & Browning (2016) have recently investigated how individual flux tubes might evolve in a 0.3 solar-mass M dwarf by effectively embedding TFTs in time-dependent flows representative of a fully convective star. We expand upon this work by initiating flux tubes at various depths in the upper 50-75% of the star in order to sample the differing convective flow pattern and differential rotation across this region. Specifically, we comment on the role of differential rotation and time-varying flows in both the suppression and promotion of the magnetic flux emergence process.

Many fully convective stars exhibit a wide variety of surface magnetism, including starspots and chromospheric activity. The manner by which bundles of magnetic field traverse portions of the convection zone to emerge at the stellar surface is not especially well understood. In the Solar context, some insight into this process has been gleaned by regarding the magnetism as consisting partly of idealized thin flux tubes (TFT). Here, we present the results of a large set of TFT simulations in a rotating spherical domain of convective flows representative of a 0.3 solar-mass, main-sequence star. This is the first study to investigate how individual flux tubes in such a star might rise under the combined influence of buoyancy, convection, and differential rotation. A time-dependent hydrodynamic convective flow field, taken from separate 3D simulations calculated with the anelastic equations, impacts the flux tube as it rises. Convective motions modulate the shape of the initially buoyant flux ring, promoting localized rising loops. Flux tubes in fully convective stars have a tendency to rise nearly parallel to the rotation axis. However, the presence of strong differential rotation allows some initially low latitude flux tubes of moderate strength to develop rising loops that emerge in the near-equatorial region. Magnetic pumping suppresses the global rise of the flux tube most efficiently in the deeper interior and at lower latitudes. The results of these simulations aim to provide a link between dynamo-generated magnetic fields, fluid motions, and observations of starspots for fully convective stars.

Observations have suggested that some low-mass stars have larger radii than predicted by 1-D structure models. Some theoretical models have invoked very strong interior magnetic fields (of order 1 MG or more) as a possible cause of such large radii. Whether fields of that strength could in principle by generated by dynamo action in these objects is unclear, and we do not address the matter directly. Instead, we examine whether such fields could remain in the interior of a low mass object for a significant time, and whether they would have any other obvious signatures. First, we estimate timescales for the loss of strong fields by magnetic buoyancy instabilities. We consider a range of field strengths and simple morphologies, including both idealized flux tubes and smooth layers of field. We confirm some of our analytical estimates using thin flux tube magnetohydrodynamic (MHD) simulations of the rise of buoyant fields in a fully-convective M-dwarf. Separately, we consider the Ohmic dissipation of such fields. We find that dissipation provides a complementary constraint to buoyancy: while small-scale, fibril fields might be regenerated faster than they rise, the dissipative heating associated with such fields would in some cases greatly exceed the luminosity of the star. We show how these constraints combine to yield limits on the internal field strength and morphology in low-mass stars. In particular, we find that for stars of 0.3 solar masses, no fields in flux tubes stronger than about 800 kG are simultaneously consistent with both constraints.

We discuss here a purely hydrodynamical mechanism to invert the sign of the kinetic helicity, which plays a key role in determining the direction of propagation of cyclical magnetism in most models of dynamo action by rotating convection. Such propagation provides a prominent, and puzzling constraint on dynamo models. In the Sun, active regions emerge first at mid-latitudes, then appear nearer the equator over the course of a cycle, but most previous global-scale dynamo simulations have exhibited poleward propagation (if they were cyclical at all). Here, we highlight some simulations in which the direction of propagation of dynamo waves is altered primarily by an inversion of the kinetic helicity throughout much of the interior, rather than by changes in the differential rotation. This tends to occur in cases with a low Prandtl number and internal heating, in regions where the local density gradient is relatively small. We analyse how this inversion arises, and contrast it to the case of convection that is either highly columnar (i.e., rapidly rotating) or locally very stratified; in both of those situations, the typical profile of kinetic helicity (negative throughout most of the northern hemisphere) instead prevails.

We are adapting the Global Circulation Model (GCM) of the UK Met Office, the so--called Unified Model (UM), for the study of hot Jupiters. In this work we demonstrate the successful adaptation of the most sophisticated dynamical core, the component of the GCM which solves the equations of motion for the atmosphere, available within the UM, ENDGame (Even Newer Dynamics for General atmospheric modelling of the environment). Within the same numerical scheme ENDGame supports solution to the dynamical equations under varying degrees of simplification. We present results from a simple, shallow (in atmospheric domain) hot Jupiter model (SHJ), and a more realistic (with a deeper atmosphere) HD 209458b test case. For both test cases we find that the large--scale, time--averaged (over the 1200 days prescribed test period), dynamical state of the atmosphere is relatively insensitive to the level of simplification of the dynamical equations. However, problems exist when attempting to reproduce the results for these test cases derived from other models. For the SHJ case the lower (and upper) boundary intersects the dominant dynamical features of the atmosphere meaning the results are heavily dependent on the boundary conditions. For the HD 209458b test case, when using the more complete dynamical models, the atmosphere is still clearly evolving after 1200 days, and in a transient state. Solving the complete (deep atmosphere and non--hydrostatic) dynamical equations allows exchange between the vertical and horizontal momentum of the atmosphere, via Coriolis and metric terms. Subsequently, interaction between the upper atmosphere and the deeper more slowly evolving (radiatively inactive) atmosphere significantly alters the results, and acts over timescales longer than 1200 days.

Young solar-type stars rotate rapidly and many are magnetically active; some undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic spherical harmonic (ASH) code to explore dynamo action achieved in the convective envelope of a solar-type star rotating at 5 times the current solar rotation rate. Striking global-scale magnetic wreaths appear in the midst of the turbulent convection zone and show rich time-dependence. The dynamo exhibits cyclic activity and undergoes quasi-periodic polarity reversals where both the global-scale poloidal and toroidal fields change in sense on a roughly 1500 day time scale. These magnetic activity patterns emerge spontaneously from the turbulent flow and are more organized temporally and spatially than those realized in our previous simulations of the solar dynamo. We assess in detail the competing processes of magnetic field creation and destruction within our simulations that contribute to the global-scale reversals. We find that the mean toroidal fields are built primarily through an $\Omega$-effect, while the mean poloidal fields are built by turbulent correlations which are not necessarily well represented by a simple $\alpha$-effect. During a reversal the magnetic wreaths propagate towards the polar regions, and this appears to arise from a poleward propagating dynamo wave. The primary response in the convective flows involves the axisymmetric differential rotation which shows variations associated with the poleward propagating magnetic wreaths. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent poleward propagating magnetic fields deep below the solar surface. [abridged]

Young solar-type stars rotate rapidly and are very magnetically active. The magnetic fields at their surfaces likely originate in their convective envelopes where convection and rotation can drive strong dynamo action. Here we explore simulations of global-scale stellar convection in rapidly rotating suns using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. We explore one case rotates five times faster than the Sun in detail. This dynamo simulation, called case D5, has repeated quasi-cyclic reversals of global-scale polarity. We compare this case D5 to the broader family of simulations we have been able to explore and discuss how future simulations and observations can advance our understanding of stellar dynamos and magnetism.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for a system rotating at three times the current solar rotation rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in this system. Striking wreaths of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a surprise, for it has been widely believed that such magnetic structures should be disrupted by magnetic buoyancy or turbulent pumping. Thus, many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. We examine how these persistent magnetic wreaths are maintained by dynamo processes and explore whether a classical mean-field $\alpha$-effect explains the regeneration of poloidal field.

When stars like our Sun are young they rotate rapidly and are very magnetically active. We explore dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in these dynamos are organized on global-scales into wreath-like structures that span the convection zone. Wreath-building dynamos can undergo quasi-cyclic reversals of polarity and such behavior is common in the parameter space we have been able to explore. These dynamos do not appear to require tachoclines to achieve their spatial or temporal organization. Wreath-building dynamos are present to some degree at all rotation rates, but are most evident in the more rapidly rotating simulations.

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.

Stars of sufficiently low mass are convective throughout their interiors, and so do not possess an internal boundary layer akin to the solar tachocline. Because that interface figures so prominently in many theories of the solar magnetic dynamo, a widespread expectation had been that fully convective stars would exhibit surface magnetic behavior very different from that realized in more massive stars. Here I describe how recent observations and theoretical models of dynamo action in low-mass stars are partly confirming, and partly confounding, this basic expectation. In particular, I present the results of 3--D MHD simulations of dynamo action by convection in rotating spherical shells that approximate the interiors of 0.3 solar-mass stars at a range of rotation rates. The simulated stars can establish latitudinal differential rotation at their surfaces which is solar-like at ``rapid'' rotation rates (defined within) and anti-solar at slower rotation rates; the differential rotation is greatly reduced by feedback from strong dynamo-generated magnetic fields in some parameter regimes. I argue that this ``flip'' in the sense of differential rotation may be observable in the near future. I also briefly describe how the strength and morphology of the magnetic fields varies with the rotation rate of the simulated star, and show that the maximum magnetic energies attained are compatible with simple scaling arguments.

We review recent insights into the dynamics of the solar convection zone obtained from global numerical simulations, focusing on two recent developments in particular. The first is quasi-cyclic magnetic activity in a long-duration dynamo simulation. Although mean fields comprise only a few percent of the total magnetic energy they exhibit remarkable order, with multiple polarity reversals and systematic variability on time scales of 6-15 years. The second development concerns the maintenance of the meridional circulation. Recent high-resolution simulations have captured the subtle nonlinear dynamical balances with more fidelity than previous, more laminar models, yielding more coherent circulation patterns. These patterns are dominated by a single cell in each hemisphere, with poleward and equatorward flow in the upper and lower convection zone respectively. We briefly address the implications of and future of these modeling efforts.

When our Sun was young it rotated much more rapidly than now. Observations of young, rapidly rotating stars indicate that many possess substantial magnetic activity and strong axisymmetric magnetic fields. We conduct simulations of dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical harmonic (ASH) code to explore the complex coupling between rotation, convection and magnetism. Here we study dynamo action realized in the bulk of the convection zone for two systems, rotating at three and five times the current solar rate. We find that substantial organized global-scale magnetic fields are achieved by dynamo action in these systems. Striking wreathes of magnetism are built in the midst of the convection zone, coexisting with the turbulent convection. This is a great surprise, for many solar dynamo theories have suggested that a tachocline of penetration and shear at the base of the convection zone is a crucial ingredient for organized dynamo action, whereas these simulations do not include such tachoclines. Some dynamos achieved in these rapidly rotating states build persistent global-scale fields which maintain amplitude and polarity for thousands of days. In the case at five times the solar rate, the dynamo can undergo cycles of activity, with fields varying in strength and even changing polarity. As the magnetic fields wax and wane in strength, the primary response in the convective flows involves the axisymmetric differential rotation, which begins to vary on similar time scales. Bands of relatively fast and slow fluid propagate toward the poles on time scales of roughly 500 days. In the Sun, similar patterns are observed in the poleward branch of the torsional oscillations, and these may represent a response to poleward propagating magnetic field deep below the solar surface.

The NSF's Astronomy and Astrophysics Postdoctoral Fellowship (AAPF) is exceptional among the available postdoctoral awards in Astronomy and Astrophysics. The fellowship is one of the few that allows postdoctoral researchers to pursue an original research program, of their own design, at the U.S. institution of their choice. However, what makes this fellowship truly unique is the ability of Fellows to lead an equally challenging, original educational program simultaneously. The legacy of this singular fellowship has been to encourage and advance leaders in the field who are equally as passionate about their own research as they are about sharing that research and their passion for astronomy with students and the public. In this positional paper we address the importance of fellowships like the AAPF to the astronomical profession by identifying the science and educational contributions that Fellows have made to the community. Further, we recommend that fellowships that encourage leading postdoctoral researchers to also become leaders in Astronomy education be continued and expanded.

We review recent advances in modeling global-scale convection and dynamo processes with the Anelastic Spherical Harmonic (ASH) code. In particular, we have recently achieved the first global-scale solar convection simulations that exhibit turbulent pumping of magnetic flux into a simulated tachocline and the subsequent organization and amplification of toroidal field structures by rotational shear. The presence of a tachocline not only promotes the generation of mean toroidal flux, but it also enhances and stabilizes the mean poloidal field throughout the convection zone, promoting dipolar structure with less frequent polarity reversals. The magnetic field generated by a convective dynamo with a tachocline and overshoot region is also more helical overall, with a sign reversal in the northern and southern hemispheres. Toroidal tachocline fields exhibit little indication of magnetic buoyancy instabilities but may be undergoing magneto-shear instabilities.

We present magnetic flux measurements in seven rapidly rotating M dwarfs. Our sample stars have X-ray and H-alpha emission indicative of saturated emission, i.e., emission at a high level independent of rotation rate. Our measurements are made using near-infrared FeH molecular spectra observed with HIRES at Keck. Because of their large convective overturn times, the rotation velocity of M stars with small Rossby numbers is relatively slow and does not hamper the measurement of Zeeman splitting. The Rossby numbers of our sample stars are as small as 0.01. All our sample stars exhibit magnetic flux of kilo-Gauss strength. We find that the magnetic flux saturates in the same regime as saturation of coronal and chromospheric emission, at a critical Rossby number of around 0.1. The filling factors of both field and emission are near unity by then. We conclude that the strength of surface magnetic fields remains independent of rotation rate below that; making the Rossby number yet smaller by a factor of ten has little effect. These saturated M-star dynamos generate an integrated magnetic flux of roughly 3 kG, with a scatter of about 1 kG. The relation between emission and flux also has substantial scatter.

Today, the generation of magnetic fields in solar-type stars and its relation to activity and rotation can coherently be explained, although it is certainly not understood in its entirety. Rotation facilitates the generation of magnetic flux that couples to the stellar wind, slowing down the star. There are still many open questions, particularly at early phases (young age), and at very low mass. It is vexing that rotational braking becomes inefficient at the threshold to fully convective interiors, although no threshold in magnetic activity is seen, and the generation of large scale magnetic fields is still possible for fully convective stars. This article briefly outlines our current understanding of the rotation-magnetic field relation.

In the solar convection zone, rotation couples with intensely turbulent convection to drive a strong differential rotation and achieve complex magnetic dynamo action. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the global-scale patterns of convection in such stars and the flows of differential rotation and meridional circulation which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of localized convection emerges. Convection near the equator in these systems is dominated by one or two nests in longitude of locally enhanced convection, with quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation Delta Omega grows with more rapid rotation while the relative shear Delta Omega/Omega_0 decreases. In contrast, at more rapid rotation the meridional circulations decrease in energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

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.

In the solar convection zone, rotation couples with intensely turbulent convection to build global-scale flows of differential rotation and meridional circulation. Our sun must have rotated more rapidly in its past, as is suggested by observations of many rapidly rotating young solar-type stars. Here we explore the effects of more rapid rotation on the patterns of convection in such stars and the global-scale flows which are self-consistently established. The convection in these systems is richly time dependent and in our most rapidly rotating suns a striking pattern of spatially localized convection emerges. Convection near the equator in these systems is dominated by one or two patches of locally enhanced convection, with nearly quiescent streaming flow in between at the highest rotation rates. These active nests of convection maintain a strong differential rotation despite their small size. The structure of differential rotation is similar in all of our more rapidly rotating suns, with fast equators and slower poles. We find that the total shear in differential rotation, as measured by latitudinal angular velocity contrast, Delta_Omega, increases with more rapid rotation while the relative shear, Delta_Omega/Omega, decreases. In contrast, at more rapid rotation the meridional circulations decrease in both energy and peak velocities and break into multiple cells of circulation in both radius and latitude.

We present three-dimensional nonlinear magnetohydrodynamic simulations of the interiors of fully convective M-dwarfs. Our models consider 0.3 solar-mass stars using the Anelastic Spherical Harmonic code, with the spherical computational domain extending from 0.08-0.96 times the overall stellar radius. Like previous authors, we find that fully convective stars can generate kG-strength magnetic fields (in rough equipartition with the convective flows) without the aid of a tachocline of shear. Although our model stars are everywhere unstably stratified, the amplitudes and typical pattern sizes of the convective flows vary strongly with radius, with the outer regions of the stars hosting vigorous convection and field amplification while the deep interiors are more quiescent. Modest differential rotation is established in hydrodynamic calculations, but -- unlike in some prior work --strongly quenched in MHD simulations because of the Maxwell stresses exerted by the dynamo-generated magnetic fields. Despite the lack of strong differential rotation, the magnetic fields realized in the simulations possess significant mean (axisymmetric) components, which we attribute partly to the strong influence of rotation upon the slowly overturning flows.

Recent observations indicate that fully convective stars can effectively build magnetic fields without the aid of a tachocline of shear, that those fields can possess large-scale components, and that they may sense the effects of rotation. Motivated by these puzzles, we present global three-dimensional simulations of convection and dynamo action in the interiors of fully convective M-dwarfs of 0.3 solar masses. We use the Anelastic Spherical Harmonic (ASH) code, adopting a spherical computational domain that extends from 0.08-0.96 times the overall stellar radius. We find that such fully convective stars can generate magnetic fields of several kG strength, roughly in equipartition with the convective flows. Differential rotation is established in hydrodynamic progenitor calculations, but strongly quenched in MHD simulations because of strong Maxwell stresses exerted by the magnetic fields. Despite the absence of interior angular velocity contrasts, the magnetic fields possess strong mean (axisymmetric) components, which we attribute partly to the very strong influence of rotation upon the slowly overturning flows.

Core convection and dynamo activity deep within rotating A-type stars of 2 solar masses are studied with 3--D nonlinear simulations. Our modeling considers the inner 30% by radius of such stars, thus capturing within a spherical domain the convective core and a modest portion of the surrounding radiative envelope. The MHD equations are solved using the ASH code to examine turbulent flows and magnetic fields, both of which exhibit intricate time dependence. By introducing small seed magnetic fields into our progenitor hydrodynamic models rotating at one and four times the solar rate, we assess here how the vigorous convection can amplify those fields and sustain them against ohmic decay. Dynamo action is indeed realized, ultimately yielding magnetic fields that are in energy equipartion with the flow. Such magnetism reduces the differential rotation obtained in the progenitors, partly by Maxwell stresses that transport angular momentum poleward and oppose the Reynolds stresses in the latitudinal balance. In contrast, in the radial direction we find that the Maxwell and Reynolds stresses may act together to transport angular momentum. The central columns of slow rotation established in the progenitors are weakened, with the differential rotation waxing and waning in strength as the simulations evolve. We assess the morphology of the flows and magnetic fields, their complex temporal variations, and the manner in which dynamo action is sustained. Differential rotation and helical convection are both found to play roles in giving rise to the magnetic fields. The magnetism is dominated by strong fluctuating fields throughout the core, with the axisymmetric (mean) fields there relatively weak.

We present the first results from three-dimensional spherical shell simulations of magnetic dynamo action realized by turbulent convection penetrating downward into a tachocline of rotational shear. This permits us to assess several dynamical elements believed to be crucial to the operation of the solar global dynamo, variously involving differential rotation resulting from convection, magnetic pumping, and amplification of fields by stretching within the tachocline. The simulations reveal that strong axisymmetric toroidal magnetic fields (about 3000 G in strength) are realized within the lower stable layer, unlike in the convection zone where fluctuating fields are predominant. The toroidal fields in the stable layer possess a striking persistent antisymmetric parity, with fields in the northern hemisphere largely of opposite polarity to those in the southern hemisphere. The associated mean poloidal magnetic fields there have a clear dipolar geometry, but we have not yet observed any distinctive reversals or latitudinal propagation. The presence of these deep magnetic fields appears to stabilize the sense of mean fields produced by vigorous dynamo action in the bulk of the convection zone.

We present the results of 3--D nonlinear simulations of magnetic dynamo action by core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, with the spherical domain thereby encompassing the convective core and a portion of the surrounding radiative envelope. The compressible Navier-Stokes equations, subject to the anelastic approximation, are solved to examine highly nonlinear flows that span multiple scale heights, exhibit intricate time dependence, and admit magnetic dynamo action. Small initial seed magnetic fields are found to be amplified greatly by the convective and zonal flows. The central columns of strikingly slow rotation found in some of our progenitor hydrodynamic simulations continue to be realized in some simulations to a lesser degree, with such differential rotation arising from the redistribution of angular momentum by the nonlinear convection and magnetic fields. We assess the properties of the magnetic fields thus generated, the extent of convective penetration, the magnitude of the differential rotation, and the excitation of gravity waves within the radiative envelope.

We present the results of 3--D simulations of core convection within A-type stars of 2 solar masses, at a range of rotation rates. We consider the inner 30% by radius of such stars, thereby encompassing the convective core and some of the surrounding radiative envelope. We utilize our anelastic spherical harmonic (ASH) code, which solves the compressible Navier-Stokes equations in the anelastic approximation, to examine highly nonlinear flows that can span multiple scale heights. The cores of these stars are found to rotate differentially, with central cylindrical regions of strikingly slow rotation achieved in our simulations of stars whose convective Rossby number (R_{oc}) is less than unity. Such differential rotation results from the redistribution of angular momentum by the nonlinear convection that strongly senses the overall rotation of the star. Penetrative convective motions extend into the overlying radiative zone, yielding a prolate shape (aligned with the rotation axis) to the central region in which nearly adiabatic stratification is achieved. This is further surrounded by a region of overshooting motions, the extent of which is greater at the equator than at the poles, yielding an overall spherical shape to the domain experiencing at least some convective mixing. We assess the overshooting achieved as the stability of the radiative exterior is varied, and the weak circulations that result in that exterior. The convective plumes serve to excite gravity waves in the radiative envelope, ranging from localized ripples of many scales to some remarkable global resonances.

The advent of massively parallel supercomputing has begun to permit explicit 3--D simulations of turbulent convection occurring within the cores of early-type main sequence stars. Such studies should complement the stellar structure and evolution efforts that have so far largely employed 1--D nonlocal mixing length descriptions for the transport, mixing and overshooting achieved by core convection. We have turned to A-type stars as representative of many of the dynamical challenges raised by core convection within rotating stars. The differential rotation and meridional circulations achieved deep within the star by the convection, the likelihood of sustained magnetic dynamo action there, and the bringing of fresh fuel into the core by overshooting motions, thereby influencing main sequence lifetimes, all constitute interesting dynamical questions that require detailed modelling of global-scale convection. Using our anelastic spherical harmonic (ASH) code tested on the solar differential rotation problem, we have conducted a series of 3--D spherical domain simulations that deal with a simplified description of the central regions of rotating A-type stars, i.e a convectively unstable core is surrounded by a stable radiative envelope. A sequence of 3--D simulations are used to assess the properties of the convection (its global patterns, differential rotation, meridional circulations, extent and latitudinal variation of the overshooting) as transitions are made between laminar and turbulent states by changing the effective diffusivities, rotation rates, and subadiabaticity of the radiative exterior. We report on the properties deduced from these models for both the extent of penetration and the profile of rotation sustained by the convection.

We have developed a code that models the formation, destruction, radiative transfer, and vibrational/rotational excitation of H_2 in a detailed fashion. We discuss how such codes, together with FUSE observations of H_2 in diffuse and translucent lines of sight, may be used to infer various physical parameters. We illustrate the effects of changes in the major physical parameters (UV radiation field, gas density, metallicity), and we point out the extent to which changes in one parameter may be mirrored by changes in another. We provide an analytic formula for the molecular fraction, f_H2, as a function of cloud column density, radiation field, and grain formation rate of H_2. Some diffuse and translucent lines of sight may be concatenations of multiple distinct clouds viewed together. Such situations can give rise to observables that agree with the data, complicating the problem of uniquely identifying one set of physical parameters with a line of sight. Finally, we illustrate the application of our code to an ensemble of data, such as the FUSE survey of H_2 in the Large and Small Magellanic Clouds (LMC/SMC), in order to constrain the elevated UV radiation field intensity and reduced grain formation rate of H_2 in those low- metallicity environments.

We describe a moderate-resolution FUSE survey of H2 along 70 sight lines to the Small and Large Magellanic Clouds, using hot stars as background sources. FUSE spectra of 67% of observed Magellanic Cloud sources (52% of LMC and 92% of SMC) exhibit absorption lines from the H2 Lyman and Werner bands between 912 and 1120 A. Our survey is sensitive to N(H2) >= 10^14 cm^-2; the highest column densities are log N(H2) = 19.9 in the LMC and 20.6 in the SMC. We find reduced H2 abundances in the Magellanic Clouds relative to the Milky Way, with average molecular fractions <f_H2> = 0.010 (+0.005, -0.002) for the SMC and <f_H2> = 0.012 (+0.006, -0.003) for the LMC, compared with <f_H2> = 0.095 for the Galactic disk over a similar range of reddening. The dominant uncertainty in this measurement results from the systematic differences between 21 cm radio emission and Lya in pencil-beam sight lines as measures of N(HI). These results imply that the diffuse H2 masses of the LMC and SMC are 8 x 10^6 Msun and 2 x 10^6 Msun, respectively, 2% and 0.5% of the H I masses derived from 21 cm emission measurements. The LMC and SMC abundance patterns can be reproduced in ensembles of model clouds with a reduced H2 formation rate coefficient, R ~ 3 x 10^-18 cm^3 s^-1, and incident radiation fields ranging from 10 - 100 times the Galactic mean value. We find that these high-radiation, low-formation-rate models can also explain the enhanced N(4)/N(2) and N(5)/N(3) rotational excitation ratios in the Clouds. We use H2 column densities in low rotational states (J = 0 and 1) to derive a mean kinetic and/or rotational temperature <T_01> = 82 +/- 21 K for clouds with N(H2) >= 10^16 cm^-2, similar to Galactic gas. We discuss the implications of this work for theories of star formation in low-metallicity environments. [Abstract abridged]