Optimal iterative schemes for computing the Moore-Penrose matrix inverse

Among the iterative schemes for computing the Moore — Penrose inverse of a woll-conditioned matrix, only those which have an order of convergence three or two are computationally efficient. A Fortran programme for these schemes is provided.

Moore meets maxwell

Moore's Law has driven the semiconductor revolution enabling over four decades of scaling in frequency, size, complexity, and power. However, the limits of physics are preventing further scaling of speed, forcing a paradigm shift towards multicore computing and parallelization. In effect, the system is taking over the role that the single CPU was playing: high-speed signals running through chips but also packages and boards connect ever more complex systems. High-speed signals making their way through the entire system cause new challenges in the design of computing hardware. Inductance, phase shifts and velocity of light effects, material resonances, and wave behavior become not only prevalent but need to be calculated accurately and rapidly to enable short design cycle times. In essence, to continue scaling with Moore's Law requires the incorporation of Maxwell's equations in the design process. Incorporating Maxwell's equations into the design flow is only possible through the combined power that new algorithms, parallelization and high-speed computing provide. At the same time, incorporation of Maxwell-based models into circuit and system-level simulation presents a massive accuracy, passivity, and scalability challenge. In this tutorial, we navigate through the often confusing terminology and concepts behind field solvers, show how advances in field solvers enable integration into EDA flows, present novel methods for model generation and passivity assurance in large systems, and demonstrate the power of cloud computing in enabling the next generation of scalable Maxwell solvers and the next generation of Moore's Law scaling of systems. We intend to show the truly symbiotic growing relationship between Maxwell and Moore!

Evidence for a new metastable crystalline compound in Ge---Te system

Pressure dependence of the electrical resistivity of bulk, melt quenched GexTe100−x glasses (15 less-than-or-equals, slant x less-than-or-equals, slant 28) has been studied up to 8GPa pressure. All the glasses exhibit a sharp, discontinuous glass to crystal transition under pressure. The high pressure crystalline phases are identified to have a face centered cubic structure. The value of the cell constant is 0.779nm for 15 less-than-or-equals, slant x less-than-or-equals, slant 17, 0.642nm for x=20 and 0.55lnm for 22 ≤ x ≤ 28 samples respectively. The cell constants of the high pressure crystalline phases suggest the possible existance of a new metastable crystalline compound in the Ge---Te system with F.C.C. structure and cell constant equal to 1.109nm as reported by Moore et al.

Dynamics after a sweep through a quantum critical point

The coherent quantum evolution of a one-dimensional many-particle system after slowly sweeping the Hamiltonian through a critical point is studied using a generalized quantum Ising model containing both integrable and nonintegrable regimes. It is known from previous work that universal power laws of the sweep rate appear in such quantities as the mean number of excitations created by the sweep. Several other phenomena are found that are not reflected by such averages: there are two different scaling behaviors of the entanglement entropy and a relaxation that is power law in time rather than exponential. The final state of evolution after the quench is not characterized by any effective temperature, and the Loschmidt echo converges algebraically for long times, with cusplike singularities in the integrable case that are dynamically broadened by nonintegrable perturbations.

Ring-Opening-Induced Toughening of a Low-Permittivity Polymer-Metal Interface

Integrating low dielectric permittivity (low-k) polymers to metals is an exacting fundamental challenge because poor bonding between low-polarizability moieties and metals precludes good interfacial adhesion. Conventional adhesion-enhancing methods such as using intermediary layers are unsuitable for engineering polymer/metal interfaces for many applications because of the collateral increase in dielectric permittivity. Here, we demonstrate a completely new approach without surface treatments or intermediary layers to obtain an excellent interfacial fracture toughness of > 13 J/m(2) in a model system comprising copper. and a cross-linked polycarbosilane with k similar to 2.7 obtained by curing a cyclolinear polycarbosilane in air.Our results suggest that interfacial oxygen catalyzed molecularring-opening and anchoring of the opened ring moieties of the polymer to copper is the main toughening mechanism. This novel approach of realizing adherent low-k polymer/metal structures without intermediary layers by activating metal-anchoring polymer moieties at the interface could be adapted for applications such as device wiring and packaging, and laminates and composites.

Ordering near the percolation threshold in models of two-dimensional interacting bosons with quenched dilution

Randomly diluted quantum boson and spin models in two dimensions combine the physics of classical percolation with the well-known dimensionality dependence of ordering in quantum lattice models. This combination is rather subtle for models that order in two dimensions but have no true order in one dimension, as the percolation cluster near threshold is a fractal of dimension between 1 and 2: two experimentally relevant examples are the O(2) quantum rotor and the Heisenberg antiferromagnet. We study two analytic descriptions of the O(2) quantum rotor near the percolation threshold. First a spin-wave expansion is shown to predict long-ranged order, but there are statistically rare points on the cluster that violate the standard assumptions of spin-wave theory. A real-space renormalization group (RSRG) approach is then used to understand how these rare points modify ordering of the O(2) rotor. A new class of fixed points of the RSRG equations for disordered one-dimensional bosons is identified and shown to support the existence of long-range order on the percolation backbone in two dimensions. These results are relevant to experiments on bosons in optical lattices and superconducting arrays, and also (qualitatively) for the diluted Heisenberg antiferromagnet La-2(Zn,Mg)(x)Cu1-xO4.

Tuning the electronic effective mass in double-doped SrTiO3

We elucidate the relationship between effective mass and carrier concentration in an oxide semiconductor controlled by a double-doping mechanism. In this model oxide system, Sr1-xLaxTiO3-delta, we can tune the effective mass ranging from 6 to 20m(e) as a function of filling (carrier concentration) and the scattering mechanism, which are dependent on the chosen lanthanum-and oxygen-vacancy concentrations. The effective mass values were calculated from the Boltzmann transport equation using the measured transport properties of thin films of Sr1-xLaxTiO3-delta. We show that the effective mass decreases with carrier concentration in this large-band-gap, low-mobility oxide, and this behavior is contrary to the traditional high-mobility, small-effective-mass semiconductors.

Concise row-pruning algorithm to invert a matrix

Presented here, in a vector formulation, is an O(mn2) direct concise algorithm that prunes/identifies the linearly dependent (ld) rows of an arbitrary m X n matrix A and computes its reflexive type minimum norm inverse A(mr)-, which will be the true inverse A-1 if A is nonsingular and the Moore-Penrose inverse A+ if A is full row-rank. The algorithm, without any additional computation, produces the projection operator P = (I - A(mr)- A) that provides a means to compute any of the solutions of the consistent linear equation Ax = b since the general solution may be expressed as x = A(mr)+b + Pz, where z is an arbitrary vector. The rank r of A will also be produced in the process. Some of the salient features of this algorithm are that (i) the algorithm is concise, (ii) the minimum norm least squares solution for consistent/inconsistent equations is readily computable when A is full row-rank (else, a minimum norm solution for consistent equations is obtainable), (iii) the algorithm identifies ld rows, if any, and reduces concerned computation and improves accuracy of the result, (iv) error-bounds for the inverse as well as the solution x for Ax = b are readily computable, (v) error-free computation of the inverse, solution vector, rank, and projection operator and its inherent parallel implementation are straightforward, (vi) it is suitable for vector (pipeline) machines, and (vii) the inverse produced by the algorithm can be used to solve under-/overdetermined linear systems.

Climate forcing and response to idealized changes in surface latent and sensible heat

Land use and land cover changes affect the partitioning of latent and sensible heat, which impacts the broader climate system. Increased latent heat flux to the atmosphere has a local cooling influence known as `evaporative cooling', but this energy will be released back to the atmosphere wherever the water condenses. However, the extent to which local evaporative cooling provides a global cooling influence has not been well characterized. Here, we perform a highly idealized set of climate model simulations aimed at understanding the effects that changes in the balance between surface sensible and latent heating have on the global climate system. We find that globally adding a uniform 1 W m(-2) source of latent heat flux along with a uniform 1 W m(-2) sink of sensible heat leads to a decrease in global mean surface air temperature of 0.54 +/- 0.04 K. This occurs largely as a consequence of planetary albedo increases associated with an increase in low elevation cloudiness caused by increased evaporation. Thus, our model results indicate that, on average, when latent heating replaces sensible heating, global, and not merely local, surface temperatures decrease.

Modeling of wheeled mobile robots using dextrous manipulation kinematics

Abstract—This document introduces a new kinematic simulation of a wheeled mobile robot operating on uneven terrain. Our modeling method borrows concepts from dextrous manipulation. This allows for an accurate simulation of the way 3-dimensional wheels roll over a smooth ground surface. The purpose of the simulation is to validate a new concept for design of off-road wheel suspensions, called Passive Variable Camber (PVC). We show that PVC eliminates kinematic slip for an outdoor robot. Both forward and inverse kinematics are discussed and simulation results are presented.

Structure and properties of zinc alloy based metal matrix composites

Zinc-aluminium cast alloys (ZA alloys) exhibit good castability and mechanical properties but these alloys lack creep resistance and high temperature stability. One solution to improve these properties is to reinforce with ceramic particles or fibres, to result in MMCs. MMCs can be produced using casting technique involving infiltration. A systematic investigation was taken and this paper discusses the salient findings of the study on the ZA-27 alloy based MMCs produced through squeeze casting. (Reinforcing fibers: SAFFIL (chopped alumina) or mullite.)

Optimal Linear Glauber Model

Contrary to the actual nonlinear Glauber model, the linear Glauber model (LGM) is exactly solvable, although the detailed balance condition is not generally satisfied. This motivates us to address the issue of writing the transition rate () in a best possible linear form such that the mean squared error in satisfying the detailed balance condition is least. The advantage of this work is that, by studying the LGM analytically, we will be able to anticipate how the kinetic properties of an arbitrary Ising system depend on the temperature and the coupling constants. The analytical expressions for the optimal values of the parameters involved in the linear are obtained using a simple Moore-Penrose pseudoinverse matrix. This approach is quite general, in principle applicable to any system and can reproduce the exact results for one dimensional Ising system. In the continuum limit, we get a linear time-dependent Ginzburg-Landau equation from the Glauber's microscopic model of non-conservative dynamics. We analyze the critical and dynamic properties of the model, and show that most of the important results obtained in different studies can be reproduced by our new mathematical approach. We will also show in this paper that the effect of magnetic field can easily be studied within our approach; in particular, we show that the inverse of relaxation time changes quadratically with (weak) magnetic field and that the fluctuation-dissipation theorem is valid for our model.

Positive isotropic curvature and self-duality in dimension 4

We study a positivity condition for the curvature of oriented Riemannian 4-manifolds: the half-PIC condition. It is a slight weakening of the positive isotropic curvature (PIC) condition introduced by M. Micallef and J. Moore. We observe that the half-PIC condition is preserved by the Ricci flow and satisfies a maximality property among all Ricci flow invariant positivity conditions on the curvature of oriented 4-manifolds. We also study some geometric and topological aspects of half-PIC manifolds.