Madelung model

It is sometimes desirable to include the effect of the rest of the system, outside of the QM and MM regions. One way to do this is using periodic boundary conditions, as is done in liquid-state simulations. Some researchers have defined a potential that is intended to reproduce the effect of the bulk solvent. This solvent potential may be defined just for this type of calculation, or it may be a continuum solvation model as described in the next chapter. For solids, a set of point charges, called a Madelung potential, is often used. 

Fig. 3.3 shows a similar ionic-model construction of the electronic energies for the isostructural transition-metal oxide MnO. In this case, the Madelung energy lifts the Mn 3d level above the top of the... 

The results for the three models are summarized in Table 9.2, which lists the effective Madelung constant n (t, defined by eiectrosuiic = — peff/a the effective Born coefficient, defined by U = — NAncff/a(l — l/neff) and the lattice energy U. 

In Chapter 2 it was shown that the Madelung field of a crystal is equivalent to a capacitive electric circuit which can be solved using a set of Kirchhoff equations. In Sections 3.1 and 3.2 it was shown that for unstrained structures the capacitances are all equal and that there is a simple relationship between the bond flux (or experimental bond valence) and the bond length. These ideas are brought together here in a summary of the three basic rules of the bond valence model, Rules 3.3, 3.4, and 3.5. 

We see that the simple rectangular d band model reproduces the behaviour found by experiment and predicted by Miedema s semi-empirical scheme. However, we must stress that the model does not give credence to any theory that bases the heat of formation of transition-metal alloys on ionic Madelung contributions that arise from electronegativity differences between the constituent atoms because in the metallic state the atoms are perfectly screened and, hence, locally charge neutral. Instead, the model supports... 

The crystal lattice energy can be estimated from a simple electrostatic model When this model is applied to an ionic crystal only the electrostatic charges and the shortest anion-cation intermiclear distance need be considered. The summation of all the geometrical interactions be/Kveeti the ions is called the Madelung constant. From this model an equatitWjor the crystal lattice energy is derived ... 

In many cases more refined molecular orbital models give a better agreement between theory and experiment. Self-consistent field- 80 81> or CNDO-82> calculations as well as other ab initio calculations 83>84) were performed and the results of several different approaches for phosphorous compounds were critically evaluated by M. Pelavin 77). For the nitrogen compounds two different linear relationships, one for cations and for neutral molecules, the other for anoins, were observed 82), a phenomenon which might be explained with a Madelung contribution. 

---