UKCP - United Kingdom Car-Parrinello Consortium

 
   Site Navigation
  UKCP Home Page
  The CASTEP Code
  Technical Support
  Organisation
  UKCP and Accelrys
  Useful Links
  Conference Calendar
  Job Opportunities
  Further Information
  Publications
   Related Sites
Accelrys

 

 

First-principles simulation

Calculating the properties of materials fron scratch

First-principles calculations involve no adjustable parameters, are based on an explicit quantum treatment of the electrons in a model system, which means solving Schroedinger's equation to find the electronic ground state There are still approximations in the solution, but these are of an extremely general kind. The calculations yield the total energy of assemblies of atoms, and forces on the atoms. This knowledge enables us to calculate a vast range of properties with near chemical accuracy.

Density Functional Theory (DFT) is a powerful, general theory relating the total energy of a system of interacting electrons in an external potential to the electron density. Almost all our calculations employ DFT to describe the valence electrons. The expansion in the scope of such calculations has been breathtaking in the last fifteen years or so, and they are now routinely used to study fundamental aspects of the physics and chemistry of condensed matter in every conceivable guise, ranging from catalysts to protiens, minerals to liquid metals. The great importance of DFT and its application in simulation was recognised by the award of the 1998 Nobel Prize in chemistry to Walter Kohn and John Pople.

Car-Parrinello Methods

The bedrock of UKCP's science is the DFT plane-wave pseudopotential method, pioneered in the 80's by Roberto Car and Michele Parrinello. This approach has transformed the first-principles landscape, making possible calculations far beyond previous approaches. Throughout this time the UK has played a leading role in the development and application of the methods, an effort founded on Mike Payne's work with the CASTEP code.

The methods rely on a plane-wave basis, pseudopotentials and the use of density-functional theory to describe the valence electrons in a model system. Other ingredients include fast-Fourier transforms and minimization of the total energy rather than matrix diagonalization.

The main approximations

In Car-Parrinello calculations there are two main approximations. The first is the way in which the exchange-and-correlation energy of the electrons is handled. DFT provides a very general theory involving functionals of the electron density, but does not say what these functionals are, only that they must exist. In practise it has been found that the Local-Density Approximation (LDA), sometimes extended to take account of density gradients (the Generalised Gradient Approximation or GGA) are remarkably accurate, despite their simplicity. The LDA has been in widespread use for well over 20 years, and is known from vast numbers of calculations to be very reliable for large classes of systems.

The second approximation is the pseudopotential approximation. This is basically a way of avoiding the need to include all the electrons in the calculations. Only the valence electrons are treated explicitly. The core electrons are assumed to be in exactly the same state as in the free atom. The effective interaction between the valence electrons and the ion cores is described by a pseudopotential for each kind of atom. This pseudopotential is generated by first-principles calculations on free atoms. In the early days of the pseudopotential method, only a limited range of chemical elements could easily be handled. This ceased to be true a long time ago, and nowadays the method is routinely used for transition metals and first-row elements, which used to be regarded as `difficult'.

What you can calculate

The total energy and forces supplied by first-principles calculations can be used to do many things. One of the simplest and most important is the determination of equilibrium structures. The LDA-pseudopotential method was already in routine use for studying crystal structures ten years ago. Now, it is widely used for more complex problems, such as the equilibrium structure of point defects, grain boundaries, surfaces and molecules at surfaces. Closely related to the treatment of equilibrium structures is the calculation of vibrational properties, which is also widely practiced.

We also use first-principles calculations for dynamical simulations, where the forces on the atoms are used to generate the time evolution of the atomic positions. This is the idea that has long been used in classical molecular dynamics simulation, but now the forces are calculated from first principles, instead of being represented by empirical models. This means that dynamical processes involving the making and breaking of chemical bonds can be studied. Dynamical simulations allow the study of solids and liquids in thermal equilibrium, through the calculation of quantities such as thermodynamic functions, radial distribution functions, diffusion coefficients and dynamical structure factors. Even the first-principles calculation of phase diagrams is becoming possible. Another use of dynamical simulation is for the study of non-equilibrium processes like the dissociative chemisorption of molecules at surfaces.

Technical Description and History

Beginnings

The Cambridge Serial Total Energy Package was developed in the mid-'80s by Mike Payne. During that time and subsequently, many people in numerous research groups contributed to the code. CASTEP calculates the total energy, forces and stresses in a 3D-periodic system. The review article Rev. Mod. Phys. 64 (1992) 1045 sets out in detail the theory and technical implementation of CASTEP. By the end of the '80s it was clear the recipe of Density-Functional Theory (DFT), plane-waves, pseudopotentials and Fast-Fourier Transforms (FFT's) was revolutionising first-principles calculations and endowing them with unprecedented power.

CETEP: Meeting the Grand Challenges

UKCP was formed in 1990 with the aim of exploiting parallel supercomputers for first-principles simulation. It was recognised that realistic and accurate simulations would inherently be on a grand scale, and that parallel machines were essential. Out of this came CETEP, a highly efficient parallel version of CASTEP in which the real- and reciprocal-space grids were decomposed and distributed among the nodes of the parallel machine. UKCP's ability to exploit parallel computers has consistently enabled us to tackle and solve the most difficult and demanding of problems - real Grand Challenges. This relies heavily on the UK's investment in supercomputing, from the Meiko Computing Surface in the early days through to the latest CSAR service.

CASTEP Now

Since 1993 Accelrys Ltd (formely MSI) have developed and consolidated CASTEP as a commercial product, while CETEP remained the workhorse of UKCP academics. However, within the new UKCP-Accelys agreement these strengths are combined. CASTEP is now a state-of-the-art plane-wave code, incorporating all the latest and best technology. Once parallel development is complete it will derive the benefits of several years' of CETEP parallel experience.

CASTEP 4.2 Key Features

  • Ultrasoft or norm-conserving pseudopotentials with non-linear core corrections
  • Complete pseudopotential library
  • Comprehensive range of minimisation methods: Density Mixing, RM-DIIS, Conjugate Gradients band-by-band and all-bands
  • Full structural relaxation and molecular dynamics capabilities
  • Local Density and Generalised Gradient approximations, spin-polarisation
  • Supported by 100's of successful and published applications
  • Compatable with Accelrys Cerius2 and Materials Studio user interfaces giving automatic setup and full visualisation facilities

Distribution

The UKCP-Accelrys agreement secures the right to the CASTEP source code for all UK academics. UKCP believes this is a remarkable benefit, offering a high-quality code which is the product of many years' effort by academic researchers and software professionals. Even more important though is the opportunity for UK researchers to contribute to the continued development and enhancement of CASTEP. Welcome to CASTEP and UKCP!

Obtaining the code

Daresbury Laboratory distributes CASTEP to UK academics. To obtain the code you will be required to sign a license agreement and provide some details about yourself and your institution. After we receive these the software will be sent to you. Our address is:
Dr Adrian Wander
UKCP/CASTEP
CLRC Daresbury Laboratory
Keckwick Lane
Warrington WA4 4AD, UK.
E_mail: ukcp@dl.ac.uk

Academic release 4.2 is now available for distribution. Please write to the above address, clearly marking your letter or e_mail "UKCP/CASTEP", to obtain the requisite forms. Once we receive these completed and in order the software will be posted to you. Please allow reasonable time for delivery!

Please note: at present this agreement covers academic researchers in the UK only, and we regret that we cannot distribute CASTEP to researchers in other countries.

Electronic versions of the CASTEP paperwork are also available via e-mail.

The Accelrys Cerius 2 interface

Accelrys provide a full graphical user interface to CASTEP as part of the Cerius 2 package. However, this software is not covered by the UKCP-Accelrys agreement, and must be bought from Accelrys. UKCP has, though, agreed special rates with Accelrys for a Cerius 2 CASTEP package available to UK academics. Accelrys's contact details are as follows.
Accelrys Ltd.
334 Cambridge Science Park
Cambridge
CB4 0WN
Tel: +44 1223 228500
Fax: +44 1223 228501

 
 
   
 
 
back to top