Open access

Classification of Mass-Produced Carbon Nanotubes and Their Physico-Chemical Properties

Written By

Heon Sang Lee

Submitted: 20 March 2012 Published: 09 May 2013

DOI: 10.5772/52613

From the Edited Volume

Syntheses and Applications of Carbon Nanotubes and Their Composites

Edited by Satoru Suzuki

Chapter metrics overview

2,946 Chapter Downloads

View Full Metrics

1. Introduction

Carbon nanotube (CNT) may be classified into single walled (SWCNT), double walled (DWCNT), and multiwalled carbon nanotube (MWCNT) according to the number of graphene layers. In some cases, bamboo-shaped multiwalled carbon nanotubes were also synthesized. Among these carbon nanotubes, multiwalled carbon nanotubes have been mass-produced in hundreds metric tons level. Many researchs on multiwalled carbon nanotubes point to an electrode, polymer composites, coating, and others. The number of graphene layers, purity, and crystallinity are the main features of multiwalled carbon nanotubes, which need to be characterized. We have proposed one important characteristic of multiwalled carbon nanotubes, the mesoscopic shape of MWCNT, of which many industrial applications may be comprised. According to our suggestion, one can determine the degree of tortuousness of MWCNT, quantitatively. (see ref.1-5 and sections 1-6 in this chapter)In this chapter, we will describe the mesoscopic shape factor of MWCNT in detail. Various physical properties as well as toxicity may strongly depend on the mesoscopic shape factor of MWCNT. Our suggestion has also been published as an international standard ISO/TS11888 by international organization for standardization (ISO) in 2011.

I hope readers enjoy the concepts and expressions shown in this chapter. Especially, this chapter shall be helpful to whom may want to develop a commercial application by selecting a proper CNT.

Advertisement

2. Static bending persistence length (SBPL, lsp)

If MWCNTs have no defect along their axis, their appearance would be straight to several hundred micro meter. Persistence length is the maximum straight length that is not bent by thermal energy. The persistence length of MWCNT is expected to be several hundred micro meter due to its exceptional high modulus. Static bending persistence length (SBPL) has been proposed in our earlier work to quantify the mesoscopic shape of MWCNT. SBPL is the maximum straight length that is not bent by permanent deformation. Fig. 1 shows the concept of SBPL. When a length considered is longer than SBPL, the shape of MWCNT looks tortuous. On the contrary, the shape of MWCNT looks straight as a length considered is shorter than SBPL.

Figure 1.

The concept of static bending persistence length of MWCNT (lsp).

If length considered is longer than SBPL, the shape of MWCNT looks tortuous. On the contrary, the shape of MWCNT looks straight as a length considered is shorter than SBPL.

Advertisement

3. Mathematical expression of SBPL (lsp)

The end-to-end vector can be obtained such as eq 1 when the distribution of bending points ({φ}(φ1,φ2,...,φk)) is given.

R=Ni=1kφiriE1

The spatial average of end-to-end distance R should be zero, since probability to bend to one direction is the same as that to the opposite direction. Then spatial average of square end-to-end vector is obtained as eq 2

R2=N2i=1kj=1k(φiri)(φjrj)=N2b2DbE2
DbR2/L2i=1kφi2E3

where Db is a bending ratio, φi=Ni/N , Ni is the number of unit segment in i-direction segment, N is the total number of unit segment, k=m+1, m is the number of static bending points on a coil, and ri is i-direction segment vector with the length of b. The expression shown in eq 3 is significant. This indicates that we can obtain the distribution function when we have enough data. This is often called as ill-posed problem. Regularization method in applied mathematics gives us the solution for solving the problems. Equation 3 holds only if a probability of the fold-back conformation is the same as that of the straight conformation.By using the scaling law, the coil expressed in eq 2 and 3 can be renormalized into the coil that has constant segment length, 2lp0. Then we can obtain eq 3 with φi=2lp0/L and k=L/2lp0. We can also consider a case where the bent angle (θ) between the ith and (i+1)th segments is a fixed small angle. The spatial average of the square end-to-end vector is obtained as following

rnrm=N(i=1kφi2)b2(cosθ)|nm|E4
R2=n=1km=1krnrm=n=1kp=n+1knrnrn+pn=1kp=rnrn+pE5
p=rnrn+p=N(i=1kφi2)b2(1+2p=1cospθ)=N(i=1kφi2)b2(1+cosθ1cosθ)E6
R2=(N2b2)(i=1kϕi2)(1+cos(θ)1cos(θ))=L2DbE7
DbR2L2(i=1kφi2)(1+cos(θ)1cos(θ))E8

Figure 2.

Tortuous MWCNT; bent points are distributed randomly along MWCNT axis.

Equation 7 can also be renormalized into the coil that has a constant segment length, 2lsp. The bending ratio (Db) is expressed as eq 9

DbR2L2(2lp0L)(1+cos(θ)1cos(θ))=C(2lp0L)=2lspLE9

where lsp=Clp0 is the static bending persistence length and Cshould be a constant for a fixed bent angle. The static bending persistence length is a statistical quantity, representing the maximum straight length that is not bent by static bending. In the case of continuous curvature, a more accurate statement is that the static bending persistence length is the mean radius of curvature of the rigid random-coil due to static bending. The same quantity arising from dynamic bending instead of static bending is dynamic bending persistence length (lp). The dynamic bending persistence length represents the stiffness of the molecules as determined by the effective bending modulus against thermal energy in Brownian motion. Equation 5 is valid when L>>lsp, the coil limit. Db=R2/L2=1 when L<lsp, the rod limit.If we know the values of end-to-end distance and contour length, the bending ratio can be obtained from the mean-squared end-to-end distance divided by the mean-squared contour length. The end-to-end distance of RRC varies with the change of bending angle. The difference can be compromised by using an arbitrary unit segment length which is similar to the scaling of polymer chain. The mean-squared end-to-end distance by the Kratky-Porod (KP) expression is given by eq 10 when the dynamic bending persistence length (lp) is replaced by the static bending persistence length (lsp) and the twice lsp equals to Kuhn length.

R2=2lspL+2lsp2(eL/lsp1)E10
Advertisement

4. Measurement methods for SBPL

The plot of eq 10 is presented in Fig. 3. Given data, the SBPL can be obtained by eq 11.

lsp=limL12dDbdlnLE11

Figure 3.

Bending ratio (Db) with respect to reciprocal contour length.

In this method, one need to have experimental data for R2 and L. In order to obtain these data, one have to cut MWCNTs into pieces with various L. Acid cutting or mechanical cutting method may be applied to obtain pieces of MWCNT. It is worth to note that R2 are Gaussian, given contour length (L). That is, various end-to-end distances may be measured for a constant L. This method is exact, but hard to obtain the experimental data.

Figure 4.

Approximation method to determine SBPL; the mean radius of curvature approximate SBPL.

The mean radius of curvature approximates the SBPL. One can easily obtain the mean value of the radius of curvatures of MWCNTs from any SEM images as seen in Fig. 4. The approximation method is convenient because SEM images of as-synthesized or as-received MWCNT can be directly used. The SBPL obtained by the approximation could have an error up to 200% compared to those obtained by exact method. However, the approximated value of SBPL still has physical significant in many applications, since many applied properties depend on the order of magnitude of SBPL.

Advertisement

5. Intrinsic viscosity of MWCNTs

From the molecular weight, the contour length, and the persistence length, the intrinsic viscosity of MWCNTs can be calculated. If we apply the intrinsic viscosity model of a worm-like coil to the rigid random-coil, the following expressions are obtained,

[η]=2.20x1021R23/2MfE12
f[1+0.926θ(Db)1/2]1E13
θ=ln(2lspe)2.431+(e/a)E14

where M is molecular weight, eis spacing between frictional elements along the contour, a=ζ/3πηs, ζis the friction factor for a single frictional element, and ηs is the solvent viscosity. For the non-draining limit for the random coil, f=1 , giving the maximum value of intrinsic viscosity in the model. When we take the static bending persistence length (lsp) as the length of a single frictional element, the friction factor of the element in eq 14 may follow the rigid-rod model such that ζT=3πηslsp/(ln(lsp/d)+0.3) for the translational motion and may be ζr=πηslsp3/(3(ln(lsp/d)0.8)) for the end-over-end rotational motion. Translational-rotational coupling and hydrodynamic shielding may also be considered for the evaluation of friction factor in eq 14.In this case, we can surmise that friction factor in eq 14 is scaled with lsps, where s is larger than unit value. We can reasonably neglect e/a in eq 14. The measurement of intrinsic viscosity assumes the deformation rate is slow enough. The intrinsic viscosity is determined by the competition of tendency of orientation toward flow direction and tendency to random orientation due to thermal motion (Brownian motion). The measurement often performed at shear rate of several hundreds reciprocal second. At this regime, CNTs may be extended to the static shape by shear force where peclet number (Pe=γ˙(2Rh)2/DT) is over 10.It is worth noting that the static bending persistence length determined from intrinsic viscosity is consistent with that determined from 3-D SEM analysis in dried state.

Advertisement

6. Diffusions of MWCNTs

Not only the toxicological issues but also researches on novel hybrid materials or nano-scale devices points to the need for the understanding of overall shape and mobility of carbon nanotube particles in a solution or in atmosphere.The degree of flexibility of carbon nanotubes is the major ingredient for the shape and mobility, however it is also puzzling.The persistence lengths of single-walled carbon nanotubes are expected to be in the order of tens to hundreds of micrometers due to their exceptionally large modulusand to have longer persistence lengths for muliwalled nanotubes, indicating currently prepared several-micrometer long nanotubes behave like rigid rods. Elastic fluctuations of semi-rigid particles by thermal energy have been described exactly by the worm-like coil model proposed more than 50 years ago by Kratky and Porod. The model describes the stiffness of molecules by dynamic bending persistence lengths (mean radius of curvatures) which are determined by effective bending modulus (Eeff) against thermal energy (kT) in a solution. Theoretical calculationshave shown that the dynamic bending persistence lengths (lps) of carbon nanotubes (CNTs) are up to several millimetersdue to their exceptionally large Young’s modulus of about 1.5TPa. Real-time visualization technique revealed that lp s of singlewalled carbon nanotubes (SWCNTs) are between 32 and 174μm, indicating SWCNTs shorter than lp(=32μm) may be rigid around room temperature in a solution. However, rippling developed on the compressive side of the tube leading to a remarkable reduction of the effective bending modulus, which is more pronounced for multiwalled carbon nanotubes (MWCNTs). Theoretical calculationshave shown that the effective bending moduli of MWCNTs are around 0.5 nNnm2 when the radii of curvatures are around 150 ~ 500nm. This indicates MWCNTs longer than 0.5 μmmight be flexible in a solution around room temperature, since thermal energy is about 4.1 x 10-3nNnm. It seems not likely that van der Waals interaction between graphene layers is the only reason that makes the effective bending moduli of MWCNTs more than 100 times smaller than SWCNT.

Both MWCNTs and SWCNTs discussed above are no more than worm-like coils (WLCs) where ensemble average of overall size (end-to-end distance) scales with the square root of molecular weight (contour length) in asymptotic limit. Our recent work has revealed that the spatial average of overall size of MWCNTs also follows the same scaling as WLCs in spite of their static bent points. We designated these MWCNTs as rigid random-coils (RRCs).The only difference between RRCs and WLCs is whether the bending points are static or dynamic by thermal energy.The relationship between the shape and size of RRCs has been characterized by static bending persistence lengths (lsps). Because both RRCs and WLCs are Gaussian, the models for the mobility of WLCs may also work for RRCs.

Translational diffusion coefficient is defined by the mobility of particle against thermal energy as Einstein relation, eq 15.

D=kTζE15

where k is Boltzman constant, T is temperature, and 1/ζ is the mobility. By analogy to macromolecule, a MWCNT with static bend points can also be considered to be made up N identical structural elements with a frictional factor ζe per unit element and a spacing e between elements along the contour of the coil.In this case, the mobility in eq 15 may be expressed as the sum of free-draining contribution (1/Nζe) and hydrodynamic interaction contribution which is called non-draining term.

1ζ=1Nζe+12!N2ij1ζij+E16

where ζij is the frictional factor by interaction between ith and jth element and ij.

Figure 5.

CNT made up N identical frictional element.

When we choose a spherical bead having diameter of a as a frictional element, the frictional factor of each element follows Stokes-Einstein relation, ζe=3aπηs where ηs is the viscosity of solvent. The frictional factor by interaction between ith and jth element may also follow Stokes-Einstein relation since a mean value of distance between elements i and jis small. Then, Kirkwood expression is obtained such as eq 17.

1ζ=1N3aπηs[1+a2Neij(1δij)rij1]E17
1R¯ij=1erij¯=(1e)ij(1δij)rij1E18

where rijRij/e and Rij is the distance between element i and j,e is solvated diameter of molecule, and N is number of frictional element, N=L/e. Here, we can see that mathematical expression for the mobility of RRCs is similar to that of WLCs. Equation 17 is widely used for the estimation of translational diffusion coefficient of macromolecules. Equation 18 has been solved by Hearst and Stockmayer using Riseman-Kirkwood theory and Daniel distribution. We notice that r¯ij is no more than a mean value of distance between the frictional element i and j. Then, r¯ij must depend on the conformation of the carbon nanotubes. Hearst and Stockmayer obtained r¯ij by the Kirkwood-Riseman theory as eq 19.

12xn(1δxn)rn1x=1Ndl1ldn0F(r,n)rdr=1Ndl1ldn0[F(r,n)f(r,n)]rdr+1Ndl1ldn0f(r,n)rdrE19

where F(r,n) is the unknown distribution for all n, f(r,n) is the known distribution, x is the contour distance of the point of interest from one end of the carbon nanotube, nis the contour distance from the point of interest to the frictional element n, and r is the displacement of frictional element n from the point of interest. Hearst and Stockmayer chose the Daniels distribution which includes a first-order correction to a Gaussian distribution as f(r,n), and obtained r¯ij as eq 20.

r¯ij=[2N(ln(2Lp)2.431+1.843(N/2Lp)1/2+0.138(N/2Lp)1/20.305(N/2Lp)1)]1E20

where N=L/e, L is contour length, eis spacing between frictional elements along the contour, Lp=lp/e, and lp is persistence length. The translational diffusion coefficient of worm-like coil can be estimated by eqs 15, 17, and 20. Rotational diffusion coefficientis expressed as eq 21.

Dr=(kTηs)(2lpL2)[0.253(L4lp)1/2+0.159ln(2Lp)0.387+0.160]E21

Equations 20 and 21 are valid for a semi-flexible rod when the contour length of rod is much longer than its persistence length such that the mean squared end-to-end distance follows random-coil scaling, R2=Nb2 where b=2lp.The semi-flexible rods in this coil limit, L>>lp , are so-called worm-like coils (WLCs). We see that the mobility is determined solely by the average conformation of particle with a given solvent viscosity and a contour length in eqs 13 and 14. We can reasonably surmise that the diffusion coefficients of RRCs are similar to those of WLCs with a given contour length, if the values of static bending persistence lengths (lsps) of RRCs are the same as those of the dynamic bending persistence lengths (lps) of WLCs.When hydrodynamic shielding effect is taken into account, the diffusion coefficients of RRCs might be slightly larger than those of WLCs due to the static bent points. The root mean-squared end-to-end distance of RRCs are given by eq 22.

R2=(N2b2)(i=1kφi2)(1+cos(θ)1cos(θ))=L2DbE22
DbR2L2(i=1kφi2)(1+cos(θ)1cos(θ))(2lp0L)(1+cos(θ)1cos(θ))=C(2lp0L)=2lspLE23

where Db is bending ratio, lsp is static bending persistence length, lp0 is an arbitrary constant segment length, θis static bent angle from the MWCNT axis, φi=Ni/N, Niis the number of unit segment in i-direction segment, N is the total number of unit segment, k=m+1, m is the number of static bending points on a coil. When RRCs have semi-flexibility by thermal energy, the ensemble average of bent angle (θ2) always becomes larger in amount of Δθ than the static bent angle θ;θ2=θ+Δθ. This is due to the fact that the effective bending modulus toward the bent direction is smaller than that toward the opposite direction.This indicates that the overall size of RRCs may be decreased when they are fluctuated by thermal energy. Because frictional elements of RRCs have Gaussian distribution by the definition of RRCs, those of semi-flexible RRCs also have Gaussian distrbution. Therefore, eqs 20 and 21 are also valid for semi-flexible RRCs when the persistence length is replaced by an apparent persistence length. The apparent persistence length (lap) is determined by the static bent angle (θ) and dynamic bent angle due to thermal energy (Δθ) as following.

lap=lsp(1+cos(θ+Δθ)1cos(θ+Δθ))(1cos(θ)1+cos(θ))E24

Expression for the translational diffusion coefficient of RRCs can be obtained from eqs 15, 17,18, and 24.

DT=kT3πηsL[1+ln(2Lap)2.431+1.843(N/2Lap)1/2+0.138(N/2Lap)1/20.305(N/2Lap)1]E25

Similary, expression for the rotational diffusion coefficient of RRCs can be obtained as following.

Dr=(kTηs)(2lapL2)[0.253(L4lap)1/2+0.159ln(2Lap)0.387+0.160]E26

where Lap=lap/e.Equations 25 and 26 are promising for the estimation of diffusion coefficients of MWCNTs synthesized by a CVD method. In other words, eqs 25 and 26 give us the information of the shape and size of MWCNTs if we have the measured values of diffusion coefficients.

Figure 6.

Size of carbon nanotube (R) is decreased to R2 at the elevated temperature.

Advertisement

7. Dynamic light scattering

The translational and rotational Brownian motions lead to fluctuation in the intensity of scattered light.The velocity of particles in Brownian motion can be directly measured by using dynamic light scattering (DLS) method, since time is correlated to obtain the intensity autocorrelation function (g(2)(t)). The intensity autocorrelation function is connected to the electric field autocorrelation function (g(1)(t)) which is given by eq 27 for a monodisperse solution.

|g(1)(t)|=exp(Γt)E27

where Γ=Dq2 with D, the translational diffusion coefficient of the molecules, and q, the scattering vector magnitude (q=4πnsin(θ/2)/λ0 where n is the solution refractive index, θ is the scattering angle, and λ0 is the incident light wavelength in vacuo). For polydisperse solutions, the electric field correlation function is given by a sum or distribution of exponentials,

|g(1)(t)|=0G(Γ)exp(Γt)dΓE28

where0G(Γ)dΓ=1

In many cases, a single exponent obtained from an average translational diffusion coefficient value fits the decay rate of the electric field correlation function such as eq 29

K1=Γav=0ΓG(Γ)dΓE29

where K1 is the first cumulant. When a solution is dilute enough neglecting the interaction between particles, the effective diffusion coefficient can be obtained by the intensity average. The intensity of light scattered by macromolecule species i is often proportional to the molecular weight (Mi) times the weight concentration (ci). In this case, the intensity average diffusion coefficient equals to z average diffusion coefficient.

The average decay rate (Γ) of the electric field autocorrelation function can be obtained by using conventional DLS.The first cumulant generally fits the data well for carbon nanotube solutions. When the incident light and detector are both vertical, Vv, translational diffusion are characterized by eq 30.

ΓVv=q2DTE30

When the incident light vertical and detector horizontal, Hv , the diffusion of anisotropic particle are characterized by eq 31

ΓHv=q2DT+6DRE31

where DTis translational diffusion coefficient and DR is rotational diffusion coefficient. This equation is valid if the particle rotates many times while diffusing a distance comparable to q1 or if there is little anisotropy in particle dimension. MWCNTs solution meets with the former case in this work. Now we have three independent mathematical model eqs 25, 26, and 12. And two equations for DLS measurement.Three unknown shape factors of static bending persistence length, contour length, and thermal fluctuation angle can be determined from the measured diffusion coefficients and intrinsic viscosity using eqs 25, 26, and 12.This is uniqueness of carbon nanotubes compared to macromolecules, since macromolecules have only two unknown shape factors of persistence length and contour length.

Advertisement

8. Micro rheology

The terminology of “microrheology” is used, to distinguish the technique from conventional (macro) rheology. In the microrheology, colloidal particles are used for probing the rheology of material of interest. The starting point is the Stokes-Einstein relation.

DT=kT6aπηsE32

If we have measured values of translational diffusion coefficient, the viscosity of material of interest can be easily obtained by

ηs=kT6aπDTE33

where a is the radius of spherical colloidal particle; colloidal particles have usually average diameter between 1 nm and 1000 nm. It is comparable the ISO definition of nanoparticles those have average diameter between 1 nm and 100 nm. In this sense, nanoparticles are just some kinds of colloidal particles. When we have measured value of mean-squared displacement of probe particle, the following equation can be utilized.

r2(t)=6DTt=kTπaηsE34

This seemingly simple idea has done a great impact on various research fields, indeed. One example is the nanoparticles dispersed in a polymer melt. It is often reported that nanoparticles seems diffuse faster than expected. The origin of this phenomenon lies in the “Nano” size. The viscosity of polymer melt is well described by integral constitutive equations such as reptation model. In this model, the viscosity is determined by the stress relaxation time of polymer chain from the constraint of entanglement. When the observation time is much shorter than any relaxation time of polymer in rheometry of frequency sweep, the polymers behave like a crosslinked rubber, exhibiting a plateau modulus. The plateau modulus of polymers is determined from the entanglement lengths of polymer such as

GoN=ρRTMeE35

The plateau modulus of polymer is usually reported in the order of 10 6~7Pa. Entanglement molecular weight of polymer is about 1000~2000 g/mole. The entanglement length is about 10~100 nm.The particles having comparable size to the entanglement length of a polymer would feel less frictional force than expected from the melt viscosity in macrorheology. Therefore, viscosity of polymer melt is much lower for the nanoparticles. This may lead to the faster thermal motion of nanoparticle compared to a larger particles.

Advertisement

9. Applications

In our previous works, we demonstrated that MWCNT having shorter SBPL have a certain merit in a polymer composite for electrical conductive application. When MWCNTs are needle-like, polymer composites comprised of them exhibit higher electrical conductivity compared to those comprised of tortuous MWCNTs. The situation changed drastically when the composites were molded into a specimen by injection molding machine. The needle like CNTs aligned to the flow direction, which broke the electrical conductive networks, then the composites lose the electrical conductivity. However, this problem was not observed when the composites contained tortuous MWCNTs which have a short SBPL. We also showed that the electrical percolation threshold depends on the length of MWCNT when MWCNT are needle-like. But, the electrical percolation threshold depends on the SBPL for a tortuous MWCNT. Thermal conductivity and linear thermal expansivity are also strongly dependant properties on the SBPL of MWCNT. Especially, thermal shrinkable material can be fabricated as well as thermal expansive material by controlling SBPL of MWCNT.

References

  1. 1. Lee, H. S.; Yun, C. H.; Kim, S. K.; Choi, J. H.; Lee, C. J.; Jin, H. J.; Lee, H.; Park, S. J.; Park, M. Appl. Phys. Lett. 2009, 95, 134104.
  2. 2. Lee, H. S.; Yun, C. H. J. Phys. Chem. C 2008, 112, 10653-10658.
  3. 3. Lee, H. S.; Yun, C. H.; Kim, H. M.; Lee, C. J. J. Phys. Chem. C 2007, 111, 18882- 18887.
  4. 4. Han, M. S.; Lee, Y. K.; Yun, C. H.; Lee, H. S.; Lee, C. J.; Kim, W. N. Synthetic Metals 2011, 161, 1629-1634.
  5. 5. Heo, Y. J.; Yun, C. H.; Kim, W. N.; Lee, H. S. Curr. Appl. Phys. 2011, 11, 1144-1148.
  6. 6. Hearst, J. E.; Stockmayer, W. H. J. Chem. Phys. 1962, 37, 1425-1433.
  7. 7. Squires, T. M.; Mason, T. G. Ann. Rev. Fluid Mech. 2010, 42, 413-438.
  8. 8. Hearst, J. E. J. Chem. Phys. 1963, 38, 1062-1065.
  9. 9. Hearst, J. E. J. Chem. Phys. 1964, 40, 1506-1509.
  10. 10. Kratky, O; Porod, G. Recl. Trav. Chim. Pays-Bas 1949, 68, 1106-1122.
  11. 11. Doi, M.; Edwards, S. F. The theory of polymer dynamics; Oxford University Press: New York, U.S., 1986.
  12. 12. Berne, B. &Pecora, R. Dynamic Light Scattering; Wiley: New York, 1976.
  13. 13. Tanford, C. Physical Chemistry of Macromolecules; John Wiley & Sons: New York, 1961.
  14. 14. Koppel, D. E. J. Chem. Phys. 1972, 57, 4814-4820.
  15. 15. Cush, R.; Russo, P. S.; Kucukyavuz, Z.; Bu, Z.; Neau, D.; Shih, D.; Kucukyavuz, S.; Ricks, H. Macromolecules1997, 30, 4920-4926.

Written By

Heon Sang Lee

Submitted: 20 March 2012 Published: 09 May 2013