Important plasma problems in astrophysics

In astrophysics, plasmas occur under very extreme conditions. For example, there are ultrastrong magnetic fields in neutron stars, relativistic plasmas around black holes and in jets, extremely energetic particles such as cosmic rays in the interstellar medium, extremely dense plasmas in accretion disks, and extremely large magnetic Reynolds numbers in the interstellar medium. These extreme limits for astrophysical plasmas make plasma phenomena much simpler to analyze in astrophysics than in the laboratory. An understanding of such phenomena often results in an interesting way, by simply taking the extreme limiting case of a known plasma theory. The author will describe one of the more exciting examples and will attempt to convey the excitement he felt when he was first exposed to it. However, not all plasma astrophysical phenomena are so simple. There are certain important plasma phenomena in astrophysics that have not been so easily resolved. In fact, a resolution of them is blocking significant progress in astrophysical research. They have not yet yielded to attacks by theoretical astrophysicists nor to extensive numerical simulation. The author will attempt to describe one of the more important of these plasma–astrophysical problems, and discuss why its resolution is so important to astrophysics. This significant example is fast, magnetic reconnection. Another significant example is the large-magnetic-Reynolds number magnetohydrodynamics (MHD) dynamos.

The graviton-to-photon conversion effect in a magnetized relativistic lepton plasma is considered. This effect can be important for the possible generation of electromagnetic radiation accompanying coalescence of relativistic compact neutron star – black hole binaries. The relativistic electron-positron plasma can be generated near the surface of a rotating magnetized neutron star (a radio pulsar). The formation of a relativistic compact binary containing a pulsar and a black hole is predicted by the evolution of massive binary stars. Prior to the coalescence of such a binary due to gravitational wave (GW) emission, a fraction of the GW power can be converted in the plasma outflow into a low-frequency electromagnetic (EM) waves, which can lead to additional radio power prior to the coalescence. Using the graviton- to-photon conversion mechanism in an external magnetic field, we calculate the fraction of GW power converted into the EM radiation in the relativistic plasma. The result is found to depend on the neutron star spin period P, plasma Lorentz factor γ and the cascade multiplicity λ, but independent of the neutron star magnetic field: K ≃ 10−35(P/1 s)2(γ/105)2(λ/105)−2. The possibility of the detection of the non-thermal EM counterparts from neutron star – black hole coalescences in the forthcoming GW observations by aLIGO/Virgo detectors is briefly discussed.

High-energy astrophysics focuses on understanding the origins of the most energetic radiation and particles observed or theoretically predicted in the universe. This area of research probes the most mysterious and most difficult-to-observe astrophysical phenomena. Understanding these phenomena is worth the effort, however, because they may be key to understanding the origin of the universe and the fundamentals of particle physics. Recent observations of extreme events such as ultrahigh-energy cosmic rays and gamma ray bursts have led to significant progress in determining their sources and modeling the mechanisms that produced them. Relativistic jets accelerated from massive compact objects, such as neutron stars or black holes, have been implicated as the primary culprits for producing high- energy events. Here we present three Reviews of the more rapidly advancing areas of research, ultrahigh-energy cosmic rays, gamma ray bursts, and relativistic jets. The highest energy cosmic rays are thought to be extragalactic. However, their source remains a mystery, as discussed in a Review by Guenter Sigl (p. [73][1]). About 12 events have been detected by large-area, ground-based detectors that sample the secondary particles of the cosmic rays. These observations, combined with theoretical models, suggest that ultrahigh-energy cosmic rays come from the earliest processes in the formation of the universe, nonstandard particle physics, or acceleration of particles by energetic processes such as gamma ray bursts. Gamma ray bursts—which last for only a few seconds and generate about 1053 ergs per second of energy output, greater than any other astrophysical phenomenon—are now being detected at a rate of about one per day. Peter Meszaros (p. [79][2]) reviews our current understanding of these curious events. The Italian-Dutch satellite BeppoSAX in particular has helped to locate gamma ray bursts, and accurate locations make it much easier to observe the afterglows with ground-based telescopes and satellites. These observations have led to a model that favors the collapse of a massive object into a neutron star or black hole, generating jets that create gamma ray bursts. Relativistic jets, collimated flows of plasma that travel at almost the speed of light, may produce gamma ray bursts, ultrahigh-energy cosmic rays, and many other high-energy phenomena. David Meier, Shinji Koide, and Yutaka Uchida (p. [84][3]) review the latest observations of jet structures, mainly from radio and x-ray wavelength observations, and focus on magnetohydrodynamic simulations of relativistic jets that reproduce these structures. A neutron star or black hole, in a binary system or in the act of formation, will accrete mass and angular momentum. Differential rotation of this mass and momentum may twist the magnetic field lines of the system, producing bipolar jets, and simulations can provide the appropriate conditions to accelerate particles to produce ultra high-energy cosmic rays or create gamma ray bursts in the jet stream. Relativistic jets also make up the centerpiece of the first of three News stories. Mark Sincell (p. [66][4]) describes how the jets of recently discovered objects known as microquasars may shed light on the mechanisms behind quasars—galaxy-spanning structures powered by supermassive black holes. Next, Govert Schilling (p. [68][5]) tells of another attempt to link two similar-looking phenomena: faint high- energy blips known as soft gamma repeaters and anomalous x-ray pulsars. Finally, Richard Stone (p. [70][6]) travels to Tatarstan to report on a project intended to rival much larger Western ground- and space-based attempts to detect gravitational waves from binary pulsars Thus, the combination of numerical simulations, observations, and theory is improving our understanding of the origins of high-energy astrophysical phenomena. The future looks bright for resolving some of these great mysteries in the early 21st century. [1]: /lookup/doi/10.1126/science.291.5501.73 [2]: /lookup/doi/10.1126/science.291.5501.79 [3]: /lookup/doi/10.1126/science.291.5501.84 [4]: /lookup/doi/10.1126/science.10.1126/SCIENCE.291.5501.66 [5]: /lookup/doi/10.1126/science.10.1126/SCIENCE.291.5501.68 [6]: /lookup/doi/10.1126/science.10.1126/SCIENCE.291.5501.70

Coalescing neutron star binaries are believed to be the most reliable sources for ground-based detectors of gravitational waves and likely progenitors of short gamma-ray bursts. In the process of coalescence, magnetic fields of neutron stars can induce interesting observational manifestations and affect the form of gravitational wave signal. In this papaer we use the population synthesis method to model the expected distribution of neutron star magnetic fields during the coalescence under different assumptions on the initial parameters of neutron stars and their magnetic field evolution. We discuss possible elecotrmagnetic phenomena preceding the coalescence of magnetized neutron star binaries and the effect of magnetic field on the gravitational wave signal. We find that a log-normal (Gaussian in logarithms) distribution of the initial magnetic fields of neutron stars, which agrees with observed properties of radio pulsars, produces the distribution of the magnetic field energy during the coalescence that adequately describes the observed luminosity function of short gamma-ray bursts under different assumptions on the field evolution and initial parameters of neutron stars. This agreement lends further support to the model of coalescing neutron star binaries as progenitors of gamma-ray bursts.

We implement a semi-analytic approach for stability analysis, addressing the ongoing uncertainty about stability and structure of neutron star magnetic fields. Applying the energy variational principle, a model system is displaced from its equilibrium state. The related energy density variation is set up analytically, whereas its volume integration is carried out numerically. This facilitates the consideration of more realistic neutron star characteristics within the model compared to analytical treatments. At the same time, our method retains the possibility to yield general information about neutron star magnetic field and composition structures that are likely to be stable. In contrast to numerical studies, classes of parametrized systems can be studied at once, finally constraining realistic configurations for interior neutron star magnetic fields. We apply the stability analysis scheme on polytropic and non-barotropic neutron stars with toroidal, poloidal and mixed fields testing their stability in a Newtonian framework. Furthermore, we provide the analytical scheme for dropping the Cowling approximation in an axisymmetric system and investigate its impact. Our results confirm the instability of simple magnetised neutron star models as well as a stabilisation tendency in the case of mixed fields and stratification. These findings agree with analytical studies whose spectrum of model systems we extend by lifting former simplifications.

The possibility of observationally discriminating between various types of neutron stars, described by different equations of state of the nuclear matter, as well as differentiating neutron stars from other types of exotic objects, like, for example, quark stars, is one of the fundamental problems in contemporary astrophysics. We consider and investigate carefully the possibility that different types of rapidly rotating neutron stars, as well as other type of compact general relativistic objects, can be differentiated from the study of the emission properties of the accretion disks around them. We obtain the energy flux, the temperature distribution and the emission spectrum from the accretion disks around several classes of rapidly rotating neutron stars, described by different equations of state of the neutron matter, and for quark stars, described by the MIT bag model equation of state and in the CFL (Color-Flavor-Locked) phase, respectively. Particular signatures appear in the electromagnetic spectrum, thus leading to the possibility of directly testing the equation of state of the dense matter by using astrophysical observations of the emission spectra from accretion disks.

When a star dies, it leaves a mark on its surrounding environment. The energy from the supernova explosion forms an expanding shock wave that interacts with interstellar and circumstellar material, creating what we know as a supernova remnant (SNR). If the original star has a mass that is greater than or equal to 8 solar masses, this can also lead to the formation of a rapidly rotating neutron star called a pulsar. As these objects evolve, they interact with the surrounding environment, producing non-thermal and thermal emission. For an SNR, its non-thermal emission arises from a population of relativistic particles being accelerated at the shock front of the SNR, while its thermal emission arises from the shock front heating ejecta and and swept-up interstellar medium to X-ray emitting temperatures. For pulsars, their non-thermal emission arises from relativistic particles being accelerated at the termination shock of a pulsar wind. These particles interact with surrounding magnetic fields and ambient photon fields producing synchrotron and inverse Compton emission which we observe as a pulsar wind nebula (PWN), while its thermal emission arises from the surface of the neutron star. These properties of SNRs and pulsars provide a unique window into studying the acceleration, injection, propagation and interaction of highly energetic particles called cosmic rays with the interstellar medium. In addition, they providing information about the evolution, and dynamics of these objects; properties of the shock fronts; details about the original progenitor star; and the impact that these objects have on their surroundings. The research presented here focuses on analysing the intimate connection between cosmic rays, the non-thermal emission arising from SNRs interacting with molecular clouds, and pulsar wind nebulae; as well as analysing the observational and evolutionary properties of these objects. In this thesis we model the propagation of cosmic rays through the Galaxy in an attempt to characterise a standard cosmic ray background with uncertainties, to reveal the origin of the cosmic ray electron positron anomaly. Furthermore, we analyse the gamma-ray emission from SNRs Kes 79 and MSH 11-61A, which are known to be interacting with molecular clouds, as well as the non-thermal X-ray emission arising from the PWN of PSR J1741-2054. We find that the emission from both SNRs most likely arises from the decay of neutral pions that resulted from the interaction of relativistic ions which are accelerated at the shock-front of a SNR, with ambient material. For PSR J1741-2054, we characterise the properties, minimum magnetic field and minimum energy of the particle population that produces the observed diffuse synchrotron emission that surrounds and trails the pulsar. In addition, we characterise the X-ray emission arising from Kes 79, MSH 11-61A and PSR J1741-2054, in an attempt to shed light on the origin and nature of these objects and their emission. Using X-ray data from XMM-Newton and Suzaku respectively, we probe the temperature, ionisation state, and elemental abundance of the shocked gas of each SNR. This allows us to determine their evolutionary properties, properties of the shock, and mass of the original progenitor; and constrain the density of the X-ray emitting plasma. Using Chandra, we determined the temperature of PSR J1741-2054, as well as characterised its proper motion, velocity, direction of motion, and presence of small scale structure immediately surrounding the pulsar.

The purpose of this chapter is to lay out the basic physical foundations of cosmicray transport in space. We will discuss the basic physics of the cosmic-ray transport equation, first written down by Parker (1965). This equation is remarkably robust and is widely used to study cosmic-ray transport in the solar system and the interstellar medium. The chapter starts with a general background on cosmic rays in the solar system. The transport equation itself is not formally derived, but the basic physical processes that lead to the various terms in this equation are discussed in detail. We will also address what assumptions are made about this equation and under what conditions it is applicable. At the end of this chapter, we provide a few applications related to specific heliophysics phenomena, such as the propagation of solar-energetic particles in the inner heliosphere, the modulation of galactic cosmic rays, and the drift motions of cosmic rays in the solar system. Cosmic rays in the solar system A fundamental and outstanding problem in astrophysics concerns the origin of high-energy charged particles in space. This problem has been known since the early 1900s when Victor Hess performed his famous electroscope experiments on balloon flights to demonstrate that the excess atmospheric radiation increased with increasing distance from the Earth's surface. This led to the discovery of cosmic rays. In the 1920s Chapman and Ferraro considered the effect on Earth of localized and intermittent streams of corpuscular radiation. It was later determined that the solar corpuscular radiation had two components: one that was steady, now known as the solar wind, and the other intermittent, that we now know to be solar energetic particles.

In the presence of a magnetic field, rotational energy can be extracted from black holes via the Blandford-Znajek mechanism. We use self-similar advection-dominated accretion flow (ADAF) models to estimate the efficiency of this mechanism for black holes accreting from geometrically thick disks, in the light of recent magnetohydrodynamic disk simulations, and show that the power from electromagnetic energy extraction exceeds the accretion luminosity for ADAFs at sufficiently low accretion rates. We consider the detectability of isolated stellar mass black holes accreting from the interstellar medium, and show that for any rapidly rotating holes, the efficiency of energy extraction, BZ ≡ LBZ/c2, could reach 10-2. The estimated total luminosity would be consistent with the tentative identification of some EGRET sources as accreting isolated black holes, if that energy is radiated primarily as gamma rays. We discuss the importance of emission from the Blandford-Znajek mechanism for the spectra of other ADAFs, especially those in low-luminosity galactic nuclei.

The origin of PeV cosmic rays (CRs) is a long-standing mystery, and ultrahigh-energy gamma-ray observations would play a crucial role in identifying it. Recently, LHAASO reported the discovery of “dark” gamma-ray sources that were detected above 100 TeV without any GeV–TeV gamma-ray counterparts. The origins of these dark gamma-ray sources are unknown. We propose isolated black holes (IBHs) wandering in molecular clouds as the origins of PeV CRs and LHAASO dark sources. An IBH accretes surrounding dense gas, which forms a magnetically arrested disk (MAD) around the IBH. Magnetic reconnection in the MAD can accelerate CR protons up to PeV energies. CR protons of GeV–TeV energies fall to the IBH, whereas CR protons at sub-PeV energies can escape from the MAD, providing PeV CRs into the interstellar medium. The sub-PeV CR protons interact with the surrounding molecular clouds, producing TeV–PeV gamma rays without emitting GeV–TeV gamma rays. This scenario can explain the dark sources detected by LHAASO. Taking into account the IBH and molecular cloud distributions in our Galaxy, we demonstrate that IBHs can provide a significant contribution to the PeV CRs observed on Earth. Future gamma-ray detectors in the southern sky and neutrino detectors would provide a concrete test to our scenario.

“Cyclotron lines” provide the only direct way to measure the magnetic field of neutron stars. Detected as strong absorption features in the spectra of accreting X‐ray pulsars, they provide a powerful tool to probe the line‐forming region above the neutron star magnetic poles. The lines result from resonant scattering of the X‐ray photons with magnetically quantized electrons in the plasma of accreted material. The interpretation of their positions and observationally clearly resolved structure can be used to constrain the physical parameters of accretion.

Context. In many astrophysical phenomena, and especially in those that involve the high-energy regimes that always accompany the astronomical phenomenology of black holes and neutron stars, physical conditions that are achieved are extreme in terms of speeds, temperatures, and gravitational fields. In such relativistic regimes, numerical calculations are the only tool to accurately model the dynamics of the flows and the transport of radiation in the accreting matter. Aims. We here continue our effort of modelling the behaviour of matter when it orbits or is accreted onto a generic black hole by developing a new numerical code that employs advanced techniques geared towards solving the equations of general-relativistic hydrodynamics. Methods. More specifically, the new code employs a number of high-resolution shock-capturing Riemann solvers and reconstruction algorithms, exploiting the enhanced accuracy and the reduced computational cost of adaptive mesh-refinement (AMR) techniques. In addition, the code makes use of sophisticated ray-tracing libraries that, coupled with general-relativistic radiation-transfer calculations, allow us to accurately compute the electromagnetic emissions from such accretion flows. Results. We validate the new code by presenting an extensive series of stationary accretion flows either in spherical or axial symmetry that are performed either in two or three spatial dimensions. In addition, we consider the highly nonlinear scenario of a recoiling black hole produced in the merger of a supermassive black-hole binary interacting with the surrounding circumbinary disc. In this way, we can present for the first time ray-traced images of the shocked fluid and the light curve resulting from consistent general-relativistic radiation-transport calculations from this process. Conclusions. The work presented here lays the ground for the development of a generic computational infrastructure employing AMR techniques to accurately and self-consistently calculate general-relativistic accretion flows onto compact objects. In addition to the accurate handling of the matter, we provide a self-consistent electromagnetic emission from these scenarios by solving the associated radiative-transfer problem. While magnetic fields are currently excluded from our analysis, the tools presented here can have a number of applications to study accretion flows onto black holes or neutron stars.

The utilization of a conical target irradiated by a high power laser is proposed to study fast magnetic reconnection in relativistic plasma interactions. Such target, placed in front of the near critical density gas jet, splits the laser pulse, forming two parallel laser pulses in the 2D case and a donut shaped pulse in the 3D case. The magnetic annihilation and reconnection occur in the density downramp region of the subsequent gas jet. The magnetic field energy is converted into the particle kinetic energy. As a result, a backward accelerated electron beam is obtained as a signature of reconnection. The above mechanisms are demonstrated using particle-in-cell simulations in both 2D and 3D cases. Facilitating the synchronization of two laser beams, the proposed approach can be used in designing the corresponding experiments on studying fundamental problems of relativistic plasma physics.

The Supernova (SN) is one of the most important and most complex phenomena in astrophysics. Detailed observations of SN require advanced techniques of astronomy and high energy astrophysics, but the theoretical explanation of SN involves virtually every branch of physics. Supernovae, however, offer more than a challenging physics problem because SN are involved in the origin of most of the heavy elements, are the birthplaces of neutron stars, pulsars and probably black holes, control the structure of the interstellar medium, may be responsible for the birth of new stars (and possibly our own solar system), and, of greatest concern in this paper, are involved either directly or indirectly in the origin of the galactic cosmic rays.

The author draws on his memories to review the decisive contributions of Ya B Zeldovich to the formation and development of the theory of disc accretion onto black holes and neutron stars in binaries. A theory developed by N I Shakura and R A Sunyaev in the early 1970s under the guidance of Ya B Zeldovich predicted these objects to be the brightest X-ray sources in the sky and defined the prospects for research in X-ray astronomy and high-energy astrophysics for decades ahead.