OFW-Gassensorik · Diss. 2002
Dissertation · RWTH Aachen · 21. Juni 2002

Oberflächenwellen​-Bauelemente als Gassensoren

Untersuchungen zur Modellierung und Optimierung von Oberflächenwellen-Bauelementen und Auswerteschaltungen für die Gassensorik — mit einem erweiterten Coupling-of-Modes-Modell, das erstmals den Einfluss einer gassensitiven Polymerbeschichtung auf das komplette Übertragungsverhalten berechenbar macht.

Verfasser
Dipl.-Ing. Peter Oechtering
Fakultät
Elektrotechnik & Informationstechnik
Berichter
Prof. Dr. W. Mokwa · Prof. Dr. D. Kohl
Betreuung
Prof. Dr. P. Kaul
0 MHz Arbeitsfrequenz der Quarz-Resonatoren auf ST-X-Schnitt
0 Resonanzgüte Q des unbeschichteten Resonators
0 µm Elektrodenperiode 2d = λ der Interdigitalwandler
0 ppm Toluol-Testkonzentration an der Gasmischeranlage
Kurzfassung

Worum es geht

Lagern sich Gasmoleküle an einer gassensitiven Polymerschicht an, so nimmt deren Masse zu — und je nach Beschichtung ändern sich zusätzlich die elektrischen und mechanischen Randbedingungen an der Substratoberfläche. Beides beeinflusst Phasengeschwindigkeit und Dämpfung einer piezoelektrisch angeregten Oberflächenwelle. Damit ändert sich das elektrische Übertragungsverhalten des Bauelements, das messtechnisch sehr genau bestimmt werden kann.

Ziel der Arbeit ist ein Modell für OFW-Gassensoren, das eine Simulation des Sensors gemeinsam mit der Auswerteelektronik erlaubt. Der Schwerpunkt liegt auf der Charakterisierung und Modellierung von Oberflächenwellen-Resonatoren; zum Vergleich werden Verzögerungsleitungen und Dickenscherschwinger untersucht.

Der eigentliche Beitrag ist das erweiterte Coupling-of-Modes-Modell: Es verbindet die aus der Literatur bekannten Modelle für den Einfluss viskoelastischer Schichten mit einem COM-Ansatz, der Substrateigenschaften und Elektrodengeometrie berücksichtigt. Ergebnis ist das komplette Übertragungsverhalten in Form von Streuparametern — und damit die Möglichkeit, Sensor und Oszillatorschaltung applikationsspezifisch zu optimieren.

Die Kernaussagen in Kürze

  • Die Beschichtungstechnik dominiert das Sensorverhalten stärker als das Beschichtungsmaterial: Pyrolyse und Sublimation liefern homogene Schichten, Drop- und Spray-Coating streuen stark.
  • Eine Temperaturkompensation über ein unbeschichtetes Referenzelement ist nicht sinnvoll — beschichtete Sensoren zeigen einen anderen, teils im Vorzeichen umgekehrten Temperaturgang.
  • Die Phasenlage, bei der ein Sensor im Oszillator betrieben wird, verändert die Sensitivität um mehr als 50 %. Weil Luftfeuchte und Toluol unterschiedlich stark davon abhängen, lassen sie sich durch Modulation der Phasenlage mit einem einzelnen Sensor unterscheiden.
  • Bei niedrigviskosen Schichten erreicht der Resonator die niedrigere Nachweisgrenze, bei hochviskosen und dicken Schichten die Verzögerungsleitung.
  • Ein gezielt „verstimmter“ Resonator, der nur im Resonanzraum beschichtet wird, erreicht in der Simulation eine bis zu 2-fach niedrigere Nachweisgrenze.

Zu dieser Webseite. Sie gibt den Inhalt der Dissertation strukturiert wieder und ergänzt ihn um interaktive 3D-Modelle und Diagramme. Alle Kurven in den Diagrammen werden live aus den Gleichungen der Arbeit gerechnet (die jeweilige Gleichungsnummer steht dabei); Messwerte stammen aus den Tabellen 2.1–5.2 und sind als solche gekennzeichnet. Die 3D-Darstellungen sind physikalisch motivierte Visualisierungen der Abbildungen 2.2, 2.5–2.13 und 2.7–2.9 — keine Finite-Elemente-Rechnungen.

Kapitel 1

Einleitung — warum künstliche Nasen schwierig sind

Der menschliche Geruchssinn ist ein außerordentlich komplexes biologisches System: Über 100 Millionen Proteinrezeptoren erlauben es, chemische Verbindungen schon in sehr geringer Konzentration wahrzunehmen und kleine Unterschiede in der Zusammensetzung komplexer Stoffgemische zu erkennen. Dennoch kann der Mensch viele toxische Substanzen und Umweltgifte gar nicht riechen.

Genau dort setzen chemische Sensoren an: Onlinemessung und Überwachung maximal zulässiger Schadstoffkonzentrationen in der Umwelt, am Arbeitsplatz oder zuhause. Die Beispiele reichen von der flächendeckenden Ozonmessung über die Überwachung von Benzolkonzentrationen an Tankstellen, Dioxin bei Müllverbrennungsanlagen, Kohlendioxid in der Lebensmittelindustrie bis zu narkotischen Gasen in Krankenhäusern — und weiter zur Detektion von Sprengstoffen an Flughäfen oder dem Aufspüren von Landminen.

Auch die Industrie nutzt Gassensoren zur Qualitätsüberwachung. In der Automobilindustrie wird häufig speziell geschultes Personal eingesetzt, um Ausdünstungen von Kunststoffen in Neufahrzeugen zu finden; bei Kaffee, Wein, Bier oder Käse ist der menschliche Geruchssinn bislang unverzichtbar. Chemische Sensoren sind hier wegen zu geringer Sensitivität und Selektivität nur bedingt geeignet — obwohl sie deutlich günstiger wären und eine kontinuierliche, objektive Überwachung ohne Erkältung oder Rauchgewöhnung erlauben würden.

Das Messprinzip dieser Arbeit

Die hier untersuchten Sensoren beruhen auf der Beeinflussung akustischer Wellen, die in piezoelektrischen Festkörpern elektrisch angeregt werden. Die Bauelementoberfläche wird mit gassensitiven Polymeren beschichtet. Lagern sich Gasmoleküle an, ändern sich Ausbreitungsgeschwindigkeit und Amplitude der akustischen Welle. In einer Oszillatorschaltung werden diese Größen in elektrische Signale umgesetzt.

Beim Dickenscherschwinger (BAW, QCM) sind Änderungen der elastischen oder elektrischen Schichteigenschaften gegenüber der Massenbelegung vernachlässigbar — er ist im Wortsinn ein Massensensor. Beim Oberflächenwellen-Sensor dagegen müssen zusätzlich Leitfähigkeit, Dielektrizität, Elastizität und Viskosität der Schicht berücksichtigt werden. Diese hohe Empfindlichkeit gegenüber allen Schichteigenschaften ermöglicht sehr niedrige Nachweisgrenzen — sie ist aber zugleich das Problem: Ohne genaue Kenntnis der Schichteigenschaften und ihrer Abhängigkeit von Gaskonzentration, Temperatur und Luftfeuchte lässt sich das Sensorverhalten nicht vorhersagen.

Und genau deshalb braucht es Modelle. Nur wenn Sensitivität und Störgrößen berechenbar sind, lässt sich ein Sensorsystem optimal auf eine Anwendung auslegen — und ein Sensorarray so entwerfen, dass die Einzelsensoren möglichst wenig redundante Information liefern.

Kapitel 2.1

Ausbreitung von Oberflächenwellen

Zur Erzeugung akustischer Wellen muss der Festkörper durch eine mechanische Spannung deformiert werden — bei piezoelektrischen Materialien geschieht das durch ein elektrisches Feld.

Im eindimensionalen Fall gilt der einfache Zusammenhang S = Δu/Δx zwischen Dehnung und Auslenkung. Im anisotropen Festkörper beeinflussen jedoch auch Kräfte senkrecht zur betrachteten Richtung die Dehnung; sie muss deshalb als Dehnungstensor geschrieben werden:

Skl = ½ ( ∂uk/∂xl + ∂ul/∂xk )   mit k, l = 1, 2, 3 Dehnungstensor. Über den Elastizitätstensor cijkl ist er mit der mechanischen Spannung Tij verknüpft. Von den 81 Elementen des Elastizitätstensors sind wegen der Symmetriebedingungen nur 21 voneinander unabhängig.
Gl. 2.1

Piezoelektrizität

Voraussetzung für piezoelektrisches Verhalten eines Kristalls ist eine polare Achse — eine Achse, bezüglich derer die Strukturzelle nicht faltungssymmetrisch ist. Drückt man in Richtung dieser Achse auf die Quarz-Strukturzelle, verschieben sich die positiv geladenen Si-Ionen gegenüber den negativ geladenen Sauerstoffionen des SiO₂-Kristalls: Es entsteht eine Polarisierung (direkter Piezoeffekt). Der dazu reziproke inverse Effekt — Deformation beim Anlegen eines elektrischen Feldes — wird zur Anregung der Wellen genutzt.

3D-Modell · Quarz-Strukturzelle und direkter Piezoeffekt interaktiv Abb. 2.2
0 %
Si⁴⁺ O²⁻ Elektroden
Ziehen zum Drehen, Mausrad zum Zoomen. Die Stauchung entlang der polaren Achse verschiebt Si- und Sauerstoffionen gegeneinander — an den Elektroden entsteht eine messbare Spannung. Den elastischen Rückstellkräften stehen dabei die Coulombkräfte gegenüber.

Für ideal elastische piezoelektrische Materialien erweitern sich die Materialgleichungen um den Tensor der piezoelektrischen Kopplungskoeffizienten eijk:

Tij = Σk,l cijkl·Skl − Σk ekij·Ek    ·    Di = Σj εij·Ej + Σj,k eijk·Sjk Gekoppeltes Gleichungssystem für mechanische Spannung T und elektrische Verschiebungsdichte D. Setzt man es in die Bewegungsgleichung ein und drückt das elektrische Feld als Gradient des Potentials Φ aus, erhält man ein lineares System gekoppelter Differentialgleichungen (Gl. 2.15 und 2.18), mit dem sich prinzipiell alle ausbreitungsfähigen Wellenmodi beschreiben lassen.
Gl. 2.7 / 2.8

Piezoelektrische Substrate

Die beiden wichtigsten Auswahlkriterien sind die piezoelektrische Kopplung KP² — sie gibt an, welcher Anteil der akustischen Wellenenergie als elektrische Energie vorliegt, und damit die Effizienz der Wandlung — und der Temperaturgang der Frequenz.

KP² = 2 · | Δv / v₀ |elektr. = 2 · | (v₀ − vE) / v₀ | Definition über die durch eine leitfähige Oberflächenbeschichtung verursachte Änderung der Phasengeschwindigkeit — dadurch messtechnisch bestimmbar. vE ist die Phasengeschwindigkeit bei elektrisch kurzgeschlossener Oberfläche.
Gl. 2.9
Tabelle 2.1 Materialparameter verschiedener piezoelektrischer Substrate
Substrat / SchnittKP² v₀ [m/s]Temperaturkoeffizient TKfRolle in dieser Arbeit
Quarz ST-X0,00163158 parabolisch, ≈ 0 bei 25 °CSubstrat der 433,9-MHz-Resonatoren
LiNbO₃ YZ0,0453488 linear, −94 ppm/Khohe Kopplung, aber starker Temperaturgang
LiTaO₃ X, 112° rot-Y0,00753300 linear, −18 ppm/KVergleichsmaterial
ZnO Z-X (gesputtert)0,024 … 0,0322630 … 2811 linear, −15 ppm/K (gemessen −22 … −37)Dünnschicht-Verzögerungsleitungen

Quarz hat zwar die niedrigste piezoelektrische Kopplung, wird aber wegen des bei Raumtemperatur nahezu verschwindenden Temperaturgangs bevorzugt eingesetzt — für die Gassensorik ist das entscheidend, weil Temperaturschwankungen hier die stärkste Störgröße darstellen. ZnO hat den strategischen Vorteil, als piezoelektrische Dünnschicht auf Silizium aufgesputtert werden zu können: Damit lassen sich Sensor, Temperaturmesswiderstände, Heizstrukturen und prinzipiell die komplette Auswerteelektronik auf einem Substrat integrieren.

Rayleighwellen

Eine mögliche Lösung des Gleichungssystems ist die nach Lord Rayleigh benannte Oberflächenwelle — der dominierende Mode der meisten kommerziell erhältlichen OFW-Bauelemente. Sie lässt sich anschaulich als Kombination einer Longitudinalwelle und einer vertikal polarisierten Scherwelle auffassen: Die Partikel nahe der Oberfläche folgen einer elliptischen, retrograden Bahn. Bei einer idealen Rayleighwelle gibt es keine Auslenkung in x-Richtung, und die Schwingungsamplitude klingt von der Oberfläche aus exponentiell in die Tiefe ab.

3D-Modell · Rayleighwelle im Substrat interaktiv Abb. 2.6
Die gelbe Ellipse zeigt die Bahn eines Oberflächenteilchens — sie wird gegen die Ausbreitungsrichtung durchlaufen (retrograd). Beachten Sie den exponentiellen Abfall der Amplitude mit der Tiefe: Die insgesamt bewegte Masse ist sehr gering, die Eindringtiefe liegt unter einer Wellenlänge. Genau daher rührt die hohe Massenempfindlichkeit des Sensors.
ui = Ai · eγRκiy · ej(ωt − γRz)    mit    γR = α + jk ,   k₀ = ω / v₀ Lösungsansatz als Linearkombination von vier Partialwellen. Der Dämpfungsbelag α beschreibt die Amplitudenabnahme längs der Laufstrecke (Einheit m⁻¹); bei verlustloser Ausbreitung wäre α = 0. Wegen der komplizierten Struktur der Gleichungen erfolgt die Berechnung in der Regel iterativ, nicht analytisch.
Gl. 2.19 / 2.21 / 2.22
Kapitel 2.2

Oberflächenwellen-Bauelemente

Zur Anregung wird ein Array gegenüberliegender, ineinandergreifender Fingerstrukturen verwendet: der Interdigitalwandler (IDT). Die akustische Wellenlänge ist ausschließlich durch den geometrischen Elektrodenabstand bestimmt.

f0 = v₀ / λ = v₀ / (2d) Für die in dieser Arbeit verwendeten Bauelemente auf Quarz (v₀ = 3158 m/s) und eine Resonanzfrequenz von 433 MHz folgt ein Elektrodenabstand von 2d = 7,29 µm. Angeregt werden können nur ungeradzahlige Vielfache dieser Frequenz — nur dann überlagern sich elektrisches Feld und Oberflächenwelle gleichphasig.
Gl. 2.23
3D-Modell · Zweitor-Resonator und Verzögerungsleitung interaktiv Abb. 2.7 – 2.9
Resonator: Reflektorstrukturen an beiden Substratenden werfen die zur Kante laufenden Wellen zurück; es bildet sich eine stehende Welle mit maximaler Energiedichte in der Substratmitte („Cavity“). Einfügungsdämpfung typisch 10 dB. Verzögerungsleitung: Keine Reflektoren — die zur Kante laufende Welle geht verloren (Kanten angeschrägt bzw. akustische Absorber). Selbst bei vollständiger Wandlung wird deshalb nur die Hälfte der Energie genutzt; die Einfügungsdämpfung liegt deutlich über 20 dB. Schalten Sie die Beschichtung zu, um die zusätzliche Bedämpfung zu sehen.
Darstellung schematisch: Die Elektrodenzahl ist stark reduziert — real besitzt jeder Wandler N = 35 Elektrodenpaare, jeder Reflektor M = 220 Elektroden, und der Wandlerabstand beträgt 180 µm ≙ 25 λ.
Resonator

Kleine Fläche, hohe Güte — empfindlich gegen Dämpfung

Da nahezu die ganze Substratfläche mit Elektroden belegt ist, sind die Resonanzpeaks wesentlich schmalbandiger als bei der Verzögerungsleitung. Die Welle wird zwischen den Reflektoren mehrfach reflektiert; man kann eine virtuelle Laufzeit definieren, die weit größer ist als der geometrische Wandlerabstand.

Der Preis: Nimmt die Dämpfung durch eine Beschichtung zu, wird die Welle seltener reflektiert, die virtuelle Laufzeit sinkt — und die Güte bricht überproportional ein.

Verzögerungsleitung

Robuste Güte, aber hohe Grunddämpfung

Die Phasenverschiebung ergibt sich aus dem geometrisch festgelegten Wandlerabstand (φ = D/v · 2πf). Die Güte hängt daher kaum von der Beschichtung ab. Zusätzlich ist die große unmetallisierte Fläche zwischen den Wandlern wesentlich einfacher und reproduzierbarer zu beschichten.

Rechnungen in der Literatur zeigen: Ab einer beschichtungsbedingten Zusatzdämpfung von etwa 10 dB hat die Verzögerungsleitung die niedrigere Nachweisgrenze. Da sie in dieser Arbeit meist darunter blieb, wurden überwiegend Zweitor-Resonatoren eingesetzt.

Alternative resonante Sensoren

Neben der Rayleighwelle lassen sich weitere Wellentypen sensorisch nutzen. Der Vergleich zeigt, warum die Rayleighwelle für Gassensoren so empfindlich, für Flüssigkeiten aber ungeeignet ist.

3D-Modell · Wellentypen im Vergleich interaktiv Abb. 2.5 · 2.6 · 2.10 – 2.12
Kapitel 2.3

Wie die Beschichtung die Welle beeinflusst

Formal lässt sich die Änderung der Phasengeschwindigkeit als totales Differential über alle Einflussgrößen schreiben — Massenbelegung m, elastische Konstante c, dielektrische Konstante ε, Oberflächenleitfähigkeit σ, Temperatur T und Druck p:

Δvv₀ = 1v₀ [ ∂v/∂m · Δm + ∂v/∂c · Δc + ∂v/∂ε · Δε + ∂v/∂σ · Δσ + ∂v/∂T · ΔT + ∂v/∂p · Δp ] Nach Ricco et al. Für einige dieser Größen lassen sich mit Näherungsmethoden geschlossene Gleichungen finden — die folgenden Abschnitte behandeln sie einzeln.
Gl. 2.25

2.3.1 Ideal elastische, nichtleitfähige Beschichtung

Aus einer Störungsanalyse folgt die von Wohltjen angegebene Gleichung. Sie ist die in der Gassensorik am häufigsten verwendete Beziehung überhaupt:

Δvv₀ = (k₁ + k₂) · f₀ · ρf · hf  −  k₂ · f₀ · hf · ff + µf)v₀² (λf + 2µf) Erster Term: Massenzunahme — weil k₁, k₂ negativ sind, verlangsamt zusätzliche Masse die Welle. Zweiter Term: Versteifung durch das Schermodul µf. Bei weichen Polymeren (µf ≈ 10⁶…10⁷ Pa) liegt sein Beitrag unter 1 %, bei kristallinen Polymeren (10⁹ Pa) bei ca. 10 %, bei einer Kupferphtalocyanin-Dünnschicht (≈ 7·10⁹ Pa) bei rund 60 %.
Gl. 2.31
Tabelle 2.2 Materialkonstanten für die Wohltjen-Gleichung
Substratv₀ [m/s]k₁ [m²s/kg] k₂ [m²s/kg]k₁+k₂ [m²s/kg]
ST-X Quarz3158−8,7·10⁻⁸ −3,9·10⁻⁸−12,6·10⁻⁸
ZnO Z-X und Y2639−5,47·10⁻⁸ −2,06·10⁻⁸−7,53·10⁻⁸

Für den Dickenscherschwinger gilt stattdessen die klassische Sauerbrey-Gleichung, die ausschließlich die Massenzunahme berücksichtigt — mit kBAW = −2/(ρQ·vQ) = −22,6·10⁻⁸ m²s/kg (ρQ = 2650 kg/m³, vQ = 3340 m/s für AT-Quarz).

ΔfBAWf0 BAW = kBAW · f₀ · ρf · hf    ·    ΔfOFWf0 OFW = (k₁ + k₂) · f₀ · ρf · hf Direkt vergleichbar: Bei gleicher Resonanzfrequenz ist die Massensensitivität des Dickenscherschwingers um den Faktor 1,8 höher als die des Rayleighwellen-Sensors — Folge der unterschiedlichen Partikelauslenkungen an der Oberfläche. Entscheidend ist aber die erreichbare Frequenz: Beim BAW bestimmt die Substratdicke die Frequenz (d = 1660/f₀ mm — bei 30 MHz nur noch 55 µm und damit zu zerbrechlich), beim OFW allein die Lithographie.
Gl. 2.35 / 2.37
Diagramm · Massensensitivität OFW vs. Dickenscherschwinger gerechnet · Gl. 2.35 / 2.37
Absolute Frequenzänderung |Δf| = f₀²·k·ρf·hf — sie wächst quadratisch mit der Resonanzfrequenz. Da das Oszillatorrauschen jedoch linear mit der Frequenz zunimmt, ist für den Sensitivitätsvergleich die relative Änderung Δf/f₀ die sinnvollere Größe; sie steigt linear.

2.3.2 Viskoelastische Beschichtung — akustisch dünn und akustisch dick

Bei Polymerbeschichtungen beobachtet man neben der Frequenzabnahme auch eine Zunahme der Einfügungsdämpfung. Zur Beschreibung dieser Verluste werden das komplexe Kompressionsmodul K = K′ + jK″ und das Schermodul G = G′ + jG″ verwendet; die Imaginärteile repräsentieren die um 90° phasenverschobenen Spannungsanteile und damit die Verluste.

Entscheidend ist, ob sich der Film synchron mit der Substratoberfläche bewegt. Als Kriterium dient die Kennzahl

RA = A · f₀ · v₀ · ρf · hf|G|    mit    A = 2π · C₂C₁ + 2C₂   →   AQuarz = 1,9 RA ≪ 1 → akustisch dünn (Gradienten über die Schichtdicke vernachlässigbar). RA ≥ 1 → akustisch dick. Für die 433-MHz-Quarz-Resonatoren liegt die Grenze bei 0,38 nm für weiche Polymere (|G| = 10⁶ Pa) und bei 380 nm für kristalline (10⁹ Pa) — bei weichen Polymeren müssen also selbst Schichten weit unterhalb der Wellenlänge (λR = 7,3 µm) als akustisch dick behandelt werden.
Gl. 2.40 / 2.41
3D-Modell · Deformation des Polymerfilms durch die Rayleighwelle interaktiv Abb. 2.13
Substrat Film synchron Scherung im Film

Für den akustisch dünnen Film ergeben sich lineare Zusammenhänge — die Phasengeschwindigkeit hängt vom Realteil des Schermoduls ab, die Dämpfung vom Imaginärteil:

Δvv₀ ≅ (k₁+k₂)·f₀·ρf·hf − k₂·f₀·hf· 4G′v₀²    ·    Δαk₀ ≅ − k₂·f₀·hf·4G″v₀² Akustisch dünner Film. Vereinfachung Ezz ≈ 4G, weil bei den meisten Polymeren |K| ≫ |G| gilt.
Gl. 2.44 / 2.45
Δvv₀ + j Δαk₀ = k₁·f₀·ρf·hf + k₂·f₀ · βRz·Gj·ω² · tanh( j·βRz·hf )   ,   βRz = ω f − 4G/v₀²) / G Akustisch dicker Film. Der Term βRz·hf bestimmt die Phasenverschiebung zwischen Unter- und Oberseite des Films. Erreicht sie π/2, tritt Filmresonanz auf: Die in den Film eingekoppelte akustische Energie wird maximal, die Oberflächenwelle stark bedämpft. Oberhalb der Filmresonanz wechselt sogar das Vorzeichen der Geschwindigkeitsänderung — eine Massenzunahme kann dann theoretisch zu einer Erhöhung der Phasengeschwindigkeit führen. Für kleine Argumente geht die Gleichung in Gl. 2.44/2.45 über.
Gl. 2.51
Diagramm · Phasengeschwindigkeit und Dämpfung über der Schichtdicke gerechnet · Gl. 2.44 / 2.45 / 2.51
Stellen Sie G′ = 107 Pa ein: Die Filmresonanz liegt bei rund 60 nm — genau wie in Abbildung 2.15 der Arbeit. Mit G′ = 108 Pa verschiebt sie sich auf etwa 190 nm (Abbildung 2.16), die Grenze RA = 1 auf 38 nm. Die gestrichelte Kurve zeigt, ab wann das einfache Modell für akustisch dünne Filme unbrauchbar wird.
Hinweis: Die absoluten Werte folgen den Konstanten aus Tabelle 2.2 und stimmen mit Tabelle 3.3 überein (23,7 kHz/nm für ρf = 1 g/cm³ bei 433,9 MHz); die Ordinatenskalierung der Abbildungen 2.15/2.16 im Original weicht davon ab.

2.3.3 Elektrisch leitfähige Beschichtung

Jede Oberflächenwelle auf piezoelektrischem Substrat wird von einem elektrischen Feld begleitet. Eine leitfähige oder halbleitende Schicht wechselwirkt damit und trägt Oberflächenladungen — bei sehr guter Leitfähigkeit als rein kapazitive Kopplung, sonst zusätzlich mit resistivem Anteil.

Δvv₀ ≅ − KP²2 · 11 + (v₀·COb.Ob.   ·   Δαk₀KP²2 · (v₀·COb.Ob.1 + (v₀·COb.Ob. Die stärkste Änderung und die maximale Dämpfung treten bei σOb. = v₀·COb. auf. Der Kapazitätsbelag ergibt sich als Summe der Dielektrizitätszahlen von Substrat und Film (Gl. 2.54); für Quarz beträgt er ca. 50 pF/m, vollständige Wasserbedeckung entspräche ≈ 760 pF/m.
Gl. 2.52 / 2.53
Diagramm · Einfluss der Oberflächenleitfähigkeit gerechnet · Gl. 2.52 / 2.53
In der halblogarithmischen Darstellung ist gut zu erkennen: Bei sehr kleinen und sehr hohen Leitfähigkeiten haben Änderungen kaum Auswirkungen. Im mittleren Bereich (10⁻⁸ … 10⁻⁶ S) fällt die Phasengeschwindigkeit stark ab, die Dämpfung wird maximal. Genau deshalb eignen sich halbleitende Schichten wie Phtalocyanine — sie ändern ihre Leitfähigkeit bei Adsorption oxidierender Gase (NO, NO₂, SO₂, Ozon).

2.3.4 Temperatureinfluss

Bei ST-X-Quarz ist der Temperaturgang des Substrats nahezu vernachlässigbar — der dominierende Effekt sind daher die Materialeigenschaften der Polymerschicht. Oberhalb der dynamischen Glastemperatur werden die Kettensegmente beweglich; Scher- und Kompressionsmodul ändern sich sprunghaft.

Die dynamische Glastemperatur hängt stark von der Anregungsfrequenz ab und liegt deutlich über der in Torsionspendelmessungen bestimmten statischen. Literaturwerte für G und K sind deshalb nur eingeschränkt verwendbar — umgekehrt sind polymerbeschichtete OFW-Bauelemente eine gute Methode, diese Parameter bei hohen Frequenzen zu messen.

Praktische Konsequenz: Man sollte Polymere wählen, deren dynamische Glastemperatur für die gewählte Anregungsfrequenz deutlich unter der Betriebstemperatur liegt. Dann ist das Polymer weich, im thermodynamischen Gleichgewicht, und Temperaturschwankungen wirken sich schwächer aus. Zusätzlich diffundieren Gasmoleküle schneller ein — die Ansprechzeit wird besser.

2.3.5 Adsorption und Absorption

Adsorption — an der Oberfläche

Im Gleichgewicht ist die adsorbierte Konzentration proportional zur Konzentration im Trägergas; der Partitionskoeffizient KAds = eEa/RT hängt exponentiell von der Bindungsenergie ab.

Die Bindungsenergie sollte hoch sein für hohe Sensitivität — aber die mittlere Verweilzeit steigt mit dem gleichen Faktor. Physisorption (Van-der-Waals, 0,5…5 kJ/mol) ist bei Raumtemperatur reversibel, Chemisorption (80…400 kJ/mol) nicht; sie ist dafür selektiver.

Absorption — im Schichtvolumen

Bei Polymerschichten dominierend. Der Partitionskoeffizient lässt sich über die LSER-Theorie (Linear Solvation Energy Relationship) berechnen:

log KAbs = c + r·R₂ + s·πH₂ + a·αH₂ + b·βH₂ + l·log L16

Die Faktoren stehen für Polarisierbarkeit, Dipolarität, Acidität und Basizität von Wasserstoffbrücken sowie für dispersive Kräfte. Ein so berechneter Koeffizient erlaubt eine Abschätzung der Sensitivität.

Der Effekt, der oft übersehen wird. Absorbierte Lösungsmitteldämpfe weichen das Polymer auf — das Schermodul sinkt. Die daraus folgende Frequenzänderung hat dasselbe Vorzeichen wie die durch die Massenzunahme und addiert sich. Bei Polyisobutylen ließ sich die beobachtete Geschwindigkeitsabnahme bei Trichlorethylen nur zu 67 %, bei Pentan sogar nur zu 40 % auf die Masse zurückführen; bei manchen Polymeren wurden vier- bis sechsmal höhere Signale gemessen als eine reine Massenzunahme erklären kann. Umgekehrt kann ein quellendes, gut haftendes Polymer durch die Grenzflächenspannung versteifen — dieser Effekt wirkt der Frequenzabnahme entgegen. Beide Mechanismen haben unterschiedliche Zeitkonstanten und sind im Signalverlauf unterscheidbar.

Kapitel 2.4

Modelle für Oberflächenwellen-Elemente

Ersatzschaltbild nach Mason

Ein infinitesimal kurzes Stück des Bauelements wird als Reihenschaltung von Induktivitäten mit parallelen Leitwerten und Kapazitäten beschrieben — Leitungstheorie in Analogie zur elektromagnetischen Welle. Betrachtet wird eine „Einheitszelle“ der Länge d = λ/2 mit einer Elektrode in der Mitte.

Sie hat drei Tore: zwei akustische (U₁,I₁ · U₂,I₂) und ein elektrisches (U₃,I₃). Die Transformation der akustischen Größen erfolgt über das Übersetzungsverhältnis Φ: U = F/Φ und I = v·Φ. Damit lässt sich eine charakteristische Impedanz definieren:

Z₀ = Zmech/Φ² = 2π / (ω₀ Cs KP²) Für ein Metallisierungsverhältnis η = 0,5 vereinfacht sich die statische Kapazität eines Elektrodenpaars zu Cs = W · εeff, also direkt proportional zur Überlappungsweite W.
Gl. 2.61 / 2.66

Coupling-of-Modes (COM)

Der COM-Ansatz beschreibt Systeme mit räumlich verteilten Parametern und erlaubt geschlossene Lösungen für komplette Interdigitalwandler und Reflektorstrukturen. Angenommen wird, dass ein Wandler vorwärts- und rückwärtslaufende Wellen a(z), b(z) mit langsam veränderlichen Amplituden R(z), S(z) erzeugt.

Weil die Phasengeschwindigkeit im metallisierten Bereich niedriger ist, wird die Welle an den Elektrodenkanten teilweise reflektiert — beide Wellen koppeln:

dR/dz = −j k₁₁ R(z) − j k₁₂ e2jδz S(z)
dS/dz =   j k₁₂* e−2jδz R(z) + j k₁₁ S(z) δ = ω/v₀ − π/d ist die Verstimmung. k₁₁ beschreibt die Änderung der Phasengeschwindigkeit im Elektrodenbereich, k₁₂ die akustische Fehlanpassung (Reflexion) an der Elektrodenkante.
Gl. 2.86 / 2.87

Sowohl die Geschwindigkeitsänderung als auch der Impedanzsprung setzen sich aus einem elektrischen Anteil (Kurzschluss der tangentialen Feldkomponente durch die Metallisierung) und einem mechanischen Anteil (Masse und Steifigkeit der Elektrode) zusammen:

ΔZ/Z₀ = Relektr. + Rmech. · Hm/(2d)   ·   Δv/v₀ = (Δv/v₀)elektr. + (Δv/v₀)mech. · Hm/(2d) Mit Relektr. = −0,42·KP²/2, Rmech. = −0,5, (Δv/v₀)elektr. = −0,7·KP²/2 und (Δv/v₀)mech. = −0,02.
Gl. 2.90 / 2.91 / 2.92
Metallisierung Hm = 60 nm Aluminium Elektrodenabstand d = 3,635 µm Wellenzahl k₀ = π/d = 0,864 µm⁻¹ k₁₁ = 0,00109 µm⁻¹ k₁₂ = 0,00123 µm⁻¹ KP² = 0,0016
Γ(0) = S(0)R(0) = − j k₁₂γR coth(γR M d) + j(δ + k₁₁)   ,   γR = k₁₂² − (δ + k₁₁)² Reflexionsfaktor eines Reflektors aus M Elektroden mit Abstand d. Randbedingungen: R(0) = 1 (Amplitude der vorwärtslaufenden Welle maximal), S(M·d) = 0 (am Ende keine Reflexionen mehr).
Gl. 2.99
Diagramm · Reflexionsfaktor eines Reflektors nach Betrag und Phase gerechnet · Gl. 2.99 Abb. 2.22
Deutlich erkennbar ist das Bandpassverhalten des Reflektors: Nur Oberflächenwellen in einem schmalen Wellenlängenbereich werden zurückgeworfen. Die Stoppbandbreite wird von k₁₂ bestimmt, die Mittenfrequenz durch k₁₁ gegenüber v₀/(2d) nach unten verschoben. Mehr Elektroden erhöhen den Reflexionsgrad, verändern die Bandbreite aber kaum.

Für den Interdigitalwandler kommt eine dritte Differentialgleichung für den eingeprägten Strom sowie der Anregungskoeffizient k₁₃ hinzu (Gl. 2.100–2.102). Der Wandler wird elektrisch als Parallelschaltung von akustischem Leitwert Ga(ω), Blindleitwert Ba(ω) und statischer Kapazität N·Cs modelliert; beide Leitwerte sind über die Hilbert-Transformation verknüpft. Daraus folgen geschlossene Ausdrücke für sämtliche Streuparameter S₁₁, S₁₂, S₁₃, S₃₃ (Gl. 2.112–2.117).

Das ist der entscheidende Punkt für alles Weitere: Durch Zusammenschalten der S-Parameter-Blöcke von Wandlern und Reflektoren lassen sich beliebige Anordnungen als elektrisches Modell im Frequenzbereich simulieren — und weil S-Parameter die Standardbeschreibung der Hochfrequenztechnik sind, direkt in Schaltungssimulatoren einbinden.

Kapitel 3

Experimentelles

3.1.1 · Hauptsensor

Quarz-Resonatoren, 433,9 MHz

Von Siemens auf gesägtem Wafer bezogen, in Leadless-Chip-Carrier gebondet. ST-X-Quarz, Rayleigh-Typ, Aluminiumelektroden mit ca. 60 nm Metallisierungshöhe. 35 Elektrodenpaare je Wandler, 220 Elektroden je Reflektor, Apertur 755 µm.

3.1.2 · Vergleich

ZnO-Verzögerungsleitungen

Am Institut für Werkstoffe der Elektrotechnik I der RWTH Aachen entwickelt. Gesputterte ZnO-Dünnschicht auf Silizium, Platinelektroden (Korrosionsschutz) mit Titan-Haftvermittler, Lift-off-Strukturierung. Heizstrukturen und Temperaturmess­widerstände auf demselben Substrat integriert. Untersucht: 154,9 · 178,2 · 227,2 MHz.

3.1.3 · Referenz

Quarz-Dickenscherschwinger

10-MHz- und 30-MHz-Grundtonquarze von Telefilter (Teltow) in unverschweißten HC45-Gehäusen, Goldelektroden gegen Alterung durch Korrosion. Bei 30 MHz beträgt die Scheibendicke nur noch 55 µm.

3.2 Gassensitive Schichten

Die Polymere wurden im Hinblick auf eine möglichst hohe Sensitivität für Benzol ausgewählt. Da Benzol stark krebserregend ist, konnte an der vorhandenen Gasmischeranlage aus Arbeitsschutzgründen nicht damit gemessen werden — als Analyt diente durchgehend Toluol, dessen Molekülstruktur der von Benzol ähnelt.

CH₂ CH₂ * * n
PPX — Poly(p-xylylen)
Parylene N · Tg 70 °C, teilkristallin
Mn > 100 000 g/mol · Pyrolyse · 90–180 nm
Si CH₃ CH₃ O * * n
PDMS — Polydimethylsiloxan
Stationäre Phase der GC · Tg ≈ −127 °C, amorph
Spray-Coating · 5–70 nm
NN NN Cu
CuPc — Kupferphtalocyanin
organischer p-Halbleiter · vier benzolähnliche Ringe
Sublimation bei ≈ 400 °C / 10⁻⁶ mbar · ~80 nm
Tabelle 3.1 / 3.2 Verwendete gassensitive Schichten
KürzelSubstanzVerfahrenDicke TgStruktur
PPXPoly(p-xylylen)Pyrolyse (IMC Marburg) 90–180 nm70 °Cteilkristallin
PPASPoly(phenylen-alt-styrol)Pyrolyse (IMC Marburg) 20–90 nm120 °Camorph, ≈ 300 000 g/mol
PPBPXPoly(phenylbutyl-p-xylylen)Drop-Coating (Sican) 20–70 nm−3 °Camorph, > 130 000 g/mol
PDMSPolydimethylsiloxanSpray-Coating (Bürkert) 5–70 nm−127 °Camorph
CuPcKupferphtalocyaninSublimation (IAP Gießen) 76–84 nmkristallin, halbleitend

Beschichtungsverfahren

A · aus Lösung

Dip- / Drop-Coating

Eintauchen bzw. Auftropfen. Sehr einfach, aber Schlieren und Tropfen — Schichtdicke kaum reproduzierbar.

B · aus Lösung

Spin- / Spray-Coating

Rotation bzw. Airbrush im Trägergasstrom. Deutlich homogener als A, für gute Reproduzierbarkeit ist aber ein automatisiertes System nötig.

C · aus Gasphase

Pyrolyse

Sublimation des Dimers bei 120 °C → Zerlegung bei 650 °C → Polymerisation bei Raumtemperatur direkt auf dem Substrat. Sehr homogen, keine Risse oder Poren.

D · aus Gasphase

Thermische Sublimation

CuPc-Pulver im Sublimationsofen bei ≈ 400 °C, Druck 10⁻⁶ mbar. Ebenfalls sehr homogen; Strukturierung über Masken möglich.

Zur Schichtdickenbestimmung. Da die Polymere sehr weich sind, ist eine Profilometermessung kaum möglich. Die Dicken wurden daher rechnerisch aus der Frequenzverschiebung gegenüber dem unbeschichteten Sensor bestimmt — über die Wohltjen-Gleichung (OFW) bzw. Sauerbrey (BAW). Dabei musste für die Dichte durchweg 1 g/cm³ angenommen werden (für CuPc 1,63 g/cm³). Bei den per Sublimation erzeugten CuPc-Schichten, deren Dicke bekannt war, stimmen gemessene und gerechnete Werte gut überein — bei den inhomogenen Drop- und Spray-Schichten sind die Werte nur Näherungen.

Tabelle 3.3 Schichtdickenabhängige Frequenzverschiebung [kHz/nm]
SchichtOFW Quarz 433 MHzOFW ZnO 227 MHz OFW ZnO 178 MHzOFW ZnO 154,9 MHz BAW 30 MHzBAW 10 MHz
PPX210,2
PPAS23,50,023
PDMS23,5
PPBPX23,542,51,80,023
CuPc23

3.3 Messelektronik

Oszillator mit elektrisch einstellbarer Phasenlage

Die Dimensionierung erfolgt nach dem Konzept der Schleifenverstärkung: Bei der Resonanzfrequenz muss |S₂₁| > 1 und die Phasendrehung ein Vielfaches von 2π sein. Ein Test-Port mit 50 Ω an beiden Enden der aufgetrennten Schleife erlaubt die Messung genau dieser Größe mit dem Netzwerkanalysator.

Kern der Schaltung ist ein zweistufiger elektronischer Phasenschieber mit Varaktordioden: über 100 el.° Phasenänderung bei gleichzeitig verbesserter 50-Ω-Anpassung; bei 80° Phasenhub ändert sich die Dämpfung nur um 0,6 dB. Eine einstufige Variante schaffte lediglich 20–30 el.°.

Als Verstärker dienten HP MSA0611 (rauscharm, 18 dB) bzw. für stark bedämpfte Sensoren INA-02184 (35 dB). Ein antiparalleles Schottkydioden-Paar begrenzt die Amplitude, damit der Verstärker im linearen Bereich bleibt — sonst würden Dämpfungsänderungen über die Verstärkerphase auf die Frequenz durchschlagen.

Dämpfungsmessung und Auswerteeinheit

Statt einer Verstärkungsregelung wird die Dämpfung über Richtkoppler vor und hinter dem Sensor erfasst: je ein logarithmischer Verstärker, dahinter ein rauscharmer Differenzverstärker mit Verstärkung 100. Ein zusätzliches Π-Dämpfungsglied gleicht die Sensordämpfung aus, damit beide Eingangspegel etwa gleich sind.

Die Platine ist vierlagig ausgeführt: Versorgungen auf den Innenlagen, HF oben, Niederfrequenz unten, getrennte Masseflächen für HF, Digital und Gleichspannung.

Die Auswerteeinheit misst nach dem Prinzip der Multiperiodendauermessung: Ein Time-to-Digital-Converter-ASIC (Acam TDC-GP1) löst Zeiten zwischen zwei Triggersignalen bis 110 ps auf. Ein Mikrocontroller 80C517A steuert Multiplexer, programmierbaren Teiler und TDC, rechnet in Frequenzen um und schreibt auf PCMCIA-SRAM-Karte oder über RS-232 zum PC. Acht Kanäle sequentiell.

Messaufbau. Gasmischeranlage nach der Sättigungsmethode: trockene synthetische Luft wird durch die flüssige Phase des Analyten bzw. durch destilliertes Wasser geleitet; Massendurchflussregler mischen rechnergesteuert definierte Toluolkonzentrationen und Luftfeuchten. Die Waschflaschen werden auf konstant 25 °C geregelt (Heizmatten, Heizdraht an den Zuleitungen, Regelung über Multimeter und programmierbares Netzteil am IEEE-Bus). In der temperaturgeregelten Aluminium-Messkammer sitzen gleichzeitig 4 OFW-Sensoren sowie 6 Dickenscherschwinger mit integriertem Temperatur- und Luftfeuchtesensor. Die Streuparameter wurden alle 20–25 s über den IEEE-Bus vom Netzwerkanalysator geholt.

Kapitel 4

Ergebnisse

4.1 Elektrische Charakterisierung

Der unbeschichtete Resonator zeigt seine Hauptresonanz bei 433,9 MHz mit 10 dB Dämpfung; die am schwächsten bedämpfte Nebenresonanz liegt bei 431,9 MHz mit 20 dB. Durch die Beschichtung verschieben sich beide zu niedrigeren Frequenzen, und die Dämpfung der Hauptresonanz nimmt zu — die der Nebenresonanz bleibt dagegen nahezu unverändert oder sinkt sogar.

Mit zunehmender Schichtdicke nimmt die Phasensteilheit und damit die Resonanzgüte Q ab (Ausgangswert Q = 9200). Die Sensoren wurden ausschließlich im niedriger bedämpften Hauptresonanzpunkt betrieben.

Messdaten · Einfluss der Beschichtung auf Dämpfung und Güte Tabelle 4.2 · gemessen
Der Befund ist eindeutig: Die per Pyrolyse (PPX, PPAS) und Sublimation (CuPc) erzeugten Schichten bedämpfen den Sensor um eine Größenordnung schwächer als die per Drop-Coating (PPBPX) und Spray-Coating (PDMS) aufgetragenen — bei vergleichbarer Frequenzverschiebung. Ursache ist nicht in erster Linie die Viskosität des Materials, sondern die Inhomogenität der Schicht: Schlieren und ungleichmäßige Dicken führen zu zusätzlichen Wellenreflexionen.
Tabelle 4.1 / 4.2 Zunahme der Dämpfung und Abnahme der Resonanzgüte (Ausgangsgüte des unbeschichteten Sensors: Q = 9200)
BeschichtungVerfahren Dämpfung Hauptres.Dämpfung Nebenres. ΔQ Hauptres.ΔQ Nebenres.
PPXPyrolyse1,8 dB/MHz · 0,038 dB/nm 0 dB/MHz1200 /MHz · 25,2 /nm210 /MHz
PPASPyrolyse1,2 dB/MHz · 0,028 dB/nm 0 dB/MHz1800 /MHz · 42,3 /nm−100 /MHz
CuPcSublimation0,8 dB/MHz · 0,018 dB/nm −0,8 dB/MHz−45 /MHz · −1 /nm−710 /MHz
PPBPXDrop-Coating12 dB/MHz · 0,28 dB/nm 1,6 dB/MHz7300 /MHz · 172 /nm180 /MHz
PDMSSpray-Coating20 dB/MHz · 0,47 dB/nm 14 dB/MHz12600 /MHz · 296 /nm1650 /MHz

Bemerkenswert: Bei CuPc ist die Nebenresonanz nach der Beschichtung sogar weniger stark bedämpft und weist eine höhere Güte auf als beim unbeschichteten Sensor (negative Werte in der Tabelle). Ein Betrieb auf der Nebenresonanz brächte trotzdem nichts, weil deren absolute Dämpfung höher bleibt. Ein Resonator lässt sich aber durch gezieltes „Verstimmen“ so auslegen, dass auch die Hauptresonanz weniger stark bedämpft wird — siehe Kapitel 5.

Verzögerungsleitungen

Tabelle 4.3 PPBPX-beschichtete ZnO-Verzögerungsleitungen im Punkt niedrigster Einfügungsdämpfung. Grunddämpfung unbeschichtet: ca. 32 dB.
VerzögerungsleitungDämpfungszunahme bezogen auf SchichtdickeGüteabnahme bezogen auf Schichtdicke
154,9 MHz9,4 dB/MHz0,017 dB/nm 1040 /MHz1,9 /nm
178,2 MHz10 dB/MHz0,025 dB/nm 260 /MHz0,64 /nm
227,2 MHz14,2 dB/MHz0,056 dB/nm 2100 /MHz8,4 /nm

Der Einfluss der Schichtdicke auf Dämpfung und besonders auf die Güte nimmt mit steigender Frequenz stark zu — hochfrequente Sensoren dürfen nur sehr dünn beschichtet werden. Verglichen mit den PPBPX-beschichteten Resonatoren (ΔQ = 7300 /MHz) fällt die Güteabnahme bei den Verzögerungsleitungen (1040 bzw. 2100 /MHz) aber deutlich geringer aus: Ihre Güte hängt von der geometrisch festgelegten Laufzeit ab.

4.2 Mikroskopie · 4.3 Temperatureinfluss

Auflichtmikroskopie mit Interferenzkontrast

Was man auf den Sensoren sieht

PPX (Pyrolyse): recht homogen; beim dünnsten Sensor einzelne Tropfen, die auf nur partielle Bedeckung schließen lassen.
PPAS (Pyrolyse): sehr homogen, keine Schlierenbildung, nur Tröpfchen am Rand.
PDMS (Spray): bei den dünnsten Sensoren keine durchgängige Schicht, nur vereinzelte Tröpfchen auf den Bondpads.
PPBPX (Drop): beim dünnsten Sensor nur Einzeltröpfchen; dickere Schichten geschlossen, aber mit starken Schlieren — bei einem Sensor führte das zur Anregung eines zusätzlichen Schwingungsmodes.

Temperaturgang, 20–45 °C

Der Referenzsensor hilft nicht

Im untersuchten Bereich hängen Resonanzfrequenz und Dämpfung aller Quarz-Resonatoren linear von der Temperatur ab. Entscheidend ist aber: Bei den mit PPBPX und PPX beschichteten Sensoren unterscheidet sich die Frequenzänderung mit der Temperatur sogar im Vorzeichen vom unbeschichteten Element.

Die Änderungen der Polymereigenschaften mit der Temperatur wirken sich also stärker aus als der Temperaturgang des Substrats. Eine Kompensation über ein unbeschichtetes Referenzelement ist damit nicht möglich. Bei den ZnO-Verzögerungsleitungen ist der Temperaturgang ohnehin um Größenordnungen höher: 8,5 kHz/°C (227 MHz, ≙ 37 ppm/K) bzw. 4 kHz/°C (178 MHz, ≙ 22 ppm/K) gegenüber 0,5 kHz/°C beim Quarz-Resonator.

4.4 Gassensitive Eigenschaften

Das Messprofil: Bei zwei verschiedenen Luftfeuchten je drei Toluolpulse unterschiedlicher Konzentration (100 / 200 / 300 ppm), anschließend ausgehend von ca. 33 % r. F. eine niedrige (20 %) und eine hohe Luftfeuchte (60 %). Die Kammertemperatur wurde über eine Regelung und eine Klimaanlage im Raum auf ±0,5 °C konstant gehalten.

Fast alle Sensoren reagieren auf Toluolpulse mit Frequenzabnahme und Dämpfungszunahme. Die höchste Toluolsensitivität — sowohl in der Frequenz als auch in der Dämpfung — zeigen die PPBPX-beschichteten Sensoren; starke Dämpfungszunahmen auch bei PDMS.

Die Ausnahme: Sensor PPX-c. Als einziger untersuchter Sensor reagiert er auf steigende Toluolkonzentration mit einer Zunahme der Frequenz und einer Abnahme der Dämpfung — und das ausgeprägt nichtlinear: Bei 300 ppm ist die Frequenzverschiebung sogar kleiner als bei den niedrigeren Konzentrationen. Eine Frequenzzunahme lässt sich mit einer Massenzunahme nicht erklären; verantwortlich ist eine Überlagerung zusätzlicher Effekte — Änderungen von Elastizität und Viskosität, die der massenbedingten Frequenzabnahme entgegenwirken (Quelleffekt, vgl. Kapitel 2.3.5).

Beim Luftfeuchteeinfluss liegt der unbeschichtete Quarz-Resonator bei ca. 25 Hz je 1 % r. F. Bei PPAS und PPX ist der Einfluss geringer als beim unbeschichteten Element, bei PDMS etwa gleich, bei PPBPX stark schichtdickenabhängig. Grund für die große Streuung: Die Quarzoberfläche ist stark hydrophil — unvollständige Benetzung, Risse oder Poren im Polymerfilm wirken sich unmittelbar auf die Feuchteempfindlichkeit aus.

Messdaten · Toluolsensitivität über der Schichtdicke Tabelle 5.1 · gemessen
Alle Sensoren wurden per Drop-Coating mit PPBPX beschichtet. Die absolute Sensitivität steigt mit der Schichtdicke — die schichtdickenbezogene Sensitivität ist bei Schichten unter 50 nm jedoch deutlich größer. Erklärung: Der gassensitive Effekt beruht zum einen auf Adsorption an der Oberfläche (unabhängig von der Dicke), zum anderen auf Absorption im Volumen. Für den Vergleich wurden die Verzögerungsleitungswerte über den Faktor aus Gl. 5.3 auf 434 MHz und Quarz umgerechnet: Bei gleicher Beschichtung, Dicke, Substrat und Frequenz ist die Massensensitivität von Resonator und Verzögerungsleitung gleich.

4.4.2 / 4.4.5 Der Einfluss der Phasenlage

Aus den kontinuierlich aufgezeichneten Streuparametern lässt sich mehr herausholen als nur die Resonanzfrequenz. Trägt man |S₂₁| und die Frequenz über der Phase auf, entspricht das genau dem Verhalten in einer Oszillatorschaltung: Der Sensor wird stets bei derselben Phasenlage betrieben; variabel sind Frequenz und Dämpfung.

Für jeden Sensor wurden 7 Arbeitspunkte ausgewertet — der Punkt minimaler Einfügungsdämpfung sowie je drei Punkte in 10-el.°-Schritten darüber und darunter. Bestätigt wurde das Ergebnis durch Oszillatormessungen, bei denen die Steuerspannung des Phasenschiebers über 12 Stunden von 0,75 V bis 11 V durchgefahren wurde (≙ ca. 80 el.°).

Messdaten · Sensitivität in Abhängigkeit der Phasenlage Abschnitt 4.4.5 · gemessen
Dargestellt sind die im Text genannten Eckwerte bei −5 el.° und −65 el.°, linear interpoliert und auf das jeweilige Maximum normiert. Bei PDMS-b fällt die Toluolsensitivität über den Phasenhub auf 67 %, die Luftfeuchtesensitivität dagegen auf 38 % — die Kurven laufen deutlich auseinander. In der Oszillatormessung entsprach das einer Variation der Toluol-Sensitivität um 60 % (1500 → 600 Hz bei 300 ppm) gegenüber 33 % bei der Luftfeuchte. Genau diese Differenz erlaubt die Unterscheidung mit einem einzigen Sensor.

Zweiter Auswerteweg: das Verhältnis Frequenz zu Dämpfung. Bei den meisten Sensoren hat die Luftfeuchte einen deutlich größeren Einfluss auf die Dämpfung als auf die Frequenz — bei Toluol ist es genau umgekehrt. Sensor PPBPX-d etwa: 0,2 mdB/ppm Toluol (≙ 60 mdB bei 300 ppm) gegenüber 0,6 mdB/% Luftfeuchte (≙ 18 mdB bei 30 %). Die gleichzeitige Messung beider Größen ermöglicht also ebenfalls eine Unterscheidung. Die praktische Hürde: Dämpfungsänderungen von wenigen mdB liefern mit der hier entwickelten Auswerteschaltung stark verrauschte Signale, während Frequenzen mit einer Auflösung von unter 0,1 Hz gemessen werden können.

Kapitel 5

Diskussion — das erweiterte COM-Modell

Die Modelle aus Kapitel 2.3 liefern nur die Änderung der Ausbreitungskonstante ΔγR/k₀ = Δα/k₀ − j·Δv/v₀. Sie betrachten ein unendlich ausgedehntes Substrat mit homogener Beschichtung — Anregung, Elektrodeneinfluss und Randeffekte bleiben außen vor.

Für die Frequenzverschiebung reicht das: Δf ist in erster Näherung proportional zu Δv. Für die Einfügungsdämpfung jedoch nicht: Δα ist ein längenbezogener Dämpfungsbelag; man müsste ihn mit der Laufstrecke multiplizieren. Bei einer Verzögerungsleitung ist das der Wandlerabstand — bei einem Resonator mit seinen Mehrfachreflexionen lässt sich aus der Geometrie überhaupt keine Laufstrecke ableiten.

Die Idee. Im COM-Ansatz wird die Ausbreitung ohnehin über die Ausbreitungskonstante γR und die Verstimmung δ beschrieben. Man kann also die aus dem viskoelastischen Modell berechnete relative Änderung einfach dort einsetzen — und erhält das komplette Übertragungsverhalten des beschichteten Sensors.

γSchicht = j ( ω/v₀ − k₀ · ( Δv/v₀ + j·Δα/k₀ ) )    ·    δSchicht = ω/v₀ − π/d + k₀ · ( Δv/v₀ + j·Δα/k₀ ) Die erweiterten Parameter werden in die COM-Gleichungen für Reflexionsfaktor (Gl. 2.99) und Wandler-Streuparameter (Gl. 2.112–2.115) eingesetzt. Die unmetallisierten Bereiche zwischen den Wandlern bzw. zwischen Wandlern und Reflektoren werden als Verzögerungsleitung mit der S-Matrix [[0, e−γRD],[e−γRD, 0]] modelliert.
Gl. 5.1 / 5.2 — der Kern dieser Arbeit
Validierung

Modell vs. Messung

Der Vergleich des simulierten mit dem gemessenen |S₂₁|-Verlauf für den unbeschichteten Resonator zeigt: Das elektrische Verhalten inklusive aller Nebenresonanzen im Durchlassbereich wird gut wiedergegeben. Deutliche Abweichungen treten erst am Rand des Durchlassbereichs auf — für die Oszillatorsimulation unerheblich.

Erklärt: Nichtlinearität

Warum die Dämpfung nicht linear steigt

Ein Polymer geringer Viskosität führt bei homogener Beschichtung zu keiner erkennbaren Dämpfungszunahme. Mit steigender Viskosität nimmt sie zu — aber nichtlinear: schon bei dünnen Schichten stark, bei dickeren nur noch wenig. Grund: Mit zunehmender Dämpfung sinkt die Zahl der Reflexionen und damit die virtuelle Laufzeit der Welle. Genau dieses Verhalten wurde bei PDMS und PPBPX gemessen.

Erklärt: Nebenresonanzen

Warum sie schwächer bedämpft werden

Die Reflektoren bestehen aus M = 220 Elektroden, die Wandler nur aus N = 35 Paaren — die Frequenzcharakteristik der Reflektoren ist also viel schmalbandiger. Das Bandpassverhalten des Resonators entsteht aus dem Zusammenspiel von Reflektorcharakteristik und der Bedingung, dass der Reflektorabstand ein ganzzahliges Vielfaches von λ/2 sein muss. Simulation und Messung (PPX, PPBPX) stimmen überein.

5.1.5 Resonator oder Verzögerungsleitung?

Bei gleicher Beschichtung, Dicke, Substrat und Resonanzfrequenz ist die Massensensitivität beider Bauformen gleich. Der Unterschied liegt im Rauschen — und damit in der Nachweisgrenze. Als Vergleichsgröße dient die rauschäquivalente Gaskonzentration: die Konzentration, die zu einer Frequenzverschiebung entsprechend der dreifachen Standardabweichung des Frequenzrauschens führt.

SφO(f) = FV·GV·k·T2 · ( 1 + fcVf − fT ) · ( 1 + ( fT2(f − fT)·Q )² ) Leistungsspektrum des Phasenrauschens bei geschlossener Schleife. Das breitbandige Verstärkerrauschen wird über den schmalbandigen Resonator gefiltert — die Güte Q geht direkt ein. Multiplikation mit (2π(f−fT))² liefert das Frequenzrauschen, Integration über die Allan-Varianz die Standardabweichung. Verwendete Werte: Rauschzahl FV = 2 dB, Eckfrequenz fcV = 1 kHz (INA-02184), Torzeit τ = 1 s; GV muss die Einfügungsdämpfung kompensieren.
Gl. 5.6 – 5.8

Was die Simulation zeigt

  • Die Dämpfung der Verzögerungsleitung steigt linear mit der Schichtdicke — beim Resonator nichtlinear.
  • Die Güte des Resonators nimmt deutlich stärker ab, weil sie von der stark beschichtungsabhängigen virtuellen Laufzeit abhängt.
  • Bei beiden Typen sinkt die Nachweisgrenze zunächst stark mit der Schichtdicke, dann kaum noch. Beim Resonator steigt sie ab G″ > 0,4·10⁸ Pa jenseits eines Optimums wieder an.
  • Niedrigviskose Schichten → Resonator gewinnt. Hochviskose, dicke Schichten → Verzögerungsleitung gewinnt.

Optimierung der Elektrodengeometrie

  • Verzögerungsleitung: Die Länge hat kaum Einfluss auf die Nachweisgrenze — nur bei sehr dicken Schichten sind kurze Leitungen besser. Entscheidend ist die Elektrodenzahl: 70 statt 35 Elektrodenpaare halbieren die Einfügungsdämpfung und senken die Nachweisgrenze deutlich.
  • Aber: Mehr Elektroden machen die Eingangsimpedanz niederohmiger und stärker kapazitiv — das muss die Auswerteschaltung auffangen. Zu breite Elektroden regen zudem störende Transversalwellen an (im 1D-COM-Modell nicht erfasst).
  • Elektrodenmaterial: Aluminium hat die geringere akustische Fehlanpassung, Gold und Platin die bessere Langzeitstabilität — für Sensoren im offenen Gehäuse entscheidend.

Der stärkste Optimierungsansatz: den Resonator absichtlich verstimmen. Beim Resonator ist die akustische Energie im Bereich zwischen den beiden Wandlern konzentriert. Beschichtet man nur diesen Resonanzraum und verkürzt gleichzeitig den Wandlerabstand um λ/4 (D = 25λ − λ/4), so kompensiert die durch die Beschichtung verursachte niedrigere Ausbreitungsgeschwindigkeit genau diese Verstimmung: Das optimale Dämpfungsverhalten stellt sich erst im beschichteten Zustand ein. Ergebnis der Simulation: eine bis zu zweifach niedrigere Nachweisgrenze — und das trotz geringerer Massensensitivität (Frequenzverschiebung 9 kHz/nm statt 24 kHz/nm bei vollständiger Beschichtung). Der Ansatz funktioniert allerdings nur für relativ dünne Schichten; bei dickeren wandert die Hauptresonanz zu weit in Richtung der Nebenresonanz.

5.2 Betrieb und Auswertung

5.2.1 Warum die Phasenlage die Sensitivität verändert

Zwei Effekte überlagern sich im Phasenverlauf von S₂₁: Die Abnahme der Phasengeschwindigkeit verschiebt den Verlauf zu kleineren Frequenzen; die zunehmende Dämpfung verringert die Phasensteilheit, der Verlauf fällt flacher ab. Betrachtet man den Sensor bei konstanter Phasenlage — genau das tut ein Oszillator —, hängt die resultierende Frequenzänderung von beiden Effekten ab.

Die Simulation reproduziert das gemessene Verhalten: Bei stark bedämpften Sensoren nimmt die Frequenzsensitivität mit zunehmend positiver Phasenlage zu; bei ca. −40 el.° bleibt die Frequenz konstant, darunter wechselt das Vorzeichen.

Da Luftfeuchte bei vielen Sensoren viel stärker bedämpft als Toluol, ist die Luftfeuchte-Sensitivität deutlich phasenabhängiger — die Grundlage für die Unterscheidung mit einem einzelnen Sensor.

Eine mögliche Erklärung für ein bekanntes Problem: In der Literatur wird berichtet, dass mit bestimmten OFW-Sensoren und Beschichtungen erreichte Sensitivitäten von anderen Arbeitsgruppen oft nicht reproduziert werden können. Neben unterschiedlicher Schichtbeschaffenheit kann auch die nicht angegebene Phasenlage der verwendeten Oszillatorschaltung dafür verantwortlich sein.

5.2.2 Einfluss des Anpassungsnetzwerks

Für die Nachweisgrenze zählt nicht die Güte des Sensors allein, sondern die der Reihenschaltung aus Anpassnetzwerken und Sensor. Untersucht wurden sechs Parametrisierungen, die alle eine Anpassung an 50 Ω bewirken — mit deutlich unterschiedlichem Ergebnis:

Tabelle 5.2 Untersuchte Anpassnetzwerke (Sensor PPX-b)
CShunt [pF]LSerie [nH] RSerie [Ω]Bewertung
1,0480breitbandig, nicht optimal
2,237,516,5 niedrigste Nachweisgrenze
3,13126
3,92632
5,12137,5
6,31841schmalbandig, hohe Dämpfung

Je hochohmiger der Serienwiderstand, desto größer die Einfügungsdämpfung — aber desto schmalbandiger auch der Resonanzpeak. Das Optimum liegt dazwischen: weder die 0-Ω-Anpassung noch die hochohmigste. Die Frequenz-Sensitivität wird durch das Netzwerk kaum beeinflusst, die Dämpfungs-Sensitivität dagegen erheblich.

Kapitel 6

Wertung und Ausblick

Untersucht wurde die Eignung polymerbeschichteter Oberflächenwellen-Bauelemente als Gassensoren — mit Temperatur und Luftfeuchte als wichtigsten Störgrößen und Toluol als Analytgas.

Was sich gezeigt hat

  • Störeinflüsse hängen stark vom verwendeten Polymer ab — und noch stärker von der Beschaffenheit der Beschichtung. Drop-Coating erwies sich als besonders schlecht reproduzierbar. Pyrolyse und Sublimation sind für hochfrequente Sensoren, die möglichst dünn beschichtet werden müssen, wesentlich besser geeignet.
  • Temperaturkompensation über ein unbeschichtetes Referenzelement funktioniert nicht. Ein Betrieb bei konstanter Temperatur ist unerlässlich; die erreichten ±0,5 °C sind als Mindestanforderung anzusehen. Der Ansatz, Temperatursensoren und Heizelemente auf dem Sensorsubstrat zu integrieren (ZnO auf Silizium), wurde hier noch nicht genutzt.
  • Die Wahl des Arbeitspunkts durch die Auswerteschaltung spielt eine große Rolle: Bei einigen Sensoren ist die Luftfeuchte-Sensitivität bei niedrigen Phasenlagen um mehr als 50 % größer als bei hohen.
  • Für zukünftige Arbeiten müssen Rauschabstand und Reproduzierbarkeit der Dämpfungsmessung verbessert werden — Änderungen von wenigen mdB sind mit der hier entwickelten Schaltung zu verrauscht.

Das entwickelte Modell — Stand und Grenzen

Das erweiterte Coupling-of-Modes-Modell berechnet das Übertragungsverhalten eines OFW-Sensors in Abhängigkeit von Dicke, Dichte und viskoelastischen Eigenschaften der Beschichtung — für beliebige Frequenzbereiche, unter Berücksichtigung von Substrateigenschaften und Elektrodengeometrie. Die Übereinstimmung mit den Netzwerkanalysator-Messungen ist für alle Sensoren gut.

Noch nicht berücksichtigt:

  • die elektrischen Eigenschaften der Beschichtung
  • der Einfluss der Temperatur
  • chemische Modelle wie LSER — dafür müssten erst die Auswirkungen absorbierter Gasmoleküle auf die physikalischen Schichteigenschaften bekannt sein
  • zeitliche Verläufe: Massenzunahme, Änderung der viskoelastischen und der elektrischen Eigenschaften haben vermutlich unterschiedliche Zeitkonstanten. Kennt man sie, lässt sich aus dem Signalverlauf deutlich mehr Information gewinnen — etwa über neuronale Netze.

Zur experimentellen Bestätigung der simulierten Optimierungsansätze braucht es Bauelemente mit verschiedenen Elektrodenkonfigurationen auf demselben Substrat sowie sehr gut reproduzierbare Beschichtungen. Unter diesen Voraussetzungen ist es auch sinnvoll, zu höheren Resonanzfrequenzen zu gehen. Offen ist ebenfalls, wie sich eine gezielte Teilbeschichtung — etwa über Masken — praktisch umsetzen lässt.

±0,5 °C Mindestanforderung an die Temperaturregelung
> 50 % Variation der Sensitivität über die Phasenlage
niedrigere Nachweisgrenze durch gezielt verstimmten Resonator
0,1 Hz Frequenzauflösung der Multiperiodendauer-Messung
Redaktionelle Anmerkungen

Nachgerechnet — und was dabei auffiel

Beim Aufbau dieser Webseite wurden die zentralen Rechnungen der Arbeit mit ihren eigenen Konstanten unabhängig nachvollzogen. Sie tragen. Dabei sind einige Stellen aufgefallen, die redaktionell nicht sauber sind — sie sind hier der Vollständigkeit halber dokumentiert.

Methodische Einschränkung. Grundlage ist der aus dem PDF extrahierte Text, nicht das gesetzte Original. Formeln und Abbildungen extrahieren unsauber; die folgenden Befunde beschränken sich deshalb auf Fließtext sowie auf Zahlenwerte, die sich unabhängig nachrechnen ließen. Ein Abgleich mit dem Originaldruck ist bei jedem Punkt angeraten.

Die Gegenprobe

Sechs unabhängige Nachrechnungen aus den Gleichungen der Arbeit gegen die dort angegebenen Werte. Alle Diagramme dieser Seite beruhen auf denselben Rechnungen.

Verifikation Nachrechnung mit den Konstanten aus Tabelle 2.2 und den Bauelementdaten aus Abschnitt 2.4.2
PrüfungAngabe in der Arbeit nachgerechnetBefund
Frequenzverschiebung pro nm Schichtdicke (Gl. 2.31) 21 – 23,5 kHz/nm Tab. 3.3 23,7 kHz/nm✓ stimmt
Filmresonanz bei G′ = 10⁷ Pa (Gl. 2.51) „ca. 60 nm“ Abb. 2.15 58,5 nm✓ stimmt
Filmresonanz bei G′ = 10⁸ Pa (Gl. 2.51) „ca. 190 nm“ Abb. 2.16 187 nm✓ stimmt
Grenze RA = 1, weiches Polymer (Gl. 2.40) 3,8 nm Abschn. 2.3.2 3,9 nm✓ stimmt
Mittenfrequenz des Reflektors (Gl. 2.99) gemessen 433,9 MHz Abschn. 4.1.1 433,84 MHz✓ stimmt
Umrechnungsfaktoren Resonator ↔ Verzögerungsleitung (Gl. 5.3) 10 bzw. 6 Gl. 5.4 9,95 / 6,12✓ stimmt

Befunde

inhaltlich

1 Stelle

Ein Wortdreher, der die Argumentation an dieser Stelle umkehrt.

formal / widersprüchlich

6 Stellen

Formelschreibweise, Symbolverzeichnis, Zahlen- und Satzwidersprüche.

Tippfehler

9 Stellen

Reine Korrektoratsfehler ohne inhaltliche Wirkung.

Liste Auffälligkeiten, nach Gewicht geordnet
EinordnungStelleBefund
inhaltlich S. 111, Abschn. 4.4.1
Sensor PPX-c
„Eine Abnahme der Frequenz kann mit einer Massenzunahme nicht erklärt werden.“
Muss Zunahme heißen. Unmittelbar davor steht, dass dieser Sensor als einziger mit steigender Toluolkonzentration zu höheren Frequenzen wandert. Eine Abnahme ist ja gerade das, was eine Massenzunahme erklärt — in der gedruckten Fassung kehrt der Satz die Begründung um.
Formel Gl. 3.4, Abschn. 3.4.1 S11 Osz(ω) = S21 Osz(ω)
Der Satz davor fordert, dass die Eingangsimpedanz beider Enden der offenen Schleife gleich ist — das ist S₁₁ = S₂₂. Abbildung 3.24 trägt konsequenterweise auch „Parameter S11, S22“ im Titel. Ein Reflexions- mit einem Transmissionsparameter gleichzusetzen ergibt physikalisch keinen Sinn.
Verzeichnis S. 4, Symbolverzeichnis R(z) — Amplitude der vorwärtslaufenden Welle
R(z) — Amplitude der rückwärtslaufenden Welle

Der zweite Eintrag muss S(z) lauten; so wird das Symbol im gesamten Kapitel 2.4.2 auch verwendet.
Zahlen Abschn. 4.3.2
ZnO-Verzögerungsleitungen
„Die stärkste Temperaturabhängigkeit wurde hier mit 2,5 kHz/K (entspr. 5,8 ppm/K) … gemessen.“
Zwei Sätze zuvor stehen 8,5 kHz/°C — 2,5 kHz/K ist also nicht die stärkste. Und 2,5 kHz entsprechen bei 227 MHz 11 ppm/K, bei 178 MHz 14 ppm/K, keinesfalls 5,8. Der Wert 5,8 ppm/K passt zu 2,5 kHz bei ≈ 431 MHz — also zum Quarz-Resonator. Hier sind offenbar Werte vertauscht.
Skalierung Abb. 2.15 / 2.16
Ordinate Δv/v₀ [%]
Mit den eigenen Konstanten (k₁+k₂ = −12,6·10⁻⁸, f₀ = 434 MHz, ρf = 1 g/cm³) ergeben sich bei 120 nm −0,66 %; die Abbildung reicht bis −6 %Faktor ≈ 10. Dass die Filmresonanz bei 58 nm liegt, belegt zugleich, dass genau diese Parameter verwendet wurden. Die Skalierung ist damit die wahrscheinliche Fehlerquelle, nicht die Rechnung: Tabelle 3.3 (21–23,5 kHz/nm) stimmt mit Gl. 2.31 überein. Die Aussagen der Abbildungen — Lage der Filmresonanz und Nichtlinearität — bleiben unberührt.
Zahlen Abschn. 2.2.1 gegen 2.4.2 Elektrodenabstand einmal 2d = 7,29 µm (⇒ d = 3,645 µm), einmal d = 3,635 µm. Entspricht 1,4 MHz Unterschied in der Resonanzfrequenz — marginal, aber inkonsistent.
Satzbau S. 92, Abschn. 4.1.1
Sensor PPBPX-e
Erst „deutlich geringere Dämpfungszunahme“, dann Gleiches gilt für die Güte …, die durch die Beschichtung wesentlich stärker reduziert wird“. „Gleiches“ und „stärker“ widersprechen sich; gemeint ist offenbar „Umgekehrtes gilt“.
Tippfehler verteilt in eiem Sputterprozeß · unverscheißten HC45-Gehäusen · Seonsoren · Dickenschwerschwinger · Imaginaerteil · propotional · Strühstrahl · Sensitiväten · an dieser Stelle diese Stelle darauf

Einordnung. Keiner dieser Punkte berührt die Substanz der Arbeit. Das erweiterte Coupling-of-Modes-Modell — der eigentliche wissenschaftliche Beitrag — ist sauber hergeleitet, die Messreihen sind in sich konsistent, und die Schlussfolgerungen halten der unabhängigen Nachrechnung stand. Was hier steht, sind Korrektoratsfehler, wie sie in praktisch jeder Dissertation zu finden sind — erst recht in einer von 2002, entstanden ohne die heute selbstverständlichen Prüfwerkzeuge. Der einzige Punkt mit inhaltlichem Gewicht ist der Wortdreher auf Seite 111; auf dieser Webseite ist er in der Darstellung des Sensors PPX-c stillschweigend korrigiert.

Anhang

Symbole, Kennzahlen und Quellen

Auswahl Verzeichnis verwendeter Symbole
ZeichenBedeutungEinheit
αDämpfungsbelag einer Rayleighwelle1/m
δVerstimmung im COM-Modell1/m
γRAusbreitungskonstante der Rayleighwelle1/m
βRzAusbreitungskonstante der z-polarisierten Teilwelle1/s
ΓReflexionsfaktor1
λf, µfLamé-Konstanten des viskoelastischen FilmsN/m²
ρf, hfDichte / Dicke der Beschichtungkg/m³ · m
σOberfl.Oberflächenleitfähigkeit1/(Ω·m)
σ₀²Allan-Varianz1
COberfl.Oberflächenkapazität (Kapazitätsbelag)F/m
Csstatische Kapazität eines ElektrodenpaarsF
dElektrodenabstand des Interdigitalwandlersm
DVZLWandlerabstand einer Verzögerungsleitungm
Ezzlongitudinales DehnungsmodulN/m²
FV, GVRauschzahl / Verstärkung des Verstärkers1
G = G′+jG″komplexes SchermodulN/m²
K = K′+jK″komplexes KompressionsmodulN/m²
HmMetallisierungshöhem
k₁, k₂Materialkonstanten der Wohltjen-Gleichungm²·s/kg
k₁₁, k₁₂, k₁₃Koppel- und Anregungskoeffizienten im COM-Modell1/m
KP²piezoelektrischer Kopplungsfaktor1
M, NZahl der Reflektorelektroden / Elektrodenpaare1
QResonanzgüte, Q = −f · δφ/δf1
RAKennzahl akustisch dick / dünn1
SijElement der Streuparametermatrix1
TGGlastemperatur eines PolymersK
TKfTemperaturkoeffizient der Frequenzppm/K
v₀, vEPhasengeschwindigkeit frei / kurzgeschlossenm/s
WApertur (Überlappungslänge) der Elektrodenm

Kenndaten der verwendeten Bauelemente

Resonanzfrequenz 433,9 MHz Substrat ST-X Quarz λR = 7,3 µm 2d = 7,29 µm N = 35 Elektrodenpaare M = 220 Reflektorelektroden Apertur W = 755 µm Wandlerabstand 180 µm ≙ 25 λ Hm = 60 nm Aluminium Qunbeschichtet = 9200 Einfügungsdämpfung 10 dB

Zentrale Referenzen der Arbeit

  • [Aul90] Auld: Acoustic Fields and Waves in Solids — Störungsanalyse, Grundlage der Gl. 2.26
  • [Woh84] Wohltjen — Massensensitivität dünner elastischer Filme (Gl. 2.31)
  • [Sau59] Sauerbrey — Dickenscherschwinger als Massensensor (Gl. 2.33)
  • [Fry94] Frye — Modell für viskoelastische, akustisch dicke Filme (Gl. 2.46–2.51)
  • [Ric85] Ricco et al. — totales Differential (Gl. 2.25), leitfähige Schichten
  • [Tho88] · [Che85] · [Hau77] — Coupling-of-Modes-Theorie und Koppelkoeffizienten
  • [Bel73] Bell, Hartmann, Rosenfeld — Impulsmodell des Interdigitalwandlers
  • [Mau95] · [Mau96] — Vergleich Resonator/Verzögerungsleitung, Rauschmodelle nach Leeson und Kurokawa
  • [Bar71] Barnes et al. — Allan-Varianz (Gl. 5.8)
  • [Ben98] Bender, Mokwa — temperaturstabilisierte ZnO-Sensoren auf Silizium
  • [Abr94] · [Gra95] — LSER-Theorie und Polymerauswahl
  • [Rap98] — Phasenlage bei Ein-Tor-Resonatoren

Danksagung (aus der Arbeit). Prof. Dr. W. Mokwa · Prof. Dr. D. Kohl (Koreferat) · Prof. Dr. P. Kaul (Anregung und Betreuung) · den Partnern des BMBF-Verbundprojekts NoseBox, insbesondere S. Bender (IWE I, RWTH Aachen) für die Zinkoxid-Verzögerungsleitungen · M. Ishaque und F. Brink-Spalink (Institut für Makromolekulare Chemie, Universität Marburg) für die Polymerbeschichtungen · V. Uwira (Institut für Angewandte Physik, Universität Gießen) für die Kupferphtalocyanin-Schichten · Prof. F. L. Dickert (Universität Wien) · M. Lösel (Bürkert) für die Spray-Coating-Beschichtungen · Dr. G. Fischerauer (Siemens) für die Quarz-Resonatoren · den Mitarbeitern der SICAN GmbH und des NIGAN e.V., besonders Dr. H. D. Hartmann, sowie den Praktikanten und Diplomanden M. Descher, E. Boerner und H. Homeyer.