The Least Equation

The general form for a matrix polynomial equation satisfied by A is 

The least equation is the polynomial equation satisfied by A that has the smallest possible degree. There is only one least equation 

The degree of the least equation, k, is called the rank of the matrix A. The degree k is never greater than n for the least equation (although there are other equations satisfied by A for which k n). If A = n, the size of a square matrix, the inverse A exists. If the matrix is not square or k n, then A has no inverse. 

One method of finding the least equation for the simple second degree case is illustrated. Find a number r such that 

Use the method given above to find the least equation of the matrix 

---

Note that the matrix from Exercise 2-8 is the matrix of coefficients in this simultaneous equation set. Note also the similarity in method between finding the least equation and Gaussian elimination. 

While Eq. III-18 has been verified for small droplets, attempts to do so for liquids in capillaries (where Rm is negative and there should be a pressure reduction) have led to startling discrepancies. Potential problems include the presence of impurities leached from the capillary walls and allowance for the film of adsorbed vapor that should be present (see Chapter X). There is room for another real effect arising from structural peiturbations in the liquid induced by the vicinity of the solid capillary wall (see Chapter VI). Fisher and Israelachvili [19] review much of the literature on the verification of the Kelvin equation and report confirmatory measurements for liquid bridges between crossed mica cylinders. The situation is similar to that of the meniscus in a capillary since Rm is negative some of their results are shown in Fig. III-3. Studies in capillaries have been reviewed by Melrose [20] who concludes that the Kelvin equation is obeyed for radii at least down to 1 fim. 

The representation of trial fiinctions as linear combinations of fixed basis fiinctions is perhaps the most connnon approach used in variational calculations optimization of the coefficients is often said to be an application of tire linear variational principle. Altliough some very accurate work on small atoms (notably helium and lithium) has been based on complicated trial functions with several nonlinear parameters, attempts to extend tliese calculations to larger atoms and molecules quickly runs into fonnidable difficulties (not the least of which is how to choose the fomi of the trial fiinction). Basis set expansions like that given by equation (A1.1.113) are much simpler to design, and the procedures required to obtain the coefficients that minimize are all easily carried out by computers. 

The simplest extension to the DH equation that does at least allow the qualitative trends at higher concentrations to be examined is to treat the excluded volume rationally. This model, in which the ion of charge z-Cq is given an ionic radius d- is temied the primitive model. If we assume an essentially spherical equation for the u. . 

The least recognized fonns of the Porod approximation are for the anisotropic system. If we consider the cylindrical scattering expression of equation (B 1.9.61). there are two principal axes (z and r directions) to be discussed... 

Equations (169) and (171), together with Eqs. (170), fomi the basic equations that enable the calculation of the non-adiabatic coupling matrix. As is noticed, this set of equations creates a hierarchy of approximations starting with the assumption that the cross-products on the right-hand side of Eq. (171) have small values because at any point in configuration space at least one of the multipliers in the product is small [115]. 

Once the least-squares fits to Slater functions with orbital exponents e = 1.0 are available, fits to Slater function s with oth er orbital expon cn ts can be obtained by siin ply m ii Itiplyin g th e cc s in th e above three equations by It remains to be determined what Slater orbital exponents to use in electronic structure calculation s. The two possibilities may be to use the "best atom" exponents (e = 1. f) for II. for exam pie) or to opiim i/e exponents in each calculation. The "best atom expon en ts m igh t be a rather poor ch oicc for mo lecular en viron men ts, and optirn i/.at ion of non linear exponents is not practical for large molecules, where the dimension of the space to be searched is very large.. 4 com prom isc is to use a set of standard exponents where the average values of expon en ts are optirn i/ed for a set of sin all rn olecules, fh e recom -mended STO-3G exponents are... 

The standard least-squares approach provides an alternative to the Galerkin method in the development of finite element solution schemes for differential equations. However, it can also be shown to belong to the class of weighted residual techniques (Zienkiewicz and Morgan, 1983). In the least-squares finite element method the sum of the squares of the residuals, generated via the substitution of the unknown functions by finite element approximations, is formed and subsequently minimized to obtain the working equations of the scheme. The procedure can be illustrated by the following example, consider... 

The basic procedure for the derivation of a least squares finite element scheme is described in Chapter 2, Section 2.4. Using this procedure the working equations of the least-squares finite element scheme for an incompressible flow are derived as follows ... 

Working equations of the least-squares scheme in Cartesian coordinate systems... 

It is evident that application of Green s theorem cannot eliminate second-order derivatives of the shape functions in the set of working equations of the least-sc[uares scheme. Therefore, direct application of these equations should, in general, be in conjunction with C continuous Hermite elements (Petera and Nassehi, 1993 Petera and Pittman, 1994). However, various techniques are available that make the use of elements in these schemes possible. For example, Bell and Surana (1994) developed a method in which the flow model equations are cast into a set of auxiliary first-order differentia] equations. They used this approach to construct a least-sciuares scheme for non-Newtonian flow equations based on equal-order C° continuous, p-version hierarchical elements. 

The degree of the least polynomial of a square matr ix A, and henee its rank, is the number of linearly independent rows in A. A linearly independent row of A is a row that eannot be obtained from any other row in A by multiplieation by a number. If matrix A has, as its elements, the eoeffieients of a set of simultaneous nonhomo-geneous equations, the rank k is the number of independent equations. If A = , there are the same number of independent equations as unknowns A has an inverse and a unique solution set exists. If k < n, the number of independent equations is less than the number of unknowns A does not have an inverse and no unique solution set exists. The matrix A is square, henee k > n is not possible. 

The least squares derivation for quadratics is the same as it was for linear equations except that one more term is canied through the derivation and, of course, there are three normal equations rather than two. Random deviations from a quadratic are ... 

The mathematical requirements for unique determination of the two slopes mi and ni2 are satisfied by these two measurements, provided that the second equation is not a linear combination of the first. In practice, however, because of experimental error, this is a minimum requirement and may be expected to yield the least reliable solution set for the system, just as establishing the slope of a straight line through the origin by one experimental point may be expected to yield the least reliable slope, inferior in this respect to the slope obtained from 2, 3, or p experimental points. In univariate problems, accepted practice dictates that we... 

Eirst decide what the integral equation corresponding to Eq. (6-29) is for the approximate wave function (6-30), then integrate it for various values of y. Report both Y at the minimum energy and Cmin for the Gaussian approximation function. This is a least upper bound to the energy of the system. Your report should include a... 

We may now understand the nature of the change which occurs when an anhydrous salt, say copper sulphate, is shaken with a wet organic solvent, such as benzene, at about 25°. The water will first combine to form the monohydrate in accordance with equation (i), and, provided suflScient anhydrous copper sulphate is employed, the effective concentration of water in the solvent is reduced to a value equivalent to about 1 mm. of ordinary water vapour. The complete removal of water is impossible indeed, the equilibrium vapour pressures of the least hydrated tem may be taken as a rough measure of the relative efficiencies of such drying agents. If the water present is more than sufficient to convert the anhydrous copper sulphate into the monohydrate, then reaction (i) will be followed by reaction (ii), i.e., the trihydrate will be formed the water vapour then remaining will be equivalent to about 6 mm. of ordinary water vapour. Thus the monohydrate is far less effective than the anhydrous compound for the removal of water. 

In order to parameterize a QSAR equation, a quantihed activity for a set of compounds must be known. These are called lead compounds, at least in the pharmaceutical industry. Typically, test results are available for only a small number of compounds. Because of this, it can be difficult to choose a number of descriptors that will give useful results without htting to anomalies in the test set. Three to hve lead compounds per descriptor in the QSAR equation are normally considered an adequate number. If two descriptors are nearly col-linear with one another, then one should be omitted even though it may have a large correlation coefficient. 

The Least Squares or Best-fit Line. The simplest type of approximating curve is a straight line, the equation of which can be written as in form number 1 above. It is customary to employ the above definition when X is the independent variable and Y is the dependent variable. 

---