Basis null space

Siace the columns of any complete B-matrix are a basis for the null space of the dimensional matrix, it follows that any two complete B-matrices are related by a nonsingular transformation. In other words, a complete B-matrix itself contains enough information as to which linear combiaations should be formed to obtain the optimized ones. Based on this observation, an efficient algorithm for the generation of an optimized complete B-matrix has been presented (22). No attempt is made here to demonstrate the algorithm. Instead, an example is being used to illustrate the results. 

Next we want to consider the fact that a stoichiometric number matrix can be calculated from the conservation matrix, and vice versa. Since Av = 0, the A matrix can be used to calculate a basis for the stoichiometric number matrix v. The stoichiometric number matrix v is referred to as the null space of the A matrix. When the conservation matrix has been row reduced it is in the form A = [/C,Z], where /c is an identity matrix with rank C. A basis for the null space is given by... 

Mathematica is very useful for carrying out these matrix operations. The operation for row reduction is RowReduce, and the operation for calculating a basis for the null space is NullSpace. Row reduction is also used to determine whether the equations in a set of conservation equations or reaction equations are independent. Rows that are dependent come out as all zeros when this is done, and they must be deleted because they do not provide any useful information. 

A basis for the null space v of conservation matrix 5.2-5 at specified pH obtained with equation 5.1-19 or with a computer is... 

Null space. If the product of two matrices is a zero matrix (all zeros), ax = 0 is said to be a homogeneous equation. The matrix jc is said to be the null space of a. Tn Mathematica a basis for the null space of a can be calculated by use of Null Space [a]. There is a degree of arbitrariness in the null space in that it provides a basis, and alternative forms can be calculated from it, that are equivalent. See Equation 5.1-19 for a method to calculate a basis for the null space by hand. When a basis for the null space of a matrix needs to be compared with another matrix of the same dimensions, they are both row reduced. If the two matrices have the same row-reduced form, they are equivalent. 

Transformation matrix. When the conservation matrix a for a system is written in terms of elemental compositions, the elements are used as components. But we can change the choice of components (change the basis) by making a matrix multiplication that does not change the row-reduced form of the a matrix or its null space. Since components are really coordinates, we can shift to a new coordinate system by multiplying by the inverse of the transformation matrix between the two coordinate systems. A new choice of components can be made by use of a component transformation matrix m, which gives the composition of the new components (columns) in terms of the old components (rows). The following matrix multiplication yields a new a matrix in terms of the new components. 

The apparent stoichiometric number matrix v" can be obtained from the row-reduced form of A" by use of the analogue of equation 5.1-19 or by calculating a basis for the null space using a computer program. 

More generally, if R is a matrix that contains a basis for the right null space of (i.e., SR = 0), then... 

Here, the matrix S of Equation (9.4) has a one-dimensional right null space, for which the vector [1 1 1 r is a basis. Equation (9.22) corresponds to summing the reaction potentials about the closed loop formed by the reactions in Equation (9.3). 

Famili, I. and Palsson, B.O. (2003) The convex basis of the left null space of the stoichiometric matrix leads to the definition of metabolically meaningful pools. Biophys. J. 85, 16-26. 

CxN) NxR) = CxR. Equation 7.1-8 is useful because it makes it possible to calculate a stoichiometric matrix from a conservation matrix. This operation is called taking the null space of A, and the Mathematica operation for doing this is called NullSpace. The use of NullSpace yields a basis for the stoichiometric number matrix. We will see what this means and how it is handled. 

This is why calculating the null space yields a basis for the conservation matrix. The row reduced conservation matrix can be labelled as follows ... 

An alternative to the method of Lagrange multipliers for imposing the necessary constraints is sketched below. It derives a lower dimensional unconstrained problem from the original constrained problem by using an orthogonal basis for the null space of the constraint matrix. This method is well suited to the potentially rank-deficient problem at hand, where steps may be taken to... 

Of course, first we need to find an initial minimal DD pair. Following the null space approach [23, 24], we compute a basis of the kernel of the stoichiometric matrix S. More specifically, we compute a column-reduced echelon form of the basis and (after a permutation of rows) obtain... 

Clearly, the first p-t-k rows of the corresponding matrices are non-negative. This fact is the basis for the binary null space approach [23, 24] for which efficient implementations are available [20,21]. Still, the DD suffers from the combinatorial explosion of intermediate extreme rays resulting from the combination of adjacent rays W, r>. Thus, a full EFM A is only applicable to medium-scale metabolic networks. 

Here, mattix A is the stoichiometric coefficient matrix formed from the reaction stoichiometry, and C is any concentration vector that satisfies Equation 6.7. Determination of the subspace of concenttations orthogonal to the stoichiometric subspace, is in fact equivalent to computing the mil space of A. We shall denote by N the mattix whose columns form a basis for the null space of A. Linear combination of the columns in N hence generates the set of concentrations orthogonal to the stoichiometric subspace. 

Hence, Xj = -2 if X2 = 1 and X3 = 0. We have therefore found a combination of values for Xj, X2, and X3 that produce zero when multiplied by A. In the language of linear algebra, we have found a basis for the null space of A. The vector... 

Determine the matrix N—the matrix with columns that form a basis for the null space of A. ... 

Finally, one may compute the set of vectors orthogonal to the stoichiometric subspace. Again computation gives a set of two vectors that form a basis for the null space. Two representative vectors for this space are given below ... 

From Section 7.2.1.1, the dimension of the AR is three (d=3) and there are four components (n = 4). It is expected that the dimension of the subspace orthogonal to the stoichiometric subspace is (4 - 3) = 1. Therefore, for the three-dimensional Van de Vusse system, the null space is given by a one-dimensional subspace (in other words, the basis for the null space is composed of a single, linearly independent vector). This is confirmed when the null space of A is computed, giving... 

An alternative method for computing S exists that relies on computing a basis for the null space of the stoichiometric matrix A. Vectors belonging to null (A ) are normal to S, and these vectors may be related to bounding hyperplane constraints that define S. This method is described in Feinberg (1987,2000a, 2000b). 

The null space KerB is a vector subspace of the whole space of Af-vectors. Let us have M-L linearly independent vectors (transposed row vectors) oej e KerB (A = 1, —, M-L), thus a basis of KerB . The basis can be completed by some L linearly independent vectors, say Bk (A = 1, —, L) to a basis of. Then the matrix... 

LDR Decomposition of G, Computation of an Orthogonal Basis of the Null Space of G Using Gaussian Elimination... 

The ssf inverse is independent of the particular basis V of the null space of G. 

Every n vector can be represented as a point in an -dimensional coordinate space. The n elements of the vector are the coordinates along n basis vectors, such as defined in the previous section. The null vector 0 defines the origin of the coordinate space. Note that the origin together with an endpoint define a directed line segment or axis, which also represents a vector. Hence, there is an equivalence between points and axes, which can both be thought as geometrical representations of vectors in coordinate space. (The concepts discussed here are extensions of those covered previously in Sections 9.2.4 to 9.2.5.)... 

TABLE VI. Values (in a.u.) of the non-null Ijij coefficients of the Cl energy increment for the first MC iteration of the Li2 molecule computed with the same basis set as in Table IV with a Cl space of 20 configurations. The optimum computed energy increment (in a.u.) is also given. For the first two sweeps all the rotations giving energy increments greater than 1.0x1 O 6 a.u. are included. Only the 5-1 rotation is shown in the last sweeps. 

Symmetry considerations alone can teach a lot about the nature and properties of crystal orbitals. In a first approximation, each orbital in the unit cell can be considered as giving rise in k-space to a separate energy band, which may be vacant, half-filled or filled, as the original MO is, and whose width and slope can be qualitatively guessed on the basis of simple symmetry and overlap arguments. For example, in a row of n-stacked ethylene molecules the bands arising from s-a and p-a molecular orbitals must be very flat because of scarce or null overlap, while the n-n band would develop some width because of jt-overlap (Fig. 6.4). Since all the s and the p-n MO are fully occupied, so are also the s and the p-n bands, and the wide gap to the lowest unoccupied band, the p-n band, explains the insulating nature of this hypothetical material. 

---