Nonsingular matrix Parameters

In other words, if the C coefficients /iH+R+i /is are given the values determined by Eq. (19), then the total of the expressions in (17) will be a direct mechanism. Furthermore, if we go through this procedure for every selection of C columns in (17) such that the C x C matrix M is nonsingular, then we get every direct mechanism for the overall reaction (14). Altogether there are R + C undetermined coefficients fiH+15. .., fis in (17), the last C of which are determined for each direct mechanism. The remaining R parameters fiH + ,..., /iH + R are in the expression (14) for the overall reaction, which is of multiplicity R. Similarly each direct mechanism must be a function of these R parameters. 

When written with the help of the Tl matrix as in (19), from (20) the OR parameter and other linear response properties are seen to afford singularities where co = coj, just like in the SOS equation (2). Therefore, at and near resonances the solutions of the TDDFT response equations (and response equations derived for other quantum chemical methods) yield diverging results that cannot be compared directly to experimental data. In reality, the excited states are broadened, which may be incorporated in the formalism by introducing dephasing constants 1 such that o, —> ooj — iT j for the excitation frequencies. This would lead to a nonsingular behavior of (20) near the coj where the real and the imaginary part of the response function varies smoothly, as in the broadened scenario at the top of Fig. 1. 

Consider a system for which K has been subjected to a similarity transformation to give a system with a coefficient matrix P KP, where P is nonsingular. Recall that under a similarity transformation, the eigenvalues do not change. Impose on P KP all the structural constraints on K and require that the response function of the system with matrix P" KP be the same as that of the system with matrix K. If the only P that satisfies those requirements is the identity matrix, all parameters are globally identifiable. If a P I satisfies the requirements, one can work out which parameters are not identifiable and which are. 

The final step is to find the pairwise correlations for the identifiable parameters. Even if identifiable it may be difficult to estimate two parameters separately in the presence of measurement error if they are highly correlated. To check that, we reduce the matrix gj gs by eliminating rows and columns corresponding to the nonidentifiable parameters to obtain gi gi- Since gi gi is the matrix corresponding to the identifiable parameters, it is nonsingular so we can invert it to obtain (gi gi)The correlation matrix is obtained by dividing the i,j element of (g/gi) by the square root of the product of the /th and yth diagonal elements. 

The determinant for the Rank matrix is nonsingular thus, the Rank of this Dimension matrix is 4. The number of dimensionless parameters is... 

In any Newton-based optimization - which, as discussed in Section 10.8, implicitly or explicitly requires the inversion of the Hessian matrix - the inclusion of redundant parameters is not only unnecessary but also undesirable since, at stationary points, these parameters make the electronic Hessian singular. The singularity of the Hessian follows from (10.2.8), which shows that the rows and columns corresponding to redundant rotations vanish at stationary points. Away from the stationary points, however, the Hessian (10.1.30) is nonsingular since the gradient elements that couple the redundant and nonredundant operators in (10.2.5) do not vanish. Still, as the optimization approaches a stationary point, the smallest eigenvalues of the Hessian will tend to zero and may create convergence problems as the stationary point is approached. Therefore, for the optimization of a closed-shell state by a Newton-based method, we should consider only those rotations that mix occupied and virtual orbitals ... 

---