Massive stars are the drivers of star formation and galactic dynamics due to their relatively short lives and explosive demises, thus impacting all of astrophysics. Since they are so impactful on their environments, through their winds on the main sequence and their ultimate supernovae, it is crucial to understand how they evolve. Recent photometric observations with space-based platforms such as CoRoT, K2, and now TESS have permitted access to their interior dynamics through asteroseismology, while ground-based spectropolarimetric measurements such as those of ESPaDOnS have given us a glimpse at their surface magnetic fields. The dynamics of massive stars involve a vast range of scales. Extant methods can either capture the long-term structural evolution or the short-term dynamics such as convection, magnetic dynamos, and waves due to computational costs. Thus, many mysteries remain regarding the impact of such dynamics on stellar evolution, but they can have strong implications both for how they evolve and what they leave behind when they die. Some of these dynamics including rotation, tides, and magnetic fields have been addressed in recent work, which is reviewed in this paper.

The dissipation of the kinetic energy of wave-like tidal flows within the convective envelope of low-mass stars is one of the key physical mechanisms that shapes the orbital and rotational dynamics of short-period exoplanetary systems. Although low-mass stars are magnetically active objects, the question of how the star's magnetic field impacts large-scale tidal flows and the excitation, propagation and dissipation of tidal waves still remains open. Our goal is to investigate the impact of stellar magnetism on the forcing of tidal waves, and their propagation and dissipation in the convective envelope of low-mass stars as they evolve. We have estimated the amplitude of the magnetic contribution to the forcing and dissipation of tidally induced magneto-inertial waves throughout the structural and rotational evolution of low-mass stars (from M to F-type). For this purpose, we have used detailed grids of rotating stellar models computed with the stellar evolution code STAREVOL. The amplitude of dynamo-generated magnetic fields is estimated via physical scaling laws at the base and the top of the convective envelope. We find that the large-scale magnetic field of the star has little influence on the excitation of tidal waves in the case of nearly-circular orbits and coplanar hot-Jupiter planetary systems, but that it has a major impact on the way waves are dissipated. Our results therefore indicate that a full magneto-hydrodynamical treatment of the propagation and dissipation of tidal waves is needed to properly assess the impact of star-planet tidal interactions throughout the evolutionary history of low-mass stars hosting short-period massive planets.

Aims: We investigate from a theoretical perspective if space asteroseismology can be used to distinguish between different thermal structures and shapes of the near-core mixing profiles for different types of coherent oscillation modes in massive stars with convective cores, and if this capacity depends on the evolutionary stage of the models along the main sequence. Methods: We compute 1D stellar structure and evolution models for four different prescriptions of the mixing and temperature gradient in the near-core region. Their effect on the frequencies of dipole prograde gravity modes in both Slowly Pulsating B and $\beta$ Cep stars is investigated, as well as for pressure modes in $\beta$ Cep stars. Results: A comparison between the mode frequencies of the different models at various stages during the main sequence evolution reveals that they are more sensitive to a change in temperature gradient than to the exact shape of the mixing profile in the near-core region. Depending on the duration of the observed light curve, one can distinguish between either just the temperature gradient, or also between the shapes of the mixing coefficient. The relative frequency differences are in general larger for more evolved models, and are largest for the higher-frequency pressure modes in $\beta$ Cep stars. Conclusions:In order to unravel the core boundary mixing and thermal structure of the near-core region, one must have asteroseismic masses and radii with $\sim 1\%$ relative precision for hundreds of stars.

Rotational scaling relationships are examined for the degree of equipartition between magnetic and kinetic energies in stellar convection zones. These scaling relationships are approached from two paradigms, with first a glance at scaling relationship built upon an energy-balance argument and second a look at a force-based scaling. The latter implies a transition between a nearly-constant inertial scaling when in the asymptotic limit of minimal diffusion and magnetostrophy, whereas the former implies a weaker scaling with convective Rossby number. Both scaling relationships are then compared to a suite of 3D convective dynamo simulations with a wide variety of domain geometries, stratifications, and range of convective Rossby numbers.

A monomodal model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale as a function of rotation rate as well as with thermal and viscous diffusivities. The convection model is used as a boundary condition for a linearization of the equations of motion in the transition region between convectively unstable and stably-stratified regions, yielding the depth to which convection penetrates into the stable region and establishing a relationship between that depth and the local convective Rossby number, diffusivity, and pressure scale height of those flows. Upward and downward penetrative convection have a similar scaling with rotation rate and diffusivities, but they depend differently upon the pressure scale height due to the differing energetic processes occurring in convective cores of early-type stars versus convective envelopes of late-type stars.

Gravito-inertial waves can be excited at the interface of convective and radiative regions and by the Reynolds stresses in the bulk of the convection zone. The magnitude of their energy flux will therefore vary with the properties of the convection. To assess how convection changes with rotation, a simplified local monomodal model for rotating convection is presented that provides the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale as a function of the convective Rossby number as well as with thermal and viscous diffusivities. In the context of this convection model, two models for assessing the gravito-inertial wave flux are considered: an interfacial model and a full treatment of the Reynolds stress impact on the waves. It is found that there are regimes where the sub-inertial waves may carry a significant energy flux relative to pure gravity waves that depend upon the convective Rossby number, the ratio of the buoyancy time-scale in the stable region to the convective overturning time, and the wave frequency.

Differential rotation plays a key role in stellar evolution by triggering hydrodynamical instabilities and large-scale motions that induce transport of chemicals and angular momentum and by modifying the propagation and the frequency spectrum of gravito-inertial waves. It is thus crucial to investigate its effect on the propagation of gravity waves to build reliable seismic diagnostic tools, especially for fast rotating stars, where perturbative treatments of rotation fail. Generalising a previous work done in the case of uniform rotation, we derived a local dispersion relation for gravity waves in a differentially rotating star, taking the full effect of rotation (both Coriolis and centrifugal accelerations) into account. Then we modelled the propagation of axisymmetric waves as the propagation of rays. This allowed us to efficiently probe the properties of the waves in various regimes of differential rotation.

The energy dissipation of wave-like tidal flows in the convective envelope of low-mass stars is one of the key physical mechanisms that shape the orbital and rotational dynamics of short-period planetary systems. Tidal flows, and the excitation, propagation, and dissipation of tidally-induced inertial waves can be modified by stellar magnetic fields (e.g., Wei 2016, 2018, Lin and Ogilvie 2018). It is thus important to assess for which stars, at which location of their internal structure, and at which phase of their evolution, one needs to take into account the effects of magnetic fields on tidal waves. Using scaling laws that provide the amplitude of dynamo-generated magnetic fields along the rotational evolution of these stars (e.g., Christensen et al. 2009, Brun et al. 2015), combined with detailed grids of stellar rotation models (e.g., Amard et al. 2016), we examine the influence of a magnetic field on tidal forcing and dissipation near the tachocline of solar-like stars. We show that full consideration of magnetic fields is required to compute tidal dissipation, but not necessarily for tidal forcing.

The purpose of this short paper is to give a summary of the derivation of a perturbative method to assess the impact of an arbitrary 3D magnetic field configuration, whose energy is everywhere small relative to the gravitational binding energy, on the form and value of the eigenfrequencies and eigenfunctions of stellar oscillations.

A simplified model for stellar and planetary convection is derived for the magnitude of the rms velocity, degree of superadiabaticity, and characteristic length scale with Rossby number as well as with thermal and viscous diffusivities. Integrating the convection model into a linearization of the dynamics in the transition region between convectively unstable and stably-stratified region yields a Rossby number, diffusivity, and pressure scale height dependent convective penetration depth into the stable region. This may have important consequences for mixing along the evolution of rotating stars.

Asteroseismology has reached a level of accuracy that may allow us to detect the effect of a deep magnetic field on oscillation modes. We thus aim to develop an asymptotic theory for short-wavelength waves in radiative zones of rotating stars with a general magnetic field topology (toroidal and poloidal). A dispersion relation is derived for a uniformly rotating star with a low-amplitude, axisymmetric magnetic field. The parameter space of this model is explored with a Hamiltonian ray-tracing method. The features of the gravito-inertial waves modified by the magnetic field are studied with different magnetic topologies.

Context. Differential rotation has a strong influence on stellar internal dynamics and evolution, notably by triggering hydrodynamical instabilities, by interacting with the magnetic field, and more generally by inducing transport of angular momentum and chemical elements. Moreover, it modifies the way waves propagate in stellar interiors and thus the frequency spectrum of these waves, the regions they probe, and the transport they generate. Aims. We investigate the impact of a general differential rotation (both in radius and latitude) on the propagation of axisymmetric gravito-inertial waves. Methods. We use a small-wavelength approximation to obtain a local dispersion relation for these waves. We then describe the propagation of waves thanks to a ray model that follows a Hamiltonian formalism. Finally, we numerically probe the properties of these gravito-inertial rays for different regimes of radial and latitudinal differential rotation. Results. We derive a local dispersion relation that includes the effect of a general differential rotation. Subsequently, considering a polytropic stellar model, we observe that differential rotation allows for a large variety of resonant cavities that can be probed by gravito-inertial waves. We identify that for some regimes of frequency and differential rotation, the properties of gravito-inertial rays are similar to those found in the uniformly rotating case. Furthermore, we also find new regimes specific to differential rotation, where the dynamics of rays is chaotic. Conclusions. As a consequence, we expect modes to follow the same trend. Some parts of oscillation spectra corresponding to regimes similar to those of the uniformly rotating case would exhibit regular patterns, while parts corresponding to the new regimes would be mostly constituted of chaotic modes with a spectrum rather characterised by a generic statistical distribution.

We seek to characterize how the change of global rotation rate influences the overall dynamics and large scale flows arising in the convective envelopes of stars covering stellar spectral types from early G to late K. We do so through numerical simulations with the ASH code, where we consider stellar convective envelopes coupled to a radiative interior with various global properties. As solar-like stars spin down over the course of their main sequence evolution, such change must have a direct impact on their dynamics and rotation state. We indeed find that three main states of rotation may exist for a given star: anti-solar-like (fast poles, slow equator), solar-like (fast equator, slow poles), or a cylindrical rotation profile. Under increasingly strict rotational constraints, the latter profile can further evolve into a Jupiter-like profile, with alternating prograde and retrograde zonal jets. We have further assessed how far the convection and meridional flows overshoot into the radiative zone and investigated the morphology of the established tachocline. Using simple mixing length arguments, we are able to construct a scaling of the fluid Rossby number $R_{of} = \tilde{\omega}/2\Omega_* \sim \tilde{v}/2\Omega_* R_*$, which we calibrate based on our 3-D ASH simulations. We can use this scaling to map the behavior of differential rotation versus the global parameters of stellar mass and rotation rate. Finally, we isolate a region on this map ($R_{of} \gtrsim 1.5-2$) where we posit that stars with an anti-solar differential rotation may exist in order to encourage observers to hunt for such targets.

This paper provides a brief look at dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Two simple models are examined, where one that assumes magnetostrophy and another that includes the effects of inertia. These models are then compared to a suite of convective dynamo simulations of the convective core of a main-sequence B-type star and applied to its later evolutionary stages.

The recent observational, theoretical, and numerical progress made in understanding stellar magnetism is discussed. Particularly, this review will cover the physical processes thought to be at the origin of these magnetic fields and their variability, namely dynamo action arising from the interaction between convection, rotation, radiation and magnetic fields. Some care will be taken to cover recent analytical advances regarding the dynamics and magnetism of radiative interiors, including some thoughts on the role of a tachocline. Moreover, recent and rapidly advancing numerical modeling of convective dynamos will be discussed, looking at rapidly rotating convective systems, grand minima and scaling laws for magnetic field strength. These topics are linked to observations or their observational implications.

This paper provides a brief overview of dynamo scaling relationships for the degree of equipartition between magnetic and kinetic energies. Three basic approaches are adopted to explore these scaling relationships, with a first look at two simple models: one assuming magnetostrophy and another that includes the effects of inertia. Next, a third scaling relationship is derived that utilizes the assumptions that the dynamo possesses two integral spatial scales and that it is driven by the balance of buoyancy work and ohmic dissipation as studied in Davidson 2013. The results of which are then compared to a suite of convective dynamo simulations that possess a fully convective domain with a weak density stratification and that captured the behavior of the resulting dynamo for a range of convective Rossby numbers (Augustson et al. 2016).

This paper provides a brief overview of the formation of stellar fossil magnetic fields and what potential instabilities may occur given certain configurations of the magnetic field. In particular, a purely magnetic instability can occur for poloidal, toroidal, and mixed poloidal-toroidal axisymmetric magnetic field configurations as originally studied in Tayler (1973), Markey & Tayler (1973), and Tayler (1980). However, most of the magnetic field configurations observed at the surface of massive stars are non- axisymmetric. Thus, extending earlier studies of the axisymmetric Tayler instability in spherical geometry (Goossens, 1980), we introduce a formulation for the global change in the potential energy contained in a convectively-stable region given an arbitrary Lagrangian perturbation, which permits the inclusion of both axisymmetric and non-axisymmetric magnetic fields. With this tool in hand, a path is shown by which more general stability criterion can be established.

The dynamo action achieved in the convective cores of main-sequence massive stars is explored here through 3D global simulations of convective core dynamos operating within a young $\Msun{10}$ B-type star, using the anelastic spherical harmonic (ASH) code. These simulations capture the inner 65\% of this star by radius, encompassing the convective nuclear-burning core (about 23\% by radius) and a portion of the overlying radiative envelope. Eight rotation rates are considered, ranging from 0.05\% to 16\% of the surface breakup velocity, thereby capturing both convection that barely senses the effects of rotation and other situations in which the Coriolis forces are prominent. The vigorous dynamo action realized within all of these turbulent convective cores builds magnetic fields with peak strengths exceeding a megagauss, with the overall magnetic energy (ME) in the faster rotators reaching super-equipartition levels compared to the convective kinetic energy (KE). The core convection typically involves turbulent columnar velocity structures roughly aligned with the rotation axis, with magnetic fields threading through these rolls and possessing complex linkages throughout the core. The very strong fields are able to coexist with the flows without quenching them through Lorentz forces. The velocity and magnetic fields achieve such a state by being nearly co-aligned, and with peak magnetic islands being somewhat displaced from the fastest flows as the intricate evolution proceeds. As the rotation rate is increased, the primary force balance shifts from nonlinear advection balancing Lorentz forces to a magnetostrophic balance between Coriolis and Lorentz forces.

We investigate the role of magnetic helicity in promoting cyclic magnetic activity in a global, 3D, magnetohydrodynamic (MHD) simulation of a convective dynamo. This simulation is characterized by coherent bands of toroidal field that exist within the convection zone, with opposite polarities in the northern and southern hemispheres. Throughout most of the cycle, the magnetic helicity in these bands is negative in the northern hemisphere and positive in the southern hemisphere. However, during the declining phase of each cycle, this hemispheric rule reverses. We attribute this to a global restructuring of the magnetic topology that is induced by the interaction of the bands across the equator. This band interaction appears to be ultimately responsible for, or at least associated with, the decay and subsequent reversal of both the toroidal bands and the polar fields. We briefly discuss the implications of these results within the context of solar observations, which also show some potential evidence for toroidal band interactions and helicity reversals.

Several lines of evidence suggest that the velocity amplitude in global simulations of solar convection, U, may be systematically over-estimated. Motivated by these recent results, we explore the factors that determine U and we consider how these might scale to solar parameter regimes. To this end, we decrease the thermal diffusivity $\kappa$ along two paths in parameter space. If the kinematic viscosity $\nu$ is decreased proportionally with $\kappa$ (fixing the Prandtl number $P_r = \nu/\kappa$), we find that U increases but asymptotes toward a constant value, as found by Featherstone & Hindman (2016). However, if $\nu$ is held fixed while decreasing $\kappa$ (increasing $P_r$), we find that U systematically decreases. We attribute this to an enhancement of the thermal content of downflow plumes, which allows them to carry the solar luminosity with slower flow speeds. We contrast this with the case of Rayleigh-Benard convection which is not subject to this luminosity constraint. This dramatic difference in behavior for the two paths in parameter space (fixed $P_r$ or fixed $\nu$) persists whether the heat transport by unresolved, near-surface convection is modeled as a thermal conduction or as a fixed flux. The results suggest that if solar convection can operate in a high-$P_r$ regime, then this might effectively limit the velocity amplitude. Small-scale magnetism is a possible source of enhanced viscosity that may serve to achieve this high-$P_r$ regime.

The 3-D magnetohydrodynamic (MHD) Anelastic Spherical Harmonic (ASH) code, using slope-limited diffusion, is employed to capture convective and dynamo processes achieved in a global-scale stellar convection simulation for a model solar-mass star rotating at three times the solar rate. The dynamo generated magnetic fields possesses many time scales, with a prominent polarity cycle occurring roughly every 6.2 years. The magnetic field forms large-scale toroidal wreaths, whose formation is tied to the low Rossby number of the convection in this simulation. The polarity reversals are linked to the weakened differential rotation and a resistive collapse of the large-scale magnetic field. An equatorial migration of the magnetic field is seen, which is due to the strong modulation of the differential rotation rather than a dynamo wave. A poleward migration of magnetic flux from the equator eventually leads to the reversal of the polarity of the high-latitude magnetic field. This simulation also enters an interval with reduced magnetic energy at low latitudes lasting roughly 16 years (about 2.5 polarity cycles), during which the polarity cycles are disrupted and after which the dynamo recovers its regular polarity cycles. An analysis of this grand minimum reveals that it likely arises through the interplay of symmetric and antisymmetric dynamo families. This intermittent dynamo state potentially results from the simulations relatively low magnetic Prandtl number. A mean-field-based analysis of this dynamo simulation demonstrates that it is of the $\alpha$-$\Omega$ type. The time scales that appear to be relevant to the magnetic polarity reversal are also identified.

The global-scale dynamo action achieved in a simulation of a Sun-like star rotating at thrice the solar rate is assessed. The 3-D MHD Anelastic Spherical Harmonic (ASH) code, augmented with a viscosity minimization scheme, is employed to capture convection and dynamo processes in this G-type star. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. It is found that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals occurring roughly every 6.2 years. The magnetic energy also rises and falls with a regular period. The magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the processes leading to polarity reversals are more complex, appearing to arise from the interaction of convection with the mean toroidal fields. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation, and potentially to a nonlinear dynamo wave. This simulation also enters a grand minimum lasting roughly 20~years, after which the dynamo recovers its regular polarity cycles.

The Kepler mission has collected light curves for almost 4 years. The excellent quality of these data has allowed us to probe the structure and the dynamics of the stars using asteroseismology. With the length of data available, we can start to look for magnetic activity cycles. The Kepler data obtained for the F star, KIC 12009504, shows a rotation period of 9.5 days and additional variability that could be due to the magnetic activity of the star. Here we present recent and preliminary 3D global-scale dynamo simulations of this star with the ASH and STELEM codes, capturing a substantial portion of the convection and the stable radiation zone below it. These simulations reveal a multi-year activity cycle whose length tentatively depends upon the width of the tachocline present in the simulation. Furthermore, the presence of a magnetic field and the dynamo action taking place in the convection zone appears to help confine the tachocline, but longer simulations will be required to confirm this.

We assess the global-scale dynamo action achieved in a simulation of a sun-like star rotating at three times the solar rate. The 3-D MHD Anelastic Spherical Harmonic code, using slope-limited diffusion, is employed to capture convection and dynamo processes in such a young sun. The simulation is carried out in a spherical shell that encompasses 3.8 density scale heights of the solar convection zone. We find that dynamo action with a high degree of time variation occurs, with many periodic polarity reversals every 6.2 years. The magnetic energy also rises and falls with a regular period, with two magnetic energy cycles required to complete a polarity cycle. These magnetic energy cycles arise from a Lorentz-force feedback on the differential rotation, whereas the polarity reversals are present due to the spatial separation of the equatorial and polar dynamos. Moreover, an equatorial migration of toroidal field is found, which is linked to the changing differential rotation and to a near-surface shear layer. This simulation also enters a grand minimum lasting roughly 20 years, after which the dynamo recovers its regular polarity cycles.

We examine the role of small-scale granulation in helping to drive supergranulation and even larger scales of convection. The granulation is modeled as localized cooling events introduced at the upper boundary of a 3-D simulation of compressible convection in a rotating spherical shell segment. With a sufficient number of stochastic cooling events compared to uniform cooling, we find that supergranular scales are realized, along with a differential rotation that becomes increasingly solar-like.

As we approach solar convection simulations that seek to model the interaction of small-scale granulation and supergranulation and even larger scales of convection within the near-surface shear layer (NSSL), the treatment of the boundary conditions and minimization of sub-grid scale diffusive processes become increasingly crucial. We here assess changes in the dynamics and the energy flux balance of the flows established in rotating spherical shell segments that capture much of the NSSL with the Curved Spherical Segment (CSS) code using two different diffusion schemes. The CSS code is a new massively parallel modeling tool capable of simulating 3-D compressible MHD convection with a realistic solar stratification in rotating spherical shell segments.

Main-sequence massive stars possess convective cores that likely harbor strong dynamo action. To assess the role of core convection in building magnetic fields within these stars, we employ the 3-D anelastic spherical harmonic (ASH) code to model turbulent dynamics within a 10 solar mass main-sequence (MS) B-type star rotating at 4 times the solar rate. We find that strong (900 kG) magnetic fields arise within the turbulence of the core and penetrate into the stably stratified radiative zone. These fields exhibit complex, time-dependent behavior including reversals in magnetic polarity and shifts between which hemisphere dominates the total magnetic energy.