BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

David Tsiklauri

David Tsiklauri appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.SRphysics.plasm-phphysics.space-phastro-ph

Catalog footprint

What is connected

12works
4topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

12 published item(s)

preprint2022arXiv

Enhanced Phase Mixing of Torsional Alfvén Waves in Stratified and Divergent Solar Coronal Structures, Paper I: Linear Solutions

We derive a corrected analytical solution for the propagation and enhanced phase mixing of torsional Alfvén waves, in a potential magnetic field with exponentially divergent field lines, embedded in a stratified solar corona. Further we develop a code named TAWAS which calculates the analytic solution describing torsional Alfvén waves using IDL software language. We then use TAWAS to demonstrate that both our correction to the analytic solution and the inclusion of wave reflection have a significant impact on Alfvén wave damping. We continue to utilise TAWAS by performing a parameter study in order to identify the conditions under which enhanced phase mixing is strongest. We find that phase mixing is the strongest for high frequency Alfvén waves in magnetic fields with highly divergent field lines and without density stratification. We then present a finite difference solver, Wigglewave, which solves the linearised evolution equations for the system directly. Comparing solutions from TAWAS and Wigglewave we see that our analytical solution is accurate within the limits of the WKB approximation but under-reports the wave damping, caused by enhanced phase mixing, beyond the WKB limit. Both TAWAS and Wigglewave solve the linearised governing equations and not the complete nonlinear MHD equations. Paper II will consider simulations that solve the full MHD equations including important nonlinear effects.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Enhanced Phase Mixing of Torsional Alfvén Waves in Stratified and Divergent Solar Coronal Structures, Paper II: Nonlinear Simulations

We use MHD simulations to detect the nonlinear effects of torsional Alfvén wave propagation in a potential magnetic field with exponentially divergent field lines, embedded in a stratified solar corona. In Paper I we considered solutions to the linearised governing equations torsional Alfvén wave propagation and showed, using a finite difference solver we developed named WiggleWave, that in certain scenarios wave damping is stronger than what would be predicted by our analytic solutions. In this paper we consider whether damping would be further enhanced by the presence of nonlinear effects. We begin by deriving the nonlinear governing equations for torsional Alfvén wave propagation and identifying the terms that cause coupling to magnetosonic perturbations. We then compare simulation outputs from an MHD solver called Lare3d, which solves the full set of nonlinear MHD equations, to the outputs from WiggleWave to detect nonlinear effects such as: the excitation of magnetosonic waves by the Alfvén wave, self-interaction of the Alfvén wave through coupling to the induced magnetosonic waves, and the formation of shock waves higher in the atmosphere caused by the steepening of these compressive perturbations. We suggest that the presence of these nonlinear effects in the solar corona would lead to Alfvén wave heating that exceeds the expectation from the phase mixing alone.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

The Effects of Oscillations & Collisions of Emerging Bipolar Regions on the Triggering of Solar Flares

The ability to predict the occurrence of solar flares in advance is important to humankind due to the potential damage they can cause to Earth's environment and infrastructure. It has been shown in Kusano et al. (2012) that a small-scale bipolar region (BR), with its flux reversed relative to the potential component of the overlying field, appearing near the polarity inversion line (PIL) is sufficient to effectively trigger a solar flare. In this study we perform further 3D magnetohydrodynamic simulations to study the effect that the motion of these small-scale BRs has on the effectiveness of flare triggering. The effect of two small-scale BRs colliding is also simulated. The results indicate that the strength of the triggered flare is dependent on how much of the overlying field is disrupted by the BR. Simulations of linear oscillations of the BR showed that oscillations along the PIL increase the flare strength whilst oscillations across the PIL detract from the flare strength. The flare strength is affected more by larger amplitude oscillations but is relatively insensitive to the frequency of oscillations. In the most extreme case the peak kinetic energy of the flare increased more than threefold compared to a non-oscillating BR. Simulations of torsional oscillations of the BR showed a very small effect on the flare strength. Finally, simulations of colliding BRs showed the generation of much stronger flares as the flares triggered by each individual BR coalesce. These results show that significantly stronger flares can result from motion of the BR along the PIL of a sheared field or from the presence of multiple BRs in the same region.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A potential magnetic field calculator for solar physics applications using staggered grids

A program has been designed for accurately generating a potential magnetic field on a staggered grid by extrapolating the magnetic field normal to the photospheric surface. The code first calculates a magnetic potential using the Green's function method and then uses a finite differencing scheme to calculate the magnetic field from the potential. A new finite differencing formula was derived which accounts for grid staggering, it is shown that this formula gives a numerical approximation that is closest to the real potential field. It is also shown that extending the region over which normal photospheric field is specified can improve the accuracy of the potential field produced. The program is a FORTRAN 90 code that can be used to generate potential magnetic field inputs for Lare3d and other MHD solvers that use a staggered grid for magnetic field components. The program can be parallelised to run quickly over multiple computing cores. The code and supporting description are provided in the appendices and can also be found at https://github.com/calboo/Potential_Field_Calculator.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Role of Electron Inertia and Reconnection Dynamics in a Stressed X-point Collapse with a Guide-Field

In previous simulations of collisionless 2D magnetic reconnection it was consistently found that the term in the generalised Ohm's law that breaks the frozen-in condition is the divergence of the electron pressure tensor's non-gyrotropic components. A fully relativistic particle-in-cell (PIC) code was used to model $X$-point collapse with a guide-field in two and three spatial dimensions. We show that in a 2D $X$-point collapse with a guide-field close to the strength of the in-plane field, the increased induced shear flows along the diffusion region lead to a new reconnection regime in which electron inertial terms play a dominant role at the $X$-point. This transition is marked by the emergence of a magnetic island - and hence a second reconnection site - as well as electron flow vortices moving along the current sheet. The reconnection electric field at the $X$-point is shown to exceed all lower guide-field cases for a brief period, indicating a strong burst in reconnection. By extending the simulation to three spatial dimensions it is shown that the locations of vortices along the current sheet (visualised by their $Q$-value) vary in the out-of-plane direction, producing tilted vortex tubes. The vortex tubes on opposite sides of the diffusion region are tilted in opposite directions, similarly to bifurcated current sheets in oblique tearing-mode reconnection. The tilt angles of vortex tubes were compared to a theoretical estimation and were found to be a good match. Particle velocity distribution functions for different guide-field runs, for 2.5D and 3D simulations, are analysed and compared.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The effects of ion mass variation and domain size on octupolar out-of-plane magnetic field generation in collisionless magnetic reconnection

J. Graf von der Pahlen and D. Tsiklauri, Phys. Plas. 21, 060705 (2014), established that the generation of octupolar out-of-plane magnetic field structure in a stressed $X$-point collapse is due to ion currents. The field has a central region, comprising of the well-known qaudrupolar field (quadrupolar components), as well as four additional poles of reversed polarity closer to the corners of the domain (octupolar components). In this extended work, the dependence of the octupolar structure on domain size and ion mass variation is investigated. Simulations show that the strength and spatial structure of the generated octupolar magnetic field is independent of ion to electron mass ratio. Thus showing that ion currents play a significant role in out-of-plane magnetic structure generation in physically realistic scenarios. Simulations of different system sizes show that the width of the octupolar structure remains the same and has a spacial extent of the order of the ion inertial length. The width of the structure thus appears to be independent on boundary condition effects. The length of the octupolar structure however increases for greater domain sizes, prescribed by the external system size. This was found to be a consequence of the structure of the in-plane magnetic field in the outflow region halting the particle flow and thus terminating the in-plane currents that generate the out-of-plane field. The generation of octupolar magnetic field structure is also established in a tearing-mode reconnection scenario. The differences in the generation of the octupolar field and resulting qualitative differences between $X$-point collapse and tearing-mode are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Whistler wave generation by non-gyrotropic, relativistic, electron beams

Particle-in-cell code, EPOCH, is used for studying features of the wave component evident to propagate backwards from the front of the non-gyrotropic, relativistic beam of electrons injected in the Maxwellian, magnetised background plasma with decreasing density profile. According to recent findings presented in Tsiklauri (2011), Schmitz & Tsiklauri (2013) and Pechhacker & Tsiklauri (2012), in a 1.5-dimensional magnetised plasma system, the non-gyrotropic beam generates freely escaping electromagnetic radiation with properties similar to the Type-III solar radio bursts. In this study the backwards propagating wave component evident in the perpendicular components of the elecromagnetic field in such a system is presented for the first time. Background magnetic field strength in the system is varied in order to prove that the backwards propagating wave's frequency, prescribed by the whistler wave dispersion relation, is proportional to the specified magnetic field. Moreover, the identified whistlers are shown to be generated by the normal Doppler-shifted relativistic resonance. Large fraction of the energy of the perpendicular electromagnetic field components is found to be carried away by the whistler waves, while a small but sufficient fraction is going into L- and R- electromagnetic modes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Electron cyclotron maser emission mode coupling to the z-mode on a longitudinal density gradient in the context of solar type III bursts

A beam of super-thermal, hot electrons was injected into maxwellian plasma with a density gradient along a magnetic field line. 1.5D particle-in-cell simulations were carried out which established that the EM emission is produced by the perpendicular component of the beam injection momentum. The beam has a positive slope in the distribution function in perpendicular momentum phase space, which is the characteristic feature of a cyclotron maser. The cyclotron maser in the overdense plasma generates emission at the electron cyclotron frequency. The frequencies of generated waves were too low to propagate away from the injection region, hence the wavelet transform shows a pulsating wave generation and decay process. The intensity pulsation frequency is twice the relativistic cyclotron frequency. Eventually, a stable wave packet formed and could mode couple on the density gradient to reach frequencies of the order of the plasma frequency, that allowed for propagation. The emitted wave is likely to be a z-mode wave. The total electromagnetic energy generated is of the order of 0.1% of the initial beam kinetic energy. The proposed mechanism is of relevance to solar type III radio bursts, as well as other situations, when the injected electron beam has a non-zero perpendicular momentum, e.g. magnetron.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The effect of electron beam pitch angle and density gradient on solar type III radio bursts

1.5D Particle-In-Cell simulations of a hot, low density electron beam injected into magnetized, maxwellian plasma were used to further explore the alternative non-gyrotropic beam driven electromagnetic emission mechanism, first studied in Tsiklauri (2011). Variation of beam injection angle and background density gradient showed that the emission process is caused by the perpendicular component of the beam injection current, whereas the parallel component only produces Langmuir waves, which play no role in the generation of EM waves in our mechanism. Particular emphasis was put on the case, where the beam is injected perpendicularly to the background magnetic field, as this turned off any electrostatic wave generation along the field and left a purely electromagnetic signal in the perpendicular components. The simulations establish the following key findings: i) Initially waves at a few w_ce/gamma are excited, mode converted and emitted at w_pe ii) The emission intensity along the beam axis is proportional to the respective component of the kinetic energy of the beam; iii) The frequency of the escaping EM emission is independent of the injection angle; iv) A stronger background density gradient causes earlier emission; v) The beam electron distribution function in phase space shows harmonic oscillation in the perpendicular components at the relativistic gyrofrequency; vi) The requirement for cyclotron maser emission, df/dv_perp > 0, is fulfilled; vii) The degree of linear polarization of the emission is strongly dependent on the beam injection angle; viii) The generated electromagnetic emission is left-hand elliptically polarized as the pitch angle tends to 90 deg; ix) The generated electromagnetic energy is of the order of 0.1% of the initial beam kinetic energy.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Vlasov-Maxwell, self-consistent electromagnetic wave emission simulations in the solar corona

1.5D Vlasov-Maxwell simulations are employed to model electromagnetic emission generation in a fully self-consistent plasma kinetic model for the first time in the solar physics context. The simulations mimic the plasma emission mechanism and Larmor drift instability in a plasma thread that connects the Sun to Earth with the spatial scales compressed appropriately. The effects of spatial density gradients on the generation of electromagnetic radiation are investigated. It is shown that 1.5D inhomogeneous plasma with a uniform background magnetic field directed transverse to the density gradient is aperiodically unstable to Larmor-drift instability. The latter results in a novel effect of generation of electromagnetic emission at plasma frequency. When density gradient is removed (i.e. when plasma becomes stable to Larmor-drift instability) and a $low$ density, super-thermal, hot beam is injected along the domain, in the direction perpendicular to the magnetic field, plasma emission mechanism generates non-escaping Langmuir type oscillations which in turn generate escaping electromagnetic radiation. It is found that in the spatial location where the beam is injected, the standing waves, oscillating at the plasma frequency, are excited. These can be used to interpret the horizontal strips observed in some dynamical spectra. Quasilinear theory predictions: (i) the electron free streaming and (ii) the beam long relaxation time, in accord with the analytic expressions, are corroborated via direct, fully-kinetic simulation. Finally, the interplay of Larmor-drift instability and plasma emission mechanism is studied by considering $dense$ electron beam in the Larmor-drift unstable (inhomogeneous) plasma. http://www.maths.qmul.ac.uk/~tsiklauri/movie1.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie2.mpg * http://www.maths.qmul.ac.uk/~tsiklauri/movie3.mpg

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Missing bits of the solar jigsaw puzzle: small-scale, kinetic effects in coronal studies

The solar corona, anomalously hot outer atmosphere of the Sun, is traditionally described by magnetohydrodynamic, fluid-like approach. Here we review some recent developments when, instead, a full kinetic description is used. It is shown that some of the main unsolved problems of solar physics, such as coronal heating and solar flare particle acceleration can be viewed in a new light when the small-scale, kinetic plasma description methods are used.

0 citations0 reviewsNo code linkedOpen access
preprint1998arXiv

Heavy neutrino ball as a possible solution to the "blackness problem" of the Galactic center

It has been recently shown (Tsiklauri & Viollier, 1998a) that the matter concentration inferred from observed stellar motion at the galactic center (Eckart & Genzel, 1997, MNRAS, 284, 576 and Genzel et al., 1996, ApJ, 472, 153) is consistent with a supermassive object of $2.5 \times 10^6$ solar masses, composed of self-gravitating, degenerate heavy neutrinos. It has been furthermore suggested (Tsiklauri & Viollier, 1998a) that the neutrino ball scenario may have an advantage that it could possibly explain the so-called "blackness problem" of the galactic center. Here, we present a quantitative investigation of this statement, by calculating the emitted spectrum of Sgr A$^*$ in the framework of standard accretion disk theory.

0 citations0 reviewsNo code linkedOpen access