Monthly Archives: August 2026

Preliminary Raman Evidence for Boron-Rich Al–B–C Phases in Natural Topaz from Schneckenstein/Vogtland, Germany

DOI: 10.31038/GEMS.2026861

Abstract

Natural topaz crystals from the Schneckenstein deposit, Vogtland, Germany, were examined by optical microscopy and low-power Raman spectroscopy to characterize rare colorless spherical inclusions associated with previously reported diamond-bearing topaz. The inclusions, commonly up to about 25 µm in diameter and including smaller examples near 5 µm, occur within the topaz matrix and are therefore unlikely to represent surface contamination. Raman spectra obtained at low laser power show a characteristic band near 471 cm-1, together with additional bands that are broadly consistent with boron-rich phases or boron-bearing Al–B–C microdomains. Spectral similarities are noted with rhombohedral boron, AlB12-type compounds, and boron carbide components, but these assignments remain tentative because several possible phases have overlapping Raman features. Bands near 1331–1332 cm-1 and around 1450 cm-1 may indicate boron-doped diamond or related carbon-bearing phases, although alternative interpretations cannot yet be excluded. The unexpected optical transparency of the inclusions, possible pressure-related band shifts, and the use of Raman spectroscopy as the principal method make the phase identification preliminary. Overall, the observations provide preliminary Raman evidence for unusual boron-rich Al–B–C microphases in Schneckenstein topaz and may be compatible with high-pressure processes and transport by supercritical fluids or melts during Variscan mineralization. Confirmation by complementary structural and chemical methods remains essential

Keywords

Topaz, Schneckenstein, Vogtland, Raman spectroscopy, Boron-rich inclusions, Rhombohedral boron, AlB12, Boron carbide, Boron-doped diamond, High-pressure mineralization, Supercritical fluids.

Summary

In summary, the Schneckenstein topaz contains rare, water-clear spherical inclusions whose Raman spectra indicate a complex association of boron-rich phases, possible AlB12-type components, boron carbide, and boron-doped diamond. Their occurrence within unprepared natural topaz excludes preparation-related contamination and points to an unusually high-pressure mineral assemblage preserved in the topaz matrix. Although the phase assignments remain provisional, the observations strengthen the interpretation that supercritical fluids or melts played an important role in transporting mantle-derived or ultrahigh-pressure components into the Variscan crustal mineralization system.

Introduction

The Schneckenstein deposit in Vogtland, Germany, is a classical occurrence of natural topaz and has long attracted mineralogical interest because of the quality, transparency, and unusual inclusion inventory of its crystals [1,2]. Earlier work on Schneckenstein topaz included determinations of the fluorine content and studies of fluid inclusions, which revealed boric-acid-rich fluids and provided evidence for complex postmagmatic processes during Variscan mineralization [3]. Recent Raman spectroscopic investigations have expanded this picture by reporting micrometer-sized spherical diamond crystals in Schneckenstein topaz (see Table 1 ), together with additional high-pressure or high-temperature mineral relics [4]. Ti-rich topaz relics with the idealized composition Al2(TiO4)F2 were also described and interpreted as possible indicators of unusual formation conditions. These observations suggest that Schneckenstein topaz may preserve a more complex inclusion assemblage than previously recognized, although the pressure–temperature significance of individual phases remains to be tested carefully. During further microscopic screening of natural, unprepared topaz cleavages, rare water-clear spherical inclusions were observed that differ from diamond in their Raman spectra. The present contribution focuses on these inclusions. Using optical microscopy and low-power Raman spectroscopy, the study documents their occurrence, compares their spectra with reference data for boron-rich phases, and evaluates whether they may represent unusual boron-bearing Al–B–C microphases preserved in the topaz matrix.

Table 1: Results of the Raman measurement in the first-order range of diamond from Schneckenstein [4].

Diamond

FWHM G-band FWHM n

Topaz from Schneckenstein

1326.6 ± 4.2 cm-1

41.6 ± 25.4 cm-1 1597.3 ± 2.8 cm-1 61.2 ± 3.5 cm-1 14
1304.2 ± 1.9 cm-1 22.7 ± 1.7 cm-1 not present –

7

FWHM – Full Width at Half Maximum. n – number of measured crystals.

Sample Material

The Schneckenstein topaz samples were collected during a visit to Schneckenstein in October 1960, together with Peter Haupt. In this study, only centimeter-sized natural crystal cleavages, generally perpendicular to the c-axis, were used (Figure 1). Descriptions of the Schneckenstein occurrence are given by Rösler et al. (1968) [1] and Lahl (2012) [2]. The general appearance of the free topaz crystals is shown in Figure 1.

Figure 1: Topaz crystals from Schneckenstein. The Petri dish has a lower diameter of 9 cm. Most specimens have crystal cleavages perpendicular to the c-axis.

A typical Raman spectrum of topaz from Schneckenstein is shown in Figure 2. According to Thomas (2026a [6], 2026b [4], the topaz from Schneckenstein contains 19.26 ± 0.74.% F; earlier measurements yielded 19.15 ± 0.14 % F (Thomas 1982) [3]. The notation % F means [% (g/g)], because wt.% is not a correct and allowed unit (IUPAC 1988), see Ebel et al. (2004) [5]. The Schneckenstein samples are ideal for optical measurements and the study of diamonds, as no preparation is required. Therefore, we can exclude contamination from diamond-bearing tools.

Figure 2: Typical Raman spectrum of topaz from Schneckenstein (532 nm laser, 0.9 mW on sample).

The topaz crystals are generally wine-yellow and water-clear. Only a few crystals contain visible black spots, mostly composed of rutile. The topaz is also rich in fluid inclusions, with homogenization temperatures of 407 ± 25°C (n = 54) and critical behavior [3]. The fluid inclusions are rich in boric acid and are arranged along planes perpendicular to the c-axis. Based on their overall character, they are interpreted as secondary inclusions.

Figure 2 shows a typical Raman spectrum of topaz from Schneckenstein, obtained with the 532 nm laser and a typical sample power of 0.9 mW.

A strong triplet characterizes the Raman spectrum of topaz at approximately 239, 267, and 285 cm-1, and by a medium-intensity band at 923 cm-1. This topaz triplet overlaps with the low-frequency Raman range studied for the inclusions.

Microscopy and Raman Spectroscopy: Methodology

Preliminary studies of the samples were performed with a Zeiss JENALAP pol equipped with a universal stage UT 124 and various compensators (Brace-Köhler and Ehringhaus). We performed all microscopic and Raman spectroscopic studies with a petrographic polarization microscope (BX 43) with a rotating stage coupled with the EnSpectr Raman spectrometer R532 (Enhanced Spectrometry, Inc., Mountain View, CA, USA) in reflection and transmission. We used Raman spectroscopy as the main tool for characterizing the inclusion crystals (diamond, DLC, graphite, carbon), as well as the newly identified boron-bearing phases. The Raman data presented here were acquired using an EnSpectr Raman microscope (RamMics R532). Measurements covered the spectral range from 0 to 4000 cm⁻¹ and were conducted with a single-mode 532 nm laser operating at a maximum output of 50 mW. Analytical settings included a 20 µm entrance aperture, a holographic grating of 1800 g mm⁻¹, and a spectral resolution of approximately 4 cm⁻¹ (generally for the whole range). For most analyses, a long-working-distance Olympus LMPlanFL 100× objective is used. Laser power at the sample surface was continuously adjustable down to 0.02 mW. Higher powers (up to 50 mW) were applied only for overview measurements. Generally, for all Raman measurements during this study, we used 0.9 mW of the 532nm laser on the sample over the whole range (0-4000 cm-1) to prevent local heating, which is particularly critical for boron solids. Raman band positions were calibrated before and after each measurement series using the Si band of a semiconductor-grade single-crystal silicon chip. Based on 20 repeated measurements, run-to-run reproducibility was ± 0.2 cm⁻¹ for silicon (520.2 ± 0.2 cm⁻¹) in the measuring range 100 to 1450 cm-1. As a second reference, we used a water-clear diamond crystal from Brazil. For the first order diamond line, we obtained (1330.1 ± 0.6) cm-1 with a FWHM = 5.1 ± 0.1 cm-1. FWHM stands for Full Width at Half Maximum.

Attention: High Risk of Thermal Phase Alteration

The green lasers possess high energy density and are highly absorbed by dark, ceramic phases like B4C and α-AlB12 or sp2-rich

DLC:

  • Laser-Induced Oxidation: If the laser power is set too high (typically above 5–10 mW at the sample surface), the localized focal spot can rapidly heat up in That converts the sample into B2O3or α-Al2O3 right under the beam.
  • Spectral evidence of burning: The sudden, irreversible emergence of sharp peaks at 378 and 418 cm-1(α-Al2O3) or the widening of the boron peaks mid-scan indicate that the laser power is too high.
  • Fix: Drop the laser power down to 1–2 mW and use longer acquisition times to collect a clean spectrum Therefore, we generally use a laser power of 0.9 mW on the sample.

For example, the fingerprint Raman line at 471.5 cm-1 (0.9 mW) shifts to 477.6 cm-1 at 30 mW and to 477.8 cm-1 at 50 mW; at 50 mW, the characteristic lines of α-Al2O3 at 377.4 and 419.5 cm-1 appear.

Results

Careful screening of the topaz crystals revealed, in addition to spherical diamond crystals, colorless spherical crystals that differ strongly from diamond by Raman spectroscopy. Figure 3 shows a typical crystal. Such large inclusions, about 25 µm in diameter, are rare; smaller ones are relatively widespread. The inset in Figure 3 shows a small inclusion (5 µm) not much wider than the large one, with the same Raman characteristics as the larger one.

Figure 3: Inclusion in topaz. The inclusion is covered by a thin topaz (Toz) layer (36 µm thick). The inserted figure on the left shows the same inclusion, and the smaller one (13 µm depth) on the top left is more typical.

The depth in the topaz matrix shows clearly that the spherical crystals are not surface contamination.

Figure 4 shows the Raman spectrum of an inclusion part not covered with a topaz layer (this inclusion lies directly on the surface.

Figure 4: Overall view of the Raman spectrum of the large inclusion, similar to the crystal-like Figure 3. The Raman band between 1300 and 1500 cm-1, according to Zaitsev (2001) [9], is observed in boron-doped diamond films around the inclusion.

The Raman spectrum has, in addition to the broad and intense band under 200 cm-1 (of unclear origin), an intense Raman peak at 471.5 cm-1, which may serve as a fingerprint of β-rhobohedral boron [7]. Raman spectra from the inclusion in natural topaz from Schneckenstein show bands compatible with boron-rich Al–B or Al–B–C structures. Although a direct assignment to α-AlB12 or γ-AlB12 is not yet possible, the spectra suggest the presence of unusual boron-bearing inclusions or microphases that warrant further structural and chemical investigation. Table 2 (first column) shows the measured results for 0.9 mW laser power on the sample, compared with the data presented by Werheit et al. (2010) [7]. In Table 2, the Raman band corresponding to boron-doped diamond, indicated by the strong band at 1452.8 ± 2.5 cm-1, is not shown. The diamond line is from 10 different inclusions at 1331.2 ± 0.6 cm-1, the FWHM = 26.3 ± 4.2 cm-1. The Raman spectrum of the spherical crystal 5 is shown in Figure 5. A strong Raman fingerprint band at 471.6 cm-1 also characterizes C-doped β-rhombohedral boron. The bands at 261, 848, and 925 cm-1 belong to the topaz matrix.

Table 2: Results of the Raman measurements on 10 different spherical crystals in topaz (this work) and comparison with data from Werheit et al. (2010) [7] and Jay et al (2023) [8] for B4.3C. Data in cm-1. Laser: 532 nm, 0.9 mW on sample.

Thomas

Werheit et al. (2010) [7] and Jay et al. (2023) [8]

First-order Raman lines

α-rh. boron β-rh. boron a-AlB12 π-AlB12 B4.3C
    282    

270

297.4 ± 9.6

  309 294 294  
         

321

351.0 ± 1.6

  357 357 346  
   

376

     

  391      
431.7 ± 1.1   414 416 414

417

  456      
471.1 ± 0.9   480 472 476

477, 488

493.1 ± 9.8

494   495 495  
  527   516 516

529, 536

570.3 ± 3.2

552 565 561 562  
  589 594   578

575

607.9 ± 2.9

    615 615  
    630 632

633

 
 

    656 652
  694 685   680

699

714.2 ± 3.8

713   711 709 719
762.1 ± 1.8 750

773

     
 

778

      768
  795     805

815

   

813

  844 841
854.0 ± 4.4 873 885 866

867

 
     

905  
925.0 ± 3.6 934   927

932

 
   

987

    994
1044.3 ± 2.5     1039 1035

1042

1075.3 ± 4.6

  1097 1058 1058 1071
  1094      

1096

1114.9 ± 9.7

1125   1125 1125 1131
 

1160

       
 

1187

       
  1201   1219

1223

 

1251.1 ± 2.6

1238 1217 1248

1252

 

Second-order Raman lines

1452.8 ± 2.5

1464

       

1719

1710

       

a-rh. boron – a-rhombohedral boron, b-rh. boron – b-rhombohedral boron, α-AlB12 and π-AlB12 – two different crystal modifications of aluminumdodecaboride (AlB12), and B4.3C – boron carbide.

Figure 5: Raman spectrum of spherical crystal 5. The bands around 261 cm-1 are due to the matrix topaz. The strong band at 471.6 cm-1 is the fingerprint line, possibly of the α-AlB12 molecule, and the band at 1100 cm-1 can be assigned to the β-rhombohedral boron. The Raman band at 67.8 cm-1 (about 52.5 cm-1 lower) may result from the E2g mode of hexagonal h-BN.

If we take a Raman spectrum over the range 0-200 cm-1, we observe a Raman band at 51.8 cm-1, which is well fit to the E2g interlayer mode of h-BN. The bands at 1331.8 and 1450 cm-1 are due to B-doped diamond [9]. Because the spherical crystals (inclusions) are completely transparent, it is not possible to achieve high sp2 content. Because we have found more than three such oval-to-spherical crystals in one remnant of Ti-topaz (Figure 2 in Thomas 2026b [4], we can assume that the β-rhombohedral boron is doped with Ti (see also Werheit et al. 2010) (Figures 6 and 7) [7].

Figure 6: Raman spectrum of a small (5µm) colorless inclusion in topaz, mainly composed of rhombohedral boron. The presence of a diamond (1331.2 cm-1) component is only outlined.

Figure 7: Part of the Raman spectrum (Figures 4 and 5), which, according to Zaitsev (2001) [7], is interpreted as boron-doped diamond.

Another interpretation is more probable (because no phase boundary in the inclusion can be seen): the bands at 1121, 1147, and 1249 cm-1 can be assigned to a-AlB12 (see also Table 2); the 1332.3 cm-1 band corresponds to diamond, and the 1301 cm-1 band may be attributed to cubic BC2N.

Interpretation

The interpretation is only provisional because, with Raman alone, a straightforward interpretation is difficult. A large complication is the water-clear transparency of the studied inclusions. Boron is a black semimetal. Under high pressure between 19 and 89 GPa, the B-atoms in different crystal structures can rearrange [10], which is associated with changes in the optical properties. Nearby, all aluminum borides are not transparent [11]. That is also the case for boron carbides [10]. According to the same author [10], for example, e-boron became colorless and fully optically transparent at 62 GPa. A further complication is that the Raman line positions of boron carbides shift under high pressure, not all of them reversibly. At high pressure (27–40 GPa), the crystal structure changes significantly, and the material becomes increasingly transparent. A further complication arises from the possibility of orthorhombic boron oxide [12], which shows a small number of Raman lines that coincide with those of our transparent spherical crystals. Because at room temperature complete transparent spheres with Raman lines of boron and AlB12 are present, we interpret the transparent spheres as a mixture of α- and β-rhombohedral boron with small amounts of B4.3C, which moved by supercritical fluids or melts from mantle depths into the crust. Due to the rapid transport, some high-pressure Raman signatures did not relax.

Discussion

The results are very difficult to discuss because of the straightforward classification. Interpreting the obtained Raman spectroscopic results is not straightforward. Although the spherical inclusions in topaz from Schneckenstein appear homogeneous and water-clear, this poses a significant problem for boron compounds listed in Table 2. All are at room temperature, black, metallic, or very deep red translucent. The origin of this discrepancy may be phase transitions at 27–40 GPa, associated with an increase in transparency [10]. After that, the author observed that e-boron became colorless and fully optically transparent at 62 GPa. All these observations, together with the diamond spheres, demonstrate that the topaz is not a single crustal formation. The participation of mantle components via supercritical fluids (SCF) or supercritical melts (SCM) is necessary. The study of the Schneckenstein topaz shows that during the Variscan period, SCF and/or SCM contributed to mineralization to a greater extent than expected. The author has, in the last four years, shown many examples that proof there influence on the ore deposits [13-20], and references in it. Furthermore, the Schneckenstein topaz is not the only host for such spherical boron minerals. Figure 8 shows a similar Raman spectrum from an inclusion in cassiterite (Sn-58, Magdalena vein, Sauberg mine, Ehrenfriedersdorf; Sn-70; Sn-81) – see Thomas, 1982. Other cassiterites with orthorhombic cassiterite spheres (CaCl2-and cotunnite-types) often contain such B-inclusions, for example, in cassiterite (Sn-70; Thomas, 2024) [21-26]. Also, the cassiterite from Zinnwald contains rhombohedral boron spheres beside diamond, graphite, orthorhombic cassiterites and coesite.

Figure 8: Raman spectrum of a spherical boron mineral mixture in cassiterite from Ehrenfriedersdorf (Sn-58).

Generally, such boron spheres in the Variscan topaz and cassiterites are not rare. By careful screening, such boron minerals cannot be overlooked and demonstrate that the interaction between mantle domains and crust via SCF and/or SCM is a common process that must be borne in mind in any genetic discussion of the Variscan tin and related deposits. These observations are consistent with the preservation of unusual boron-rich Al–B–C microphases in Schneckenstein topaz and may point to high-pressure processes and transport by supercritical fluids or melts during Variscan mineralization; however, confirmation by complementary structural and chemical methods remains essential.

Acknowledgment

I thank all my colleagues with whom I could discuss the problems of interaction between SCF and SCM from the mantle region into the crust.

References

  1. Rösler HJ, Baumann L, Jung W (1968) Postmagmatic mineral deposits of the northern edge of the Bohemian Massif (Erzgebirge-Harz). International Geological Congress XXIII Session, Guide to Excursion 22 AC © German Democratic Republic. ZGI Berlin. Pg: 57.
  2. Lahl, B (2012) Königliche Topase vom Schneckenstein: Edelstein aus dem Chemnitzer Verlag. Pg: 141.
  3. Thomas R (1982) Ergebnisse der thermobarometrischen Untersuchungen an Flüssigkeitseinschlüssen in Mineralen der postmagmatischen Zinn-wolfram-mineralisation des Freiberger Forschungshefte (FFH) C370: 5-85. + I-XVI.
  4. Thomas R (2026b) Diamond in pycnite-rock from Altenberg and in the topaz from the famous Schneckenstein Geol Earth Mar Sci 8: 1-7.
  5. Ebel HF, Bliefert C, Russey WE (2004) The Art of Scientific Writing. Wiley-VCH-Verlag GmbH, Weinheim (second Edition). Pg: 1-595.
  6. Thomas R (2026a) Raman spectroscopic determination of fluorine in topaz. Geol Earth Mar Sci 8: 1-8.
  7. Werheit H, Filipov V, Kuhlmann U, Schwarz U, Aermbrüster M, et (2010) Raman effect in icosahedral boron-rich solids. Sci Technol Adv Mater 11: 1-27.
  8. Jay A, Dupare OH, Sjakste J, Vast N (2023) Theoretical Raman spectrum of boron carbide B43C under pressure. Elsevier, Acta Materialia. HALId: hal-04294014: 1-11.
  9. Zaitsev AM (2001) Optical Properties of Diamond -A Data Berlin, Heidelberg, New York. Springer. Pg: 502.
  10. Parakhonskiy G (2012) Synthesis and investigation of boron phases at high pressure and Dissertation, University of Bayreuth/Germany. Pg: 117.
  11. Kohn JA (1960). Crystallography of the Aluminum In: Kohn JA, Nye WF, Gaulé GK (eds) Boron Synthesis, Structure, and Properties. Springer, Boston, MA.
  12. Cherednichenko KA, Godec YL, Kalinko A, Mezouar M, Solozhenko VL (2016) Orthorhombic boron oxide under pressure: in situ study by X-ray diffraction and Raman J Appl Phys 120: 1-24.
  13. Thomas R (2026) Diamond in pycnite-rock from Altenberg and in the topaz from the famous Schneckenstein Geol Earth Mar Sci 8: 1-7.
  14. Thomas R, Rericha A (2025) Strong isotope fractionation between 13C and 12C in the supercritical fluids related to the Variscan mineralizations in Erzgebirge, Slavkovsky Les (Kaiserwald), and Lusatian Mountains, Germany, and the Czech Republic and some remarks on the low-pressure formation of diamond. Geol Earth Mar Sci 5: 1-5.
  15. Thomas R, Davidson P, Rericha A, Recknagel U (2022a) Water-rich coesite in prismatine-granulite from Waldheim/Saxony. Veröffentlichungen Museum für Naturkunde Chemnitz 45: 67-80.
  16. Thomas R, Davidson P, Rericha A, Recknagel U (2022b) Discovery of stishovite in the prismatine-bearing granulite from Waldheim, Germany: A possible role of supercritical fluids of ultrahigh-pressure origin. Geoscience 12: 1-13.
  17. Thomas R, Davidson P, Rericha A, Recknagel U (2023a) Ultrahigh-pressure mineral inclusions in a crustal granite: Evidence for a novel transcrustal transport mechanism. Geoscience 94: 1-13.
  18. Thomas R, Recknagel U, Rericha A (2023b) A moissanite-diamond-graphite paragenesis in a small beryl-quartz vein related to the Varsican tin-mineralization of the Ehrenfriedersdorf deposit, Aspects in Mining & Mineral Science 11: 1310-1319.
  19. Thomas R (2023a) Growth of SiC whiskers in beryl by a natural supercritical VLS Aspects in Mining & Mineral Science 11: 1292-1297.
  20. Thomas R (2023b) Unusual cassiterite mineralization, related to the Variscan tin-mineralization of the Ehrenfriedersdorf deposit, Aspects in Mining & Mineral Sciences 11: 1233-1236.
  21. Thomas R (2024) The CaCl2 to rutile phase transition in SnO2 from high to low pressure in Geol Earth Mar Sci 6: 1-4.
  22. Hubble OHW, Kudryashov I, Solozhenko VL, Zinin PV, Sharma SK, Ming LC (2004) Raman studies of cubic BC2N, a new superhard phase. Journal of Raman Spectroscopy 35: 822-825.
  23. Thomas R (2024a) The CaCl2-to-rutile phase transition in SnO2 from high to low pressure in Geol Earth Mar Sci 6: 1-4.
  24. Thomas R (2024b) Rhomboedric cassiterite as inclusions in tetragonal cassiterite from Slavkovský les – North Bohemia (Czech Republic). Earth Mar Sci 6: 1-6.
  25. Thomas R (2025a) Extremely 13C-rich diamond in orthorhombic cassiterites in the Variscan Erzgebirge, Saxony/Germany. Earth Mar Sci 7: 1-5.
  26. Thomas R (2025b) Diamond, diamond whisker, graphite, carbon, and coesite in a quartz crystal from Zinnwald, E-Erzgebirge. Geol Earth Mar Sci 7: 1-6.

Effect of Double-Coil Repetitive Peripheral Magnetic Stimulation on Knee Osteoarthritis: A Pilot Randomized Controlled Trial

DOI: 10.31038/IJOT.2026911

Abstract

Background: Knee osteoarthritis (KOA) is a leading cause of disability in older adults, with impaired walking ability significantly affecting independence and quality of life. Maintaining mobility is therefore a key therapeutic goal, particularly in patients awaiting total knee arthroplasty or managed conservatively. However, responses to standard non-invasive rehabilitation are often heterogeneous, highlighting the need for adjunctive therapeutic strategies.

Methods: Forty patients with KOA were recruited and randomly assigned to either a control group receiving standard rehabilitation or an experimental group receiving standard rehabilitation supplemented with double-coil repetitive peripheral magnetic stimulation (rPMS). Pain intensity, functional mobility outcomes, and knee range of motion were evaluated.

Results: Significant between-group differences in favor of the experimental group were observed for pain intensity, Timed Up and Go (TUG), and stair climb test (SCT) (p < 0.05), with the largest effect seen in pain reduction. Improvements in other outcomes consistently favored the experimental group. Responder analysis showed higher rates of clinically meaningful improvement in the experimental group.

Conclusions: This pilot study suggests that double-coil rPMS may contribute to pain reduction and improved mobility in patients with knee osteoarthritis, particularly in short-duration functional performance. Future studies should confirm these findings in larger populations and further investigate the comparative effectiveness of double-coil and single-coil rPMS approaches.

Keywords

Repetitive peripheral magnetic stimulation, Knee osteoarthritis, Total knee arthroplasty

Introduction

Osteoarthritis (OA) is a leading cause of disability in older adults, and its prevalence is expected to increase further due to population aging and rising obesity rates. The knee is the most commonly affected joint, followed by the hand and hip [1]. In addition to age and body weight, other important risk factors for knee osteoarthritis (KOA) include a history of joint injury and high bone mineral density [2]. Interestingly, unlike hip OA, the prevalence of KOA in women is strongly associated with postmenopausal age, particularly between 50 and 75 years [2]. The condition is characterized by chronic pain, reduced joint function, and progressive limitations in activities of daily living [1,2].

Impaired walking ability is one of the most clinically relevant consequences of KOA and has been associated with reduced quality of life and functional capacity compared to healthy individuals [3]. Maintaining mobility is therefore a key therapeutic goal, particularly in patients awaiting total knee arthroplasty or those managed conservatively in earlier stages of the disease. In older adults especially, the ability to walk with reduced pain and sufficient functional capacity is essential for preserving autonomy and preventing further physical decline [4].

Exercise-based, non-invasive rehabilitation is recommended as a core component of KOA management and has been shown to improve pain and function [5]. However, responses to rehabilitation are often heterogeneous, with considerable variability in outcomes and overall effects typically ranging from small to moderate [6,7]. This variability may reflect the presence of distinct patient subgroups and insufficient targeting of interventions, suggesting that current approaches may not fully address the complexity of knee joint dysfunction [8]. Additionally, treatment adherence can be limited by pain and related factors such as fear of movement, which may reduce patients’ ability to fully engage in exercise programs and contribute to inconsistent treatment effects [9-11]. These limitations highlight the need for adjunctive therapeutic strategies that are less dependent on active patient participation.

Repetitive peripheral magnetic stimulation (rPMS) is an emerging neuromodulatory technique that has been widely used in the management of spasticity and functional impairment in neurological conditions. However, its application in musculoskeletal disorders, including KOA, remains insufficiently explored [12,13]. Recently, a novel double-coil configuration with an adjustable angle between coils has been introduced, allowing improved adaptation to the treatment area, particularly in large joints such as the knee [14]. Mathematical simulations suggest that this configuration may deliver 45–121% greater energy to deep ligament tissue layers (3–8 cm) and up to 20% greater total magnetic energy in knee cartilage compared to a single-coil setup [14,15]. Preliminary clinical studies have indicated a potential benefit of double-coil rPMS in both upper and lower limb conditions, including KOA, with reported improvements in pain, functional disability, range of motion, and mobility [16,17].

However, to date, no randomized controlled trial has evaluated the effectiveness of this novel double-coil rPMS approach in combination with standard rehabilitation in patients with KOA, particularly with respect to both pain and functional mobility outcomes. Therefore, the aim of this pilot randomized controlled trial was to evaluate the feasibility and preliminary effectiveness of combining double-coil rPMS therapy with standard rehabilitation compared to standard rehabilitation alone in patients with KOA, with a focus on pain intensity and functional mobility.

Materials and Methods

Study Design

This study was designed as a randomized controlled trial conducted at the Rehamil Clinic (Milovice, Czech Republic) between March 2025 and February 2026. The study was carried out in accordance with the Declaration of Helsinki. All participants were informed about the study procedures and provided written informed consent prior to inclusion.

Participants

Participants were recruited from patients undergoing conservative treatment for KOA at the Rehamil Clinic. Eligible individuals were those diagnosed with KOA of Kellgren–Lawrence (KL) grades II–IV. No upper age limit was applied. To ensure sufficient functional capacity, only patients with at least 90° of active knee flexion and the ability to ambulate independently without assistive devices were included. The presence of osteoarthritis in other joints was allowed provided it was clinically stable and did not limit participation. Participants were required to be at least one month after intra-articular knee injection and at least six months after knee arthroscopy. Throughout the study, they were instructed to maintain their usual level of physical activity and stable analgesic medication; any changes in medication were recorded.

Exclusion criteria included pregnancy and the presence of implanted metallic or electronic devices (e.g., pacemakers, neurostimulators, or defibrillators). Additional exclusions comprised a history of seizures, active malignancy, systemic infection, or skin lesions in the treatment area. Patients with severe cardiovascular, pulmonary, or renal conditions, febrile illness, or other neurological or musculoskeletal disorders affecting lower limb function were also excluded.

Participants were randomly assigned to the control or experimental group using a computer-generated randomization algorithm. Group allocation and treatment scheduling were managed by the therapist responsible for delivering rPMS therapy, who was aware of group assignment. Therapists providing standard rehabilitation were not informed about group allocation. Outcome assessors were not involved in the randomization process.

Intervention

All participants underwent a standardized rehabilitation program consisting of 12 physiotherapy sessions. Each session lasted approximately 60 minutes and included exercise therapy (muscle strengthening and stretching), manual therapy (joint mobilization), and transcutaneous electrical nerve stimulation (TENS).

In addition to standard rehabilitation, participants in the experimental group received rPMS using a double-coil system (BTL SIS DUO, BTL Industries, Ltd., Prague, Czech Republic). Each rPMS session had a total duration of approximately 30 minutes, including preparation and positioning, with the active stimulation phase lasting approximately 13 minutes. The rehabilitation protocol was otherwise identical in both groups.

rPMS was applied with the patient in a supine position. The applicator was positioned around the knee joint with the coils placed medially and laterally, forming an approximately 90° angle to ensure coverage of the target area. A predefined treatment protocol intended for chronic musculoskeletal pain was used. Stimulation was delivered using a combination of modulation patterns with frequencies ranging approximately between 5 and 50 Hz, aiming to achieve both analgesic and neuromuscular effects. The intervention consisted of sequential phases including gradual intensity adjustment, pain-modulating stimulation, and phases supporting local circulation, followed by a gradual reduction in intensity at the end of the session. Stimulation intensity was individually adjusted based on patient tolerance and modified as needed during treatment.

Outcome Measures

Outcome measures were assessed at baseline and at a follow-up visit scheduled within 3 to 7 days after completion of the intervention.

Pain intensity was evaluated using the visual analogue scale (VAS; 0–10), where lower values indicate less pain [18]. Functional performance was assessed using the 30-second chair stand test (30SCHR; number of repetitions), Timed Up and Go test (TUG; seconds), and the stair climb test (SCT; seconds), with higher values indicating better performance for 30SCHR and lower values indicating better performance for TUG and SCT [19-21]. Walking performance was evaluated using the 40-meter walk test (40MWT; seconds) and the 6-minute walk test (6MWT; meters). For the 40MWT, shorter time indicates better performance, while for the 6MWT, longer distance reflects better walking capacity [22,23]. Knee range of motion (ROM; degrees) was assessed using standard goniometric measurement, with higher values indicating greater joint mobility. All assessments were performed under standardized conditions by the same examiner.

Changes between baseline and post-treatment values were analyzed for all outcomes. In addition to absolute changes, percentage changes were calculated to describe the magnitude of improvement across participants.

In addition to statistical significance, clinical relevance was evaluated using responder analysis based on minimal clinically important difference (MCID) thresholds, derived from available literature for the KOA population. The following MCID thresholds were applied: VAS ≥ 2-point reduction; 30SCHR ≥ 3 repetitions increase; 40MWT ≥ 8 s reduction; 6MWT ≥ 72 m increase; TUG ≥ 1 s reduction; SCT ≥ 2 s reduction; ROM ≥ 5° increase [22,24-28].

Sample Size and Statistical Analysis

Due to the pilot nature of the study, no a priori sample size calculation was performed. The sample size was determined pragmatically based on the capacity of the clinical setting and was considered sufficient to provide preliminary estimates of treatment effects.

Data processing, visualization, and statistical analyses were performed using R software (version 4.5.2). The distribution of all variables was assessed using normality testing. Variables with a normal distribution are presented as mean ± standard deviation (SD), while non-normally distributed variables are reported as median and interquartile range (IQR).

Within-group comparisons were performed using paired t-tests for normally distributed data and the Wilcoxon signed-rank test for non-normally distributed data. Between-group differences were analyzed using analysis of covariance (ANCOVA), with post-treatment values as the dependent variable, group as the independent factor, and baseline values included as covariates. Adjusted mean differences between groups, 95% confidence intervals (CI), and corresponding p-values were calculated.

In addition to statistical significance, clinical relevance was assessed using responder analysis based on MCID thresholds. For each outcome, the proportion of responders was compared between groups, and odds ratios (OR) with 95% confidence intervals were calculated using Fisher’s exact test. All tests were two-sided, and a p-value < 0.05 was considered statistically significant.

Results

A total of 40 patients were enrolled and evenly allocated to the control and experimental groups. Two participants, one in each group, did not complete the full treatment protocol due to discontinuation during the intervention period and were not included in the final outcome analysis. Baseline characteristics of both groups are presented in Table 1. The intervention was well tolerated, and no adverse events or other reasons for treatment discontinuation were reported. The complete patient flow, including recruitment, allocation, and analysis, is shown in Figure 1.

Table 1: Baseline characteristics of participants.

Variables

Control group (n=20)

Experimental group (n=20)

Age

63.35 ± 7.08 60.3 ± 8.68
Female sex, n (%) 15 (75%)

14 (70%)

BMI

28.85 ± 4.72 27.49 ± 6.26
KL grade, mean ± SD 2.55 ± 0.54

2.52 ± 0.64

BMI: Body Mass Index; KL grade: Kellgren–Lawrence Grade; SD: Standard Deviation.

Figure 1: Flow diagram of participant recruitment, allocation and analysis.

Within-group analysis demonstrated significant improvements across all outcome measures in the experimental group. In contrast, in the control group, changes in VAS and 30SCHR did not reach statistical significance. The magnitude of improvement, including percentage changes, consistently favored the experimental group across all outcomes (Tables 2 and 3). The largest treatment effect was observed for pain intensity. Changes in pain intensity over the study period are illustrated in Figure 2. Although the experimental group presented with higher baseline VAS values, it achieved markedly lower post-treatment scores compared to the control group.

Table 2: Within-group changes in clinical outcomes of experimental group.

Experimental group (n=19)
Outcome Baseline Post-treatment Δ Δ%

p-value

VAS

4.53 ± 2.32 1.84 ± 1.38 -2.68 ± 1.83 -53% ± 28% <0.001
30SCHR 11.37 ± 3.04 14.37 ± 3.04 3.00 ± 2.45 31% ± 27%

<0.001

40MWT (s)

40.40 ± 10.54 33.09 ± 6.34 -7.31 ± 5.51 -16% ± 11% <0.001
TUG (s) 7.05 ± 2.34 5.09 ± 1.60 -1.96 ± 1.57 -26% ± 18%

<0.001

6MWT (m)

343.96 ± 94.61 426.77 ± 96.58 82.81 ± 51.60 27% ± 18% <0.001
SCT (s) 12.37 ± 5.55 8.91 ± 2.73 -3.47 ± 3.30 -24% ± 15%

<0.001

ROM (°)

125.00 (7.50) 130.00 (7.50) 5.00 (8.50) 4% (7%)

<0.001

Values are presented as mean ± standard deviation, except for ROM, which is reported as median (interquartile range). Δ represents the absolute change between baseline and post-treatment values, and Δ% represents the percentage change relative to baseline. Negative values indicate improvement for VAS, 40MWT, TUG, and SCT, while positive values indicate improvement for 30SCHR, 6MWT, and ROM. Within-group differences were assessed utilizing paired t-tests, with the exception of ROM, which was evaluated using the Wilcoxon signed-rank test. VAS: Visual Analogue Scale; 30SCHR: 30-Second Chair Stand Test; 40MWT: 40-Meter Walk Test; TUG: Timed Up and Go; 6MWT: 6-Minute Walk Test; SCT: Stair Climb Test; ROM: Range of Motion.

Table 3: Within-group changes in clinical outcomes of control group.

Control group (n=19)
Outcome Baseline Post-treatment Δ Δ%

p-value

VAS

3.05 ± 1.84 2.68 ± 2.08 -0.37 ± 2.24 14% ± 133% 0.483
30SCHR 13.63 ± 3.64 14.79 ± 3.81 1.16 ± 2.75 10% ± 19%

0.08

40MWT (s)

37.32 ± 9.19 33.20 ± 8.88 -4.12 ± 4.74 -10% ± 12% <0.001
TUG (s) 6.53 ± 2.02 5.62 ± 1.77 -0.90 ± 1.35 -13% ± 18%

0.009

6MWT (m)

386.66 ± 112.41 448.90 ± 101.66 62.25 ± 51.69 18% ± 16% <0.001
SCT (s) 10.23 ± 4.38 8.76 ± 3.59 -1.47 ± 1.46 -13% ± 11%

<0.001

ROM (°)

125.00 (10.00) 130.00 (5.00) 0.00 (5.00) 0% (4%)

0.006

Values are presented as mean ± standard deviation, except for ROM, which is reported as median (interquartile range). Δ represents the absolute change between baseline and post-treatment values, and Δ% represents the percentage change relative to baseline. Negative values indicate improvement for VAS, 40MWT, TUG, and SCT, while positive values indicate improvement for 30SCHR, 6MWT, and ROM. Within-group differences were assessed utilizing paired t-tests, with the exception of ROM, which was evaluated using the Wilcoxon signed-rank test. VAS: Visual Analogue Scale; 30SCHR: 30-Second Chair Stand Test; 40MWT: 40-Meter Walk Test; TUG: Timed Up and Go; 6MWT: 6-Minute Walk Test; SCT: Stair Climb Test; ROM: Range of : Motion.

Figure 2: Changes in pain intensity (VAS) over time in the control and experimental groups. Values represent mean ± 95% confidence interval at baseline and post-treatment.

These between-group differences were confirmed by ANCOVA, which demonstrated statistically significant effects in favor of the experimental group for VAS, TUG, and SCT (Table 4). The forest plot (Figure 3) further illustrates that adjusted between-group effects consistently favored the experimental group across all outcomes, although not all reached statistical significance. This pattern was supported by responder analysis, which showed higher MCID responder rates in the experimental group across all outcomes. Statistically significant differences were observed for VAS and 40MWT (Table 5). Individual changes in pain intensity are presented in Figure 4.

Table 4: Between-group comparison of outcomes (ANCOVA).

Outcome

Control mean Experimental mean Difference (95% CI)* p-value
VAS 2.96 1.56 1.40 (0.28 to 2.52)

0.016

30SCHR

14.0 15.2 -1.22 (-2.95 to 0.50) 0.159
40MWT (s) 34.2 32.1 2.17 (-0.51 to 4.86)

0.109

TUG (s)

5.77 4.94 0.83 (0.08 to 1.59) 0.032
6MWT (m) 431 445.0 -13.4 (-46.56 to 19.80)

0.418

SCT

9.37 8.29 1.08 (0.13 to 2.04) 0.028
ROM (°) 127.0 129.0 -2.41 (-5.85 to 1.04)

0.165

Adjusted means were estimated using analysis of covariance (ANCOVA) with baseline values included as a covariate. Differences are presented as control minus experimental group. Positive values indicate better outcomes in the experimental group for VAS, TUG, 40MWT, and SCT, whereas negative values indicate better outcomes for 30SCHR, 6MWT, and ROM. CI denotes 95% confidence interval. VAS: Visual Analogue Scale; 30SCHR: 30-Second Chair Stand Test; 40MWT: 40-Meter Walk Test; TUG: Timed Up and Go; 6MWT: 6-Minute Walk Test; SCT: Stair Climb Test; ROM: Range of Motion.

Figure 3: Adjusted between-group differences in clinical outcomes. Points represent adjusted mean differences estimated using ANCOVA, and horizontal lines indicate 95% confidence intervals. Positive values favor the experimental group. VAS: Visual Analogue Scale; 30SCHR: 30-Second Chair Stand Test; 40MWT: 40-Meter Walk Test; TUG: Timed Up and Go; 6MWT: 6-Minute Walk Test; SCT: Stair Climb Test ROM: Range of Motion.

Table 5: Proportion of patients achieving clinically meaningful improvement (MCID).

Outcome

Control, n (%) Experimental, n (%) Odds ratio (95% CI) p-value
VAS 4/19 (21.1%) 14/19 (73.7%) 9.72 (1.93 to 62.58)

0.003

30SCHR

6/19 (31.6%) 10/19 (52.6%) 2.35 (0.54 to 11.20) 0.325
40MWT 2/19 (10.5%) 10/19 (52.6%) 8.86 (1.44 to 100.00)

0.013

6MWT

7/19 (36.8%) 13/19 (68.4%) 3.58 (0.81 to 17.70) 0.103
TUG 9/19 (47.4%) 13/19 (68.4%) 2.35 (0.54 to 11.20)

0.325

SCT

6/19 (31.6%) 12/19 (63.2%) 3.58 (0.81 to 17.70) 0.103
ROM 9/19 (47.4%) 13/19 (68.4%) 2.35 (0.54 to 11.20)

0.325

Responders were defined as patients achieving a clinically meaningful improvement based on predefined minimal clinically important difference (MCID) thresholds. The following thresholds were used: VAS ≥ 2-point reduction; 30SCHR ≥ 3 repetitions increase; 40MWT ≥ 8 s reduction; 6MWT ≥ 72 m increase; TUG ≥ 1 s reduction; SCT ≥ 2 s reduction; ROM ≥ 5° increase. Odds ratios (OR) and 95% confidence intervals (CI) were calculated using Fisher’s exact test. VAS: Visual Analogue Scale; 30SCHR: 30-Second Chair Stand Test; 40MWT: 40-Meter Walk Test; TUG: Timed Up and Go; 6MWT: 6-Minute Walk Test; SCT: Stair Climb Test; ROM: Range of Motion.

Figure 4: Individual changes in pain intensity (VAS) in the control and experimental groups. Each point represents one patient plotted according to baseline and post-treatment values. The solid diagonal line indicates no change (y=x), while the dashed line represents the threshold for clinically meaningful improvement (≥2-point reduction in VAS). Points below the dashed line indicate patients achieving clinically meaningful pain reduction.

Discussion

This pilot randomized controlled trial demonstrated the added value of double-coil rPMS as an adjunct to standard rehabilitation in patients with KOA. The combined intervention resulted in the greatest effect in pain reduction compared to standard rehabilitation alone. Statistically significant between-group differences were also observed for functional outcomes, specifically the TUG and the SCT.

Overall, treatment outcomes consistently favored the experimental group across all measured parameters, both in terms of percentage improvement and MCID responder rates. Patients receiving the combined intervention appeared to benefit across multiple domains, with the most pronounced differences observed in pain intensity and functional mobility, particularly in TUG, SCT, and 40MWT. These findings suggest that the addition of rPMS may enhance pain modulation and improve short-duration functional performance, which are likely more sensitive to changes in neuromuscular activation and reduced pain-related inhibition during movement. This is consistent with previous evidence showing that neuromuscular stimulation preferentially improves muscle activation and short-term functional tests such as the TUG, while having less consistent effects on longer-duration performance measures such as the 6MWT [29-32].

Although the available evidence on rPMS in the treatment of KOA remains limited, the findings of the only existing study evaluating double-coil rPMS are consistent with the results of the present study. Bednar et al. reported comparable percentage improvements in pain and range of motion, reaching approximately 50% and 4%, respectively, which closely aligns with the outcomes observed in this trial [17]. A greater improvement in TUG was observed in the present study (26% vs. 16%), which may be explained by the addition of a standardized rehabilitation protocol, potentially enabling greater functional gains. In contrast, Bednar et al. investigated the effect of double-coil rPMS as a standalone intervention. Direct comparison of MCID outcomes is limited due to differences in threshold definitions. However, a similar trend toward higher responder rates, approaching 70%, was observed in both studies. To date, no other studies have evaluated this specific intervention in patients with KOA, limiting the possibility of broader comparison [17].

The limited number of randomized controlled trials investigating rPMS in patients with KOA may be partly explained by the technical limitations of conventional single-coil systems, which have been predominantly used to date. These systems primarily target superficial tissues, typically within a depth of approximately 3 cm, and may therefore be less effective for deeper joint structures such as the knee. The recently introduced double-coil design allows for partial overlap of the generated magnetic fields when appropriately positioned, potentially enabling higher stimulation intensities at greater depths (approximately 3–8 cm). To date, this concept has been supported mainly by experimental simulations and preliminary clinical studies. Direct clinical comparison between single-coil and double-coil approaches in patients with KOA remains lacking, and further research is needed to clarify the potential advantages of this technology.

Several limitations of this study should be acknowledged. First, the pilot design and relatively small sample size limit the generalizability of the findings and increase the risk of overestimating treatment effects. In addition, no a priori sample size calculation was performed, and the sample size was determined pragmatically based on the capacity of the clinical setting. Second, this was a single-center study, which may further limit the external validity of the results. The follow-up period was restricted to the immediate post-treatment assessment, and therefore no conclusions can be drawn regarding the long-term sustainability of the observed effects. Third, blinding was not fully implemented, as the therapist responsible for delivering rPMS was aware of group allocation. Although therapists providing standard rehabilitation were not informed, the potential for performance bias cannot be excluded.

Despite its pilot nature, this study provides preliminary evidence supporting the potential clinical value of double-coil rPMS in the treatment of KOA. Future research should confirm these findings in larger populations and investigate the comparative effectiveness of double-coil and single-coil rPMS approaches.

Conclusions

This pilot study suggests that double-coil rPMS may contribute to pain reduction and improved mobility in patients with knee osteoarthritis, particularly in short-duration functional performance. Supporting mobility in this population is essential for maintaining independence in daily activities and may help prevent further physical and psychological decline, especially during the period preceding total knee arthroplasty. Future studies should confirm these findings in larger populations and further investigate the comparative effectiveness of double-coil and single-coil rPMS approaches.

References

  1. Hunter DJ, Bierma-Zeinstra S (2019) Osteoarthritis Lancet.
  2. Duong V, Abdel Shaheed C, Ferreira ML, Narayan SW, Venkatesha V, Hunter DJ, Zhu J, Atukorala I, Kobayashi S, Goh SL, Briggs AM, Cross M, Espinosa-Morales R, Fu K, Guillemin F, Keefe F, Stefan Lohmander L, March L, Milne GJ, Mei Y, et al. (2025) Risk factors for the development of knee osteoarthritis across the lifespan: A systematic review and meta-analysis. Osteoarthritis and Cartilage. [crossref]
  3. Pavel RMS, Purza AL, Tit DM, Radu AF, Iovanovici DC, Vasileva D, Uivaraseanu B, Bungau G, Nistor-Cseppento CD (2025) Functional Burden and Quality of Life in Hip and Knee Osteoarthritis: A Cross-Sectional Study. Medicina. [crossref]
  4. Shafrin J, Sullivan J, Goldman DP, Gill TM (2017) The association between observed mobility and quality of life in the near elderly. PLoS One. [crossref]
  5. Sánchez J, Álvarez-Barbosa F, Rodríguez M, Ramos-Munell J (2025) Physical exercise to manage pain in patients with knee osteoarthritis: a systematic review. Archivos de Medicina del Deporte
  6. Quicke JG, Runhaar J, Windt D, Healey E, Foster NE, Holden M (2020) Moderators of the effects of exercise for hip and knee osteoarthritis: a systematic review of subgroup analyses from randomised controlled trials. Osteoarthritis and Cartilage. [crossref]
  7. Yazdi F, Salajegheh A, Yahyaabadi F, Salajegheh P (2025) The Impact of Strength/ Resistance Training on Related Outcomes of Patients with Knee Osteoarthritis: A Systematic Review and Meta-Analysis Study. Galen Medical Journal. [crossref]
  8. Holden MA, Hattle M, Runhaar J, Riley RD, Healey EL, Quicke J, van der Windt DA, Dziedzic K, van Middelkoop M, Burke D, Corp N, Legha A, Bierma-Zeinstra S, Foster NE, STEER OA Patient Advisory Group, OA Trial Bank Exercise Collaborative (2023) Moderators of the effect of therapeutic exercise for knee and hip osteoarthritis: a systematic review and individual participant data meta-analysis. The Lancet Rheumatology. [crossref]
  9. Jinks C, Moore A, Holden M, Foster N (2016) Identifying barriers and facilitators to exercise and physical activity in older adults with knee pain: Findings from the BEEP trial longitudinal qualitative study. Osteoarthritis and Cartilage
  10. Hinman R, Jones S, Nelligan R, Campbell P, Hall M, Foster N, Russell T, Bennell K (2023) Absence of Improvement With Exercise in Some Patients With Knee Osteoarthritis: A Qualitative Study of Responders and Arthritis Care & Research. [crossref]
  11. Bennell K, Nelligan R, Schwartz S, Kasza J, Kimp A, Crofts S, Hinman R (2020) Behavior Change Text Messages for Home Exercise Adherence in Knee Osteoarthritis: Randomized Trial. Journal of Medical Internet Research. [crossref]
  12. Sakai K, Yasufuku Y, Kamo T, Ota E, Momosaki R (2019) Repetitive peripheral magnetic stimulation for impairment and disability in people after stroke. The Cochrane Database of Systematic Reviews. [crossref]
  13. Fawaz SI, Izumi SI, Zaki AS, et (2023) Repetitive peripheral magnetic stimulation for improving upper limb function in post-stroke hemiparesis. Egypt Rheumatol Rehabil
  14. Vrbová B, Malá J, Barna M (2025) Comparison of single- and double-coil rPMS configurations in physiotherapy: Perception and field distribution in commonly treated J Orthop Sports Med
  15. Vrbová B, Malá J (2026) Single versus double-coil repetitive peripheral magnetic stimulation targeting knee cartilage and ligaments: A finite element analysis. J Clin Diagn Res
  16. Modica V, Sciarrone G, Barna M (2026) A Novel Double-Coil rPMS Approach for Treating Joint Disorders of the Upper and Lower Extremities: A Pilot Open Journal of Therapy and Rehabilitation. [crossref]
  17. Bednár R, Flašková M, Fejková N (2026) A Pilot Study Investigating Clinical and Functional Outcomes of Novel Double-Coil rPMS in Knee Osteoarthritis. Biomedicines. [crossref]
  18. Price DD, McGrath PA, Rafii A, Buckingham B (1983) The validation of visual analogue scales as ratio scale measures for chronic and experimental pain. Pain. [crossref]
  19. Jones CJ, Rikli RE, Beam WC (1999) A 30-s Chair-Stand Test as a Measure of Lower Body Strength in Community-Residing Older Adults. Research Quarterly for Exercise and Sport. [crossref]
  20. Podsiadlo D, Richardson S (1991) The Timed “Up & Go”: A Test of Basic Functional Mobility for Frail Elderly J Am Geriatr Soc. [crossref]
  21. Sharma S, Wilson R, Pryymachenko Y, Pathak A, Chua J, Gwynne-Jones D, Bennell KL, Metcalf B, Hinman RS, Hunt MA, Lim BW, Abbott JH (2023) Reliability, Validity, Responsiveness, and Minimum Important Change of the Stair Climb Test in Adults with Hip and Knee Arthritis Care Res. [crossref]
  22. Ramalho RB, Chaves TC, Terluin B, Selistre LFA (2025) Minimal Important Changes of Common Outcome Measures of Physical Function in Individuals With Knee Osteoarthritis: A Prospective Clinical Archives of Physical Medicine and Rehabilitation. [crossref]
  23. Wang TJ, Belza B, Thompson FE, Whitney JD, Bennett K (2007) Effects of Aquatic Exercise on Flexibility, Strength and Aerobic Fitness in Adults with Osteoarthritis of the Hip or J Adv Nurs. [crossref]
  24. Franceschini M, Boffa A, Pignotti E, Andriolo L, Zaffagnini S, Filardo G (2023) The Minimal Clinically Important Difference Changes Greatly Based on the Different Calculation The American Journal of Sports Medicine. [crossref]
  25. King L, Waugh E, Woodhouse L, Hawker G, Jones A (2020) What change in six-minute walk test indicates a clinically meaningful improvement in mobility to patients with knee osteoarthritis after total knee arthroplasty? Osteoarthritis and Cartilage. [crossref]
  26. Alghadir A, Anwer S, Brismée JM (2015) The reliability and minimal detectable change of Timed Up and Go test in individuals with grade 1-3 knee BMC Musculoskeletal Disorders. [crossref]
  27. Pajaniaye JB, Cheret EP, Langvad CH, Holm PM, Skou ST, Larsen JB, Mechkenburg I (2026) Knee osteoarthritis and minimal important change for the nine-step stair climb Danish Medical Journal. [crossref]
  28. Wright AA, Cook CE, Baxter GD, Dockerty JD, Abbott JH (2011) A comparison of methodological approaches to defining clinically important improvement of performance measures in hip osteoarthritis. J Orthop Sports Phys Ther. [crossref]
  29. Peng L, Wang K, Zeng Y, Wu Y, Si H, Shen B (2021) Effect of Neuromuscular Electrical Stimulation After Total Knee Arthroplasty: A Systematic Review and Meta-Analysis of Randomized Controlled Frontiers in Medicine. [crossref]
  30. Kourek C, Kanellopoulos M, Raidou V, Antonopoulos M, Karatzanos E, Patsaki I, Dimopoulos S (2024) Safety and effectiveness of neuromuscular electrical stimulation in cardiac surgery: A systematic review. World Journal of Cardiology. [crossref]
  31. Capin J, Bade M, Jennings J, Snyder-Mackler L, Stevens-Lapsley J (2022) Total Knee Arthroplasty Assessments Should Include Strength and Performance-Based Functional Tests to Complement Range-of-Motion and Patient-Reported Outcome Physical Therapy. [crossref]
  32. Nussbaum E, Houghton P, Anthony J, Rennie S, Shay B, Hoens A (2017) Neuromuscular Electrical Stimulation for Treatment of Muscle Impairment: Critical Review and Recommendations for Clinical Physiotherapy Canada. [crossref]

Diamond and Carbon Whiskers in Hydrothermal Quartz from Zinnwald, Germany: Evidence for Disequilibrium Fluid-controlled Growth

DOI: 10.31038/GEMS.2026853

Abstract

This study documents the presence of diamond, diamond-like carbon (DLC), graphite, and carbon whiskers in the outer growth zone of a smoky quard crystal from Zinnwald, eastern Erzgebirge, Germany. Optical microscopy and Raman spectroscopy show that these carbon phases occur as needle- and whisker-like crystals concentrated in the late stage of quard growth. Coesite associated with the assemblage indicates pressures of about 2 GPa, which are too low for the equilibrium graphite-to-diamond transformation. In addition, fluid inclusions containing CO2, CH4, and H2 point to a carbon-bearing C–O–H fluid system. The coexistence of multiple carbon phases, their unusual morphologies, and their restriction to the final growth zone are best explained by rapid disequilibrium precipitation from a supercritical fluid or melt, rather than by stable low-pressure diamond growth. The Zinnwald occurrence, therefore, supports fluid-controlled carbon deposition and may reflect mantle–crust interaction during Variscan mineralization.

Keywords

Diamonds, Carbon whiskers, Hydrothermal quartz, Raman spectroscopy, C-O-H fluids, Disequilibrium crystallization, Zinnwald

Introduction

This contribution examines diamond, diamond-like carbon (DLC), graphite, and carbon whiskers in the outer growth zone of a smoky quartz crystal from Zinnwald, eastern Erzgebirge, Germany. The central question is whether these carbon phases record equilibrium diamond growth or, instead, rapid disequilibrium precipitation from a carbon-bearing fluid or melt during late quartz crystallization. A concise description of the pre-Mesozoic geology of the Saxo-Thuringia zone can be taken from Linneman and Romer (2010) [1]. Previous work from Variscan mineralizations in the Erzgebirge and related regions has reported high-pressure mineral associations, including diamond, lonsdaleite, and orthorhombic cassiterite, and has been interpreted as evidence for mantle–crust interaction mediated by supercritical fluids or melts. So the found CaCl2-type and cotunnite-type cassiterite give pressures of 12 and 54 GPa, respectively [2-4]. Against that background, the Zinnwald quartz sample is important because it contains mixed carbon phases in a late-growth zone, together with coesite and C–O–H-bearing fluid inclusions (Table 1).

Table 1: Diamonds in orthorhombic cassiterite from the Sauberg mine near Ehrenfriedersdorf.

Sample

First-order Raman band FWHM n
Sn-58 1334.3 ± 4.1 cm⁻¹ 33.8 ± 4.8 cm⁻¹

6

Sn-70

1326.4 ± 4.8 cm⁻¹ 64.4 ± 15.2 cm⁻¹

11

FWHM – Full-Width at Half Maximum.

The Zinnwald occurrence is therefore significant because diamond, graphite, DLC, and carbon whiskers appear in the final growth zone of a hydrothermal quartz crystal, whereas the coesite-constrained pressure of about 2 GPa is too low for equilibrium graphite-to-diamond transformation. This mismatch motivates the interpretation developed here: the carbon phases most likely formed by disequilibrium, fluid-controlled precipitation rather than by stable low-pressure diamond growth.

Sample Material

The sample material is detailed described in Thomas (2025b) [5]. It is a 6 cm-long smoky quartz crystal from a pegmatitic vein. From this crystal, we prepared a 500 µm-thick, both-side-polished section for the microscopic and Raman spectroscopic studies (Figure 1).

Figure 1: Thick section of the used quartz sample (ZQ-2) from Zinnwald, E-Erzgebirge. Qtz-o is the older core of the a-quartz crystal. D-C-rich: region with many whisker-like crystals of diamond, carbon, graphite, and DLC. The quartz, apart from the older core, shows a pronounced mosaic structure.

Methods: Microscopy and Raman Spectroscopy

We performed all microscopic and Raman spectroscopic studies with a petrographic polarization microscope (BX 43) with a rotating stage coupled with the EnSpectr Raman spectrometer R532 (Enhanced Spectrometry, Inc., Mountain View, CA, USA) in reflection and transmission. We used Raman spectroscopy as the main tool for characterizing the carbon materials (diamond, DLC, graphite, and carbon). According to Zaitsev (2001) [6], the more correct term for DLC is “diamond-like materials.” The data presented here were acquired using an EnSpectr Raman microscope (RamMics R532). Measurements covered the spectral range from 0 to 4000 cm⁻¹ and were conducted with a single-mode 532 nm laser operating at a maximum output of 50 mW. Analytical settings included a 20 µm entrance aperture, a holographic grating of 1800 g mm⁻¹, and a spectral resolution of approximately 4 cm⁻¹ (generally for the whole range). For most analyses, a long-working-distance Olympus LMPlanFL 100× objective is used. Laser power at the sample surface was continuously adjustable down to 0.02 mW. Higher powers (up to 50 mW) were applied only for overview measurements. Generally, for measurements, we used 0.9 mW on the sample over the range 750-1450 cm⁻¹ to prevent local heating. Raman band positions were calibrated before and after each measurement series using the Si band of a semiconductor-grade single-crystal silicon chip. Based on 20 repeated measurements, run-to-run reproducibility was ± 0.2 cm⁻¹ for silicon (520.2 ± 0.2 cm⁻¹) in the measuring range 100 to 1450 cm⁻¹. As a second reference, we used a water-clear diamond crystal from Brazil. For the first order diamond line, we obtained (1330.1 ± 0.6) cm⁻¹ with a FWHM = 5.1 ± 0.1 cm⁻¹.

Results

The outer growth zone of the quartz crystal contains abundant needle- and whisker-like crystals of carbon, graphite, DLC, and diamond (Figures 2-6). These phases are concentrated in the final growth zone and are absent or rare toward the crystal interior, except for a single colorless diamond whisker reported near the center [5]. Coesite identified at the end of quartz growth indicates pressures of about 2 GPa, which are insufficient for equilibrium diamond formation from graphite. Most fluid inclusions are secondary and contain a vapor phase dominated by CO2. A smaller number of high-temperature inclusions also contain H2, identified from the 354 cm⁻¹ band, and CH4, identified from the 2917 cm⁻¹ band. Together, these inclusions indicate a carbon-bearing C–O–H fluid system associated with late quartz growth. Next, we show some examples of whisker-like diamond, DLC, graphite, and carbon. Such an appearance is exceptional. Also, the often observed bending is unusual.

Figure 2: Carbon needles in rand zone of the quartz crystal. Some needles are bent.

Figure 3: Graphite (Gr) needles in the same growth zone as Figure 2.

Figure 4: DLC-needles, often bent, in the outer growth zone of the quartz crystal. The most needles have a brownish shade.

Figure 5: Diamond whiskers (DW), often colorless to brownish, in the outer growth zone of the quartz.

Using the needle in Figure 6, we have performed Raman measurements in the first-order diamond range. The results are in Table 2. To compare the data, the measured values from Thomas (2025b) [5] are included. Conspicuously, the FWHM values for all carbon types shown here are relatively large. That can be traced back to higher portions of sp2 (one carbon atom forms three bonds in one plane, additionally one p-orbital for one π-bond remaining). Pure diamond material consists mainly of sp3 (four equivalent bonds in a tetrahedral arrangement). A high proportion of sp2 leads to an increase in the FWHM.

Figure 6: A brownish translucent diamond-carbon needle in the outer quartz growth zone.

Table 2: Raman measurements in the first-order diamond range.

Sample

First order D band (cm⁻¹) FWHM (cm⁻¹) G-band (cm⁻¹) FWHM (cm⁻¹) n
Needle Fig. 6 1330.1 ± 4.3 73.6 ± 12.0 1582.8 ± 4.8 67.7 ± 7.7

20

Thomas 2025b

1321.6 ± 3.8 51.6 ± 6.7 1552.9 ± 7.1 47.9 ± 4.3 12
Thomas 2025b 1321.6 ± 7.0 14.3 ± 0.8 1575.6 ± 6.5 51.7 ± 13.2

14

As already mentioned, the estimated pressure (2 GPa) from the presence of coesite is insufficient for the diamond formation at the equilibrium graphite-diamond. Figure 7 shows a typical high-temperature fluid inclusion.

Figure 7: Complex high-temperature fluid inclusion in the Zinnwald quartz consists of carbon dioxide (CO2), methane (CH4), and hydrogen (H2). The water-rich phase (H2O) consists of CO2 and higher hydrocarbons. Note the low contrast between the quartz host and the water-rich phase.

Interpretation

Disequilibrium Growth from C–O–H Fluids Rather than Direct Graphite-to-Diamond Conversion

The occurrence of diamond in the Zinnwald quartz is unlikely to reflect equilibrium crystallization along the graphite–diamond boundary. The coesite-bearing assemblage indicates pressures of about 2 GPa, which are too low for direct graphite-to-diamond conversion under equilibrium conditions. A more consistent explanation is precipitation from a carbon-rich supercritical C–O–H fluid or melt, in which redox state, carbon speciation, and rapid disequilibrium changes controlled carbon deposition. This interpretation is supported by four observations: mixed carbon phases, whisker-like morphologies, evidence for rapid growth, and fluid inclusions rich in C–O–H species.

  • fluid inclusions containing CH4, H2, and CO2,
  • needle- and whisker-like crystal morphologies,
  • the coexistence of diamond, DLC, graphite, and carbon,
  • and their concentration in the final growth zone of the quartz

Taken together, these features argue against stable low-pressure diamond growth. They are more consistent with a strongly disequilibrium fluid system in which carbon precipitated rapidly from a supercritical fluid or melt during late quartz growth. Transient local overpressure at microsites, or the transport of pre-existing, deeper-formed diamond by an ascending fluid, may also have contributed. Thus, the Zinnwald occurrence is best explained by fluid-controlled, disequilibrium carbon deposition rather than by simple equilibrium graphite-to-diamond transformation.

Discussion

The key result of this study is that the Zinnwald carbon assemblage is best explained by disequilibrium precipitation during late quartz growth, not by equilibrium graphite-to-diamond transformation. This conclusion follows from three linked observations: coesite constrains pressure to about 2 GPa, which is too low for equilibrium diamond formation; the carbon phases occur as mixed needle- and whisker-like morphologies concentrated in the final growth zone; and associated inclusions record a carbon-bearing C–O–H fluid containing CO2, CH4, and H2. Together, these features support rapid, heterogeneous carbon deposition from a supercritical fluid or melt in which redox fluctuations, carbon speciation, and kinetic effects control phase formation. Transient local overpressure or transport of deeper-formed diamond by an ascending fluid remain possible contributing mechanism, but they do not alter the main inference that fluid-controlled disequilibrium processes dominated carbon deposition in the Zinnwald quartz. This interpretation also fits the broader regional context. Nearby major tin deposits, including Ehrenfriedersdorf, Altenberg, Krupka, and Slavkovský les, as well as the Schneckenstein deposit, have yielded spherical diamonds [7-13], and related Variscan tin mineralizations also contain orthorhombic cassiterite [2,3] and diamond enclosed in topaz [14,15]. Taken together, these occurrences strengthen the view that some Variscan mineral assemblages record transcrustal carbon transport and mantle–crust interaction rather than simple closed-system crystallization.

Limitations

This interpretation should nevertheless be regarded as provisional because the present study is based on a limited sample set and primarily on optical observations and Raman spectroscopy. Although the coexistence of diamond, DLC, graphite, carbon whiskers, coesite, and C–O–H-bearing fluid inclusions is consistent with disequilibrium, fluid-controlled growth, the data do not yet fully exclude alternative explanations, such as the transport of deeper-formed diamond or a transient local overpressure during quartz growth. Additional work, including broader sampling, quantitative Raman mapping, microstructural analysis, and isotopic measurements, would be needed to test the generality of the Zinnwald occurrence and further constrain the origin of the carbon phases.

Acknowledgment

I thank Michael Leh in D-02699 Neschwitz, Germany for his longstanding interest and support of my home office work.

References

  1. Linnemann U, Romer RL (2010) Pre-Mesozoic Geology of Saxo-Thuringia. Schweizerbart Science Pg: 488.
  2. Thomas R (2024a) The CaCl2-to-rutile phase transition in SnO2 from high to low pressure in Geol Earth Mar Sci 6: 1-4.
  3. Thomas R (2024b) Rhomboedric cassiterite as inclusions in tetragonal cassiterite from Slavkovský les – North Bohemia (Czech Republic). Geol Earth Mar Sci 6: 1-6.
  4. Yehdego DT (2009) Pressure induced phase transformation of SnO2 an ab initio constant pressure Master of Science Thesis, University of Texas at el Paso. Pg: 54.
  5. Thomas R (2025b) Diamond, diamond whisker, graphite, carbon, and coesite in a quartz crystal from Zinnwald, E-Erzgebirge. Geol Earth Mar Sci 7: 1-6.
  6. Zaitsev AM (2001) Optical Properties of Diamond -A Data Berlin, Heidelberg, New York, Springer. Pg: 502.
  7. Thomas R, Davidson P, Rericha A, Recknagel U (2022a) Water-rich coesite in prismatine-granulite from Waldheim/Saxony. Veröffentlichungen Museum für Naturkunde Chemnitz 45: 67-80.
  8. Thomas R, Davidson P, Rericha A, Recknagel U (2022b) Discovery of stishovite in the prismatine-bearing granulite from Waldheim, Germany: A possible role of supercritical fluids of ultrahigh-pressure origin. Geoscience 12: 1-13.
  9. Thomas R (2023a) Growth of SiC whiskers in beryl by a natural supercritical VLS Aspects in Mining & Mineral Science 11: 1292-1297.
  10. Thomas R, Davidson P, Rericha A, Recknagel U (2023a) Ultrahigh-pressure mineral inclusions in a crustal granite: Evidence for a novel transcrustal transport mechanism. Geoscience 94: 1-13.
  11. Thomas R (2023b) Unusual cassiterite mineralization, related to the Variscan tin-mineralization of the Ehrenfriedersdorf deposit, Aspects in Mining & Mineral Sciences 11: 1233-1236.
  12. Thomas R, Recknagel U, Rericha A (2023b) A moissanite-diamond-graphite paragenesis in a small beryl-quartz vein related to the Varsican tin-mineralization of the Ehrenfriedersdorf deposit, Aspects in Mining & Mineral Science 11: 1310-1319.
  13. Thomas R, Rericha A (2025) Strong isotope fractionation between 13C and 12C in the supercritical fluids related to the Variscan mineralizations in Erzgebirge, Slavkovsky Les (Kaiserwald), and Lusatian Mountains, Germany, and the Czech Republic and some remarks on the low-pressure formation of diamond. Geol Earth Mar Sci 5: 1-5.
  14. Thomas R (2026) Diamond in pycnite-rock from Altenberg and in the topaz from the famous Schneckenstein Geol Earth Mar Sci 8: 1-7.
  15. Thomas R (2025a) Extremely 13C-rich diamond in orthorhombic cassiterites in the Variscan Erzgebirge, Saxony/Germany. Earth Mar Sci 7: 1-5.