N-body simulation: Difference between revisions

From formulasearchengine
Jump to navigation Jump to search
en>Sjzaslaw
The Fourier transform of the poisson equation in the particle mesh method was not compiling properly. Fixed
 
(One intermediate revision by one other user not shown)
Line 1: Line 1:
In [[condensed matter physics]] and [[crystallography]], the '''static structure factor''' (or '''structure factor''' for short) is a mathematical description of how a material scatters incident radiation. The structure factor is a particularly useful tool in the interpretation of [[interference pattern]]s obtained in [[X-ray diffraction|X-ray]], [[Electron diffraction|electron]] and [[neutron diffraction|neutron]] [[diffraction]] experiments.
I would like to introduce myself to you, I am Andrew and my spouse doesn't like it at all. Ohio is where her house is. My day job is an information officer but I've currently utilized for an additional 1. What I adore doing is football but I don't have the time recently.<br><br>Also visit my blog post :: [http://www.indosfriends.com/profile-253/info/ spirit messages]
 
The static structure factor is measured without resolving the energy of scattered photons/electrons/neutrons. Energy-resolved measurements yield the [[dynamic structure factor]].
 
==Derivation==
Let us consider a scalar (real) quantity <math>\phi(\mathbf{r})</math> defined in a volume <math>V</math>; it may correspond, for instance, to a mass or charge distribution or to the refractive index of an inhomogeneous medium. If the scalar function is assumed to be integrable, we can define its [[Fourier transform]] <math>\textstyle \phi(\mathbf{q}) = \int_{V} \phi(\mathbf{r}) \exp (-i \mathbf{q} \mathbf{r}) \, \mathrm{d} \mathbf{r}</math>. Expressing the field <math>\phi</math> in terms of the [[spatial frequency]] <math>\mathbf{q}</math> instead of the point position <math>\mathbf{r}</math> is very useful, for instance, when interpreting [[scattering]] experiments. Indeed, in the [[Born approximation]] (weak interaction between the field and the medium), the amplitude of the signal corresponding to the scattering vector <math>\mathbf{q}</math> is proportional to <math>\textstyle \phi(\mathbf{q})</math>. Very often, only the intensity of the scattered signal <math>\textstyle I(\mathbf{q})</math> is detectable, so that <math>\textstyle I(\mathbf{q}) \sim \left | \phi(\mathbf{q}) \right |^2</math>.
 
If the system under study is composed of a number <math>N</math> of identical constituents (atoms, molecules, colloidal particles, etc.) it is very convenient to explicitly capture the variation in <math>\phi</math> due to the morphology of the individual particles using an auxiliary function <math>f(\mathbf{r})</math>, such that:
 
{{NumBlk|:|<math>\phi(\mathbf{r}) = \sum_{j=1}^{N} f(\mathbf{r} - \mathbf{R}_{j}) = f(\mathbf{r}) \ast \sum_{j=1}^{N} \delta(\mathbf{r} - \mathbf{R}_{j})</math>,|{{EquationRef|1}}}}
 
with <math>\textstyle \mathbf{R}_{j}, j = 1, \, \ldots, \, N</math> the particle positions. In the second equality, the field is decomposed as the [[convolution]] product <math>\ast</math> of the function <math>f</math>, describing the "form" of the particles, with a sum of [[Dirac delta function]]s depending only on their positions. Using the property that the Fourier transform of a convolution product is simply the product of the Fourier transforms of the two factors, we have <math>\textstyle \phi(\mathbf{q}) = f(\mathbf{q}) \sum_{j=1}^{N} \exp (-i \mathbf{q} \mathbf{R}_{j})</math>, such that:
 
{{NumBlk|:|<math>I(\mathbf{q}) \sim \left | \phi(\mathbf{q}) \right |^2 = \left | f(\mathbf{q}) \right |^2 \times \left ( \sum_{j=1}^{N} \mathrm{e}^{-i \mathbf{q} \mathbf{R}_{j}} \right ) \times \left ( \sum_{k=1}^{N} \mathrm{e}^{i \mathbf{q} \mathbf{R}_{k}} \right )= \left | f(\mathbf{q}) \right |^2 \sum_{jk} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_j - \mathbf{R}_k)}</math>.|{{EquationRef|2}}}}
 
In general, the particle positions are not fixed and the measurement takes place over a finite exposure time and with a macroscopic sample (much larger than the interparticle distance). The experimentally accessible intensity is thus an averaged one <math>\textstyle \langle I(\mathbf{q}) \rangle</math>; we need not specify whether <math>\langle \cdot \rangle</math> denotes a time or [[ensemble average]]. We can finally write:
 
{{NumBlk|:|<math> \langle I(\mathbf{q}) \rangle \sim \langle \left | \phi(\mathbf{q}) \right |^2 \rangle = N \left | f(\mathbf{q}) \right |^2 S(\mathbf{q})</math>,|{{EquationRef|3}}}}
 
thus defining the '''structure factor'''
{{NumBlk|:|<math>S(\mathbf{q}) = \frac{1}{N} \left \langle \sum_{jk} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_j - \mathbf{R}_k)} \right \rangle</math>.|{{EquationRef|4}}}}
 
==Perfect crystals==
In a [[crystal]], the constitutive particles are arranged periodically, forming a [[crystal lattice|lattice]]. In the following, we will consider that all particles are identical (so the above separation in factor and structure factors ({{EquationNote|3}}) holds). We also assume that all atoms have an identical environment (i.e. they form a [[Bravais lattice]]). The general case of lattice with a basis (see below) is not fundamentally different.
 
If the lattice is infinite and completely regular, the system is a [[perfect crystal]]. In addition, we will neglect all thermal motion, so that there is no need for averaging in ({{EquationNote|4}}). As in ({{EquationNote|2}}), we can write:
: <math>S(\mathbf{q}) = \frac{1}{N} \left | \sum_{j=1}^{N} \mathrm{e}^{-i \mathbf{q} \mathbf{R}_{j}} \right | ^2</math>.
 
The structure factor is simply the squared modulus of the [[Fourier transform]] of the lattice, and it is itself a periodic arrangement of points, known as the [[reciprocal lattice]].
 
===One dimension===
 
[[File:Sq linear.svg|thumb|Structure factor of a periodic chain, for different particle numbers <math>N</math>.]]
 
The reciprocal lattice is easily constructed in one dimension: for particles on a line with a period <math>a</math>, the atom positions <math>\textstyle R_j =  a (j - (N-1)/2)</math> (for simplicity, we consider that <math>N</math> is odd). The sum of the phase factors is a simple geometric series and the structure factor becomes:
 
: <math>S(q) = \frac{1}{N} \left | \frac{1 - \mathrm{e}^{-i N q a}}{1 - \mathrm{e}^{-i q a}} \right | ^2 =
\frac{1}{N} \left [ \frac{\sin(N q a/2)}{\sin(q a/2)} \right ] ^2 </math>.
 
This function is shown in the Figure below for different values of <math>N</math>.
 
Based on this expression for <math>S(q)</math>, one can draw several conclusions: the reciprocal lattice has a spacing <math>2\pi/a</math>; the intensity of the maxima increases with the number of particles <math>S(q = 2 k \pi/a) = N</math> (this is apparent from the Figure and can be shown by estimating the limit <math>S(q \to 0)</math> using, for instance, [[L'Hôpital's rule]]); the intensity at the midpoint <math>S(q = (2 k +1) \pi/a) = 1/N</math> (by direct evaluation); the peak width also decreases like <math>1/N</math>. In the large <math>N</math> limit, the peaks become infinitely sharp Dirac delta functions.
 
===Two dimensions===
 
[[File:square lattice scattering.png|thumb|Diagram of scattering by a square (planar) lattice. The incident and outgoing beam are shown, as well as the relation between their wave vectors <math>\mathbf{k}_i</math>, <math>\mathbf{k}_o</math> and the scattering vector <math>\mathbf{q}</math>.]]
 
In two dimensions, there are only five [[Bravais lattice]]s. The corresponding reciprocal lattices have the same symmetry as the direct lattice. The Figure shows the construction of one vector of the reciprocal lattice and its relation with a scattering experiment.
 
A parallel beam, with wave vector <math>\mathbf{k}_i</math> is incident on a square lattice of parameter <math>a</math>. The scattered wave is detected at a certain angle, which defines the wave vector of the outgoing beam, <math>\mathbf{k}_o</math> (under the assumption of [[elastic scattering]], <math>|\mathbf{k}_o| = |\mathbf{k}_i|</math>). One can equally define the scattering vector <math>\mathbf{q}=\mathbf{k}_o - \mathbf{k}_i</math> and construct the harmonic pattern <math>\exp (i \mathbf{q}\mathbf{r})</math>. In the depicted example, the spacing of this pattern coincides to the distance between particle rows: <math>q = 2\pi /a</math>, so that contributions to the scattering from all particles are in phase (constructive interference). Thus, the total signal in direction <math>\mathbf{k}_o</math> is strong, and <math>\mathbf{q}</math> belongs to the reciprocal lattice. It is easily shown that this configuration fulfills [[Bragg's law]].
 
===Three dimensions===
 
===Lattice with a basis===
{{Confusing|date=May 2013}}
To compute structure factors for a specific lattice, compute the sum above over the atoms in the unit cell. Since crystals are often described in terms of their [[Miller index|Miller indices]], it is useful to examine a specific structure factor in terms of these.
 
'''Body-centered cubic (BCC) '''
 
As a convention, the [[body-centered cubic]] system is described in terms of a simple cubic lattice with primitive vectors <math>a\hat{x}, a\hat{y}, a\hat{z}</math>, with a basis consisting of <math>\mathbf{r}_0 = \vec{0}</math> and <math>\mathbf{r}_1 = (a/2)(\hat{x} + \hat{y} + \hat{z})</math>. The corresponding reciprocal lattice is also simple cubic with side <math>2\pi/a</math>.
 
In a monatomic crystal, all the form factors <math>f</math> are the same.    The intensity of a diffracted beam scattered with a vector <math>\mathbf{K}=h\hat{x}^* + k\hat{y}^* + l\hat{z}^*=(2\pi/a)(h\hat{x} + k\hat{y} + l\hat{z})</math> by a crystal plane with Miller indices <math>(hkl)</math> is then given by:
 
:<math>\begin{matrix}
F_{\mathbf{K}} & = & f \left[ e^{-i\mathbf{K}\cdot\vec{0}} + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{y} + \hat{z})} \right] \\
& = & f \left[ 1 + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{y} + \hat{z})} \right] \\
& = & f \left[ 1 + e^{-i\pi(h + k + l)} \right]\\
& = & f \left[ 1 + (-1)^{h + k + l} \right] \\
\end{matrix}</math>
 
We then arrive at the following result for the structure factor for scattering from a plane <math>(hkl)</math>:
 
<math>
F_{hkl} =  \begin{cases} 2f, & h + k + l \ \ \mbox{even}\\
                    0, & h + k + l \ \  \mbox{odd} \end{cases}
</math>
 
This result tells us that for a reflection to appear in a diffraction experiment involving a body-centered crystal, the sum of the Miller indices of the scattering plane must be even. If the sum of the Miller indices is odd, the intensity of the diffracted beam is reduced to zero due to destructive interference. This zero intensity for a group of diffracted beams is called a systematic absence. Since atomic form factors fall off with increasing diffraction angle corresponding to higher Miller indices, the most intense diffraction peak from a material with a BCC structure is typically the (110).    The (110) plane is the most densely packed of BCC crystal structures and is therefore the lowest energy surface for a thin film to grow.  Films of BCC materials like [[iron]] and [[tungsten]] therefore grow in  a characteristic (110) orientation.
 
'''Face-centered cubic (FCC) '''
 
In the case of a monatomic FCC crystal, the atoms in the basis are at the origin <math>\mathbf{r}_0 = \vec{0}</math> with indices (0,0,0) and at the three face centers <math>\mathbf{r}_1 = (a/2)(\hat{x} + \hat{y})</math>,  <math>\mathbf{r}_2 = (a/2)(\hat{y} + \hat{z})</math>,  <math>\mathbf{r}_3 = (a/2)(\hat{x} + \hat{z})</math> with indices given by (1/2,1/2,0), (0,1/2,1/2), (1/2,0,1/2).   An argument similar to the one above gives the expression
 
:<math>\begin{matrix}
F_{\mathbf{K}} & = & f \left[ e^{-i\mathbf{K}\cdot\vec{0}} + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{y})} + e^{-i\mathbf{K}\cdot(a/2)(\hat{y} + \hat{z})} + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{z})} \right] \\
& = & f \left[ 1 + (-1)^{h + k} + (-1)^{k + l} + (-1)^{h + l} \right] \\
\end{matrix}</math>
 
with the result
 
<math>
F_{hkl} =  \begin{cases} 4f, & h,k,l  \ \ \mbox{all even or all odd}\\
                    0, & h,k,l \ \  \mbox{mixed parity} \end{cases}
</math>
 
The most intense diffraction peak from a material that crystallizes in the FCC structure is typically the (111).  Films of FCC materials like [[gold]] tend to grow in a (111) orientation with a triangular surface symmetry.
 
'''Diamond Crystal Structure'''
 
The [[Diamond cubic]] crystal structure occurs in [[diamond]] ([[carbon]]), most [[semiconductors]] and [[tin]]. The basis cell contains 8 atoms located at cell positions:
 
<math>\mathbf{r}_0 = \vec{0}</math>
 
<math>\mathbf{r}_1 = (a/4)(\hat{x} + \hat{y} + \hat{z})</math>
 
<math>\mathbf{r}_2 = (a/4)(2\hat{x} + 2\hat{y})</math>
 
<math>\mathbf{r}_3 = (a/4)(3\hat{x} + 3\hat{y} + \hat{z})</math>
 
<math>\mathbf{r}_4 = (a/4)(2\hat{x} + 2\hat{z})</math>
 
<math>\mathbf{r}_5 = (a/4)(2\hat{y} + 2\hat{z})</math>
 
<math>\mathbf{r}_6 = (a/4)(3\hat{x} + \hat{y} + 3\hat{z})</math>
 
<math>\mathbf{r}_7 = (a/4)(\hat{x} + 3\hat{y} + 3\hat{z})</math>
 
The Structure factor then takes on a form like this:
 
:<math>\begin{matrix}
F_{\mathbf{K}} & = & f \left[
\begin{matrix}
e^{-i\mathbf{K}\cdot\vec{0}} + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{y})} + e^{-i\mathbf{K}\cdot(a/2)(\hat{y} + \hat{z})} + e^{-i\mathbf{K}\cdot(a/2)(\hat{x} + \hat{z})} + \\
e^{-i\mathbf{K}\cdot(a/4)(\hat{x} + \hat{y} + \hat{z})} +
e^{-i\mathbf{K}\cdot(a/4)(3\hat{x} + \hat{y} + 3\hat{z})} +
e^{-i\mathbf{K}\cdot(a/4)(3\hat{x} + 3\hat{y} + \hat{z})} +
e^{-i\mathbf{K}\cdot(a/4)(\hat{x} + 3\hat{y} + 3\hat{z})}
\end{matrix}
\right] \\
& = & f \left[
\begin{matrix}
1 + (-1)^{h + k} + (-1)^{k + l} + (-1)^{h + l} + \\
(-i)^{h + k + l} + (-i)^{3h + k + 3l} + (-i)^{3h + 3k + l} + (-i)^{h + 3k + 3l}
\end{matrix}
\right] \\
& = & f \left[ 1 + (-1)^{h + k} + (-1)^{k + l} + (-1)^{h + l} \right] \cdot \left[ 1 + (-i)^{h + k + l} \right]\\
\end{matrix}
</math>
 
with the result
 
* for mixed values (odds and even values combined) of h, k, and l, F<sup>2</sup> will be 0
* if the values are unmixed and...
**h+k+l is odd then F=4f(1+i) or 4f(1-i), FF<sup>*</sup>=32f<sup>2</sup>
**h+k+l is even and exactly divisible by 4 (satisfies h+k+l=4n) then ''F'' = 8''f''
**h+k+l is even but not exactly divisible by 4(doesn't satisfy h+k+l=4n) then ''F'' = 0
 
==Imperfect crystals==
Although the perfect lattice is an extremely useful model, real crystals always exhibit imperfections, which can have profound effects on the structure and properties of the material. [[André Guinier]] <ref>See Guinier, chapters 6-9</ref> proposed a widely employed distinction between imperfections that preserve the [[long-range order]] of the crystal (disorder of the first kind) and those that destroy it (disorder of the second kind).
 
===Disorder of the first kind===
 
===Disorder of the second kind===
 
==Liquids==
In contrast with crystals, liquids have no [[long-range order]] (in particular, there is no regular lattice), so the structure factor does not exhibit sharp peaks. They do however show a certain degree of [[short-range order]], depending on their density and on the strength of the interaction between particles. Liquids are isotropic, so that, after the averaging operation in Equation ({{EquationNote|4}}), the structure factor only depends on the absolute magnitude of the scattering vector <math>q = \left |\mathbf{q} \right |</math>. For further evaluation, it is convenient to separate the diagonal terms <math>j = k</math> in the double sum, whose phase is identically zero, and therefore each contribute a unit constant:
 
{{NumBlk|:|<math>S(q) = 1 + \frac{1}{N} \left \langle \sum_{j \neq k} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_j - \mathbf{R}_k)} \right \rangle</math>.|{{EquationRef|5}}}}
 
One  can obtain an alternative expression for <math>S(q)</math> in terms of the [[radial distribution function]] <math>g(r)</math>:<ref>See Chandler, section 7.5.</ref>
 
{{NumBlk|:| <math>S(q) = 1 + \rho \int_V \mathrm{d} \mathbf{r} \, \mathrm{e}^{-i  \mathbf{q}\mathbf{r}} g(r)</math>.|{{EquationRef|6}}}}
 
===Ideal gas===
In the limiting case of no interaction, the system is an [[ideal gas]] and the structure factor is completely featureless: <math>S(q) = 1</math>, because there is no correlation between the positions <math>\mathbf{R}_j</math> and <math>\mathbf{R}_k</math> of different particles (they are [[independent random variables]]), so the off-diagonal terms in Equation ({{EquationNote|5}}) average to zero: <math>\langle \exp [-i \mathbf{q} (\mathbf{R}_j - \mathbf{R}_k)]\rangle = \langle \exp (-i \mathbf{q} \mathbf{R}_j) \rangle \langle \exp (i \mathbf{q} \mathbf{R}_k) \rangle = 0</math>.
 
===High-<math>q</math> limit===
Even for interacting particles, at high scattering vector the structure factor goes to 1. This result follows from Equation ({{EquationNote|6}}), since <math>S(q)-1</math> is the [[Fourier transform]] of the "regular" function <math>g(r)</math> and thus goes to zero for high values of the argument <math>q</math>. This reasoning does not hold for a perfect crystal, where the distribution function exhibits infinitely sharp peaks.
 
===Low-<math>q</math> limit===
In the low-<math>q</math> limit, as the system is probed over large length scales, the structure factor contains thermodynamic information, being related to the [[isothermal compressibility]] <math>\chi _T</math> of the liquid by the [[compressibility equation]]:
: <math>\lim _{q \rightarrow 0} S(q) = \rho kT \chi _T= kT\left(\frac{\partial \rho}{\partial p}\right)</math>.
 
===Hard-sphere liquids===
[[File:HS structure factor PY.svg|thumb|Structure factor of a hard-sphere fluid, calculated using the Percus-Yevick approximation, for volume fractions <math>\Phi</math> from 1% to 40%.]]
In the [[hard sphere]] model, the particles are described as impenetrable spheres with radius <math>R</math>; thus, their center-to-center distance <math>r \geq 2R</math> and they experience no interaction beyond this distance. Their interaction potential can be written as:
: <math> V(r) = \left \lbrace
  \begin{array}{l l}
    \infty \, &\text{for} \,\, r < 2 R\\
    0\, &\text{for} \,\, r \geq 2 R\\
  \end{array}
\right .</math>
 
This model has an analytical solution<ref>{{Cite doi|10.1103/PhysRevLett.10.321}}</ref> in the [[Percus–Yevick approximation]]. Although highly simplified, it provides a good description for systems ranging from liquid metals<ref>{{Cite doi|10.1103/PhysRev.145.83}}</ref> to colloidal suspensions.<ref>{{Cite doi|10.1038/320340a0}}</ref> In an illustration, the structure factor for a hard-sphere fluid is shown in the Figure, for volume fractions <math>\Phi</math> from 1% to 40%.
 
==Polymers==
In [[polymer]] systems, the general definition ({{EquationNote|4}}) holds; the elementary constituents are now the [[monomer]]s making up the chains. However, the structure factor being a measure of the correlation between particle positions, one can reasonably expect that this correlation will be different for monomers belonging to the same chain or to different chains.
 
Let us assume that the volume <math>V</math> contains <math>N_c</math> identical molecules, each composed of <math>N_p</math> monomers, such that <math>N_c N_p = N</math> (<math>N_p</math> is also known as the [[degree of polymerization]]). We can rewrite ({{EquationNote|4}}) as:
{{NumBlk|:|<math>S(\mathbf{q}) = \frac{1}{N_c N_p} \left \langle \sum_{\alpha \beta = 1}^{N_c} \sum_{jk = 1}^{N_p} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_{\alpha j} - \mathbf{R}_{\beta k})} \right \rangle = \frac{1}{N_c N_p} \left \langle \sum_{\alpha = 1}^{N_c} \sum_{jk = 1}^{N_p} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_{\alpha j} - \mathbf{R}_{\alpha k})} \right \rangle + \frac{1}{N_c N_p} \left \langle \sum_{\alpha \neq \beta = 1}^{N_c} \sum_{jk = 1}^{N_p} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_{\alpha j} - \mathbf{R}_{\beta k})} \right \rangle</math>,|{{EquationRef|7}}}}
where indices <math>\alpha , \beta</math> label the different molecules and <math>j, k</math> the different monomers along each molecule. On the right-hand side we separated ''intramolecular'' (<math>\alpha = \beta</math>) and ''intermolecular'' (<math>\alpha \neq \beta</math>) terms. Using the equivalence of the chains, ({{EquationNote|7}}) can be simplified:<ref>See Teraoka, Section 2.4.4.</ref>
{{NumBlk|:|<math>S(\mathbf{q}) = \underbrace{\frac{1}{N_p} \left \langle \sum_{jk = 1}^{N_p} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_{j} - \mathbf{R}_{k})} \right \rangle}_{S_1(q)} + \frac{N_c - 1}{N_p} \left \langle \sum_{jk = 1}^{N_p} \mathrm{e}^{-i \mathbf{q} (\mathbf{R}_{1 j} - \mathbf{R}_{2 k})} \right \rangle</math>,|{{EquationRef|8}}}}
where <math>S_1 (q)</math> is the single-chain structure factor.
 
==See also==
*[[R-factor (crystallography)]]
*[[Patterson function]]
 
==Notes==
{{Reflist}}
 
==References==
# Als-Nielsen, N. and McMorrow, D. (2011). Elements of Modern X-ray Physics (2nd edition). John Wiley & Sons.
#[[André Guinier|Guinier, A.]] (1963). X-ray Diffraction. In Crystals, Imperfect Crystals, and Amorphous Bodies. W. H. Freeman and Co.
#[[David Chandler (chemist)|Chandler, D.]] (1987). [[Introduction to Modern Statistical Mechanics]]. Oxford University Press.
#[[Jean-Pierre Hansen|Hansen, J. P.]] and McDonald, I. R. (2005). Theory of Simple Liquids (3rd edition). Academic Press.
#Teraoka, I. (2002). Polymer Solutions: An Introduction to Physical Properties. John Wiley & Sons.
 
==External links==
*[http://www.ysbl.york.ac.uk/~cowtan/sfapplet/sfintro.html Structure Factor Tutorial] located at the [[University of York]].
 
[[Category:Crystallography]]

Latest revision as of 01:59, 21 December 2014

I would like to introduce myself to you, I am Andrew and my spouse doesn't like it at all. Ohio is where her house is. My day job is an information officer but I've currently utilized for an additional 1. What I adore doing is football but I don't have the time recently.

Also visit my blog post :: spirit messages