Computational Methodologies

One should remember that the critical properties for effective high-energetic materials are large exothermicities and high dissociation barriers [2-4], This fact emphasizes that accurate calculations of the potential energy surface of polynitrogens 

Tetrahedral N4 is expected to dissociate into two N2 molecules, but this reaction is forbidden by orbital symmetry. Dunn and Morokuma [33] characterized a transition state for the exothermic dissociation of tetrahedral N4 into two N2 and estimated the activation barrier to be 63 kcal/mol at the CASSCF(12e,12o) level, which indicates that N4 is a metastable species with significant kinetic stability. The calculated potential energy surface of N4 suggests that the low-lying triplet state might cross with the singlet surface (Fig. 3), which could reduce the activation energy barrier to about 30 kcal/mol [29,31,32], 

Recent experimental studies supports the existence of tetrahedral N4 in the nitrogen plasma as well as in liquid and solid nitrogen [9,10], whereas neutralization-reionization mass spectrometry study of N4 suggests that observed species is the open-chain triplet structure [11]. 

Previous theoretical calculations [22,35,48-55] have shown that the potential energy surface (PES) of N6 isomers is highly dependent on the level of theory and basis sets. For instance, cyclic N6 with D6h symmetry, isoelectronic with benzene, was calculated to be a minimum at the HF level [49], but to become a higher-order saddle point at the MP2 level on the PES [35,49]. Higher-level calculations indicate that the lowest-energy form of cyclic N6 is not the benzene-like D6h structure but the twisted-boat geometry with D6h symmetry [22,35,50]. This D6h structure is probably not stable at room temperature because the dissociation frequency mode is only 73.6 cm-i at the CCSD(T) level [22], 

The global minimum of N6has an open-chain diazide structure with C2h symmetry 

So far we have indicated how once the scattering matrix of a reactive system is known, the scattering amplitudes and cross sections can be obtained. In this section we outline how that matrix can be calculated once the system s electronically adiabatic potential energy function V Rx,rx, x) is known. 

As in the collinear case, several approaches are possible. Around 1976, three coupled-channel methods had been used in cross-section calculations for 3-PD systems. One of them, developed by Elkowitz and Wyatt [50, 51], used natural collision coordinates (NCC) and local hindered asymmetric-top-vibrator basis sets [41]. Another, developed by Kuppermann and Schatz [106], used asymptotic free rotor and local vibrator basis sets, and different coordinates in different regions of configuration space, similar to those described 

These are simultaneous eigenfunctions of the operators Jop Jzop J defined at the beginning of an earlier section. They span the 4-MD subspace defined by the angular coordinates 0a, 7a, We now expand, in this basis set, the wave functions which are simultaneous eigenfunctions of 

Replacement of this expansion in the Schrodinger equation furnishes the following set of coupled-channel partial differential equations for the functions 

terms are centrifugal coupling terms that are independent 

---

Gasteiger J, C Rudolph and J Sadowski 1990. Automatic Generation of 3D Atomic Coordinates fo Organic Molecules. Tetrahedron Computer Methodology 3 537-547. 

A number of books and journal articles reviewing computational methods relevant to biophysical problems have been published in the last decade. Two of the most popular texts, however, were published more than ten years ago those of McCammon and Harvey in 1987 and Brooks, Karplus, and Pettitt in 1988. There has been significant progress in theoretical and computational methodologies since the publication of these books. Therefore, we feel that there is a need for an updated, comprehensive text including the most recent developments and applications in the field. 

More recently, two-state E/B models have been proposed by Chujo and Doi (9.10) for the analysis of polypropylene. Similar E/B models were proposed by Cheng(11) and Asakura, et al(12) for polybutylene. For copolymers, two-state B/B models have been proposed for ethylene-propylene copolymers,(11,13-15) and propylene-butylene copolymers.(11,13) Recently, Cheng(11) generalized these multi-state models and developed computer methodology for the general analysis of such systems. A number of polymer systems were treated. 

This chapter summarizes the computational methodologies used for conformational analysis. Specifically, Section 8.1 gives a theoretical outUne of the problem and presents details of various implementations of computer codes to perform conformational analysis. Section 8.2 describes calculations illustrative of the current accuracy in generating the conformation of a ligand when bound to proteins (the bioactive conformer) by comparisons to crystaUographically observed data. Finally, Section 8.3 concludes by presenting some practical... 

Certain computational methodologies such as some approaches to quantitative structure-activity relationship (QSAR) studies use 3D ligand structures [37, 38]. These methods generally assume that a bioactive conformation has been estab-Hshed for a set of molecules and that these conformers can be ahgned in a maimer that reflects the relative orientation they would adopt in a binding site. It is thus... 

In more recent years, additional progress and new computational methodologies in macromolecular quantum chemistry have placed further emphasis on studies in transferability. Motivated by studies on molecular similarity [69-115] and electron density representations of molecular shapes [116-130], the transferability, adjustability, and additivity of local density fragments have been analyzed within the framework of an Additive Fuzzy Density Fragmentation (AFDF) approach [114, 131, 132], This AFDF approach, motivated by the early charge assignment approach of Mulliken [1, 2], is the basis of the first technique for the computation of ab initio quality electron densities of macromolecules such as proteins [133-141],... 

Kearsley, S.K. and Smith, G.M. (1990). An alternative method for the alignment of molecular structures Maximizing electrostatic and steric overlap. Tetrahedron Computer Methodology 3 615-633. 

In this chapter the basic theory of the structurally coupled QM/MM is summarized. This is followed by some technical points important in the practical use of the method. In particular, details about the treatment of the QM/MM boundary are discussed. The thermodynamically coupled quantum mechanical/ free energy (QM/FE) method is then introduced. Some representative applications of QM/MM methods are then described. The examples are selected to provide a representative picture of the potential applications of QM/MM methods on studies of reaction mechanisms. Here there is special emphasis on recent advances in the computational methodologies and in the future developments needed to improve the applicability of the methods. 

Rapid advances are taking place in the application of DFT to describe complex chemical reactions. Researchers in different fields working in the domain of quantum chemistry tend to have different perspectives and to use different computational approaches. DFT owes its popularity to recent developments in predictive powers for physical and chemical properties and its ability to accurately treat large systems. Both theoretical content and computational methodology are developing at a pace, which offers scientists working in diverse fields of quantum chemistry, cluster science, and solid state physics. 

We start by discussing various means of incorporating surface roughness into the model systems in order to perform more realistic simulations. The means of subjecting the system to shear and load are discussed below. Thermostats are then discussed. Finally, we consider cases in which one can neglect the walls and treat the system as a bulk fluid. We finish with a discussion of different computational methodologies that are used in tribological simulations. 

An important characteristic of ab initio computational methodology is the ability to approach the exact description - that is, the focal point [11] - of the molecular electronic structure in a systematic manner. In the standard approach, approximate wavefunctions are constructed as linear combinations of antisymmetrized products (determinants) of one-electron functions, the molecular orbitals (MOs). The quality of the description then depends on the basis of atomic orbitals (AOs) in terms of which the MOs are expanded (the one-electron space), and on how linear combinations of determinants of these MOs are formed (the n-electron space). Within the one- and n-electron spaces, hierarchies exist of increasing flexibility and accuracy. To understand the requirements for accurate calculations of thermochemical data, we shall in this section consider the one- and n-electron hierarchies in some detail [12]. 

G. W. Robinson, S. Singh, and M. W. Evans, Water in Biology, Chemistry and Physics. Experimental Overviews and Computational Methodologies, World Scientific, Singapore, 1996. 

Academia has been rich in producing theoretical computational methodology that underpins molecular modeling. The following software arose from universities or private and publicly funded institutes AMBER [10], INSIGHT [11,12], CHARMM [13], SYBYL [14], GRID [15], DOCK [16] and HINT (Hydropathic INTeractions) [39]. All except AMBER were commercialized. 

Tristan Gerard Alfred Youngs obtained his PhD in computational chemistry from the University of Reading in 2004, and moved to Queen s University Belfast to take up a research fellowship working on ionic liquids shortly afterward. His interests focus on the use of computational methodology to examine ionic liquids at the atomic and molecular levels, and works closely with experimentalists and theorists alike. 

A number of computational methodologies have been used in this experiment. Reagent selections for library designs were exported... 

In order to provide accurate static values for this purpose, we have recently undertaken a series of calculations of the polarizability and hyperpolarizability of the rare gases. An extensive basis set investigation was performed for Ne, and we shall consider these results in detail. We shall also discuss aspects of the correlation treatment and computational methodology. We begin by considering methods for the calculation of the polarizability and hyperpolarizability. 

---