Truncated Newton method for sparse unconstrained optimization using automatic differentiation

When solving large complex optimization problems, the user is faced with three major problems. These are (i) the cost in human time in obtaining accurate expressions for the derivatives involved; (ii) the need to store second derivative information; and (iii), of lessening importance, the time taken to solve the problem on the computer. For many problems, a significant part of the latter can be attributed to solving Newton-like equations. In the algorithm described, the equations are solved using a conjugate direction method that only needs the Hessian at the current point when it is multiplied by a trial vector. In this paper, we present a method that finds this product using automatic differentiation while only requiring vector storage. The method takes advantage of any sparsity in the Hessian matrix and computes exact derivatives. It avoids the complexity of symbolic differentiation, the inaccuracy of numerical differentiation, the labor of finding analytic derivatives, and the need for matrix store. When far from a minimum, an accurate solution to the Newton equations is not justified, so an approximate solution is obtained by using a version of Dembo and Steihaug's truncated Newton algorithm (Ref. 1).This paper was presented at the SIAM National Meeting, Boston, Massachusetts, 1986.     

3D-TOFSIMS characterization of black spots in polymer light emitting diodes

The occurrence and formation of black spots areas in PolyLED devices has been studied by time-of-flight SIMS (TOFSIMS). The composition, shape and position of the black spots is visualised by three-dimensional (3D)-TOFSIMS depth-profiling. It has been established that the formation of non-emissive spots is due to the growth of aluminium oxide clusters at the AlBa/polymer interface. Electron injection in the black spots is lost by the resulting local increase of the resistivity of the cathode.     

A model for adatom structures

A model concerning adatom structures is proposed. Attractive nearest neighbour interactions, which may be of electronic nature lead to 2-dimensional condensation. Every pair bond causes and elastic dipole. The elastic dipoles interact via substrate strains with an anisotropics ^–3 power law. Different types of adatoms or sites are permitted and many-body effects result, from the assumptions. Electric dipole interactions of adatoms are included for comparison. The model is applied to the W(110) surface and compared with superstructures experimentally found in the W(110)-0 system. It is found that there is still lack for an additional next-nearest neighbour interaction.     

Timing and detuning studies of a synchronously pumped mode-locked dye laser

The timing relationship between a synchronously pumped mode-locked dye laser and its pump laser is explored experimentally and theoretically. The experimental data is in quantitative agreement with numerical solutions of a semiclassical model that includes no free parameters. The data also agrees qualitatively with an approximate, analytic distillation of this model which is presented here for the first time. Our study has a direct bearing on the jitter and the long term stability of these important ultrashort pulse lasers.     

Theory of interfacial phase transitions in surfactant systems

The spin-1 Ising model, which is equivalent to the three-component lattice gas model, is used to study wetting transitions in three-component surfactant systems consisting of an oil, water, and a nonionic surfactant. Phase equilibria, interfacial profiles, and interfacial tensions for three-phase equilibrium are determined in mean field approximation, for a wide range of temperature and interaction parameters. Surfactant interaction parameters are found to strongly influence interfacial tensions, reducing them in some cases to ultralow values. Interfacial tensions are used to determine whether the middle phase, rich in surfactant, wets or does not wet the interface between the oil-rich and water-rich phases. By varying temperature and interaction parameters, a wetting transition is located and found to be of the first order. Comparison is made with recent experimental results on wetting transitions in ternary surfactant systems.This paper is dedicated to J. K. Percus in honor of his 65th birthday.     

An applied cell mapping method for optimal control problems

From the application point of view, a series of modifications are proposed for the cell mapping method discussed in Ref. 1 for the optimal control analysis of dynamical systems. The cell order around the target set is rearranged. A set of common discriminate principles is used for the selection of the optimal one among competing control strategies of the same cost. Inequality constraints of the system are taken into account. The number of elements in the set of allowable time intervals is not prescribed, but left open. These modifications seem to make the cell mapping method more efficient for analyzing feedback systems and for obtaining their global optimal control information. The algorithms presented in this paper could broaden the application of the cell mapping approach of Ref. 1 to a wider class of engineering problems.     

Rehabilitation of the Gauss-Jordan algorithm

Summary In this paper a Gauss-Jordan algorithm with column interchanges is presented and analysed. We show that, in contrast with Gaussian elimination, the Gauss-Jordan algorithm has essentially differing properties when using column interchanges instead of row interchanges for improving the numerical stability. For solutions obtained by Gauss-Jordan with column interchanges, a more satisfactory bound for the residual norm can be given. The analysis gives theoretical evidence that the algorithm yields numerical solutions as good as those obtained by Gaussian elimination and that, in most practical situations, the residuals are equally small. This is confirmed by numerical experiments. Moreover, timing experiments on a Cyber 205 vector computer show that the algorithm presented has good vectorisation properties.     

An expedient route to diaza-spirocycles utilizing a sequential multicomponent α-aminoallylation/ring-closing metathesis strategy

A general method for the preparation of diaza-spirocycles is reported. This method used an olefin metathesis in order to construct the desired spirocyclic framework. Beginning with commercially available protected amino ketones, this strategy ultimately produced pharmacologically relevant diaza-scaffolds in an efficient and high-yielding process.