Tutaj prezentujemy protokół do wykonywania czułej, przestrzennie rozdzielczej spektroskopii gazów w obszarze średniej podczerwieni, wykorzystując zdegenerowane mieszanie czterofalowe w połączeniu z detekcją konwersji w górę.
Artykuł metodologiczny
Tutaj prezentujemy protokół do wykonywania czułej, przestrzennie rozdzielczej spektroskopii gazów w obszarze średniej podczerwieni, wykorzystując zdegenerowane mieszanie czterofalowe w połączeniu z detekcją konwersji w górę.
Przedstawiamy protokół do wykonywania spektroskopii gazowej za pomocą czterofalowego mieszania w podczerwieni (IR-DFWM), do ilościowego wykrywania gatunków gazów w zakresie ppm do pojedynczego procentu. Głównym celem metody jest przestrzennie rozdzielcza detekcja gatunków o niskim stężeniu, które nie mają przejść w zakresie widma widzialnym lub bliskiej podczerwieni, które można by wykorzystać do detekcji. IR-DFWM jest metodą nieinwazyjną, co jest wielką zaletą w badaniach spalania, ponieważ włożenie sondy do płomienia może go drastycznie zmienić. IR-DFWM jest połączony z detekcją konwersji w górę. Ten schemat detekcji wykorzystuje generowanie częstotliwości sumarycznej do przenoszenia sygnału IR-DFWM ze środkowej podczerwieni do obszaru bliskiej podczerwieni, aby wykorzystać doskonałą charakterystykę szumów detektorów na bazie krzemu. Proces ten odrzuca również większość termicznego promieniowania tła. Przedstawiony tutaj protokół koncentruje się na prawidłowym ustawieniu optyki IR-DFWM oraz na tym, jak ustawić system detekcji konwersji w górę wewnątrz wnęki.
IR-DFWM zapewnia możliwość pomiaru stężeń aktywnych substancji IR aż do poziomu ppm1, z rozdzielczością przestrzenną. IR-DFWM ma kilka zalet, które czynią go atrakcyjną techniką do badań spalania. Płomienie można drastycznie zmienić poprzez włożenie sond, ale IR-DFWM jest nieinwazyjny. Ma rozdzielczość przestrzenną, dzięki czemu można mierzyć stężenia gatunków w różnych punktach struktury płomienia. Zapewnia spójny sygnał, który można odizolować od emisji termicznej płomienia. Dodatkowo, DFWM jest mniej wrażliwy na środowisko zderzenia niż na przykład fluorescencja indukowana laserem (LIF), która może być trudna do określenia w płomieniu. Technika ta zapewnia również dostęp do gatunków molekularnych, które są aktywne w podczerwieni, ale nie mają widocznych lub prawie widocznych przejść, które można wykorzystać do ich pomiaru za pomocą innych technik.
Chociaż DFWM ma wiele zalet, alternatywne techniki mogą być preferowane, jeśli jedna lub więcej z tych zalet nie jest potrzebna. Jeśli rozdzielczość przestrzenna nie jest potrzebna, techniki oparte na absorpcji będą zarówno prostsze, jak i dokładniejsze. Jeśli dane gatunki molekularne mają przejścia w obszarze widzialnym lub bliskiej podczerwieni, LIF może być preferowany, ponieważ LIF może dostarczyć przestrzennie rozdzielone informacje z płaszczyzny, a nie tylko z jednego punktu. W odpowiednich warunkach metody nieliniowe, takie jak DFWM i PS, mogą być również używane do pojedynczych pomiarów 2D2. Sygnał tych metod nieliniowych jest proporcjonalny do natężenia wiązki sondy w kostce, a ponieważ wiązka pompy musi zostać rozszerzona, aby pokryć obszar pomiaru 2D, wymaga to do działania albo bardzo wysokich energii impulsów, albo kombinacji wysokiej podatności trzeciego rzędu, wysokich stężeń i niskiego szumu tła. Dlatego w większości zależy to od gatunku molekularnego, czy jest to możliwe.
W bardziej bezpośredniej konkurencji z DFWM, istnieją inne techniki spektroskopowe oparte na mieszaniu czterech fal: koherentna spektroskopia Ramana anty-Stokesa (CARS), spektroskopia siatki indukowanej laserem (LIGS) i spektroskopia polaryzacyjna (PS). CARS to dobrze znana technika pomiaru temperatury i głównych gatunków w środowisku spalania. Brakuje mu jednak czułości do wykrywania mniejszych gatunków, ponieważ granica wykrywalności wynosi zwykle około 1%2. Wcześniej wykazano, że PS i DFWM mają podobne granice czułości i wykrywalności3; jednak wykazano, że stosunek sygnału do szumu DFWM wzrasta o współczynnik 500 w połączeniu z wykrywaniem konwersji w górę4, podczas gdy PS wykazał tylko 64-krotny wzrost5. LIGS ma tę zaletę, że indukuje siatkę za pomocą światła średniej podczerwieni, ale mierzy efekt przez załamanie lasera sondy tej siatki, a długość fali tego lasera sondy można dowolnie wybierać6. Długość fali lasera sondy może zatem znajdować się w obszarze widzialnym, gdzie dostępne są szybkie, niskoszumowe detektory na bazie krzemu. Jest to ta sama zaleta, którą osiąga się dzięki zastosowaniu konwersji w górę. LIGS ma tę wadę, że jest bardzo wrażliwy na kolizje2, co oznacza, że stężenie głównych gatunków gazów musi być znane do precyzyjnych pomiarów stężenia lub temperatury za pomocą LIGS. Jeśli ten problem zostanie rozwiązany, LIGS ma podobną czułość zarówno dla DFWM, jak i PS przy ciśnieniu atmosferycznym3, ale tam, gdzie sygnał LIGS wzrasta wraz ze wzrostem ciśnienia, sygnał z DFWM i PS wzrasta przy niższych ciśnieniach, co oznacza, że preferowana technika będzie zależeć od środowiska ciśnienia.
Wykrywanie konwersji w górę to technika konwersji sygnału z długich fal na krótsze za pomocą generowania sumy częstotliwości. Zaletą tego rozwiązania jest to, że detektory w zakresie widzialnym lub bliskiej podczerwieni mają niższy szum i wyższą czułość niż ich odpowiedniki w zakresie średniej podczerwieni. Po raz pierwszy zbadano to pięć dekad temu7, ale od tego czasu spotkało się to z bardzo małym zainteresowaniem i zastosowaniem, ze względu na niską wydajność konwersji. Jednak dzięki postępom w technikach produkcji okresowo spolaryzowanego niobianu litu (PPLN) i innych materiałów o wysokich współczynnikach nieliniowych, a także zwiększonej dostępności diod laserowych dużej mocy (LD), technika ta przyciągnęła coraz większą uwagę w ostatniej dekadzie, a jej zastosowania obejmują takie obszary, jak detekcja pojedynczych fotonów w średniej podczerwieni8,9,10,11, IR lidar12,13, oraz obrazowanie hiperspektralne14,15 i microscopy16. Główną zaletą połączenia detekcji konwersji w górę z IR-DFWM jest to, że warunek dopasowania fazowego ma wąskie pasmo akceptacji kątowej i spektralnej, które silnie dyskryminuje tło termiczne, umożliwiając wykrywanie słabszych sygnałów.
Konfiguracja detektora konwersji w górę jest pokazana w Rysunek 1; lustra, soczewki lub inne optyki, do których odwołuje się protokół, są zidentyfikowane tutaj lub na schemacie konfiguracji IR-DFWM pokazanym w Rysunek 2. Sekcja protokołu zajmuje się głównie dostosowaniem konfiguracji optycznej używanej w tej metodzie, a proces można wstrzymać w dowolnym momencie, wyłączając wszystkie uruchomione urządzenia. Wszystkie lusterka są regulowane ręcznie. Oprogramowanie służące do sterowania kamerą i LD zostało dostarczone wraz z detektorem konwersji w górę. Sposób korzystania z oprogramowania jest opisany na końcu protokołu.
1. Konwersja w górę
2. Wyrównanie IR-DFWM
UWAGA: Zobacz Rysunek 2 dla diagramu konfiguracji DFWM.
3. Oprogramowanie diod laserowych
4. Systemy rozwoju obrazowania IDS
Rysunek 3 przedstawia sygnał dla różnych stężeń HCN w N2, uśredniony z trzech skanów dla każdego stężenia. Mieszaninę przygotowano poprzez zmieszanie 300 ppm HCN w N2 z czystym N2 przy użyciu sterowników przepływu masowego i podgrzanie jej do 843 K. Centralny pik to linia P(20) pasma wibracyjnego ν1 HCN. Wstawka na Rysunku 3 pokazuje wartość szczytową sygnału z tej linii dla każdego stężenia, wraz z dopasowaniem wielomianowym drugiego stopnia. Zależność sygnału od stężenia można opisać wzorem S = ax2 + b, gdzie S to sygnał, a a i b to stałe dopasowania17. Bezwzględne pomiary stężenia w płomieniu wymagają pomiaru kalibracyjnego, jak pokazano tutaj, przy znanej temperaturze, w celu wyznaczenia stałej a. Należy również zmierzyć temperaturę w objętości pomiarowej w płomieniu, ponieważ stała a zmienia się wraz z temperaturą; pełna dyskusja na ten temat została już opublikowana17. Okres polaryzacji zastosowany w tym pomiarze wynosił 21,5 µm, przy temperaturze kryształu 104,5 °C.
Rysunek 4 przedstawia surowe dane z płomienia premiksowanego. Pokazuje on pięć kolejnych skanów w zakresie 3229,5-3232 cm-1, przy czym każdy skan trwał około 65 s. Obejmują one trzy grupy linii wody wykorzystywanych do pomiarów temperatury. W idealnych warunkach, przy pracy ze stabilnym układem, każdy skan w tym samym zakresie powinien być identyczny, ponieważ stężenie, ciśnienie i temperatura powinny pozostać niezmienne. Intensywność widocznych tutaj linii zmienia się gwałtownie z jednego skanu na drugi, co wynika z braku stabilności trybu i energii impulsu lasera pomiędzy skanami. Takie wyniki są nieużyteczne, chyba że zarejestrowano energię impulsu lasera, co pozwala oddzielić pomiary o wystarczającej energii od pozostałych. Okres polaryzacji zastosowany w tym pomiarze wynosił 21,5 µm, przy temperaturze kryształu 123 °C.
Na Rycynie 4 rozproszenie tła nie jest widoczne, ponieważ w celu zmniejszenia sygnału i uniknięcia nasycenia detektora zastosowano filtr ND2. W przypadku słabszych sygnałów ustalono, że rozproszenie tła wynosi około 5 pJ na impuls, co odpowiada sygnałowi generowanemu przez linię P(20) pasma wibracyjnego ν1 dla 100 ppm HCN w temperaturze pokojowej.

Rycina 1: Schemat detektora upkonwersyjnego. U1-U7 i UH to lustra z powłoką wysokoodblaskową (HR) dla 1,064 nm. Wszystkie lustra są płaskie z wyjątkiem U3, które ma promień krzywizny 200 mm. Lustra U1-U5 zostały wykonane tak, aby były przepuszczalne dla długości fali diody laserowej, co zapobiega przedostawaniu się światła LD do detektora. U6 jest przepuszczalne dla sygnału po upkonwersji w zakresie 650-1,050 nm. U7 jest przepuszczalne dla sygnału w średniej podczerwieni. Lustro UH wykazuje 95% odbicia dla 1,064 nm i 5% przejścia. Długość drogi od U1 do U3 wynosi 156 mm, a długość drogi od U3 do U7 wynosi 202 mm. L4 i L5 to soczewki achromatyczne o ogniskowych odpowiednio 60 mm i 75 mm. Obie są przezroczyste dla zakresu 650-1,050 nm. Kamera pełniąca rolę detektora jest umieszczona 75 mm od L5. Pole wnęki jest spolaryzowane pionowo. Zastosowany tutaj kryształ PPLN ma okresy polaryzacji 21,0 µm, 21,5 µm, 22,0 µm, 22,5 mm i 23,0 µm, a długość kryształu wynosi 20 mm. Jako detektor światła widzialnego i bliskiej podczerwieni wykorzystano kamerę UI-5240CP-NIR-GL firmy IDS Imaging Development Systems. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 2: Schemat układu DFWM. M1 to lustro dielektryczne o wysokim współczynniku odbicia (HR) dla długości fali lasera impulsowego. M2 to lustro dielektryczne z powłoką zapewniającą wysokie odbicie (HR) dla długości fali lasera impulsowego oraz przepuszczalność dla wiązki prowadzącej HeNe. M3-M6 to chronione lustra złote. B.C.1 i B.C.2 to płytki BOXCARS 1 i 2. L1 to soczewka CaF2 o ogniskowej 500 mm i średnicy 5,1 cm. L2 to soczewka CaF2 o ogniskowej 500 mm i średnicy 2,54 cm. L3 to soczewka CaF2 o ogniskowej 100 mm. Laser impulsowy jest spolaryzowany pionowo. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rycina 3: Sygnał dla różnych stężeń HCN w N2. Centralny pik to linia P(20) pasma wibracyjnego ν1 HCN. Wstawka przedstawia sygnał szczytowy dla każdego stężenia (znaczniki w kształcie rombów) wraz z dopasowaniem wielomianowym drugiego stopnia. Kliknij tutaj, aby wyświetlić powiększoną wersję tej ryciny.

Rysunek 4: Pięć kolejnych skanów o czasie trwania ok. 65 s każdy, wykonanych w płomieniu wymieszanym. Laser był skanowany w zakresie 3229,5-3232 cm-1. Widoczne tutaj piki stanowią sygnał z kilku grup linii przejść H2O. Sygnał został zredukowany za pomocą filtrów ND1 i ND0,6, aby uniknąć nasycenia detektora. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 5: Płytka BOXCARS widziana z boku. Jest to blok wykonany z przezroczystego materiału. Po stronie wejściowej połowa powierzchni pokryta jest powłoką antyrefleksyjną. Tutaj wpada wiązka lasera i dociera do strony wyjściowej, gdzie połowa powierzchni jest powlekana w celu uzyskania 50% transmisji. Światło odbite wewnętrznie w płytce ulega załamaniu do części strony wejściowej pokrytej powłoką o wysokim współczynniku odbicia, a następnie jest odbijane przez górną połowę strony wyjściowej. Powoduje to rozdzielenie jednej wiązki na dwie równoległe wiązki. Ten sam efekt można by osiągnąć za pomocą dzielnika wiązki i lustra, jednak dzielnik wiązki powodowałby pewne odbicia od tylnej powierzchni, co mogłoby zwiększyć szum tła. Ponadto płytka BOXCARS nie wymaga żadnej kalibracji, aby zapewnić równoległość dwóch wytworzonych wiązek. Prosimy kliknąć tutaj, aby wyświetlić powiększoną wersję tego rysunku.

Rysunek 6: Moc wewnątrzwnękowa w funkcji prądu pompowego diody laserowej dla modułu upkonwersji. Każdy punkt stanowi średnią mocy zmierzonej dla trzech oddzielnych ustawień wnęki, a słupki błędów wskazują rozrzut między poszczególnymi ustawieniami. Odchylenie od idealnej charakterystyki lasera jest spowodowane efektami termicznymi w krysztale laserowym oraz krysztale PPLN. Kliknij tutaj, aby wyświetlić powiększoną wersję tego rysunku.
Precyzja wyrównania impulsowej wiązki laserowej ma kluczowe znaczenie dla czułości metody. Należy zwrócić szczególną uwagę, aby upewnić się, że belki są oddzielone od siebie w równej odległości za tabliczkami BOXCARS i aby belki były równomiernie rozmieszczone wokół środka L1. Odchylenie od tego doprowadzi do znacznego spadku intensywności sygnału, a tym samym czułości. Podobnie należy zwrócić uwagę na to, aby wnęka modułu konwersji w górę pracowała w trybie podstawowym i aby wiązka sygnału była wyrównana w celu optymalnego nakładania się z pompą konwersji w górę. Sygnał można łatwo zmniejszyć o jeden lub dwa rzędy wielkości, jeśli wnęka konwersji w górę działa w niewłaściwym trybie lub nakładanie się wiązki sygnału na pole wnęki jest nieoptymalne. Obejmuje to umieszczenie L3 z milimetrową precyzją, tak aby ognisko wiązki sygnału znajdowało się w środku kryształu PPLN. Dzięki optymalnemu nakładaniu się i mocy 80 W wnęki możliwa jest 6% sprawność kwantowa stopnia SFG. Dzięki zastosowanemu tutaj detektorowi i długości fali całkowita skuteczność wykrywania wynosi 3%. Maksymalna moc wewnątrz wnęki, jaką można osiągnąć, wynosi 120 W, ale 80 W można osiągnąć niezawodnie. Wydajność konwersji jest proporcjonalna do mocy wewnątrzwnękowej, więc sygnały zarejestrowane z inną mocą wewnątrzwnękową można porównać, jeśli zarejestrowana jest moc wewnątrzwnękowa.
Głównym czynnikiem ograniczającym czułość tej metody jest rozpraszanie tła, które zagłusza słabe sygnały. Aby ograniczyć to rozpraszanie, bardzo ważne jest, aby optyka była wolna od kurzu, zwłaszcza soczewka L1. Należy również zwrócić uwagę, aby położenie bloku belki minimalizowało hałas tła. Blok belki powinien być umieszczony na stoliku xy, tak aby mógł być przesuwany w sposób kontrolowany zarówno w płaszczyźnie poziomej, jak i pionowej, prostopadłej do kierunku belek.
Omawiane tutaj skanowanie odbywa się za pomocą PPLN w stałej temperaturze. Wydajność konwersji jest proporcjonalna do(ΔkL/2π)2, gdzie Δk to niedopasowanie fazowe, a L to długość kryształu. Maksimum połowy pełnej szerokości (FWHM) tej funkcji to szerokość pasma detektora przy stałej temperaturze kryształu PPLN. FWHM tej funkcji zmienia się wraz z temperaturą kryształu i długością fali, ale na ogół jest rzędu 5 cm-1 w średniej podczerwieni, dla kryształu o długości 20 mm. Wyjątkiem jest okolica 4200 nm, gdzie szerokość znacznie wzrasta18.
Na schemacie konfiguracji na rysunku 2 nie uwzględniono optyki skalowania, ponieważ istnieje szereg kwestii, które należy rozważyć przed podjęciem decyzji, jakie skalowanie, jeśli w ogóle, jest potrzebne. W opisanej tutaj konfiguracji impulsowa wiązka laserowa jest kolimowana przy średnicy wiązki ok. 2 mm po osiągnięciu L1. Daje to talię wiązki w ognisku około 400 μm, przy długości fali 3 μm. Podczas stosowania tej techniki może być pożądana zmiana ogniskowej L1, ponieważ ze względów praktycznych potrzebna jest większa przestrzeń między L1 a punktem ogniskowym, albo w celu skrócenia objętości pomiaru poprzez zwiększenie kątów zbieżności, co można osiągnąć przy użyciu krótszej ogniskowej. W takim przypadku pas wiązki w punkcie ogniskowym powinien być utrzymany na poziomie ok. 400 μm, a skolimowana wiązka powinna być odpowiednio wyskalowana. Należy jednak wziąć pod uwagę, że zwiększenie średnicy belki bez zwiększania rozstawu belek zwiększy rozpraszanie od krawędzi bloku belki. Rozdzielczość przestrzenna jest określona przez zachodzenie na siebie belek pompy. W opisanej tutaj konfiguracji zakładka ma długość 6 mm, więc objętość pomiarowa to cylinder o długości 6 mm i promieniu 0,4 mm.
Aby osiągnąć quasi-fazowe dopasowanie w krysztale PPLN, zarówno sygnał średniej podczerwieni, jak i pole wewnątrzwnękowe wnęki konwersji w górę muszą być wyjątkowo spolaryzowane w krysztale PPLN. Wnęka konwersji w górę powinna być zbudowana tak, aby polaryzacja pola wewnątrzwnękowego była automatycznie prawidłowa. Jeśli laser średniej podczerwieni jeszcze do tego nie pasuje, można włożyć płytkę falową do wyjścia lasera w średniej podczerwieni, aby zmienić polaryzację.
IR-DFWM wymaga impulsów o stosunkowo wysokiej energii, 1-4 mJ, w połączeniu z wystarczająco wąską szerokością linii lasera, aby rozdzielić linie molekularne, które są rzędu 0,1 cm-1. Lasery, które spełniają te kryteria, mają na ogół niską częstotliwość powtarzania, a ponieważ pozyskiwanie danych za pomocą DFWM odbywa się zwykle poprzez skanowanie długości fali lasera, ogranicza to szybkość pomiarów. Oznacza to, że metoda ta jest najchętniej stosowana do pomiarów, w których obiekt nie zmienia się w czasie, chociaż została również zastosowana do pomiarów rozdzielczychw czasie 17. Innym ograniczeniem jest to, że ze względu na wrażliwość na rozproszone światło, cząstki znajdujące się w objętości pomiarowej lub w jej pobliżu będą tworzyć zdarzenia rozpraszające, które całkowicie zatopią sygnał17. Warunek dopasowania fazowego procesu konwersji w górę jest wąski spektralnie, co pomaga wyeliminować szum z termicznego promieniowania tła, ale sprawia, że skanowanie w szerokim zakresie długości fal jest bardziej czasochłonne, ponieważ temperatura PPLN musi być dostrojona, aby utrzymać długość fali sygnału dopasowaną do fazy.
Planowane są przyszłe zastosowania IR-DFWM w celu wykrywania NH3 w płomieniach lub kontynuowania prac z HCN w bardziej praktycznych środowiskach. Najbardziej oczywistym sposobem ulepszenia metody jest dalsza redukcja tła ze światła rozproszonego. Można to zrobić za pomocą przestrzennego filtrowania wiązki sygnału po odebraniu sygnału przez L2.
Autorzy nie mają nic do ujawnienia.
Wysoko cenione jest finansowanie otrzymane przez autorów w ramach programu Horyzont 2020 przez Unię Europejską. Prace te prowadzono w ramach innowacyjnej sieci szkoleniowej Mid-TECH Marie Curie [H2020-MSCA-ITN-2014-642661].
| Nazwa | Firma | Numer katalogowy | Komentarze |
|---|---|---|---|
| Laser Nd:YAG, pulsacyjny | Spectra Physics | Quantarau Pro-290-10 | Ilość: 1 (Do pompowania lasera mid-IR) |
| Laser Nd:YAG, system wysiewu iniekcyjnego | Spectra Physics | 6350 | Ilość: 1 |
| Laser barwnikowy NIR - system OPA | Sirah | OPANIR | Ilość: 1 |
| HeNe laser | Thorlabs | HNL100LB | Ilość: 1 |
| Lustro dichroiczne | LASEROPTIK | NA | Ilość: 1, na zamówienie: HR dla średniej podczerwieni, przezroczysty dla 632 nm |
| Chronione złote lustra | Thorlabs | PF10-03-M01 | Ilość: 5 |
| BoxCars Plate | LASEROPTIK | NA | Ilość: 2, Na zamówienie |
| xy-stage | Thorlabs | DTS25/M Ilość | : 2 |
| 500 mm ogniskowa obiektywu CaF2, Ø 2'' | Eksmaoptics | 110-5523E | Ilość: |
| obiektyw CaF2 o ogniskowej 1 500 mm, Ø 1' | 'Thorlabs | LA5464 | Ilość: 1 |
| obiektyw o ogniskowej 100 mm CaF2, Ø 1'' | Thorlabs | LA5817 | Ilość: 1 |
| Irys i Oslash; 50 mm | Thorlabs | ID50/M | Ilość: 1 |
| Filtr IR ND1 (ZnSe) | Thorlabs | NDIR10B | Ilość: 1 |
| Filtr IR ND2 (ZnSe) | Thorlabs | NDIR20B | Ilość: 1 |
| Filtr IR ND3 (ZnSe) | Thorlabs | NDIR30B | Ilość: 2 |
| Detektor konwersji w górę | NLIR | NA (najbliższy odpowiednik: U3055 3.0-5.5µ m) | Ilość: 1, na zamówienie niestandardowe |
| Karta detektora VIS/NIR | Thorlabs | VRC2 | Ilość: 1, (niska intensywność) |
| Karta detektora NIR | Thorlabs | VRC4 | Ilość: 1, (wysoka intensywność) |
| Karta detektora MIR | Thorlabs | VRC6S | Ilość: 1 |
| Głowica czujnika mocy cieplnej | Thorlabs | S302C | Ilość: 1 |
| Konsola miernika mocy | Thorlabs | PM100D | Ilość: 1 |
Poproś o pozwolenie na ponowne wykorzystanie tekstu lub ilustracji tego artykułu JoVE
Poproś o pozwolenie