Fundamental organometallic reactions: Applications on the CYBER 205

Two of the most challenging problems of Organometallic chemistry (loosely defined) are pollution control with the large space velocities needed and nitrogen fixation, a process so capably done by nature and so relatively poorly done by man (industry). For a computational chemist these problems are on the fringe of what is possible with conventional computers (large models needed and accurate energetics required). A summary of the algorithmic modification needed to address these problems on a vector processor such as the CYBER 205 and a sketch of findings to date on deNOx catalysis and nitrogen fixation are presented

Structural and electronic properties of the metal-metal intramolecular junctions of single-walled carbon nanotubes

Several intramolecular junctions (IMJs) connecting two metallic (11, 8) and (9, 6) carbon nanotubes along their common axis have been realized by using a layer-divided technique to the nanotubes and introducing the topological defects. Atomic structure of each IMJ configuration is optimized with a combination of density-functional theory (DFT) and the universal force field (UFF) method, based upon which a four-orbital tight-binding calculation is made on its electronic properties. Different topological defect structures and their distributions on the IMJ interfaces have been found, showing decisive effects on the localized density of states, while the sigma-pi coupling effect is negligible near Fermi energy (EF). Finally, a new IMJ model has been proposed, which probably reflects a real atomic structure of the M-M IMJ observed in the experiment [Science 291, 97 (2001)].Comment: 11 pages and 3 figure

A mixed ultrasoft/normconserved pseudopotential scheme

A variant of the Vanderbilt ultrasoft pseudopotential scheme, where the normconservation is released for only one or a few angular channels, is presented. Within this scheme some difficulties of the truly ultrasoft pseudopotentials are overcome without sacrificing the pseudopotential softness. i) Ghost states are easily avoided without including semicore shells. ii) The ultrasoft pseudo-charge-augmentation functions can be made more soft. iii) The number of nonlocal operators is reduced. The scheme will be most useful for transition metals, and the feasibility and accuracy of the scheme is demonstrated for the 4d transition metal rhodium.Comment: 4 pages, 2 figure

Ab Initio Study of Screw Dislocations in Mo and Ta: A new picture of plasticity in bcc transition metals

We report the first ab initio density-functional study of screw dislocations cores in the bcc transition metals Mo and Ta. Our results suggest a new picture of bcc plasticity with symmetric and compact dislocation cores, contrary to the presently accepted picture based on continuum and interatomic potentials. Core energy scales in this new picture are in much better agreement with the Peierls energy barriers to dislocation motion suggested by experiments.Comment: 3 figures, 3 table

Structure and Vibrations of the Vicinal Copper (211) Surface

We report a first principles theoretical study of the surface relaxation and lattice dynamics of the Cu(211) surface using the plane wave pseudopotential method. We find large atomic relaxations for the first several atomic layers near the step edges on this surface, and a substantial step-induced renormalization of the surface harmonic force constants. We use the results to study the harmonic fluctuations around the equilibrium structure and find three new step-derived features in the zone center vibrational spectrum. Comparison of these results with previous theoretical work and weith experimental studies using inelastic He scattering are reported.Comment: 6 Pages RevTex, 7 Figures in Postscrip

Angular Forces Around Transition Metals in Biomolecules

Quantum-mechanical analysis based on an exact sum rule is used to extract an semiclassical angle-dependent energy function for transition metal ions in biomolecules. The angular dependence is simple but different from existing classical potentials. Comparison of predicted energies with a computer-generated database shows that the semiclassical energy function is remarkably accurate, and that its angular dependence is optimal.Comment: Tex file plus 4 postscript figure

Surface Enhanced Second Harmonic Generation from Macrocycle, Catenane, and Rotaxane Thin Films: Experiments and Theory

Surface enhanced second harmonic generation (SE SHG) experiments on molecular structures, macrocycles, catenanes, and rotaxanes, deposited as monolayers and multilayers by vacuum sublimation on silver, are reported. The measurements show that the molecules form ordered thin films, where the highest degree of order is observed in the case of macrocycle monolayers and the lowest in the case of rotaxane multilayers. The second harmonic generation activity is interpreted in terms of electric field induced second harmonic (EFISH) generation where the electric field is created by the substrate silver atoms. The measured second order nonlinear optical susceptibility for a rotaxane thin film is compared with that obtained by considering only EFISH contribution to SHG intensity. The electric field on the surface of a silver layer is calculated by using the Delphi4 program for structures obtained with TINKER molecular mechanics/dynamics simulations. An excellent agreement is observed between the calculated and the measured SHG susceptibilities.

Supercell technique for total-energy calculations of finite charged and polar systems

We study the behavior of total-energy supercell calculations for dipolar molecules and charged clusters. Using a cutoff Coulomb interaction within the framework of a plane-wave basis set formalism, with all other aspects of the method (pseudopotentials, basis set, exchange-correlation functional) unchanged, we are able to assess directly the interaction effects present in the supercell technique. We find that the supercell method gives structures and energies in almost total agreement with the results of calculations for finite systems, even for molecules with large dipole moments. We also show that the performance of finite-grid calculations can be improved by allowing a degree of aliasing in the Hartree energy, and by using a reciprocal space definition of the cutoff Coulomb interaction

A Geometric Formulation of Quantum Stress Fields

We present a derivation of the stress field for an interacting quantum system within the framework of local density functional theory. The formulation is geometric in nature and exploits the relationship between the strain tensor field and Riemannian metric tensor field. Within this formulation, we demonstrate that the stress field is unique up to a single ambiguous parameter. The ambiguity is due to the non-unique dependence of the kinetic energy on the metric tensor. To illustrate this formalism, we compute the pressure field for two phases of solid molecular hydrogen. Furthermore, we demonstrate that qualitative results obtained by interpreting the hydrogen pressure field are not influenced by the presence of the kinetic ambiguity.Comment: 22 pages, 2 figures. Submitted to Physical Review B. This paper supersedes cond-mat/000627