Solving large-scale dynamic systems using band Lanczos method in Rockwell NASTRAN on CRAY X-MP

The improved cost effectiveness using better models, more accurate and faster algorithms and large scale computing offers more representative dynamic analyses. The band Lanczos eigen-solution method was implemented in Rockwell's version of 1984 COSMIC-released NASTRAN finite element structural analysis computer program to effectively solve for structural vibration modes including those of large complex systems exceeding 10,000 degrees of freedom. The Lanczos vectors were re-orthogonalized locally using the Lanczos Method and globally using the modified Gram-Schmidt method for sweeping rigid-body modes and previously generated modes and Lanczos vectors. The truncated band matrix was solved for vibration frequencies and mode shapes using Givens rotations. Numerical examples are included to demonstrate the cost effectiveness and accuracy of the method as implemented in ROCKWELL NASTRAN. The CRAY version is based on RPK's COSMIC/NASTRAN. The band Lanczos method was more reliable and accurate and converged faster than the single vector Lanczos Method. The band Lanczos method was comparable to the subspace iteration method which was a block version of the inverse power method. However, the subspace matrix tended to be fully populated in the case of subspace iteration and not as sparse as a band matrix

Charge redistribution in the formation of one-dimensional lithium wires on Cu(001)

We describe the formation of one-dimensional lithium wires on a Cu(001) substrate, providing an atomic-scale description of the onset of metallization in this prototypical adsorption system. A combination of helium atom scattering and density-functional theory reveals pronounced changes in the electronic charge distribution on the formation of the c(5√2×√2)R45° Li/Cu(001) structure, as in-plane bonds are created. Charge donation from Li-substrate bonds is found to facilitate the formation of stable, bonded, and depolarized chains of Li adatoms that coexist with an interleaved phase of independent adatoms. The resultant overlayer has a commensurate charge distribution and lattice modulations but differs fundamentally from structurally similar charge-density wave systems

Ab Initio Liquid Hydrogen Muon Cooling Simulations with ELMS in ICOOL

This paper presents new theoretical results on the passage of muons through liquid hydrogen which have been confirmed in a recent experiment. These are used to demonstrate that muon bunches may be compressed by ionisation cooling more effectively than suggested by previous calculations. Muon cooling depends on the differential cross section for energy loss and scattering of muons. We have calculated this cross section for liquid H2 from first principles and atomic data, avoiding traditional assumptions. Thence, 2-D probability maps of energy loss and scattering in mm-scale thicknesses are derived by folding, and stored in a database. Large first-order correlations between energy loss and scattering are found for H2, which are absent in other simulations. This code is named ELMS, Energy Loss & Multiple Scattering. Single particle trajectories may then be tracked by Monte Carlo sampling from this database on a scale of 1 mm or less. This processor has been inserted into the cooling code ICOOL. Significant improvements in 6-D muon cooling are predicted compared with previous predictions based on GEANT. This is examined in various geometries. The large correlation effect is found to have only a small effect on cooling. The experimental scattering observed for liquid H2 in the MUSCAT experiment has recently been reported to be in good agreement with the ELMS prediction, but in poor agreement with GEANT simulation.Comment: 6 pages, 3 figure

Domino tilings and the six-vertex model at its free fermion point

At the free-fermion point, the six-vertex model with domain wall boundary conditions (DWBC) can be related to the Aztec diamond, a domino tiling problem. We study the mapping on the level of complete statistics for general domains and boundary conditions. This is obtained by associating to both models a set of non-intersecting lines in the Lindstroem-Gessel-Viennot (LGV) scheme. One of the consequence for DWBC is that the boundaries of the ordered phases are described by the Airy process in the thermodynamic limit.Comment: 14 pages, 8 figure

Minimal Brownian Ratchet: An Exactly Solvable Model

We develop an exactly-solvable three-state discrete-time minimal Brownian ratchet (MBR), where the transition probabilities between states are asymmetric. By solving the master equations we obtain the steady-state probabilities. Generally the steady-state solution does not display detailed balance, giving rise to an induced directional motion in the MBR. For a reduced two-dimensional parameter space we find the null-curve on which the net current vanishes and detailed balance holds. A system on this curve is said to be balanced. On the null-curve, an additional source of external random noise is introduced to show that a directional motion can be induced under the zero overall driving force. We also indicate the off-balance behavior with biased random noise.Comment: 4 pages, 4 figures, RevTex source, General solution added. To be appeared in Phys. Rev. Let

Algebraic arctic curves in the domain-wall six-vertex model

The arctic curve, i.e. the spatial curve separating ordered (or `frozen') and disordered (or `temperate) regions, of the six-vertex model with domain wall boundary conditions is discussed for the root-of-unity vertex weights. In these cases the curve is described by algebraic equations which can be worked out explicitly from the parametric solution for this curve. Some interesting examples are discussed in detail. The upper bound on the maximal degree of the equation in a generic root-of-unity case is obtained.Comment: 15 pages, no figures; v2: metadata correcte