Construction of optimal multi-level supersaturated designs

A supersaturated design is a design whose run size is not large enough for estimating all the main effects. The goodness of multi-level supersaturated designs can be judged by the generalized minimum aberration criterion proposed by Xu and Wu [Ann. Statist. 29 (2001) 1066--1077]. A new lower bound is derived and general construction methods are proposed for multi-level supersaturated designs. Inspired by the Addelman--Kempthorne construction of orthogonal arrays, several classes of optimal multi-level supersaturated designs are given in explicit form: Columns are labeled with linear or quadratic polynomials and rows are points over a finite field. Additive characters are used to study the properties of resulting designs. Some small optimal supersaturated designs of 3, 4 and 5 levels are listed with their properties.Comment: Published at http://dx.doi.org/10.1214/009053605000000688 in the Annals of Statistics (http://www.imstat.org/aos/) by the Institute of Mathematical Statistics (http://www.imstat.org

Runup and rundown generated by three-dimensional sliding masses

To study the waves and runup/rundown generated by a sliding mass, a numerical simulation model, based on the large-eddy-simulation (LES) approach, was developed. The Smagorinsky subgrid scale model was employed to provide turbulence dissipation and the volume of fluid (VOF) method was used to track the free surface and shoreline movements. A numerical algorithm for describing the motion of the sliding mass was also implemented. To validate the numerical model, we conducted a set of large-scale experiments in a wave tank of 104m long, 3.7m wide and 4.6m deep with a plane slope (1:2) located at one end of the tank. A freely sliding wedge with two orientations and a hemisphere were used to represent landslides. Their initial positions ranged from totally aerial to fully submerged, and the slide mass was also varied over a wide range. The slides were instrumented to provide position and velocity time histories. The time-histories of water surface and the runup at a number of locations were measured. Comparisons between the numerical results and experimental data are presented only for wedge shape slides. Very good agreement is shown for the time histories of runup and generated waves. The detailed three-dimensional complex flow patterns, free surface and shoreline deformations are further illustrated by the numerical results. The maximum runup heights are presented as a function of the initial elevation and the specific weight of the slide. The effects of the wave tank width on the maximum runup are also discussed

Possible Weyl fermions in the magnetic Kondo system CeSb

Materials where the electronic bands have unusual topologies allow for the realization of novel physics and have a wide range of potential applications. When two electronic bands with linear dispersions intersect at a point, the excitations could be described as Weyl fermions which are massless particles with a particular chirality. Here we report evidence for the presence of Weyl fermions in the ferromagnetic state of the low-carrier density, strongly correlated Kondo lattice system CeSb, from electronic structure calculations and angle-dependent magnetoresistance measurements. When the applied magnetic field is parallel to the electric current, a pronounced negative magnetoresistance is observed within the ferromagnetic state, which is destroyed upon slightly rotating the field away. These results give evidence for CeSb belonging to a new class of Kondo lattice materials with Weyl fermions in the ferromagnetic state.Comment: 18 pages, 4 figures, Supplementary Information available from journal link (open access

Critical manifold of the kagome-lattice Potts model

Any two-dimensional infinite regular lattice G can be produced by tiling the plane with a finite subgraph B of G; we call B a basis of G. We introduce a two-parameter graph polynomial P_B(q,v) that depends on B and its embedding in G. The algebraic curve P_B(q,v) = 0 is shown to provide an approximation to the critical manifold of the q-state Potts model, with coupling v = exp(K)-1, defined on G. This curve predicts the phase diagram both in the ferromagnetic (v>0) and antiferromagnetic (v<0) regions. For larger bases B the approximations become increasingly accurate, and we conjecture that P_B(q,v) = 0 provides the exact critical manifold in the limit of infinite B. Furthermore, for some lattices G, or for the Ising model (q=2) on any G, P_B(q,v) factorises for any choice of B: the zero set of the recurrent factor then provides the exact critical manifold. In this sense, the computation of P_B(q,v) can be used to detect exact solvability of the Potts model on G. We illustrate the method for the square lattice, where the Potts model has been exactly solved, and the kagome lattice, where it has not. For the square lattice we correctly reproduce the known phase diagram, including the antiferromagnetic transition and the singularities in the Berker-Kadanoff phase. For the kagome lattice, taking the smallest basis with six edges we recover a well-known (but now refuted) conjecture of F.Y. Wu. Larger bases provide successive improvements on this formula, giving a natural extension of Wu's approach. The polynomial predictions are in excellent agreement with numerical computations. For v>0 the accuracy of the predicted critical coupling v_c is of the order 10^{-4} or 10^{-5} for the 6-edge basis, and improves to 10^{-6} or 10^{-7} for the largest basis studied (with 36 edges).Comment: 31 pages, 12 figure

Applications of numerical codes to space plasma problems

Solar wind, earth's bowshock, and magnetospheric convection and substorms were investigated. Topics discussed include computational physics, multifluid codes, ionospheric irregularities, and modeling laser plasmas