Electron density optimization

Fig. 2 Atomic shell structure as it emerges from electron-density optimization on a golden spiral. The variable convergence angle of An/ In — 1) manifests in the appearance of 2n — 1 additional cycles s, p, d, /) in each interval between Bohr levels n and n — I, shown here as elementary ripples. In contrast to the Bohr-Schrodinger (BS) model, closed shells in the Ford-circle simulation (FC) invariably coincide with noble-gas configurations...

Davidsou-Fletcher-Powell (DFP) a geometry optimization algorithm De Novo algorithms algorithms that apply artificial intelligence or rational techniques to solving chemical problems density functional theory (DFT) a computational method based on the total electron density... 

Once you have calculated an ab initio or a semi-empirical wave function via a single point calculation, geometry optimization, molecular dynamics or vibrations, you can plot the electrostatic potential surrounding the molecule, the total electronic density, the spin density, one or more molecular orbitals /i, and the electron densities of individual orbitals You can examine orbital energies and select orbitals for plotting from an orbital energy level diagram. 

However, it is also impirrtant to emphasize that orbitals are actually mathematical conv enience.s and not physical quantities (despite how real models may make them seem). While the energy, electron density, and optimized geonipctry are physical obseivables, the orbitals are not. In fact, several different sets of orbitals can lead to the same energy. Nevertheless, orbitals are very irseful in qualitative descriptions of bonding and reactivity. 

Optimize the structure of acetyl radical using the 6-31G(d) basis set at the HF, MP2, B3LYP and QCISD levels of theory. We chose to perform an Opt Freq calculation at the Flartree-Fock level in order to produce initial force constants for the later optimizations (retrieved from the checkpoint file via OptsReadFC). Compare the predicted spin polarizations (listed as part of the population analysis output) for the carbon and oxygen atoms for the various methods to one another and to the experimental values of 0.7 for the C2 carbon atom and 0.2 for the oxygen atom. Note that for the MP2 and QCISD calculations you will need to include the keyword Density=Current in the job s route section, which specifies that the population analysis be performed using the electron density computed by the current theoretical method (the default is to use the Hartree-Fock density). 

The aromatic portion of the molecules discussed in this chapter is frequently, if not always, an essential contributor to the intensity of their pharmacological action. It is, however, usually the aliphatic portion that determines the nature of that action. Thus it is a common observation in the practice Ilf medicinal chemistry that optimization of potency in these drug classes requires careful attention to the correct spatial orientation of the functional groups, their overall electronic densities, and the contribution that they make to the molecule s solubility in biological fluids. These factors are most conveniently adjusted by altering the substituents on the aromatic ring. 

Energy minimization methods that exploit information about the second derivative of the potential are quite effective in the structural refinement of proteins. That is, in the process of X-ray structural determination one sometimes obtains bad steric interactions that can easily be relaxed by a small number of energy minimization cycles. The type of relaxation that can be obtained by energy minimization procedures is illustrated in Fig. 4.4. In fact, one can combine the potential U r) with the function which is usually optimized in X-ray structure determination (the R factor ) and minimize the sum of these functions (Ref. 4) by a conjugated gradient method, thus satisfying both the X-ray electron density constraints and steric constraint dictated by the molecular potential surface. 

Two models of practical interest using quantum chemical parameters were developed by Clark et al. [26, 27]. Both studies were based on 1085 molecules and 36 descriptors calculated with the AMI method following structure optimization and electron density calculation. An initial set of descriptors was selected with a multiple linear regression model and further optimized by trial-and-error variation. The second study calculated a standard error of 0.56 for 1085 compounds and it also estimated the reliability of neural network prediction by analysis of the standard deviation error for an ensemble of 11 networks trained on different randomly selected subsets of the initial training set [27]. 

Petersilka, M., Gossmann, U. J., Gross, E. K. U., 1998, Time Dependent Optimized Effective Potential in the Linear Response Regime in Electronic Density Functional Theory. Recent Progress and New Directions, Dobson, J. F., Vignale, G., Das, M. P. (eds.), Plenum Press, New York. 

BSSE also opposes the tendency of the Hartree-Fock model to keep the interacting closed shell fragments too far apart. So, when optimized geometries are considered for the complex, BSSE is found to mimic some of those effects on the electron density distribution which would be induced by the interfragment dispersion contributions. 

Figure 13-12. The snapshots of the S2 (green) and S3 (cyan) binding states from the PMF calculation (top) and the results from the QM/MM X-ray optimization (bottom) superimposed on top of the crystal structure 1ZB6 (grey) and the electron density contour at 0.5a...

The first step in our procedure is to compute an optimized structure for each molecule and then to use this geometry to compute the electronic density and the electrostatic potential. A large portion of our work in this area has been carried out at the SCF/STO-5G //SCF/STO-3G level, although some other basis sets have also been used. We then compute V(r) on 0.28 bohr grids over molecular surfaces defined as the 0.001 au contour of the electronic density (Bader et al. 1987). The numbers of points on these grids are converted to surface areas (A2), and the and Fs min are determined. Our statistically based interaction in-... 

In contrast, the NBO and NRT methods make no use of molecular geometry information (experimental or theoretical), but instead provide optimal descriptions of orbital composition or electron-density distributions based directly on the first-order density operator. For this reason the NBO/NRT indices have predictive utility for a broad range of chemical phenomena, without bias toward geometry or other particular empirical properties. 

---