Share to: share facebook share twitter share wa share telegram print page

Neutron backscattering


Neutron backscattering spectroscopy (NBS) is a neutron scattering technique used to study atomic and molecular motion. It is distinguished by its ultra-high energy resolution, typically in the micro-electronvolt (μeV) range, which allows for the observation of processes on a nanosecond timescale.[1] The technique achieves this high resolution by using Bragg diffraction from crystals at a scattering angle very close to 180° (backscattering), a geometry where the wavelength spread of the diffracted neutrons is minimized.[2]

Introduction

Neutron backscattering spectroscopy is an inelastic neutron scattering technique that improves the energy resolution of neutron spectrometers by approximately two orders of magnitude compared to conventional methods.[1] The method exploits the reflection of neutrons from a near-perfect crystal back toward their source direction to achieve very high energy resolution for studying slow dynamical processes in condensed matter.[3] This technique is particularly valuable for investigating phenomena such as slow motions in complex liquids, jump diffusion, and quantum rotational tunneling.[2]

The accessible energy transfer range typically extends from -15 to +15 μeV, which corresponds to timescales from approximately 0.1 to 10 nanoseconds.[4] This makes backscattering spectroscopy uniquely suited for probing slow relaxation processes that are too fast for neutron spin echo but too slow for a conventional triple-axis spectrometer.[5]

Historical Background

The development of neutron backscattering spectroscopy emerged from the recognition that operating spectrometers at extreme scattering angles could significantly improve energy resolution. The theoretical foundation was established by Heinz Maier-Leibnitz and his colleagues in the 1960s. They demonstrated that the wavelength spread (Δλ) of a Bragg-diffracted neutron beam decreases as the scattering angle (2θ) approaches 180°.[6]

The first operational backscattering spectrometer was constructed at the Jülich Research Centre in 1966 by B. Alefeld and colleagues.[7] This prototype instrument demonstrated the feasibility of achieving μeV energy resolution but suffered from low neutron flux. Subsequent developments at high-flux reactor facilities such as the Institut Laue-Langevin (ILL) led to the construction of IN10 in 1972 and IN13 in 1978, which established backscattering as a routine technique for the scientific community.[8]

Fundamental Principles

Neutron Properties

Neutrons used in backscattering spectroscopy typically have wavelengths of 6–7 Å, corresponding to energies of approximately 2 meV.[2] The neutron's properties make it an ideal probe for this technique:

  • Zero electric charge, allowing deep penetration into materials without interacting with electron shells.
  • Magnetic moment, enabling studies of magnetic dynamics.
  • Mass comparable to atomic nuclei, facilitating efficient energy transfer during scattering events.
  • Wavelength comparable to interatomic distances, making it sensitive to atomic-scale structure and motion.[9]

The de Broglie wavelength (λ) of a neutron is related to its energy (E) by the equation:

where h is the Planck constant and mn is the neutron mass.[10]

Scattering Kinematics

The principle of backscattering relies on Bragg's law: λ = 2d sin θ, where d is the lattice spacing of the monochromator or analyzer crystal and θ is half the scattering angle.[11]

For inelastic scattering, the energy transfer (ℏω) and momentum transfer (Q) are related to the incident (ki) and final (kf) wavevectors:

In a backscattering geometry, the momentum transfer is maximized for elastic scattering, where Q = 2k.[10]

Elastic, Inelastic, and Quasielastic Scattering

The dynamic structure factor S(Q,ω) measured in backscattering experiments typically contains three components:[9]

  • Elastic scattering: A delta function, δ(ω), arising from the static (time-averaged) structure of the sample.
  • Quasielastic scattering: A broadened component centered at ω = 0 that results from slow, diffusive, or random motions.
  • Inelastic scattering: Discrete peaks at finite energy transfer (ω ≠ 0) that correspond to vibrational excitations or de-excitations (phonons, magnons, etc.).

The quasielastic component often follows a Lorentzian profile:

where Γ(Q) is the half-width at half-maximum (HWHM), which is related to the characteristic relaxation time (τ) of the motion by Γ = ℏ/τ.[9]

Energy and Momentum Resolution

The energy resolution in backscattering is determined by three main contributions:[12]

1. **Primary extinction**: Intrinsic width of the Bragg reflection from perfect crystals, with δE/E ≈ 10-6. 2. **Mosaic spread**: Spread in crystal plane orientations, where δE/E ≈ η·cot θ, with η being the mosaic width. 3. **Beam divergence**: Angular spread of the neutron beam, where δE/E ≈ α²/8 for a divergence angle α.

The total resolution is the quadratic sum of these contributions.[12]

Instrumentation & Techniques

Neutron Sources

Reactor Sources

Continuous neutron beams from research reactors provide a stable flux for backscattering measurements. Notable reactor-based instruments include:[13]

  • IN16B at ILL (France) - flux: 1.3×10⁴ n/cm²/s, resolution: 0.75 μeV
  • HFBS at NIST (USA) - flux: 1.0×10⁴ n/cm²/s, resolution: 0.8 μeV
  • SPHERES at FRM II (Germany) - flux: 8×10³ n/cm²/s, resolution: 0.65 μeV

Spallation Sources

Pulsed neutron beams from spallation sources enable time-of-flight methods for energy discrimination. Key instruments include:[4]

  • BASIS at SNS (USA) - flux: 2×10⁴ n/cm²/s, resolution: 3.5 μeV
  • IRIS at ISIS (UK) - flux: 5×10³ n/cm²/s, resolution: 1 μeV
  • DNA at J-PARC (Japan) - flux: 1×10⁴ n/cm²/s, resolution: 1.8 μeV

Backscattering Spectrometers

Fixed-Focus Instruments

Traditional backscattering spectrometers employ fixed monochromator and analyzer crystals arranged in an exact backscattering geometry.[13] The IN16B spectrometer at ILL exemplifies this design with:

  • Si(111) monochromator cooled to 200 K
  • A large, spherically arranged Si(111) analyzer array covering 2 steradians
  • A phase space transformation chopper for flux enhancement[14]

Time-of-Flight Backscattering

Modern instruments at spallation sources combine backscattering analyzers with time-of-flight techniques.[15] The BASIS spectrometer employs:

  • 84 Si(111) analyzer crystals at a scattering angle of 2θ = 176° (near-backscattering)
  • Bandwidth choppers for elastic line discrimination
  • A dynamic range of ±100 μeV with 3.5 μeV resolution[4]

Monochromators and Analyzers

Silicon crystals are predominantly used in backscattering applications due to their crystalline perfection and suitable lattice parameters.[6] Common crystal reflections include:

  • Si(111): d = 3.135 Å, λ = 6.27 Å, E = 2.08 meV
  • Si(311): d = 1.637 Å, λ = 3.27 Å, E = 7.64 meV
  • Ge(111): d = 3.266 Å, λ = 6.53 Å, E = 1.92 meV

The neutron energy can be varied by scanning the temperature of the monochromator, which changes the lattice spacing d through thermal expansion:[13]

where α ≈ 2.6×10-6 K-1 for silicon.

Sample Environments

Modern backscattering spectrometers can accommodate diverse sample environments to study materials under various conditions:[16]

  • Cryostats: Temperature range of 0.05–800 K
  • Pressure cells: Up to 10 kbar for studies of liquids and soft matter
  • Humidity chambers: 0–100% relative humidity control
  • Magnetic fields: Up to 14 T for quantum phenomena
  • Sample changers: Automated multi-position systems for high-throughput experiments

Data Acquisition & Reduction

Data Collection Modes

Backscattering experiments typically employ three measurement protocols:[2]

1. **Fixed window scans**: Energy-integrated intensity is measured as a function of temperature. 2. **Elastic scans**: Elastic intensity is measured as a function of Q at a fixed temperature. 3. **Inelastic scans**: The full dynamic structure factor, S(Q,ω), is measured at selected temperatures.

Time-of-flight discrimination is used at spallation sources to separate elastically scattered neutrons from the background.[4] The arrival time is:

where L₁ and L₂ are the primary and secondary flight paths, and v₀ is the incident neutron velocity.

Background Subtraction and Normalization

Standard data reduction procedures include:[17]

1. **Empty cell subtraction**: To remove scattering from the sample container. 2. **Vanadium normalization**: To correct for detector efficiency, using a sample of vanadium which scatters neutrons almost perfectly incoherently. 3. **Multiple scattering correction**: To account for neutrons that scatter more than once within the sample. 4. **Absorption correction**: To compensate for the attenuation of the neutron beam by the sample.

The corrected intensity is given by:

Resolution Function

The instrumental resolution function, R(Q,ω), is typically measured using a vanadium standard, which scatters purely elastically at room temperature.[17] The measured spectrum is a convolution of the true scattering function and the resolution function:

Deconvolution methods include Fourier deconvolution, maximum entropy methods, and Bayesian inference.[17]

Data Analysis & Modeling

Lineshape Fitting

The dynamic structure factor S(Q,ω) is analyzed by fitting it with model functions:[9]

  • **Lorentzian model** for simple diffusion:
  • **Stretched exponential** for complex relaxation, typically fitted in the time domain:

where β < 1 indicates a non-exponential relaxation process.[18]

Diffusion and Relaxation Models

Common models used to analyze the Q-dependence of quasielastic scattering include:[9]

  • **Jump diffusion (Chudley-Elliott model)**:

where Djump is the jump diffusion coefficient and τ₀ is the residence time between jumps.

  • **Continuous diffusion (Fick's law)**:

where D is the translational diffusion coefficient.

  • **Rotational diffusion**:

where jl are spherical Bessel functions and Drot is the rotational diffusion coefficient.[9]

Computational Tools and Software

Standard analysis packages for backscattering data include:[17]

  • DAVE (NIST): A comprehensive suite for data reduction and analysis.
  • LAMP (ILL): Large Array Manipulation Program.
  • Mantid (ISIS/SNS): A framework for neutron and muon data analysis.
  • QENS (JCNS): A specialized library for quasielastic neutron scattering analysis.

Molecular dynamics simulations provide complementary insight by calculating intermediate scattering functions that can be directly compared with experimental data:[19]

Applications

Condensed-Matter Physics

Quantum Tunneling

Backscattering spectroscopy has revealed quantum tunneling in numerous systems:[20]

  • Methyl group rotation in molecular crystals (tunnel splitting: 0.1–10 μeV)
  • NH₄⁺ rotation in ammonium salts (splitting: 1–100 μeV)
  • Hydrogen tunneling in metals (activation energy: 10–100 meV)

Magnetic Excitations

Studies of spin dynamics in frustrated magnets and quantum spin systems benefit from the μeV resolution:[21]

Soft Matter and Biology

Protein Dynamics

Backscattering reveals protein motions that are essential for biological function:[22]

  • The protein dynamical transition at ~200–230 K
  • Methyl group rotations (τ ~ 10-10 s)
  • Side chain fluctuations (τ ~ 10-9 s)
  • Domain motions (τ ~ 10-8 s)

The mean square displacement ⟨u²⟩, extracted from elastic scans, provides insight into protein flexibility:[23]

Polymer Dynamics

Backscattering probes polymer chain dynamics across multiple length and time scales:[24]

  • Segmental relaxation (τ ~ 10-10–10-8 s)
  • Side group rotation (τ ~ 10-11–10-9 s)
  • Methyl group dynamics (τ ~ 10-12–10-10 s)

Membranes and Lipids

Studies of lipid bilayers reveal:[25]

  • Lateral diffusion of lipids (D ~ 10-7 cm²/s)
  • Conformational dynamics of lipid tails (τ ~ 10-10 s)
  • Collective membrane undulations

Materials Science

Hydrogen in Metals

Backscattering is uniquely suited to probe hydrogen dynamics in metal hydrides:[26]

  • Jump diffusion between interstitial sites
  • Quantum tunneling at low temperatures
  • Trapping at defects and interfaces

Ionic Conductors

Studies of battery materials and solid electrolytes reveal:[27]

  • Li⁺ ion hopping mechanisms (τ ~ 10-10–10-8 s)
  • Correlation effects in ion transport
  • Coupling of framework dynamics to ion motion

Confined Systems

The dynamics of molecules in nanoporous materials can be investigated, including:[28]

Advantages and Limitations

Strengths

  • Ultra-high energy resolution: A full width at half maximum (FWHM) of 0.1–1 μeV enables observation of slow dynamics on nanosecond timescales.[1]
  • Bulk sensitivity: Neutrons penetrate deeply into materials, probing the entire sample volume.[29]
  • Isotope contrast: Substituting deuterium (D) for hydrogen (H) allows for selective labeling and contrast variation to highlight specific components.[9]
  • Simultaneous Q-range: Modern instruments with large detector arrays cover a wide range of momentum transfer (e.g., 0.2–2.0 Å⁻¹) simultaneously.[4]
  • Non-destructive: Low-energy neutrons typically do not cause radiation damage to the sample.[29]
  • Magnetic sensitivity: The neutron's magnetic moment directly probes spin dynamics.[10]

Challenges

  • Limited neutron flux: The high resolution is achieved at the cost of flux, leading to low count rates and long acquisition times (hours to days).[6]
  • Sample size: Relatively large sample quantities (0.1–10 g) are often required, depending on the scattering cross-section.[4]
  • Fixed geometry: The backscattering condition limits flexibility in accessing the full Q-ω space compared to other techniques.[1]
  • Limited energy range: The dynamic range is typically limited to ±15–100 μeV, which restricts studies of faster dynamics.[29]
  • Background signals: Weak signals can be dominated by instrumental background, requiring careful measurement and subtraction.[2]
  • Data analysis complexity: The need to deconvolve the measured spectrum from the instrumental resolution function complicates interpretation.[17]

Comparison with Other Spectroscopies

Neutron Spin Echo

Neutron spin echo (NSE) achieves even higher energy resolution (neV) but differs fundamentally:[29]

  • NSE directly measures the intermediate scattering function I(Q,t), whereas backscattering measures the dynamic structure factor S(Q,ω).
  • NSE is primarily sensitive to coherent scattering, while backscattering can measure both coherent and incoherent scattering.
  • NSE typically covers a lower Q-range (0.02–2.0 Å⁻¹) compared to backscattering (0.2–2.0 Å⁻¹).

Triple-Axis Spectroscopy

Triple-axis spectroscopy (TAS) offers greater flexibility but has lower resolution:[30]

  • Energy resolution is typically 10–100 μeV for TAS versus 0.1–1 μeV for backscattering.
  • TAS allows for flexible scans along any trajectory in Q-ω space, unlike the fixed geometry of backscattering.
  • TAS is better suited for studying dispersive excitations like phonons and magnons, while backscattering excels at quasielastic scattering.

Inelastic X-ray Scattering

Inelastic X-ray scattering (IXS) provides complementary information:[31]

  • IXS probes dynamics linked to electron density, whereas neutrons probe nuclear positions.
  • IXS is limited to coherent scattering, while neutron scattering can be dominated by incoherent scattering (especially from hydrogen).
  • IXS requires much smaller samples (μg) compared to neutrons (grams).

Raman and Brillouin Spectroscopy

Optical spectroscopies have much higher flux but different selection rules and sensitivities:[32]

  • They are limited to optically active modes, whereas neutrons are sensitive to all modes of motion.
  • Optical techniques are often more surface-sensitive, while neutrons provide bulk information.
  • The Q-range accessible with light scattering is very small (near the Brillouin zone center), whereas neutrons can cover the entire zone.

Safety Considerations

Neutron backscattering experiments are conducted at large-scale facilities and require adherence to strict radiation safety protocols:[33]

  • **Biological shielding**: Instruments are enclosed in concrete and polyethylene walls to minimize radiation dose rates in experimental halls.
  • **Access control**: Interlocked enclosures (hutches) prevent personnel from being exposed to the direct beam during operation.
  • **Sample activation**: Samples can become radioactive after exposure to the beam. This activation is monitored, and cooling periods may be required for materials containing certain elements.
  • **Dosimetry**: All personnel wear dosimeters to track their accumulated radiation dose, which is kept below regulatory limits (e.g., 20 mSv/year for occupational exposure).
  • **Training**: Users must complete radiation safety training and certification before being allowed to operate the instrument.

Next-Generation Spectrometers

The European Spallation Source (ESS) is expected to host advanced backscattering instruments with significantly higher performance:[34]

  • MIRACLES at ESS (under construction): This instrument aims for an energy resolution of 0.3 μeV, a flux gain of up to 30 times over existing instruments, a large dynamic range (±600 μeV), and a new analyzer geometry with 360° coverage.

Advanced concepts currently under development include polarized neutron backscattering for magnetic studies, multi-energy analyzers for simultaneous resolution options, and the use of machine learning for real-time data analysis.[35]

Integration with Complementary Techniques

A growing trend is the simultaneous or in situ measurement combining backscattering with other techniques:[29]

  • **Diffraction**: To establish direct structure-dynamics correlations.
  • **Small-angle neutron scattering (SANS)**: To probe multi-scale dynamics from nanometers to micrometers.
  • **Optical spectroscopy**: Combining neutron data with IR/Raman spectroscopy for a complete dynamical picture.
  • **Computation**: Validating and refining molecular dynamics simulations against high-resolution neutron data.

Emerging Applications

New frontiers for backscattering applications include:[36]

  • **Quantum materials**: Probing subtle excitations in topological insulators and unconventional superconductors.
  • **Energy materials**: Investigating ion transport in solid-state batteries and hydrogen dynamics in storage materials.
  • **Biological systems**: Studying dynamics in live cells, drug delivery mechanisms, and intrinsically disordered proteins.
  • **Extreme conditions**: Exploring dynamics at high pressures (>100 GPa) and high magnetic fields (>20 T).

See Also

References

  1. ^ a b c d "Back-scattering spectroscopy - Dynamics - Techniques". NMI3.
  2. ^ a b c d e A Meyer (2003). "The high-flux backscattering spectrometer at the NIST Center for Neutron Research". Review of Scientific Instruments. 74 (6): 2759–2777. doi:10.1063/1.1568557.
  3. ^ "Neutron Spectroscopy". Spectroscopy Online.
  4. ^ a b c d e f E Mamontov; K Herwig (2011). "A time-of-flight backscattering spectrometer at the Spallation Neutron Source, BASIS". Review of Scientific Instruments. 82: 085109. doi:10.1063/1.3626214.
  5. ^ CJ Carlile; MA Adams (1988). "The design of the IRIS inelastic neutron spectrometer". Physica B. 182: 431–440. doi:10.1016/0921-4526(92)90047-V.
  6. ^ a b c "REVIEW ON THE APPLICATIONS". ILL.
  7. ^ B Alefeld (1969). "A Neutron Backscattering Spectrometer with Very High Energy Resolution". Naturwissenschaften. 56: 410. doi:10.1007/BF00593617.
  8. ^ M Birr; A Heidemann; B Alefeld (1971). "A neutron crystal spectrometer with extremely high energy resolution". Nuclear Instruments and Methods. 95: 435–439. doi:10.1016/0029-554X(71)90542-8.
  9. ^ a b c d e f g M Bée (1988). Quasielastic Neutron Scattering. Adam Hilger. ISBN 0-85274-543-2. {{cite book}}: Check |isbn= value: checksum (help)
  10. ^ a b c GL Squires (1978). Introduction to the Theory of Thermal Neutron Scattering. Cambridge University Press. ISBN 0-521-21884-8. {{cite book}}: Check |isbn= value: checksum (help)
  11. ^ DA Neumann; JRD Copley; RL Cappelletti; WA Kamitakahara (1993). "Ultra-High Resolution Inelastic Neutron Scattering". Journal of Research of NIST. 98: 1–15. doi:10.6028/jres.098.001.
  12. ^ a b J Schelten; B Alefeld (1984). "Backscattering Spectroscopy". Proceedings of Workshop on Neutron Scattering Instrumentation: 378–398.
  13. ^ a b c B Frick; M Gonzalez (2001). "Five years operation of the second generation backscattering spectrometer IN16". Physica B. 301: 8–19. doi:10.1016/S0921-4526(01)00492-6.
  14. ^ F Demmel; A Fleischmann; W Häußler (2010). "The phase space transformation chopper for the ESS backscattering spectrometer MIRACLES". Nuclear Instruments and Methods A. 620: 524–529. doi:10.1016/j.nima.2010.03.150.
  15. ^ MTF Telling; KH Andersen (2005). "Spectroscopic characteristics of the OSIRIS near-backscattering crystal analyser spectrometer". Physical Chemistry Chemical Physics. 7: 1255–1261. doi:10.1039/B413934H.
  16. ^ E Mamontov (2016). "Fast macromolecular dynamics studied by inelastic neutron scattering". Biochimica et Biophysica Acta. 1861: 3513–3519. doi:10.1016/j.bbagen.2016.04.024.
  17. ^ a b c d e RT Azuah; LR Kneller; Y Qiu (2009). "DAVE: A comprehensive software suite for the reduction, visualization, and analysis of low energy neutron spectroscopic data". Journal of Research of NIST. 114: 341–358. doi:10.6028/jres.114.025.
  18. ^ A Arbe; J Colmenero (1998). "Characterization of the "simple liquid" state in a polymeric system". Physical Review Letters. 81: 590–593. doi:10.1103/PhysRevLett.81.590.
  19. ^ JH Roh; JE Curtis; S Azzam (2006). "Dynamics of hydration water in protein characterized by quasielastic neutron scattering". Biophysical Journal. 91: 2573–2588. doi:10.1529/biophysj.106.085688.
  20. ^ M Prager; A Heidemann (1998). "Rotational tunneling and neutron spectroscopy". Chemical Reviews. 97: 2933–2966. doi:10.1021/cr9500848.
  21. ^ B Lake; DA Tennant; CD Frost (2005). "Quantum criticality and universal scaling of a quantum antiferromagnet". Nature Materials. 4: 329–334. doi:10.1038/nmat1327.
  22. ^ W Doster (2010). "The protein-solvent glass transition". Biochimica et Biophysica Acta. 1804: 3–14. doi:10.1016/j.bbapap.2009.06.019.
  23. ^ G Zaccai (2000). "How soft is a protein?". Science. 288: 1604–1607. doi:10.1126/science.288.5471.1604.
  24. ^ D Richter; M Monkenbusch; A Arbe; J Colmenero (2005). "Neutron spin echo in polymer systems". Advances in Polymer Science. 174: 1–221. doi:10.1007/b106578.
  25. ^ MC Rheinstädter; T Seydel; F Demmel; T Salditt (2006). "Molecular motions in lipid bilayers". Physical Review E. 71: 061908. doi:10.1103/PhysRevE.71.061908.{{cite journal}}: CS1 maint: article number as page number (link)
  26. ^ M Karlsson (2015). "Proton dynamics in oxides". Dalton Transactions. 44: 10524–10538. doi:10.1039/C4DT03948C.
  27. ^ S Molto; M Kannan; JC Orain (2014). "Lithium diffusion in Li-ion battery materials". Journal of Power Sources. 245: 501–509. doi:10.1016/j.jpowsour.2013.06.153.
  28. ^ H Jobic; DI Kolokolov (2016). "Dynamics of ethane and propane in ZIF-8". Microporous and Mesoporous Materials. 225: 552–558. doi:10.1016/j.micromeso.2016.01.033.
  29. ^ a b c d e JS Gardner; G Ehlers; A Faraone; VM García Sakai (2020). "High-resolution neutron spectroscopy using backscattering and neutron spin-echo spectrometers". Nature Reviews Physics. 2: 103–116. doi:10.1038/s42254-019-0128-1.
  30. ^ G Shirane; SM Shapiro; JS Tranquada (2002). Neutron Scattering with a Triple-Axis Spectrometer. Cambridge University Press. ISBN 0-521-41126-8. {{cite book}}: Check |isbn= value: checksum (help)
  31. ^ AQR Baron (2020). "High-Resolution Inelastic X-Ray Scattering". Synchrotron Light Sources and Free-Electron Lasers: 1643–1719. doi:10.1007/978-3-030-23201-6_41.
  32. ^ M Cardona; G Güntherodt, eds. (1982). Light Scattering in Solids. Springer. ISBN 3-540-11513-7.
  33. ^ IAEA (2018). Radiation Protection and Safety of Radiation Sources. International Atomic Energy Agency. ISBN 978-92-0-102418-6. {{cite book}}: Check |isbn= value: checksum (help)
  34. ^ KH Andersen; DJ Argyriou; AJ Jackson (2020). "The instrument suite of the European Spallation Source". Nuclear Instruments and Methods A. 957: 163402. doi:10.1016/j.nima.2020.163402.{{cite journal}}: CS1 maint: article number as page number (link)
  35. ^ N Tsapatsaris; PK Willendrup; HN Bordallo (2016). "From BASIS to MIRACLES". EPJ Web of Conferences. 83: 03015. doi:10.1051/epjconf/20158303015.
  36. ^ VG Sakai; JW Arbe (2020). "Quasielastic neutron scattering in soft matter". Current Opinion in Colloid & Interface Science. 49: 94–106. doi:10.1016/j.cocis.2020.05.002.
Prefix: a b c d e f g h i j k l m n o p q r s t u v w x y z 0 1 2 3 4 5 6 7 8 9

Portal di Ensiklopedia Dunia

Kembali kehalaman sebelumnya