Relativistic dispersion effects on SS433 jets

We here investigate the dispersion properties of radiation in the SS433 relativistic jets. We assume that the jet is composed of cold electron-proton plasma immersed in a predominantly parallel magnetic field to the jet axis. We find that for the mildly relativistic source SS433 (for which similar or equal to 79 degrees), the bulk velocity is too small (v similar or equal to 0.26c) to produce significant changes in the dispersion properties of the medium. Nonetheless, in the rarefied outer regions of the jets, where radio emission dominates, even a weak magnetic field has some influence on the dispersion properties and there appear two different electromagnetic branches that are slightly sensitive to the bulk relativistic motion. In the inner, X-ray region, the magnetic field is much stronger, but in this region the high electron density preserves the isotropic character of the local plasma and no branch separation occurs. In the region of the jet where the IR and optical emission dominates, the cold plasma may be also considered isotropic, i.e., neither the magnetic field nor the bulk velocity is able to affect the propagation of the radiation. Finally, we find that the Doppler line displacement in SS433 is affected by plasma dispersion only in a narrow frequency range in the far IR. As a consequence, although the shift (z) modulation due to precession of the SS433 jets is well described by previous work, it has to be corrected by plasma dispersion effects in the far-IR range.

SIMULATIONS OF WINDS OF WEAK-LINED T TAURI STARS. II. THE EFFECTS OF A TILTED MAGNETOSPHERE AND PLANETARY INTERACTIONS

Based on our previous work, we investigate here the effects on the wind and magnetospheric structures of weak-lined T Tauri stars due to a misalignment between the axis of rotation of the star and its magnetic dipole moment vector. In such a configuration, the system loses the axisymmetry presented in the aligned case, requiring a fully three-dimensional (3D) approach. We perform 3D numerical magnetohydrodynamic simulations of stellar winds and study the effects caused by different model parameters, namely the misalignment angle theta(t), the stellar period of rotation, the plasma-beta, and the heating index.. Our simulations take into account the interplay between the wind and the stellar magnetic field during the time evolution. The system reaches a periodic behavior with the same rotational period of the star. We show that the magnetic field lines present an oscillatory pattern. Furthermore, we obtain that by increasing theta(t), the wind velocity increases, especially in the case of strong magnetic field and relatively rapid stellar rotation. Our 3D, time-dependent wind models allow us to study the interaction of a magnetized wind with a magnetized extrasolar planet. Such interaction gives rise to reconnection, generating electrons that propagate along the planet`s magnetic field lines and produce electron cyclotron radiation at radio wavelengths. The power released in the interaction depends on the planet`s magnetic field intensity, its orbital radius, and on the stellar wind local characteristics. We find that a close-in Jupiter-like planet orbiting at 0.05 AU presents a radio power that is similar to 5 orders of magnitude larger than the one observed in Jupiter, which suggests that the stellar wind from a young star has the potential to generate strong planetary radio emission that could be detected in the near future with LOFAR. This radio power varies according to the phase of rotation of the star. For three selected simulations, we find a variation of the radio power of a factor 1.3-3.7, depending on theta(t). Moreover, we extend the investigation done in Vidotto et al. and analyze whether winds from misaligned stellar magnetospheres could cause a significant effect on planetary migration. Compared to the aligned case, we show that the timescale tau(w) for an appreciable radial motion of the planet is shorter for larger misalignment angles. While for the aligned case tau(w) similar or equal to 100 Myr, for a stellar magnetosphere tilted by theta(t) = 30 degrees, tau(w) ranges from similar to 40 to 70 Myr for a planet located at a radius of 0.05 AU. Further reduction on tau(w) might occur for even larger misalignment angles and/or different wind parameters.

SURFACE ALFVEN WAVE DAMPING IN A THREE-DIMENSIONAL SIMULATION OF THE SOLAR WIND

Here we investigate the contribution of surface Alfven wave damping to the heating of the solar wind in minima conditions. These waves are present in the regions of strong inhomogeneities in density or magnetic field (e.g., the border between open and closed magnetic field lines). Using a three-dimensional (3D) magnetohydrodynamics (MHD) model, we calculate the surface Alfven wave damping contribution between 1 and 4 R(circle dot) (solar radii), the region of interest for both acceleration and coronal heating. We consider waves with frequencies lower than those that are damped in the chromosphere and on the order of those dominating the heliosphere: 3 x 10(-6) to 10(-1) Hz. In the region between open and closed field lines, within a few R(circle dot) of the surface, no other major source of damping has been suggested for the low frequency waves we consider here. This work is the first to study surface Alfven waves in a 3D environment without assuming a priori a geometry of field lines or magnetic and density profiles. We demonstrate that projection effects from the plane of the sky to 3D are significant in the calculation of field line expansion. We determine that waves with frequencies >2.8 x 10(-4) Hz are damped between 1 and 4 R(circle dot). In quiet-Sun regions, surface Alfven waves are damped at further distances compared to active regions, thus carrying additional wave energy into the corona. We compare the surface Alfven wave contribution to the heating by a variable polytropic index and find it as an order of magnitude larger than needed for quiet-Sun regions. For active regions, the contribution to the heating is 20%. As it has been argued that a variable gamma acts as turbulence, our results indicate that surface Alfven wave damping is comparable to turbulence in the lower corona. This damping mechanism should be included self-consistently as an energy driver for the wind in global MHD models.

THREE-DIMENSIONAL NUMERICAL SIMULATIONS OF MAGNETIZED WINDS OF SOLAR-LIKE STARS

By means of self-consistent three-dimensional magnetohydrodynamics (MHD) numerical simulations, we analyze magnetized solar-like stellar winds and their dependence on the plasma-beta parameter (the ratio between thermal and magnetic energy densities). This is the first study to perform such analysis solving the fully ideal three-dimensional MHD equations. We adopt in our simulations a heating parameter described by gamma, which is responsible for the thermal acceleration of the wind. We analyze winds with polar magnetic field intensities ranging from 1 to 20 G. We show that the wind structure presents characteristics that are similar to the solar coronal wind. The steady-state magnetic field topology for all cases is similar, presenting a configuration of helmet streamer-type, with zones of closed field lines and open field lines coexisting. Higher magnetic field intensities lead to faster and hotter winds. For the maximum magnetic intensity simulated of 20 G and solar coronal base density, the wind velocity reaches values of similar to 1000 km s(-1) at r similar to 20r(0) and a maximum temperature of similar to 6 x 10(6) K at r similar to 6r(0). The increase of the field intensity generates a larger ""dead zone"" in the wind, i.e., the closed loops that inhibit matter to escape from latitudes lower than similar to 45 degrees extend farther away from the star. The Lorentz force leads naturally to a latitude-dependent wind. We show that by increasing the density and maintaining B(0) = 20 G the system recover back to slower and cooler winds. For a fixed gamma, we show that the key parameter in determining the wind velocity profile is the beta-parameter at the coronal base. Therefore, there is a group of magnetized flows that would present the same terminal velocity despite its thermal and magnetic energy densities, as long as the plasma-beta parameter is the same. This degeneracy, however, can be removed if we compare other physical parameters of the wind, such as the mass-loss rate. We analyze the influence of gamma in our results and we show that it is also important in determining the wind structure.

Local star formation triggered by supernova shocks in magnetized diffuse neutral clouds

In this work, considering the impact of a supernova remnant (SNR) with a neutral magnetized cloud we derived analytically a set of conditions that are favourable for driving gravitational instability in the cloud and thus star formation. Using these conditions, we have built diagrams of the SNR radius, R(SNR), versus the initial cloud density, n(c), that constrain a domain in the parameter space where star formation is allowed. This work is an extension to previous study performed without considering magnetic fields (Melioli et al. 2006, hereafter Paper I). The diagrams are also tested with fully three-dimensional MHD radiative cooling simulations involving a SNR and a self-gravitating cloud and we find that the numerical analysis is consistent with the results predicted by the diagrams. While the inclusion of a homogeneous magnetic field approximately perpendicular to the impact velocity of the SNR with an intensity similar to 1 mu G within the cloud results only a small shrinking of the star formation zone in the diagram relative to that without magnetic field, a larger magnetic field (similar to 10 mu G) causes a significant shrinking, as expected. Though derived from simple analytical considerations these diagrams provide a useful tool for identifying sites where star formation could be triggered by the impact of a supernova blast wave. Applications of them to a few regions of our own Galaxy (e.g. the large CO shell in the direction of Cassiopeia, and the Edge Cloud 2 in the direction of the Scorpious constellation) have revealed that star formation in those sites could have been triggered by shock waves from SNRs for specific values of the initial neutral cloud density and the SNR radius. Finally, we have evaluated the effective star formation efficiency for this sort of interaction and found that it is generally smaller than the observed values in our own Galaxy (SFE similar to 0.01-0.3). This result is consistent with previous work in the literature and also suggests that the mechanism presently investigated, though very powerful to drive structure formation, supersonic turbulence and eventually, local star formation, does not seem to be sufficient to drive global star formation in normal star-forming galaxies, not even when the magnetic field in the neutral clouds is neglected.

Paleointensity data from Early Cretaceous Ponta Grossa dikes (Brazil) using a multisample method

Definition of the long-term variation of the geomagnetic virtual dipole moment requires more reliable paleointensity results. Here, we applied a multisample protocol to the study of the 130.5 Ma Ponta Grossa basaltic dikes (southern Brazil) that carry a very stable dual-polarity magnetic component. The magnetic stability of the samples wits checked using thermomagnetic curves and by monitoring the magnetic Susceptibility evolution through the paleointensity experiments. Twelve sites containing the least alterable samples were chosen for the paleointensity measurements. Although these rocks failed stepwise double-heating experiments, they yielded coherent results in the multisample method for all sites but one. The coherent sites show low to moderate field intensities between 5.7 +/- 0.2 and 26.4 +/- 0.7 mu T (average 13.4 +/- 1.9 mu T). Virtual dipole moments for these sites range from 1.3 +/- 0.04 to 6.0 +/- 0.2 x 10(22) A m(2) (average 2.9 +/- 0.5 x 10(22) A m(2)). Our results agree with the tendency for low dipole moments during the Early Cretaceous, immediately prior to the Cretaceous Normal Superchron (CNS). The available paleointensity database shows a strong variability of the field between 80 and 160 Ma. There seems to be no firm evidence for a Mesozoic Dipole Low, but a long-term tendency does emerge from the data with the highest dipole moments Occurring at the middle of the CNS.

Fractal structures in nonlinear plasma physics

Fractal structures appear in many situations related to the dynamics of conservative as well as dissipative dynamical systems, being a manifestation of chaotic behaviour. In open area-preserving discrete dynamical systems we can find fractal structures in the form of fractal boundaries, associated to escape basins, and even possessing the more general property of Wada. Such systems appear in certain applications in plasma physics, like the magnetic field line behaviour in tokamaks with ergodic limiters. The main purpose of this paper is to show how such fractal structures have observable consequences in terms of the transport properties in the plasma edge of tokamaks, some of which have been experimentally verified. We emphasize the role of the fractal structures in the understanding of mesoscale phenomena in plasmas, such as electromagnetic turbulence.

The non-twist standard map with robust tori

The non-twist standard map occurs frequently in many fields of science specially in modelling the dynamics of the magnetic field lines in tokamaks. Robust tori, dynamical barriers that impede the radial transport among different regions of the phase space, are introduced in the non-twist standard map in a conservative fashion. The resulting non-twist standard map with robust tori is an improved model to study transport barriers in plasmas confined in tokamaks.

Integrable maps with non-trivial topology: application to divertor configurations

We explore a method for constructing two-dimensional area-preserving, integrable maps associated with Hamiltonian systems, with a given set of fixed points and given invariant curves. The method is used to find an integrable Poincare map for the field lines in a large aspect ratio tokamak with a poloidal single-null divertor. The divertor field is a superposition of a magnetohydrodynamic equilibrium with an arbitrarily chosen safety factor profile, with a wire carrying an electric current to create an X-point. This integrable map is perturbed by an impulsive perturbation that describes non-axisymmetric magnetic resonances at the plasma edge. The non-integrable perturbed map is applied to study the structure of the open field lines in the scrape-off layer, reproducing the main transport features obtained by integrating numerically the magnetic field line equations, such as the connection lengths and magnetic footprints on the divertor plate.

High-frequency extensions of magnetorotational instability in astrophysical plasmas

High-frequency extensions of magnetorotational instability driven by the Velikhov effect beyond the standard magnetohydrodynamic (MHD) regime are studied. The existence of the well-known Hall regime and a new electron inertia regime is demonstrated. The electron inertia regime is realized for a lesser plasma magnetization of rotating plasma than that in the Hall regime. It includes the subregime of nonmagnetized electrons. It is shown that, in contrast to the standard MHD regime and the Hall regime, magnetorotational instability in this subregime can be driven only at positive values of dln Omega/dlnr, where Omega is the plasma rotation frequency and r is the radial coordinate. The permittivity of rotating plasma beyond the standard MHD regime, including both the Hall regime and the electron inertia regime, is calculated.

Ideal internal kink modes in a differentially rotating cylindrical plasma

The Velikhov effect leading to magnetorotational instability (MRI) is incorporated into the theory of ideal internal kink modes in a differentially rotating cylindrical plasma column. It is shown that this effect can play a stabilizing role for suitably organized plasma rotation profiles, leading to suppression of MHD (magnetohydrodynamic) instabilities in magnetic confinement systems. The role of this effect in the problem of the Suydam and the m = 1 internal kink modes is elucidated, where m is the poloidal mode number.

Chaotic transport in reversed shear tokamaks

For tokamak models using simplified geometries and reversed shear plasma profiles, we have numerically investigated how the onset of Lagrangian chaos at the plasma edge may affect the plasma confinement in two distinct but closely related problems. Firstly, we have considered the motion of particles in drift waves in the presence of an equilibrium radial electric field with shear. We have shown that the radial particle transport caused by this motion is selective in phase space, being determined by the resonant drift waves and depending on the parameters of both the resonant waves and the electric field profile. Moreover, we have shown that an additional transport barrier may be created at the plasma edge by increasing the electric field. In the second place, we have studied escape patterns and magnetic footprints of chaotic magnetic field lines in the region near a tokamak wall, when there are resonant modes due to the action of an ergodic magnetic limiter. A non-monotonic safety factor profile has been used in the analysis of field line topology in a region of negative magnetic shear. We have observed that, if internal modes are perturbed, the distributions of field line connection lengths and magnetic footprints exhibit spatially localized escape channels. For typical physical parameters of a fusion plasma, the two Lagrangian chaotic processes considered in this work can be effective in usual conditions so as to influence plasma confinement. The reversed shear effects discussed in this work may also contribute to evaluate the transport barrier relevance in advanced confinement scenarios in future tokamak experiments.

Contributions to the theory of magnetorotational instability and waves in a rotating plasma

The one-fluid magnetohydrodynamic (MHD) theory of magnetorotational instability (MRI) in an ideal plasma is presented. The theory predicts the possibility of MRI for arbitrary 0, where 0 is the ratio of the plasma pressure to the magnetic field pressure. The kinetic theory of MRI in a collisionless plasma is developed. It is demonstrated that as in the ideal MHD, MRI can occur in such a plasma for arbitrary P. The mechanism of MRI is discussed; it is shown that the instability appears because of a perturbed parallel electric field. The electrodynamic description of MRI is formulated under the assumption that the dispersion relation is expressed in terms of the permittivity tensor; general properties of this tensor are analyzed. It is shown to be separated into the nonrotational and rotational parts. With this in mind, the first step for incorporation of MRI into the general theory of plasma instabilities is taken. The rotation effects on Alfven waves are considered.

Plasma confinement in tokamaks with robust torus

We present a non-linear symplectic map that describes the alterations of the magnetic field lines inside the tokamak plasma due to the presence of a robust torus (RT) at the plasma edge. This RT prevents the magnetic field lines from reaching the tokamak wall and reduces, in its vicinity, the islands and invariant curve destruction due to resonant perturbations. The map describes the equilibrium magnetic field lines perturbed by resonances created by ergodic magnetic limiters (EMLs). We present the results obtained for twist and non-twist mappings derived for monotonic and non-monotonic plasma current density radial profiles, respectively. Our results indicate that the RT implementation would decrease the field line transport at the tokamak plasma edge. (C) 2010 Elsevier B.V. All rights reserved.

Observation in the MINOS far detector of the shadowing of cosmic rays by the sun and moon

The shadowing of cosmic ray primaries by the moon and sun was observed by the MINOS far detector at a depth of 2070 mwe using 83.54 million cosmic ray muons accumulated over 1857.91 live-days. The shadow of the moon was detected at the 5.6 sigma level and the shadow of the sun at the 3.8 sigma level using a log-likelihood search in celestial coordinates. The moon shadow was used to quantify the absolute astrophysical pointing of the detector to be 0.17 +/- 0.12 degrees. Hints of interplanetary magnetic field effects were observed in both the sun and moon shadow. Published by Elsevier B.V.