@phdthesis{Henn2004,
  author    = {Henn, Julian},
  title     = {The electron density : a bridge between exact quantum mechanics and fuzzy chemical concepts},
  url       = {http://nbn-resolving.de/urn:nbn:de:bvb:20-opus-9003},
  school      = {Universit{\"a}t W{\"u}rzburg},
  year      = {2004},
  abstract  = {Summary The nature of the chemical bond is a topic under constant debate. What is known about individual molecular properties and functional groups is often taught and rationalized by explaining Lewis structures, which, in turn, make extensive use of the valence concept. The valence concept distinguishes between electrons, which do not participate in chemical interactions (core electrons) and those, which do (single, double, triple bonds, lone-pair electrons, etc.). Additionally, individual electrons are assigned to atomic centers. The valence concept is of paramount success: It allows the successful planning of chemical syntheses and analyses, it explains the behavior of individual functional groups, and, moreover, it provides the "language" to think of and talk about molecular structure and chemical interactions. The resounding success of the valence concept may be misleading to forget its approximative character. On the other hand, quantum mechanics provide in principle a quantitative description of all chemical phenomena, but there is no discrimination between electrons in quantum mechanics. From the quantum mechanical point of view there are only indistinguishable electrons in the field of the nuclei, i.e., it is impossible to assign a given electron to a particular center or to ascribe a particular purpose to individual electrons. The concept of indistinguishability of micro particles is founded on the Heisenberg uncertainty relation, which states, that wavepackets diverge in the 6N dimensional phase space, such that individual trajectories can not be identified. Hence it is a deep-rooted and approved physical concept. As an introduction to the present work density partitioning schemes were discussed, which divide the total molecular density into chemically meaningful areas. These partitioning schemes are intimately related to either the concepts of bound atoms in a molecule (as in the Atoms In Molecules theory (AIM) according to Bader or as in the Hirshfeld partitioning scheme) or to the concept of chemical structure in the sense of Lewis structures, which divide the total molecular density into core and valence density, where the valence density is split up again into bonding and non-bonding electron densities. Examples are early and recent loge theories, the topological analysis by means of the Electron Localization Function (ELF), and the Natural Bond Orbital (NBO) approach. Of these partitioning schemes, the theories according to Bader (AIM), to Becke and Edgecomb (ELF) and according to Weinhold (NBO and Natural Resonance Theory, NRT), respectively, were reviewed in detail critically. Points of criticism were explicated for each of the mentioned theories. Since theoretically derived electron densities are to be compared to experimentally derived densities, a brief introduction into the theory of X-ray di®raction experiments was given and the multipole formalism was introduced. The procedure of density refinement was briefly discussed. Various suggestions for improvements were developed: One strategy would be the employment of model parameters, which are to a maximum degree mutually orthogonal, with the object of minimizing correlations among the model parameters, e.g., to introduce nodal planes into the radial functions of the multipole model. A further suggestion involves the guidance of the iterative refinement procedure by an extremum principle, which states, that when di®erent solutions to the least squares minimization problem are available with about the same statistical measures of quality and with about the same residual density, then the solution is to prefer, which yields a minimum density at the bond critical point (BCP) and a maximum polarity in terms of the ratio of distances between the BCP and the nuclei. This suggestion is based on the well known fact, that the bond polarity (in terms of the ratio of distances between the BCP and the respective nuclei) is underestimated in the experiment. Another suggestion for including physical constraints is the explicit consideration of the virial theorem, e.g., by evaluating the integration of the Laplacian over the entire atomic basins and comparing this value to zero and to the value obtained from the integration of the electron gradient field over the atomic surface. The next suggestion was to explicitly use the electrostatic theorem of Feynman (often also denoted as Hellmann-Feynman theorem), which states, that the forces onto the nuclei can be calculated from the purely classical electrostatic forces of the electron distribution and the nuclei distribution. For a stationary system, these forces must add to zero. This also provides an internal quality criterion of the density model. This can be performed in an iterative way during the refinement procedure or as a test of the final result. The use of the electrostatic theorem is expected to reduce significantly correlations among static density parameters and parameters describing vibrations, since it is a valuable tool to discriminate between physically reasonable and artificial static electron densities. All of these mentioned suggestions can be applied as internal quality criteria. The last suggestion is based on the idea to initiate the experimental refinement with a set of model parameters, which is, as much as possible close to the final solution. This can be achieved by performing periodic boundary conditions calculations, from which theoretically created files are obtained, which contain the Miller indices (h, k, l) and the respective intensity I. This file is used for a model parameter estimation (refinement), which excludes vibrations. The resulting parameters can be used for the experimental refinement, where, in a first step, the density parameters are fixed to determine the parameters describing vibrations. For a fine tuning, again the electrostatic theorem and the other above mentioned suggestions could be applied. Theoretical predictions should not be biased by the method of computation. Therefore the dependence of the density analyzing tools on the level of calculation (method of calculation/basis set) and on the substituents in complex chemical bonding situations were evaluated in the second part of the present work. A number of compounds containing formal single and double sulfur nitrogen bonds was investigated. For these compounds, experimental data were also available. The calculated data were compared internally and with the experimental results. The internal comparison was drawn with regard to questions of convergency as well as with regard to questions of consistency: The resulting molecular properties from NBO/NRT analyses were found to be very stable, when the geometries were optimized at the respective level of theory. This stability is valid for variations in the methods of calculation as well as for variations in the basis set. Only the individual resonance weights of the contributing Natural Lewis Structures differed considerably depending on the level of calculation and depending on the substituents. However, the deviations were in both cases to a large extent within a limit which preserves the descending order of the leading resonance structure weights. The resulting bond orders, i.e., the total, covalent and ionic bond order from NRT calculations, were not affected by the shift in the resonance weights. The analysis of the bond topological parameters resulted in a discrimination between insensitive parameters and sensitive parameters. The stable parameters do neither depend strongly on the method of calculation nor on the basis set. Only minor variation occurs in the numerical values of these parameters, when the level of calculation is changed or even when other functional groups (H, Me, or tBu) are employed, as long as the methods of calculation do not drop considerably below a standard level. The bond descriptors of the sulfur nitrogen bonds were found to be also stable with respect to the functional groups R = H, R = Me, and R = tBu. Stable parameters are the bond distance, the density at the bond critical point (BCP) and the ratio of distances between the BCP and the nuclei A and B, which varies clearly when considering the formal bond type. For very small basis sets like the 3-21G basis set, this characteristic stability collapses. The sensitive parameters are based on the second derivatives of the density with respect to the coordinates. This is in accordance with the well known fact, that the total second derivative of the density with respect to the coordinates is a strongly oscillating function with positive as well as negative values. A profound deviation has to be anticipated as a consequence of strong oscillations. lambda3, which describes the local charge depletion in the direction of the interaction line, is the most varying parameter. A detailed analysis revealed that the position of the BCP in the rampant edge of the Laplacian distribution is responsible for the sensitivity of the numerical value of lambda3 in formal double bonds. Since the slope of the Laplacian assumes very high values in its rampant edge, a tiny displacement of the BCP leads already to a considerable change in lambda3. This instability is not a failure of the underlying theory, but it yields de facto to a considerable dependence of sensitive bond topological properties on the method of calculation and on the applied basis sets. Since the total second derivative is important to judge on the nature of the bond in the AIM theory (closed shell interactions versus shared interactions), the changes in lambda3 can lead to differing chemical interpretations. The comparison of theoretically derived bond topological properties of various sulfur nitrogen bonds provides the possibility to measure the self consistency of this data set. All data sets clearly exhibit a linear correlation between the bond distances and the density at the BCP on one hand and between the bond distances and the Laplacian values at the BCP on the other hand. These correlations were almost independent of the basis set size. In this context, the linear regression has to be regarded exclusively as a descriptive statistics tool. There is no correlation anticipated a priori. The formal bond type was found to be readily deducible from the theoretically obtained bond topological descriptors of the model systems. In this sense, the bond topological properties are self consistent despite of the numerical sensitivity of the derivatives, as exemplified above. Often, calculations are performed with the experimentally derived equilibrium geometries and not with optimized ones. Applying this approach, the computationally costly geometry optimizations are saved. Following this approach the bond topological properties were calculated using very flexible basis sets and employing the fixed experimental geometry (which, of course, includes the application of tBu groups). Regression coe±cients similar to those from optimized geometries were obtained for correlations between bond distances and the densities at the BCP as well as for the correlation between bond distances and the Laplacian at the BCP, i.e. the approach is valid. However, the data points scattered less and the coe±cient of correlation was clearly increased when geometry optimizations were performed beforehand. The comparison between data obtained from theory and experiment revealed fundamental discrepancies: In the data set of bond topological parameters from the experiment, the behavior of only 2 out of 3 insensitive parameters was comparable to the behavior of the theoretically obtained values, i.e. theoretical and experimental bond distances as well as theoretical and experimental densities at the BCP correlate. From the theoretically obtained data it was easy to deduce the formal bond type from the position of the BCP, since it changed in a systematic manner. The respective experimentally obtained values were almost constant and did not change systematically. For the SN bonds containing compounds, the total second derivative assumes exclusively negative values in the experiment. Due to the different internal behavior, experimentally and theoretically sensitive bond topological values could not be compared directly. The qualitative agreement in the Laplacian distribution, however, was excellent. In the third and last part of this work, the application to chemical systems follows. Formal hypervalent molecules, i.e. molecules where some atoms are considered to hold more than 8 electrons in their valence shell, were investigated. These were compounds containing sulfur nitrogen bonds (H(NtBu)2SMe, H2C{S(NtBu)2(NHtBu)}2, S(NtBu)2 and S(NtBu)3) and a highly coordinated silicon compound. The set of sulfur nitrogen compounds also contained a textbook example for valence expansion, the sulfur triimide. For these molecules, experimental reference values were available from high resolution X-ray experiments. The experimental results were in the case of the sulfur triimide not unique. Furthermore, from the experimental bond topological data no definite conclusion about the formal bonding type could be drawn. The situation of sulfur nitrogen bonds in the above mentioned set of molecules was analyzed in terms of a geometry discussion and by means of a topological analysis. The methyl-substituted isolated molecules served as model compounds. For the interpretation of the bonding situation additional NBO/NRT calculations were preformed for the sulfur nitrogen compounds and an ELF calculation and analysis was performed for the silicon compound. The ELF analysis included not only the presentation and discussion of the ELF-isosurfaces (eta = 0.85), but also the investigation of populations of disynaptic valence basins and the percentage contributions to these populations of the individual atoms when the disynaptic valence basins are split into atomic contributions according to Bader's partitioning scheme. The question of chemical interest was whether hypervalency is present in the set of molecules or not. In the first case the octet rule would be violated, in the second case Pauling's verdict would be violated. While the concept of hypervalency is well established in chemistry, the violation of Pauling's verdict is not. The quantitative numbers of the sensitive bond topological values from theory and experiment were not comparable, since no systematic relationship between the experimentally and theoretically determined sensitive bond descriptors was found. However, the insensitive parameters are in good agreement and the qualitative Laplacian distribution is, with few exceptions, in excellent agreement. The formal bonding type was deduced from experimental and theoretical topological data by considering the number and shape of valence shell charge concentrations in proximity to the sulfur and nitrogen centers. The results from NBO/NRT calculations confirmed the findings. All employed density analyzing tools AIM, ELF and NBO/NRT coincided in describing the bonding situation in the formally hypervalent molecules as highly polar. A comparison and analysis of experimentally and theoretically derived electron densities led consistently to the result, that regarding this set of molecules, hypervalency has to be excluded unequivocally.},
  subject      = {Elektronendichtebestimmung},
  language  = {en}
}