Compressional waves propagating in the partially ionised solar lower atmospheric plasmas can easily steepen into nonlinear waves, including shocks. Here we investigate the effect of weak dispersion generated by Hall currents perpendicular to the ambient magnetic field on the characteristics of shock waves. Our study will also focus on the interplay between weak dispersion and partial ionisation of the plasma. Using a multiple scale technique we derive the governing equation in the form of a Korteweg-de Vries-Burgers equation. The effect of weak dispersion on shock waves is obtained using a perturbation technique. The secular behaviour of second order terms is addressed with the help of a renormalisation technique. Our results show that dispersion modifies the characteristics of shock waves and this change is dependent also on the ionisation degree of the plasma. Dispersion can create short lived oscillations in the shocked plasma. The shock fronts become wider with the increase in the number of neutrals in the plasma.

We present a strongly interacting quadruple system associated with the K2 target EPIC 220204960. The K2 target itself is a Kp = 12.7 magnitude star at Teff ~ 6100 K which we designate as "B-N" (blue northerly image). The host of the quadruple system, however, is a Kp = 17 magnitude star with a composite M-star spectrum, which we designate as "R-S" (red southerly image). With a 3.2" separation and similar radial velocities and photometric distances, 'B-N' is likely physically associated with 'R-S', making this a quintuple system, but that is incidental to our main claim of a strongly interacting quadruple system in 'R-S'. The two binaries in 'R-S' have orbital periods of 13.27 d and 14.41 d, respectively, and each has an inclination angle of >89 degrees. From our analysis of radial velocity measurements, and of the photometric lightcurve, we conclude that all four stars are very similar with masses close to 0.4 Msun. Both of the binaries exhibit significant ETVs where those of the primary and secondary eclipses 'diverge' by 0.05 days over the course of the 80-day observations. Via a systematic set of numerical simulations of quadruple systems consisting of two interacting binaries, we conclude that the outer orbital period is very likely to be between 300 and 500 days. If sufficient time is devoted to RV studies of this faint target, the outer orbit should be measurable within a year.

KIC 7177553 was observed by the Kepler satellite to be an eclipsing eccentric binary star system with an 18-day orbital period. Recently, an eclipse timing study of the Kepler binaries has revealed eclipse timing variations in this object with an amplitude of about 100 sec, and an outer period of 529 days. The implied mass of the third body is that of a superJupiter, but below the mass of a brown dwarf. We therefore embarked on a radial velocity study of this binary to determine its system configuration and to check the hypothesis that it hosts a giant planet. From the radial velocity measurements, it became immediately obvious that the same Kepler target contains another eccentric binary, this one with a 16.5-day orbital period. Direct imaging using adaptive optics reveals that the two binaries are separated by 0.4 arcsec (about 167 AU), and have nearly the same magnitude (to within 2%). The close angular proximity of the two binaries, and very similar Gamma velocities, strongly suggest that KIC 7177553 is one of the rare SB4 systems consisting of two eccentric binaries where at least one system is eclipsing. Both systems consist of slowly rotating, non-evolved, solar-like stars of comparable masses. From the orbital separation and the small difference in Gamma velocity, we infer that the period of the outer orbit most likely lies in the range 1000 to 3000 years. New images taken over the next few years, as well as the high-precision astrometry of the Gaia satellite mission, will allow us to set much narrower constraints on the system geometry. Finally, we note that the observed eclipse timing variations in the Kepler data cannot be produced by the second binary. Further spectroscopic observations on a longer time scale will be required to prove the existence of the massive planet.

Babcock-Leighton type solar dynamo models with single-celled meridional circulation are successful in reproducing many solar cycle features. Recent observations and theoretical models of meridional circulation do not indicate a single-celled flow pattern. We examine the role of complex multi-cellular circulation patterns in a Babcock-Leighton solar dynamo in advection- and diffusion-dominated regimes. We show from simulations that presence of a weak, second, high-latitude reverse cell speeds up the cycle and slightly enhances the poleward branch in butterfly diagram, whereas the presence of a second cell in depth reverses the tilt of butterfly wing to an anti-solar type. A butterfly diagram constructed from middle of convection zone yields a solar-like pattern, but this may be difficult to realize in the Sun because of magnetic buoyancy effects. Each of the above cases behaves similarly in higher and lower magnetic diffusivity regimes. However, our dynamo with a meridional circulation containing four cells in latitude behaves distinctly differently in the two regimes, producing solar-like butterfly diagrams with fast cycles in the higher diffusivity regime, and complex branches in butterfly diagrams in the lower diffusivity regime. We also find that dynamo solutions for a four-celled pattern, two in radius and two in latitude, prefer to quickly relax to quadrupolar parity if the bottom flow-speed is strong enough, of similar order of magnitude as the surface flow-speed.

The chemistry of molecular clouds has been studied for decades, with an increasingly general and sophisticated treatment of the reactions involved. Yet the treatment of turbulent diffusion has remained extremely sketchy, assuming simple Fickian diffusion with a scalar diffusivity D. However, turbulent flows similar to those in the interstellar medium are known to give rise to anomalous diffusion phenomena, more specifically superdiffusion (increase of the diffusivity with the spatial scales involved). This paper considers to what extent and in what sense superdiffusion modifies molecular abundances in interstellar clouds. For this first exploration of the subject we employ a very rough treatment of the chemistry and the effect of non-unifom cloud density on the diffusion equation is also treated in a simplified way. The results nevertheless clearly demonstrate that the effect of superdiffusion is quite significant, abundance values at a given radius being modified by order of unity factors. The sense and character of this influence is highly nontrivial.

We present detailed numerical calculations of the fast solar tachocline based on the assumption that the dynamo field dominates over the dynamics of the tachocline. In the present paper of the series, we focus on three shortfalls of the earlier models. First, instead of the simple oscillating dipole poloidal field we study the more general magnetic field structures reminiscent of the butterfly diagram. The migrating field is prescribed as the observed axisymmetric radial magnetic field Stenflo (1988, 1994). Our results are in good agreement with our analitical estimate and our previous works in Forgacs-Dajka & Petrovay (2001,2002), but the polar "dip" in isorotational surfaces is strongly reduced in this case. On the other hand, a more realistic model should have a magnetic diffusivity decreasing significantly inside the radiative interior, so we also explore the effect of diffusivity and magnetic Prandtl number varying with depth. We found that the downwards decreasing magnetic diffusivity and Prandtl number have no significant effect on the solution, although the temporal variation of the tachocline thickness has decreased.

We studied the relation between the distribution of sunspot groups and the Gleissberg cycle. As the magnetic field is related to the area of the sunspot groups, we used area-weighted sunspot group data. On the one hand, we confirm the previously reported long-term cyclic behaviour of the sum of the northern and southern sunspot group mean latitudes, although we found a somewhat longer period (P~104 years). We introduced the difference between the ensemble average area of sunspot groups for the two hemispheres, which turns out to show similar behaviour. We also investigated a further aspect of the Gleissberg cycle where while in the 19th century the consecutive Schwabe cycles are sharply separated from each other, one century later the cycles overlap each other more and more.

One possible scenario for the origin of the solar tachocline, known as the "fast tachocline", assumes that the turbulent diffusivity exceeds eta>10^9 cm^2/s. In this case the dynamics will be governed by the dynamo-generated oscillatory magnetic field on relatively short timescales. Here, for the first time, we present detailed numerical models for the fast solar tachocline with all components of the magnetic field calculated explicitly, assuming axial symmetry and a constant turbulent diffusivity eta and viscosity nu. We find that a sufficiently strong oscillatory poloidal field with dipolar latitude dependence at the tachocline-convective zone boundary is able to confine the tachocline. Exploring the three-dimensional parameter space defined by the viscosity in the range log(nu)=9-11, the magnetic Prandtl number in the range Prm=0.1-10, and the meridional flow amplitude (-3 to +3 cm/s), we also find that the confining field strength B_conf, necessary to reproduce the observed thickness of the tachocline, increases with viscosity nu, with magnetic Prandtl number nu/eta, and with equatorward meridional flow speed. Nevertheless, the resulting B_conf values remain quite reasonable, in the range 10^3-10^4 G, for all parameter combinations considered here. The thickness of the tachocline shows a marked dependence on both time and latitude. A comparison with seismic constraints suggests that best agreement with our models is achieved for the highest values of nu and Prm considered here.

Helioseismic measurements indicate that the solar tachocline is very thin, its full thickness not exceeding 4% of the solar radius. The mechanism that inhibits differential rotation to propagate from the convective zone to deeper into the radiative zone is not known, though several propositions have been made. In this paper we demonstrate by numerical models and analytic estimates that the tachocline can be confined to its observed thickness by a poloidal magnetic field B_p of about one kilogauss, penetrating below the convective zone and oscillating with a period of 22 years, if the tachocline region is turbulent with a diffusivity of eta 10^10 cm^2/s (for a turbulent magnetic Prandtl number of unity). We also show that a similar confinement may be produced for other pairs of the parameter values (B_p, eta). The assumption of the dynamo field penetrating into the tachocline is consistent whenever eta>10^9 cm^2/s.

We present a model for torsional oscillations where the inhibiting effect of active region magnetic fields on turbulence locally reduces turbulent viscous torques, leading to a cycle- and latitude-dependent modulation of the differential rotation. The observed depth dependence of torsional oscillations as well as their phase relationship with the sunspot butterfly diagram are reproduced quite naturally in this model. The resulting oscillation amplitudes are significantly smaller than observed, though they depend rather sensitively on model details. Meridional circulation is found to have only a weak effect on the oscillation pattern.