Long-term variations of turbulent transport coefficients in a solar-like convective dynamo simulation

The Sun, aside from its eleven year sunspot cycle is additionally subject to long term variation in its activity. In this work we analyse a solar-like convective dynamo simulation, containing approximately 60 magnetic cycles, exhibiting equatorward propagation of the magnetic field, multiple frequencies, and irregular variability, including a missed cycle and complex parity transitions between dipolar and quadrupolar modes. We compute the turbulent transport coefficients, describing the effects of the turbulent velocity field on the mean magnetic field, using the test-field method. The test-field analysis provides a plausible explanation of the missing cycle in terms of the reduction of $\alpha_{\phi\phi}$ in advance of the reduced surface activity, and enhanced downward turbulent pumping during the event to confine some of the magnetic field at the bottom of the convection zone, where local maximum of magnetic energy is observed during the event. At the same time, however, a quenching of the turbulent magnetic diffusivities is observed, albeit differently distributed in depth compared to the other transport coefficients. Therefore, dedicated mean-field modelling is required for verification.Comment: 11 pages, 12 figures, accepted by AN for 14th Potsdam Thinksho

The supernova-regulated ISM. I. The multi-phase structure

We simulate the multi-phase interstellar medium randomly heated and stirred by supernovae, with gravity, differential rotation and other parameters of the solar neighbourhood. Here we describe in detail both numerical and physical aspects of the model, including injection of thermal and kinetic energy by SN explosions, radiative cooling, photoelectric heating and various transport processes. With 3D domain extending 1 kpc^2 horizontally and 2 kpc vertically, the model routinely spans gas number densities 10^-5 - 10^2 cm^-3, temperatures 10-10^8 K, local velocities up to 10^3 km s^-1 (with Mach number up to 25). The thermal structure of the modelled ISM is classified by inspection of the joint probability density of the gas number density and temperature. We confirm that most of the complexity can be captured in terms of just three phases, separated by temperature borderlines at about 10^3 K and 5x10^5 K. The probability distribution of gas density within each phase is approximately lognormal. We clarify the connection between the fractional volume of a phase and its various proxies, and derive an exact relation between the fractional volume and the filling factors defined in terms of the volume and probabilistic averages. These results are discussed in both observational and computational contexts. The correlation scale of the random flows is calculated from the velocity autocorrelation function; it is of order 100 pc and tends to grow with distance from the mid-plane. We use two distinct parameterizations of radiative cooling to show that the multi-phase structure of the gas is robust, as it does not depend significantly on this choice.Comment: 28 pages, 22 figures and 8 table

The supernova-regulated ISM. II. The mean magnetic field

The origin and structure of the magnetic fields in the interstellar medium of spiral galaxies is investigated with 3D, non-ideal, compressible MHD simulations, including stratification in the galactic gravity field, differential rotation and radiative cooling. A rectangular domain, 1x1x2 kpc^{3} in size, spans both sides of the galactic mid-plane. Supernova explosions drive transonic turbulence. A seed magnetic field grows exponentially to reach a statistically steady state within 1.6 Gyr. Following Germano (1992) we use volume averaging with a Gaussian kernel to separate magnetic field into a mean field and fluctuations. Such averaging does not satisfy all Reynolds rules, yet allows a formulation of mean-field theory. The mean field thus obtained varies in both space and time. Growth rates differ for the mean-field and fluctuating field and there is clear scale separation between the two elements, whose integral scales are about 0.7 kpc and 0.3 kpc, respectively.Comment: 5 pages, 10 figures, submitted to Monthly Notices Letter

Backbone Fragility and the Local Search Cost Peak

The local search algorithm WSat is one of the most successful algorithms for solving the satisfiability (SAT) problem. It is notably effective at solving hard Random 3-SAT instances near the so-called `satisfiability threshold', but still shows a peak in search cost near the threshold and large variations in cost over different instances. We make a number of significant contributions to the analysis of WSat on high-cost random instances, using the recently-introduced concept of the backbone of a SAT instance. The backbone is the set of literals which are entailed by an instance. We find that the number of solutions predicts the cost well for small-backbone instances but is much less relevant for the large-backbone instances which appear near the threshold and dominate in the overconstrained region. We show a very strong correlation between search cost and the Hamming distance to the nearest solution early in WSat's search. This pattern leads us to introduce a measure of the backbone fragility of an instance, which indicates how persistent the backbone is as clauses are removed. We propose that high-cost random instances for local search are those with very large backbones which are also backbone-fragile. We suggest that the decay in cost beyond the satisfiability threshold is due to increasing backbone robustness (the opposite of backbone fragility). Our hypothesis makes three correct predictions. First, that the backbone robustness of an instance is negatively correlated with the local search cost when other factors are controlled for. Second, that backbone-minimal instances (which are 3-SAT instances altered so as to be more backbone-fragile) are unusually hard for WSat. Third, that the clauses most often unsatisfied during search are those whose deletion has the most effect on the backbone. In understanding the pathologies of local search methods, we hope to contribute to the development of new and better techniques

Introduction

Non peer reviewe

Oscillating Fracture in Rubber

We have found an oscillating instability of fast-running cracks in thin rubber sheets. A well-defined transition from straight to oscillating cracks occurs as the amount of biaxial strain increases. Measurements of the amplitude and wavelength of the oscillation near the onset of this instability indicate that the instability is a Hopf bifurcation