Simulations of the axisymmetric magnetospheres of neutron stars

In this paper we present the results of time-dependent simulations of the dipolar axisymmetric magnetospheres of neutron stars carried out within the frameworks of both relativistic magnetohydrodynamics (MHD) and resistive force-free electrodynamics. The results of force-free simulations reveal the inability of our numerical method to accommodate the equatorial current sheets of pulsar magnetospheres, and raise a question mark about the robustness of this approach. On the other hand, the MHD approach allows us to make significant progress. We start with a non-rotating magnetically dominated dipolar magnetosphere and follow its evolution as the stellar rotation is switched on. We find that the time-dependent solution gradually approaches a steady state that is very close to the stationary solution of the pulsar equation found in 1999 by Contopoulos, Kazanas & Fendt. This result suggests that other stationary solutions that have the Y-point located well inside the light cylinder are unstable. The role of particle inertia and pressure on the structure and dynamics of MHD magnetospheres is studied in detail, as well as the potential implications of dissipative processes in the equatorial current sheet. We argue that pulsars may have differentially rotating magnetospheres which develop noticeable structural oscillations, and that this may help to explain the nature of the subpulse phenomena.     

Dynamic evolution of a quasi-spherical general polytropic magnetofluid with self-gravity

In various astrophysical contexts, we analyze self-similar behaviours of magnetohydrodynamic (MHD) evolution of a quasi-spherical polytropic magnetized gas under self-gravity with the specific entropy conserved along streamlines. In particular, this MHD model analysis frees the scaling parameter n in the conventional polytropic self-similar transformation from the constraint of n+γ=2 with γ being the polytropic index and therefore substantially generalizes earlier analysis results on polytropic gas dynamics that has a constant specific entropy everywhere in space at all time. On the basis of the self-similar nonlinear MHD ordinary differential equations, we examine behaviours of the magnetosonic critical curves, the MHD shock conditions, and various asymptotic solutions. We then construct global semi-complete self-similar MHD solutions using a combination of analytical and numerical means and indicate plausible astrophysical applications of these magnetized flow solutions with or without MHD shocks.     

Self-Collimated Stationary Disk Winds

In this contribution, we first review the theory of self-collimated jets launched from magnetized accretion disks (disk-winds originating from the first AUs). We show why it is crucial to solve in a self-consistent way the interplay between the resistive accretion disk and the ideal MHD jets. Indeed, this is the only way to get exact values for the disk ejection efficiency ξ (the jet mass load issue). Then, we show self-similar calculations of such accretion-ejection structures: first cold jets, then warm jets obtained in the presence of a hot disk chromosphere. Finally, we present for the first time an accretion-ejection flow crossing all three critical points.     

Interfacing MHD Single Fluid and Kinetic Exospheric Solar Wind Models and Comparing Their Energetics

An exospheric kinetic solar wind model is interfaced with an observation-driven single-fluid magnetohydrodynamic (MHD) model. Initially, a photospheric magnetogram serves as observational input in the fluid approach to extrapolate the heliospheric magnetic field. Then semi-empirical coronal models are used for estimating the plasma characteristics up to a heliocentric distance of 0.1 AU. From there on, a full MHD model that computes the three-dimensional time-dependent evolution of the solar wind macroscopic variables up to the orbit of Earth is used. After interfacing the density and velocity at the inner MHD boundary, we compare our results with those of a kinetic exospheric solar wind model based on the assumption of Maxwell and Kappa velocity distribution functions for protons and electrons, respectively, as well as with in situ observations at 1 AU. This provides insight into more physically detailed processes, such as coronal heating and solar wind acceleration, which naturally arise from including suprathermal electrons in the model. We are interested in the profile of the solar wind speed and density at 1 AU, in characterizing the slow and fast source regions of the wind, and in comparing MHD with exospheric models in similar conditions. We calculate the energetics of both models from low to high heliocentric distances.     

Presupernova structure of massive stars

Issues concerning the structure and evolution of core collapse progenitor stars, and stellar evolution in general, are discussed with an emphasis on interior evolution. We discuss some recent results that address quantifying the uncertainties inherent in modern stellar evolution calculations, and we describe a research effort aimed at investigating the transport and mixing processes associated with stellar turbulence, which is arguably the greatest source of uncertainty in supernova progenitor structure, besides mass loss, at the time of core collapse. We highlight the important role played by precision observations of stellar parameters in constraining theoretical models, as well as the physical insight that can be garnered from three-dimensional hydrodynamic simulation.     

From Forbidden Coronal Lines to Meaningful Coronal Magnetic Fields

We review methods to measure magnetic fields within the corona using the polarized light in magnetic-dipole (M1) lines. We are particularly interested in both the global magnetic-field evolution over a solar cycle, and the local storage of magnetic free energy within coronal plasmas. We address commonly held skepticisms concerning angular ambiguities and line-of-sight confusion. We argue that ambiguities are, in principle, no worse than more familiar remotely sensed photospheric vector fields, and that the diagnosis of M1 line data would benefit from simultaneous observations of EUV lines. Based on calculations and data from eclipses, we discuss the most promising lines and different approaches that might be used. We point to the S-like [Fe xi] line (J=2 to J=1) at 789.2 nm as a prime target line (for the Advanced Technology Solar Telescope (ATST) for example) to augment the hotter 1074.7 and 1079.8 nm Si-like lines of [Fe xiii] currently observed by the Coronal Multi-channel Polarimeter (CoMP). Significant breakthroughs will be made possible with the new generation of coronagraphs, in three distinct ways: i) through single-point inversions (which encompasses also the analysis of MHD wave modes), ii) using direct comparisons of synthetic MHD or force-free models with polarization data, and iii) using tomographic techniques.     

Magnetohydrodynamic behaviour of a finite amplitude disturbance in the solar magnetic field

In this paper, dynamic processes in the solar atmosphere are studied numerically from a complete set of MHD equations. Dynamic evolution of the non-linear magnetic field is produced by the finite amplitude of the azimuthai magnetic field at the base of the flux tube of the solar atmosphere. It is assumed that the initial configuration of the magnetic field is a force-free and potential field, the magnetic field is disturbed at the base, the plasma is driven and a part of the magnetic energy is transformed into the kinetic energy of the plasma.The compressed flow of the plasma has the features of fast MHD waves. The computation results give quantitatively the non-linear evolution of strong magnetic fields. These results could be used in an explanation of coronal transients, surge, spray and eruptive prominence events in the solar atmosphere, as well as in a modelling of plasma behaviour in high-β structure experiments in the laboratory.     

The impact of correlated noise on SuperWASP detection rates for transiting extrasolar planets

The evolution of the Alfvén turbulence due to three-wave interactions is discussed using kinetic theory for a collisionless, thermal plasma. There are three low-frequency modes, analogous to the three modes of compressible magnetohydrodynamics (MHD). When only Alfvén waves are considered, the known anisotropy of turbulence in incompressible MHD theory is reproduced. Inclusion of a fast mode wave leads to the separation of turbulence into two regimes: small wave numbers where three-wave processes involving a fast mode are dominant, and large wave numbers where the three Alfvén wave process is dominant. Possible application of the anisotropic Alfvén turbulence to the interstellar medium and dissipation of magnetic energy in magnetars are discussed.     

Collapse of primordial clouds — II. The role of dark matter

In this article we extend the study performed in our previous article of the collapse of primordial objects. We here analyse the behaviour of the physical parameters for clouds ranging from 10⁷ to 10¹⁵ M_⊙. We study the dynamical evolution of these clouds in two ways: as purely baryonic clouds and as clouds with non-baryonic dark matter included. We start the calculations at the beginning of the recombination era, following the evolution of the structure until the collapse (which we defined as the time when the density contrast of the baryonic matter is greater than 10⁴). We analyse the behaviour of several physical parameters of the clouds (e.g. the density contrast and the velocities of the baryonic matter and the dark matter) as a function of time and radial position in the cloud. In this study all physical processes that are relevant to the dynamical evolution of the primordial clouds, such as for example photon drag (due to the cosmic background radiation) and hydrogen molecular production, besides the expansion of the Universe, are included in the calculations. In particular we find that the clouds with dark matter collapse at higher redshift when we compare the results with the purely baryonic models. As a general result we find that the distribution of the non-baryonic dark matter is more concentrated than the baryonic one. It is important to stress that we do not take into account the putative virialization of the non-baryonic dark matter; we just follow the time and spatial evolution of the cloud, solving its hydrodynamical equations. We also studied the role of cooling–heating processes in the purely baryonic clouds.     

On the MHD Acceleration of Astrophysical Jets

We present a 2.5D magnetohydrodynamic (MHD) simulation of the acceleration of a collimated jet from a magnetized accretion disk. We employ a MHD Adaptive Mesh Refinement (AMR) code (FLASH—University of Chicago). Thanks to this tool we can follow the evolution of the system for many dynamical timescales with a high-spatial resolution. Assuming an initial condition in which a Keplerian disk, thus with no accretion motions, is threaded by a uniform poloidal magnetic field, we show how both the accretion flow and the acceleration of the outflow occur, and we present in detail which are the forces responsible for the jet launching and collimation. Our simulation also shows how the collimating forces due to the self-generated toroidal magnetic field can produce some peculiar knotty features.     

Damping of Longitudinal Magneto–Acoustic Oscillations in Slowly Varying Coronal Plasma

We investigate the propagation of MHD waves in a magnetised plasma in a weakly stratified atmosphere, representative of hot coronal loops. In most earlier studies, a time-independent equilibrium was considered. Here we abandon this restriction and allow the equilibrium to develop as a function of time. In particular, the background plasma is assumed to be cooling due to thermal conduction. The cooling is assumed to occur on a time scale greater than the characteristic travel times of the perturbations. We investigate the influence of cooling of the background plasma on the properties of magneto–acoustic waves. The MHD equations are reduced to a 1D system modelling magneto–acoustic modes propagating along a dynamically cooling coronal loop. A time-dependent dispersion relation that describes the propagation of the magneto–acoustic waves is derived using the WKB theory. An analytic solution for the time-dependent amplitude of waves is obtained, and the method of characteristics is used to find an approximate analytical solution. Numerical calculations of the analytically derived solutions are obtained to give further insight into the behaviour of the MHD waves in a system with a variable, time-dependent background. The results show that there is a strong damping of MHD waves and the damping also appears to be independent of the position along the loop. Studies of MHD wave behaviour in a time-dependent backgrounds seem to be a fundamental and very important next step in the development of MHD wave theory that is applicable to a wide range of situations in solar physics.     

Waves and Oscillations in the Corona – (Invited Review)

It has long been suggested on theoretical grounds that MHD waves must occur in the solar corona, and have important implications for coronal physics. An unequivocal identification of such waves has however proved elusive, though a number of events were consistent with an interpretation in terms of MHD waves. Recent detailed observations of waves in events observed by SOHO and TRACE removes that uncertainty, and raises the importance of MHD waves in the corona to a higher level. Here we review theoretical aspects of how MHD waves and oscillations may occur in a coronal medium. Detailed observations of waves and oscillations in coronal loops, plumes and prominences make feasible the development of coronal seismology, whereby parameters of the coronal plasma (notably the Alfvén speed and through this the magnetic field strength) may be determined from properties of the oscillations. MHD fast waves are refracted by regions of low Alfvén speed and slow waves are closely field-guided, making regions of dense coronal plasma (such as coronal loops and plumes) natural wave guides for MHD waves. There are analogies with sound waves in ocean layers and with elastic waves in the Earth's crust. Recent observations also indicate that coronal oscillations are damped. We consider the various ways this may be brought about, and its implications for coronal heating.     

Resonant Damping of Propagating Kink Waves in Time-Dependent Magnetic Flux Tube

We explore the notion of resonant absorption in a dynamic time-dependent magnetised plasma background. Very many works have investigated resonance in the Alfvén and slow MHD continua under both ideal and dissipative MHD regimes. Jump conditions in static and steady systems have been found in previous works, connecting solutions at both sides of the resonant layer. Here, we derive the jump conditions in a temporally dependent, magnetised, inhomogeneous plasma background to leading order in the Wentzel–Kramers–Billouin (WKB) approximation. Next, we exploit the results found in Williamson and Erdélyi (Solar Phys. 289, 899, 2014) to describe the evolution of the jump condition in the dynamic model considered. The jump across the resonant point is shown to increase exponentially in time. We determined the damping as a result of the resonance over the same time period and investigated the temporal evolution of the damping itself. We found that the damping coefficient, as a result of the evolution of the resonance, decreases as the density gradient across the transitional layer decreases. This has the consequence that in such time-dependent systems resonant absorption may not be as efficient as time progresses.     

Escape from a Crisis in Fokker-Planck Models

Recent N-body simulations have shown that there is a serious discrepancy between the results of N-body simulations and the results of Fokker-Planck simulations for the evolution of globular and rich open clusters under the influence of the galactic tidal field. In some cases, the lifetime obtained from Fokker-Planck calculations is more than an order of magnitude smaller than those from N-body simulations. In this paper we show that the principal cause for this discrepancy is an over-simplified treatment of the tidal field used in previous Fokker-Planck simulations. We performed new Fokker-Planck calculations using a more appropriate implementation for the boundary condition of the tidal field. The implementation is only possible with anisotropic Fokker-Planck models, while all previous Fokker-Planck calculations rely on the assumption of isotropy. Our new Fokker-Planck results agree well with N-body results. Comparison of the two types of simulations gives a better understanding of the evolution of such clusters. This revised version was published online in July 2006 with corrections to the Cover Date.     

Turbulent molecular clouds

Stars form within molecular clouds but our understanding of this fundamental process remains hampered by the complexity of the physics that drives their evolution. We review our observational and theoretical knowledge of molecular clouds trying to confront the two approaches wherever possible. After a broad presentation of the cold interstellar medium and molecular clouds, we emphasize the dynamical processes with special focus to turbulence and its impact on cloud evolution. We then review our knowledge of the velocity, density and magnetic fields. We end by openings towards new chemistry models and the links between molecular cloud structure and star-formation rates.     

Mapping Solar Wind Streams from the Sun to 1 AU: A Comparison of Techniques

A variety of techniques exist for mapping solar wind plasma and magnetic field measurements from one location to another in the heliosphere. Such methods are either applied to extrapolate solar data or coronal model results from near the Sun to 1 AU (or elsewhere), or to map in-situ observations back to the Sun. In this study, we estimate the sensitivity of four models for evolving solar wind streams from the Sun to 1 AU. In order of increasing complexity, these are: i) ballistic extrapolation; ii) ad hoc kinematic mapping; iii) 1-D upwinding propagation; and iv) global heliospheric MHD modeling. We also consider the effects of the interplanetary magnetic field on the evolution of the stream structure. The upwinding technique is a new, simplified method that bridges the extremes of ballistic extrapolation and global heliospheric MHD modeling. It can match the dynamical evolution captured by global models, but is almost as simple to implement and as fast to run as the ballistic approximation.     

Populating the Galaxy with low-mass X-ray binaries

We perform binary population-synthesis calculations to investigate the incidence of low-mass X-ray binaries (LMXBs) and their birth rate in the Galaxy. We use a binary-evolution algorithm that models all the relevant processes including tidal circularization and synchronization. Parameters in the evolution algorithm that are uncertain and may affect X-ray binary formation are allowed to vary during the investigation. We agree with previous studies that under standard assumptions of binary evolution the formation rate and number of black hole (BH) LMXBs predicted by the model are more than an order of magnitude less than what is indicated by observations. We find that the common-envelope process cannot be manipulated to produce significant numbers of BH LMXBs. However, by simply reducing the mass-loss rate from helium stars adopted in the standard model, to a rate that agrees with the latest data, we produce a good match to the observations. Including LMXBs that evolve from intermediate-mass systems also leads to favourable results. We stress that constraints on the X-ray binary population provided by observations are used here merely as a guide as surveys suffer from incompleteness and much uncertainty is involved in the interpretation of results.     

Basic N-body modelling of the evolution of globular clusters – I. Time scaling

We consider the use of N -body simulations for studying the evolution of rich star clusters (i.e. globular clusters).The dynamical processes included in this study are restricted to gravitational (point-mass) interactions, the steady tidal field of a galaxy, and instantaneous mass loss resulting from stellar evolution. With evolution driven by these mechanisms, it is known that clusters fall roughly into two broad classes: those that dissipate promptly in the tidal field, as a result of mass loss; and those that survive long enough for their evolution to become dominated by two-body relaxation.
The time-scales of the processes we consider scale in different ways with the number of stars in the simulation, and the main aim of the paper is to suggest how the scaling of a simulation should be done so that the results are representative of the evolution of a 'real' cluster. We investigate three different ways of scaling time. One of these is appropriate to the first type of cluster, i.e. those that dissipate rapidly; similarly, a second scaling is appropriate only to the second (relaxation-dominated) type. We also develop a hybrid scaling, which is a satisfactory compromise for both types of cluster. Finally we present evidence that the widely used Fokker–Planck method produces models that are in good agreement with N -body models of those clusters that are relaxation-dominated, at least for N -body models with several thousand particles, but that the Fokker–Planck models evolve too fast for clusters that dissipate promptly.     

Self-similar dynamics of a magnetized polytropic gas

In broad astrophysical contexts of large-scale gravitational collapses and outflows and as a basis for various further astrophysical applications, we formulate and investigate a theoretical problem of self-similar magnetohydrodynamics (MHD) for a non-rotating polytropic gas of quasi-spherical symmetry permeated by a completely random magnetic field. Within this framework, we derive two coupled nonlinear MHD ordinary differential equations (ODEs), examine properties of the magnetosonic critical curve, obtain various asymptotic and global semi-complete similarity MHD solutions, and qualify the applicability of our results. Unique to a magnetized gas cloud, a novel asymptotic MHD solution for a collapsing core is established. Physically, the similarity MHD inflow towards the central dense core proceeds in characteristic manners before the gas material eventually encounters a strong radiating MHD shock upon impact onto the central compact object. Sufficiently far away from the central core region enshrouded by such an MHD shock, we derive regular asymptotic behaviours. We study asymptotic solution behaviours in the vicinity of the magnetosonic critical curve and determine smooth MHD eigensolutions across this curve. Numerically, we construct global semi-complete similarity MHD solutions that cross the magnetosonic critical curve zero, one, and two times. For comparison, counterpart solutions in the case of an isothermal unmagnetized and magnetized gas flows are demonstrated in the present MHD framework at nearly isothermal and weakly magnetized conditions. For a polytropic index γ=1.25 or a strong magnetic field, different solution behaviours emerge. With a strong magnetic field, there exist semi-complete similarity solutions crossing the magnetosonic critical curve only once, and the MHD counterpart of expansion-wave collapse solution disappears. Also in the polytropic case of γ=1.25, we no longer observe the trend in the speed-density phase diagram of finding infinitely many matches to establish global MHD solutions that cross the magnetosonic critical curve twice.      

The heating of the solar corona

The heating of the solar corona has been a fundamental astrophysical issue for over sixty years. Over the last decade in particular, space-based solar observatories (Yohkoh, SOHO and TRACE) have revealed the complex and often subtle magnetic-field and plasma interactions throughout the solar atmosphere in unprecedented detail. It is now established that any energy release mechanism is magnetic in origin - the challenge posed is to determine what specific heat input is dominating in a given coronal feature throughout the solar cycle. This review outlines a range of possible magnetohydrodynamic (MHD) coronal heating theories, including MHD wave dissipation and MHD reconnection as well as the accumulating observational evidence for quasi-periodic oscillations and small-scale energy bursts occurring in the corona. Also, we describe current attempts to interpret plasma temperature, density and velocity diagnostics in the light of specific localised energy release. The progress in these investigations expected from future solar missions (Solar-B, STEREO, SDO and Solar Orbiter) is also assessed.Received: 6 February 2003, Published online: 14 November 2003 Correspondence to: R. W. Walsh