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.     

不完全信息群体决策专家权重的集结

本文对于属性权重信息和属性效用信息都不完全的群体多属性决策问题，通过构造属性值区间和运用系统聚类分析法，对群体决策中的专家进行分类，并确定每位专家的权重．     

Non‐empirical calculations of NMR indirect carbon–carbon coupling constants. Part 4: Bicycloalkanes

A systematic study of the one‐bond and long‐range J(C,C), J(C,H) and J(H,H) in the series of nine bicycloalkanes was performed at the SOPPA level with special emphasis on the coupling transmission mechanisms at bridgeheads. Many unknown couplings were predicted with high reliability. Further refinement of SOPPA computational scheme adjusted for better performance was carried out using bicyclo[1.1.1]pentane as a benchmark to investigate the influence of geometry, basis set and electronic correlation. The calculations performed demonstrated that classical ab initio SOPPA applied with the locally dense Dunning's sets augmented with inner core s‐functions used for coupled carbons and Sauer's sets augmented with tight s‐functions used for coupled hydrogens performs perfectly well in reproducing experimental values of different types of coupling constants (the estimated reliability is ca 1–2 Hz) in relatively large organic molecules of up to 11 carbon atoms. Additive coupling increments were derived for J(C,C), J(C,H) and J(H,H) based on the calculated values of coupling constants within SOPPA in the model bicycloalkanes, in reasonably good agreement with the known values obtained earlier on pure empirical grounds. Most of the bridgehead couplings in all but one bicycloalkane appeared to be essentially additive within ca 2–3 Hz while bicyclo[1.1.1]pentane demonstrated dramatic non‐additivity of ?14.5 Hz for J(C,C), +16.6 Hz for J(H,H) and ?5.5 Hz for J(C,H), in line with previous findings. Non‐additivity effects in the latter compound established at the SOPPA level should be attributed to the through‐space non‐bonded interactions at bridgeheads due to the essential overlapping of the bridgehead rear lobes which provides an additional and effective non‐bonding coupling path for the bridgehead carbons and their protons in the bicyclopentane framework. Copyright © 2003 John Wiley & Sons, Ltd.     

Simulation of viscous water column collapse using adapting hierarchical grids

An adaptive hierarchical grid‐based method for predicting complex free surface flows is used to simulate collapse of a water column. Adapting quadtree grids are combined with a high‐resolution interface‐capturing approach and pressure‐based coupling of the Navier–Stokes equations. The Navier–Stokes flow solution scheme is verified for simulation of flow in a lid‐driven cavity at Re=1000. Two approaches to the coupling of the Navier–Stokes equations are investigated as are alternative face velocity and hanging node interpolations. Collapse of a water column as well as collapse of a water column and its subsequent interaction with an obstacle are simulated. The calculations are made on uniform and adapting quadtree grids, and the accuracy of the quadtree calculations is shown to be the same as those made on the equivalent uniform grids. Results are in excellent agreement with experimental and other numerical data. A sharp interface is maintained at the free surface. The new adapting quadtree‐based method achieves a considerable saving in the size of the computational grid and CPU time in comparison with calculations made on equivalent uniform grids. Copyright © 2005 John Wiley & Sons, Ltd.     

The Gram Matrix of a PT-Symmetric Quantum System

The eigenstates of a diagonalizable PT-symmetric Hamiltonian satisfy unconventional completeness and orthonormality relations. These relations reflect the properties of a pair of bi-orthonormal bases associated with non-hermitean diagonalizable operators. In a similar vein, such a dual pair of bases is shown to possess, in the presence of PT symmetry, a Gram matrix of a particular structure: its inverse is obtained by simply swapping the signs of some its matrix elements.     

Stereological simulations of asymmetric polymerization in channels of bile acid inclusion compounds: A detailed analysis of the monomer arrangements in the channels

We studied simulations by computer graphics to estimate the steric mechanism of the asymmetric polymerization of prochiral diene monomers in channels of inclusion compounds of steroidal bile acids, such as deoxycholic acid (DCA) and cholic acid. We applied a hierarchization method to interpret the crystal structures of bile acids, clarifying that the chiral host molecules associated to form characteristic 2₁-helical assemblies with uneven surfaces. A detailed analysis of the uneven channels in a close-packing state indicated that there were many possible arrangements of the monomers in the channels. The plausible arrangements in the channel could explain a previous study, which showed that the polymerization in the DCA channel yielded chiral polymers with a predominant configuration from prochiral diene monomers, such as 2-methyl-trans-1,3-pentadiene. On the basis of such simulation studies of the arrangements of guest monomers in the channel, we examined a plausible steric mechanism for asymmetric inclusion polymerization. © 2004 Wiley Periodicals, Inc. J Polym Sci Part A: Polym Chem 42: 4648–4655, 2004     

多小波子空间上的单小波表示

本文在较弱的条件下,建立了2重多小波子空间与单小波子空间的关系.即由2重多小波构造出单小波.一方面,这种单小波的平移伸缩与2重多小波的平移伸缩生成的子空间是完全相同的;另一方面,它具有插值性.因此通过构造出的单小波建立了多小波子空间上的Shannon型采样定理.     

Banach空间上的框架与Riesz基

本文讨论Banach空间上框架、无冗框架与Riesz基之间的关系及它们的稳定性．     

反恐怖安全检测技术的物理基础

反恐怖安全检测技术已成为人们的热门话题。安全检测技术的诞生和发展离不开物理学家的贡献。物理学为反恐怖安全检测孕育了多种手段和方法。文章简要介绍了X射线探测、γ射线探测、四极共振分析、脉冲式快中子分析、离子迁移率光谱测量、毫米波探测、声学探测、磁场探测、电磁脉冲探测等安全检测技术的物理基础。     

The Lascoux, Leclerc and Thibon algorithm and Soergel's tilting algorithm

We generalize Soergel's tilting algorithm to singular weights and deduce from this the validity of the Lascoux-Leclerc-Thibon conjecture on the connection between the canonical basis of the basic submodule of the Fock module and the representation theory of the Hecke-algebras at root of unity. Supported in part by Programa Reticulados y Ecuaciones and by FONDECYT grant 1051024.