A heliospheric hybrid model: hydrodynamic plasma flow and kinetic cosmic ray transport

International audienceIn this paper we present a new five particle species hybrid model for calculating cosmic ray particle transport and acceleration in a dynamic heliospheric environment. In particular the effects of solar cycle related changes in the solar wind speed on the heliospheric geometry, solar wind flow and cosmic ray distribution are discussed, when a polar-ecliptic asymmetry at the inner boundary is modeled. It is shown that the disappearance of the fast solar wind over the solar poles toward solar maximum influences the geometry of the termination shock which is an important structure for cosmic ray acceleration. For solar maximum conditions, the shock radius is smaller in the polar regions and in the heliospheric tail compared to solar minimum. These changes influence cosmic ray transport and acceleration in these regions, especially for the polarity cycle where positive particles drift in along the heliospheric current sheet. For this polarity cycle, and for both the anomalous and galactic cosmic ray protons, an increase in particle intensities at the shock in the heliospheric tail is computed as the shock moves inward toward the Sun. For the heliospheric nose, it is also shown that both the plasma speed and cosmic ray intensities are relative insensitive to changes in the latitudinal profile of the solar wind speed. Therefore toward solar maximum conditions there is a decrease in the nose-tail asymmetry of the computed cosmic ray distribution compared to solar minimum conditions

Unidentified Galactic High-Energy Sources as Ancient Pulsar Wind Nebulae in the light of new high energy observations and the new code

In a Pulsar Wind Nebula (PWN), the lifetime of inverse Compton (IC) emitting electrons exceeds the lifetime of its progenitor pulsar (as well as its shell-type remnant), but it also exceeds the age of those that emit via synchrotron radiation. Therefore, during its evolution, the PWN can remain bright in IC so that its GeV-TeV gamma-ray flux remains high for timescales much larger (for 10^5 - 10^6 yrs) than the pulsar lifetime and the X-ray PWN lifetime. In this scenario, the magnetic field in the cavity induced by the wind of the progenitor star plays a crucial role. This scenario is in line with the discovery of several unidentified or "dark" sources in the TeV gamma-ray band without X-ray counterparts; and it is also finding confirmation in the recent discoveries at GeV gamma rays. Moreover, these consequences could be also important for reinterpreting the detection of starburst galaxies in the TeV gamma-ray band when considering a leptonic origin of the gamma-ray signal. Both theoretical aspects and their observational proofs will be discussed, as well as the first results of our new modeling code.Comment: Proceedings of the 5th International Symposium on High-Energy Gamma-Ray Astronomy (Gamma2012

Non-collinear coupling between magnetic adatoms in carbon nanotubes

The long range character of the exchange coupling between localized magnetic moments indirectly mediated by the conduction electrons of metallic hosts often plays a significant role in determining the magnetic order of low-dimensional structures. In addition to this indirect coupling, here we show that the direct exchange interaction that arises when the moments are not too far apart may induce a non-collinear magnetic order that cannot be characterized by a Heisenberg-like interaction between the magnetic moments. We argue that this effect can be manipulated to control the magnetization alignment of magnetic dimers adsorbed to the walls of carbon nanotubes.Comment: 13 pages, 5 figures, submitted to PR

Outflows and Jets from Collapsing Magnetized Cloud Cores

Star formation is usually accompanied by outflow phenomena. There is strong evidence that these outflows and jets are launched from the protostellar disk by magneto-rotational processes. Here, we report on our three dimensional, adaptive mesh, magneto-hydrodynamic simulations of collapsing, rotating, magnetized Bonnor-Ebert-Spheres whose properties are taken directly from observations. In contrast to the pure hydro case where no outflows are seen, our present simulations show an outflow from the protodisk surface at ~ AU and a jet at ~ 0.07 AU after a strong toroidal magnetic field build up. The large scale outflow, which extends up to ~ AU at the end of our simulation, is driven by toroidal magnetic pressure (spring), whereas the jet is powered by magneto-centrifugal force (fling). At the final stage of our simulation these winds are still confined within two respective shock fronts. Furthermore, we find that the jet-wind and the disk-anchored magnetic field extracts a considerable amount of angular momentum from the protostellar disk. The initial spin of our cloud core was chosen high enough to produce a binary system. We indeed find a close binary system (separation ~ 3 R_sol) which results from the fragmentation of an earlier formed ring structure. The magnetic field strength in these protostars reaches ~ 3 kG and becomes about 3 G at 1 AU from the center in agreement with recent observational results.Comment: revised version, accepted for publication in ApJ, a higher resolution version of this paper can be downloaded at http://www.physics.mcmaster.ca/~banerjee/outflows.pd

Two-component jet simulations: I. Topological stability of analytical MHD outflow solutions

Observations of collimated outflows in young stellar objects indicate that several features of the jets can be understood by adopting the picture of a two-component outflow, wherein a central stellar component around the jet axis is surrounded by an extended disk-wind. The precise contribution of each component may depend on the intrinsic physical properties of the YSO-disk system as well as its evolutionary stage. In this context, the present article starts a systematic investigation of two-component jet models via time-dependent simulations of two prototypical and complementary analytical solutions, each closely related to the properties of stellar-outflows and disk-winds. These models describe a meridionally and a radially self-similar exact solution of the steady-state, ideal hydromagnetic equations, respectively. By using the PLUTO code to carry out the simulations, the study focuses on the topological stability of each of the two analytical solutions, which are successfully extended to all space by removing their singularities. In addition, their behavior and robustness over several physical and numerical modifications is extensively examined. It is found that radially self-similar solutions (disk-winds) always reach a final steady-state while maintaining all their well-defined properties. The different ways to replace the singular part of the solution around the symmetry axis, being a first approximation towards a two-component outflow, lead to the appearance of a shock at the super-fast domain corresponding to the fast magnetosonic separatrix surface. Conversely, the asymptotic configuration and the stability of meridionally self-similar models (stellar-winds) is related to the heating processes at the base of the wind.Comment: Accepted for publication in A&

FEATURE SELECTION APPLIED TO THE TIME-FREQUENCY REPRESENTATION OF MUSCLE NEAR-INFRARED SPECTROSCOPY (NIRS) SIGNALS: CHARACTERIZATION OF DIABETIC OXYGENATION PATTERNS

Diabetic patients might present peripheral microcirculation impairment and might benefit from physical training. Thirty-nine diabetic patients underwent the monitoring of the tibialis anterior muscle oxygenation during a series of voluntary ankle flexo-extensions by near-infrared spectroscopy (NIRS). NIRS signals were acquired before and after training protocols. Sixteen control subjects were tested with the same protocol. Time-frequency distributions of the Cohen's class were used to process the NIRS signals relative to the concentration changes of oxygenated and reduced hemoglobin. A total of 24 variables were measured for each subject and the most discriminative were selected by using four feature selection algorithms: QuickReduct, Genetic Rough-Set Attribute Reduction, Ant Rough-Set Attribute Reduction, and traditional ANOVA. Artificial neural networks were used to validate the discriminative power of the selected features. Results showed that different algorithms extracted different sets of variables, but all the combinations were discriminative. The best classification accuracy was about 70%. The oxygenation variables were selected when comparing controls to diabetic patients or diabetic patients before and after training. This preliminary study showed the importance of feature selection techniques in NIRS assessment of diabetic peripheral vascular impairmen

Transition to complete synchronization in phase coupled oscillators with nearest neighbours coupling

We investigate synchronization in a Kuramoto-like model with nearest neighbour coupling. Upon analyzing the behaviour of individual oscillators at the onset of complete synchronization, we show that the time interval between bursts in the time dependence of the frequencies of the oscillators exhibits universal scaling and blows up at the critical coupling strength. We also bring out a key mechanism that leads to phase locking. Finally, we deduce forms for the phases and frequencies at the onset of complete synchronization.Comment: 6 pages, 4 figures, to appear in CHAO