Turbulent magnetic pumping in a Babcock-Leighton solar dynamo model

Context. The turbulent pumping effect corresponds to the transport of magnetic flux due to the presence of density and turbulence gradients in convectively unstable layers. In the induction equation it appears as an advective term and for this reason it is expected to be important in the solar and stellar dynamo processes. Aims. We explore the effects of turbulent pumping in a flux-dominated Babcock-Leighton solar dynamo model with a solar-like rotation law. Methods. As a first step, only vertical pumping has been considered through the inclusion of a radial diamagnetic term in the induction equation. In the second step, a latitudinal pumping term was included and then, a near-surface shear was included. Results. The results reveal the importance of the pumping mechanism in solving current limitations in mean field dynamo modeling, such as the storage of the magnetic flux and the latitudinal distribution of the sunspots. If a meridional flow is assumed to be present only in the upper part of the convective zone, it is the full turbulent pumping that regulates both the period of the solar cycle and the latitudinal distribution of the sunspot activity. In models that consider shear near the surface, a second shell of toroidal field is generated above r = 0.95 R(circle dot) at all latitudes. If the full pumping is also included, the polar toroidal fields are efficiently advected inwards, and the toroidal magnetic activity survives only at the observed latitudes near the equator. With regard to the parity of the magnetic field, only models that combine turbulent pumping with near-surface shear always converge to the dipolar parity. Conclusions. This result suggests that, under the Babcock-Leighton approach, the equartorward motion of the observed magnetic activity is governed by the latitudinal pumping of the toroidal magnetic field rather than by a large scale coherent meridional flow. Our results support the idea that the parity problem is related to the quadrupolar imprint of the meridional flow on the poloidal component of the magnetic field and the turbulent pumping positively contributes to wash out this imprint.

CONTRASTING PHYLOGEOGRAPHIC PATTERNS IN MITOCHONDRIAL DNA AND MICROSATELLITES: EVIDENCE OF FEMALE PHILOPATRY AND MALE-BIASED GENE FLOW AMONG REGIONAL POPULATIONS OF THE BLUE-AND-YELLOW MACAW (PSITTACIFORMES: ARA ARARAUNA) IN BRAZIL

Comparing the patterns of population differentiation among genetic markers with different modes of inheritance call provide insights into patterns of sex-biased dispersal and gene flow. The blue-and-yellow Macaw (Ara ararauna) is a Neotropical parrot with a broad geographic distribution ill South America. However, little is known about the natural history and current status Of remaining wild populations, including levels of genetic variability. The progressive decline and possible fragmentation of populations may endanger this species in the near future. We analyzed mitochondrial DNA (mtDNA) control-region sequences and six microsatellite 106 Of Blue-and-yellow Macaws sampled throughout their geographic range ill Brazil to describe population genetic Structure, to make inferences about historical demography and dispersal behavior, and to provide insight for conservation efforts. Analyses of population genetic structure based on mtDNA showed evidence of two major populations ill western and eastern Brazil that share a few low-frequency haplotypes. This phylogeographic pattern seems to have originated by the historical isolation of Blue-and-yellow Macaw populations similar to 374,000 years ago and has been maintained by restricted gene flow and female philopatry. By contrast, variation ill biparentally inherited microsatellites was not structured geographically, Male-biased dispersal and female philopatry best explain the different patterns observed in these two markers. Because females disperse less than males, the two regional populations with well-differentiated mtDNA haplogroups should be considered two different management units for conservation purposes. Received 4 November 2007 accepted 10 December 2008.

PARTITIONED ANALYSIS FOR DIMENSIONALLY-HETEROGENEOUS HYDRAULIC NETWORKS

In this work an iterative strategy is developed to tackle the problem of coupling dimensionally-heterogeneous models in the context of fluid mechanics. The procedure proposed here makes use of a reinterpretation of the original problem as a nonlinear interface problem for which classical nonlinear solvers can be applied. Strong coupling of the partitions is achieved while dealing with different codes for each partition, each code in black-box mode. The main application for which this procedure is envisaged arises when modeling hydraulic networks in which complex and simple subsystems are treated using detailed and simplified models, correspondingly. The potentialities and the performance of the strategy are assessed through several examples involving transient flows and complex network configurations.

Escape patterns of chaotic magnetic field lines in a tokamak with reversed magnetic shear and an ergodic limiter

The existence of a reversed magnetic shear in tokamaks improves the plasma confinement through the formation of internal transport barriers that reduce radial particle and heat transport. However, the transport poloidal profile is much influenced by the presence of chaotic magnetic field lines at the plasma edge caused by external perturbations. Contrary to many expectations, it has been observed that such a chaotic region does not uniformize heat and particle deposition on the inner tokamak wall. The deposition is characterized instead by structured patterns called magnetic footprints, here investigated for a nonmonotonic analytical plasma equilibrium perturbed by an ergodic limiter. The magnetic footprints appear due to the underlying mathematical skeleton of chaotic magnetic field lines determined by the manifold tangles. For the investigated edge safety factor ranges, these effects on the wall are associated with the field line stickiness and escape channels due to internal island chains near the flux surfaces. Comparisons between magnetic footprints and escape basins from different equilibrium and ergodic limiter characteristic parameters show that highly concentrated magnetic footprints can be avoided by properly choosing these parameters. (c) 2008 American Institute of Physics.

Nonaxisymmetric magnetorotational instability in ideal and viscous plasmas

The excitation of magnetorotational instability (MRI) in rotating laboratory plasmas is investigated. In contrast to astrophysical plasmas, in which gravitation plays an important role, in laboratory plasmas it can be neglected and the plasma rotation is equilibrated by the pressure gradient. The analysis is restricted to the simple model of a magnetic confinement configuration with cylindrical symmetry, in which nonaxisymmetric perturbations are investigated using the local approximation. Starting from the simplest case of an ideal plasma, the corresponding dispersion relations are derived for more complicated models including the physical effects of parallel and perpendicular viscosities. The Friemann-Rotenberg approach used for ideal plasmas is generalized for the viscous model and an analytical expression for the instability boundary is obtained. It is shown that, in addition to the standard effect of radial derivative of the rotation frequency (the Velikhov effect), which can be destabilizing or stabilizing depending on the sign of this derivative in the ideal plasma, there is a destabilizing effect proportional to the fourth power of the rotation frequency, or, what is the same, to the square of the plasma pressure gradient, and to the square of the azimuthal mode number of the perturbations. It is shown that the instability boundary also depends on the product of the plasma pressure and density gradients, which has a destabilizing effect when it is negative. In the case of parallel viscosity, the MRI looks like an ideal instability independent of viscosity, while, in the case of strong perpendicular viscosity, it is a dissipative instability with the growth rate inversely proportional to the characteristic viscous decay rate. We point out, however, that the modes of the continuous range of the magnetohydrodynamics spectrum are not taken into account in this paper, and they can be more dangerous than those that are considered. (c) 2008 American Institute of Physics.

Free-standing urethane/urea elastomer films undoped and doped with ferro-nano-particles

We report on an experimental study of the structures presented by urethane/urea elastomeric films without and with ferromagnetic nanoparticles incorporated. The study is made by using the X-ray diffraction, nuclear magnetic resonance (NMR), optical, atomic and magnetic force (MFM) microscopy techniques, and mechanical assays. The structure of the elastomeric matrix is characterized by a distance of 0.46 nm between neighboring molecular segments, almost independent on the stretching applied. The shear casting performed in order to obtain the elastomeric films tends to orient the molecules parallel to the flow direction thus introducing anisotropy in the molecular network which is reflected on the values obtained for the orientational order parameter and its increase for the stretched films. In the case of nanoparticles-doped samples, the structure remains nearly unchanged although the local order parameter is clearly larger for the undoped films. NMR experiments evidence modifications in the molecular network local ordering. Micrometer size clusters were observed by MFM for even small concentration of magnetic particles.

Two-dimensional supersonic nonlinear Schrodinger flow past an extended obstacle

Supersonic flow of a superfluid past a slender impenetrable macroscopic obstacle is studied in the framework of the two-dimensional (2D) defocusing nonlinear Schroumldinger (NLS) equation. This problem is of fundamental importance as a dispersive analog of the corresponding classical gas-dynamics problem. Assuming the oncoming flow speed is sufficiently high, we asymptotically reduce the original boundary-value problem for a steady flow past a slender body to the one-dimensional dispersive piston problem described by the nonstationary NLS equation, in which the role of time is played by the stretched x coordinate and the piston motion curve is defined by the spatial body profile. Two steady oblique spatial dispersive shock waves (DSWs) spreading from the pointed ends of the body are generated in both half planes. These are described analytically by constructing appropriate exact solutions of the Whitham modulation equations for the front DSW and by using a generalized Bohr-Sommerfeld quantization rule for the oblique dark soliton fan in the rear DSW. We propose an extension of the traditional modulation description of DSWs to include the linear ""ship-wave"" pattern forming outside the nonlinear modulation region of the front DSW. Our analytic results are supported by direct 2D unsteady numerical simulations and are relevant to recent experiments on Bose-Einstein condensates freely expanding past obstacles.

Elliptic and Hexadecapole Flow of Charged Hadrons in Au plus Au Collisions at root s(NN)=200 GeV

Differential measurements of the elliptic (upsilon(2)) and hexadecapole (upsilon(4)) Fourier flow coefficients are reported for charged hadrons as a function of transverse momentum (p(T)) and collision centrality or number of participant nucleons (N(part)) for Au + Au collisions at root s(NN) = 200 GeV/ The upsilon(2,4) measurements at pseudorapidity vertical bar eta vertical bar <= 0.35, obtained with four separate reaction-plane detectors positioned in the range 1.0 < vertical bar eta vertical bar < 3.9, show good agreement, indicating the absence of significant Delta eta-dependent nonflow correlations. Sizable values for upsilon(4)(p(T)) are observed with a ratio upsilon(4)(p(T), N(part))/upsilon(2)(2)(p(T), N(part)) approximate to 0.8 for 50 less than or similar to N(part) less than or similar to 200, which is compatible with the combined effects of a finite viscosity and initial eccentricity fluctuations. For N(part) greater than or similar to 200 this ratio increases up to 1.7 in the most central collisions.

Systematic studies of elliptic flow measurements in Au plus Au collisions at s(NN)=200 GeV

We present inclusive charged hadron elliptic flow (v(2)) measured over the pseudorapidity range vertical bar eta vertical bar < 0.35 in Au+Au collisions at s(NN)=200 GeV. Results for v(2) are presented over a broad range of transverse momentum (p(T)=0.2-8.0 GeV/c) and centrality (0-60%). To study nonflow effects that are correlations other than collective flow, as well as the fluctuations of v(2), we compare two different analysis methods: (1) the event-plane method from two independent subdetectors at forward (vertical bar eta vertical bar=3.1-3.9) and beam (vertical bar eta vertical bar>6.5) pseudorapidities and (2) the two-particle cumulant method extracted using correlations between particles detected at midrapidity. The two event-plane results are consistent within systematic uncertainties over the measured p(T) and in centrality 0-40%. There is at most a 20% difference in the v(2) between the two event-plane methods in peripheral (40-60%) collisions. The comparisons between the two-particle cumulant results and the standard event-plane measurements are discussed.

Directed flow at midrapidity in event-by-event hydrodynamics

Fluctuations in the initial geometry of a nucleus-nucleus collision have been recently shown to result in a new type of directed flow (v(1)) that, unlike the usual directed flow, is also present at midrapidity. We compute this new v(1) versus transverse momentum and centrality for Au-Au collisions at RHIC using the hydrodynamic code NeXSPheRIO. We find that the event plane of v(1) is correlated with the angle of the initial dipole of the distribution, as predicted, though with a large dispersion. It is uncorrelated with the reaction plane. Our results are in excellent agreement with results inferred from STAR correlation data.

Importance of granular structure in the initial conditions for the elliptic flow

We show the effects of the granular structure of the initial conditions of a hydrodynamic description of high-energy nucleus-nucleus collisions on some observables, especially on the elliptic-flow parameter upsilon(2). Such a structure enhances production of isotropically distributed high-p(T) particles, making upsilon(2) smaller there. Also, it reduces upsilon(2) in the forward and backward regions where the global matter density is smaller and, therefore, where such effects become more efficacious.

Entanglement and classical instabilities: Fingerprints of electron-hole-to-exciton phase transition in a simple model

We propose a schematic model to study the formation of excitons in bilayer electron systems. The phase transition is signalized both in the quantum and classical versions of the model. In the present contribution we show that not only the quantum ground state but also higher energy states, up to the energy of the corresponding classical separatrix orbit, ""sense"" the transition. We also show two types of one-to-one correspondences in this system: On the one hand, between the changes in the degree of entanglement for these low-lying quantum states and the changes in the density of energy levels; on the other hand, between the variation in the expected number of excitons for a given quantum state and the behavior of the corresponding classical orbit.

Simple analytic model for astrophysical S factors

We propose a physically transparent analytic model of astrophysical S factors as a function of a center-of-mass energy E of colliding nuclei (below and above the Coulomb barrier) for nonresonant fusion reactions. For any given reaction, the S(E) model contains four parameters [two of which approximate the barrier potential, U(r)]. They are easily interpolated along many reactions involving isotopes of the same elements; they give accurate practical expressions for S(E) with only several input parameters for many reactions. The model reproduces the suppression of S(E) at low energies (of astrophysical importance) due to the shape of the low-r wing of U(r). The model can be used to reconstruct U(r) from computed or measured S(E). For illustration, we parametrize our recent calculations of S(E) (using the Sao Paulo potential and the barrier penetration formalism) for 946 reactions involving stable and unstable isotopes of C, O, Ne, and Mg (with nine parameters for all reactions involving many isotopes of the same elements, e. g., C+O). In addition, we analyze astrophysically important (12)C+(12)C reaction, compare theoretical models with experimental data, and discuss the problem of interpolating reliably known S(E) values to low energies (E less than or similar to 2-3 MeV).

Charged and strange hadron elliptic flow in Cu plus Cu collisions at root s(NN)=62.4 and 200 GeV

We present the results of an elliptic flow, v(2), analysis of Cu + Cu collisions recorded with the solenoidal tracker detector (STAR) at the BNL Relativistic Heavy Ion Collider at root s(NN) = 62.4 and 200 GeV. Elliptic flow as a function of transverse momentum, v(2)(p(T)), is reported for different collision centralities for charged hadrons h(+/-) and strangeness-ontaining hadrons K(S)(0), Lambda, Xi, and phi in the midrapidity region vertical bar eta vertical bar < 1.0. Significant reduction in systematic uncertainty of the measurement due to nonflow effects has been achieved by correlating particles at midrapidity, vertical bar eta vertical bar < 1.0, with those at forward rapidity, 2.5 < vertical bar eta vertical bar < 4.0. We also present azimuthal correlations in p + p collisions at root s = 200 GeV to help in estimating nonflow effects. To study the system-size dependence of elliptic flow, we present a detailed comparison with previously published results from Au + Au collisions at root s(NN) = 200 GeV. We observe that v(2)(p(T)) of strange hadrons has similar scaling properties as were first observed in Au + Au collisions, that is, (i) at low transverse momenta, p(T) < 2 GeV/c, v(2) scales with transverse kinetic energy, m(T) - m, and (ii) at intermediate p(T), 2 < p(T) < 4 GeV/c, it scales with the number of constituent quarks, n(q.) We have found that ideal hydrodynamic calculations fail to reproduce the centrality dependence of v(2)(p(T)) for K(S)(0) and Lambda. Eccentricity scaled v(2) values, v(2)/epsilon, are larger in more central collisions, suggesting stronger collective flow develops in more central collisions. The comparison with Au + Au collisions, which go further in density, shows that v(2)/epsilon depends on the system size, that is, the number of participants N(part). This indicates that the ideal hydrodynamic limit is not reached in Cu + Cu collisions, presumably because the assumption of thermalization is not attained.

Eddy Formation in the Tropical Atlantic Induced by Abrupt Changes in the Meridional Overturning Circulation

The variability of the meridional overturning circulation (MOC) in the upper tropical Atlantic basin is investigated using a reduced-gravity model in a simplified domain. Four sets of idealized numerical experiments are performed: (i) switch-on of the MOC until a fixed value when a constant northward flow is applied along the western boundary; (ii) MOC with a variable flow; (iii) MOC in a quasi-steady flow; and (iv) shutdown of the MOC in the Northern Hemisphere. Results from experiments (i) show that eddies are generated at the equatorial region by shear instability and detached northward; eddies are responsible for an enhancement of the mean flow and the variability of the MOC. Results from experiments (ii) show a transitional behavior of the MOC related to the eddy generation in interannual-decadal time scales as the Reynolds number varies due to the variations in the MOC. In experiments (iii), a critical Reynolds number Re(c) around 30 is found, above which eddies are generated. Experiments (iv) demonstrate that even after the collapse of MOC in the Northern Hemisphere, eddies can still be generated and carry energy across the equator into the Northern Hemisphere; these eddies act to attenuate the impact of the MOC shutdown on short time scales. The results described here may be particularly pertinent to ocean general circulation models in which the Reynolds number lies close to the bifurcation point separating the laminar and turbulent regimes.