The non-linear evolution of magnetic flux ropes: 3. effects of dissipation

International audienceWe study the evolution (expansion or oscillation) of cylindrically symmetric magnetic flux ropes when the energy dissipation is due to a drag force proportional to the product of the plasma density and the radial speed of expansion. The problem is reduced to a single, second-order, ordinary differential equation for a damped, non-linear oscillator. Motivated by recent work on the interplanetary medium and the solar corona, we consider polytropes whose index, ?, may be less than unity. Numerical analysis shows that, in contrast to the small-amplitude case, large-amplitude oscillations are quasi-periodic with frequencies substantially higher than those of undamped oscillators. The asymptotic behaviour described by the momentum equation is determined by a balance between the drag force and the gradient of the gas pressure, leading to a velocity of expansion of the flux rope which may be expressed as (1/2?)r/t, where r is the radial coordinate and t is the time. In the absence of a drag force, we found in earlier work that the evolution depends both on the polytropic index and on a dimensionless parameter, ?. Parameter ? was found to have a critical value above which oscillations are impossible, and below which they can exist only for energies less than a certain energy threshold. In the presence of a drag force, the concept of a critical ? remains valid, and when ? is above critical, the oscillatory mode disappears altogether. Furthermore, critical ? remains dependent only on ? and is, in particular, independent of the normalized drag coefficient, ?*. Below critical ?, however, the energy required for the flux rope to escape to infinity depends not only on ? (as in the conservative force case) but also on ?*. This work indicates how under certain conditions a small change in the viscous drag coefficient or the initial energy may alter the evolution drastically. It is thus important to determine ?* and ? from observations

Kilohertz QPOs in Neutron Star Binaries modeled as Keplerian Oscillations in a Rotating Frame of Reference

Since the discovery of kHz quasi-periodic oscillations (QPO) in neutron star binaries, the difference between peak frequencies of two modes in the upper part of the spectrum, i.e. Delta (omega)=omega_h-omega_K has been studied extensively. The idea that the difference Delta(omega) is constant and (as a beat frequency) is related to the rotational frequency of the neutron star has been tested previously. The observed decrease of Delta(omega) when omega_h and omega_k increase has weakened the beat frequency interpretation. We put forward a different paradigm: a Keplerian oscillator under the influence of the Coriolis force. For such an oscillator, omega_h and the assumed Keplerian frequency omega_k hold an upper hybrid frequency relation: omega^2_h-omega^2_K=4*Omega^2, where Omega is the rotational frequency of the star's magnetosphere near the equatorial plane. For three sources (Sco X-1, 4U 1608-52 and 4U 1702-429), we demonstrate that the solid body rotation Omega=Omega_0=const. is a good first order approximation. Within the second order approximation, the slow variation of Omega as a function of omega_K reveals the structure of the magnetospheric differential rotation. For Sco X-1, the QPO have frequencies approximately 45 and 90 Hz which we interpret as the 1st and 2nd harmonics of the lower branch of the Keplerian oscillations for the rotator with vector Omega not aligned with the normal of the disk: omega_L/2 pi=(Omega/pi)(omega_K/omega_h)sin(delta) where delta is the angle between vector Omega and the vector normal to the disk.Comment: 13 pages, 3 figures, accepted for publications in ApJ Letter

Experimental and theoretical lifetimes and transition probabilities in Sb I

We present experimental atomic lifetimes for 12 levels in Sb I, out of which seven are reported for the first time. The levels belong to the 5p$^2$($^3$P)6s $^{2}$P, $^{4}$P and 5p$^2$($^3$P)5d $^{4}$P, $^{4}$F and $^{2}$F terms. The lifetimes were measured using time-resolved laser-induced fluorescence. In addition, we report new calculations of transition probabilities in Sb I using a Multiconfigurational Dirac-Hartree-Fock method. The physical model being tested through comparisons between theoretical and experimental lifetimes for 5d and 6s levels. The lifetimes of the 5d $^4$F$_{3/2, 5/2, 7/2}$ levels (19.5, 7.8 and 54 ns, respectively) depend strongly on the $J$-value. This is explained by different degrees of level mixing for the different levels in the $^4$F term.Comment: 10 page

On The Low Frequency Quasi Periodic Oscillations of X-ray Sources

Based on the interpretation of the twin kilohertz Quasi Periodic Oscillations (kHz QPOs) of X-ray spectra of Low Mass X-Ray Binaries (LMXBs) to the Keplerian and the periastron precession frequencies at the magnetosphere-disk of X-ray neutron star (NS) respectively, we ascribe the low frequency Quasi Periodic Oscillations (LFQPO) and HBO (15-60 Hz QPO for Z sources or Atoll sources) to the periastron precession at some outer disk radius. The obtained conclusions include: all QPO frequencies increase with increasing the accretion rate. The obtained theoretical relations between HBO (LFQPO) frequency and the kHz QPO frequency are similar to the measured empirical formula. Further, the possible dynamical mechanism for QPO production is discussed.Comment: 6 pages, 2 figures, accepted by APSS, 200

Modeling the Subsurface Structure of Sunspots

While sunspots are easily observed at the solar surface, determining their subsurface structure is not trivial. There are two main hypotheses for the subsurface structure of sunspots: the monolithic model and the cluster model. Local helioseismology is the only means by which we can investigate subphotospheric structure. However, as current linear inversion techniques do not yet allow helioseismology to probe the internal structure with sufficient confidence to distinguish between the monolith and cluster models, the development of physically realistic sunspot models are a priority for helioseismologists. This is because they are not only important indicators of the variety of physical effects that may influence helioseismic inferences in active regions, but they also enable detailed assessments of the validity of helioseismic interpretations through numerical forward modeling. In this paper, we provide a critical review of the existing sunspot models and an overview of numerical methods employed to model wave propagation through model sunspots. We then carry out an helioseismic analysis of the sunspot in Active Region 9787 and address the serious inconsistencies uncovered by \citeauthor{gizonetal2009}~(\citeyear{gizonetal2009,gizonetal2009a}). We find that this sunspot is most probably associated with a shallow, positive wave-speed perturbation (unlike the traditional two-layer model) and that travel-time measurements are consistent with a horizontal outflow in the surrounding moat.Comment: 73 pages, 19 figures, accepted by Solar Physic

Molecular Interpretation of ACTH-β-Endorphin Coaggregation: Relevance to Secretory Granule Biogenesis

Peptide/protein hormones could be stored as non-toxic amyloid-like structures in pituitary secretory granules. ACTH and β-endorphin are two of the important peptide hormones that get co-stored in the pituitary secretory granules. Here, we study molecular interactions between ACTH and β-endorphin and their colocalization in the form of amyloid aggregates. Although ACTH is known to be a part of ACTH-β-endorphin aggregate, ACTH alone cannot aggregate into amyloid under various plausible conditions. Using all atom molecular dynamics simulation we investigate the early molecular interaction events in the ACTH-β-endorphin system, β-endorphin-only system and ACTH-only system. We find that β-endorphin and ACTH formed an interacting unit, whereas negligible interactions were observed between ACTH molecules in ACTH-only system. Our data suggest that ACTH is not only involved in interaction with β-endorphin but also enhances the stability of mixed oligomers of the entire system