Fitting complex potential energy surfaces to simple model potentials: application of the simplex-annealing method

A stochastic method of optimization, which combines simulated annealing with simplex, is implemented to fit the parameters of a simple model potential. The main characteristic of the method is that it explores the whole space of the parameters of the model potential, and therefore it is very efficient in locating the global minimum of the cost function, in addition to being independent of the initial guess of the parameters. The method is employed to fit the complex intermolecular potential energy surface of the dimer of water, using as a reference the spectroscopic quality anisotropic site-site potential of Feller et al. The simple model potential chosen for its reparameterization is the MCY model potential of Clementi et al. The quality of the fit is assessed by comparing the geometry of the minimum, the harmonic frequencies, and the second virial coefficients of the parameterized potential with the reference one. Finally, to prove more rigorously the robustness of this method, it is compared with standard nonstochastic methods of optimization.     

Counterpoise-corrected potential energy surfaces of simple H-bonded systems

Geometries and stabilization energies of various simple H-bonded complexes (water dimer, hydrogen fluoride dimer, formamide dimer, formic acid dimer) have been determined by a gradient optimization that eliminates the basis set superposition error (BSSE) by the counterpoise (CP) method in each gradient cycle as well as by the standard gradient optimization. Both optimization methods lead to different potential energy surfaces (PES). The difference depends on the theoretical level used and is larger if correlation energy is considered. Intermolecular distances from the CP-corrected PES are consistently longer, and this difference might be significant (∼0.1 ?); also angular characteristics determined from both surfaces differ significantly. Different geometries were obtained even when passing to larger basis sets (aug-cc-pVDZ). The standard optimization procedure can result in a completely wrong structure. For example, the “quasi-linear” structure of the (HF)₂ (global minimum) does not exist at the standard MP2/ 6-31G** PES (where only cyclic structure was detected) and is found only at the CP-corrected PES. Stabilization energies obtained from the CP-corrected PES are always larger than these from the standard PES where the BSSE is added only a posteriori for the final optimized structure; both energies converge only when passing to a larger basis set (aug-cc-pVDZ). Received: 11 March 1998 / Accepted: 19 June 1998 / Published online: 4 September 1998 RID=" ID=" <E6>Acknowledgements.</E6> The project was supported by the Grant Agency of the Czech Republic (Grant No. 203/98/1166). RID=" ID=" <E5>Correspondence to</E5>: P. Hobza     

应用扩展的重叠分辨率法优化液相色谱条件分离二茂铁类化合物

本文对优化液相色谱分离条件的重叠分辨率法进行了改进与扩展：首先利用全范围的二元线性梯度淋洗进行溶剂强度的优化，再通过重叠分辨率法进行了溶剂选择性的优化，并且建立了计算机系统，使此方法程序化。利用此方法建立了分离Ｎ－二茂铁甲酰基－Ｎ＇－芳基硫脲类衍生物的最佳流动相体系。     

Parameter identification of thermophilic anaerobic degradation of valerate

The considered mathematical model of the decomposition of valerate presents three unknown kinetic parameters, two unknown stoichiometric coefficients, and three unknown initial concentrations for biomass. Applying a structural identifiability study, we concluded that it is necessary to perform simultaneous batch experiments with differenitial conditions for estimating these parameters. Four simultaneous batch experiments were conducted at 55°C, characterized by four different initial acetate concentrations. Product inhibition of valerate degradation by acetate was considered. Practical identification was done optimizing the sum of the multiple determination coefficients for all measured state viariables and for all experiments simultaneously. The estimated values of kinetic parameters and stoichiometric coefficients were characterized by the parameter correlation matrix, the confidence interval, and the student's t-test at 9% significance level with positive results except for the saturation constant, for which more eperiments for improving its identifiability should be conducted. In this article, we discussekinetic parameter estimation methods.     

Data loci in algebraic optimization

We consider parametric optimization problems from an algebraic viewpoint. The idea is to find all of the critical points of an objective function thereby determining a global optimum. For generic parameters (data) in the objective function the number of critical points remains constant. This number is known as the algebraic degree of an optimization problem. In this article, we go further by considering the inverse problem of finding parameters of the objective function so it gives rise to critical points exhibiting a special structure. For example if the critical point is in the singular locus, has some symmetry, or satisfies some other algebraic property. Our main result is a theorem describing such parameters.     

A greedy algorithm for finding maximum spanning trees in infinite graphs

In finite graphs, greedy algorithms are used to find minimum spanning trees (MinST) and maximum spanning trees (MaxST). In infinite graphs, we illustrate a general class of problems where a greedy approach discovers a MaxST while a MinST may be unreachable. Our algorithm is a natural extension of Prim's to infinite graphs with summable and strictly positive edge weights, producing a sequence of finite trees that converge to a MaxST.     

Finding normal solutions in piecewise linear programming

Letf: ⁿ (–, ] be a convex polyhedral function. We show that if any standard active set method for quadratic programming (QP) findsx(t)= arg min _x ¦x¦²/2+t f(x) for somet> 0, then its final working set defines a simple equality QP subproblem, whose Lagrange multiplier can be used both for testing ift is large enough forx(t) to coincide with the normal minimizer off, and for increasingt otherwise. The QP subproblem may easily be solved via the matrix factorizations used for findingx(t). This opens up the way for efficient implementations. We also give finite methods for computing the whole trajectory {x(t)} _t 0, minimizingf over an ellipsoid, and choosing penalty parameters inL ₁QP methods for strictly convex QP.This research was supported by the State Committee for Scientific Research under Grant 8S50502206.     

A two-phase algorithm for the partial accessibility constrained vehicle routing problem

In the partial accessibility constrained vehicle routing problem, a route can be covered by two types of vehicles, i.e. truck or truck + trailer. Some customers are accessible by both vehicle types, whereas others solely by trucks. After introducing an integer programming formulation for the problem, we describe a two-phase heuristic method which extends a classical vehicle routing algorithm. Since it is necessary to solve a combinatorial problem that has some similarities with the generalized assignment problem, we propose an enumerative procedure in which bounds are obtained from a Lagrangian relaxation. The routine provides very encouraging results on a set of test problems.     

Computing global minima to polynomial optimization problems using Gröbner bases

The local optimality conditions to polynomial optimization problems are a set of polynomial equations (plus some inequality conditions). With the recent techniques of Gröbner bases one can find all solutions to such systems, and hence also find global optima. We give a short survey of these methods. We also apply them to a set of problems termed with exact solutions unknown in the problem sets of Hock and Schittkowski. To these problems we give exact solutions.     

Approximate necessary conditions for locally weak Pareto optimality

Necessary conditions for a given pointx ⁰ to be a locally weak solution to the Pareto minimization problem of a vector-valued functionF=(f ₁,...,f _m ),F:XR ^m,XR ^m, are presented. As noted in Ref. 1, the classical necessary condition-conv {Df ₁(x ⁰)|i=1,...,m}T ^*(X, x ⁰) need not hold when the contingent coneT is used. We have proven, however, that a properly adjusted approximate version of this classical condition always holds. Strangely enough, the approximation form>2 must be weaker than form=2.The authors would like to thank the anonymous referee for the suggestions which led to an improved presentation of the paper.