2-D R-matrix propagation: A large scale electron scattering simulation dominated by the multiplication of dynamically changing matrices

We examine the computational aspects of propagating a global R-matrix, R, across sub-regions in a 2-D plane. This problem originates in the large scale simulation of electron collisions with atoms and ions at intermediate energies. The propagation is dominated by matrix multiplications which are complicated because of the dynamic nature of R, which changes the designations of its rows and columns and grows in size as the propagation proceeds. The use of PBLAS to solve this problem on distributed memory HPC machines is the main focus of the paper.

Full-Scale Marine Exposure Tests on Treated and Untreated Concretes Initial 7-Year Results

An extensive chloride profiling program was undertaken on concrete pier stems erected in the vicinity of the Dornoch Bridge located at the Dornoch Firth in Northeast Scotland. The pier stems were 2 m (6.562 ft) high and octagonal in plan with 0.66 m (2.165 ft) wide faces. The piers were constructed in sets of three with the lowest of each set in the tidal zone and the highest in the atmospheric zone. The pier stems were placed in such a way that they would represent the exposure conditions of the actual bridge piers of the Dornoch Bridge. In all, six of the pier stems were made using plain ordinary portland cement (OPC) concrete (with three of these having the surface treated with silane); the remaining three pier stems had a concrete containing caltite as an additive. Three exposurezones were studied: the tidal zone, the splash zone, and the atmospheric zone. The tidal zone was further subdivided into two levels defined as low-level and high-level. Chloride profiles were obtained from the different regimes over a period of 7 years for all nine pier stems. This paper describes the nature of chloride ingress and the usefulness of diffusion parameters in classifying each exposure regimes. Furthermore, the effectiveness of silane and caltite in protecting concrete from chloride ingress in different exposure zones was studied.

Late Holocene linkages between decade-century scale climate variability and productivity at Lake Tanganyika, Africa

Microlaminated sediment cores from the Kalya slope region of Lake Tanganyika provide a near-annually resolved paleoclimate record between similar to 2,840 and 1,420 cal. yr B.P. demonstrating strong linkages between climate variability and lacustrine productivity. Laminae couplets comprise dark, terrigenous-dominated half couplets, interpreted as low density underflows deposited from riverine sources during the rainy season, alternating with light, planktonic diatomaceous ooze, with little terrigenous component, interpreted as windy/dry season deposits. Laminated portions of the studied cores consist of conspicuous dark and light colored bundles of laminae couplets. Light and dark bundles alternate at decadal time scales. Within dark bundles, both light and dark half couplets are significantly thinner than within light bundles, implying slower sediment accumulation rates during both seasons over those intervals.

Large-scale density functional calculations of the surface properties of the Wigner crystal

The surface properties of the jellium model have been investigated by large supercell computations in the density functional theory-local spin-density (DFT-LSD) approach for planar slabs with up to 1000 electrons. A wide interval of densities has been explored, extending into the stability range of the Wigner crystal. Most computations have been carried out on nominally paramagnetic samples with an equal number of spin-up and spin-down electrons. The results show that within DFT-LSD spontaneous spin polarization and charge localization start nearly simultaneously at the surface for r(s) similar to 20, then, with decreasing density, they progress toward the center of the slab. Electrons are fully localized and spin polarized at r(s) = 30. At this density the charge distribution is the superposition of disjoint charge blobs, each corresponding to one electron. The distribution of blobs displays both regularities and disorder, the first being represented by well-defined planes and simple in-plane geometries, and the latter by a variety of surface defects. The surface energy, surface dipole, electric polarisability, and magnetization pattern have been determined as a function of density. All these quantities display characteristic anomalies at the density of the localization transition. The analysis of the low-frequency electric conductivity shows that in the fluid paramagnetic regime the in-plane current preferentially flows in the central region of the slab and the two spin channels are equally conducting. In the charge localized, spin-polarized regime, conductivity is primarily a surface effect, and an apparent asymmetry is observed in the two spin currents.

Nonlinear Dynamics of Rotating Multi-Component Pair Plasmas and e-p-i Plasmas

The propagation of small amplitude stationary profile nonlinear electrostatic excitations in a pair plasma is investigated, mainly drawing inspiration from experiments on fullerene pair-ion plasmas. Two distinct pair ion species are considered of opposite polarity and same mass, in addition to a massive charged background species, which is assumed to be stationary, given the frequency scale of interest. In the pair-ion context, the third species is thought of as a background defect (e.g. charged dust) component. On the other hand, the model also applies formally to electron-positron-ion (e-p-i) plasmas, if one neglects electron-positron annihilation. A two-fluid plasma model is employed, incorporating both Lorentz and Coriolis forces, thus taking into account the interplay between the gyroscopic (Larmor) frequency ?c and the (intrinsic) plasma rotation frequency O0. By employing a multi-dimensional reductive perturbation technique, a Zakharov-Kuznetsov (ZK) type equation is derived for the evolution of the electric potential perturbation. Assuming an arbitrary direction of propagation, with respect to the magnetic field, we derive the exact form of nonlinear solutions, and study their characteristics. A parametric analysis is carried out, as regards the effect of the dusty plasma composition (background number density), species temperature(s) and the relative strength of rotation to Larmor frequencies. It is shown that the Larmor and mechanical rotation affect the pulse dynamics via a parallel-to-transverse mode coupling diffusion term, which in fact diverges at ?c ? ±2O0. Pulses collapse at this limit, as nonlinearity fails to balance dispersion. The analysis is complemented by investigating critical plasma compositions, in fact near-symmetric (T- ˜ T+) “pure” (n- ˜ n+) pair plasmas, i.e. when the concentration of the 3rd background species is negligible, case in which the (quadratic) nonlinearity vanishes, so one needs to resort to higher order nonlinear theory. A modified ZK equation is derived and analyzed. Our results are of relevance in pair-ion (fullerene) experiments and also potentially in astrophysical environments, e.g. in pulsars.

Spatially Resolved Measurements of Laser Filamentation in Long Scale Length Underdense Plasmas with and without Beam Smoothing

The onset of filamentation, following the interaction of a relatively long (tau(L) similar or equal to 1 ns) and intense (I-L similar or equal to 5 x 10(14) W/cm(2)) laser pulse with a neopentane filled gas bag target, has been experimentally studied via the proton radiography technique, in conditions of direct relevance to the indirect drive inertial confinement fusion scheme. The density gradients associated with filamentation onset have been spatially resolved yielding direct and unambiguous evidence of filament formation and quantitative information about the filamentation mechanism in agreement with previous theoretical modelings. Experimental data confirm that, once spatially smoothed laser beams are used, filamentation is not a relevant phenomenon during the heating laser beams propagation through typical hohlraum gas fills.