Physics 250-06 “Advanced Electronic Structure”
Download
1 / 25

Physics 250-06 “Advanced Electronic Structure” Pseudopotentials Contents: - PowerPoint PPT Presentation


  • 194 Views
  • Uploaded on

Physics 250-06 “Advanced Electronic Structure” Pseudopotentials Contents: 1. Plane Wave Representation 2. Solution for Weak Periodic Potential 3. Solution for a single scattering problem 4. Pseudopotential Approximation. Solving Schroedingers Equation via Fourier Transoforms.

loader
I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.
capcha
Download Presentation

PowerPoint Slideshow about ' Physics 250-06 “Advanced Electronic Structure” Pseudopotentials Contents:' - vina


An Image/Link below is provided (as is) to download presentation

Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.


- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

Physics 250-06 “Advanced Electronic Structure”

Pseudopotentials

Contents:

1. Plane Wave Representation

2. Solution for Weak Periodic Potential

3. Solution for a single scattering problem

4. Pseudopotential Approximation


Solving Schroedingers Equation via Fourier Transoforms.

Periodic array of potential wells placed at distance a between the wells.

Basic property of the potential V(x+l)=V(x), and of the Hamiltonian H(x+l)=H(x), i.e

periodicity.

Periodic boundary condition is imposed

Properties of solutions: The solutions are traveling (or Bloch waves). There are infinite

number of solutions which can be labeled by wave vector K, i.e.

where t=m*l, m is any integer. Vectors k can take the values following from the

periodic boundary condition:

when n is any integer.


Non trivial values of K are those for which n runs from 0 to N, where N is number of wells considered,

i.e. L=N*l.Important properties of the Bloch waves: Orthogonality


Example of periodic potential

This potential is periodic function as for x=l*m, m=0,…so that it can be a model to

analyze the solutions of Schroedingers equation in periodic potential.

For V1=0 the solutions can be written as plane waves

The eigenvalues E are


Consider now the equation

Let us use a single plane wave as a basis:

The expansion coefficient is simply equal 1 (assuming plane waves are normalized to 1

within the volume). Let us find the eigenvalues. The matrix eigenvalue problem is collapsed to

which means that there is no correction to the free electrons eigenstates linear to V1.


We need to improve our basis! With the plane waves it is pretty easy.

For this recall first what is

where n is any integer. Let us agree on n running from 0 to N and this will be called k. Then,

where n is an integer from 0 to N and n’ is any integer. In other words,

so that

and note also that the latter is periodic function on the lattice

We thus separated the space of all wave vectors K on the subspace of so called irreducible

wave vectors k and the subspace of so called reciprocal lattice vectors G. (concept of Brillouin Zone)

This can always be done as soon as we introduce a period set by length l, and the

periodic boundary condition set by length L=Nl. Now we have a powerful basis set since we can expand

In other words, for each given k, we use the subspace G which delivers us basis functions. We use

those basis functions to represent the wave function for given wave vector k.


Let us return back to our example. Restrict the expansion by two nearest G, i.e. pick n’=-1,0,1, or

G=-2pi/l,0,+2pi/l. We have

where G0=2pi/l.

Compute matrix eigenvalue problem in this basis

Now we realize that plane waves with different K’s or different G’s are orthogonal and therefore

We also realize that diagonal element

does not depend on the oscillating part of the potential!We realize that off diagonal elements depend

only on oscillating part of the potential.because


The matrix elements are trivially evaluated two nearest G, i.e. pick n’=-1,0,1, or

We finally obtain

The roots can be found by looking at the determinant

Opening the determinant produces

This quation delivers us three roots E1(k), E2(k), E3(k).


Solution for Periodic Potential two nearest G, i.e. pick n’=-1,0,1, or


Scattering by a single potential two nearest G, i.e. pick n’=-1,0,1, or

Consider scattering by a potential assumed to be spherically

symmetric inside a sphere

Sphere


Solve radial Schroedinger equation inside the sphere two nearest G, i.e. pick n’=-1,0,1, or

Solve Helmholtz equation outside

the sphere

Sphere


Solution of Helmholtz equation outside two nearest G, i.e. pick n’=-1,0,1, or

the sphere

where coefficients provide smooth

matching with


A single L-partial wave solves single scattering problem two nearest G, i.e. pick n’=-1,0,1, or

Any linear combination solves it as well


A potential information is hidden inside the coefficients providing

smooth matching

One can replace real potential by a pseudopotential

which would provide the same scattering property!


Replace real potential for each L by its by its pseudopotential

should produce such solution that

In other words, (i) Solve Equation for , find

(ii) Adjust so that it produces pseudo

with the same scattering properties.


The Pseudopotential Approximation for Solids pseudopotential

Solving solid state problem using plane wave basis set

Unfortunately a plane wave basis set is usually very poorly suited to expanding the electronic wavefunctions because a very large number are required to accurately describe the rapidly oscillating wavefunctions of electrons in the core region.

Also we have a divergency problem in 3D ionic potential


It is well known that most physical properties of solids are dependent on the valence electrons to a much greater degree than that of the tightly bound core

electrons. It is for this reason that the pseudopotential approximation is introduced.

This approximation uses this fact to remove the core electrons and the strong nuclear potential and replace them with a weaker pseudopotential which acts on a set of pseudowavefunctions rather than the true valence wavefunctions.

So, for each atom in periodic table ionic potential can be replaced by pseudopotential which would produce the same scattering process:


A pseudopotential is not unique, therefore several methods of

generation exist. However they must obey several criteria.

The core charge produced by the pseudo wavefunctions

must be the same as that produced by the atomic wavefunctions.

This ensures that the pseudo atom produces the same

scattering properties as the ionic core.

Pseudopotentials of this type are known as non-local

norm-conserving pseudopotentials and are the most

transferable since they are capable of describing the

scattering properties of an ion in a variety of atomic environments.

Pseudo-electron eigenvalues must be the same as the

valence eigenvalues obtained from the real wavefunctions.


Mathematically, the angular dependent pseudopotential can be written as:

is an angular momentum projection operator. In computation, this can be evaluated

Here =

Analytical formula exists to write down the matrix elements. For large systems

he Kleinman-Bylander implementation of the nonlocal potential is used. Basically, we first take

one l (s, or, p, or d) as the local potential,

Then we can define a nonlocal part as

Then the above nonlocal pseudopotential is approximated as:

are the atomic pseudowavefunctions for angular l.

in G-space), or real-space. In the real-space representation, the projection function can be cut off beyond a radius rcut (which is used as a parameter in PEtot).

As a result, for large system calculations, real-space implementation of the Kleinman-Bylander form is faster than G-space implementation


Plane Wave Grids written as:

The sphere of Gc2 being inside the Fourier box of n1*n2*n3 grid.  


Gygi’s adaptive grid: ideas from the theory of gravity written as:

Francois Gygi (UC Davis)


ad