We solve for waves in an isothermal, stratified medium with a magnetic field that points along a direction perpendicular to that of gravity and varies exponentially in the direction of gravity. For waves propagating along the magnetic field, we calculate approximate dispersion relation as a function of the magnetic field strength. We also find exact solutions for two different cases: (a) waves propagating along the direction of the magnetic field and (b) waves propagating along the direction of gravity. In each of these cases, we find solutions in terms of either confluent hypergeometric functions or Gauss hypergeometric functions depending on whether the ratio of the scale height of the magnetic field over the density scale height is equal to two or not.

Long-term solar magnetic activity reconstructions indicate the solar dynamo operates in two distinct -- grand minimum and regular activity -- modes. By employing bifurcation analysis of a reduced dynamo model we establish this to be a consequence of dynamo hysteresis. We reproduce the observed bimodal distribution of sunspots, but only for subcritical dynamos. A theoretical framework consistent with this finding explains observations of an abrupt midlife transition in stellar activity, characterized by reduced angular momentum loss rates and breakdown of gyrochronology relations. Our study indicates that an evolving dynamo bridges a diversity of phenomena in solar-like stars across their lifetime.

We explore a reduced Babcock-Leighton (BL) dynamo model based on delay differential equations using numerical bifurcation analysis. This model reveals hysteresis, seen in the recent mean-field dynamo model and the direct numerical simulations of turbulent dynamos. The BL model with 'magnetic noise' as an additional weak-source of the poloidal field recovers the solar cycle every time from grand minima, which BL source alone cannot do. The noise-incorporated model exhibits a bimodal distribution of toroidal field energy confirming two modes of solar activity. It also shows intermittency and reproduces phase space collapse, an experimental signature of the Maunder Minimum. The occurrence statistics of grand minima in our model agree reasonably well with the observed statistics in the reconstructed sunspot number. Finally, we demonstrate that the level of magnetic noise controls the duration of grand minima and even has a handle over its waiting period, suggesting a triggering effect of grand minima by the noise and thus shutting down the global dynamo. Therefore, we conclude that the 'magnetic noise' due to small-scale turbulent dynamo action (or other sources) plays a vital role even in Babcock-Leighton dynamo models.