Showing 1–50 of 60 results for author: Stepanov, R

Numerical simulation of tidal synchronization of the large-scale circulation in Rayleigh-Bénard convection with aspect ratio 1

Authors: S. Röhrborn, P. Jüstel, V. Galindo, F. Stefani, R. Stepanov

Abstract: A possible explanation for the apparent phase stability of the 11.07-year Schwabe cycle of the solar dynamo was the subject of a series of recent papers. The synchronization of the helicity of an instability with azimuthal wavenumber m=1 by a tidal m=2 perturbation played a key role here. To analyze this type of interaction in a paradigmatic set-up, we study a thermally driven Rayleigh-Bénard Conv… ▽ More A possible explanation for the apparent phase stability of the 11.07-year Schwabe cycle of the solar dynamo was the subject of a series of recent papers. The synchronization of the helicity of an instability with azimuthal wavenumber m=1 by a tidal m=2 perturbation played a key role here. To analyze this type of interaction in a paradigmatic set-up, we study a thermally driven Rayleigh-Bénard Convection (RBC) of a liquid metal under the influence of a tide-like electromagnetic forcing. As shown previously, the time-modulation of this forcing emerges as a peak frequency in the m=2 mode of the radial flow velocity component. In this paper we present new numerical results on the interplay between the Large Scale Circulation (LSC) of a RBC flow and the time modulated electromagnetic forcing. △ Less

Submitted 15 September, 2022; originally announced September 2022.

Comments: 6 pages, 7 figures

Journal ref: Magnetohydrodynamics 58 (2022), pp. 285-292

Inverse cascades of kinetic energy and thermal variance in three-dimensional horizontally extended turbulent convection

Authors: Philipp P. Vieweg, Janet D. Scheel, Rodion Stepanov, Jörg Schumacher

Abstract: Inverse cascades of kinetic energy and thermal variance in the subset of vertically homogeneous modes in spectral space are found to cause a slow aggregation to a pair of convective supergranules that eventually fill the whole horizontally extended, three-dimensional, turbulent Rayleigh-Bénard convection layer when a heat flux is prescribed at the top and bottom. An additional weak rotation of the… ▽ More Inverse cascades of kinetic energy and thermal variance in the subset of vertically homogeneous modes in spectral space are found to cause a slow aggregation to a pair of convective supergranules that eventually fill the whole horizontally extended, three-dimensional, turbulent Rayleigh-Bénard convection layer when a heat flux is prescribed at the top and bottom. An additional weak rotation of the layer around the vertical axis stops this aggregation at a scale that is smaller than the lateral domain extension and ceases the inverse cascade for the thermal variance. The inverse cascade for the kinetic energy remains intact, even for times at which the root mean square values of temperature and velocity have reached the statistically steady state. This kinetic energy inverse cascade sustains the horizontally extended convection patterns which are best visible in the temperature field. The resulting characteristic length of the aggregated convection patterns depends on the thermal driving and linearly on the strength of rotation. Our study demonstrates the importance of inverse energy cascades beyond the two-dimensional turbulence case in a three-dimensional convection flow that is subject to a multi-scale energy injection by thermal plumes and driven by boundary heat fluxes as typically present in natural geo- and astrophysical systems, such as solar convection. △ Less

Submitted 20 October, 2022; v1 submitted 25 July, 2022; originally announced July 2022.

Comments: 15 pages, 10 figures

Journal ref: Phys. Rev. Research 4, 043098 (2022)

Synchronizing the helicity of Rayleigh-Bénard convection by a tide-like electromagnetic forcing

Authors: Peter Jüstel, Sebastian Röhrborn, Sven Eckert, Vladimir Galindo, Thomas Gundrum, Rodion Stepanov, Frank Stefani

Abstract: We present results on the synchronization of the helicity in a liquid-metal Rayleigh-Bénard (RB) experiment under the influence of a tide-like electromagnetic forcing with azimuthal wavenumber m=2. We show that for a critical forcing strength the typical Large Scale Circulation (LSC) in the cylindrical vessel of aspect ratio unity is entrained by the period of the tide-like forcing, leading to syn… ▽ More We present results on the synchronization of the helicity in a liquid-metal Rayleigh-Bénard (RB) experiment under the influence of a tide-like electromagnetic forcing with azimuthal wavenumber m=2. We show that for a critical forcing strength the typical Large Scale Circulation (LSC) in the cylindrical vessel of aspect ratio unity is entrained by the period of the tide-like forcing, leading to synchronized helicity oscillations with opposite signs in two half-spaces. The obtained experimental results are consistent with and supported by numerical simulations. A similar entrainment mechanism for the helicity in the solar tachocline may be responsible for the astonishing synchronization of the solar dynamo by the 11.07-year triple synodic alignment cycle of the tidally dominant planets Venus, Earth and Jupiter. △ Less

Submitted 25 July, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Extragalactic Magnetism with SOFIA (Legacy Program) -- I: The magnetic field in the multi-phase interstellar medium of M51

Authors: Alejandro S. Borlaff, Enrique Lopez-Rodriguez, Rainer Beck, Rodion Stepanov, Eva Ntormousi, Annie Hughes, Konstantinos Tassis, Pamela M. Marcum, Lucas Grosset, John E. Beckman, Leslie Proudfit, Susan E. Clark, Tanio Díaz-Santos, Sui Ann Mao, William T. Reach, Julia Roman-Duval, Kandaswamy Subramanian, Le Ngoc Tram, Ellen G. Zweibel, SOFIA Legacy Team

Abstract: The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio… ▽ More The recent availability of high-resolution far-infrared (FIR) polarization observations of galaxies using HAWC+/SOFIA has facilitated studies of extragalactic magnetic fields in the cold and dense molecular disks.We investigate if any significant structural differences are detectable in the kpc-scale magnetic field of the grand design face-on spiral galaxy M51 when traced within the diffuse (radio) and the dense and cold (FIR) interstellar medium (ISM). Our analysis reveals a complex scenario where radio and FIR polarization observations do not necessarily trace the same magnetic field structure. We find that the magnetic field in the arms is wrapped tighter at 154um than at 3 and 6 cm; statistically significant lower values for the magnetic pitch angle are measured at FIR in the outskirts (R > 7 kpc) of the galaxy. This difference is not detected in the interarm region. We find strong correlations of the polarization fraction and total intensity at FIR and radio with the gas column density and 12CO(1-0) velocity dispersion. We conclude that the arms show a relative increase of small-scale turbulent B-fields at regions with increasing column density and dispersion velocities of the molecular gas. No correlations are found with HI neutral gas. The star formation rate shows a clear correlation with the radio polarized intensity, which is not found in FIR, pointing to a small-scale dynamo-driven B-field amplification scenario. This work shows that multi-wavelength polarization observations are key to disentangling the interlocked relation between star formation, magnetic fields, and gas kinematics in the multi-phase ISM. △ Less

Submitted 8 September, 2021; v1 submitted 19 May, 2021; originally announced May 2021.

Comments: 35 pages, 25 figures, 5 tables. Accepted for publication in ApJ

Generation of zonal flows in convective systems by travelling thermal waves

Authors: Philipp Reiter, Xuan Zhang, Rodion Stepanov, Olga Shishkina

Abstract: This work addresses the effect of travelling thermal waves applied at the fluid layer surface, on the formation of global flow structures in 2D and 3D convective systems. For a broad range of Rayleigh numbers ($10^3\leq Ra \leq 10^7$) and thermal wave frequencies ($10^{-4}\leq Ω\leq 10^{0}$), we investigate flows with and without imposed mean temperature gradients. Our results confirm that the tra… ▽ More This work addresses the effect of travelling thermal waves applied at the fluid layer surface, on the formation of global flow structures in 2D and 3D convective systems. For a broad range of Rayleigh numbers ($10^3\leq Ra \leq 10^7$) and thermal wave frequencies ($10^{-4}\leq Ω\leq 10^{0}$), we investigate flows with and without imposed mean temperature gradients. Our results confirm that the travelling thermal waves can cause zonal flows, i.e. strong mean horizontal flows. We show that the zonal flows in diffusion dominated regimes are driven purely by the Reynolds stresses, always travelling retrograde, while in convection dominated regimes, mean flow advection, caused by tilted convection cells, becomes dominant, which generally leads to prograde mean zonal flows. By means of direct numerical simulations we validate theoretical predictions made for the diffusion dominated regime. Furthermore, we make use of the linear stability analysis and explain the existence of the tilted convection cell mode. Our extensive 3D simulations support the results for 2D flows and thus confirm the relevance of the findings for geopyhsical and astrophysical systems. △ Less

Submitted 20 December, 2020; v1 submitted 3 September, 2020; originally announced September 2020.

Comments: 24 pages, 13 figures, submitted to the journal of fluid mechanics (jfm)

Shaken and stirred: When Bond meets Suess-de Vries and Gnevyshev-Ohl

Authors: Frank Stefani, Rodion Stepanov, Tom Weier

Abstract: We argue that the most prominent temporal features of the solar dynamo, in particular the Hale cycle, the Suess-de Vries cycle (associated with variations of the Gnevyshev-Ohl rule), Gleissberg-type cycles, and grand minima can be self-consistently explained by double synchronization with the 11.07-years periodic tidal forcing of the Venus-Earth-Jupiter system and the (mainly) 19.86-years periodic… ▽ More We argue that the most prominent temporal features of the solar dynamo, in particular the Hale cycle, the Suess-de Vries cycle (associated with variations of the Gnevyshev-Ohl rule), Gleissberg-type cycles, and grand minima can be self-consistently explained by double synchronization with the 11.07-years periodic tidal forcing of the Venus-Earth-Jupiter system and the (mainly) 19.86-years periodic motion of the Sun around the barycenter of the solar system. In our numerical simulation, grand minima, and clusters thereof, emerge as intermittent and non-periodic events on millennial time scales, very similar to the series of Bond events which were observed throughout the Holocene and the last glacial period. If confirmed, such an intermittent transition to chaos would prevent any long-term prediction of solar activity, notwithstanding the fact that the shorter-term Hale and Suess-de Vries cycles are clocked by planetary motion. △ Less

Submitted 19 June, 2020; v1 submitted 15 June, 2020; originally announced June 2020.

Comments: 25 pages, 13 figures

Generating a tide-like flow in a cylindrical vessel by electromagnetic forcing

Authors: Peter Jüstel, Sebastian Röhrborn, Peter Frick, Vladimir Galindo, Thomas Gundrum, Felix Schindler, Frank Stefani, Rodion Stepanov, Tobias Vogt

Abstract: We show and compare numerical and experimental results on the electromagnetic generation of a tide-like flow structure in a cylindrical vessel which is filled with the eutectic liquid metal alloy GaInSn. Fields of various strengths and frequencies are applied to drive liquid metal flows. The impact of the field variations on amplitude and structure of the flows is investigated. The results represe… ▽ More We show and compare numerical and experimental results on the electromagnetic generation of a tide-like flow structure in a cylindrical vessel which is filled with the eutectic liquid metal alloy GaInSn. Fields of various strengths and frequencies are applied to drive liquid metal flows. The impact of the field variations on amplitude and structure of the flows is investigated. The results represent the basis for a future Rayleigh-Bénard experiment, in which a modulated tide-like flow perturbation is expected to synchronize the typical sloshing mode of the large-scale circulation. A similar entrainment mechanism for the helicity in the Sun may be responsible for the synchronization of the solar dynamo with the alignment cycle of the tidally dominant planets Venus, Earth and Jupiter. △ Less

Submitted 5 June, 2020; originally announced June 2020.

Comments: 11 pages, 19 figure

Wavelet analysis of the long-term activity of V833 Tau

Authors: R. Stepanov, N. I. Bondar', M. M. Katsova, D. Sokoloff, P. Frick

Abstract: The bulk of available stellar activity observations is frequently checked for the manifestation of signs in comparison with the known characteristic of solar magnetic modulation. The problem is that stellar activity records are usually an order of magnitude shorter than available observations of solar activity variation. Therefore, the resolved time scales of stellar activity are insufficient to d… ▽ More The bulk of available stellar activity observations is frequently checked for the manifestation of signs in comparison with the known characteristic of solar magnetic modulation. The problem is that stellar activity records are usually an order of magnitude shorter than available observations of solar activity variation. Therefore, the resolved time scales of stellar activity are insufficient to decide reliably that a cyclic variation for a particular star is similar to the well-known 11-yr sunspot cycles. As a result, recent studies report several stars with double or multiple cycles which serve to challenge the underlying theoretical understanding. This is why a consistent method to separate 'true' cycles from stochastic variations is required. In this paper, we suggest that a conservative method, based on the best practice of wavelet analysis previously applied to the study of solar activity, for studying and interpreting the longest available stellar activity record - photometric monitoring of V833 Tau for more than 100 years. We find that the observed variations of V833 Tau with timescales of 2-50 yr should be comparable with the known quasi-periodic solar mid-term variations, whereas the true cycle of V833 Tau, if it exists, should be of about a century or even longer. We argue that this conclusion does not contradict the expectations from stellar dynamo theory. △ Less

Submitted 22 May, 2020; originally announced May 2020.

Comments: accepted in MNRAS

Inverse cascade of energy in helical turbulence

Authors: Franck Plunian, Andrei Teimurazov, Rodion Stepanov, Mahendra Kumar Verma

Abstract: Using direct numerical simulation of hydrodynamic turbulence with helicity forcing applied at all scales, a near-maximum helical turbulent state is obtained, with an inverse energy cascade at scales larger than the energy forcing scale and a forward helicity cascade at scales smaller than the energy forcing scale. In contrast to previous studies using decimated triads, our simulations contain all… ▽ More Using direct numerical simulation of hydrodynamic turbulence with helicity forcing applied at all scales, a near-maximum helical turbulent state is obtained, with an inverse energy cascade at scales larger than the energy forcing scale and a forward helicity cascade at scales smaller than the energy forcing scale. In contrast to previous studies using decimated triads, our simulations contain all possible triads. By computing the shell-to-shell energy fluxes, we show that the inverse energy cascade results from weakly non-local interactions among homochiral triads. Varying the helicity injection range of scales leads to necessary conditions to obtain an inverse energy cascade. △ Less

Submitted 10 May, 2020; originally announced May 2020.

Comments: 13 pages, 10 figures

Journal ref: J. Fluid Mech. 895 (2020) A13

Enstrophy transfers in helical turbulence

Authors: Shubhadeep Sadhukhan, Roshan Samuel, Franck Plunian, Rodion Stepanov, Ravi Samtaney, Mahendra Kumar Verma

Abstract: In this paper we study the enstrophy transers in helical turbulence using direct numerical simulation. We observe that the helicity injection does not have significant effects on the inertial-range energy and helicity spectra ($\sim k^{-5/3}$) and fluxes (constants). We also calculate the separate contributions to enstrophy transfers via velocity to vorticity and vorticity to vorticity channels. T… ▽ More In this paper we study the enstrophy transers in helical turbulence using direct numerical simulation. We observe that the helicity injection does not have significant effects on the inertial-range energy and helicity spectra ($\sim k^{-5/3}$) and fluxes (constants). We also calculate the separate contributions to enstrophy transfers via velocity to vorticity and vorticity to vorticity channels. There are four different enstrophy fluxes associated with the former channel or vorticity stretching, and one flux associated with the latter channel or vorticity advection. In the inertial range, the fluxes due to vorticity stretching are larger than that due to advection. These transfers too are insensitive to helicity injection. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: 7 pages, 4 figures

Journal ref: Physical Review Fluids, 4, 084607 (2019)

On uniqueness of transfer rates in magnetohydrodynamic turbulence

Authors: Franck Plunian, Rodion Stepanov, Mahendra Kumar Verma

Abstract: In hydrodynamic and magnetohydrodynamic turbulence, expressions for the transfer rates rely on integrals over wavenumber triads $(\textbf{k,p,q})$ satisfying $\textbf{k+p+q=0}$. As an example $S_E^{uu}(\textbf{k|p,q})$ denotes the kinetic energy transfer rate to the mode $\textbf{k}$, from the two other modes in the triad, $\textbf{p}$ and $\textbf{q}$. As noted by Kraichnan (1958), in… ▽ More In hydrodynamic and magnetohydrodynamic turbulence, expressions for the transfer rates rely on integrals over wavenumber triads $(\textbf{k,p,q})$ satisfying $\textbf{k+p+q=0}$. As an example $S_E^{uu}(\textbf{k|p,q})$ denotes the kinetic energy transfer rate to the mode $\textbf{k}$, from the two other modes in the triad, $\textbf{p}$ and $\textbf{q}$. As noted by Kraichnan (1958), in $S_E^{uu}(\textbf{k|p,q})$, what fraction of the energy transferred to the mode $\textbf{k}$ originated from $\textbf{p}$ and which from $\textbf{q}$ is unknown . Such an expression is thus incongruent with the customary description of turbulence in terms of two-scale energy exchange. Notwithstanding this issue, Dar etal. (2001) further decomposed these transfers into separate contributions from $\textbf{p}$-to-$\textbf{k}$ and $\textbf{q}$-to-$\textbf{k}$, thus introducing the concept of mode-to-mode transfers that they applied to MHD turbulence. Doing so, they had to set aside additional transfers circulating within each triad, but failed to calculate them. In the present paper we explain how to derive the complete expressions of the mode-to-mode transfers, including the circulating transfers. We do it for kinetic energy and kinetic helicity in hydrodynamic turbulence, for kinetic energy, magnetic energy and magnetic helicity in MHD turbulence. Separating the contribution of magnetic advection from magnetic stretching, the energy mode-to-mode transfer rates involving the magnetic field become uniquely defined, in striking contrast to the hydrodynamic case. The magnetic helicity mode-to-mode transfer rate is also found to be uniquely defined, contrary to kinetic helicity in hydrodynamics. We find that shell-to-shell transfer rates have the same properties as mode-to-mode transfer rates. Finally calculating the fluxes, we show how the circulating transfers cancel in accordance with conservation laws. △ Less

Submitted 21 April, 2020; originally announced April 2020.

Comments: 19 pages, 3 figures

Journal ref: Journal of Plasma Physics, Volume 85, Issue 5 905850507 (2019)

Phase coherence and phase jumps in the Schwabe cycle

Authors: F. Stefani, J. Beer, A. Giesecke, T. Gloaguen, M. Seilmayer, R. Stepanov, T. Weier

Abstract: Guided by the working hypothesis that the Schwabe cycle of solar activity is synchronized by the 11.07 years alignment cycle of the tidally dominant planets Venus, Earth and Jupiter, we reconsider the phase diagrams of sediment accumulation rates in Lake Holzmaar, and of methanesulfonate (MSA) data in the Greenland ice core GISP2, which are available for the period 10000-9000 cal. BP. Since some h… ▽ More Guided by the working hypothesis that the Schwabe cycle of solar activity is synchronized by the 11.07 years alignment cycle of the tidally dominant planets Venus, Earth and Jupiter, we reconsider the phase diagrams of sediment accumulation rates in Lake Holzmaar, and of methanesulfonate (MSA) data in the Greenland ice core GISP2, which are available for the period 10000-9000 cal. BP. Since some half-cycle phase jumps appearing in the output signals are, very likely, artifacts of applying a biologically substantiated transfer function, the underlying solar input signal with a dominant 11.04 years periodicity can be considered as mainly phase-coherent over the 1000 years period in the early Holocene. For more recent times, we show that the re-introduction of a hypothesized "lost cycle" at the beginning of the Dalton minimum would lead to a real phase jump. Similarly, by analyzing various series of $^{14}$C and $^{10}$Be data and comparing them with Schove's historical cycle maxima, we support the existence of another "lost cycle" around 1565, also connected with a real phase jump. Viewed synoptically, our results lend greater plausibility to the starting hypothesis of a tidally synchronized solar cycle, which at times can undergo phase jumps, although the competing explanation in terms of a non-linear solar dynamo with increased coherence cannot be completely ruled out. △ Less

Submitted 10 July, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: 14 pages, 16 figures; to be published in Astronomische Nachrichten

Spectral characteristic of mid-term quasi-periodicities in sunspots data

Authors: P. Frick, D. Sokoloff, R. Stepanov, V. Pipin, I. Usoskin

Abstract: Numerous analyses suggest the existence of various quasi-periodicities in solar activity. The power spectrum of solar activity recorded in sunspot data is dominated by the $\sim$11-year quasi-periodicity, known as the Schwabe cycle. In the mid-term range (1 month -- 11 years) a pronounced variability known as a quasi-biennial oscillation (QBO) is widely discussed. In the shorter time scale a prono… ▽ More Numerous analyses suggest the existence of various quasi-periodicities in solar activity. The power spectrum of solar activity recorded in sunspot data is dominated by the $\sim$11-year quasi-periodicity, known as the Schwabe cycle. In the mid-term range (1 month -- 11 years) a pronounced variability known as a quasi-biennial oscillation (QBO) is widely discussed. In the shorter time scale a pronounced peak, corresponding to the synodic solar rotation period ($\sim$ 27 days) is observed. Here we revisited the mid-term solar variability in terms of statistical dynamic of fully turbulent systems, where solid arguments are required to accept an isolated dominant frequency in a continuous (smooth) spectrum. For that, we first undertook an unbiased analysis of the standard solar data, sunspot numbers and the F10.7 solar radioflux index, by applying a wavelet tool, which allows one to perform a frequency-time analysis of the signal. Considering the spectral dynamics of solar activity cycle by cycle, we showed that no single periodicity can be separated, in a statistically significant manner, in the specified range of periods. We examine whether a model of solar dynamo can reproduce the mid-term oscillation pattern observed in solar data. We found that a realistically observed spectrum can be explained if small spatial (but not temporal) scales are effectively smoothed. This result is important because solar activity is a it global feature, although monitored via small-scale tracers like sunspots. △ Less

Submitted 15 November, 2019; originally announced November 2019.

Comments: accepted in MNRAS

Schwabe, Gleissberg, Suess-de Vries: Towards a consistent model of planetary synchronization of solar cycles

Authors: F. Stefani, A. Giesecke, M. Seilmayer, R. Stepanov, T. Weier

Abstract: Aiming at a consistent planetary synchronization model of both short-term and long-term solar cycles, we start with an analysis of Schove's historical data of cycle maxima. Their deviations (residuals) from the average cycle duration of 11.07 years show a high degree of regularity, comprising a dominant 200-year period (Suess-de Vries cycle), and a few periods around 100 years (Gleissberg cycle).… ▽ More Aiming at a consistent planetary synchronization model of both short-term and long-term solar cycles, we start with an analysis of Schove's historical data of cycle maxima. Their deviations (residuals) from the average cycle duration of 11.07 years show a high degree of regularity, comprising a dominant 200-year period (Suess-de Vries cycle), and a few periods around 100 years (Gleissberg cycle). Encouraged by their robustness, we support previous forecasts of an upcoming grand minimum in the 21st century. To explain the long-term cycles, we enhance our tidally synchronized solar dynamo model by a modulation of the field storage capacity of the tachocline with the orbital angular momentum of the Sun, which is dominated by the 19.86-year periodicity of the Jupiter-Saturn synodes. This modulation of the 22.14 years Hale cycle leads to a 193-year beat period of dynamo activity which is indeed close to the Suess-de Vries cycle. For stronger dynamo modulation, the model produces additional peaks at typical Gleissberg frequencies, which seem to be explainable by the non-linearities of the basic beat process, leading to a bi-modality of the Schwabe cycle. However, a complementary role of beat periods between the Schwabe cycle and the Jupiter-Uranus/Neptune synodic cycles cannot be completely excluded. △ Less

Submitted 30 September, 2020; v1 submitted 23 October, 2019; originally announced October 2019.

Comments: 11 pages, 4 figures, final version as published

Journal ref: Magnetohydrodynamics, Vol. 56, No. 2/3, pp. 269-280, 2020

Cross helicity sign reversals in the dissipative scales of magnetohydrodynamic turbulence

Authors: V. Titov, R. Stepanov, N. Yokoi, M. Verma, R. Samtaney

Abstract: We perform direct numerical simulations of magnetohydrodynamic (MHD) turbulence with kinetic energy and cross helicity injections at large scales. We observe that the cross helicity changes sign as we go from large and intermediate scales to small scales. In addition, the magnetic reconnections are strongest at the regions where the cross helicity changes sign and becomes smallest in magnitude. Th… ▽ More We perform direct numerical simulations of magnetohydrodynamic (MHD) turbulence with kinetic energy and cross helicity injections at large scales. We observe that the cross helicity changes sign as we go from large and intermediate scales to small scales. In addition, the magnetic reconnections are strongest at the regions where the cross helicity changes sign and becomes smallest in magnitude. Thus, our simulations provide an important window to explore the regions of magnetic reconnections in nonlinear MHD. △ Less

Submitted 21 February, 2019; originally announced February 2019.

Comments: appear in MAGNETOHYDRODYNAMICS

Energy transfers in MHD turbulence and its applications to dynamo

Authors: Mahendra K. Verma, Rodion Stepanov, Franck Plunian

Abstract: In this paper we describe the mode-to-mode energy transfers and energy fluxes of MHD turbulence. These energy transfers are very useful for understanding the turbulence dynamics, as well as for applications such as dynamo. We illustrate how the energy fluxes provide valuable insights into the mechanism of growth of the large-scale magnetic energy in dynamo. In this paper we describe the mode-to-mode energy transfers and energy fluxes of MHD turbulence. These energy transfers are very useful for understanding the turbulence dynamics, as well as for applications such as dynamo. We illustrate how the energy fluxes provide valuable insights into the mechanism of growth of the large-scale magnetic energy in dynamo. △ Less

Submitted 18 February, 2019; originally announced February 2019.

Comments: to appear in Magnetohydrodynamics

Measuring the filamentary structure of interstellar clouds through wavelets

Authors: Volker Ossenkopf-Okada, Rodion Stepanov

Abstract: The ubiquitous presence of filamentary structures in the interstellar medium asks for an unbiased characterization of their properties including a stability analysis. We propose a novel technique to measure the spectrum of filaments in any two-dimensional data set. Using anisotropic wavelets we can quantify and distinguish local and global anisotropies and measure the size distribution of filament… ▽ More The ubiquitous presence of filamentary structures in the interstellar medium asks for an unbiased characterization of their properties including a stability analysis. We propose a novel technique to measure the spectrum of filaments in any two-dimensional data set. Using anisotropic wavelets we can quantify and distinguish local and global anisotropies and measure the size distribution of filaments. The wavelet analysis does not need any assumptions on the alignment or shape of filaments in the maps, but directly measures their typical spatial dimensions. In a rigorous test program, we calibrate the scale-dependence of the method and test the angular and spatial sensitivity. We apply the method to molecular line maps from magneto-hydrodynamic (MHD) simulations and observed column density maps from Herschel observations. When applying the anisotropic wavelet analysis to the MHD data, we find that the observed filament sizes depend on the combination of magnetic-field dominated density-velocity correlations with radiative transfer effects. This can be exploited by observing tracers with different optical depth to measure the transition from a globally ordered large-scale structure to small-scale filaments with entangled field lines. The unbiased view to Herschel column density maps does not confirm a universal characteristic filament width. The map of the Polaris Flare shows an almost scale-free filamentary spectrum up to the size of the dominating filament of about 0.4pc. For the Aquila molecular cloud the range of filament widths is limited to 0.05-0.2pc. The filaments in Polaris show no preferential direction in contrast to the global alignment that we trace in Aquila. By comparing the power in isotropic and anisotropic structures we can measure the relative importance of spherical and cylindrical collapse modes and their spatial distribution. △ Less

Submitted 5 November, 2018; originally announced November 2018.

Comments: Accepted for publication by A&A

Journal ref: A&A 621, A5 (2019)

Magnetic arms of NGC6946 traced in the Faraday cubes at low radio frequencies

Authors: Anton Chupin, Rainer Beck, Peter Frick, George Heald, Dmitry Sokoloff, Rodion Stepanov

Abstract: Magnetic fields in galaxies exist on various spatial scales. Large-scale magnetic fields are thought to be generated by the $α-Ω$ dynamo. Small-scale galactic magnetic fields (1 kpc and below) can be generated by tangling the large-scale field or by the small-scale turbulent dynamo. The analysis of field structures with the help of polarized radio continuum emission is hampered by the effect of Fa… ▽ More Magnetic fields in galaxies exist on various spatial scales. Large-scale magnetic fields are thought to be generated by the $α-Ω$ dynamo. Small-scale galactic magnetic fields (1 kpc and below) can be generated by tangling the large-scale field or by the small-scale turbulent dynamo. The analysis of field structures with the help of polarized radio continuum emission is hampered by the effect of Faraday dispersion (due to fluctuations in magnetic field and/or thermal electron density) that shifts signals from large to small scales. At long observation wavelengths large-scale magnetic fields may become invisible, as in the case of spectro-polarimetric data cube of the spiral galaxy NGC~6946 observed with the Westerbork Radio Synthesis Telescope in the wavelength range 17-23 cm. The application of RM Synthesis alone does not overcome this problem. We propose to decompose the Faraday data cube into data cubes at different spatial scales by a wavelet transform. Signatures of the `magnetic arms' observed in NGC~6946 at shorter wavelengths become visible. Our method allows us to search for large-scale field patterns in data cubes at long wavelengths, as provided by new-generation radio telescopes. △ Less

Submitted 16 August, 2018; originally announced August 2018.

Comments: 7 pages, 7 figures, 1 table

Transient flows and reorientations of large-scale convection in a cubic cell

Authors: Andrei Vasiliev, Peter Frick, Abhishek Kumar, Rodion Stepanov, Andrei Sukhanovskii, Mahendra K. Verma

Abstract: The transient processes of a turbulent large-scale convective circulation (LSC) in a cubic cell are investigated using large-eddy simulations for Rayleigh number $\Ray=10^8$ and Prandtl number $\Pran=0.7$. For the first time, we have explicitly shown that LSC is accompanied by large-scale azimuthal flows with non-zero total angular momentum. It is also shown that solid-body rotation of the entire… ▽ More The transient processes of a turbulent large-scale convective circulation (LSC) in a cubic cell are investigated using large-eddy simulations for Rayleigh number $\Ray=10^8$ and Prandtl number $\Pran=0.7$. For the first time, we have explicitly shown that LSC is accompanied by large-scale azimuthal flows with non-zero total angular momentum. It is also shown that solid-body rotation of the entire fluid is not realized. It is found that correlation between rotation of LSC plane and the mean azimuthal motion is high during quasiperiodic oscillations of LSC near the diagonal plane and relatively weak during LSC reorientations. We propose a new plausible scenario for the reorientations of the LSC in a cube that does not involve a mean azimuthal flow. Instead of a single-roll, we introduce the superposition of a pair of large-scale orthogonal quasi-two-dimensional (Q2D) rolls and the reorientation of the LSC occurs as a result of the cessation of one of the Q2D rolls. This scenario is consistent with all known experimental and numerical data. △ Less

Submitted 27 July, 2019; v1 submitted 17 May, 2018; originally announced May 2018.

Comments: submitted to International Communications in Heat and Mass Transfer

Energy Spectra and Fluxes in Dissipation Range of Turbulent and Laminar Flows

Authors: Mahendra K. Verma, Abhishek Kumar, Praveen Kumar, Satyajit Barman, Anando G. Chatterjee, Ravi Samtaney, Rodion Stepanov

Abstract: Two well-known turbulence models that describe the energy spectrum in the inertial and dissipative ranges simultaneously are by Pao~(1965) and Pope~(2000). In this paper, we compute the energy spectrum $E(k)$ and energy flux $Π(k)$ using direct numerical simulations on grids up to $4096^3$, and show consistency between the numerical results and the predictions by the aforementioned models. We also… ▽ More Two well-known turbulence models that describe the energy spectrum in the inertial and dissipative ranges simultaneously are by Pao~(1965) and Pope~(2000). In this paper, we compute the energy spectrum $E(k)$ and energy flux $Π(k)$ using direct numerical simulations on grids up to $4096^3$, and show consistency between the numerical results and the predictions by the aforementioned models. We also construct a model for laminar flows that predicts $E(k)\sim k^{-1} \exp(-k)$ and $Π(k)\sim k \exp(-k)$. Our model predictions match with the numerical results. We emphasize differences on the energy transfers in the two flows---they are {\em local} in the turbulent flows, and {\em nonlocal} in laminar flows. △ Less

Submitted 18 February, 2019; v1 submitted 14 May, 2017; originally announced May 2017.

Journal ref: Fluid Dynamics, 2018, Vol. 53, No. 6, pp. 862-873

Magnetic and gaseous spiral arms in M83

Authors: P. Frick, R. Stepanov, R. Beck, D. Sokoloff, A. Shukurov, M. Ehle, A. Lundgren

Abstract: Isotropic and anisotropic wavelet transforms are used to decompose the images of the spiral galaxy M83 in various tracers to quantify structures in a range of scales from 0.2 to 10 kpc. We used radio polarization observations at λ6 cm and 13 cm obtained with the VLA, Effelsberg and ATCA telescopes and APEX sub-mm observations at 870 μm, which are first published here, together with maps of the emi… ▽ More Isotropic and anisotropic wavelet transforms are used to decompose the images of the spiral galaxy M83 in various tracers to quantify structures in a range of scales from 0.2 to 10 kpc. We used radio polarization observations at λ6 cm and 13 cm obtained with the VLA, Effelsberg and ATCA telescopes and APEX sub-mm observations at 870 μm, which are first published here, together with maps of the emission of warm dust, ionized gas, molecular gas, and atomic gas. The spatial power spectra are similar for the tracers of dust, gas, and total magnetic field, while the spectra of the ordered magnetic field are significantly different. The wavelet cross-correlation between all material tracers and total magnetic field is high, while the structures of the ordered magnetic field are poorly correlated with those of other tracers. -- The magnetic field configuration in M83 contains pronounced magnetic arms. Some of them are displaced from the corresponding material arms, while others overlap with the material arms. The magnetic field vectors at λ6 cm are aligned with the outer material arms, while significant deviations occur in the inner arms and in the bar region, possibly due to non-axisymmetric gas flows. Outside the bar region, the typical pitch angles of the material and magnetic spiral arms are very close to each other at about 10°. The typical pitch angle of the magnetic field vectors is about 20° larger than that of the material spiral arms. One of the main magnetic arms in M83 is displaced from the gaseous arms, while the other main arm overlaps a gaseous arm. We propose that a regular spiral magnetic field generated by a mean-field dynamo is compressed in material arms and partly aligned with them. The interaction of galactic dynamo action with a transient spiral pattern is a promising mechanism for producing such complicated spiral patterns as in M83. △ Less

Submitted 19 November, 2015; v1 submitted 2 October, 2015; originally announced October 2015.

Comments: Accepted for publication in Astronomy & Astrophysics

MSC Class: 85-05

Journal ref: A&A 585, A21 (2016)

Current helicity and magnetic field anisotropy in solar active regions

Authors: H. Xu, R. Stepanov, K. Kuzanyan, D. Sokoloff, H. Zhang, Y. Gao

Abstract: The electric current helicity density $\displaystyle χ=\langleε_{ijk}b_i\frac{\partial b_k}{\partial x_j}\rangle$ contains six terms, where $b_i$ are components of the magnetic field. Due to the observational limitations, only four of the above six terms can be inferred from solar photospheric vector magnetograms. By comparing the results for simulation we distinguished the statistical difference… ▽ More The electric current helicity density $\displaystyle χ=\langleε_{ijk}b_i\frac{\partial b_k}{\partial x_j}\rangle$ contains six terms, where $b_i$ are components of the magnetic field. Due to the observational limitations, only four of the above six terms can be inferred from solar photospheric vector magnetograms. By comparing the results for simulation we distinguished the statistical difference of above six terms for isotropic and anisotropic cases. We estimated the relative degree of anisotropy for three typical active regions and found that it is of order 0.8 which means the assumption of local isotropy for the observable current helicity density terms is generally not satisfied for solar active regions. Upon studies of the statistical properties of the anisotropy of magnetic field of solar active regions with latitudes and with evolution in the solar cycle, we conclude that the consistency of that assumption of local homogeneity and isotropy requires further analysis in the light of our findings. △ Less

Submitted 9 September, 2015; originally announced September 2015.

Comments: MNRAS submitted

Non-Kolmogorov cascade of helicity driven turbulence

Authors: Mouloud Kessar, Franck Plunian, Rodion Stepanov, Guillaume Balarac

Abstract: We solve the Navier-Stokes equations with two simultaneous forcings. One forcing is applied at a given large-scale and it injects energy. The other forcing is applied at all scales belonging to the inertial range and it injects helicity. In this way we can vary the degree of turbulence helicity from non helical to maximally helical. We find that increasing the rate of helicity injection does not c… ▽ More We solve the Navier-Stokes equations with two simultaneous forcings. One forcing is applied at a given large-scale and it injects energy. The other forcing is applied at all scales belonging to the inertial range and it injects helicity. In this way we can vary the degree of turbulence helicity from non helical to maximally helical. We find that increasing the rate of helicity injection does not change the energy flux. On the other hand the level of total energy is strongly increased and the energy spectrum gets steeper. The energy spectrum spans from a Kolmogorov scaling law $k^{-5/3}$ for a non-helical turbulence, to a non-Kolmogorov scaling law $k^{-7/3}$ for a maximally helical turbulence. In the later case we find that the characteristic time of the turbulence is not the turnover time but a time based on the helicity injection rate. We also analyse the results in terms of helical modes decomposition. For a maximally helical turbulence one type of helical mode is found to be much more energetic than the other one, by several orders of magnitude. The energy cascade of the most energetic type of helical mode results from the sum of two fluxes. One flux is negative and can be understood in terms of a decimated model. This negative flux is however not sufficient to lead an inverse energy cascade. Indeed the other flux involving the least energetic type of helical mode is positive and the largest. The least energetic type of helical mode is then essential and cannot be neglected. △ Less

Submitted 9 September, 2015; originally announced September 2015.

Hindered energy cascade in highly helical isotropic turbulence

Authors: Rodion Stepanov, Ephim Golbraikh, Peter Frick, Alexander Shestakov

Abstract: The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade to sma… ▽ More The conventional approach to the turbulent energy cascade, based on Richardson-Kolmogorov phenomenology, ignores the topology of emerging vortices, which is related to the helicity of the turbulent flow. It is generally believed that helicity can play a significant role in turbulent systems, e.g., supporting the generation of large-scale magnetic fields, but its impact on the energy cascade to small scales has never been observed. We suggest for the first time a generalized phenomenology for isotropic turbulence with an arbitrary spectral distribution of the helicity. We discuss various scenarios of direct turbulent cascades with new helicity effect, which can be interpreted as a hindering of the spectral energy transfer. Therefore the energy is accumulated and redistributed so that the efficiency of non-linear interactions will be sufficient to provide a constant energy flux. We confirm our phenomenology by high Reynolds number numerical simulations based on a shell model of helical turbulence. The energy in our model is injected at a certain large scale only, whereas the source of helicity is distributed over all scales. In particular, we found that the helical bottleneck effect can appear in the inertial interval of the energy spectrum. △ Less

Submitted 10 November, 2015; v1 submitted 28 August, 2015; originally announced August 2015.

Comments: 5 pages, 2 figures, accepted in PRL

Journal ref: Phys. Rev. Lett. 115, 234501 (2015)

The formation of regular interarm magnetic fields in spiral galaxies

Authors: David Moss, Rodion Stepanov, Marita Krause, Rainer Beck, Dmitri Sokoloff

Abstract: Observations of regular magnetic fields in several nearby galaxies reveal magnetic arms situated between the material arms. The nature of these magnetic arms is a topic of active debate. Previously we found a hint that taking into account the effects of injections of small-scale magnetic fields generated, e.g., by turbulent dynamo action, into the large-scale galactic dynamo can result in magnetic… ▽ More Observations of regular magnetic fields in several nearby galaxies reveal magnetic arms situated between the material arms. The nature of these magnetic arms is a topic of active debate. Previously we found a hint that taking into account the effects of injections of small-scale magnetic fields generated, e.g., by turbulent dynamo action, into the large-scale galactic dynamo can result in magnetic arm formation. We now investigate the joint roles of an arm/interarm turbulent diffusivity contrast and injections of small-scale magnetic field on the formation of large-scale magnetic field ("magnetic arms") in the interarm region. We use the relatively simple "no-$z$" model for the galactic dynamo. This involves projection on to the galactic equatorial plane of the azimuthal and radial magnetic field components; the field component orthogonal to the galactic plane is estimated from the solenoidality condition. We find that addition of diffusivity gradients to the effect of magnetic field injections makes the magnetic arms much more pronounced. In particular, the regular magnetic field component becomes larger in the interarm space compared to that within the material arms.The joint action of the turbulent diffusivity contrast and small-scale magnetic field injections (with the possible participation of other effects previously suggested) appears to be a plausible explanation for the phenomenon of magnetic arms. △ Less

Submitted 18 May, 2015; v1 submitted 29 April, 2015; originally announced April 2015.

Comments: 9 pages, 9 figures

Journal ref: A&A 578, A94 (2015)

Measuring magnetism in the Milky Way with the Square Kilometre Array

Authors: Marijke Haverkorn, Takuya Akahori, Ettore Carretti, Katia Ferriere, Peter Frick, Bryan Gaensler, George Heald, Melanie Johnston-Hollitt, David Jones, Tom Landecker, Sui Ann Mao, Aris Noutsos, Niels Oppermann, Wolfgang Reich, Timothy Robishaw, Anna Scaife, Dominic Schnitzeler, Rodion Stepanov, Xiaohui Sun, Russ Taylor

Abstract: Magnetic fields in the Milky Way are present on a wide variety of sizes and strengths, influencing many processes in the Galactic ecosystem such as star formation, gas dynamics, jets, and evolution of supernova remnants or pulsar wind nebulae. Observation methods are complex and indirect; the most used of these are a grid of rotation measures of unresolved polarized extragalactic sources, and broa… ▽ More Magnetic fields in the Milky Way are present on a wide variety of sizes and strengths, influencing many processes in the Galactic ecosystem such as star formation, gas dynamics, jets, and evolution of supernova remnants or pulsar wind nebulae. Observation methods are complex and indirect; the most used of these are a grid of rotation measures of unresolved polarized extragalactic sources, and broadband polarimetry of diffuse emission. Current studies of magnetic fields in the Milky Way reveal a global spiral magnetic field with a significant turbulent component; the limited sample of magnetic field measurements in discrete objects such as supernova remnants and HII regions shows a wide variety in field configurations; a few detections of magnetic fields in Young Stellar Object jets have been published; and the magnetic field structure in the Galactic Center is still under debate. The SKA will unravel the 3D structure and configurations of magnetic fields in the Milky Way on sub-parsec to galaxy scales, including field structure in the Galactic Center. The global configuration of the Milky Way disk magnetic field, probed through pulsar RMs, will resolve controversy about reversals in the Galactic plane. Characteristics of interstellar turbulence can be determined from the grid of background RMs. We expect to learn to understand magnetic field structures in protostellar jets, supernova remnants, and other discrete sources, due to the vast increase in sample sizes possible with the SKA. This knowledge of magnetic fields in the Milky Way will not only be crucial in understanding of the evolution and interaction of Galactic structures, but will also help to define and remove Galactic foregrounds for a multitude of extragalactic and cosmological studies. △ Less

Submitted 2 January, 2015; originally announced January 2015.

Comments: 19 pages, 2 figures; to appear as part of 'Cosmic Magnetism' in Proceedings 'Advancing Astrophysics with the SKA (AASKA14)', PoS(AASKA14)096

Joint inverse cascade of magnetic energy and magnetic helicity in MHD turbulence

Authors: Rodion Stepanov, Peter Frick, Irina Mizeva

Abstract: We show that oppositely directed fluxes of energy and magnetic helicity coexist in the inertial range in fully developed magnetohydrodynamic (MHD) turbulence with small-scale sources of magnetic helicity. Using a helical shell model of MHD turbulence, we study the high Reynolds number magnetohydrodynamic turbulence for helicity injection at a scale that is much smaller than the scale of energy inj… ▽ More We show that oppositely directed fluxes of energy and magnetic helicity coexist in the inertial range in fully developed magnetohydrodynamic (MHD) turbulence with small-scale sources of magnetic helicity. Using a helical shell model of MHD turbulence, we study the high Reynolds number magnetohydrodynamic turbulence for helicity injection at a scale that is much smaller than the scale of energy injection. In a short range of scales larger than the forcing scale of magnetic helicity, a bottleneck-like effect appears, which results in a local reduction of the spectral slope. The slope changes in a domain with a high level of relative magnetic helicity, which determines that part of the magnetic energy related to the helical modes at a given scale. If the relative helicity approaches unity, the spectral slope tends to $-3/2$. We show that this energy pileup is caused by an inverse cascade of magnetic energy associated with the magnetic helicity. This negative energy flux is the contribution of the pure magnetic-to-magnetic energy transfer, which vanishes in the non-helical limit. In the context of astrophysical dynamos, our results indicate that a large-scale dynamo can be affected by the magnetic helicity generated at small scales. The kinetic helicity, in particular, is not involved in the process at all. An interesting finding is that an inverse cascade of magnetic energy can be provided by a small-scale source of magnetic helicity fluctuations without a mean injection of magnetic helicity. △ Less

Submitted 10 December, 2014; originally announced December 2014.

Comments: 5 pages, accepted by ApJL

A systematic bias in the calculation of spectral density from a 3D spatial grid

Authors: Rodion Stepanov, Franck Plunian, Mouloud Kessar, Guillaume Balarac

Abstract: The energy spectral density $E(k)$, where $k$ is the spatial wave number, is a well-known diagnostic of homogeneous turbulence and magnetohydrodynamic turbulence. However in most of the curves plotted by different authors, some systematic kinks can be observed at $k=9$, $k=15$ and $k=19$. We claim that these kinks have no physical meaning, and are in fact the signature of the method which is used… ▽ More The energy spectral density $E(k)$, where $k$ is the spatial wave number, is a well-known diagnostic of homogeneous turbulence and magnetohydrodynamic turbulence. However in most of the curves plotted by different authors, some systematic kinks can be observed at $k=9$, $k=15$ and $k=19$. We claim that these kinks have no physical meaning, and are in fact the signature of the method which is used to estimate $E(k)$ from a 3D spatial grid. In this paper we give another method, in order to get rid of the spurious kinks and to estimate $E(k)$ much more accurately. △ Less

Submitted 30 October, 2014; originally announced October 2014.

Comments: 4 pages, in PRE

Comparison of algorithms for determination of rotation measure and Faraday structure I. 1100 - 1400 MHz

Authors: X. H. Sun, L. Rudnick, Takuya Akahori, C. S. Anderson, M. R. Bell, J. D. Bray, J. S. Farnes, S. Ideguchi, K. Kumazaki, T. O'Brien, S. P. O'Sullivan, A. M. M. Scaife, R. Stepanov, J. Stil, K. Takahashi, R. J. van Weeren, M. Wolleben

Abstract: (abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying sc… ▽ More (abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying science: a) an average RM weighted by polarized intensity, RMwtd, b) the separation $Δφ$ of two Faraday components and c) the reduced chi-squared. Based on the current test data of signal to noise ratio of about 32, we find that: (1) When only one Faraday thin component is present, most methods perform as expected, with occasional failures where two components are incorrectly found; (2) For two Faraday thin components, QU-fitting routines perform the best, with errors close to the theoretical ones for RMwtd, but with significantly higher errors for $Δφ$. All other methods including standard Faraday synthesis frequently identify only one component when $Δφ$ is below or near the width of the Faraday point spread function; (3) No methods, as currently implemented, work well for Faraday thick components due to the narrow bandwidth; (4) There exist combinations of two Faraday components which produce a large range of acceptable fits and hence large uncertainties in the derived single RMs; in these cases, different RMs lead to the same Q, U behavior, so no method can recover a unique input model. △ Less

Submitted 15 September, 2014; originally announced September 2014.

Comments: 15 pages, 8 figures, 5 tables, accepted by AJ. Some figures have been degraded, for a full resolution version, see http://physics.usyd.edu.au/~xhsun/datachallengeI.pdf

Faraday signature of magnetic helicity from reduced depolarization

Authors: Axel Brandenburg, Rodion Stepanov

Abstract: Using one-dimensional models, we show that a helical magnetic field with an appropriate sign of helicity can compensate the Faraday depolarization resulting from the superposition of Faraday-rotated polarization planes from a spatially extended source. For radio emission from a helical magnetic field, the polarization as a function of the square of the wavelength becomes asymmetric with respect to… ▽ More Using one-dimensional models, we show that a helical magnetic field with an appropriate sign of helicity can compensate the Faraday depolarization resulting from the superposition of Faraday-rotated polarization planes from a spatially extended source. For radio emission from a helical magnetic field, the polarization as a function of the square of the wavelength becomes asymmetric with respect to zero. Mathematically speaking, the resulting emission occurs then either at observable or at unobservable (imaginary) wavelengths. We demonstrate that rotation measure (RM) synthesis allows for the reconstruction of the underlying Faraday dispersion function in the former case, but not in the latter. The presence of positive magnetic helicity can thus be detected by observing positive RM in highly polarized regions in the sky and negative RM in weakly polarized regions. Conversely, negative magnetic helicity can be detected by observing negative RM in highly polarized regions and positive RM in weakly polarized regions. The simultaneous presence of two magnetic constituents with opposite signs of helicity is shown to possess signatures that can be quantified through polarization peaks at specific wavelengths and the gradient of the phase of the Faraday dispersion function. Similar polarization peaks can tentatively also be identified for the bi-helical magnetic fields that are generated self-consistently by a dynamo from helically forced turbulence, even though the magnetic energy spectrum is then continuous. Finally, we discuss the possibility of detecting magnetic fields with helical and non-helical properties in external galaxies using the Square Kilometre Array. △ Less

Submitted 13 April, 2014; v1 submitted 16 January, 2014; originally announced January 2014.

Comments: 12 pages, 12 figures, ApJ, in press (with 3-D turbulence results now included)

Report number: NORDITA-2014-1

Journal ref: Astrophys. J. 786, 91 (2014)

The relation between magnetic and material arms in models for spiral galaxies

Authors: D. Moss, R. Beck, D. Sokoloff, R. Stepanov, M. Krause, T. G. Arshakian

Abstract: Context. Observations of polarized radio emission show that large-scale (regular) magnetic fields in spiral galaxies are not axisymmetric, but generally stronger in interarm regions. In some nearby galaxies such as NGC 6946 they are organized in narrow magnetic arms situated between the material spiral arms. Aims. The phenomenon of magnetic arms and their relation to the optical spiral arms (the m… ▽ More Context. Observations of polarized radio emission show that large-scale (regular) magnetic fields in spiral galaxies are not axisymmetric, but generally stronger in interarm regions. In some nearby galaxies such as NGC 6946 they are organized in narrow magnetic arms situated between the material spiral arms. Aims. The phenomenon of magnetic arms and their relation to the optical spiral arms (the material arms) call for an explanation in the framework of galactic dynamo theory. Several possibilities have been suggested but are not completely satisfactory; here we attempt a consistent investigation. Methods. We use a 2D mean-field dynamo model in the no-z approximation and add injections of small-scale magnetic field, taken to result from supernova explosions, to represent the effects of dynamo action on smaller scales. This injection of small scale field is situated along the spiral arms, where star-formation mostly occurs. Results. A straightforward explanation of magnetic arms as a result of modulation of the dynamo mechanism by material arms struggles to produce pronounced magnetic arms, at least with realistic parameters, without introducing new effects such as a time lag between Coriolis force and α-effect. In contrast, by taking into account explicitly the small-scale magnetic field that is injected into the arms by the action of the star forming regions that are concentrated there, we can obtain dynamo models with magnetic structures of various forms that can be compared with magnetic arms. (abbrev). Conclusions. We conclude that magnetic arms can be considered as coherent magnetic structures generated by large-scale dynamo action, and associated with spatially modulated small-scale magnetic fluctuations, caused by enhanced star formation rates within the material arms. △ Less

Submitted 4 July, 2013; originally announced July 2013.

Comments: 13 pages, 8 figures, accepted for publication to A&A

Shell Models of Magnetohydrodynamic Turbulence

Authors: Franck Plunian, Rodion Stepanov, Peter Frick

Abstract: Shell models of hydrodynamic turbulence originated in the seventies. Their main aim was to describe the statistics of homogeneous and isotropic turbulence in spectral space, using a simple set of ordinary differential equations. In the eighties, shell models of magnetohydrodynamic (MHD) turbulence emerged based on the same principles as their hydrodynamic counter-part but also incorporating intera… ▽ More Shell models of hydrodynamic turbulence originated in the seventies. Their main aim was to describe the statistics of homogeneous and isotropic turbulence in spectral space, using a simple set of ordinary differential equations. In the eighties, shell models of magnetohydrodynamic (MHD) turbulence emerged based on the same principles as their hydrodynamic counter-part but also incorporating interactions between magnetic and velocity fields. In recent years, significant improvements have been made such as the inclusion of non-local interactions and appropriate definitions for helicities. Though shell models cannot account for the spatial complexity of MHD turbulence, their dynamics are not over simplified and do reflect those of real MHD turbulence including intermittency or chaotic reversals of large-scale modes. Furthermore, these models use realistic values for dimensionless parameters (high kinetic and magnetic Reynolds numbers, low or high magnetic Prandtl number) allowing extended inertial range and accurate dissipation rate. Using modern computers it is difficult to attain an inertial range of three decades with direct numerical simulations, whereas eight are possible using shell models. In this review we set up a general mathematical framework allowing the description of any MHD shell model. The variety of the latter, with their advantages and weaknesses, is introduced. Finally we consider a number of applications, dealing with free-decaying MHD turbulence, dynamo action, Alfven waves and the Hall effect. △ Less

Submitted 1 April, 2013; v1 submitted 18 September, 2012; originally announced September 2012.

Comments: published in Physics Reports

Journal ref: Physics Reports, Volume 523, Issue 1, February 2013, Pages 1-60

An observational test for correlations between cosmic rays and magnetic fields

Authors: Rodion Stepanov, Anvar Shukurov, Andrew Fletcher, Rainer Beck, Laura La Porta, Fatemeh Tabatabaei

Abstract: We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way and M33, from both observations and theory under various assumption about the relation between cosmic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in the Galactic magnetic field (the ratio of the rms fluctuations to the mean magnetic field strength) suggested by Farada… ▽ More We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way and M33, from both observations and theory under various assumption about the relation between cosmic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in the Galactic magnetic field (the ratio of the rms fluctuations to the mean magnetic field strength) suggested by Faraday rotation and synchrotron polarization, the observations are inconsistent with local energy equipartition between cosmic rays and magnetic fields. Our analysis of relative synchrotron intensity fluctuations indicates that the distribution of cosmic rays is nearly uniform at the scales of the order of and exceeding $100\p$, in contrast to strong fluctuations in the interstellar magnetic field at those scales. A conservative upper limit on the ratio of the the fluctuation magnitude in the cosmic ray number density to its mean value is 0.2--0.4 at scales of order 100\,pc. Our results are consistent with a mild anticorrelation between cosmic-ray and magnetic energy densities at these scales, in both the Milky Way and M33. Energy equipartition between cosmic rays and magnetic fields may still hold, but at scales exceeding 1\,kpc. Therefore, we suggest that equipartition estimates be applied to the observed synchrotron intensity smoothed to a linear scale of kiloparsec order (in spiral galaxies) to obtain the cosmic ray distribution and a large-scale magnetic field. Then the resulting cosmic ray distribution can be used to derive the fluctuating magnetic field strength from the data at the original resolution. The resulting random magnetic field is likely to be significantly stronger than existing estimates. △ Less

Submitted 19 September, 2013; v1 submitted 2 May, 2012; originally announced May 2012.

Comments: submitted to MNRAS

Journal ref: MNRAS 437 (2014) 2201-2216

Recognizing magnetic structures by present and future radio telescopes with Faraday rotation measure synthesis

Authors: R. Beck, P. Frick, R. Stepanov, D. Sokoloff

Abstract: We investigate whether the method of wavelet-based Faraday rotation measure (RM) Synthesis can help us to identify structures of regular and turbulent magnetic fields in extended magnetized objects, such as galaxies and galaxy clusters. Wavelets allow us to reformulate the RM synthesis method in a scale-dependent way and to visualize the data as a function of Faraday depth and scale. We present ob… ▽ More We investigate whether the method of wavelet-based Faraday rotation measure (RM) Synthesis can help us to identify structures of regular and turbulent magnetic fields in extended magnetized objects, such as galaxies and galaxy clusters. Wavelets allow us to reformulate the RM synthesis method in a scale-dependent way and to visualize the data as a function of Faraday depth and scale. We present observational tests to recognize magnetic field structures. A region with a regular magnetic field generates a broad "disk" in Faraday space ("Faraday spectrum"), with two "horns" when the distribution of cosmic-ray electrons is broader than that of the thermal electrons. Each magnetic field reversal generates one asymmetric "horn" on top of the "disk". A region with a turbulent field can be recognized as a "Faraday forest" of many components. These tests are applied to the spectral ranges of various synthesis radio telescopes. We argue that the ratio of maximum to minimum wavelengths determines the range of scales that can be identified in Faraday space. A reliable recognition of magnetic field structures requires the analysis of data cubes in position-position-Faraday depth space ("PPF cubes"), observed over a wide and continuous wavelength range, allowing the recognition of a wide range of scales as well as high resolution in Faraday space. The planned Square Kilometre Array (SKA) will fulfill this condition and will be close to representing a perfect "Faraday telescope". The combination of data from the Low Frequency Array (LOFAR) and the Expanded Very Large Array (EVLA) appears to be a promising approach for the recognition of magnetic structures on all scales. The addition of data at intermediate frequencies from the Westerbork Synthesis Radio Telescope (WSRT) or the Giant Meterwave Radio Telescope (GMRT) would fill the gap between the LOFAR and EVLA frequency ranges. △ Less

Submitted 21 June, 2012; v1 submitted 25 April, 2012; originally announced April 2012.

MSC Class: 85 ACM Class: J.2

Multiscale magnetic fields in spiral galaxies: evolution and reversals

Authors: David Moss, Rodion Stepanov, Tigran Arshakian, Rainer Beck, Marita Krause, Dmitri Sokoloff

Abstract: Magnetic fields in nearby, star-forming galaxies reveal both large-scale patterns and small-scale structures. A large-scale field reversal may exist in the Milky Way but no such reversals have been observed so far in external galaxies. The effects of star-forming regions of galaxies need to be included when modelling the evolution of their magnetic fields, which can then be compared to future radi… ▽ More Magnetic fields in nearby, star-forming galaxies reveal both large-scale patterns and small-scale structures. A large-scale field reversal may exist in the Milky Way but no such reversals have been observed so far in external galaxies. The effects of star-forming regions of galaxies need to be included when modelling the evolution of their magnetic fields, which can then be compared to future radio polarization observations. The causes of large-scale field reversals also need clarification. Our model of field evolution in isolated disc galaxies includes a standard mean-field dynamo and continuous injection of turbulent fields (the effect of supernova explosions) in discrete star forming regions by implicit small-scale dynamo action. Synthetic maps of radio synchrotronemission and Faraday rotation measures are computed. A large-scale dynamo is essential to obtain regular large-scale spiral magnetic fields, observed in many galaxies. These appear, on kpc scales in near energy equilibrium with the turbulence, after 1-2 Gyr (redshift 4-3). Turbulent field injection generates small-scale field structures. Depending on model parameters, large-scale field reversals may persist over many Gyrs and can survive until the present epoch. Significant polarized radio synchrotron emission from young galaxies is expected at redshifts less than 4. Faraday rotation measures (RM) are crucial to detect field reversals. Large-scale patterns ofrotation measures can be observed at redshifts less than 3. Our model can explain the general form of axisymmetric spiral fields with many local distortions, as observed in nearby galaxies. For a slightly different choice of parameters, large-scale field reversals can persist over the lifetime of a galaxy. Comparison of our synthetic maps with future observations of distant galaxies with the planned Square Kilometre Array (SKA) will allow refinement of models. △ Less

Submitted 19 September, 2011; originally announced September 2011.

Faraday Rotation Measure Synthesis for Magnetic Fields of Galaxies

Authors: P. Frick, D. Sokoloff, R. Stepanov, R. Beck

Abstract: RM Synthesis was recently developed as a new tool for the interpretation of polarized emission data in order to separate the contributions of different sources lying on the same line of sight. Until now the method was mainly applied to discrete sources in Faraday space (Faraday screens). Here we consider how to apply RM Synthesis to reconstruct the Faraday dispersion function, aiming at the furthe… ▽ More RM Synthesis was recently developed as a new tool for the interpretation of polarized emission data in order to separate the contributions of different sources lying on the same line of sight. Until now the method was mainly applied to discrete sources in Faraday space (Faraday screens). Here we consider how to apply RM Synthesis to reconstruct the Faraday dispersion function, aiming at the further extraction of information concerning the magnetic fields of extended sources, e.g. galaxies. The main attention is given to two related novelties in the method, i.e. the symmetry argument in Faraday space and the wavelet technique. We give a relation between our method and the previous applications of RM Synthesis to point-like sources. We demonstrate that the traditional RM Synthesis for a point-like source indirectly implies a symmetry argument and, in this sense, can be considered as a particular case of the method presented here. Investigating the applications of RM Synthesis to polarization details associated with small-scale magnetic fields, we isolate an option which was not covered by the ideas of the Burn theory, i.e. using quantities averaged over small-scale fluctuations of magnetic field and electron density. We describe the contribution of small-scale fields in terms of Faraday dispersion and beam depolarization. We consider the complex polarization for RM Synthesis without any averaging over small-scale fluctuations of magnetic field and electron density and demonstrate that it allows us to isolate the contribution from small-scale field. △ Less

Submitted 21 February, 2011; originally announced February 2011.

Comments: 10 pages, 12 figures, accepted in MNRAS

Modeling the total and polarized emission in evolving galaxies: "spotty" magnetic structures

Authors: T. G. Arshakian, R. Stepanov, R. Beck, M. Krause, D. Sokoloff

Abstract: Future radio observations with the SKA and its precursors will be sensitive to trace spiral galaxies and their magnetic field configurations up to redshift $z\approx3$. We suggest an evolutionary model for the magnetic configuration in star-forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron emission, and the Faraday rotation measures for disk ga… ▽ More Future radio observations with the SKA and its precursors will be sensitive to trace spiral galaxies and their magnetic field configurations up to redshift $z\approx3$. We suggest an evolutionary model for the magnetic configuration in star-forming disk galaxies and simulate the magnetic field distribution, the total and polarized synchrotron emission, and the Faraday rotation measures for disk galaxies at $z\la 3$. Since details of dynamo action in young galaxies are quite uncertain, we model the dynamo action heuristically relying only on well-established ideas of the form and evolution of magnetic fields produced by the mean-field dynamo in a thin disk. We assume a small-scale seed field which is then amplified by the small-scale turbulent dynamo up to energy equipartition with kinetic energy of turbulence. The large-scale galactic dynamo starts from seed fields of 100 pc and an averaged regular field strength of 0.02\,$μ$G, which then evolves to a "spotty" magnetic field configuration in about 0.8\,Gyr with scales of about one kpc and an averaged regular field strength of 0.6\,$μ$G. The evolution of these magnetic spots is simulated under the influence of star formation, dynamo action, stretching by differential rotation of the disk, and turbulent diffusion. The evolution of the regular magnetic field in a disk of a spiral galaxy, as well as the expected total intensity, linear polarization and Faraday rotation are simulated in the rest frame of a galaxy at 5\,GHz and 150\,MHz and in the rest frame of the observer at 150\,MHz. We present the corresponding maps for several epochs after disk formation. (abridged) △ Less

Submitted 17 March, 2011; v1 submitted 9 February, 2011; originally announced February 2011.

Comments: 13 pages, 6 figures, 2 tables, revised version is accepted for publication in Astronomische Nachrichten

Dissipation scales of kinetic helicities in turbulence

Authors: T. Lessinnes, F. Plunian, R. Stepanov, D. Carati

Abstract: A systematic study of the influence of the viscous effect on both the spectra and the nonlinear fluxes of conserved as well as non conserved quantities in Navier-Stokes turbulence is proposed. This analysis is used to estimate the helicity dissipation scale which is shown to coincide with the energy dissipation scale. However, it is shown using the decomposition of helicity into eigen modes of the… ▽ More A systematic study of the influence of the viscous effect on both the spectra and the nonlinear fluxes of conserved as well as non conserved quantities in Navier-Stokes turbulence is proposed. This analysis is used to estimate the helicity dissipation scale which is shown to coincide with the energy dissipation scale. However, it is shown using the decomposition of helicity into eigen modes of the curl operator, that viscous effects have to be taken into account for wave vector smaller than the Kolomogorov wave number in the evolution of these eigen components of the helicity. △ Less

Submitted 28 September, 2010; originally announced September 2010.

Comments: 6 pages, 2 figures, submited to PoF

Long-term decaying evolution of MHD turbulence

Authors: Peter Frick, Rodion Stepanov

Abstract: The free decay of MHD turbulence at large Reynolds numbers is studied numerically using a shell model. We study the statistical properties based on representative sample of realisations (128 realisations for each type of initial conditions) over the period of $10^5$ large-scale turnover times. The performed simulations show that the force-free non-helical MHD turbulence can demonstrate two differe… ▽ More The free decay of MHD turbulence at large Reynolds numbers is studied numerically using a shell model. We study the statistical properties based on representative sample of realisations (128 realisations for each type of initial conditions) over the period of $10^5$ large-scale turnover times. The performed simulations show that the force-free non-helical MHD turbulence can demonstrate two different scenarios of evolution in spite of similar initial conditions. Within the first scenario, the cross-helicity accumulation is so fast that the energy cascade vanishes before significant magnetic energy dissipates. Then the system approaches the state of maximal cross-helicity. Within the second scenario, the cascade process continues to remain active until time $10^4$ in units of large-scale turnover time. Then the magnetic field becomes vastly helical due to magnetic helicity conservation. Thus the magnetic energy does not dissipate with kinetic energy. △ Less

Submitted 22 September, 2010; originally announced September 2010.

Comments: 6 pages, 6 figures, submitted to EPL

Direct Measurement of Effective Magnetic Diffusivity in Turbulent Flow of Liquid Sodium

Authors: Peter Frick, Vitaliy Noskov, Sergey Denisov, Rodion Stepanov

Abstract: The first direct measurements of effective magnetic diffusivity in turbulent flow of electro-conductive fluids (the so-called beta-effect) under magnetic Reynolds number Rm >> 1 are reported. The measurements are performed in a nonstationary turbulent flow of liquid sodium, generated in a closed toroidal channel. The peak level of the Reynolds number reached Re \approx 3 10^6, which corresponds to… ▽ More The first direct measurements of effective magnetic diffusivity in turbulent flow of electro-conductive fluids (the so-called beta-effect) under magnetic Reynolds number Rm >> 1 are reported. The measurements are performed in a nonstationary turbulent flow of liquid sodium, generated in a closed toroidal channel. The peak level of the Reynolds number reached Re \approx 3 10^6, which corresponds to the magnetic Reynolds number Rm \approx 30. The magnetic diffusivity of the liquid metal was determined by measuring the phase shift between the induced and the applied magnetic fields. The maximal deviation of magnetic diffusivity from its basic (laminar) value reaches about 50% . △ Less

Submitted 22 September, 2010; originally announced September 2010.

Comments: 5 pages, 6 figuser, accepted in PRL

Cascades and dissipation ratio in rotating MHD turbulence at low magnetic Prandtl number

Authors: Franck Plunian, Rodion Stepanov

Abstract: A phenomenology of isotropic magnetohydrodynamic turbulence subject to both rotation and applied magnetic field is presented. It is assumed that the triple correlations decay-time is the shortest between the eddy turn-over time and the ones associated to the rotating frequency and Alfvén wave period. For $Pm=1$ it leads to four kinds of piecewise spectra, depending on the four parameters, injectio… ▽ More A phenomenology of isotropic magnetohydrodynamic turbulence subject to both rotation and applied magnetic field is presented. It is assumed that the triple correlations decay-time is the shortest between the eddy turn-over time and the ones associated to the rotating frequency and Alfvén wave period. For $Pm=1$ it leads to four kinds of piecewise spectra, depending on the four parameters, injection rate of energy, magnetic diffusivity, rotation rate and applied field. With a shell model of MHD turbulence (including rotation and applied magnetic field), spectra for $Pm \le 1$ are presented, together with the ratio between magnetic and viscous dissipation. △ Less

Submitted 15 October, 2010; v1 submitted 18 September, 2010; originally announced September 2010.

Comments: 5 figures, 1 table, appear in PRE

Journal ref: Phys.Rev.E82:046311,2010

Helicity detection of the astrophysical magnetic fields from radio emission statistics

Authors: A. Volegova, R. Stepanov

Abstract: We discuss inverse problem of detection turbulence magnetic field helical properties using radio survey observations statistics. In this paper, we present principal solution which connects magnetic helicity and correlation between Faraday rotation measure and polarization degree of radio synchrotron emission. The effect of depolarization plays the main role in this problem and allows to detect m… ▽ More We discuss inverse problem of detection turbulence magnetic field helical properties using radio survey observations statistics. In this paper, we present principal solution which connects magnetic helicity and correlation between Faraday rotation measure and polarization degree of radio synchrotron emission. The effect of depolarization plays the main role in this problem and allows to detect magnetic helicity for certain frequency range of observable radio emission. We show that the proposed method is mainly sensitive to a large-scale magnetic field component. △ Less

Submitted 18 January, 2010; originally announced January 2010.

Comments: 5 pages, 6 figures

Journal ref: JETP Letters, 2009, V. 90, N. 10, pp. 637-641

Exploring the magnetic fields in local and distant galaxies

Authors: T. G. Arshakian, R. Stepanov, R. Beck, M. Krause, D. Sokoloff, P. Frick

Abstract: (abridged) We investigate the possibility to recognize the magnetic field structures in nearby galaxies and to test the cosmological evolution of their large- and small-scale magnetic fields with the SKA and its precursors. We estimate the required density of the background polarized sources detected with the SKA for reliable reconstruction and reconstruction of magnetic field structures in near… ▽ More (abridged) We investigate the possibility to recognize the magnetic field structures in nearby galaxies and to test the cosmological evolution of their large- and small-scale magnetic fields with the SKA and its precursors. We estimate the required density of the background polarized sources detected with the SKA for reliable reconstruction and reconstruction of magnetic field structures in nearby spiral galaxies. The dynamo theory is applied to distant galaxies to explore the evolution of magnetic fields in distant galaxies in the context of a hierarchical dark matter cosmology. Under favorite conditions, a \emph{recognition} of large-scale magnetic structures in local star-forming disk galaxies (at a distance $\la 100$ Mpc) is possible from $\ga 10$ RMs towards background polarized sources. Galaxies with strong turbulence or small inclination need more polarized sources for a statistically reliable recognition. A reliable \emph{reconstruction} of the field structure without precognition needs at least 20 RM values on a cut along the projected minor axis which translates to $\approx1200$ sources towards the galaxy. We demonstrate that early regular fields are already in place at $z \sim 4$ (approximately 1.5 Gyr after the disk formation) in massive gas-rich galaxies ($>10^9$ M$_{\sun}$) which then evolve to Milky-Way type galaxies. Major and minor mergers influence the star formation rate and geometry of the disk which has an effect of shifting the generation of regular fields in disks to later epochs. Predictions of the evolutionary model of regular fields, simulations of the evolution of turbulent and large-scale regular fields, total and polarized radio emission of disk galaxies, as well as future observational tests with the SKA are discussed. △ Less

Submitted 18 February, 2010; v1 submitted 8 December, 2009; originally announced December 2009.

Comments: 4 pages, 3 figures, contribution to the proceedings of the conference "Wide Field Science and Technology for the Square Kilometre Array", Final Meeting of the EC FP6 Programme SKADS 4-6 November 2009, Chateau de Limelette, Belgium (final vesrsion)

Wavelet-based Faraday Rotation Measure Synthesis

Authors: P. Frick, D. Sokoloff, R. Stepanov, R. Beck

Abstract: Faraday Rotation Measure (RM) Synthesis, as a method for analyzing multi-channel observations of polarized radio emission to investigate galactic magnetic fields structures, requires the definition of complex polarized intensity in the range of the negative lambda square. We introduce a simple method for continuation of the observed complex polarized intensity into this domain using symmetry arg… ▽ More Faraday Rotation Measure (RM) Synthesis, as a method for analyzing multi-channel observations of polarized radio emission to investigate galactic magnetic fields structures, requires the definition of complex polarized intensity in the range of the negative lambda square. We introduce a simple method for continuation of the observed complex polarized intensity into this domain using symmetry arguments. The method is suggested in context of magnetic field recognition in galactic disks where the magnetic field is supposed to have a maximum in the equatorial plane. The method is quite simple when applied to a single Faraday-rotating structure on the line of sight. Recognition of several structures on the same line of sight requires a more sophisticated technique. We also introduce a wavelet-based algorithm which allows us to consider a set of isolated structures. The method essentially improves the possibilities for reconstruction of complicated Faraday structures using the capabilities of modern radio telescopes. △ Less

Submitted 2 November, 2009; originally announced November 2009.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS Letters

Deciphering solar turbulence from sunspots records

Authors: Franck Plunian, Graeme Sarson, Rodion Stepanov

Abstract: It is generally believed that sunspots are the emergent part of magnetic flux tubes in the solar interior. These tubes are created at the base of the convection zone and rise to the surface due to their magnetic buoyancy. The motion of plasma in the convection zone being highly turbulent, the surface manifestation of sunspots may retain the signature of this turbulence, including its intermitten… ▽ More It is generally believed that sunspots are the emergent part of magnetic flux tubes in the solar interior. These tubes are created at the base of the convection zone and rise to the surface due to their magnetic buoyancy. The motion of plasma in the convection zone being highly turbulent, the surface manifestation of sunspots may retain the signature of this turbulence, including its intermittency. From direct observations of sunspots, and indirect observations of the concentration of cosmogenic isotopes $^{14}$C in tree rings or $^{10}$Be in polar ice, power spectral densities in frequency are plotted. Two different frequency scalings emerge, depending on whether the Sun is quiescent or active. %magnetic activity is maximum or minimum. From direct observations we can also calculate scaling exponents. These testify to a strong intermittency, comparable with that observed in the solar wind. △ Less

Submitted 22 September, 2009; originally announced September 2009.

Comments: 5 pages, 6 figures, accepted for publication in MNRAS Letters

Evolution of magnetic fields in galaxies in the frame of hierarchical structure formation cosmology: future tests with the SKA

Authors: Tigran G. Arshakian, Rainer Beck, Marita Krause, Dmitry Sokolff, Rodion Stepanov

Abstract: Results from simulations of hierarchical structure formation cosmology provide a tool to develop an evolutionary model of regular magnetic fields coupled to galaxy formation and evolution. We use the dynamo theory to derive the timescales of amplification and ordering of magnetic fields in disk and puffy galaxies. Galaxies similar to the Milky Way formed their disks at $z\approx10$ and regular f… ▽ More Results from simulations of hierarchical structure formation cosmology provide a tool to develop an evolutionary model of regular magnetic fields coupled to galaxy formation and evolution. We use the dynamo theory to derive the timescales of amplification and ordering of magnetic fields in disk and puffy galaxies. Galaxies similar to the Milky Way formed their disks at $z\approx10$ and regular fields of $μ$G strength and a few kpc coherence length were generated within 2 Gyr (at $z\approx3$), but field ordering up to the coherence scale of the galaxy size took another 6 Gyr (at $z\approx0.5$). Giant galaxies formed their disk already at $z\approx10$, allowing more efficient dynamo generation of strong regular fields (with kpc coherence length) already at $z\approx4$. Dwarf galaxies should have hosted fully coherent fields at $z\approx1$. This evolutionary scenario and number of predictions of the model can be tested by measurements of polarized synchrotron emission and Faraday rotation with the planned Square Kilometre Array. This model is used to simulate the evolution of regular fields in disk galaxies and the polarized radio sky as part of the Square Kilometer Array Design Studies (SKADS). △ Less

Submitted 16 September, 2009; originally announced September 2009.

Comments: 6 pages, 3 figures. To be published in "Panoramic Radio Astronomy: Wide-field 1-2 GHz research on galaxy evolution", Proceedings of Science (PRA2009)

Role of Cross Helicity in Cascade Processes of MHD turbulence

Authors: Irina Mizeva, Rodion Stepanov, Peter Frick

Abstract: The purpose of this work is to investigate the spectral properties of the developed isotropic (non-Alfven) MHD turbulence stationary excited by an external force, which injects the cross helicity into the flow simultaneously with the energy. It is shown that the cross helicity blocks the spectral energy transfer in MHD turbulence and results in energy accumulation in the system. This accumulatio… ▽ More The purpose of this work is to investigate the spectral properties of the developed isotropic (non-Alfven) MHD turbulence stationary excited by an external force, which injects the cross helicity into the flow simultaneously with the energy. It is shown that the cross helicity blocks the spectral energy transfer in MHD turbulence and results in energy accumulation in the system. This accumulation proceeds until the vortex intensification compensates the decreasing efficiency of nonlinear interactions. The formula for estimating the average turbulence energy is obtained for the set ratio between the injected helicity and energy. It is remarkable that the turbulence accumulates the injected cross helicity at its low rate injection -- the integral correlation coefficient significantly exceeds the ratio between the injected helicity and the energy. It is shown that the spectrum slope gradually increases from "5/3" to "2" with the cross helicity level. △ Less

Submitted 10 April, 2009; originally announced April 2009.

Comments: 4 pages

Journal ref: Doklady Physics, 2009, Vol. 54, No. 2, pp. 93-97

Full perturbation solution for the flow in a rotated torus

Authors: A. Chupin, R. Stepanov

Abstract: We present a perturbation solution for a pressure-driven fluid flow in a rotating toroidal channel. The analysis shows the difference between the solutions of full and simplified equations studied earlier. The result is found to be reliable for {\it low} Reynolds number ($\R$) as was the case for a previously studied solution for high $\R$. The convergence conditions are defined for the whole ra… ▽ More We present a perturbation solution for a pressure-driven fluid flow in a rotating toroidal channel. The analysis shows the difference between the solutions of full and simplified equations studied earlier. The result is found to be reliable for {\it low} Reynolds number ($\R$) as was the case for a previously studied solution for high $\R$. The convergence conditions are defined for the whole range of governing parameters. The viscous flow exhibits some interesting features in flow pattern and hydrodynamic characteristics. △ Less

Submitted 25 February, 2008; v1 submitted 28 December, 2007; originally announced December 2007.

Comments: 4 pages, 6 figures

Testing the magnetic field models of galaxies with the SKA

Authors: T. G. Arshakian, R. Stepanov, R. Beck, P. Frick, M. Krause

Abstract: The future new-generation radio telescope SKA (Square Kilometre Array) and its precursors will provide a rapidly growing number of polarized radio sources. Hundred and thousands polarized background sources can be measured towards nearby galaxies thus allowing their detailed magnetic field mapping by means of Faraday rotation measures (RM). We aim to estimate the required density of the backgrou… ▽ More The future new-generation radio telescope SKA (Square Kilometre Array) and its precursors will provide a rapidly growing number of polarized radio sources. Hundred and thousands polarized background sources can be measured towards nearby galaxies thus allowing their detailed magnetic field mapping by means of Faraday rotation measures (RM). We aim to estimate the required density of the background polarized sources detected with the SKA for reliable recognition and reconstruction of the magnetic field structure in nearby spiral galaxies. We construct a galaxy model which includes the ionized gas and magnetic field patterns of different azimuthal symmetry (axisymmetric (ASS), bisymmetric (BSS) and quadrisymmetric spiral (QSS), and superpositions) plus a halo magnetic field. RM fluctuations with a Kolmogorov spectrum due to turbulent fields and/or fluctuations in ionized gas density are superimposed. Recognition of magnetic structures is possible from RM towards background sources behind galaxies or a continuous RM map obtained from the diffuse polarized emission from the galaxy itself. Under favourite conditions, about a few dozens of polarized sources are sufficient for a reliable recognition. Reconstruction of the field structure without precognition becomes possible for a large number of background sources. A reliable reconstruction of the field structure needs at least 20 RM values on a cut along the projected minor axis which translates to approximately 1200 sources towards the galaxy. Radio telescopes operating at low frequencies (LOFAR, ASKAP and the low-frequency SKA array) may also be useful instruments for field recognition or reconstruction with the help of RM, if background sources are still significantly polarized at low frequencies (abriged). △ Less

Submitted 3 December, 2007; v1 submitted 30 November, 2007; originally announced November 2007.

Comments: 5 pages, 2 figures; contribution to the proceedings of the meeting "From planets to dark energy: the modern radio universe", 1-5 October, Manchester, UK (minor changes are added in the replaced version)

Magnetic field structures of galaxies derived from analysis of Faraday rotation measures, and perspectives for the SKA

Authors: Rodion Stepanov, Tigran G. Arshakian, Rainer Beck, Peter Frick, Marita Krause

Abstract: The forthcoming new-generation radio telescope SKA (Square Kilometre Array) and its precursors will provide a rapidly growing number of polarized radio sources. Our analysis aims on what can be learned from these sources concerning the structure and evolution of magnetic fields of external galaxies. Recognition of magnetic structures is possible from Faraday rotation measures (${\rm RM}$) toward… ▽ More The forthcoming new-generation radio telescope SKA (Square Kilometre Array) and its precursors will provide a rapidly growing number of polarized radio sources. Our analysis aims on what can be learned from these sources concerning the structure and evolution of magnetic fields of external galaxies. Recognition of magnetic structures is possible from Faraday rotation measures (${\rm RM}$) towards background sources behind galaxies. We construct models for the ionized gas and magnetic field patterns of different azimuthal symmetry (axisymmetric, bisymmetric and quadrisymmetric spiral, and superpositions) plus a halo magnetic field. $\RM$ fluctuations with a Kolmogorov spectrum due to turbulent fields and/or fluctuations in ionized gas density are superimposed. Assuming extrapolated number density counts of polarized sources, we generate a sample of $\RM$ values within the solid angle of the galaxy. Applying various templates, we derive the minimum number of background sources and the minimum quality of the observations. For a large number of sources, reconstruction of the field structure without precognition becomes possible. Any large-scale regular component of the magnetic field can be clearly recognized from $\RM$ data with help of the $χ^2$ criterium. Under favourite conditions, about a few dozens of polarized sources are sufficient for a reliable result. △ Less

Submitted 21 January, 2008; v1 submitted 8 November, 2007; originally announced November 2007.

Comments: 16 pages, 18 figures, accepted for publication in A&A