Optimized strong stability preserving IMEX Runge-Kutta methods

Esta es la versión no revisada del artículo: Inmaculada Higueras, Natalie Happenhofer, Othmar Koch, and Friedrich Kupka. 2014. Optimized strong stability preserving IMEX Runge-Kutta methods. J. Comput. Appl. Math. 272 (December 2014), 116-140. Se puede consultar la versión final en https://doi.org/10.1016/j.cam.2014.05.011

Similar prevalence and magnitude of auditory-evoked and visually evoked activity in the frontal eye fields: implications for multisensory motor control.

Saccadic eye movements can be elicited by more than one type of sensory stimulus. This implies substantial transformations of signals originating in different sense organs as they reach a common motor output pathway. In this study, we compared the prevalence and magnitude of auditory- and visually evoked activity in a structure implicated in oculomotor processing, the primate frontal eye fields (FEF). We recorded from 324 single neurons while 2 monkeys performed delayed saccades to visual or auditory targets. We found that 64% of FEF neurons were active on presentation of auditory targets and 87% were active during auditory-guided saccades, compared with 75 and 84% for visual targets and saccades. As saccade onset approached, the average level of population activity in the FEF became indistinguishable on visual and auditory trials. FEF activity was better correlated with the movement vector than with the target location for both modalities. In summary, the large proportion of auditory-responsive neurons in the FEF, the similarity between visual and auditory activity levels at the time of the saccade, and the strong correlation between the activity and the saccade vector suggest that auditory signals undergo tailoring to match roughly the strength of visual signals present in the FEF, facilitating accessing of a common motor output pathway.

Cold crucible melting of reactive metals using combined DC and AC magnetic fields

Cold crucible furnace is widely used for melting reactive metals for high quality castings. Although the water cooled copper crucible avoids contamination, it produces a low superheat of the melt. Experimental and theoretical investigations of the process showed that the increase of the supplied power to the furnace leads to a saturation in the temperature rise of the melt, and no significant increase of the melt superheat can be obtained. The computer model of theprocess has been developed to simulate the time dependent turbulent flow, heat transfer with phase change, and AC and DC magnetohydrodynamics in a time varying liquid metal envelope. The model predicts that the supermimposition of a strong DC field on top of the normal AC field reduces the level of turbulience and stirring in the liquid metal, thereby reducing the heat loss through the base of the crucible and increasing the superheat. The direct measurements of the temperature in the commercial size cold crucbile has confirmed the computer redictions and showed that the addition of a DC field increased the superheat in molten TiAl from ~45C (AC field only) to ~81C (DC+AC fields). The present paper reports further predictions of the effect of a dDC field on top of the AC field and compares these with experimental data.

The study of flow and temperature fields in conducting droplets suspended in a DC/AC combination field

Electromagnetic Levitation (EML) is a valuable method for measuring the thermo-physical properties of metals - surface tensions, viscosity, thermal/electrical conductivity, specific heat, hemispherical emissivity, etc. – beyond their melting temperature. In EML, a small amount of the test specimen is melted by Joule heating in a suspended AC coil. Once in liquid state, a small perturbation causes the liquid envelope to oscillate and the frequency of oscillation is then used to compute its surface tension by the well know Rayleigh formula. Similarly, the rate at which the oscillation is dampened relates to the viscosity. To measure thermal conductivity, a sinusoidally varying laser source may be used to heat the polar axis of the droplet and the temperature response measured at the polar opposite – the resulting phase shift yields thermal conductivity. All these theoretical methods assume that convective effects due to flow within the droplet are negligible compared to conduction, and similarly that the flow conditions are laminar; a situation that can only be realised under microgravity conditions. Hence the EML experiment is the method favoured for Spacelab experiments (viz. TEMPUS). Under terrestrial conditions, the full gravity force has to be countered by a much larger induced magnetic field. The magnetic field generates strong flow within the droplet, which for droplets of practical size becomes irrotational and turbulent. At the same time the droplet oscillation envelope is no longer ellipsoidal. Both these conditions invalidate simple theoretical models and prevent widespread EML use in terrestrial laboratories. The authors have shown in earlier publications that it is possible to suppress most of the turbulent convection generated in the droplet skin layer, through use of a static magnetic field. Using a pseudo-spectral discretisation method it is possible compute very accurately the dynamic variation in the suspended fluid envelope and simultaneously compute the time-varying electromagnetic, flow and thermal fields. The use of a DC field as a dampening agent was also demonstrated in cold crucible melting, where suppression of turbulence was achieved in a much larger liquid metal volume and led to increased superheat in the melt and reduction of heat losses to the water-cooled walls. In this paper, the authors describe the pseudo-spectral technique as applied to EML to compute the combined effects of AC and DC fields, accounting for all the flow-induced forces acting on the liquid volume (Lorentz, Maragoni, surface tension, gravity) and show example simulations.

Effects of magnetic fields on crystal growth

The effects of a constant uniform magnetic field on thermoelectric currents during dendritic solidification were investigated using a two-dimensional enthalpy based numerical model. Using an approximation for three-dimensional unconstricted growth, the resulting Lorentz forces generate a circulating flow influencing the solidification pattern. Under the presence of a strong magnetic field secondary growth on the clockwise side of the primary arm of the dendrite was encouraged, whereas the anticlockwise side is suppressed due to a reduction in local free energy. The preferred direction of growth rotated in the clockwise sense under an anticlockwise flow. The tip velocity is significantly increased compared with growth in stagnant flow. This is due to a small recirculation at the tip of the dendrite; bringing in colder liquid and lowering the concentration of solute.

Strong parallel magnetic field effects on the hydrogen molecular ion

Equilibrium distances, binding energies and dissociation energies for the ground and low-lying states of the hydrogen molecular ion in a strong magnetic field parallel to the internuclear axis are calculated and refined, by using the two- dimensional pseudospectral method. High-precision results are presented for the binding energies over a wider field regime than already given in the literature (Kravchenko and Liberman 1997 Phys. Rev. A 55 2701). The present work removes a long- standing discrepancy for the R-eq value in the 1sigma(u) state at a field strength of 1.0 x 10(6) T. The dissociation energies of the antibonding 1pi(g) state induced by magnetic fields are determined accurately. We have also observed that the antibonding 1pi(g) potential energy curve develops a minimum if the field is sufficiently strong. Some unreliable results in the literature are pointed out and discussed. A way to efficiently treat vibrational processes and coupling between the nuclear and the electronic motions in magnetic fields is also suggested within a three-dimensional pseudospectral scheme.

Guiding of Relativistic Electron Beams in Solid Targets by Resistively Controlled Magnetic Fields

Guided transport of a relativistic electron beam in solid is achieved experimentally by exploiting the strong magnetic fields created at the interface of two metals of different electrical resistivities. This is of substantial relevance to the Fast Ignitor approach to fusion energy production [M. Tabak et al., Phys. Plasmas 12, 057305 (2005)], since it allows the electron deposition to be spatially tailored-thus adding substantial design flexibility and preventing inefficiencies due to electron beam spreading. In the experiment, optical transition radiation and thermal emission from the target rear surface provide a clear signature of the electron confinement within a high resistivity tin layer sandwiched transversely between two low resistivity aluminum slabs. The experimental data are found to agree well with numerical simulations.

SIMLA: Simulating particle dynamics in intense laser and other electromagnetic fields via classical and quantum electrodynamics

We present the Fortran program SIMLA, which is designed for the study of charged particle dynamics in laser and other background fields. The dynamics can be determined classically via the Lorentz force and Landau–Lifshitz equations or, alternatively, via the simulation of photon emission events determined by strong-field quantum-electrodynamics amplitudes and implemented using Monte-Carlo routines. Multiple background fields can be included in the simulation and, where applicable, the propagation direction, field type (plane wave, focussed paraxial, constant crossed, or constant magnetic), and time envelope of each can be independently specified.

Local conductance: A means to extract polarization and depolarizing fields near domain walls in ferroelectrics

Conducting atomic force microscopy images of bulk semiconducting BaTiO3 surfaces show clear stripe domain contrast. High local conductance correlates with strong out-of-plane polarization (mapped independently using piezoresponse force microscopy), and current- voltage characteristics are consistent with dipole-induced alterations in Schottky barriers at the metallic tip-ferroelectric interface. Indeed, analyzing current-voltage data in terms of established Schottky barrier models allows relative variations in the surface polarization, and hence the local domain structure, to be determined. Fitting also reveals the signature of surface-related depolarizing fields concentrated near domain walls. Domain information obtained from mapping local conductance appears to be more surface-sensitive than that from piezoresponse force microscopy. In the right materials systems, local current mapping could therefore represent a useful complementary technique for evaluating polarization and local electric fields with nanoscale resolution.

Explaining the Coexistence of Large-scale and Small-scale Magnetic Fields in Fully Convective Stars

Despite the lack of a shear-rich tachocline region, low-mass fully convective (FC) stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three-dimensional model of magnetic field generation in low-mass FC stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman–Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass FC stars.

Transverse Spreading of Electrons in High-Intensity Laser Fields

We show that for collisions of electrons with a high-intensity laser, discrete photon emissions introduce a transverse beam spread that is distinct from that due to classical (or beam shape) effects. Via numerical simulations, we show that this quantum induced transverse momentum gain of the electron is manifest in collisions with a realistic laser pulse of intensity within reach of current technology, and we propose it as a measurable signature of strong-field quantum electrodynamics.

Multiple edges of a quantum Hall system in a strong electric

In this paper we show that if the electrons in a quantum Hall sample are subjected to a constant electric field in the plane of the material, comparable in magnitude to the background magnetic field on the system of electrons, a multiplicity of edge states localized at different regions of space is produced in the sample. The actions governing the dynamics of these edge states are obtained starting from the well-known Schrödinger field theory for a system of nonrelativistic electrons, where on top of the constant background electric and magnetic fields, the electrons are further subject to slowly varying weak electromagnetic fields. In the regions between the edges, dubbed as the "bulk," the fermions can be integrated out entirely and the dynamics expressed in terms of a local effective action involving the slowly varying electromagnetic potentials. It is further shown how the bulk action is gauge noninvariant in a particular way, and how the edge states conspire to restore the U(1) electromagnetic gauge invariance of the system. In the edge action we obtain a heretofore unnoticed gauge-invariant term that depends on the particular edge. We argue that this term may be detected experimentally as different edges respond differently to a monochromatic probe due to this term

Exchange-correlation effects on quantum wires with spin-orbit interactions under the influence of in-plane magnetic fields.

Within the noncollinear local spin-density approximation, we have studied the ground state structure of a parabolically confined quantum wire submitted to an in-plane magnetic field, including both Rashba and Dresselhaus spin-orbit interactions. We have explored a wide range of linear electronic densities in the weak (strong) coupling regimes that appear when the ratio of spin-orbit to confining energy is small (large). These results are used to obtain the conductance of the wire. In the strong coupling limit, the interplay between the applied magnetic field¿irrespective of the in-plane direction, the exchange-correlation energy, and the spin-orbit energy-produces anomalous plateaus in the conductance vs linear density plots that are otherwise absent, or washes out plateaus that appear when the exchange-correlation energy is not taken into account.

Organic inputs and farmers' management strategies in millet fields of western Niger

Research on soil fertility management in sub-Saharan Africa was criticized lately for largely ignoring farmers’ management strategies and the underlying principles. To fill this gap of knowledge, detailed interviews were conducted with 108 farm households about their rationale in managing the soil fertility of 307 individual fields in the agro-pastoral village territory of Chikal in western Niger. To amplify the farmers’ information on manuring and corralling practices, repeated measurements of applied amounts of manure were carried out within six 1-km^2 monitoring areas from February to October 1998. The interviews revealed that only 2% of the fields were completely fallowed for a period of 1–15 years, but 40% of the fields were at least partially fallowed. Mulching of crop residues was mainly practiced to fight wind erosion but was restricted to 36% of the surveyed fields given the alternative use of straw as livestock feed. Manure application and livestock corralling, as most effective tools to enhance soil fertility, were targeted to less than 30% of the surveyed fields. The application of complete fallow and manuring and corralling practices were strongly related to the households’ endowment with resources, especially with land and livestock. Within particular fields, measures were mainly applied to spots of poor soil fertility, while the restoration of the productivity of hard pans was of secondary importance. Given the limited spatial coverage of indigenous soil fertility measures and their strong dependence on farmers’ wealth, supplementary strategies to restrict the decline of soil fertility in the drought prone areas of Niger with their heavily weathered soils are needed.

Intrinsic channel closing in strong-field single ionization of H2

Die Ionisation von H2 in intensiven Laserpulsen wird mit Hilfe der numerischen Integration der zeitabhängigen Schrödingergleichung für ein Einelektronenmodell untersucht, das die Vibrationsbewegung berücksichtigt. Die Spektren der kinetischen Elektronenenergie hängen stark von der Vibrationsquantenzahl des erzeugten H2+ Ions ab. Für bestimmte Vibrationszustände ist die Ausbeute der Elektronen in der Mitte des Plateaus stark erhöht. Der Effekt wird "channel closings" zugeschrieben, die in Atomen durch Variation der Laserintensität beobachtet wurden. The ionization of H2 in intense laser pulses is studied by numerical integration of the time-dependent Schrödinger equation for a single-active-electron model including the vibrational motion. The electron kinetic energy spectra in high-order above-threshold ionization are strongly dependent on the vibrational quantum number of the created H2+ ion. For certain vibrational states, the electron yield in the mid-plateau region is strongly enhanced. The effect is attributed to channel closings, which were previously observed in atoms by varying the laser intensity.