Review Article | Open Access
Topological Excitations in Quantum Spin Systems
The origin and significance of topological excitations in quantum spin models in low dimensions are presented in detail. Besides a general review, our own work in this area is described in great depth. Apart from theoretical analysis of the existence and properties of spin vortices and antivortices, the possible experimental consequences and signatures are also highlighted. In particular, the distinguishing features between the even and odd charged topological excitations are brought out through a detailed analysis of the topological term in the quantum action. Moreover, an interesting symmetry property is predicted between the excitations from a ferromagnetic model and an antiferromagnetic model. Through a novel approach of ours, a bridge is established between field theoretical formalism and the well-known statistical mechanical treatment of Berezinskii-Kosterlitz-Thouless (BKT) transition involving these topological excitations. Furthermore, a detailed phenomenological analysis of the experimentally observed static and dynamic magnetic properties of the layered magnetic materials, possessing XY anisotropy in the in-plane spin-spin couplings, is undertaken to test the theoretical predictions regarding the behaviour of these excitations. The importance and the crucial role of quantum spin fluctuations in these studies are also brought out very clearly by our analysis.
The question of existence of nonconventional excitations in spin systems has been probed for the last four decades. The seminal papers of Kosterlitz and Thouless (1973) and by Berezinskii (1971, 1972) laid a foundation for this investigation [1–3]. The models considered there involved classical spins situated on two dimensional (2D) spatial lattice. Since then, developments in various branches of condensed matter physics have prompted investigation with quantum spin models in low-dimensional systems as well.
Recent progress in the field of low-dimensional strongly correlated electronic systems and in particular their magnetic and transport properties both from theoretical and experimental standpoints has stimulated and motivated the urgent need for the search of non-conventional spin excitations in purely quantum systems [4–8].
Many of these quantum systems exhibit topological character in their excitations, single particle and/or collective. It is therefore very important that a proper understanding of their origin and properties is taken up.
The collective excitations in quantum spin systems can broadly be divided into two categories: (i) conventional ones, namely, spin waves (quantised as magnons) and (ii) topological ones, namely, vortices and merons. The spin wave modes exist in the magnetically long range ordered phases of a Heisenberg magnet and are very stable and well defined in the low temperature regime. They are low energy excitations and represent a coherent precessional motion of all the spins around the direction of spontaneous magnetization on a lattice. In the short range ordered phase, these spin wave like modes become highly damped and fragile. In contrast, vortices and merons are high energy excitations and generally occur at higher temperatures. They exhibit large deformations and twists of the planar components and slow tilting of the out of plane component of the spins over a finite region of the lattice. Besides, the spin configurations contained in these topological excitations have the property that they cannot be shrunk to a point.
As a first step towards the understanding of these topological excitations, it would be appropriate to investigate the two-dimensional quantum XY or XY-anisotropic Heisenberg models. With the results of the analysis for the corresponding classical models being known, it is very instructive as well as natural to probe the quantum effects arising particularly in the low spin cases [9–11]. A pertinent step in this direction would be a semiclassical approach based on phenomenology. Parallel to this, a full fledged quantum treatment involving combination of field theoretic approach and conventional many body techniques can also be pursued.
We had originally adopted a phenomenological approach to quantum Heisenberg ferromagnet and antiferromagnet on square lattice in two dimensions . For K2CuF4, a spin- ferromagnet with small XY anisotropy, the static results from experiments and classical numerical calculations agree quite well with the Berezinski-Kosterlitz-Thouless (BKT) proposed phenomenology [1–3]. For La2CuO4, however, a spin- antiferromagnet, the static results from the experiments and classical renormalization group calculations exhibit a substantial disagreement with BKT-like phenomenology. The spin dynamics display even more interesting features in the above two systems. The inelastic neutron scattering experiments performed on the above two materials show a “central peak”, that is, a peak at in the constant -scan for the dynamical structure function for both of them [13–16]. In the first system it occurs in the low -regime, whereas for the second material, it takes place somewhat in the midway towards the zone boundary.
The application of the BKT inspired phenomenological treatment reproduces the central peak as well for the first system; however, the dynamical structure factor becomes negative in a large regime of the -space and hence exhibits unphysical behaviour . This signals the failure and limitation of the semiclassical treatment involving BKT (see Figure 1) and calls for a fully fledged quantum treatment . In the second system, the results from BKT-like phenomenology for the experimental observations from inelastic neutron scattering experiment involving staggered spin field and the RG calculational results available for the very weak XY-anisotropy limit disagree quantitatively by a large factor (about 4) . This calls for a thorough investigation. The experimental results for static and dynamic correlations are presented in Figures 2 and 3, respectively . The origin of this discrepancy might as well be in the presence of enormous quantum fluctuations in spin quantum antiferromagnet on 2D lattice and the failure of the classical/semiclassical treatment. This is more vigorous than that expected for a ferromagnet on a 2D lattice. However, at zero temperature, the quantum simulation performed by several theoretical groups brings out the existence of topological excitations (bound vortex-antivortex pair) even for an XY-like ferromagnet on 2D lattice [18–20]. On the other hand, even classical Monte-Carlo simulation allowing tunnelling between different topological sectors seems to wash out the BKT transition .
In view of the above developments and confusion, a long programme was taken up by the present authors to systematically analyse the interacting quantum spin models, namely, XY-anisotropic quantum Heisenberg spin systems, on low-dimensional lattices. Our method is based on the application of the path integral technique to derive an effective action corresponding to the spin Hamiltonian. The major aim was to explore whether a nonvanishing “topological term” arises in this effective action and if so, to what extent it indicates a possibility for the existence of topological excitations in the system. The presence of such a term in fact indicates a very strong possibility for the existence of topological excitations in the system and furthermore the numerical magnitude of this term (always an integer) represents the “winding number” for such configurations.
The approach we adopted is the following [12, 17, 22–26].(1)Calculate the partition function of the quantum spin model by making use of Trotter's formula and working in the spin coherent state basis.(2)Generate the corresponding “action” by taking a quasicontinuum approximation.(3)Rearrange and combine the various terms and contributions to extract out an effective nonlinear sigma model type part and a so-called topological part.(4)The sigma model part essentially represents spin wave modes, whereas the topological part implies the breaking up of the configurations into distinct topological sectors. For antiferromagnetic models, the above decomposition into two parts can be done easily. For ferromagnetic models, however, it is found that there can be another additional 3rd piece which is nontopological but has an interesting geometrical interpretation.(5)Evaluate the topological part both in the long wavelength and medium wavelength limits. In the case of nonvanishing contribution, the topological term has been explicitly calculated for vortices (antivortices) in an XY-anisotropic quantum Heisenberg ferromagnet in the low-q regime and antiferromagnet in the medium wavelength limit, on a square lattice.(6)When the net topological term is nonzero, we expect to see its signature in the partition function of the corresponding spin model. This also leads and motivates us to study the static and the dynamic spin-spin correlation functions in this model theoretically and compare them with the experimental results from materials which are very good realizations of these models.(7)Besides, various important physical and mathematical properties of the topological term and the topological excitations were studied. In fact, it was found that they carry enormous physical significance like magnitude of the vortex excitation energy, the thermodynamic feasibility of existence of these excitations, and last but not the least the characterisation of the even and odd winding number sectors. This evaluation of the properties backed by the suitable calculations also threw light on many more dramatic and fundamental issues, as will be discussed later.(8)For the time being for pragmatic reasons, the calculations for correlation functions have so far been based on combinations of both analytical and numerical techniques. Moreover, the phenomenological as well as the microscopic approaches have been made use of in this investigation. (9)The temperature dependence of the experimentally obtained static correlation functions has been studied thoroughly. The materials which are candidates for exhibiting BKT transitions are in reality very often layered or quasi-two dimensional. Moreover, the spin-spin couplings in the real materials often have additional contributions with the theoretical model being investigated. It is therefore very important to look into the right temperature regime to extract the pure spatial two-dimensional effects and those appropriate to the theoretical spin model.(10)The above procedure is followed for dynamical correlation function as well. In this case, the detailed energy dependence of dynamical structure factor in constant -scans from experiments is analysed very thoroughly. The instrumental resolution function is also taken into account while making comparison with the experimental data from inelastic neutron scattering experiments.
3. Analytical Description of the Quantum Action
It is algebraically quite complicated to write down the partition function of a quantum spin system in low dimensions containing huge amount of quantum fluctuations, using conventional angular momentum basis. However, making use of coherent spin state basis, one can get around these difficulties. Furthermore, one can formulate a path integral representation of the quantum partition function by introducing coherent spin states at each lattice point . It is only this formulation where the topological term, that is, the Wess-Zumino (WZ) type term, becomes manifested explicitly in the action corresponding to the spin system. This formulation was first established by Fradkin and Stone in the long wave length limit of an isotropic quantum Heisenberg antiferromagnetic chain . In this limit, the system takes the form of a nonlinear sigma model with a Hopf term, that is, the WZ term in 1 + 1 dimensions, and distinguishes the striking difference between the spectra corresponding to integer and half-integer spins [28, 29].
Making use of Trotter formula and inserting resolution of identity , one can write down the partition function as (chapter 5 in  and references therein).
, where the quantum Euclidean action for the coherent spin fields is given by where “” is the magnitude of the spin ( in the present case), is the position vector of a lattice site, and is the expectation value of the spin Hamiltonian () ((3) and (6) below) in a coherent state on the lattice [23, 24].
The variable denotes the pseudo-time (Euclidean time) appropriate to the above fields and has the dimension of inverse temperature. The parameter stands for as usual; is the real thermodynamic temperature of the spin system. The quantity is the Wess-Zumino term (WZ term) (see (4) below), corresponding to a single spatial lattice point at a fixed time . At this stage, let us mention that the WZ term in the action is analytically obtained in this way and thus it predicts the possibility of topological excitations in the quantum spin system. The spin operator is connected with the coherent spin variable in the following fashion : where denotes the direct product of all coherent spin states over the spatial lattice and stands for a coherent state on the lattice. The variable denotes the pseudo-time having the dimension of inverse temperature with as usual; is the real thermodynamic temperature of the spin system. The WZ term corresponding to a single spatial lattice point at a fixed time “” is given as follows: with , , and , , , and is the usual inverse temperature.
The expression in (4) is the area of the cap bounded by the trajectory parametrized by on the sphere: Furthermore, the fields are the fields in the higher dimensional -space and the boundary values are the coherent spin fields. The field is the fixed point on the above sphere and the state vector appearing on the right-hand side of (3) is the spin coherent state at a single lattice point . For the two-dimensional lattice, we express as for brevity, with being the lattice spacing. The spin Hamiltonian is given by where corresponding to anisotropic Heisenberg spin system of type with ferromagnetic coupling and corresponding to that with antiferromagnetic coupling. We represent here the explicit calculational procedure for the quantity given by (4).
It is possible to evaluate only the “difference” of terms from two neighbouring lattice sites in terms of the coherent spin fields . Thus, to extract the topological-like contribution from , we use the following expressions for the above “difference” : where with , and .
Now, the total Wess-Zumino term on the lattice is as given in the right-hand side of (2) and has to be evaluated in the case of ferro- or antiferromagnet.
It is in the case of antiferromagnet that we can calculate the term (see (9)) exactly since it appears with difference of two terms in the neighbouring lattice points owing to staggering. This is not so in the case of ferromagnet since we will not get alternating positive and negative signs in the the lattice sum of the WZ term as there is no staggering. However, we can extract the topological information content from the Wess-Zumino term as demonstrated by (10) and (11) (for more details see [23, 24]) In the case of ferromagnet, we have the following expression for the WZ term:
Now, one can extract the quantity which can be written in terms of the coherent spin fields obtained from difference of WZ terms at two neighbouring lattice sites. In the case of one-dimensional antiferromagnetic system, this turns out to be exactly the winding number in the long wavelength limit [22, 27–29]. For a unified description of Heisenberg ferromagnetic and antiferromagnetic spin systems in one and two spatial dimensions, we refer to . In more precise terms the difference between WZ terms from two neighbouring lattice sites is analogous to the surface density of the above winding number in terms of coherent spin field n expressed in a discretized form. where The term in the right-hand side of the first equation in (12) above stands for the “area of the cap” at the th vertex within a vortex of topological charge and is the number of such vortices. By topological charge , we mean that if we go around the vortex (antivortex) once, the spin at any vertex undergoes a rotation of ().
It should be kept in mind, however, that the above arguments involving the low temperature approximation to the nontopological terms do not interfere with our general formalism for studying the topological excitations, which is valid at all temperatures. The term represents the topological content in the WZ term in the case of ferromagnet, whereas in the case of antiferromagnet (9) represents the topological content in the entire WZ term. The existence of a nonvanishing topological term on the lattice implies that the entire spin configuration is fragmented into distinct topological sectors characterised by winding numbers or topological charges.
4. The Quantum Action
The quantum action corresponding to antiferromagnet is given by  where we have assumed the spin configurations to have correlations exhibiting a bipartite symmetry when the temperature is not too high .
The quantum action for ferromagnet is given by  The above actions can be derived by making use of the equations (2) to (6) (and chapter 5 in ). The integral containing the expression on the right-hand sides of both the actions corresponds to the constraint given by (5). The term in (14) acts as a topological charge measuring quantity for vortices (antivortices) which are constructed using usual coherent states for ferromagnets . The object in (13) acts as a topological charge measuring quantity for vortices (antivortices) which are constructed using staggered coherent states for antiferromagnets . One important result which came out from this investigation is that the vortices and antivortices in the case of Heisenberg ferromagnet behave as antivortices and vortices, respectively, in the case of Heisenberg antiferromagnet.
It can be shown that the above actions lead to existence of topological objects, even when with , and the WZ term represents the topological charge [23–25]. The topological objects are solely vortices and antivortices in the limit when the 3rd component of the coherent spin field is vanishingly small. The WZ term is evaluated on a vortex on the lattice by equations of motion and local periodic boundary conditions (see appendix in [23, 24]). For illustration, we present features of 2-vortices and 3-vortices to bring out essential properties of vortices and antivortices with or odd valued topological charges. These are essentially frozen (not mobile) vortices assuming that the system is below . For detailed descriptions, we refer to the above.
5. Analysis of 2-Vortex
For a typical 2-vortex, we refer to Figure 4. Here , , , and , each carry topological charge (for a 2-antivortex, we have the spins at the outer boundary oriented in the opposite direction and consequently each of the corresponding , , , and will carry a charge ). The horizontal and vertical spins are precisely the expectation values of the operators , in the coherent states at the lattice point (coherent state at a lattice point for in the present case), with the expectation being small.
However, the spin at the central lattice point of the vortex (the point O in Figure 4) becomes nonunique, as is clear from the construction of the subvortices. Thus, the central point situated on a lattice site turns out to be a defect or a singular point, and cancellation of the contribution of along the common bonds does not take place and therefore one does not have a consistent spin field configuration. The same statements hold for in the case of antiferromagnet, where the spin fields are now staggered spin fields and vortices (antivortices) of ferromagnet become antivortices (vortices). In other words, the contributions of the WZ term on the boundary of the 2-vortex are not well defined and the construction of the 2-vortex in this scheme becomes problematic. The scenario persists in all the vortices (antivortices) possessing even-valued topological charges.
6. Analysis of 3-Vortex
For 3-vortex, we have a consistent spin field configuration (Figure 5). We have now an elementary anti-vortex plaquette at the central region with well-defined usual or staggered spin field configurations, and the contributions of the WZ term (i.e., of or ) along the common bonds exactly cancel each other, giving rise to consistent spin field configurations, while along the boundary, it has a well defined finite value. It may be pointed out that in contrast to the case with even-valued charge, we now have a subvortex with opposite charge occupying the central region.
7. Phenomenology Based Analysis
So far, we have carried out the analysis, more of a mathematical nature, assuming the physical existence of topological excitations in these models. In order to substantiate the feasibility of this physical scenario, we now look at (i) the energetics and (ii) dynamics of these excitations.
Simple calculations based on the Hamiltonian show that the excitation energy of a 1-vortex, as measured from the ground state, is approximately both for ferromagnet and antiferromagnet [23, 24]. Please note that in our case. It is worthwhile to mention that the ground state spin configuration is obtained by taking the coherent spin vectors () aligned in the “same direction” (as appropriate to ferromagnet or antiferromagnet). Typically, these excitations are of much higher energy than the average thermal energy. Therefore, these topological excitations can be produced at very low temperature predominantly by means of quantum fluctuations, thereby maintaining a steady vortex antivortex pair density even at temperature close to zero.
The next important step in this direction would be to identify signature of vortex dynamics (a) in the theoretical calculation for the dynamical structure factor corresponding to our quantum spin models and (b) in the observed data from inelastic neutron scattering experiment corresponding to the appropriate materials. We have made a theoretical study analyzing available results from inelastic neutron scattering experiment performed on a quasi-two-dimensional spin- ferromagnetic material K2CuF4 and spin- antiferromagnetic material La2CuO4 [12, 17].
8. Concluding Remarks
(i)Our field theoretical approach will have to be extended to the case of finite , keeping in mind the relevant materials. Making use of the theoretical results one, can in principle evaluate the static and dynamic spin correlations for spin anisotropic quantum Heisenberg magnets at any temperature in two dimensions. These would contain contributions from both merons (antimerons) and magnons . The merons (antimerons) will emerge from the equations of motion satisfied by the quantum mechanical action. The spin-spin (in our case this spin is the coherent spin field, we have introduced) correlations would contain contributions from both merons (antimerons) and magnons. Some of the characteristics, distinguishing even and odd topological charge sectors, are likely to be reflected in our results for spin correlation functions in the case of nonvanishing .(ii)The aim of our analysis has been to formulate a consistent theory which can explain the occurrence of central peak in the neutron scattering experiments. The evaluation of the static and dynamic spin correlations will play an important role in the theoretical determination of the dynamical structure factor for low-dimensional quantum ferromagnets. The occurrences of “central peak” in the experimentally observed strongly indicate possible existence and dynamics of topological excitations of the kind described above. The “central peak” refers to the peak occurring at in the plot of versus in the constant -scan. Generally, this is an important signature for the spin dynamics driven also by the translational motion of the topological excitations and defects. These will play a crucial role in the quantitative verification of quantum BKT scenerio by making detailed comparison with the available experimental results.
The authors would like to acknowledge the help of Subhajit Sarkar in the preparation of the manuscript.
- J. M. Kosterlitz and D. J. Thouless, “Ordering, metastability and phase transitions in two-dimensional systems,” Journal of Physics C, vol. 6, no. 7, pp. 1181–1203, 1973.
- V. L. Berezinskii, “Destruction of long-range order in one-dimensional and two-dimensional systems having a continuous symmetry group—I. Classical Systems,” Journal of Experimental and Theoretical Physics, vol. 32, no. 3, pp. 493–500, 1971.
- V. L. Berezinskii, “Destruction of long-range order in one-dimensional and two-dimensional systems with a continuous symmetry group—II. Quantum systems,” Journal of Experimental and Theoretical Physics, vol. 34, no. 3, pp. 610–616, 1972.
- T. K. Ng, “Topological spin excitations of heisenberg antiferromagnets in two dimensions,” Physical Review Letters, vol. 82, pp. 3504–3507, 1999.
- A. A. Katanin and V. Y. Irkhin, “Magnetic order and spin fluctuations in low-dimensional insulating systems,” Physics-Uspekhi, vol. 50, no. 6, pp. 613–635, 2007.
- Y. il Shin, “Creation of 2D skyrmions in spin-1 polar Bose-Einstein condensates,” in Proceedings of the ULT, KAIST, Daejeon, Republic of Korea, 2011.
- C. Timm and K. H. Benneman, “Doping dependence of the Néel temperature in Mott-Hubbard antiferromagnets: effect of vortices,” Physical Review Letters, vol. 84, pp. 4994–4997, 2000.
- J. Des Cloizeaux and J. J. Pearson, “Spin-wave spectrum of the antiferromagnetic linear chain,” Physical Review, vol. 128, no. 5, pp. 2131–2135, 1962.
- H. J. Mikeska, “Solitons in a one-dimensional magnet with an easy plane,” Journal of Physics C, vol. 11, article L29, 1978.
- G. M. Wysin and A. R. Bishop, “Quantum Monte Carlo and transfer-matrix calculations for one-dimensional easy-plane ferromagnets,” Physical Review B, vol. 34, no. 5, pp. 3377–3387, 1986.
- C. Kawabata and A. R. Bishop, “Monte Carlo simulation of the two-dimensional classical Heisenberg model with easy-plane anisotropy,” Solid State Communications, vol. 42, no. 8, pp. 595–600, 1982.
- R. Chaudhury, “High temperature superconductivity—current status, our theoretical and experimental work,” Indian Journal of Physics A: Special Issue on High Temperature Superconductivity, vol. 66, pp. 159–180, 1992.
- K. Hirakawa, H. Yoshizawa, and K. Ubukoshi, “Neutron scattering study of the phase transition in two-dimensional planar ferromagnet K2CuF4,” Journal of the Physical Society of Japan, vol. 51, no. 7, pp. 2151–2158, 1982.
- K. Hirakawa, “Neutron scattering study of spin dynamics at the magnetic phase transition in two-dimensional planar ferromagnet K2CuF4,” Journal of the Physical Society of Japan, vol. 52, pp. 4220–4230, 1983.
- Y. Endoh, “Static and dynamic spin correlations in pure and doped La2CuO4,” Physical Review B, vol. 37, pp. 7443–7453, 1988.
- K. Yamada, “Spin dynamics in the two-dimensional quantum antiferromagnet La2CuO4,” Physical Review B, vol. 40, pp. 4557–4565, 1989.
- S. Sarkar, S. K. Paul, and R. Chaudhury, “Theoretical analysis of neutron scattering results for quasi-two dimensional ferromagnets,” European Physical Journal, vol. 85, no. 11, p. 380, 2012.
- E. Loh Jr., D. J. Scalapino, and P. M. Grant, “Monte Carlo studies of the quantum XY model in two dimensions,” Physical Review B, vol. 31, no. 7, pp. 4712–4714, 1985.
- R. H. Swendsen, “First- and second-order phase transitions in the d = 2 XY model,” Physical Review Letters, vol. 49, article 1302, 1982.
- D. D. Betts, F. C. Salevsky, and J. Rogiers, “Vortices in the two-dimensional s = 1/2 XY model,” Journal of Physics A, vol. 14, no. 2, pp. 531–538, 1981.
- R. Gupta, “Is Kosterlitz Thouless transition possible in 2D XY model?” Los Alamos National Laboratory preprint, 1985.
- E. Fradkin, Field Theories of Condensed Matter Systems, Addision-Wesley, San Francisco, Calif, USA, 1991.
- R. Chaudhury and S. K. Paul, “Physical realization and possible identification of topological excitations in quantum Heisenberg anti-ferromagnet on a two dimensional lattice,” European Physical Journal B, vol. 69, no. 4, p. 491, 2009.
- R. Chaudhury and S. K. Paul, “Physical realization of topological excitations in quantum Heisenberg ferromagnet on lattice,” European Physical Journal B, vol. 76, no. 4, p. 391, 2010.
- R. Chaudhury and S. K. Paul, “Possible existence of topological excitations in quantum spin models in low dimensions,” Physical Review B, vol. 60, pp. 6234–6237, 1999.
- R. Chaudhury and S. K. Paul, “Effective theory for a quantum spin system in low dimensions -beyond the long wavelength limit,” Modern Physics Letters B, vol. 16, p. 251, 2002.
- E. Fradkin and M. Stone, “Topological terms in one- and two-dimensional quantum Heisenberg antiferromagnets,” Physical Review B, vol. 38, pp. 7215–7218, 1988.
- F. D. M. Haldane, “O(3) Nonlinear σ model and the topological distinction between integer- and half-integer-spin antiferromagnets in two simensions,” Physical Review Letters, vol. 61, pp. 1029–1032, 1988.
- T. Senthil and M. P. A. Fisher, “Competing orders, nonlinear sigma models, and topological terms in quantum magnets,” Physical Review B, vol. 74, Article ID 064405, 11 pages, 2006.
- P. W. Anderson, S. John, G. Baskaran, B. Doucot, and S. D. Liang, Princeton University Preprint, 1988.
- A. M. M. Pruisken, “On localization in the theory of the quantized hall effect: a two-dimensional realization of the θ-vacuum,” Nuclear Physics B, vol. 235, pp. 277–298, 1984.
- R. B. Laughlin, “The relationship between high-temperature superconductivity and the fractional quantum hall effect,” Science, vol. 242, no. 4878, pp. 525–533, 1988.
- A. M. Finkelstein and A. I. Larkin, “Two coupled chains with Tomonaga-Luttinger interactions,” Physical Review B, vol. 47, pp. 10461–10473, 1993.
- R. Chaudhury and D. Gangopadhyay, “Interacting fermions, scaling and possible departure from fermi liquid behavior,” Modern Physics Letters B, vol. 9, no. 25, p. 1657, 1995.
Copyright © 2013 Ranjan Chaudhury and Samir K. Paul. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.