Theoretical Results

In contrast to the few reports of attempts to synthesize polynitrogens, there are himdreds of computational papers on these species. Fortunately, a fraction of these will suffice to give a fair perspective on om theoretical understanding of these (almost all experimentally imknown) compoimds. We will first report on some sitrveys of a variety of polynitrogens, then look in some detail at a few of those of particttlar irrterest, considering cyclic and then acyclic molecules. 

Hexaazabenzene (12), D6hMP2/6-31G geom Hexaazabenzene (18), D2 MP2/6-31G geom Acyclic Ng, C2 (17a) MP2/6-31 G geom Acyclic Ng, C2h (17b) B3LYP/6-31 G geom 

Most numerical approaches have considered the adsorption of protons and hydroxyl groups on surfaces in the limit of zero coverage. By modelling the surfaces by small clusters of a few atoms, on the basis of static energetic considerations, they have studied whether a dissociative adsorption occurs, and have calculated the charge transfers between H or OH and the surface and OH stretching frequencies. They have focused on the (110) faces of Ti02 and the (001) face of SrTiOs (Tsukada et al, 1983), several 

Due to the development of advanced numerical methods in the last decades, quantum approaches are now able to accurately describe the chemical bonds formed between two reactants. Nevertheless, when a surface is involved, the actual systems met in practice, for example a dense polymeric layer adsorbed on a rough surface, cannot yet be simulated, because this would require too large a memory size or too long a computation time. Quantum calculations, thus, cannot compete with empirical models in the prediction of adhesion strengths. However, they may allow one to check their validity in model cases, for example small molecules adsorbed on a substrate, or large molecules adsorbed on a cluster of a few atoms which simulates the substrate. This has been done in a number of cases but, to the author s knowledge, mostly for adsorption processes on metallic surfaces. Numerical results for the adsorption of molecules on oxide surfaces may be found in the literature (Henrich and Cox, 1994), but there exists no systematic discussion in the framework of acid-base interactions. 

When a proton adsorbs on a surface oxygen, an electron transfer h+ takes place from the surface towards the proton. The sign of the transfer is reversed between an OH group and a surface cation (Awoh- 0). In the latter case, there is also a modification of the O-H bond, the terminal hydrogen bearing less electrons when the transfer to the surface is larger. By its very definition, the importance of the adsorbate-substrate charge transfer is a measure of the Lewis acidity or basicity. 

As far as electron transfers are concerned, i.e. in the framework of Lewis acidity or basicity, it may, thus, be concluded that the oxides from BaO 

On the three MgO faces, the charge transfer per bond does not depend much upon the surface orientation in the limit of zero coverage. When the density of adsorbates becomes large, the conclusion is modified. It is found that surfaces become more and more basic as the coordination number of the surface atoms decreases. 

---

We do not wish to go into the details of Figure 10. As an illustration of the reliability of the present results we compare, however, in Figure 11, the structure of the measured HCCS spectrum published by Tang and Saito (Fig. 3 of [139]) with the results of the theoretical study. Taking into account the very complex and unusual structure of this kind of spectra, we find the agreement between our ab initio theoretical results and those following from the interpretation of experimental spectra more than satisfactory. While strongly... 

We now compare the results calculated for the fundamental frequency of the symmetric stretching mode with the only available experimental datum [78] of 326 cm . The theoretical result is seen to exceed experiment by only 8.3%. It should be recalled that the Li3 and Li3 tiimers have for lowest J the values 0 and respectively. Thus, the istopic species Li3 cannot contribute to the nuclear spin weight in Eq. (64), since the calculations for half-integer J should employ different nuclear spin weights. Note that atomic masses have been used... 

Finally, in Sec. 5, the theoretical results are illustrated by applying two adaptive schemes to the collinear photo dissociation of ArHCl. 

We have derived time-reversible, symplectic, and second-order multiple-time-stepping methods for the finite-dimensional QCMD model. Theoretical results for general symplectic methods imply that the methods conserve energy over exponentially long periods of time up to small fluctuations. Furthermore, in the limit m —> 0, the adiabatic invariants corresponding to the underlying Born-Oppenheimer approximation will be preserved as well. Finally, the phase shift observed for symmetric methods with a single update of the classical momenta p per macro-time-step At should be avoided by... 

Table 1 describes the timing results (in seconds) for a system of 4000 atoms on 4, 8 and 16 nodes. The average CPU seconds for 10 time steps per processor is calculated. In the case of the force-stripped row and force-row interleaving algorithms the CPU time is reduced by half each time the number of processors is doubled. This indicates a perfect speedup and efficiency as described in Table 2. Tables 3, refibm table3 and 5 describe the timing results, speedups and efficiencies for larger systems. In particular. Table 4 shows the scaling in the CPU time with increase in the system size. These results are very close to predicted theoretical results.

The system of atomic units was developed to simplify mathematical equations by setting many fundamental constants equal to 1. This is a means for theorists to save on pencil lead and thus possible errors. It also reduces the amount of computer time necessary to perform chemical computations, which can be considerable. The third advantage is that any changes in the measured values of physical constants do not affect the theoretical results. Some theorists work entirely in atomic units, but many researchers convert the theoretical results into more familiar unit systems. Table 2.1 gives some conversion factors for atomic units. 

Because mesoscale methods are so new, it is very important to validate the results as much as possible. One of the best forms of validation is to compare the computational results to experimental results. Often, experimental results are not available for the system of interest, so an initial validation calculation is done for a similar system for which experimental results are available. Results may also be compared to any other applicable theoretical results. The researcher can verify that a sulficiently long simulation was run by seeing that the same end results are obtained after starting from several different initial configurations. 

An analytical model of the process has been developed to expedite process improvements and to aid in scaling the reactor to larger capacities. The theoretical results compare favorably with the experimental data, thereby lending vahdity to the appHcation of the model to predicting directions for process improvement. The model can predict temperature and compositional changes within the reactor as functions of time, power, coal feed, gas flows, and reaction kinetics. It therefore can be used to project optimum residence time, reactor si2e, power level, gas and soHd flow rates, and the nature, composition, and position of the reactor quench stream. 

The dominant crystal size, is most often used as a representation of the product size, because it represents the size about which most of the mass in the distribution is clustered. If the mass density function defined in equation 33 is plotted for a set of hypothetical data as shown in Figure 10, it would typically be observed to have a maximum at the dominant crystal size. In other words, the dominant crystal size is that characteristic crystal dimension at which drajdL = 0. Also shown in Figure 10 is the theoretical result obtained when the mass density is determined for a perfectiy mixed, continuous crystallizer within which invariant crystal growth occurs. That is, mass density is found for such systems to foUow a relationship of the form m = aL exp —bL where a and b are system-dependent parameters. 

Having considered how solvents can affect the reactivities of molecules in solution, let us consider some of the special features that arise in the gas phase, where solvation effects are totally eliminated. Although the majority of organic preparative reactions and mechanistic studies have been conducted in solution, some important reactions are carried out in the gas phase. Also, because most theoretical calculations do not treat solvent effects, experimental data from the gas phase are the most appropriate basis for comparison with theoretical results. Frequently, quite different trends in substituent effects are seen when systems in the gas phase are compared to similar systems in solution. 

The calculated energy differences give a good correlation with The p parameter (p = —17) is larger than that observed experimentally for proton exchange (p — 8). A physical interpretation of this is that the theoretical results pertain to the gas phase, where... 

Much of the experimental observations on carbon nanotubes thus far have been made on multi-wall tu-bules[15-19]. This has inspired a number of theoretical calculations to extend the theoretical results initially obtained for single-wall nanotubes to observations in multilayer tubules. These calculations for multi-wall tubules have been informative for the interpretation of experiments, and influential for suggesting new re-... 

When using these theoretical results to analyze onion-like particles, we must take into account that calculations are performed for single graphitic shells, which are subsequently arranged concentrically and, then, conclusions are obtained about the minimal energy configuration. This fact arises from the limited number of atoms that may be included in a calculation due to present computational capabilities (the smallest onion-like particles are formed by in a C240 and this system represents 300 atoms). 

The experimental and theoretical results for E., are shown in Figure 3-41 for a resin content by weight ranging from 10% to 100%. Because E. is not a function of C, only k was varied — two values were chosen k = 1 and k =. 9. Some experimental results in Figure 3-41 lie above the curve for k=1 (i.e., above the upper bound ) some results lie below k =. 9. However, most results lie between k =. 9 and k = 1 with k =. 9 being a conservative estimate of the behavior. The actual specimens were handmade, so the resin content might not be precise, and fiber misalignment is not unexpected. Thus, the results above the upper bound are not unusual nor is the basic fact of variation in E. ... 

The theoretical and measured results for E, are shown in Figure 3-41 as a function of resin content by weight. Theoretical results from Equation (3.64) are shown for C = 0,. 2,. 4, and 1, and the data are bounded by the curves for C = 0 and C =. 4. The theoretical curve labeled glass-resin connected in series is a lower, lower bound than the C = 0 curve and is an overly conservative estimate of the stiffness. 

The experimental results for of a glass-epoxy composite material are shown along with the theoretical prediction from Equation (3.66) as a function of resin content by weight in Figure 3-44. Theoretical results are shown for contiguity factors of C = 0,. 2,. 4, and 1. Apparently, C = 0 is the upper limit of the data whereas C =. 4 is the lower limit. Thus, the concept of contiguity factor is further reinforced. 

The measured stiffnesses for two- and three-layered special cross-ply laminates are shown with symbols in Figure 4-28, and the theoretical results are shown with solid lines. In all cases, the load was kept so low that no strain exceeded SOOp. Thus, the behavior was linear and elastic. The agreement between theory and experiment is quite good. Both the qualitative and the quantitative aspects of the theory are verified. Thus, the capability to predict cross-ply laminate stiffnesses exists and is quite accurate. 

Michael F. Card and Robert M. Jones, Experimental and Theoretical Results for Buckling of Eccentrically Stiffened Cylinders, NASA TN D-3639, October 1966. 

Finally, let us discuss the adsorption isotherms. The chemical potential is more difficult to evaluate adequately from integral equations than the structural properties. It appears, however, that the ROZ-PY theory reflects trends observed in simulation perfectly well. The results for the adsorption isotherms for a hard sphere fluid in permeable multiple membranes, following from the ROZ-PY theory and simulations for a matrix at p = 0.6, are shown in Fig. 4. The agreement between the theoretical results and compu-... 

Theoretical results of similar quality have been obtained for thermodynamics and the structure of adsorbed fluid in matrices with m = M = 8, see Figs. 8 and 9, respectively. However, at a high matrix density = 0.273) we observe that the fluid structure, in spite of qualitatively similar behavior to simulations, is described inaccurately (Fig. 10(a)). On the other hand, the fluid-matrix correlations from the theory agree better with simulations in the case m = M = S (Fig. 10(b)). Very similar conclusions have been obtained in the case of matrices made of 16 hard sphere beads. As an example, we present the distribution functions from the theory and simulation in Fig. 11. It is worth mentioning that the fluid density obtained via GCMC simulations has been used as an input for all theoretical calculations. 

Equation (5-40) gives the transition state theoretical result for the rate constant. 

Ethyl acetoacetate reacts with ehaleone in the presence of boron fluoride etherate affording 3-carbethoxy-2-methyl-4,6-diphenylpyryl-ium this can be hydrolyzed and decarboxylated to 2-methyl-4,6-diphenylpyrylium which should theoretically result from acetone and ehaleone. 

---