Structural properties, calculating

Figure 12.11 The class diagram of structure property calculator.

In spite of the great success of the computer simulation methods in the determination of the microscopic properties of the solutions, the capacity of the traditional MD and MC simulations is always limited by the choice of the suitable potential functions to describe the interatomic interactions. The potentials are most often checked by comparison of the structural properties calculated from the simulation with those determined experimentally. The reverse Monte Carlo (RMC) method, developed by McGreevy and Pusztai [41] does not rely upon knowledge of any interaction potential, instead it generates a large set of atomic configurations on the condition that the difference between the experimental and calculated structure functions (or pair-distribution functions) should be minimum. The same structural... 

We have developed a model to study the basic structural properties of solid C<,q. The model consists of two distinct types of intermolecular interactions. The dominant one is the van der Waals-type interactions between carbon atoms on different Cm molecules. A secondary short-range Coulomb interaction is modeled by a small charge transfer between the two types of bonds in the C60 molecule. In contrast to early calculations [6] which include the van der Waals interactions only, our model predicts correctly the observed cubic ground-state structure Pa3. Many structural properties calculated, such as the compressibility, cohesive energy, and specific heat, are in good agreement with experiments l7l. 

In conclusion, we have constructed a model for the in-termolecular interactions in solid C o. The basic structural properties calculated are in good agreement with experiments. It is shown that around T, 270 K the structure undergoes a first-order transition from the high-... 

Computational solid-state physics and chemistry are vibrant areas of research. The all-electron methods for high-accuracy electronic stnicture calculations mentioned in section B3.2.3.2 are in active development, and with PAW, an efficient new all-electron method has recently been introduced. Ever more powerfiil computers enable more detailed predictions on systems of increasing size. At the same time, new, more complex materials require methods that are able to describe their large unit cells and diverse atomic make-up. Here, the new orbital-free DFT method may lead the way. More powerful teclmiques are also necessary for the accurate treatment of surfaces and their interaction with atoms and, possibly complex, molecules. Combined with recent progress in embedding theory, these developments make possible increasingly sophisticated predictions of the quantum structural properties of solids and solid surfaces. 

For long-term simulations, it generally proves advantageous to consider numerical integrators which pass the structural properties of the model onto the calculated solutions. Hence, a careful analysis of the conservation properties of QCMD model is required. A particularly relevant constant of motion of the QCMD model is the total energy of the system... 

All the techniques described above can be used to calculate molecular structures and energies. Which other properties are important for chemoinformatics Most applications have used semi-empirical theory to calculate properties or descriptors, but ab-initio and DFT are equally applicable. In the following, we describe some typical properties and descriptors that have been used in quantitative structure-activity (QSAR) and structure-property (QSPR) relationships. 

A challenging task in material science as well as in pharmaceutical research is to custom tailor a compound s properties. George S. Hammond stated that the most fundamental and lasting objective of synthesis is not production of new compounds, but production of properties (Norris Award Lecture, 1968). The molecular structure of an organic or inorganic compound determines its properties. Nevertheless, methods for the direct prediction of a compound s properties based on its molecular structure are usually not available (Figure 8-1). Therefore, the establishment of Quantitative Structure-Property Relationships (QSPRs) and Quantitative Structure-Activity Relationships (QSARs) uses an indirect approach in order to tackle this problem. In the first step, numerical descriptors encoding information about the molecular structure are calculated for a set of compounds. Secondly, statistical and artificial neural network models are used to predict the property or activity of interest based on these descriptors or a suitable subset. 

The chirality code of a molecule is based on atomic properties and on the 3D structure. Examples of atomic properties arc partial atomic charges and polarizabilities, which are easily accessible by fast empirical methods contained in the PETRA package. Other atomic properties, calculated by other methods, can in principle be used. It is convenient, however, if the chosen atomic property discriminates as much as possible between non-equivalent atoms. 3D molecular structures are easily generated by the GORINA software package (see Section 2.13), but other sources of 3D structures can be used as well. 

Mujica A and R J Needs 1993. First-principles Calculations of the Structural Properties, Stability, aind Band Structure of Complex Tetrahedral Phases of Germanium ST12 and BC8. Physical Review B48 17010-17017. 

With all-atom simulations the locations of the hydrogen atoms are known and so the order parameters can be calculated directly. Another structural property of interest is the ratio of trans conformations to gauche conformations for the CH2—CH2 bonds in the hydrocarbon tail. The trans gauche ratio can be estimated using a variety of experimental techniques such as Raman, infrared and NMR spectroscopy. 

PW91 (Perdew, Wang 1991) a gradient corrected DFT method QCI (quadratic conhguration interaction) a correlated ah initio method QMC (quantum Monte Carlo) an explicitly correlated ah initio method QM/MM a technique in which orbital-based calculations and molecular mechanics calculations are combined into one calculation QSAR (quantitative structure-activity relationship) a technique for computing chemical properties, particularly as applied to biological activity QSPR (quantitative structure-property relationship) a technique for computing chemical properties... 

To calculate the properties of a molecule, you need to generate a well-defined structure. A calculation often requires a structure that represents a minimum on a potential energy surface. HyperChem contains several geometry optimizers to do this. You can then calculate single point properties of a molecule or use the optimized structure as a starting point for subsequent calculations, such as molecular dynamics simulations. 

However, before proceeding with the description of simulation data, we would like to comment the theoretical background. Similarly to the previous example, in order to obtain the pair correlation function of matrix spheres we solve the common Ornstein-Zernike equation complemented by the PY closure. Next, we would like to consider the adsorption of a hard sphere fluid in a microporous environment provided by a disordered matrix of permeable species. The fluid to be adsorbed is considered at density pj = pj-Of. The equilibrium between an adsorbed fluid and its bulk counterpart (i.e., in the absence of the matrix) occurs at constant chemical potential. However, in the theoretical procedure we need to choose the value for the fluid density first, and calculate the chemical potential afterwards. The ROZ equations, (22) and (23), are applied to decribe the fluid-matrix and fluid-fluid correlations. These correlations are considered by using the PY closure, such that the ROZ equations take the Madden-Glandt form as in the previous example. The structural properties in terms of the pair correlation functions (the fluid-matrix function is of special interest for models with permeabihty) cannot represent the only issue to investigate. Moreover, to perform comparisons of the structure under different conditions we need to calculate the adsorption isotherms pf jSpf). The chemical potential of a... 

In this situation computer simulation is useful, since the conditions of the simulation can be chosen such that full equihbrium is established, and one can test the theoretical concepts more stringently than by experiment. Also, it is possible to deal with ideal and perfectly flat surfaces, very suitable for testing the general mechanisms alluded to above, and to disregard in a first step all the complications that real substrate surfaces have (corrugation on the atomistic scale, roughness on the mesoscopic scale, surface steps, adsorbed impurities, etc.). Of course, it may be desirable to add such complications at a later stage, but this will not be considered here. In fact, computer simulations, i.e., molecular dynamics (MD) and Monte Carlo (MC) calculations, have been extensively used to study both static and dynamic properties [11] in particular, structural properties at interfaces have been considered in detail [12]. 

We observe that for the Fe-Co system a sim le spin polarized canonical model is able to reproduce qualitatively the results obtained by LMTO-CPA calculations. Despite the simplicity of this model the structural properties of the Fe-Co alloy are explained from simple band-filling arguments. 

In a previous work we showed that we could reproduce qualitativlely the LMTO-CPA results for the Fe-Co system within a simple spin polarized canonical band model. The structural properties of the Fe-Co alloy can thus be explained from the filling of the d-band. In that work we presented the results in canonical units and we could of course not do any quantitative comparisons. To proceed that work we have here done calculations based on the virtual crystal approximation (VGA). In this approximation each atom in the alloy has the same surrounding neighbours, it is thus not possible to distinguish between random and ordered alloys, but one may analyze the energy difference between different crystal structures. 

To summarize we have reproduced the intricate structural properties of the Fe-Co, Fe-Ni and the Fe-Cu alloys by means of LMTO-ASA-CPA theory. We conclude that the phase diagram of especially the Fe-Ni alloys is heavily influenced by short range order effects. The general trend of a bcc-fcc phase transition at lower Fe concentrations is in accordance with simple band Ailing effects from canonical band theory. Due to this the structural stability of the Fe-Co alloys may be understood from VGA and canonical band calculations, since the common band model is appropriate below the Fermi energy for this system. However, for the Fe-Ni and the Fe-Cu system this simple picture breaks down. 

From the optical absorption of two different hexaphenyl films, one with its chains predominantly standing upright on the substrate, the other with the chains randomly distributed in all orientations, similar structure property relations can be concluded [139]. By comparing the calculated absorption coefficient [139J perpendicular to the chains with the observed optical absorption spectra of both films we see that the optical absorption, plotted in Figure 9-9, in the visible and... 

---