# What are the methods that can be used to ensure that an optimized geometry is a local minimum?

Usually it is suggested to do a phonon calculation. What are the other ways?

• +1. Welcome to the site! How about just perturbing the geometry and seeing if the energy goes down? How about running a geometry optimization calculation with the current geometry as a starting point? Since you say "phonon" instead of "vibrational frequency", I guess you are dealing with a periodic lattice? – Nike Dattani Aug 1 '20 at 15:20
• If this in in the context of DFT packages like Quantum ESPRESSO or VASP, geometry optimization is always done to a local minimum. You can't find the global minimum in regular Molecular dynamics-optimization techniques. – Xivi76 Aug 1 '20 at 18:36
• @Xivi76 -- strictly speaking, that is not correct. There is no guarantee that the optimized geometry is a local minimum. One can hope that it is often the case, but the only way to know for sure is to confirm all real vibrational modes. – Andrew Rosen Aug 2 '20 at 7:35
• @AndrewRosen Poor choice of words on my part. What I meant to say is, when one uses regular MD optimization, you can only expect to get a local minimum (this need not be achieved at all, as you said) but you can't guarantee a global minimum because it is only as good as the initial guess - if your initial guess is poor and you land up in a 'smaller' valley, there's no way to climb out of it and hence you can't find a global minimum. Please correct me if I'm wrong with this evaluation. – Xivi76 Aug 2 '20 at 16:16
• @Xivi76 but you can also converge to a saddle point geometry, in which case the gradient is still zero but you have one or more negative Hessian eigenvalues. – Susi Lehtola Aug 3 '20 at 8:48

This depends on what you are studying. For molecular systems without periodicity, the simplest approach is to carry out a vibrational frequency analysis and confirm that there are no imaginary modes. It is considered standard to carry out a vibrational frequency analysis for all investigated structures, provided the number of investigated systems is not prohibitively large.

For extended solids, which it sounds like you're modeling, it is much less common to confirm the result of a structure relaxation is indeed a minimum, although it can't hurt and is always a good idea if feasible. A phonon calculation can indeed be used to confirm your location in the potential energy surface.

There are other options though for extended solids. If you are modeling the adsorption of a molecule onto a surface, you can likely assume that the surface atoms are near or at their minimum energy positions (if you have pre-relaxed the adsorbate-free structure). In this case, you could perform a vibrational frequency analysis of just the adsorbate atoms (or, even better, the adsorbate atoms and a few nearby surface atoms) to confirm there are no imaginary modes. ASE has an excellent vibrational analysis function that works with many of the most popular DFT packages, which I recommend if you're using a code like VASP. Beyond this, you can always try slightly modifying the atomic positions of the converged structure in several different ways, re-optimizing each structure to see if the same low-energy configuration is found. Of course, this does not guarantee you have arrived at a local minimum, but it can provide further support.

• What if you get a imaginary frequency from phonon calculation with a perfectly optimized structure (with large basis sets and tight criterions). Does this mean that this structure not stable at all? If so how all different initial geometries end up at the same geometry. – Alfred Aug 3 '20 at 23:31
• @alfred yes, imaginary frequencies means the system is at a maximum of the local parabolic energy surface with respect to at least one perturbation -- so it is mechanically unstable. – Phil Hasnip Nov 13 '20 at 14:16

In order to check if a geometry is a local minimum, it is a necessary and sufficient condition that the Hessian is positive (semi)definite, i.e. that the lowest eigenvalue of the nuclear Hessian is non-negative.

Namely, expanding the energy $$E({\bf R})$$ around the reference point $${\bf R}_0$$ you have the Taylor expansion $$E({\bf R}) = E({\bf R}_0) - {\bf g} \cdot ({\bf R}-{\bf R}_0) + \frac 1 2 ({\bf R}-{\bf R}_0)^{\rm T} \cdot {\bf H} \cdot ({\bf R}-{\bf R}_0) + \mathcal{O}(|({\bf R}-{\bf R}_0)|^3)$$ where $${\bf g} = -[\nabla E({\bf R})]_{{\bf R}={\bf R}_0}$$ is the gradient and $$H_{ij} = -\partial^2 E / \partial R_i \partial R_j$$ is the nuclear Hessian, which has the size $$3N_{\rm atoms} \times 3N_{\rm atoms}$$.

At an extremum $${\bf g}={\bf 0}$$; this is what the optimizer finds for you. However, in some cases you can find that the Hessian has negative eigenvalues $$\lambda_i$$ with corresponding eigenvectors $$\boldsymbol{\rho}_i$$. Now, if you set $${\bf R}={\bf R}_0+\epsilon \boldsymbol{\rho}_i$$, you will find that $$E(\epsilon)=E({\bf R})=E({\bf R}_0) + \frac 1 2 \epsilon^2 \lambda_i + \mathcal{O}(|\epsilon|^3)$$. Since $$\lambda_i<0$$, this means you can find a lower energy by moving the atoms slightly along $$\boldsymbol{\rho}_i$$.

A frequency calculation means computing all eigenvalues of the Hessian, which becomes costly in large systems. (I'm not sure what is done in phonon calculations.) However, since you only need the lowest eigenvalue of the Hessian to check if you are at a local minimum, iterative diagonalizers such as the Davidson method can be used to focus only on the lowest roots of the nuclear Hessian.

Alternatively, there are also some methods that need only gradient information to estimate the lowest eigenvalue of the nuclear Hessian; a "dimer method" has been suggested e.g. by Henkelman and Jónsson in J. Chem. Phys. 111, 7010 (1999).