Fallstudie: Transversale Modeninstabilität in einem Kilowatt-Faserverstärker

Wir betrachten einen mantelgepumpten Faserverstärker mit einer Ytterbium-dotierten Faser mit großer Modenfläche. Diese hat einen Stufenindex-Kern mit 25 µm Durchmesser und einer numerischen Apertur von 0,06, einen Pumpmantel mit 400 µm Durchmesser und eine Länge von 4 m. Der Kern ist gleichmäßig mit Yb dotiert. Ein Signal mit 10 W bei 1064 nm wird mit 1 kW Pumpleistung bei 976 nm verstärkt, die in Vorwärtsrichtung eingekoppelt wird (Vorwärtspumpen). Verstärkte spontane Emission (ASE) wird berücksichtigt, erweist sich aber als schwach.

Bei 1064 nm führt der Faserkern (mit der niedrigen numerischen Apertur von 0,06) nicht nur die Grundmode LP₀₁, sondern auch die LP₁₁-Mode. Der Verstärker wird mit dem Signal in LP₀₁ betrieben, und dort möchten wir es auch behalten. Bei hohen Leistungen kann jedoch TMI Leistung in die LP₁₁-Mode übertragen. Dies verschlechtert die Strahlqualität und lässt das Strahlprofil am Ausgang fluktuieren.

Für die Simulationen verwenden wir die Software RP Fiber Power. Sie bietet eine Abschätzung der TMI-Schwelle auf der Basis einer linearen Stabilitätsanalyse des stationären Zustands des Verstärkers, die nur wenige Sekunden dauert. (Dieselbe Abschätzung ist auch in den Power Forms für kontinuierlich betriebene Faserverstärker und für Faserlaser mit Verstärkern verfügbar, dort für die Verstärkerstufen.) Zum Vergleich: Eine vollständige Simulation von TMI, die auch zeigt, was im Bereich starker TMI geschieht, könnte für eine einzige Konfiguration Stunden dauern.

Die Signal-Ausgangsleistung für 1 kW Pumpleistung beträgt 794 W. Die dissipierte Leistung ist 71 W, mit einer maximalen Wärmelast von 28 W/m nahe dem Pumpeingang. Der größte Teil der Wärme entsteht durch den Quantendefekt, d. h. die Energiedifferenz zwischen Pump- und Signalphotonen.

Der Mechanismus wird im Enzyklopädieartikel über transversale Modeninstabilität ausführlich erklärt (siehe auch [1, 5]); hier eine kurze Zusammenfassung:

• Eine schwache LP₁₁-Welle interferiert mit der starken LP₀₁-Welle. Da die beiden Moden unterschiedliche Ausbreitungskonstanten haben, bildet ihr Schwebungsmuster ein Intensitätsgitter entlang der Faser, dessen Periode die Schwebungslänge ist – hier 3 mm.
• Durch die Wärmeerzeugung im dotierten Kern und den thermooptischen Effekt erzeugt dieses Intensitätsgitter ein Gitter des Brechungsindex, das Licht von einer Mode in die andere streuen kann.
• Eine Netto-Leistungsübertragung erfordert eine Phasenverschiebung zwischen dem Brechungsindexgitter und dem Intensitätsmuster. Bei einem stationären Gitter ist das nicht möglich. Es geschieht jedoch, wenn die LP₁₁-Welle eine etwas niedrigere optische Frequenz hat (Versatz im kHz-Bereich): Das Intensitätsmuster bewegt sich dann langsam entlang der Faser, und das Temperaturgitter hinkt hinterher, weil die Wärmeleitung etwas Zeit benötigt. Diese Verzögerung erzeugt eine Phasenverschiebung, die eine Energieübertragung in die LP₁₁-Mode ermöglicht. Bei genügend hoher Pumpleistung führt dies zu einem exponentiellen Wachstum der LP₁₁-Leistung, sodass die Strahlqualität verdorben und der Betrieb instabil wird.

Das folgende Diagramm zeigt dieses Brechungsindexgitter für unseren Verstärker, als xz-Schnitt durch die Faser an der Position, wo die Kopplung am stärksten ist:

Man sieht, dass die Maxima des Brechungsindex tatsächlich um etwa 0,3 mm hinter den Intensitätsmaxima zurückbleiben – etwa ein Zehntel der Schwebungslänge. Diese kleine Verschiebung ist, wie oben erklärt, der Schlüssel zur Instabilität. Beachten Sie auch, dass das Temperaturgitter etwas in den Mantel hineinreicht und dass seine Fronten gekrümmt sind: Je weiter von der Faserachse entfernt, desto stärker hinken die Temperaturmaxima hinterher. Der Grund ist, dass die Wärme Zeit braucht, um nach außen zu diffundieren; am Rand des dargestellten Bereichs, etwa 6 µm außerhalb des Kerns, beträgt die Verzögerung bereits etwa ein Viertel der Schwebungslänge.

Die LP₁₁-Welle wächst entlang der Faser exponentiell an, mit einem Exponenten ($G$), der vom Frequenzversatz abhängt. Die Frage ist, von welchem Anfangsniveau aus sie wächst. Selbst bei perfekter Einkopplung und ohne technisches Rauschen liefert das Quantenrauschen einen minimalen LP₁₁-Seed, der mit einer Seed-Leistung in der Größenordnung von 10⁻¹⁶ W innerhalb der relevanten Bandbreite simuliert werden kann. Wenn wir die Schwelle als den Punkt definieren, an dem der LP₁₁-Anteil am Ausgang 1 % erreicht, ist der für 10 W Eingangsleistung erforderliche Exponent

Das folgende Diagramm zeigt den Exponenten als Funktion des Frequenzversatzes:

Für 1 kW Pumpleistung beträgt der maximale Exponent 25,7, deutlich unter 34,5. Indem wir die Pumpleistung erhöhen, finden wir die Schwelle bei 1,32 kW Pumpleistung, wo die Signal-Ausgangsleistung 1,04 kW beträgt. Der am stärksten wachsende Frequenzversatz ist dann 3,6 kHz. Diese Frequenz ist im Wesentlichen durch die thermische Diffusionszeit über den Bereich der Mode bestimmt – einen Bruchteil einer Millisekunde für den Abstand zwischen den beiden Maxima der LP₁₁-Mode.

Wenn nur das Quantenrauschen die Instabilität anstößt, können wir diesen Verstärker also mit 1 kW Pumpleistung und einer gewissen Sicherheitsmarge betreiben, aber nicht viel mehr als etwa 1 kW Ausgangsleistung erwarten. Diese Marge entspricht allerdings nur einem Faktor von etwa 7000 in der Seed-Leistung: Ein LP₁₁-Seed von etwa 10⁻¹² W nahe der TMI-Frequenz würde genügen, um die Schwelle bereits bei 1 kW Pumpleistung zu erreichen. Technisches Rauschen kann solche Seeds liefern (siehe unten).

Man könnte erwarten, dass TMI hauptsächlich dort entsteht, wo die Wärmelast am höchsten ist, oder nahe dem Ausgang, wo die Signalleistung am höchsten ist. Beides trifft nicht ganz zu:

Der lokale TMI-Verstärkungskoeffizient hat sein Maximum von etwa 20 m⁻¹ nur 0,4 m nach dem Eingang. Nahe dem Eingang ist er klein, weil die Kopplung proportional zur Signalleistung ist, die dort nur 10 W beträgt. Weiter stromabwärts ist die Signalleistung hoch, aber die TMI-Verstärkung nimmt wieder ab. Ein Grund dafür ist, dass die Wärmelast mit der Pumpleistung abnimmt. Wichtiger ist jedoch die Sättigung der Verstärkung: Für TMI zählt nicht die Wärmelast selbst, sondern wie stark sich die Wärmeerzeugung ändert, wenn die lokale Signalintensität durch das Schwebungsmuster moduliert wird. Bei starker Sättigung erhöht eine höhere Signalintensität die extrahierte Leistung kaum, sodass die Modulation der Wärmeerzeugung schwach wird [2].

Die Konsequenz für die LP₁₁-Leistung zeigt das nächste Diagramm:

An der Schwelle wächst die LP₁₁-Leistung von 10⁻¹⁶ W auf etwa 10 W am Ausgang – um 17 Größenordnungen, davon 13 innerhalb der ersten 1,5 m. Die restlichen 2,5 m, wo die Signalleistung am höchsten ist, tragen vergleichsweise wenig bei.

Ein großer Vorteil von Modellierung und Simulation ist, dass wir die Abhängigkeit von verschiedenen Parametern untersuchen können. Dies unterstützt die Optimierung von Geräten und ist weit schneller und günstiger, als all dies experimentell zu erkunden. Da unsere Abschätzung der Schwelle nur wenige Sekunden dauert, können wir leicht untersuchen, wie die TMI-Schwelle von verschiedenen Parametern abhängt.

Signal-Eingangsleistung

Eine höhere Signal-Eingangsleistung (von einer leistungsstärkeren Vorverstärkerstufe) erhöht die TMI-Schwelle: von 0,92 kW für 1 W Eingangsleistung auf 1,42 kW für 100 W. Dafür gibt es zwei Gründe:

• Ein stärkeres Eingangssignal sättigt die Verstärkung von Anfang an stärker, was die Modulation der Wärmeerzeugung im Eingangsbereich abschwächt, wo sich die Instabilität aufbaut.
• Das relative LP₁₁-Seed-Niveau ist niedriger, sodass ein etwas größerer Exponent erforderlich ist.

Die mögliche Verbesserung ist deutlich, wenn auch nicht dramatisch:

Pumprichtung

Mit Rückwärtspumpen (Pumplicht am Signalausgang eingekoppelt) steigt die Schwelle von 1,04 kW auf 1,22 kW. Hier wird der Eingangsteil der Faser, wo das Signal noch schwach ist, nur schwach gepumpt und erzeugt wenig Wärme; die meiste Wärme entsteht nahe dem Ausgangsende, wo die Verstärkung durch die hohe Signalleistung stark gesättigt ist. Interessanterweise wird die höchste Schwelle (1,26 kW) nicht für reines Rückwärtspumpen erreicht, sondern wenn 70 % bis 80 % der Pumpleistung in Rückwärtsrichtung eingekoppelt werden. Bei reinem Rückwärtspumpen ist die Wärme am Ausgangsende konzentriert, wo die Signalleistung – und damit die Kopplung – am höchsten ist.

Beachten Sie, dass das Rückwärtspumpen die maximale Wärmelast pro Längeneinheit erheblich erhöht (an der Schwelle 67 W/m statt 37 W/m), was eine effizientere Kühlung erfordern kann. Letztlich könnte man sich für beidseitiges Pumpen entscheiden, z. B. mit gleichen Leistungen von beiden Seiten.

Verluste höherer Moden

Zusätzliche Verluste für LP₁₁, wie sie durch Aufwickeln der Faser entstehen können (Biegeverluste), erhöhen die Schwelle – aber weniger, als man erwarten könnte: 10 dB/m zusätzlicher Verlust erhöhen sie nur um 25 %. Der Grund ist, dass der TMI-Verstärkungskoeffizient im Eingangsbereich weit höher ist (bis zu 20 m⁻¹) als ein solcher Verlust (2,3 m⁻¹).

Eine weitere interessante Richtung ist die Optimierung des Faserdesigns für erhöhte LP₁₁-Verluste.

Vermeidung von Hintergrundabsorption

Nun kommt eine weniger erfreuliche Überraschung: Hintergrundabsorption des Signallichts, wie sie durch Photodarkening verursacht werden kann, senkt die Schwelle drastisch. Bereits 0,01 dB/m senken sie von 1,04 kW auf 0,84 kW, und 0,1 dB/m auf 0,28 kW. Die absorbierte Leistung trägt zwar nur mäßig zur Wärmeerzeugung bei, aber diese Wärme folgt jeder Modulation der Signalintensität vollständig, während die Modulation der Wärme aus dem Quantendefekt durch die Sättigung der Verstärkung stark reduziert ist. Daher sind Fasern mit geringem Photodarkening für hohe TMI-Schwellen unerlässlich [4].

Man könnte meinen, dass es hilft, die Wärme über eine längere Faser zu verteilen – zum Beispiel mit einer Faser mit halber Dotierungskonzentration und doppelter Länge, sodass die Pumpabsorption gleich bleibt. Die maximale Wärmelast pro Längeneinheit wird dann halbiert. Das Modell liefert jedoch exakt dieselbe Schwelle.

Eine geringere Wärmelast pro Längeneinheit ist also für sich genommen kein Hebel – zumindest solange es keine konzentrationsabhängigen Prozesse (z. B. Quenching) und keine Hintergrundabsorption gibt, die kürzere Fasern begünstigen würden.

Absolute Werte der Schwelle aus solchen Modellen sollten als Abschätzungen betrachtet werden. Vergleiche mit veröffentlichten Modellrechnungen zeigen Abweichungen von bis zu etwa 20 %, und die Eingabewerte sind oft nicht gut bekannt:

• Seed-Niveau: Der relevante Seed ist der Teil des LP₁₁-Feldes, der relativ zu LP₀₁ im kHz-Bereich fluktuiert; das Quantenrauschen ist nur eine Untergrenze. Glücklicherweise geht das Seed-Niveau nur logarithmisch ein: Eine zehnfache Seed-Leistung senkt die Schwelle um etwa 6 %. Dennoch würde ein Seed von 10⁻⁸ W statt 10⁻¹⁶ W die Schwelle halbieren. Solch starke Seeds können in der Praxis auftreten: Intensitätsrauschen der Pumpe im kHz-Bereich moduliert die Wärmelast und damit die Phasenbeziehung der Moden; zusammen mit einem statischen LP₁₁-Anteil aufgrund einer unvollkommenen Einkopplung erzeugt es fluktuierendes LP₁₁-Licht bei den relevanten Frequenzen. In Experimenten senkten wenige Prozent Intensitätsrauschen der Pumpe die Schwelle um mehr als die Hälfte [6]. Rauscharme Treiber für die Pumpdioden und eine saubere Einkopplung sind daher wichtig. Dabei zählen nur Rauschanteile nahe der TMI-Frequenz (hier einige kHz), nicht das gesamte Intensitätsrauschen. Da die Seed-Leistung mit dem Quadrat der Rauschamplitude skaliert, ändert jeder Faktor 10 in dieser Rauschamplitude die Schwelle um grob 13 %.
• Brechungsindexprofil: Die Brechungsindexprofile kommerziell erhältlicher Fasern sind den Anwendern oft nicht bekannt. Simulationen zeigen jedoch, dass typische Abweichungen vom Stufenindexprofil, etwa ein schmaler zentraler Einbruch des Brechungsindex oder Welligkeiten, die Schwelle um weniger als 10 % verändern.
• Verluste höherer Moden und Photodarkening: Wie oben gezeigt, können diese erhebliche Auswirkungen haben, sind aber oft nicht genau bekannt. Photodarkening bedeutet mit der Zeit ansteigende Verluste – und deren Stärke kann stark von Betriebsparametern abhängen, etwa vom Grad der Yb-Anregung. Solche Parameter können ebenfalls mithilfe der Modellierung überwacht und optimiert werden.
• Signal-Linienbreite: Das Modell nimmt an, dass das Signal über die Laufzeitdifferenz der beiden Moden kohärent ist. Für diese Faser gilt das für Linienbreiten unter etwa 200 GHz – weit mehr als die wenigen GHz, die typischerweise zur Unterdrückung der Brillouin-Streuung verwendet werden. Bei sehr großen Signal-Linienbreiten kann TMI reduziert werden, weil das Interferenzmuster verwaschen wird.

Generell sind Trends, wie sie sich beim Vergleich von Design-Optionen ergeben, zuverlässiger als absolute Werte. So stimmt zum Beispiel die Abhängigkeit von der Pumpwellenlänge gut mit experimentellen Befunden überein, wie oben gezeigt. In der Praxis kann eine gemessene Schwelle eines existierenden Verstärkers verwendet werden, um unsichere Parameter zu kalibrieren; das Modell sagt dann die Auswirkungen von Design-Änderungen voraus.

Die Software RP Fiber Power ist ein prima Tool für solche Arbeiten – sehr leistungsfähig und doch einfach zu bedienen!

Sie können von dieser Studie einiges lernen:

• Die TMI-Schwelle eines Kilowatt-Faserverstärkers lässt sich innerhalb von Sekunden abschätzen; das Ergebnis ist eine Abschätzung, eignet sich aber gut für den Vergleich von Design-Optionen.
• Die Instabilität baut sich hauptsächlich in einem relativ kurzen Teil der Faser nahe dem Eingang auf, wo das Signal schon recht stark, die Verstärkung aber noch nicht stark gesättigt ist.
• Die Sättigung der Verstärkung ist der Schlüssel zum Verständnis vieler Abhängigkeiten: Eine höhere Signal-Eingangsleistung, (teilweises) Rückwärtspumpen und Pumpen abseits des Absorptionsmaximums erhöhen die Schwelle allesamt, indem sie die Sättigung verstärken.
• Hintergrundabsorption, z. B. durch Photodarkening, kann die Schwelle drastisch senken, da ihre Wärme nicht durch Sättigung reduziert wird.
• Die Wärme bei gleicher Pumpabsorption über eine längere Faser mit niedrigerer Dotierung zu verteilen, hilft nicht.
• Die berechnete Schwelle setzt voraus, dass die Instabilität nur durch Quantenrauschen angestoßen wird. Intensitätsrauschen der Pumpe im kHz-Bereich kann sie zusammen mit einem kleinen statischen LP₁₁-Anteil aufgrund der Einkopplung erheblich senken.

(Suggest additional literature!)