We introduce a simple local atomic structure optimization algorithm which is significantly faster than standard implementations of the conjugate gradient method and often competitive with more sophisticated quasi-Newton schemes typically used in ab initio calculations. It is based on conventional molecular dynamics with additional velocity modifications and adaptive time steps. The surprising efficiency and especially the robustness and versatility of the method is illustrated using a variety of test cases from nanoscience, solid state physics, materials research, and biochemistry.

Abstract. We present a program called potfit which generates an effective atomic interaction potential by matching it to a set of reference data computed in first-principles calculations. It thus allows to perform large-scale atomistic simulations of materials with physically justified potentials. We describe the fundamental principles behind the program, emphasizing its flexibility in adapting to different systems and potential models, while also discussing its limitations. The program has been used successfully in creating effective potentials for a number of complex intermetallic alloys, notably quasicrystals. Modelling

Classical effective potentials are indispensable for any large-scale atomistic simulations, and the relevance of simulation results crucially depends on the quality of the potentials used. For complex alloys like quasicrystals, however, realistic effective potentials are practically inexistent. We report here on our efforts to develop effective potentials especially for quasicrystalline alloy systems. We use the so-called force matching method, in which the potential parameters are adapted so as to optimally reproduce the forces and energies in a set of suitably chosen reference configurations. These reference data are calculated with ab-initio methods. As a first application, EAM potentials for decagonal Al-Ni-Co, icosahedral Ca-Cd, and both icosahedral and decagonal Mg-Zn quasicrystals have been constructed. The influence of the potential range and degree of specialisation on the accuracy and other properties is discussed and compared.

Abs1rad. The two main techniques for the generation of quasipcriodic tilings. de Bruijn's grid method and the projection formalism. are generalised. A vel)' broad class or quasi periodic tilings is obtained in this way. The two generalised methods are Shown to be equivalent The standard calculation of Fourier spectra is extended to the whole general class of tilings. Various upplications are discussed.

Based on a parallel scalable library for Coulomb interactions in particle systems, a comparison between the fast multipole method (FMM), multigrid-based methods, fast Fourier transform (FFT)-based methods, and a Maxwell solver is provided for the case of three-dimensional periodic boundary conditions. These methods are directly compared with respect to complexity, scalability, performance, and accuracy. To ensure comparable conditions for all methods and to cover typical applications, we tested all methods on the same set of computers using identical benchmark systems. Our findings suggest that, depending on system size and desired accuracy, the FMM- and FFT-based methods are most efficient in performance and stability.