Open Access Article

Matter–wave interference of particles selected from a molecular library with masses exceeding 10[thin space (1/6-em)]000 amu

Sandra Eibenberger a, Stefan Gerlich a, Markus Arndt *a, Marcel Mayor *bc and Jens Tüxen b
aUniversity of Vienna, Faculty of Physics, VCQ, QuNaBioS, Boltzmanngasse 5, 1090 Vienna, Austria. E-mail: markus.arndt@univie.ac.at
bDepartment of Chemistry, University of Basel, St. Johannsring 19, 4056 Basel, Switzerland. E-mail: marcel.mayor@unibas.ch
cKarlsruhe Institute of Technology (KIT), Institute of Nanotechnology, P.O. Box 3640, 76021 Karlsruhe, Germany

Received 9th April 2013 , Accepted 5th July 2013

First published on 8th July 2013


Abstract

The quantum superposition principle, a key distinction between quantum physics and classical mechanics, is often perceived as a philosophical challenge to our concepts of reality, locality or space-time since it contrasts with our intuitive expectations with experimental observations on isolated quantum systems. While we are used to associating the notion of localization with massive bodies, quantum physics teaches us that every individual object is associated with a wave function that may eventually delocalize by far more than the body's own extension. Numerous experiments have verified this concept at the microscopic scale but intuition wavers when it comes to delocalization experiments with complex objects. While quantum science is the uncontested ideal of a physical theory, one may ask if the superposition principle can persist on all complexity scales. This motivates matter–wave diffraction and interference studies with large compounds in a three-grating interferometer configuration which also necessitates the preparation of high-mass nanoparticle beams at low velocities. Here we demonstrate how synthetic chemistry allows us to prepare libraries of fluorous porphyrins which can be tailored to exhibit high mass, good thermal stability and relatively low polarizability, which allows us to form slow thermal beams of these high-mass compounds, which can be detected using electron ionization mass spectrometry. We present successful superposition experiments with selected species from these molecular libraries in a quantum interferometer, which utilizes the diffraction of matter–waves at an optical phase grating. We observe high-contrast quantum fringe patterns of molecules exceeding a mass of 10[thin space (1/6-em)]000 amu and having 810 atoms in a single particle.


Introduction

Quantum physics has long been regarded as the science of ‘small things’, but experimental progress throughout the last two decades has led to the insight that it can also be applied to mesoscopic or even macroscopic objects. This applies for instance to the superposition of macroscopic numbers of electrons in superconducting quantum devices,1 the realization of large quantum degenerate atomic clouds in Bose–Einstein condensates,2 and the cooling of micromechanical oscillators to their mechanical ground state.3

Quantum superposition studies with complex molecules4 became possible with the advent of new matter–wave interferometers5–7 and techniques for slow macromolecular beams.8,9 These interferometers were, for instance, practically used to enable quantum enhanced measurements of internal molecular properties. The quantum fringe shift of a molecular interference pattern in the presence of external electric fields provided information, for example, on electric polarizabilities,10 dipole moments,11 or configuration changes.12 Molecule interferometry can complement mass spectrometry13 and help to distinguish constitutional isomers.14

In addition to their applications in chemistry, quantum interference experiments with massive molecules currently set the strongest boundaries on certain models that challenge the linearity of quantum mechanics.15

Further exploration of the frontiers of de Broglie coherence now profits from new capabilities in tailoring molecular properties to the needs of quantum optics. Our quantum experiment dictates the design of the molecules and the challenges increase with the number of atoms involved. In order to realize a molecular beam of sufficient intensity, the model compounds must be volatile, thermally stable and accessible in quantities of several hundred milligrams. Moreover, in order to minimize absorption at the wavelength of the optical diffraction (see below) they need to feature low absorption but sufficient polarizability at 532 nm. In response to these needs, the concept of a dendritic library is particularly appealing since it can be scaled up to complex particles, once a suitable candidate has been identified.

To meet these requirements we have functionalized organic chromophores with extended perfluoroalkyl chains. Such compounds show low inter-molecular binding and relatively high vapor pressures.16,17 They possess strong intra-molecular bonds and therefore sufficient thermal stability. In addition we start with a porphyrin core which is compatible with the required optical and electronic properties.18

In the past, monodisperse fluorous porphyrins were generated by substituting the four para-fluorine substituents of the tetrakis-pentafluorophenylporphyrin (TPPF20) by dendritically branched fluorous moieties.9 Using this approach, molecules composed of 430 atoms were successfully synthesized and applied in quantum interference experiments.19

With increasing complexity it becomes more challenging to purify monodisperse particles in sufficient amounts. Here we profit from the fact that our interferometer arrangement allows us to work with compound mixtures since each molecule interferes only with itself. By substituting some of the twenty fluorine atoms of TPPF20 with a branched, terminally perfluorinated alkylthiol (1), we obtain a mixture of compounds with molecular masses that differ exactly by an integer multiple of a particular value as a molecular library. The molecular beam density is sufficiently low for the molecules not to interact with each other and the individual library compounds can be mass-specifically detected using a quadrupole mass spectrometer (QMS Extrel, 16[thin space (1/6-em)]000 amu).

Our synthetic approach is based on the fact that penta-fluorophenyl moieties can be used to attach up to five polyfluoroalkyl substituents in nucleophilic aromatic substitution reactions. Substitutions at TPPF20 with its 20 potentially reactive fluorine substituents lead to a molecular library of derivatives with a varying number of fluorous side chains.

Results and discussion

We used sodium hydride as a base, microwave radiation as a heating source and diethylene glycol dimethyl ether (diglyme) as fluorophilic solvent. TPPF20 and a large excess of thiol 1 (60 equivalents) and sodium hydride in diglyme were heated in a sealed microwave vial to 220 °C for 5 minutes. After aqueous workup the resulting mixture was analyzed by MALDI-ToF mass spectrometry (Fig. 1) and subsequently used in our quantum interference experiments without further purification. We found up to 15 substituted fluorous thiol chains reaching a molecular weight well beyond 10[thin space (1/6-em)]000 amu.
Synthetic scheme and the MALDI-ToF mass spectrum of the fluorous porphyrin library L. High-mass matter–wave experiments were performed with component L12 of the library L. This structure is composed of 810 atoms and has a nominal molecular weight of 10 123 amu.
Fig. 1 Synthetic scheme and the MALDI-ToF mass spectrum of the fluorous porphyrin library L. High-mass matter–wave experiments were performed with component L12 of the library L. This structure is composed of 810 atoms and has a nominal molecular weight of 10[thin space (1/6-em)]123 amu.

In order to study the delocalized quantum wave nature of compounds in the fluorous library we utilize a Kapitza–Dirac–Talbot–Lau interferometer (KDTLI), which has already proven to be a viable tool with good mass scalability in earlier studies.6,19 The interferometer is sketched in Fig. 2: a molecular beam is created by thermal evaporation of the entire library in a Knudsen cell. The mixture traverses three gratings G1, G2, and G3, which all have the same period of d ≃ 266 nm. The molecules first pass the transmission grating G1, a SiNx mask with a slit opening of s ≈ 110 nm, where each molecule is spatially confined to impose the required spatial coherence by virtue of Heisenberg's uncertainty principle.20 This is sufficient for the emerging quantum wavelets to cover several nodes of the optical phase grating G2, 10.5 cm further downstream. The standing light wave G2 is produced by retro-reflection of a green laser (λL = 532 nm) at a plane mirror.


KDTL interferometer setup: the molecules are evaporated in a furnace. Three height delimiters, D1–D3, define the particle velocity by selecting a flight parabola in the gravitational field. The interferometer consists of three gratings with identical periods of d = 266 nm. G1 and G3 are SiNx gratings, whereas G2 is a standing light wave which is produced by retro-reflection of a green laser at a plane mirror. A phase modulation Φ ∝ (α·P)/(v·wy) is imprinted onto the molecular matter–wave via the optical dipole force which is exerted by the light grating onto the molecular optical polarizability αopt. Here P is the laser power, v the molecular velocity, wx ≃ 18 μm and wy ≃ 945 μm the Gaussian laser beam waists. The transmitted molecules are detected using electron ionization quadrupole mass spectrometry after their passage through G3, which can be shifted along the z-axis to sample the interference pattern.
Fig. 2 KDTL interferometer setup: the molecules are evaporated in a furnace. Three height delimiters, D1–D3, define the particle velocity by selecting a flight parabola in the gravitational field. The interferometer consists of three gratings with identical periods of d = 266 nm. G1 and G3 are SiNx gratings, whereas G2 is a standing light wave which is produced by retro-reflection of a green laser at a plane mirror. A phase modulation Φ ∝ (α·P)/(v·wy) is imprinted onto the molecular matter–wave via the optical dipole force which is exerted by the light grating onto the molecular optical polarizability αopt. Here P is the laser power, v the molecular velocity, wx ≃ 18 μm and wy ≃ 945 μm the Gaussian laser beam waists. The transmitted molecules are detected using electron ionization quadrupole mass spectrometry after their passage through G3, which can be shifted along the z-axis to sample the interference pattern.

When the molecular matter–wave traverses the standing light wave, the dipole interaction between the electric field of power P and the molecular optical polarizability αopt entails a periodic phase modulation Φ = Φ0·sin2(2πz/λL) with the maximum phase shift ugraphic, filename = c3cp51500a-t1.gif. Here vx is the forward directed molecular velocity, wy ≃ 945 μm is the Gaussian laser beam waist along the grating slits and z is the coordinate along the laser beam.

Interference of the molecular wavelets behind G2 leads to a molecular density pattern of the same period d in front of the third grating. As long as we can neglect photon absorption by the molecules in the standing light wave G2 acts effectively as a pure phase grating. The expected fringe visibility V is then given by21V = 2(sin(πf)/(πf))2J2(−sgn(Φ0[thin space (1/6-em)]sin(πL/LT))Φ0[thin space (1/6-em)]sin(πL/LT)). Here f designates the grating open fraction, i.e. the ratio between the open slit width s and the grating period d. J2 is the Bessel function of second order, L the separation of the gratings, LT = d2/λdB the Talbot length and λdB = h/mv the de Broglie wavelength with h as Planck's quantum of action, m the molecular mass and v the modulus of its velocity.

G3 is again a SiNx structure and lends spatial resolution to the detector. The interferogram is sampled by tracing the transmitted particle beam intensity as a function of the lateral (z) position of G3 and the mass selection is performed in the mass spectrometer behind this stage.

Talbot–Lau interferometers22 offer the important advantage over simple grating diffraction that the required grating period d only weakly depends on the molecular de Broglie wavelength: ugraphic, filename = c3cp51500a-t2.gif. The setup accepts a wide range of velocities and low initial spatial coherence. This facilitates the use of dilute thermal molecular beams. Diffraction at the standing light wave G2 avoids the dephasing caused by the van der Waals interaction between the molecules and a dielectric wall. This is indeed present in G1 and G3 but can be neglected there since the molecular momentum distribution at G1 is wider than that caused by the grating and any phase shift in G3 is irrelevant if we are only interested in counting the number of particles that reach the mass spectrometer.

Indistinguishability in all degrees of freedom is the basis for quantum interference23 and naturally given if a single molecule interferes only with itself. We only have to make sure that every molecule contributes to the final pattern in a similar way, which is true for all members of the library with about the same mass, independent of their internal state.

Differences between various molecules, such as their isotopic distribution or the addition of a single atom, are still acceptable. The KDTLI can tolerate a wavelength distribution of ΔλdB/λdB ≤ 20% and still produce a quantum fringe visibility in excess of the classical threshold.

We present here quantum interference collected at the mass of one specific library compound, particularly for L12 = C284H190F320N4S12 which has 12 fluorous side chains, a mass of 10[thin space (1/6-em)]123 amu and 810 atoms bound in a single hot nanoparticle.

All molecules of the library evaporated at a temperature of about 600 K. We selected the velocity class around v = 85 m s−1vFWHM = 30 m s−1) – corresponding to the most probable de Broglie wavelength of approximately 500 fm. This is about four orders of magnitude smaller than the diameter of each individual molecule. We detected the signal by electron ionization quadrupole mass spectrometry. During the interference measurements the mass filter was set to the target mass of L12 and only this compound contributed to the collected interference pattern.

The molecular beam was dilute enough to prevent classical interactions between any two molecules within the interferometer. Given that 80 mg of library L molecules were evaporated in 45 minutes, we estimate a flux at the source exit of ∼2 × 1015 particles per second. Including the acceptance angle of the instrument, the velocity selection as well as the grating transmission we estimate a molecular density inside the interferometer as 30 mm−3. This corresponds to a mean particle distance of about 300 μm which is sufficient to exclude interactions with other neutral molecules in the beam.

The average flight time of a molecule through the tightly focused standing light wave amounts to about 400 ns, i.e. much longer than the time scale of molecular vibrations (10−14–10−12 s) and rotations (∼10−10 s). Therefore, the mean scalar polarizability governs the interaction with the standing light wave although the optical polarizability is generally described by a tensor. Thermal averaging also occurs for the orientation of any possibly existing molecular electric dipole moment.11 The internal molecular states are decoupled from de Broglie interference as long as we exclude effects of collisional or thermal decoherence24,25 or external force fields.10

The thermal mixture of internal states is another reason why two-particle interference, i.e. mutual coherence of two macromolecules, is excluded in our experiments. The chances of finding two of them in the same indistinguishable set of all internal states – electronic, vibrational and rotational levels, configuration, orientation and spin – are vanishingly small.

In Fig. 3(a) we show a high contrast quantum interference pattern of L12. In contrast to far-field diffraction where the fringe separation is governed by the de Broglie wave-length of the transmitted molecules,26,27 near-field interferometry of the Talbot–Lau type generates fringes of a fixed period, which are determined by the experimental geometry. Specifically, the expected interference figure in our KDTLI configuration is a sine curve whose contrast varies with the phase-shifting laser power as well as with the molecular beam velocity and polarizability.


(a) Quantum interference pattern of L12 recorded at a laser power of P ≃ 1 W. The circles represent the experimental signal s as a function of the position z of the third grating. The solid line is a sinusoidal fit to the data, with a quantum fringe visibility of V = 33(2)%. The shaded area represents the background signal of the detector. A classical picture predicts a visibility of only 8% for the same experimental parameters. (b) Measured fringe visibility V as a function of the diffracting laser power P. The expected contrasts according to the quantum and the classical model are plotted as the blue and red lines, respectively.17 The dashed blue lines correspond to the expected quantum contrast when the mean velocity is increased (reduced) by 5 m s−1.
Fig. 3 (a) Quantum interference pattern of L12 recorded at a laser power of P ≃ 1 W. The circles represent the experimental signal s as a function of the position z of the third grating. The solid line is a sinusoidal fit to the data, with a quantum fringe visibility of V = 33(2)%. The shaded area represents the background signal of the detector. A classical picture predicts a visibility of only 8% for the same experimental parameters. (b) Measured fringe visibility V as a function of the diffracting laser power P. The expected contrasts according to the quantum and the classical model are plotted as the blue and red lines, respectively.17 The dashed blue lines correspond to the expected quantum contrast when the mean velocity is increased (reduced) by 5 m s−1.

We distinguish the genuine quantum interferogram5,21 from a classical shadow image by comparing the expected and experimental interference fringe visibility (contrast) with a classical model.

The far off-resonance optical polarizability is assumed to be well approximated by the static value αoptαstat ≃ 410 Å3 × 4πε0 as estimated using Gaussian G0928 with the 6-31 G basis set. The absorption cross-section of L12 at 532 nm was estimated using the value of the pure tetraphenylporphyrin dissolved in toluene29 assuming that the perfluoroalkyl chains contribute at least an order of magnitude less to this value.30 We thus find σ532 ≃ 1.7 × 10−21 m2.

In Fig. 3(b) we show the expected classical and quantum contrast as a function of the diffracting laser power. Our experimental contrast is derived from the recorded signal curves, such as that shown in Fig. 3(a), by V = (SmaxSmin)/(Smax + Smin), where Smax and Smin are the maxima and minima of the sine curve fitted to the data. The contrast is plotted as the black dots in Fig. 3(b). These points are in good agreement with the quantum prediction (blue line). The classical effect, describing the shadow image by two material gratings and one array of dipole force lenses,21 is shown as the red line.

The experiment clearly excludes this classical picture. Since the amount of precious molecular material and thermal degradation processes in the source limited the measurement time we performed our experiments only in the parameter regime that maximizes the fringe visibility, to show the clear difference between the quantum and classical contrasts. The experimental fringe visibility reproduces well the maximally expected quantum contrast.

Conclusion

We have shown that a library approach towards stable and volatile high-mass molecules can substantially extend the complexity range of molecular de Broglie coherence experiments to masses in excess of 10[thin space (1/6-em)]000 amu. Our data confirm the fully coherent quantum delocalization of single compounds composed of about 5000 protons, 5000 neutrons and 5000 electrons. The internal complexity, the number of vibrational modes and also the internal energy of each of these particles are higher than those in any other matter–wave experiment so far.

Acknowledgements

This work was supported by the FWF programs Wittgenstein (Z149-N16) and CoQuS (W1210-2), the European Commission in the project NANOQUESTFIT (304 886), the Swiss National Science Foundation, the NCCR ‘Nanoscale Science’, and the Swiss Nanoscience Institute (SNI). The authors thank Michel Rickhaus for the artwork of Fig. 2.

Notes and references

  1. J. E. Mooij, T. P. Orlando, L. Levitov, L. Tian, C. H. van der Wal and S. Lloyd, Science, 1999, 285, 1036–1039 CrossRef CAS.
  2. M. H. Anderson, J. R. Ensher, M. R. Matthews, C. E. Wieman and E. A. Cornell, Science, 1995, 269, 198–201 CAS.
  3. A. D. O'Connell, M. Hofheinz, M. Ansmann, R. C. Bialczak, M. Lenander, E. Lucero, M. Neeley, D. Sank, H. Wang, M. Weides, J. Wenner, J. M. Martinis and A. N. Cleland, Nature, 2010, 464, 697–703 CrossRef CAS.
  4. M. Arndt, O. Nairz, J. Vos-Andreae, C. Keller, G. van der Zouw and A. Zeilinger, Nature, 1999, 401, 680–682 CrossRef CAS.
  5. B. Brezger, M. Arndt and A. Zeilinger, J. Opt. B: Quantum Semiclassical Opt., 2003, 5, 82–89 CrossRef.
  6. S. Gerlich, L. Hackermüller, K. Hornberger, A. Stibor, H. Ulbricht, M. Gring, F. Goldfarb, T. Savas, M. Müri, M. Mayor and M. Arndt, Nat. Phys., 2007, 3, 711–715 CrossRef CAS.
  7. P. Haslinger, N. Dörre, P. Geyer, J. Rodewald, S. Nimmrichter and M. Arndt, Nat. Phys., 2013, 9, 144–148 CrossRef CAS.
  8. S. Deachapunya, P. J. Fagan, A. G. Major, E. Reiger, H. Ritsch, A. Stefanov, H. Ulbricht and M. Arndt, Eur. Phys. J. D, 2008, 46, 307–313 CrossRef CAS.
  9. J. Tüxen, S. Eibenberger, S. Gerlich, M. Arndt and M. Mayor, Eur. J. Org. Chem., 2011, 4823–4833 Search PubMed.
  10. M. Berninger, A. Stefanov, S. Deachapunya and M. Arndt, Phys. Rev. A: At., Mol., Opt. Phys., 2007, 76, 013607 CrossRef.
  11. S. Eibenberger, S. Gerlich, M. Arndt, J. Tüxen and M. Mayor, New J. Phys., 2011, 13, 043033 CrossRef.
  12. M. Gring, S. Gerlich, S. Eibenberger, S. Nimmrichter, T. Berrada, M. Arndt, H. Ulbricht, K. Hornberger, M. Müri, M. Mayor, M. Böckmann and N. Doltsinis, Phys. Rev. A: At., Mol., Opt. Phys., 2010, 81, 031604 CrossRef.
  13. S. Gerlich, M. Gring, H. Ulbricht, K. Hornberger, J. Tüxen, M. Mayor and M. Arndt, Angew. Chem., Int. Ed., 2008, 47, 6195–6198 CrossRef CAS.
  14. J. Tüxen, S. Gerlich, S. Eibenberger, M. Arndt and M. Mayor, Chem. Commun., 2010, 46, 4145–4147 RSC.
  15. A. Bassi, K. Lochan, S. Satin, T. P. Singh and H. Ulbricht, Rev. Mod. Phys., 2013, 85, 471–527 CrossRef.
  16. N. L. Stock, D. A. Ellis, L. Deleebeeck, D. C. G. Muir and S. A. Mabury, Environ. Sci. Technol., 2004, 38, 1693–1699 CrossRef CAS.
  17. P. J. Krusic, A. A. Marchione, F. Davidson, M. A. Kaiser, C.-P. C. Kao, R. E. Richardson, M. Botelho, R. L. Waterland and R. C. Buck, J. Phys. Chem. A, 2005, 109, 6232–6241 CrossRef CAS.
  18. M. Meot-Ner, J. H. Green and A. D. Adler, Ann. N. Y. Acad. Sci., 1973, 206, 641–648 CrossRef CAS.
  19. S. Gerlich, S. Eibenberger, M. Tomandl, S. Nimmrichter, K. Hornberger, P. J. Fagan, J. Tüxen, M. Mayor and M. Arndt, Nat. Commun., 2011, 2, 263 CrossRef.
  20. O. Nairz, M. Arndt and A. Zeilinger, Phys. Rev. A: At., Mol., Opt. Phys., 2002, 65, 0321091–0321094 CrossRef.
  21. K. Hornberger, S. Gerlich, H. Ulbricht, L. Hackermüller, S. Nimmrichter, I. V. Goldt, O. Boltalina and M. Arndt, New J. Phys., 2009, 11, 043032 CrossRef.
  22. J. Clauser and S. Li, Generalized Talbot–Lau Atom Interferometry, 1997 Search PubMed.
  23. P. A. M. Dirac, The principles of quantum mechanics, Clarendon Pr., Oxford, 1958 Search PubMed.
  24. L. Hackermüller, K. Hornberger, B. Brezger, A. Zeilinger and M. Arndt, Appl. Phys. B: Lasers Opt., 2003, 77, 781–787 CrossRef.
  25. L. Hackermüller, K. Hornberger, B. Brezger, A. Zeilinger and M. Arndt, Nature, 2004, 427, 711–714 CrossRef.
  26. O. Nairz, M. Arndt and A. Zeilinger, Am. J. Phys., 2003, 71, 319–325 CrossRef CAS.
  27. T. Juffmann, A. Milic, M. Müllneritsch, P. Asenbaum, A. Tsukernik, J. Tüxen, M. Mayor, O. Cheshnovsky and M. Arndt, Nat. Nanotechnol., 2012, 7, 297–300 CrossRef CAS.
  28. M. J. Frisch, et al., Gaussian 09, Revision A, Gaussian Inc., Wallingford, CT, 1st edn, 2009 Search PubMed.
  29. J. M. Dixon, M. Taniguchi and J. S. Lindsey, Photochem. Photobiol., 2005, 81, 212–213 CrossRef CAS.
  30. N. Gotsche, H. Ulbricht and M. Arndt, Laser Phys., 2007, 17, 583–589 CrossRef CAS.

Footnote

Synthetic protocol and analytical data of the porphyrin libraries L: General remarks: all commercially available starting materials were of reagent grade and used as received. Microwave reactions were carried out in an Initiator 8 (400 W) obtained from Biotage. Glass coated magnetic stirring bars were used during the reactions. The solvents for the extractions were of technical grade and distilled prior to use. Matrix Assisted Laser Desorption Ionization Time of Flight (MALDI-ToF) mass spectra were recorded on an Applied Bio Systems Voyager-DePro mass spectrometer or a Bruker microflex mass spectrometer. Significant signals are given in mass units per charge (m/z) and the relative intensities are given in brackets. Porphyrin library L: thiol 1 was synthesized in seven reaction steps in an overall yield of 70% as reported elsewhere.9 5,10,15,20-Tetrakis(pentafluoro-phenyl)-porphyrin (TPPF20, 4.0 mg, 4.10 μmol, 1.0 eq.), thiol 1 (193 mg, 246 μmol, 60 eq.) and sodium hydride (60% dispersion in mineral oil, 14.8 mg, 369 μmol, 90 eq.) were added to diglyme (4 mL) in a microwave vial. The sealed tube was heated under microwave irradiation to 220 °C for 5 minutes. After cooling to room temperature the reaction mixture was quenched with water and subsequently extracted with diethylether. The organic layer was washed with brine and water, dried over sodium sulfate and evaporated to dryness. The resulting product mixture (183 mg) was analyzed by MALDI-ToF mass spectrometry. MS (MALDI-ToF, m/z): 12[thin space (1/6-em)]403 (29%), 11[thin space (1/6-em)]645 (52%), 10[thin space (1/6-em)]884 (100%), 10[thin space (1/6-em)]121 (78%), 9339 (15%), 8597 (7%).

This journal is © the Owner Societies 2013