Solid-state systems periodicity

In our opinion, the use of and calculations for one-particle Green s functions are uniquely suitable for solid-state systems periodic in any number of dimensions. When faithfully implemented, it satisfies all criteria above. Green s functions offer analytically compact and physically rich tools for representing many properties for extended, periodic systems. They satisfy powerful and elegant relations for quantities such as density of states, lifetimes for excitations, dielectric functions, photo-emission and absorption spectra, total crystal energies, and many more. 

The theoretical chemistry community developed density functional theory for finite molecular systems which involve molecules and cluster models that describe the catalytic systems. They use the same constructs used in many ab initio wavefunction methods, i.e. Gaussian or Slater basis sets. The solid-state physics community, on the other hand, developed density functional theory to describe bulk solid-state systems and infinite surfaces by using a supercell approach along with periodic basis functions, i.e. plane waves . Nearly all of our discussion has focused on finite molecular systems. In the next section we will describe in more detail infinite periodic systems. 

The approximations described above for H and P work well for performing quantum-mechanical calculations on isolated (i.e. gas-phase) molecular systems. We can extend this methodology to tackle solid-state problems, but the calculations are more complex. A molecular crystal is composed of an infinite three-dimensional periodic array and the electronic energy levels are more complicated than those of isolated molecules. Both of these factors will need to be considered when modeling solid-state systems. 

I focus here on a condensed-phase simulation problem involving the interactions of macromolecular stmctures, small molecules, and solvent. Periodic solid-state systems have their own special set of difficulties and simplifications, which 1 do not discuss. The problems of protein structure prediction, protein-ligand binding, and enzymatic catalysis, which are discussed in the next section, fall into this category. 

This produces reliable results for the chemical properties of molecules and solids. DPT calculations are also fairly computationally inexpensive [75], making this the method of choice for accurate calculations on large molecular and solid-state systems. Finally, DPT programs for periodic systems are now widely available... 

L. Trnflandier, M. Paris, C. Payen and F. Boucher. First-principles calculations within periodic bonndary conditions of the NMR shielding tensor for a transition metal nucleus in a solid state system The example of V-51 in AIVO4. J. Phys. Chem. B 110, 2006, 21403. 

Isolated gas ph ase molecules are th e sim plest to treat com pii tation -ally. Much, if not most, ch emistry lakes place in the liq iiid or solid state, however. To treat these condensed phases, you must simulate continnons, constant density, macroscopic conditions. The usual approach is to invoke periodic boundary conditions. These simulate a large system (order of 10" inoleeti les) as a contiruiotis replication in all direction s of a sm nII box, On ly th e m olceti Ics in the single small box are simulated and the other boxes arc just copies of the single box. 

Within the periodic Hartree-Fock approach it is possible to incorporate many of the variants that we have discussed, such as LFHF or RHF. Density functional theory can also be used. I his makes it possible to compare the results obtained from these variants. Whilst density functional theory is more widely used for solid-state applications, there are certain types of problem that are currently more amenable to the Hartree-Fock method. Of particular ii. Icvance here are systems containing unpaired electrons, two recent examples being the clci tronic and magnetic properties of nickel oxide and alkaline earth oxides doped with alkali metal ions (Li in CaO) [Dovesi et al. 2000]. 

The quantity x is a dimensionless quantity which is conventionally restricted to a range of —-ir < x < tt, a central Brillouin zone. For the case yj = 0 (i.e., S a pure translation), x corresponds to a normalized quasimomentum for a system with one-dimensional translational periodicity (i.e., x s kh, where k is the traditional wavevector from Bloch s theorem in solid-state band-structure theory). In the previous analysis of helical symmetry, with H the lattice vector in the graphene sheet defining the helical symmetry generator, X in the graphene model corresponds similarly to the product x = k-H where k is the two-dimensional quasimomentum vector of graphene. 

Following a decay period of at least 15 hrs, the 1,524 MeV 7-photopeak of 42K is counted for 40 min with a Ge-U solid state detector and a 4096-channel analyzer. This more sophisticated counting system is necessary to obtain the required resolution, since with a Na iodide detector, the 1,524 MeV peak of 42 K overlaps with the 1.369 MeV 7-peak of the 24Na decay spectrum... 

On the basis of these redox potentials it seems likely that direct electron release to the benzenediazonium ion takes place only with iodide. This corresponds well with experience in organic synthesis iodo-de-diazoniations are possible without catalysts, light, or other special procedures (Sec. 10.6). For bromo- and chloro-de-di-azoniations, catalysis by cuprous salts (Sandmeyer reaction, Sec. 10.5) is necessary. For fluorination the Balz-Schiemann reaction of arenediazonium tetrafluoroborates in the solid state (thermolysis) or in special solvents must be chosen (see Sec. 10.4). With astatide (211At-), the heaviest of the halide ions, Meyer et al. (1979) found higher yields for astato-de-diazoniation than for iodo-de-diazoniation, a result consistent with the position of At in the Periodic System. It has to be emphasized, however, that in investigations based on measuring yields of final products (Ar-Hal), the possibility that part of the yield may be due to heterolytic dediazoniation is very difficult to quantify. 

Stick-slip motion is another issue that has been explored using SFA. It is found that the occurrence of stick-slip depends on the sliding velocity and the stiffness of the system, and the mechanism of the phenomenon can be interpreted in terms of periodic transition between liquid and solid states of the conhned lubricant [40],... 

The SCF method for molecules has been extended into the Crystal Orbital (CO) method for systems with ID- or 3D- translational periodicityiMi). The CO method is in fact the band theory method of solid state theory applied in the spirit of molecular orbital methods. It is used to obtain the band structure as a means to explain the conductivity in these materials, and we have done so in our study of polyacetylene. There are however some difficulties associated with the use of the CO method to describe impurities or defects in polymers. The periodicity assumed in the CO formalism implies that impurities have the same periodicity. Thus the unit cell on which the translational periodicity is applied must be chosen carefully in such a way that the repeating impurities do not interact. In general this requirement implies that the unit cell be very large, a feature which results in extremely demanding computations and thus hinders the use of the CO method for the study of impurities. 

The three-dimensional symmetry is broken at the surface, but if one describes the system by a slab of 3-5 layers of atoms separated by 3-5 layers of vacuum, the periodicity has been reestablished. Adsorbed species are placed in the unit cell, which can exist of 3x3 or 4x4 metal atoms. The entire construction is repeated in three dimensions. By this trick one can again use the computational methods of solid-state physics. The slab must be thick enough that the energies calculated converge and the vertical distance between the slabs must be large enough to prevent interaction. 

Rao et al. [140] described a study using a commercial pH sensitive radio-telemetry capsule (RTC) to evaluate small bowel and colonic transit time in athletes with gastrointestinal symptoms. The RTC (type 7006 Remote Control Systems, London, UK) consists of a glass electrode with an integral reference cap and battery. RF transmissions from the capsule are detected by a solid-state receiver worn on the belt of the patient. The recorder samples the pH from the capsule at 6 second intervals for a period of 24 hours. They used pH changes as an indication of the pH capsule s movement. A sharp rise in pH from around pH 2 to pH 6 indicates that the capsule has moved into the duodenum from the stomach. Then the pH progressively rises to a plateau around pH 8, which indicates that the capsule has moved into the terminal ileum. Another commercially available wireless pH sensor (Bravo) from Medtronic Inc. has been used to measure esophageal pH for a period of 48 hours and it will be discussed later. 

The same group reported in 1986 a sensitive and selective HPLC method employing CL detection utilizing immobilized enzymes for simultaneous determination of acetylcholine and choline [187], Both compounds were separated on a reversed-phase column, passed through an immobilized enzyme column (acetylcholine esterase and choline oxidase), and converted to hydrogen peroxide, which was subsequently detected by the PO-CL reaction. In this period, other advances in this area were carried out such as the combination of solid-state PO CL detection and postcolumn chemical reaction systems in LC [188] or the development of a new low-dispersion system for narrow-bore LC [189],... 

Another approach is that of including dynamics in the calculations. A dynamical formalism of DFT was first developed by Car and Parrinello [31], and has been employed in a wide range of areas, e.g. solvation problems, reactions on surfaces, solid-state interactions, and a variety of biochemical applications. In CP-MD one normally uses a plane wave basis to reduce the computational requirements and enable easy implementation of periodic boundary conditions. Nonetheless, CP-MD simulations are rather costly, and are normally not applied to systems larger than, say, 1-200 atoms, and over relatively short time frames. 

---