Theoretical astrophysics and plasma physics at RPC

Abstract:

We examine the consequences of, and apply, the formalism developed in Terquem (2021) for calculating the rate DR at which energy is exchanged between fast tides and convection. In this previous work, DR (which is proportional to the gradient of the convective velocity) was assumed to be positive in order to dissipate the tidal energy. Here we argue that, even if energy is intermittently transferred from convection to the tides, it must ultimately return to the convective flow and transported efficiently to the stellar surface on the convective timescale. This is consistent with, but much less restrictive than, enforcing DR > 0. Our principle result is a calculation of the circularization timescale of late-type binaries, taking into account the full time evolution of the stellar structure. We find that circularization is very efficient during the PMS phase, inefficient during the MS, and once again efficient when the star approaches the RGB. These results are in much better agreement with observations than earlier theories. We also apply our formalism to hot Jupiters, and find that tidal dissipation in a Jupiter mass planet yields a circularization timescale of 1 Gyr for an orbital period of 3 d, also in good overall agreement with observations. The approach here is novel, and the apparent success of the theory in resolving longstanding timescale puzzles is compelling.

Abstract:

Magnetic reconnection, a plasma process converting magnetic energy to particle kinetic energy, is often invoked to explain magnetic energy releases powering high-energy flares in astrophysical sources including pulsar wind nebulae and black hole jets. Reconnection is usually seen as the (essentially two-dimensional) nonlinear evolution of the tearing instability disrupting a thin current sheet. To test how this process operates in three dimensions, we conduct a comprehensive particle-in-cell simulation study comparing two- and three-dimensional evolution of long, thin current sheets in moderately magnetized, collisionless, relativistically hot electron–positron plasma, and find dramatic differences. We first systematically characterize this process in two dimensions, where classic, hierarchical plasmoid-chain reconnection determines energy release, and explore a wide range of initial configurations, guide magnetic field strengths and system sizes. We then show that three-dimensional (3-D) simulations of similar configurations exhibit a diversity of behaviours, including some where energy release is determined by the nonlinear relativistic drift-kink instability. Thus, 3-D current sheet evolution is not always fundamentally classical reconnection with perturbing 3-D effects but, rather, a complex interplay of multiple linear and nonlinear instabilities whose relative importance depends sensitively on the ambient plasma, minor configuration details and even stochastic events. It often yields slower but longer-lasting and ultimately greater magnetic energy release than in two dimensions. Intriguingly, non-thermal particle acceleration is astonishingly robust, depending on the upstream magnetization and guide field, but otherwise yielding similar particle energy spectra in two and three dimensions. Although the variety of underlying current sheet behaviours is interesting, the similarities in overall energy release and particle spectra may be more remarkable.

Abstract:

The physical picture of interacting magnetic islands provides a useful paradigm for certain plasma dynamics in a variety of physical environments, such as the solar corona, the heliosheath and the Earth's magnetosphere. In this work, we derive an island kinetic equation to describe the evolution of the island distribution function (in area and in flux of islands) subject to a collisional integral designed to account for the role of magnetic reconnection during island mergers. This equation is used to study the inverse transfer of magnetic energy through the coalescence of magnetic islands in two dimensions. We solve our island kinetic equation numerically for three different types of initial distribution: Dirac delta, Gaussian and power-law distributions. The time evolution of several key quantities is found to agree well with our analytical predictions: magnetic energy decays as$\tilde {t}^{-1}$, the number of islands decreases as$\tilde {t}^{-1}$and the averaged area of islands grows as$\tilde {t}$, where$\tilde {t}$is the time normalised to the characteristic reconnection time scale of islands. General properties of the distribution function and the magnetic energy spectrum are also studied. Finally, we discuss the underlying connection of our island-merger models to the (self-similar) decay of magnetohydrodynamic turbulence.