CECAM/Psi-k Workshop,Lyon, France, 27-29 June 2007

Maximally Localized Wannier Functions:
Concepts, Applications, and Beyond

CECAM/Psi-k workshop, Lyon, France, 27-29 June 2007

Conference Dinner Menu


Nicola Marzari Massachusetts Institute of Technology, United States
Alfredo Pasquarello EPFL, Lausanne, Switzerland
Ivo Souza University of California, Berkeley, United States


Please organize your accomodation referring to the specific CECAM web page.

Scientific program

Here is the scientific program of the workshop

Confirmed Invited Participants

Ole K. Andersen Max Planck Institut, Stuttgart, Germany
Christian Brouder University Pierre et Marie Curie, Paris, France
Roberto Car Princeton University, United States
Stefano de Gironcoli Sissa, Trieste, Italy
Oswaldo Dieguez Rutgers University, United States
Andrea Ferretti University of Modena and Reggio Emilia, Italy
Feliciano Giustino University of California, Berkeley, United States
François Gygi University of California, Davis, United States
Wei Ku Brookhaven National Laboratory, United States
Frank Lechermann University of Hamburg, Germany
Young-Su Lee Massachusetts Institute of Technology, United States
Nicola Marzari Massachusetts Institute of Technology, United States
Marco Bongiorno Nardelli North Carolina State University, United States
Alfredo Pasquarello EPFL, Lausanne, Switzerland
Michel Posternak EPFL, Lausanne, Switzerland
Emil Prodan Princeton University, United States
Raffaele Resta University of Trieste, Italy
Ivo Souza University of California, Berkeley, United States
Massimiliano Stengel University of California, Santa Barbara, United States
Timo Thonhauser Rutgers University, United States
Kristian Thygesen Technical University of Denmark, Denmark
Mark E. Tuckerman New York University, United States
Paolo Umari INFM-Democritos, Trieste, Italy
David Vanderbilt Rutgers University, United States
Jonathan Yates Cambridge University, United Kingdom
Oleg Yazyev EPFL, Lausanne, Switzerland

Other Confirmed Participants

Hazem Abu-Farsakh MPI fuer Eisenforschung, Germany
Berlijn Tom Brookhaven National Laboratory, United States
Bernard Amadon CEA-DIF, Bruyères le Châtel, France
Benedetta Bonferroni Universita' di Modena e Reggio Emilia, Italy
Albert Bartok-Partay University of Cambridge, United Kingdom
Leonardo Bernasconi Vrije Universiteit Amsterdam, Netherlands
Eric Cancès CERMICS-ENPC, France
Silvia Casassa Università di Torino, Italy
Frank Freimuth Juelich Research Centre, Germany
Javad Seyed Hashemifar University of Duisburg, Germany
Deepa Kasinathan Max Planck Institute, Dresden, Germany
Amel Laref Oldenburg university, Germany
Oliver Marquardt Max Planck Institute for Iron Research, Germany
Yuriy Mokrousov University of Hamburg, Germany
Gianluca Panati University of Rome La Sapienza, Italy
Noejung Park Dankook University, Seoul, South Korea
Tonatiuh Rangel Université Caholique de Louvain, Belgium
Peter Reinhardt University Paris VI, France
Manu Sharma University of California, Davis, United States
Wenhui Xie Max Planck Institute, Stuttgart, Germany
Eva Zurek Max Planck Institute, Stuttgart, Germany


A tutorial on Wannier functions will be organized at CECAM during the day preceding our workshop, i.e. on Tuesday 26 June 2007. There will be some limited amount of support to cover the extra day and contributions to travel costs from a Marie Curie program (courtesy of the CECAM). This will only apply to PhD students or junior postdocs. For registering officially for this event, please contact Nicola Marzari.

Funding Options

We would like to point out that the NSF supports travel of junior US-based researchers to CECAM workshops through the Materials Computiation Center at the University of Illinois Urbana-Champaing. This is intended for young researchers (up to the assistant professor level, including graduate students and post-docs, but excluding senior researchers or scientists at US National Labs). We note that for participants eligible for NSF support CECAM's regulations prevent us from providing financial support for travel expenses unless it is made clear that an application to MCC was rejected. For information refer to this web page. For your convenience, there is also on-line facility for application.

Additional sources of funding for all US-based researchers are the European offices of the Naval Research and Aerospace Research. Psi-k strongly encourages US-based participants to apply to those institutions. Addresses are:
- U.S. Office of Naval Research Global
- European Office of Aerospace Research & Development

Scientific Background

The electronic structure of periodic systems is usually described in terms of the energy eigenfunctions, the Bloch functions, which are also eigenstates of the translation operator. However, these functions are spatially delocalized and offer a description which is not very intuitive. An alternative description of the electronic structure is provided by Maximally Localized Wannier Functions (MLWFs), which favor a representation of the electronic structure in terms of functions localized in real space [1,2]. These functions are obtained by unitary transformation from the energy eigenfunctions and therefore provide a fully equivalent description. Through their center, their extension and their shape, MLWFs provide a description of the electronic structure which is physically insightful. The electron orbitals can be visualized and located with respect to the positions of the atoms. Early investigations on complex crystals [2] and on disordered solids [3] and liquids [4] immediately illustrated the strength of this new concept. Essentially, a MLWF provides information on the local electronic structure in a straightforward and compact form.

Maximally localized Wannier functions are experiencing an increasing success as a monitoring tool of the electronic structure within both the electronic structure [5-9] and ab-initio molecular dynamics communities [10-12]. However, more recently, the particular properties of MLWFs have also attracted interest as a primary methodological tool within more complex theoretical developments. The most interesting aspect of MLWFs is that they allow a decomposition of the electronic structure in local components [1,2]. This property naturally allows one to locally monitor the response of the electronic structure when the investigated system is perturbed, either by the application of external (electric or magnetic) fields or by structural modification. In the study of the polarization, the usefulness of MLWFs is further enhanced by the fact that the displacement of their centers naturally connects with the modern definition of the polarization [1,2]. Hence MLWFs have become the method of choice to investigate the dielectric response in insulator-insulator [13-16] or insulator-metal interfaces [17]. Similarly, MLWFs have shown their usefulness in the study of the response of the electronic structure to a magnetic field [18-20]. These functions have been used in the calculation of NMR shifts [18], but have also turned out to be valuable in establishing the correct definition of the orbital magnetization in periodic systems [19,20]. Recent work on the interpretation of photoemission shifts at semiconductor surfaces and interfaces relies on MLWFs to highlight the role of strain fields [21].

Another exploited property of a set of MLWFs is that they constitute a complete basis set for the occupied states in which hierarchy can be established on the basis of distance to a particular point in space. This property is particularly useful for the introduction of approximations in which role of the local electronic structure is prioritized. Such basis functions have found a natural application in the development of methodologies for large-scale electronic calculations [22,23]. Recent developments also concern the research area of strongly correlated electrons, where MLWFs have recently been introduced as preferred functions for describing the local correlation [24-31].

Maximally localized Wannier functions have also found application in the construction of simplified Hamiltonians for multi-scale approaches. In this case, one uses the property that MLWFs can represent a particular energy window. This is straightforward when the relevant electronic bands are well separated from the other ones, but recent developments have extended such applications also to coupled bands introducing purposely designed disentanglement procedures [32-35]. In particular, these methodological developments are at the basis of several successful applications in the context of transport in nano-sized structures [36-38].

Although most of the early applications of MLWFs were to insulating systems, very recently the usefulness of partially-occupied MLWFs as fast and accurate interpolators for describing Fermi-surface properties of metals has been explored. For example, this has been used to reduce dramatically the computational load needed to evaluate properties such as the anomalous Hall effect of ferromagnets [39] and the electron-phonon coupling [40]. A further recent extension has shown that generalized Wannier functions can be constructed which show small spreads both in space and in energy [41].

The raise of interest in MLWFs has also prompted theoretical work in order to complete the characterization of these functions. Very recent work has been able to demonstrate the exponential nature of the localization [42], extending to three dimensions previous work for one dimension [43]. Furthermore, the real nature of MLWFs, empirically observed in many applications, could be formally demonstrated [42]. The wide-spread use of MLWFs in the ab initio simulation community has also led to the development of molecular dynamics integrators which preserve the localization properties during the course of the evolution [44-46].

A computer module for the generation of MLWFs is freely available [47] and facilitates their wide-spread use across the electronic-structure community. The success of MLWFs has brought several major software packages, such as DACAPO, CPMD and QUANTUM-ESPRESSO, to offer to their users a module for their generation. From the technical point of view, there have been important evolutions in the algorithms that are used to generate MLWFs [48-51].

Finally, it should be mentioned that the concept of MLWF has acquired generality giving rise to applications which extend beyond the description of the electronic structure. For instance, the concept has been used in the context of photonic crystals [52] and recent work shows how generalized Wannier functions can be used to study phonon Hamiltonians [41,40].

As clearly illustrated by this evolution, the use of MLWFs is rapidly expanding across diverse research areas. One of the primary aims of the present workshop is to bring together the major players in these different disciplines in order to share their experiences. These will include researchers in the electronic structure and ab initio simulation community, in physics and in chemistry. Researchers that are primarily interested in new developments and researchers that are more concerned with practical applications will be both represented. While the research areas are diverse, the optimal integration of MLWFs in broader schemes necessarily will face difficulties which are to a large extent common to all concerned research areas. The workshop represents a forum where researchers from various fields meet to expose their views, discuss their ideas, and confront their solutions. Furthermore the interaction between development and application will be beneficial to both parties leading to more focused developments on the one side, and to richer applications on the other.

The planned workshop is the very first one of this nature. The impact of the workshop is expected to be twofold. First, at this time it is not clear whether all the major researchers that have been identified in this proposal are fully aware of all the developments that are occurring in parallel in other research areas. The workshop will therefore contribute in removing potential communication barriers. This will be beneficial to several research fields. Second, the workshop will promote the interaction between development and application. Eventually this will lead to the optimal exploitation of the concept of MLWFs in various contexts.


[1] Marzari N, Vanderbilt D Maximally localized generalized Wannier functions for composite energy bands PHYSICAL REVIEW B 56, 12847 (1997).
[2] Marzari N, Souza I, Vanderbilt D An introduction to maximally-localized Wannier functions, Highlight of the Month, Psi-K Newsletter 57, 129 (2003), http://psi-k.dl.ac.uk/psi-k/newsletters.html.
[3] Silvestrelli PL, Marzari N, Vanderbilt D, Parrinello M Maximally-localized Wannier functions for disordered systems: Application to amorphous silicon SOLID STATE COMMUNICATIONS 107, 7 (1998).
[4] Silvestrelli PL, Parrinello M Structural, electronic, and bonding properties of liquid water from first principles JOURNAL OF CHEMICAL PHYSICS 111, 3572 (1999).
[5] Fernandez P, Dal Corso A, Baldereschi A Ab initio study of the dielectric properties of silicon and gallium arsenide using polarized Wannier functions PHYSICAL REVIEW B 58, R7480 (1998).
[6] Souza I, Martin RM, Marzari N, Zhao XY, Vanderbilt D Wannier-function description of the electronic polarization and infrared absorption of high-pressure hydrogen PHYSICAL REVIEW B 62, 15505 (2000).
[7] Posternak M, Baldereschi A, Massidda S, Marzari N Maximally localized Wannier functions in antiferromagnetic MnO within the FLAPW formalism PHYSICAL REVIEW B 65, 184422 (2002).
[8] Bert A, Llunell M, Dovesi R, Zicovich-Wilson CM Electronic structure characterization of six semiconductors through their localized Wannier functions PHYSICAL CHEMISTRY CHEMICAL PHYSICS 5, 5319 (2003).
[9] Cangiani G, Baldereschi A, Posternak M, Krakauer H Born charge differences of TiO2 polytypes: Multipole expansion of Wannier charge densities PHYSICAL REVIEW B 69, 121101 (2004).
[10] Kirchner B, Hutter J Solvent effects on electronic properties from Wannier functions in a dimethyl sulfoxide/water mixture JOURNAL OF CHEMICAL PHYSICS 121, 5133 (2004).
[11] Iftimie R, Tuckerman ME Decomposing total IR spectra of aqueous systems into solute and solvent contributions: A computational approach using maximally localized Wannier orbitals JOURNAL OF CHEMICAL PHYSICS 122, 214508 (2005).
[12] Sharma M, Resta R, Car R Intermolecular dynamical charge fluctuations in water: A signature of the H-bond network PHYSICAL REVIEW LETTERS 95, 187401 (2005).
[13] Giustino F, Umari P, Pasquarello A Dielectric discontinuity at interfaces in the atomic-scale limit: Permittivity of ultrathin oxide films on silicon PHYSICAL REVIEW LETTERS 91, 267601 (2003).
[14] Giustino F, Pasquarello A Theory of atomic-scale dielectric permittivity at insulator interfaces PHYSICAL REVIEW B 71, 144104 (2005).
[15] Stengel M, Spaldin NA Accurate polarization within a unified Wannier function formalism PHYSICAL REVIEW B 73, 075121 (2006).
[16] Wu X, Dieguez O, Rabe KM, Vanderbilt D Wannier-based definition of layer polarizations in perovskite superlattices PHYSICAL REVIEW LETTERS 97, 107602 (2006).
[17] Stengel M, Spaldin NA Ab-initio theory of metal-insulator interfaces in a finite electric field cond-mat/0511042.
[18] Sebastiani D, Parrinello M A new ab-initio approach for NMR chemical shifts in periodic systems JOURNAL OF PHYSICAL CHEMISTRY A 105, 1951 (2001).
[19] Thonhauser T, Ceresoli D, Vanderbilt D, Resta R Orbital magnetization in periodic insulators PHYSICAL REVIEW LETTERS 95, 137205 (2005).
[20] Ceresoli D, Thonhauser T, Vanderbilt D, Resta R Orbital magnetization in crystalline solids: Multi-band insulators, Chern insulators, and metals PHYSICAL REVIEW B 74, 024408 (2006).
[21] Yazyev OV, Pasquarello A Origin of fine structure in Si 2p photoelectron spectra at silicon surfaces and interfaces PHYSICAL REVIEW LETTERS 96, 157601 (2006).
[22] Skylaris CK, Mostofi AA, Haynes PD, Dieguez O, Payne MC Nonorthogonal generalized Wannier function pseudopotential plane-wave method PHYSICAL REVIEW B 66, 035119 (2002).
[23] Wahn M, Neugebauer J Generalized Wannier functions: An efficient way to construct ab-initio tight-binding parameters for group-III nitrides PHYSICA STATUS SOLIDI B 243, 1583 (2006).
[24] Schnell I, Czycholl G, Albers RC Hubbard-U calculations for Cu from first-principle Wannier functions PHYSICAL REVIEW B 65, 075103 (2002).
[25] Schnell I, Czycholl G, Albers RC Unscreened Hartree-Fock calculations for metallic Fe, Co, Ni, and Cu from ab initio Hamiltonians PHYSICAL REVIEW B 68, 245102 (2003).
[26] Nakamura K, Arita R, Yoshimoto Y, Tsuneyuki S Fabris S, de Gironcoli S, Baroni S, Vicario G, Balducci G Taming multiple valency with density functionals: A case study of defective ceria PHYSICAL REVIEW B 71, 041102 (2005).
[27] Zurek E, Jepsen O, Andersen OK Muffin-tin orbital Wannier-like functions for insulators and metals CHEMPHYSCHEM 6, 1934 (2005).
[28] Buth C, Birkenheuer U, Albrecht M, Fulde P Ab initio Green's function formalism for band structures PHYSICAL REVIEW B 72, 195107 (2005).
[29] Yin W-G, Volja D, Ku W Orbital Ordering in LaMnO3: Electron-Electron versus Electron-Lattice Interactions PHYSICAL REVIEW LETTERS 96, 116405 (2006).
[30] Nakamura K, Ryotaro A, Yoshimoto Y, Tsuneyuki S First-principles Calculation of Effective Onsite Coulomb Interaction of 3d Transition Metals: Constrained Local Density Functional Approach with Maximally Localized Wannier Function cond-mat/0510425.
[31] Lechermann F, Georges A, Poteryaev A, Biermann S, Posternak M, Yamasaki A, Andersen O K Dynamical mean-field theory using Wannier functions: a flexible route to electronic structure calculations of strongly correlated materials cond-mat/0605539.
[32] Souza I, Marzari N, Vanderbilt D Maximally localized Wannier functions for entangled energy bands PHYSICAL REVIEW B 65, 035109 (2002).
[33] Thygesen KS, Hansen LB, Jacobsen KW Partly occupied Wannier functions PHYSICAL REVIEW LETTERS 94, 026405 (2005).
[34] Thomas JW, Iftimie R, Tuckerman ME, Birkenheuer U, Izotov D Localization of Wannier functions for entangled energy bands PHYSICAL REVIEW B 71, 125116 (2005).
[35] Thygesen KS, Hansen LB, Jacobsen KW Partly occupied Wannier functions: Construction and applications PHYSICAL REVIEW B 72, 125119 (2005).
[36] Calzolari A, Marzari N, Souza I, Nardelli MB Ab initio transport properties of nanostructures from maximally localized Wannier functions PHYSICAL REVIEW B 69, 035108 (2004).
[37] Lee YS, Nardelli MB, Marzari N Band structure and quantum conductance of nanostructures from maximally localized wannier functions: The case of functionalized carbon nanotubes PHYSICAL REVIEW LETTERS 95, 076804 (2005).
[38] Thygesen KS, Jacobsen KW Molecular transport calculations with Wannier functions CHEMICAL PHYSICS 319, 111 (2005).
[39] Wang X, Yates JR, Souza I, Vanderbilt D Ab initio calculation of the anomalous Hall conductivity by Wannier interpolation cond-mat/0608257.
[40] Giustino F, Yates JR, Souza I, Cohen ML, Louie SG Electron-phonon interaction via electronic and lattice Wannier functions: Superconductivity in boron-doped diamond reexamined PHYSICAL REVIEW LETTERS, submitted.
[41] Giustino F, Pasquarello A Mixed Wannier-Bloch functions for electrons and phonons in periodic systems PHYSICAL REVIEW LETTERS 96, 216403 (2006).
[42] Brouder C, Panati G, Calandra M, Mourougane C, Marzari N Exponential localization of Wannier functions in insulators cond-mat/0606726.
[43] He LX, Vanderbilt D Exponential decay properties of Wannier functions and related quantities PHYSICAL REVIEW LETTERS 86, 5341 (2001).
[44] Sharma M, Wu YD, Car R Ab initio molecular dynamics with maximally localized Wannier functions INTERNATIONAL JOURNAL OF QUANTUM CHEMISTRY 95, 821 (2003).
[45] Thomas JW, Iftimie R, Tuckerman ME Field theoretic approach to dynamical orbital localization in ab initio molecular dynamics PHYSICAL REVIEW B 69, 125105 (2004).
[46] Iftimie R, Thomas JW, Tuckerman ME On-the-fly localization of electronic orbitals in Car-Parrinello molecular dynamics JOURNAL OF CHEMICAL PHYSICS 120, 2169 (2004).
[47] http://www.wannier.org/
[48] Silvestrelli PL Maximally localized Wannier functions for simulations with supercells of general symmetry PHYSICAL REVIEW B 59, 9703 (1999).
[49] Berghold G, Mundy CJ, Romero AH, Hutter J, Parrinello M General and efficient algorithms for obtaining maximally localized Wannier functions PHYSICAL REVIEW B 61, 10040 (2000).
[50] Gygi F, Fattebert JL, Schwegler E Computation of Maximally Localized Wannier Functions using a simultaneous diagonalization algorithm COMPUTER PHYSICS COMMUNICATIONS 155, 1 (2003).
[51] Maximally localized Wannier functions from PAW or ultrasoft pseudopotentials Ferretti A, Calzolari A, Bonferroni B, Di Felice R cond-mat/0603256.
[52] Whittaker DM, Croucher MP Maximally localized Wannier functions for photonic lattices PHYSICAL REVIEW B 67, 085204 (2003).