Chemical Dynamics/Electrostatics/Fourier Transforms

From testwiki
Revision as of 00:32, 12 February 2012 by imported>QUBot (Adding BookCat to uncategorised page)
(diff) ← Older revision | Latest revision (diff) | Newer revision β†’ (diff)
Jump to navigation Jump to search

The Fourier transform is a useful mathematical transformation often utilized in many scientific and engineering fields. Here we extract useful concepts of Fourier transformation and logically arrange them to form a foundation for the Ewald summation and other related methods in electrostatics. Readers could check out other more mathematically formal introduction of Fourier transform

Definition

We use the following convention in which the Fourier transform is a unitary transformation on the 3-D Cartesian space R3, the Fourier transform and its inverse transform are symmetric:

f^(𝐀)=1(2π)3/2f(𝐫)ei𝐀𝐫d3𝐫
f(𝐫)=1(2π)3/2f^(𝐀)ei𝐀𝐫d3𝐀

The translation theorem

Given a fixed position vector R0, if g(r) = Ζ’(r βˆ’ R0), then  

g^(𝐀)=ei𝐀𝐑0f^(𝐀).

Template:Hidden begin

g^(𝐀)=1(2π)3/2g(𝐫)ei𝐀𝐫d3𝐫
=1(2π)3/2f(𝐫𝐑0)ei𝐀𝐫d3𝐫
=1(2π)3/2f(𝐫𝐑0)ei𝐀𝐫d3𝐫

Now, change r to a new variable by: 𝐫=𝐫𝐑0

g^(𝐀)=1(2π)3/2f(𝐫)ei𝐀(𝐫+𝐑0)d3𝐫
=1(2π)3/2ei𝐀𝐑0f(𝐫)ei𝐀𝐫d3𝐫
=1(2π)3/2ei𝐀𝐑0f(𝐫)ei𝐀𝐫d3𝐫
=ei𝐀𝐑0f^(𝐀).

Template:Hidden end

The convolution theorem

The convolution of f and g is usually denoted as fβˆ—g, using an asterisk or star. It is defined as the integral of the product of the two functions after one is reversed and shifted:


(f*g)(t)=def f(τ)g(tτ)dτ


The convolution theorem for the Fourier transform says:

If

h(𝐫)=(f*g)(𝐫)

then

h^(𝐀)=f^(𝐀)g^(𝐀).

Template:Hidden begin

h^(𝐀)=ei𝐀𝐫h(𝐫)d𝐫
=ei𝐀𝐫f(𝐫)g(𝐫𝐫)d3𝐫d3𝐫

Template:Hidden end Template:BookCat