Atomistic calculations

Since and depend only on die valence charge densities, they can be detennined once the valence pseudo- wavefiinctions are known. Because the pseudo-wavefiinctions are nodeless, the resulting pseudopotential is well defined despite the last temi in equation Al.3.78. Once the pseudopotential has been constructed from the atom, it can be transferred to the condensed matter system of interest. For example, the ionic pseudopotential defined by equation Al.3.78 from an atomistic calculation can be transferred to condensed matter phases without any significant loss of accuracy. 

Initial atomistic calculations on nucleic acids were perfonned in the absence of an explicit solvent representation, as discussed earlier. To compensate for this omission, various... 

The most obvious feature of this mechanism to test by our theoretical calculations concerns the energetics of APB formation. Thus atomistic calculations of the type discussed in Section 2 for the shear planes, showed that the <011> APB had a considerably lower formation energy than boundaries with other orientations. The result is in line with the proposed mechanism, as interstitial capture at the <011>... 

The subscript j refers to the j member of the database, and thus Ej is the energy of that member as obtained from the relevant atomistic calculation and 4>Q ( cr )7... 

One of the central conclusions derived from the Peierls-Nabarro analysis is the role of nonlinear effects in the dislocation core. From an atomistic perspective, the far field atomic displacements could be derived just as well from linear elasticity as the full nonlinear function that results from direct atomistic calculation. By way of contrast, in the core region it is the nonlinear terms that give rise to some of the complex core rearrangements that we take up now. Our discussion will be built around two key examples cores in fee metals and the core reconstructions found in covalent semiconductors such as in Si. 

Mills M. J., High Resolution Transmission Electron Microscopy and Atomistic Calculations of Grain Boundaries in Metals and Intermetallies, Mat. Sci. Eng. A166, 35 (1993). 

Hence, a further advancement in studies of dislocations and interfaces in TiAl requires introduction of non-central forces into atomistic calculations. This character of bonding is, of course, included in ab-initio electronic structure calculations but such calculations are still not feasible on the scale needed when investigating dislocations and interfaces. Empirical attempts to include directionality of bonding have been made (Panova and Farkas 1995) but only employment of queintum mechemics based potentials can reveal the effects of non-central forces unambiguously. Such potentials have been emerging recently but have not yet been employed extensively in studies of lattice defects (Aoki and Pettifor, 1994 Memh et al., 1995 Pettlfor et al., 1995). Incorporation of... 

A particularly striking illustration that makes us question the importance of E as a measure of the energy of a GB is provided by observations of the X = 99 and X = 41 boundaries in spinel, where, by chance, the interface can facet parallel to several pairs of (different) low-index planes in both grains simultaneously. A similar situation exists for phase boundaries. It could be that the faceting onto low-index planes would lower the energy to below that of similar GBs having a lower X. There are no atomistic calculations for such GBs. 

Recently, Yamamoto and Hyodo have employed the DPD method for studying Nafion membranes [20]. The systems considered in this study were built using two distinct molecular species, denoted comb-shaped polymer ip) and water (w). The polymer was presented as a branched sequence of beads. It consisted of a main chain (backbone) of iV = 20 effective monomer units (-CF2CF2CF2CF2 ) linked with rig = 5 short side chains of = 2 units [-0CF2C(CF3)F0 and F2CF2S03H] the total number of interaction sites in the macromolecule was Np= N/, + n xn = 30. A water-like particle was modeled as the same size as the units of the Nafion fragment (<7 = 6.1 A) and represented four water molecules. The x parameters were found using an atomistic calculation. The DPD simulation was performed for water volume... 

In the last section, we discussed the use of QC calculations to elucidate reaction mechanisms. First-principle atomistic calculations offer valuable information on how reactions happen by providing detailed PES for various reaction pathways. Potential energy surfaces can also be obtained as a function of electrode potential (for example see Refs. [16, 18, 33, 38]). However, these calculations do not provide information on the complex reaction kinetics that occur on timescales and lengthscales of electrochemical experiments. Mesoscale lattice models can be used to address this issue. For example, in Refs. [25, 51, 52] kinetic Monte Carlo (KMC) simulations were used to simulate voltammetry transients in the timescale of seconds to model Pt(l 11) and Pt(lOO) surfaces containing up to 256x256 atoms. These models can be developed based on insights obtained from first-principle QC calculations and experiments. Theory and/or experiments can be used to parameterize these models. For example, rate theories [22, 24, 53, 54] can be applied on detailed potential energy surfaces from accurate QC calculations to calculate electrochemical rate constants. On the other hand, approximate rate constants for some reactions can be obtained from experiments (for example see Refs. [25, 26]). This chapter describes the later approach. 

When chromophores are covalently coupled in super/supramolecular objects (such as multichromo-phore-containing dendrimers), Monte Carlo-molecular dynamic calculations must be modified to take into account the restrictions on motion associated with covalent bond potentials. To accurately account for covalent bond potentials, atomistic Monte Carlo methods are required [68]. However because of the large number of atoms involved, fiilly atomistic calculations would be prohibitively time-consuming... 

Last, in order to overcome some shortcomings in the KMC approach—such as the requirement for defining all transitions prior to the start of the simulation—linking KMC with molecular dynamics or with inputs from atomistic calculations would provide more quantitative support for the energy barrier calculations and allow for less rigid constraints on the simulated structure [44—46]. 

To understand how metals are transported in the Earth s crust, we need to predict the speciation of metal complexes in aqueous solutions as a function of pressure, temperature and composition. To do this using atomistic calculations, we must define an adequately large system (as a function of composition) and obtain a thermodynamic average of all the possible states of the system (as a function of pressure and temperature). 

The outline of the paper is as follows. In Sect. 2 we describe the basic RISM and PRISM formalisms, and the fundamental approximations invoked that render the polymer problem tractable. The predicticms of PRISM theory for the structure of polymer melts are described in Sect. 3 for a variety of single chain models, including a comparison of atomistic calculations for polyethylene melt with diffraction experiments. The general problem of calculating thermodynamic properties, and particularly the equation-of-state, within the PRISM formalism is described in Sect. 4. A detailed application to polyethylene fluids is summarized and compared with experiment. The develojanent of a density functional theory to treat polymer crystallization is briefly discussed in Sect. 5, and numerical predictions for polyethylene and polytetrafluoroethylene are summarized. 

Day, G. M., Price, S. L., and Leslie, M. 2003. Atomistic calculations of phonon frequencies and thermodynamic quantities for crystals of rigid organic molecules. J. Phys. Chem. B 107 10919. 

---