Sensitivity analysis of the greatest eigenvalue of a symmetric matrix via the ɛ-subdifferential of the associated convex quadratic form

LetA(·) be ann × n symmetric affine matrix-valued function of a parameteruR ^m , and let (u) be the greatest eigenvalue ofA(u). Recently, there has been interest in calculating (u), the subdifferential of atu, which is useful for both the construction of efficient algorithms for the minimization of (u) and the sensitivity analysis of (u), namely, the perturbation theory of (u). In this paper, more generally, we investigate the Legendre-Fenchel conjugate function of (·) and the -subdifferential (u) of atu. Then, we discuss relations between the set (u) and some perturbation bounds for (u).The author is deeply indebted to Professor J. B. Hiriart-Urruty who suggested this study and provided helpful advice and constant encouragement. The author also thanks the referees and the editors for their substantial help in the improvement of this paper.     

Generalized Bethe ansatz solution of a one-dimensional asymmetric exclusion process on a ring with blockage

We present a model for a one-dimensional anisotropic exclusion process describing particles moving deterministically on a ring of lengthL with a single defect, across which they move with probability 0 p 1. This model is equivalent to a two-dimensional, six-vertex model in an extreme anisotropic limit with a defect line interpolating between open and periodic boundary conditions. We solve this model with a Bethe ansatz generalized to this kind of boundary condition. We discuss in detail the steady state and derive exact expressions for the currentj, the density profilen(x), and the two-point density correlation function. In the thermodynamic limitL the phase diagram shows three phases, a low-density phase, a coexistence phase, and a high-density phase related to the low-density phase by a particle-hole symmetry. In the low-density phase the density profile decays exponentially with the distance from the boundary to its bulk value on a length scale . On the phase transition line diverges and the currentj approaches its critical valuej _c = p as a power law,j _c – j ^–1/2. In the coexistence phase the width of the interface between the high-density region and the low-density region is proportional toL ^1/2 if the density f 1/2 and=0 independent ofL if = 1/2. The (connected) two-point correlation function turns out to be of a scaling form with a space-dependent amplitude n(x₁, x₂) =A(x₂)A ^Ke^–r/ withr = x ₂ –x ₁ and a critical exponent = 0.     

A generalized Dantzig—Wolfe decomposition principle for a class of nonconvex programming problems

Since Dantzig—Wolfe's pioneering contribution, the decomposition approach using a pricing mechanism has been developed for a wide class of mathematical programs. For convex programs a linear space of Lagrangean multipliers is enough to define price functions. For general mathematical programs the price functions could be defined by using a subclass of nondecreasing functions. However the space of nondecreasing functions is no longer finite dimensional. In this paper we consider a specific nonconvex optimization problem min {f(x):h _j (x)g(x),j=1, ,m, xX}, wheref(·),h _j (·) andg(·) are finite convex functions andX is a closed convex set. We generalize optimal price functions for this problem in such a way that the parameters of generalized price functions are defined in a finite dimensional space. Combining convex duality and a nonconvex duality we can develop a decomposition method to find a globally optimal solution.This paper is dedicated to Phil Wolfe on the occasion of his 65th birthday.     

AnL 1 smoothing spline algorithm with cross validation

We propose an algorithm for the computation ofL ₁ (LAD) smoothing splines in the spacesW _M (D), with . We assume one is given data of the formy _i =(f(t _i ) + _i , i=1,...,N with {itt_i} _i=1 ^{N D} , the _i are errors withE( _i )=0, andf is assumed to be inW _M . The LAD smoothing spline, for fixed smoothing parameter0, is defined as the solution,s , of the optimization problem (1/N) _i=1 ^N ¦y _i –g(t _i ¦+J _M (g), whereJ _M (g) is the seminorm consisting of the sum of the squaredL ₂ norms of theMth partial derivatives ofg. Such an LAD smoothing spline,s , would be expected to give robust smoothed estimates off in situations where the _i are from a distribution with heavy tails. The solution to such a problem is a thin plate spline of known form. An algorithm for computings is given which is based on considering a sequence of quadratic programming problems whose structure is guided by the optimality conditions for the above convex minimization problem, and which are solved readily, if a good initial point is available. The data driven selection of the smoothing parameter is achieved by minimizing aCV() score of the form .The combined LAD-CV smoothing spline algorithm is a continuation scheme in 0 taken on the above SQPs parametrized in, with the optimal smoothing parameter taken to be that value of at which theCV() score first begins to increase. The feasibility of constructing the LAD-CV smoothing spline is illustrated by an application to a problem in environment data interpretation.     

Diffusion in a one-dimensional gas of hard point particles

We simulate the classical diffusion of a particle of massM in an infinite one-dimensional system of hard point particles of massm in equilibrium. Each computer run corresponds to about 10⁸ collisions of the diffusive particle. We find that (t) 1/t fort large enough, and a crossover from an M m regime where=2 to=3 forM=m. The diffusion constant has a sharp maximum atM=m. We study moments x(t)² and x(t)⁴, and examine the behavior ofq ² (t)=x(t)⁴/3x(t)²². We find thatq(t)1 (consistent with a normal distribution) in theM limit (for all timest) and in the t limit for allM. On sabbatical leave from IVIC-Instituto Venezolano de Investigaciones Cientificas.     

Mechanism of Reppe's nickel-catalyzed ethyne tetramerization to cyclooctatetraene: a DFT study

In this B3 LYP model study, homoleptic nickel(0) ethyne complexes have been predicted as the catalyst resting state for the title reaction. Ethyne ligand coupling of Ni(C(2)H(2))(3) yields monoethyne nickelacyclopentadiene in the rate-determining step. Ethyne coordination is followed by insertion of an ethyne ligand into the Ni--C sigma bond. A highly strained monoethyne trans-nickelacycloheptatriene is formed. This trans intermediate is unable to reductively eliminate benzene without prior isomerization to a cis-structure. Instead, it rapidly collapses to a nickelacyclononatetraene. Ethyne coordination induces reductive elimination to the cyclooctatetraene complex Ni(eta(2)-C(2)H(2))(eta(2)-C(8)H(8)), followed by facile ligand exchange. Other ethyne coupling pathways have been computed to be less favored. The cyclooctatetraene ligand binds significantly weaker to nickel(0) than ethyne, both for mononuclear, and for dinuclear species. For this reason, C--C bond formation steps at Ni(2)(micro-cot) fragments have been predicted to feature prohibitively high overall reaction barriers.     

DFT study of the conductance of molecular wire:The effect of coupling geometry and intermolecular interaction on the transport properties

The density functional theory (DFT) combining with the non-equilibrium Green functions (NEGF) method is applied to the study of the electronic transport properties for a Di-thiol-benzene (DTB) molecule coupled to two Au(111) surfaces. The dependence of the transport properties on the bias, the coupling geometry of the molecule-electrode interface, and the intermolecular interaction are examined in detail. The results show that the existence of the hydrogen atom at the end of the DTB molecule would significantly decrease the transmission coefficients, and then the differential conductance (dI/dV). By changing the position of the DTB molecule located between two electrodes a maximum value of calculated current is observed. It is also found that the intermolecular interaction will strongly influence the transport properties of the system studied.     

Optimization of St John's Wort flavonoid separation by reversed phase liquid chromatography on a silica-based monolithic column

Summary An optimization procedure for the separation of flavonoids from St John's Wort by reversed-phase high-performance liquid chromatography on monolithic stationary phase is described. Three-component mobile phase systems are studied using experimental design methodology. The starting experimental domain is a triangle, each vertex of which is a pure component. From preliminary isocratic experiments, truncations in the domain are performed leading to a quadrilateral shape working domain with no symmetry. To cope with both separation and analysis time, desirability functions are used. Optimal conditions are finally given by binary systems and the four flavonoids are separated in less than seven minutes.     

A new algorithm to evaluate bielectronic integrals with 1s Slater‐type orbitals obtained by the integral transforms

This article presents a variation of the integral transform method to evaluate multicenter bielectronic integrals (12|34), with 1s Slater‐type orbitals. It is proved that it is possible to define, out of the expression of (12|34) given by the integral transform method, a function F(q) that has the property of having a unique Q, such that F(Q) = (12|34). Therefore, F(q) may be used to calculate (12|34). It is shown that the evaluation of F(Q) turns out to be simpler than the three‐dimensional integral involved in the calculation of (12|34), and an algorithm is presented to calculate Q. The results show that relative errors on the order of 10^?3 or lower are obtained very efficiently. In addition, it is shown that the proposed algorithm is very stable. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004