Penalty-Barrier Function Methods

When equality constraints are also present, some recent algorithms derived from interior point methods for linear programming use a barrier method that generates steps by solving the following minimization problem  

It is possible to use a penalty function for the equality constraints and a barrier function for the inequality constraints. The following function represents one of 

---

Penalty function methods, barrier function methods, and hybrid penalty-barrier function methods are discussed in Chapter 12 they are implemented and adopted... 

Barrier function methods are very similar to penalty function methods except that they start at an interior point of the feasible region and set a barrier against leaving the feasible region. In this case, the feasible region must have an interior, so this method is generally restricted to inequality constraints. Consider the nonlinear problem with inequality constraints. 

Murray, W., and Wright, M., Projected Lagrangian methods based on trajectories of barrier and penalty functions, SOL Report 78-23, Stanford University, Stanford, California (1978). 

The chiral discrimination in the self-association of chiral l,3a,4,6a-tetrahydroi-midazo[4,5-d]imidazoles 3 has been studied using density functional theory methods [37], (Scheme 3.20). Clusters from dimers to heptamers have been considered. The heterochiral dimers (RR SS or SS RR) are more stable than the homochiral ones (RR RR or SS SS) with energy differences up to 17.5 kJ mol-1. Besides, in larger clusters, the presence of two adjacent homochiral molecules imposes an energetic penalty when compared to alternated chiral systems (RR SS RR SS...). The differences in interaction energy within the dimers of the different derivatives have been analyzed based on the atomic energy partition carried out within the AIM framework. The mechanism of proton transfer in the homo- and heterochiral dimers shows large transition-state barriers, except in those cases where a third additional molecule is involved in the transfer. The optical rotatory power of several clusters of the parent compound has been calculated and rationalized based on the number of homochiral interactions and the number of monomers of each enantiomer within the complexes. 

---