PARTICLE ACCELERATION AND TRANSPORT IN RECONNECTING PLASMAS

We model the behaviour of particles in and around X-type magnetic configurations, a possible solar flare reconnection geometry. Particles are accelerated by a time-varying electric field close to the neutral point, and followed by integrating the equations of motion. When their motion becomes adiabatic a stochastic simulation is used to model their further transport in a collisional magnetised medium.     

Rotational fission of contact binary asteroids

The energetics and dynamics of contact binary asteroids as they approach and pass the rotational fission limit is studied. We presume that the asteroids are subject to an external torque, such as from the YORP effect, that increases their angular momentum. Furthermore, we assume the asteroids can be described by a fairly minimal model comprised of a sphere and ellipsoid resting on each other. The minimum energy configurations for contact binary asteroids at different levels of angular momentum are computed and discussed. We find distinct transitions between different configurations as the angular momentum of the system is increased. These indicate that rapidly rotating contact binary asteroids may seek out clearly different relative configurations than slowly rotating systems. We find a single end state of the systems prior to rotational fission, and distinct dynamical outcomes as a function of mass distribution and shape when the rotational fission limit is exceeded. Our theoretical results agree qualitatively with observed properties of near-Earth asteroids, and can be used to help explain the spin-rate barrier, contact binaries, and the observed morphology of most NEO binaries.     

The decay of torque free toroidal magnetic fields confined beneath the solar convection zone

Calculation similar to those of Mestel and Moss (1983) are performed to investigate the decay of a toroidal field through configurations satisfying the torque free condition, imposed by the presence of a poloidal field of dipolar form confined beneath the solar convection zone. It is found that initially stable field configurations diffuse into unstable configurations on time-scales of order a few x 10⁸ yr. The results are similar to those of Tayler (1982) for a simpler field model without any dynamical constraints.     

Experimental Design for the Laboratory Simulation of Magnetized Astrophysical Jets

Concepts of several experimental configurations for the investigation of magnetized jets and their interaction with magnetized environments are presented. In the planned experiments, the plasma jets will be created by laser ablation of shaped targets, while magnetic and electric fields with the required configurations will be produced independently by a pulsed power generator. In particular, the recently coupled Terawatt laser Tomcat and Terawatt pulsed power generator Zebra will be used for experiments.     

Basic Properties of Mutual Magnetic Helicity

We derive the magnetic helicity for configurations formed by flux tubes contained fully or only partially in the spatial domain considered (called closed and open configurations, respectively). In both cases, magnetic helicity is computed as the sum of mutual helicity over all possible pairs of magnetic flux tubes weighted by their magnetic fluxes. We emphasize that these mutual helicities have properties which are not those of mutual inductances in classical circuit theory. For closed configurations, the mutual helicity of two closed flux tubes is their relative winding around each other (known as the Gauss linkage number). For open configurations, the magnetic helicity is derived directly from the geometry of the interlaced flux tubes so it can be computed without reference to a ground state (such as a potential field). We derive the explicit expression in the case of a planar and spherical boundary. The magnetic helicity has two parts. The first one is given only by the relative positions of the flux tubes on the boundary. It is the only part if all flux tubes are arch-shaped. The second part counts the integer number of turns each pair of flux tubes wind about each other. This provides a general method to compute the magnetic helicity with discrete or continuous distributions of magnetic field. The method sets closed and open configurations on an equal level within the same theoretical framework.     

大型太阳活动区磁场的一个模型

用时间缓变的非线性无力场模拟超级活动区（弧岛式大型δ黑子）的磁场位形。这个复杂磁场包含了向量磁场的主要观测特征：正负磁流极端不平衡性（正负磁流之比为１：６）,Ｕ形磁反变线,局域磁场的二极子、四极子差异性。模拟结果厅用来解释一些观测结果：（１）大耀斑主要产生在Ｕ形中性线的磁性混杂区或四极子区（２）Ｕ形反变线的准双极性区几乎没有大耀斑很小。（３）活动区内部的大型旋转运动和磁沲运动会导致四极子场磁拓扑分     

A two-dimensional model for a solar prominence: Effect of an external magnetic field

Using analytical approximations we study the effects of different external magnetic configurations on the half-width, mass, and internal magnetic structure of a quiescent solar prominence, modelled as a thin vertical sheet of cool plasma. Firstly, we build up a zeroth-order model and analyse the effects produced by a potential coronal field or a constant- force-free field. This model allows us to obtain the half-width and mass of the prominence for different values of the external field, pressure and shear angle. Secondly, the effects of these external magnetic configurations on a two-dimensional model proposed by Ballester and Priest (1987) are studied. The main effects are a change of the half-width with height, an increase of the mass, a decrease of the magnetic field strength with height and a change in the shape of the magnetic field lines.     

A study on dynamic processes at late stages in the formation of planetary systems in gas and dust disks

The discovered exoplanetary systems have highly diverse dynamic properties, which differ from those of the Solar System. A single model including planet migration effects and their gravitational interaction is used to investigate the features of dynamic processes that lead to the formation of giant-planet systems with different orbital characteristics. It is shown for a system of four giant planets similar to the Solar System how Type I migration could lead to all the planets being captured into resonant configurations. The resonant motion can continue for a long period of time after the transition to Type II migration and after the dissipation of the gas-and-dust disk. The three-planet system of GJ 876 is used to investigate the migration of the planets inward the orbit of the most massive planet and their capture into low-order resonant configurations under the conditions of Type II migration. A system similar to the exoplanetary system of HD 102272 is used to study the capture into high-order resonances followed by an increase in the orbit’s eccentricity.     

On rotating isothermal configurations in the post-Newtonian approximation of general relativity

The equations which govern the structure of a rotating, truncated isothermal sphere in the post-Newtonian approximation of general relativity are derived and solved numerically. Each model is parameterized by both a rotation and a relativity parameter. The density inside the configurations is tabulated and graphed as a function of both distance from the center and co-latitude. Relativistic gravitational effects are found to pull the models into states which are considerably more centrally condensed than one predicts classically. Rotation tends to flatten the isothermal configurations into oblate spheroids, though for even the largest rotation parameters the degree of flattening is only a few percent. The computed models may be similar to the cores of relativistic star clusters.     

Numerical Simulations of Magnetoacoustic–Gravity Waves in the Solar Atmosphere

We investigate the excitation of magnetoacoustic–gravity waves generated from localized pulses in the gas pressure as well as in the vertical component of velocity. These pulses are initially launched at the top of the solar photosphere, which is permeated by a weak magnetic field. We investigate three different configurations of the background magnetic field lines: horizontal, vertical, and oblique to the gravitational force. We numerically model magnetoacoustic–gravity waves by implementing a realistic (VAL-C) model of the solar temperature. We solve the two-dimensional ideal magnetohydrodynamic equations numerically with the use of the FLASH code to simulate the dynamics of the lower solar atmosphere. The initial pulses result in shocks at higher altitudes. Our numerical simulations reveal that a small-amplitude initial pulse can produce magnetoacoustic–gravity waves, which are later reflected from the transition region due to the large-temperature gradient. The cavities in the lower solar atmosphere are found to have the best conditions to act as a resonator for various oscillations, including their trapping and leakage into the higher atmosphere. Our numerical simulations successfully model the excitation of such wave modes, their reflection and trapping, as well as the associated plasma dynamics.     

Slowly rotating superdense configurations in Rosen's bimetric theory of gravitation

In the limits of the Bimetric Theory of Gravitation the equations of axial-symmetric gravitational field, describing the field of rotating configurations, are obtained and examined. The analytical solutions of these equations out of configurations are found. The equations are solved numerically in the first approximation in the angular velocity for the configurations consisting of the real baryon gas (cf. Sahakian, 1972). The moments of inertia of these configurations in BTG are calculated. It is shown that the equations and the appropriate solution, obtained by Falik and Opher (1979), are incorrect.     

Expulsion of magnetic flux from the core and its dissipation in the crust of a neutron star

We construct a model for the magnetic-field evolution of an isolated neutron star by assuming that its core is a type II superconductor and that the field penetrates the core in the form of magnetic lines (fluxoids). We consider the fluxoid expulsion from the core and the field dissipation in a conducting crust. The magnetic-field evolution is calculated self-consistently by taking into account the inverse effect of crustal magnetic line bending on the fluxoid velocity in the core. We consider the evolution of two magnetic configurations, in which the bulk of the magnetic flux passes through the neutron-star core and crust. The buoyancy of fluxoids and the force from the neutron vortexes are mainly responsible for their expulsion from the core in the former and latter cases, respectively.     

Some Equilibrium of the Restricted Three-Body Problem with Axial Symmetric Primaries and a Gyrostat as Infinitesimal Mass

The restricted three-body problem is reconsidered by replacing the point-like primaries of the classical problem by a pair of axisymmetric rigid bodies which have a plane of symmetry perpendicular to their axes, and the infinitesimal mass by a gyrostat. The conditions for the circular motion of the primaries around their center of mass are stated and they yield the classification of all possible orientations of these bodies into four groups according to the value of their angular velocity. Then the equations of motion of the gyrostat are derived and solved for the equilibrium configurations of the system.     

Order and Chaos in Self-Consistent Galactic Models

We construct and compare two different self-consistent N-body equilibrium configurations of galactic models. The two systems have their origin in cosmological initial conditions selected so that the radial orbit instability appears in one model and gives an E5 type elliptical galaxy, but not in the other that gives an E1 type. We examine their phase spaces using uniformly distributed orbits of test particles in the resulting potential and compare with the distribution of the orbits of the real particles in the two systems. The main types of orbits in both cases are box, tube and chaotic orbits. One main conclusion is that the orbits of the test particles in the 3-dimensional potential are foliated in a way quite close to the foliation of invariant tori in a 2-dimensional potential. The real particles describe orbits having similar foliation. However, their distribution is far from being uniform. The difference between the two models of equilibrium is realized mainly by different balances of the populations of real particles in box and tube orbits.     

Central configurations for the planar Newtonian four-body problem

Planar central configurations of four different masses are analyzed theoretically and computed numerically. We follow Dziobek’s approach to four-body central configurations with a straightforward implicit (in the masses and distances) method of our own in which the fundamental quantities are each the quotient of a directed area divided by the corresponding mass. We apply a new simple numerical algorithm to construct general four-body central configurations. We use this tool to obtain new properties of the symmetric and non-symmetric central configurations. The explicit continuous connection between three-body and four-body central configurations where one of the four masses approaches zero is clarified. Some cases of coorbital 1+3 problems are also considered.     

The role of multipolar magnetic fields in pulsar magnetospheres

We explore the role of complex multipolar magnetic fields in determining physical processes near the surface of rotation powered pulsars. We model the actual magnetic field as the sum of global dipolar and star-centred multipolar fields. In configurations involving axisymmetric and uniform multipolar fields, 'neutral points' and 'neutral lines' exist close to the stellar surface. Also, the curvature radii of magnetic field lines near the stellar surface can never be smaller than the stellar radius, even for very high-order multipoles. Consequently, such configurations are unable to provide an efficient pair-creation process above pulsar polar caps, necessary for plasma mechanisms of generation of pulsar radiation. In configurations involving axisymmetric and non-uniform multipoles, the periphery of the pulsar polar cap becomes fragmented into symmetrically distributed narrow subregions where curvature radii of complex magnetic field lines are less than the radius of the star. The pair-production process is only possible just above these 'favourable' subregions. As a result, the pair plasma flow is confined within narrow filaments regularly distributed around the margin of the open magnetic flux tube. Such a magnetic topology allows us to model the system of 20 isolated subbeams observed in PSR B0943+10 by Deshpande & Rankin. We suggest a physical mechanism for the generation of pulsar radio emission in the ensemble of finite subbeams, based on specific instabilities. We propose an explanation for the subpulse drift phenomenon observed in some long-period pulsars.     

Remarks on the magnetic support of quiescent prominences

The development of magnetic field structures which can lead to prominence configurations of the Kuperus-Raadu type is discussed. Starting from streamer type configurations and preserving the total current in the system we find that simple two-dimensional static configurations lead to prominences which in general lie systematically much lower than the heights found from observations. We therefore conclude that either more complex field configurations are needed to explain the recent observations by Leroy et al. (1983) or the initial configurations must be very special.     

The importance of photospheric magnetic field complexity for coronal energy storage

Possibilities for the storage of energy in coronal electric currents in different magnetic background field configurations are investigated in the framework of the solar flare energy build-up model of Van Tend and Kuperus (1978). The results are compared to characteristics of filaments and X-ray loops. Empirical flare predictors are interpreted theoretically.     

Resistive instabilities in coronal conditions

Resistive instabilities in a context referring to the solar corona are rigorously investigated. Various equilibrium configurations are considered, differing, among other things, by their behaviour with respect to fast, ideal instabilities. The computations presented cover in a unified scheme all known regimes of resistive modes and allow one to determine the fastest timescale over which resistivity can play a role. Comparisons with previous work as well as possible extensions are presented.     

Building the terrestrial planets: Constrained accretion in the inner Solar System

To date, no accretion model has succeeded in reproducing all observed constraints in the inner Solar System. These constraints include: (1) the orbits, in particular the small eccentricities, and (2) the masses of the terrestrial planets - Mars’ relatively small mass in particular has not been adequately reproduced in previous simulations; (3) the formation timescales of Earth and Mars, as interpreted from Hf/W isotopes; (4) the bulk structure of the asteroid belt, in particular the lack of an imprint of planetary embryo-sized objects; and (5) Earth’s relatively large water content, assuming that it was delivered in the form of water-rich primitive asteroidal material. Here we present results of 40 high-resolution (N = 1000-2000) dynamical simulations of late-stage planetary accretion with the goal of reproducing these constraints, although neglecting the planet Mercury. We assume that Jupiter and Saturn are fully-formed at the start of each simulation, and test orbital configurations that are both consistent with and contrary to the “Nice model”. We find that a configuration with Jupiter and Saturn on circular orbits forms low-eccentricity terrestrial planets and a water-rich Earth on the correct timescale, but Mars’ mass is too large by a factor of 5-10 and embryos are often stranded in the asteroid belt. A configuration with Jupiter and Saturn in their current locations but with slightly higher initial eccentricities (e = 0.07-0.1) produces a small Mars, an embryo-free asteroid belt, and a reasonable Earth analog but rarely allows water delivery to Earth. None of the configurations we tested reproduced all the observed constraints. Our simulations leave us with a problem: we can reasonably satisfy the observed constraints (except for Earth’s water) with a configuration of Jupiter and Saturn that is at best marginally consistent with models of the outer Solar System, as it does not allow for any outer planet migration after a few Myr. Alternately, giant planet configurations which are consistent with the Nice model fail to reproduce Mars’ small size.