We present a radiative magneto-hydrodynamic simulation set-up using the PENCIL CODE to study the generation, propagation and dissipation of Alfv\'en waves in the solar atmosphere which includes a convective layer, photosphere below and chromosphere, transition region and the corona above. We prepare a setup of steady-state solar convection where the imposed external magnetic field also has reached the final value gradually starting from a very small value. From that state, we start several simulations by varying the magnetic Prandtl number and the forcing strengths. We find the propagation characteristics of waves excited in this simulation run depend strongly on the magnetic Prandtl number and the wave number of the forcing. For magnetic Prandtl number of unity, we obtain localized heating in the corona due to shock dissipation.

Here, we test the weighted horizontal magnetic gradient ($WG_M$) as a flare precursor, introduced by Korsos et. al. (2015, ApJ, 802, L21), by applying it to a magneto-hydrodynamic (MHD) simulation of solar-like flares Chatterjee et. al. (2016, Physical Review Letters, 116, 10, 101101). The pre-flare evolution of the $WG_M$ and the behavior of the distance parameter between the area-weighted barycenters of opposite polarity sunspots at various heights is investigated in the simulated $\delta$-type sunspot. Four flares emanated from this sunspot. We found the optimum heights above the photosphere where the flare precursors of the $WG_M$ method are identifiable prior to each flare. These optimum heights agree reasonably well with the heights of the occurrence of flares identified from the analysis of their thermal and Ohmic heating signatures in the simulation. We also estimated the expected time of the flare onsets from the duration of the approaching-receding motion of the barycenters of opposite polarities before each single flare. The estimated onset time and the actual time of occurrence of each flare are in good agreement at the corresponding optimum heights. This numerical experiment further supports the use of flare precursors based on the $WG_M$ method.

To explain the large-scale magnetic field of the Sun and other bodies, mean-field dynamo theory is commonly applied where one solves the averaged equations for the mean magnetic field. However, the standard approach breaks down when the scale of the turbulent eddies becomes comparable to the scale of the variations of the mean magnetic field. Models showing sharp magnetic field structures have therefore been regarded as unreliable. Our aim is to look for new effects that occur when we relax the restrictions of the standard approach, which becomes particularly important at the bottom of the convection zone where the size of the turbulent eddies is comparable to the depth of the convection zone itself. We approximate the underlying integro-differential equation by a partial differential equation corresponding to a reaction-diffusion type equation for the mean electromotive force, making an approach that is nonlocal in space and time feasible under conditions where spherical geometry and nonlinearity are included. In agreement with earlier findings, spatio-temporal nonlocality lowers the excitation conditions of the dynamo. Sharp structures are now found to be absent. However, in the surface layers the field remains similar to before.

Active regions (AR) appearing on the surface of the Sun are classified into $\alpha$, $\beta$, $\gamma$, and $\delta$ by the rules of the Mount Wilson Observatory, California on the basis of their topological complexity. Amongst these, the $\delta$-sunspots are known to be super-active and produce the most X-ray flares. Here, we present results from a simulation of the Sun by mimicking the upper layers and the corona, but starting at a more primitive stage than any earlier treatment. We find that this initial state consisting of only a thin sub-photospheric magnetic sheet breaks into multiple flux-tubes which evolve into a colliding-merging system of spots of opposite polarity upon surface emergence, similar to those often seen on the Sun. The simulation goes on to produce many exotic $\delta$-sunspot associated phenomena: repeated flaring in the range of typical solar flare energy release and ejective helical flux ropes with embedded cool-dense plasma filaments resembling solar coronal mass ejections.

The Sun's internal structure and dynamics can be studied with helioseismology, which uses the Sun's natural acoustic oscillations to build up a profile of the solar interior. We discuss how solar acoustic oscillations are affected by the Sun's magnetic field. Careful observations of these effects can be inverted to determine the variations in the structure and dynamics of the Sun's interior as the solar cycle progresses. Observed variations in the structure and dynamics can then be used to inform models of the solar dynamo, which are crucial to our understanding of how the Sun's magnetic field is generated and maintained.

We report the first results of a magnetohydrodynamic (MHD) simulation of the development of a homologous sequence of three coronal mass ejections (CMEs) and demonstrate their so-called cannibalistic behavior. These CMEs originate from the repeated formations and partial eruptions of kink unstable flux ropes as a result of continued emergence of a twisted flux rope across the lower boundary into a pre-existing coronal potential arcade field. The simulation shows that a CME erupting into the open magnetic field created by a preceding CME has a higher speed. The second of the three successive CMEs simulated is cannibalistic catching up and merging with the first into a single fast CME before exiting the domain. All the CMEs including the leading merged CME, attained speeds of about 1000 km/s as they exit the domain. The reformation of a twisted flux rope after each CME eruption during the sustained flux emergence can naturally explain the X-ray observations of repeated reformations of sigmoids and "sigmoid-under-cusp" configurations at a low-coronal source of homologous CMEs.

Evidence for a parity-breaking nature of the magnetic buoyancy instability in a stably stratified gas is reported. In the absence of rotation, no helicity is produced, but the non-helical state is found to be unstable to small helical perturbations during the development of the instability. The parity-breaking nature of an instability in magnetohydrodynamics appears to be the first of its kind and is similar to chiral symmetry breaking in biochemistry. Applications to the production of mean fields in galaxy clusters are being discussed.

(abridged) Context. Turbulent fluxes of angular momentum and heat due to rotationally affected convection play a key role in determining differential rotation of stars. Here we perform a systematic comparison between Cartesian and spherical geometries as a function of the rotation rate. Aims. We extend the earlier studies by using spherical wedges to obtain turbulent angular momentum and heat transport as functions of the rotation rate from stratified convection. We compare results from spherical and Cartesian models in the same parameter regime. In particular, we want to clarify whether the sharp equatorial profile of the horizontal Reynolds stress found in earlier Cartesian models is reproduced in spherical models. Methods. We employ direct numerical simulations of turbulent convection. In order to reach as high spatial resolution as possible in the spherical runs, we model only parts of the latitude and longitude. The rotational influence, measured by the Coriolis number, is varied from zero to roughly seven, which is the regime that is likely to be realised in the solar convection zone. Cartesian simulations are performed in overlapping parameter regimes. Results. For slow rotation we find that the radial and latitudinal turbulent angular momentum fluxes are directed inward and equatorward, respectively. In the rapid rotation regime the radial flux changes sign in contradiction with theory. The latitudinal flux remains mostly equatorward and develops a maximum close to the equator. In Cartesian simulations this peak can be explained by the strong `banana cells'. The latitudinal heat flux is mostly equatorward for slow rotation but changes sign for rapid rotation. The rotation profiles vary from anti-solar (slow equator) for slow and intermediate rotation to solar-like (fast equator) for rapid rotation. The solar-like profiles are dominated by the Taylor--Proudman balance.

A strong toroidal field can exist in form of a magnetic layer in the overshoot region below the solar convection zone. This motivates a more detailed study of the magnetic buoyancy instability with rotation. We calculate the alpha effect due to helical motions caused by a disintegrating magnetic layer in a rotating density-stratified system with angular velocity Omega making an angle theta with the vertical. We also study the dependence of the alpha effect on theta and the strength of the initial magnetic field. We carry out three-dimensional hydromagnetic simulations in Cartesian geometry. A turbulent EMF due to the correlations of the small scale velocity and magnetic field is generated. We use the test-field method to calculate the transport coefficients of the inhomogeneous turbulence produced by the layer. We show that the growth rate of the instability and the twist of the magnetic field vary monotonically with the ratio of thermal conductivity to magnetic diffusivity. The resulting alpha effect is inhomogeneous and increases with the strength of the initial magnetic field. It is thus an example of an "anti-quenched" alpha effect. The alpha effect is nonlocal, requiring around 8--16 Fourier modes to reconstruct the actual EMF based on the actual mean field.

We have calculated the effects of giant convection cells also know as sectoral rolls or banana cells, on p-mode splitting coefficients. We use the technique of quasi-degenerate perturbation theory formulated by Lavely & Ritzwoller in order to estimate the frequency shifts. A possible way of detecting giant cells is to look for even splitting coefficients of 'nearly degenerate' modes in the observational data since these modes have the largest shifts. We find that banana cells having an azimuthal wave number of 16 and maximum vertical velocity of 180 m/s cannot be ruled out from GONG data for even splitting coefficients.

Assuming that the torsional oscillation is driven by the Lorentz force of the magnetic field associated with the sunspot cycle, we use a flux transport dynamo to model it and explain its initiation at a high latitude before the beginning of the sunspot cycle.

The small-scale magnetic helicity produced as a by-product of the large-scale dynamo is believed to play a major role in dynamo saturation. In a mean-field model the generation of small-scale magnetic helicity can be modelled by using the dynamical quenching formalism. Catastrophic quenching refers to a decrease of the saturation field strength with increasing Reynolds number. It has been suggested that catastrophic quenching only affects the region of non-zero helical turbulence (i.e. where the kinematic alpha operates) and that it is possible to alleviate catastrophic quenching by separating the region of strong shear from the alpha layer. We perform a systematic study of a simple axisymmetric two-layer alpha-omega dynamo in a spherical shell for Reynolds numbers in the range 1 < Rm < 10^5. In the framework of dynamical quenching we show that this may not be the case, suggesting that magnetic helicity fluxes would be necessary.

Dynamos in the Sun and other bodies tend to produce magnetic fields that possess magnetic helicity of opposite sign at large and small scales, respectively. The build-up of magnetic helicity at small scales provides an important saturation mechanism. In order to understand the nature of the solar dynamo we need to understand the details of the saturation mechanism in spherical geometry. In particular, we want to understand the effects of magnetic helicity fluxes from turbulence and meridional circulation. We consider a model with just radial shear confined to a thin layer (tachocline) at the bottom of the convection zone. The kinetic alpha owing to helical turbulence is assumed to be localized in a region above the convection zone. The dynamical quenching formalism is used to describe the build-up of mean magnetic helicity in the model, which results in a magnetic alpha effect that feeds back on the kinetic alpha effect. In some cases we compare with results obtained using a simple algebraic alpha quenching formula. In agreement with earlier findings, the magnetic alpha effect in the dynamical alpha quenching formalism has the opposite sign compared with the kinetic alpha effect and leads to a catastrophic decrease of the saturation field strength with increasing magnetic Reynolds numbers. However, at high latitudes this quenching effect can lead to secondary dynamo waves that propagate poleward due to the opposite sign of alpha. Magnetic helicity fluxes both from turbulent mixing and from meridional circulation alleviate catastrophic quenching.

We present nonlinear mean-field alpha-Omega dynamo simulations in spherical geometry with simplified profiles of kinematic alpha effect and shear. We take magnetic helicity evolution into account by solving a dynamical equation for the magnetic alpha effect. This gives a consistent description of the quenching mechanism in mean-field dynamo models. The main goal of this work is to explore the effects of this quenching mechanism in solar-like geometry, and in particular to investigate the role of magnetic helicity fluxes, specifically diffusive and Vishniac-Cho (VC) fluxes, at large magnetic Reynolds numbers (Rm). For models with negative radial shear or positive latitudinal shear, the magnetic alpha effect has predominantly negative (positive) sign in the northern (southern) hemisphere. In the absence of fluxes, we find that the magnetic energy follows an Rm^-1 dependence, as found in previous works. This catastrophic quenching is alleviated in models with diffusive magnetic helicity fluxes resulting in magnetic fields comparable to the equipartition value even for Rm=10^7. On the other hand, models with a shear-driven Vishniac-Cho flux show an increase of the amplitude of the magnetic field with respect to models without fluxes, but only for Rm<10^4. This is mainly a consequence of assuming a vacuum outside the Sun which cannot support a significant VC flux across the boundary. However, in contrast with the diffusive flux, the VC flux modifies the distribution of the magnetic field. In addition, if an ill-determined scaling factor in the expression for the VC flux is large enough, subcritical dynamo action is possible that is driven by the action of shear and the divergence of current helicity flux.

The predictive power of mean-field theory is emphasized by comparing theory with simulations under controlled conditions. The recently developed test-field method is used to extract turbulent transport coefficients both in kinematic as well as nonlinear and quasi-kinematic cases. A striking example of the quasi-kinematic method is provided by magnetic buoyancy-driven flows that produce an alpha effect and turbulent diffusion.

We use direct numerical simulations of forced MHD turbulence with a forcing function that produces two different signs of kinetic helicity in the upper and lower parts of the domain. We show that the mean flux of magnetic helicity from the small-scale field between the two parts of the domain can be described by a Fickian diffusion law with a diffusion coefficient that is approximately independent of the magnetic Reynolds number and about one third of the estimated turbulent magnetic diffusivity. The data suggest that the turbulent diffusive magnetic helicity flux can only be expected to alleviate catastrophic quenching at Reynolds numbers of more than several thousands. We further calculate the magnetic helicity density and its flux in the domain for three different gauges. We consider the Weyl gauge, in which the electrostatic potential vanishes, the pseudo-Lorenz gauge, where the speed of light is replaced by the sound speed, and the `resistive gauge' in which the Laplacian of the magnetic vector potential acts as resistive term. We find that, in the statistically steady state, the time-averaged magnetic helicity density and the magnetic helicity flux are the same in all three gauges.

We have calculated the effects of large scale solar flows like the meridional circulation, giant convection cells and solar rotation on the helioseismic splitting coefficients using quasi-degenerate perturbation theory (QDPT). Our investigation reveals that the effect of poloidal flows like the large scale meridional circulation are difficult to detect in observational data of the global acoustic modes since the frequency shifts are much less than the errors. However, signatures of large scale convective flows may be detected if their amplitude is sufficiently large by looking for frequency shifts due to nearly degenerate modes coupled by convection. In this comprehensive study, we attempt to put limits on the magnitude of flow velocities in giant cells by comparing the splitting coefficients obtained from the QDPT treatment with observational data.

Although the Sun's torsional oscillation is believed to be driven by the Lorentz force associated with the sunspot cycle, this oscillation begins 2--3 yr before the sunspot cycle. We provide a theoretical explanation of this with the help of a solar dynamo model having a meridional circulation penetrating slightly below the bottom of the convection zone, because only in such dynamo models the strong toroidal field forms a few years before the sunspot cycle and at a higher latitude.

Using mean-field models with a dynamical quenching formalism we show that in finite domains magnetic helicity fluxes associated with small-scale magnetic fields are able to alleviate catastrophic quenching. We consider fluxes that result either from advection by a mean flow, the turbulent mixing down the gradient of mean small-scale magnetic helicity concentration, or the explicit removal which may be associated with the effects of coronal mass ejections in the Sun. In the absence of shear, all the small-scale magnetic helicity fluxes are found to be equally strong both for large-scale and small-scale fields. In the presence of shear there is also an additional magnetic helicity flux associated with the mean field, but this flux does not alleviate catastrophic quenching. Outside the dynamo-active region there are neither sources nor sinks of magnetic helicity, so in a steady state this flux must be constant. It is shown that unphysical behavior emerges if the small-scale magnetic helicity flux is forced to vanish within the computational domain.

We test the relation between relative f-mode frequency variation and Lagrangian perturbation in the solar radius obtained by Dziembowski and Goode (2004) using several pairs of solar models and show that it doesn't hold true for any of the model pairs we have used. We attempt to derive a better approximation for the kernel linking the relative frequency changes and the solar radius variation in the subsurface layers

Although systematic measurements of the solar polar magnetic field exist only from mid 1970s, other proxies can be used to infer the polar field at earlier times. The observational data indicate a strong correlation between the polar field at a sunspot minimum and the strength of the next cycle, although the strength of the cycle is not correlated well with the polar field produced at its end. This suggests that the Babcock Leighton mechanism of poloidal field generation from decaying sunspots involves randomness, whereas the other aspects of the dynamo process must be reasonably ordered and deterministic. Only if the magnetic diffusivity within the convection zone is assumed to be high, we can explain the correlation between the polar field at a minimum and the next cycle. We give several independent arguments that the diffusivity must be of this order. In a dynamo model with diffusivity like this, the poloidal field generated at the mid latitudes is advected toward the poles by the meridional circulation and simultaneously diffuses towards the tachocline, where the toroidal field for the next cycle is produced. To model actual solar cycles with a dynamo model having such high diffusivity, we have to feed the observational data of the poloidal field at the minimum into the theoretical model. We develop a method of doing this in a systematic way. Our model predicts that cycle 24 will be a very weak cycle. Hemispheric asymmetry of solar activity is also calculated with our model and compared with observational data.

Whether the upcoming cycle 24 of solar activity will be strong or not is being hotly debated. The solar cycle is produced by a complex dynamo mechanism. We model the last few solar cycles by `feeding' observational data of the Sun's polar magnetic field into our solar dynamo model. Our results fit the observed sunspot numbers of cycles 21-23 extremely well and predict that cycle~24 will be about 35% weaker than cycle~23.

Following an earlier proposal by Choudhuri (2003) for the origin of twist in the magnetic fields of solar active regions, we model the penetration of a wrapped up background poloidal field into a toroidal magnetic flux tube rising through the solar convective zone. The rise of the straight, cylindrical flux tube is followed by numerically solving the induction equation in a comoving Lagrangian frame, while an external poloidal magnetic field is assumed to be radially advected onto the tube with a speed corresponding to the rise velocity. One prediction of our model is the existence of a ring of reverse current helicity on the periphery of active regions. On the other hand, the amplitude of the resulting twist depends sensitively on the assumed structure (diffuse vs. concentrated/intermittent) of the active region magnetic field right before its emergence, and on the assumed vertical profile of the poloidal field. Nevertheless, in the model with the most plausible choice of assumptions a mean twist comparable to the observations can result rather naturally. Our results indicate that the contribution of this mechanism to the twist can be quite significant, and under favourable circumstances it can potentially account for most of the current helicity observed in active regions.

We present here our response to Dikpati et al.'s criticism of our recent solar dynamo model.

We present calculations of helicity based on our solar dynamo model and show that the results are consistent with observational data.

We explore a two-dimensional kinematic solar dynamo model in a full sphere, based on the helioseismically determined solar rotation profile and with an $\alpha$ effect concentrated near the solar surface, which captures the Babcock--Leighton idea that the poloidal field is created from the decay of tilted bipolar active regions. The meridional circulation, assumed to penetrate slightly below the tachocline, plays an important role. Some doubts have recently been raised regarding the ability of such a model to reproduce solar-like dipolar parity. We specifically address the parity issue and show that the dipolar mode is preferred when certain reasonable conditions are satisfied, the most important condition being the requirement that the poloidal field should diffuse efficiently to get coupled across the equator. Our model is shown to reproduce various aspects of observational data, including the phase relation between sunspots and the weak, diffuse field.

In response to the claim by Dziembowski et al. (2001) that the solar radius decreases with magnetic activity at the rate of 1.5 km/yr, we consider the theoretical question whether a radius variation is expected with the solar cycle. If the radius variation is caused by the magnetic pressure of toroidal flux tubes at the bottom of the convection zone, then the dynamo model of Nandy and Choudhuri (2002) would suggest a radius decrease with magnetic activity, in contrast to other dynamo models which would suggest a radius increase. However, the radius decrease is estimated to be only of the order of hundreds of metres.

We examine whether massive binary black holes in spherically symmetric bulges of galaxies can achieve coalescence through the emission of gravitational radiation under the action of stellar dynamical processes alone. In particular, we address the importance of the Brownian motion of a binary's center-of-mass to its continued interaction with stellar orbits that allows it to keep hardening: the restoring force holding the binary at the center decreases as it depletes the central stellar density, causing the amplitude of wandering to increase. We use an analytical model and N-body simulations to calculate the time required to reach the gravitational-radiation dominated stage. We find that a substantial fraction of all massive binaries in galaxies can coalesce within a Hubble time.

We develop a simple physical model to describe the dynamics of a massive point-like object, such as a black hole, near the center of a dense stellar system. It is shown that the total force on this body can be separated into two independent parts, one of which is the slowly varying influence of the aggregate stellar system, and the other being the rapidly fluctuating stochastic force due to discrete encounters with individual stars. For the particular example of a stellar system distributed according to a Plummer model, it is shown that the motion of the black hole is then similar to that of a Brownian particle in a harmonic potential, and we analyze its dynamics using an approach akin to Langevin's solution of the Brownian motion problem. The equations are solved to obtain the average values, time-autocorrelation functions, and probability distributions of the black hole's position and velocity. By comparing these results to N-body simulations, we demonstrate that this model provides a very good statistical description of the actual black hole dynamics. As an application of our model, we use our results to derive a lower limit on the mass of the black hole Sgr A* in the Galactic center.

We derive a model that describes the motion of a Brownian particle in a system which is dominated by gravitational forces. An example of such a system is a massive black hole immersed in a cluster of stars. We compute the dispersion in the position and velocity of such a black hole, and examine whether it achieves a state of equipartition of kinetic energy with the stars. This problem has been considered before only for stellar systems with an isothermal Maxwellian distribution of velocities; here we study other examples and confirm our calculations with N-body simulations. In certain cases, depending on the stellar distribution function, the black hole can acquire a steady state kinetic energy which is very far from equipartition relative to the stars.

We consider the accretion model for anomalous X-ray pulsars proposed recently by Chatterjee, Hernquist and Narayan, in which the emission is powered by accretion from a fossil disk formed by the fallback of material from a supernova explosion. We demonstrate that this model is able to account for the spin period, luminosity and age distributions of the observed population of AXPs for reasonable and broad distributions of the free parameters of the model, namely, the surface magnetic field of the neutron star, the mass of its accretion disk and its initial spin period. In particular, this model is able statistically to account for the puzzlingly narrow observed spin distribution of the AXPs. We show also that if the establishment of fallback accretion disks around isolated neutron stars is a universal phenomenon, then a fairly large minority ($\sim 20%$) of these objects become X-ray bright AXPs or X-ray faint systems spinning down by propeller action, while the rest become radio pulsars.

We present a model for the anomalous X-ray pulsars (AXPs) in which the emission is powered by accretion from a fossil disk, established from matter falling back onto the neutron star following its birth. The time-dependent accretion drives the neutron star towards a ``tracking'' solution in which the rotation period of the star increases slowly, in tandem with the declining accretion rate. For appropriate choices of disk mass, neutron star magnetic field strength and initial spin period, we demonstrate that a rapidly rotating neutron star can be spun down to periods characteristic of AXPs on timescales comparable to the estimated ages of these sources. In other cases, accretion onto the neutron star switches off after a short time, and the star becomes an ordinary radio pulsar. Thus, in our picture, radio pulsars and AXPs are drawn from the same underlying population, in contrast to models involving neutron stars with ultrastrong magnetic fields, which require a new population of stars with very different properties.