Nonlinear force-free reconstruction of the global solar magnetic field: methodology

We present a novel numerical method that allows the calculation of nonlinear force-free magnetostatic solutions above a boundary surface on which only the distribution of the normal magnetic field component is given. The method relies on the theory of force-free electrodynamics and applies directly to the reconstruction of the solar coronal magnetic field for a given distribution of the photospheric radial field component. The method works as follows: we start with any initial magnetostatic global field configuration (e.g. zero, dipole), and along the boundary surface we create an evolving distribution of tangential (horizontal) electric fields that, via Faraday's equation, give rise to a respective normal field distribution approaching asymptotically the target distribution. At the same time, these electric fields are used as boundary condition to numerically evolve the resulting electromagnetic field above the boundary surface, modelled as a thin ideal plasma with non-reflecting, perfectly absorbing outer boundaries. The simulation relaxes to a nonlinear force-free configuration that satisfies the given normal field distribution on the boundary. This is different from existing methods relying on a fixed boundary condition - the boundary evolves toward the a priori given one, at the same time evolving the three-dimensional field solution above it. Moreover, this is the first time a nonlinear force-free solution is reached by using only the normal field component on the boundary. This solution is not unique, but depends on the initial magnetic field configuration and on the evolutionary course along the boundary surface. To our knowledge, this is the first time that the formalism of force-free electrodynamics, used very successfully in other astrophysical contexts, is applied to the global solar magnetic field.Comment: 18 pages, 5 figures, Solar Physic

Nodal points and the transition from ordered to chaotic Bohmian trajectories

We explore the transition from order to chaos for the Bohmian trajectories of a simple quantum system corresponding to the superposition of three stationary states in a 2D harmonic well with incommensurable frequencies. We study in particular the role of nodal points in the transition to chaos. Our main findings are: a) A proof of the existence of bounded domains in configuration space which are devoid of nodal points, b) An analytical construction of formal series representing regular orbits in the central domain as well as a numerical investigation of its limits of applicability. c) A detailed exploration of the phase-space structure near the nodal point. In this exploration we use an adiabatic approximation and we draw the flow chart in a moving frame of reference centered at the nodal point. We demonstrate the existence of a saddle point (called X-point) in the vicinity of the nodal point which plays a key role in the manifestation of exponential sensitivity of the orbits. One of the invariant manifolds of the X-point continues as a spiral terminating at the nodal point. We find cases of Hopf bifurcation at the nodal point and explore the associated phase space structure of the nodal point - X-point complex. We finally demonstrate the mechanism by which this complex generates chaos. Numerical examples of this mechanism are given for particular chaotic orbits, and a comparison is made with previous related works in the literature.Comment: 32 pages, 13 figures, Accepted for publication in Journal of Physics

Asymptotic Orbits in Barred Spiral Galaxies

We study the formation of the spiral structure of barred spiral galaxies, using an $N$-body model. The evolution of this $N$-body model in the adiabatic approximation maintains a strong spiral pattern for more than 10 bar rotations. We find that this longevity of the spiral arms is mainly due to the phenomenon of stickiness of chaotic orbits close to the unstable asymptotic manifolds originated from the main unstable periodic orbits, both inside and outside corotation. The stickiness along the manifolds corresponding to different energy levels supports parts of the spiral structure. The loci of the disc velocity minima (where the particles spend most of their time, in the configuration space) reveal the density maxima and therefore the main morphological structures of the system. We study the relation of these loci with those of the apocentres and pericentres at different energy levels. The diffusion of the sticky chaotic orbits outwards is slow and depends on the initial conditions and the corresponding Jacobi constant.Comment: 17 pages, 24 figure

Invariant manifolds and the response of spiral arms in barred galaxies

The unstable invariant manifolds of the short-period family of periodic orbits around the unstable Lagrangian points $L_1$ and $L_2$ of a barred galaxy define loci in the configuration space which take the form of a trailing spiral pattern. In the present paper we investigate this association in the case of the self-consistent models of Kaufmann & Contopoulos (1996) which provide an approximation of real barred-spiral galaxies. We also examine the relation of `response' models of barred-spiral galaxies with the theory of the invariant manifolds. Our main results are the following: The invariant manifolds yield the correct form of the imposed spiral pattern provided that their calculation is done with the spiral potential term turned on. We provide a theoretical model explaining the form of the invariant manifolds that supports the spiral structure. The azimuthal displacement of the Lagrangian points with respect to the bar's major axis is a crucial parameter in this modeling. When this is taken into account, the manifolds necessarily develop in a spiral-like domain of the configuration space, delimited from below by the boundary of a banana-like non-permitted domain, and from above either by rotational KAM tori or by cantori forming a stickiness zone. We construct `spiral response' models on the basis of the theory of the invariant manifolds and examine the connection of the latter to the `response' models (Patsis 2006) used to fit real barred-spiral galaxies, explaining how are the manifolds related to a number of morphological features seen in such models.Comment: 16 Page

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics.Comment: 11 pages, 4 figures, Astro2020 Science White Paper submitted to the National Academies of Science

Pulsar Pair Cascades in Magnetic Fields with Offset Polar Caps

Neutron star magnetic fields may have polar caps (PC) that are offset from the dipole axis, through field-line sweepback near the light cylinder or non-symmetric currents within the star. The effects of such offsets on electron-positron pair cascades are investigated, using simple models of dipole magnetic fields with small distortions that shift the PCs by different amounts or directions. Using a Monte Carlo pair cascade simulation, we explore the changes in the pair spectrum, multiplicity and energy flux across the PC, as well as the trends in pair flux and pair energy flux with spin-down luminosity, L_{sd}. We also give an estimate of the distribution of heating flux from returning positrons on the PC for different offsets. We find that even modest offsets can produce significant increases in pair multiplicity, especially for pulsars that are near or beyond the pair death lines for centered PCs, primarily because of higher accelerating fields. Pair spectra cover several decades in energy, with the spectral range of millisecond pulsars (MSPs) two orders of magnitude higher than for normal pulsars, and PC offsets allow significant extension of all spectra to lower pair energies. We find that the total PC pair luminosity L_{pair} is proportional to L_{sd}, with L_{pair} ~ 10^{-3} L_{sd} for normal pulsars and L_{pair} ~ 10^{-2} L_{sd} for MSPs. Remarkably, the total PC heating luminosity for even large offsets increases by less than a factor of two, even though the PC area increases by much larger factors, because most of the heating occurs near the magnetic axis.Comment: 41 pages, 17 figures, accepted for publication in Ap

A Tale of Two Current Sheets

I outline a new model of particle acceleration in the current sheet separating the closed from the open field lines in the force-free model of pulsar magnetospheres, based on reconnection at the light cylinder and "auroral" acceleration occurring in the return current channel that connects the light cylinder to the neutron star surface. I discuss recent studies of Pulsar Wind Nebulae, which find that pair outflow rates in excess of those predicted by existing theories of pair creation occur, and use those results to point out that dissipation of the magnetic field in a pulsar's wind upstream of the termination shock is restored to life as a viable model for the solution of the "$\sigma$" problem as a consequence of the lower wind 4-velocity implied by the larger mass loading.Comment: 17 pages, 6 figures, Invited Review, Proceedings of the "ICREA Workshop on The High-Energy Emission from Pulsars and their Systems", Sant Cugat, Spain, April 12-16, 201

Magnetars as Astrophysical Laboratories of Extreme Quantum Electrodynamics: The Case for a Compton Telescope

A next generation of Compton and pair telescopes that improve MeV-band detection sensitivity by more than a decade beyond current instrumental capabilities will open up new insights into a variety of astrophysical source classes. Among these are magnetars, the most highly magnetic of the neutron star zoo, which will serve as a prime science target for a new mission surveying the MeV window. This paper outlines the core questions pertaining to magnetars that can be addressed by such a technology. These range from global magnetar geometry and population trends, to incisive probes of hard X-ray emission locales, to providing cosmic laboratories for spectral and polarimetric testing of exotic predictions of QED, principally the prediction of the splitting of photons and magnetic pair creation. Such fundamental physics cannot yet be discerned in terrestrial experiments. State of the art modeling of the persistent hard X-ray tail emission in magnetars is presented to outline the case for powerful diagnostics using Compton polarimeters. The case highlights an inter-disciplinary opportunity to seed discovery at the interface between astronomy and physics

Two long-term intermittent pulsars discovered in the PALFA Survey

We report the discovery of two long-term intermittent radio pulsars in the ongoing Pulsar Arecibo L-Band Feed Array survey. Following discovery with the Arecibo Telescope, extended observations of these pulsars over several years at Jodrell Bank Observatory have revealed the details of their rotation and radiation properties. PSRs J1910+0517 and J1929+1357 show long-term extreme bimodal intermittency, switching between active (ON) and inactive (OFF) emission states and indicating the presence of a large, hitherto unrecognized underlying population of such objects. For PSR J1929+1357, the initial duty cycle was fON = 0.008, but two years later, this changed quite abruptly to fON = 0.16. This is the first time that a significant evolution in the activity of an intermittent pulsar has been seen, and we show that the spin-down rate of the pulsar is proportional to the activity. The spin-down rate of PSR J1929+1357 is increased by a factor of 1.8 when it is in active mode, similar to the increase seen in the other three known long-term intermittent pulsars. These discoveries increase the number of known pulsars displaying long-term intermittency to five. These five objects display a remarkably narrow range of spin-down power ($\dot{E}\,\sim \,{10}^{32}\,\mathrm{erg}\,{{\rm{s}}}^{-1}$) and accelerating potential above their polar caps. If confirmed by further discoveries, this trend might be important for understanding the physical mechanisms that cause intermittency