Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
ISSN: 1412-0917
PROSES AKUSTIK DAN TERMAL DALAM SEL SAMPEL PADA SPEKTROSKOPI FOTOAKUSTIK Oleh: Andhy Setiawan, Asep Sutiadi, Parlindungan Sinaga Jurusan Pendidikan Fisika FPMIPA Universitas Pendidikan Indonesia ABSTRAK Sinyal fotoakustik (FA) yang dibangkitkan pada sel sampel dalam suatu sistem spektroskopi fotoakustik dapat diprediksi melalui persamaan-persamaan yang menggambarkan tekanan dan temperatur pada sel. Persamaan-persamaan yang berkaitan dengan proses akustik dan termal dalam sel FA telah diturunkan dari persamaan linear tergandeng untuk tekanan dan temperatur yang bergantung waktu dalam suatu fluida yang viskositasnya diabaikan. Dengan memasukkan syarat batas diperoleh 10 persamaan yang melibatkan ukuran sampel dan ukuran lain yang berkaiatan dengan posis sampel dalam sel, serta konstanta-konstanta yang berkaitan dengan persamaan akustik dan temperatur. Hasil yang diperoleh ini merupakan titik awal untuk analisis lebih lanjut mengenai sinyal FA. Kata kunci: akustik, termal, fotoakustik.
PENDAHULUAN Tujuan utama pengembangan model teoretis dari efek fotoakustik ialah agar dapat menginterpretasikan sinyal pada sistem spektroskopi fotoakustik (SFA) yang merupakan penyerapan optik oleh sampel. Sinyal akustik dalam suatu sel fotoakustik (FA) merupakan pemanasan periodik gas dalam ruang tertutup sebagai hasil penyerapan cahaya termodulasi oleh sampel. Variasi temperatur dalam lapisan tipis gas pada perbatasan sampel dengan gas mempengaruhi variasi tekanan (akustik). Oleh karena itu diperlukan telaah mengenai proses akustik dan termal dalam sel FA. Analisis mengenai sinyal FA yang lebih terperinci memerlukan persamaan distribusi tekanan dan temperatur. Persamaan umum untuk distribusi tersebut dapat diperoleh jika persamaan-persamaan yang melibatkan tekanan dan temperatur diturunkan dengan memasukkan syarat batas. Pada tulisan ini dijabarkan penurunan solusi harmonik untuk tekanan dan gas dengan memasukkan syarat batas yang diperoleh dari sel FA. Persamaan yang diturunkan ini diharapkan dapat digunakan untuk analisis lebih lanjut mengenai sinyal FA pada suatu sistem SFA.
50
ISSN: 1412-0917
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
EFEK FOTOAKUSTIK Efek fotoakustik (FA) dilandasi oleh adanya konversi radiasi elektromagnet (foton) menjadi gelombang bunyi (acoustic). Apabila radiasi yang dapat ditala pada frekuensi v dikenakan pada suatu sampel, maka sebagian molekul yang menghuni aras dasar akan tereksitasi. Molekul yang menyerap foton akan berubah dari keadaan dasar dengan energi E0 ke keadaan tereksitasi dengan energi E1, dengan perbedaan energi sebesar E1 E0 = hv dinyatakan dalam frekuensi foton v yang diserap. Molekul tersebut secara cepat atau lambat akan melepaskan energi dan kembali pada keadaan dasar dengan melalui cara seperti ditunjukkan oleh Gb. 1b (Rosencwaig, 1980; International School of Photonics, 2004). Sel Fotoakustik Berkas Laser Termodulasi
Gelombang Tekanan Gas Kontak (Udara)
Mikrofon (4x6x2mm) Sampel
Aras Eksitasi E1
Transisi Absorbsi Transisi Radiatif Non Radiatif Cahaya E0 Aras Dasar
(b)
(a) Gambar 1 Efek fotoakustik pada sampel
Transisi nonradiatif dapat menyebabkan terjadinya gelombang bunyi yang dapat dideteksi dengan mikrofon. Bertambahnya energi termal disebabkan molekul mengalami kenaikan energi kinetik rata-rata diikuti kenaikan suhunya. Bila intensitas radiasi foton yang datang dimodulasi dengan frekuensi tertentu maka akan menyebabkan suhu sampel termodulasi secara sepadan. Modulasi suhu sampel ini, dalam ruang tertutup berbentuk sel fotoakustik Gb. 1a (International School of Photonics, 2004), menghasilkan fluktuasi tekanan dengan frekuensi yang sama dengan frekuensi modulasi intensitas. Dalam eksperimen modulasi intensitas dilakukan secara mekanik (menggunakan chopper cahaya) dengan frekuensi chopping yang dapat diatur. Fluktuasi tekanan yang dihasilkan muncul sebagai gelombang bunyi yang dapat dideteksi dengan menggunakan sensor elektronik berupa mikrofon (Svanberg, 1990, h. 252). Proses yang dijelaskan di atas secara bagan ditunjukkan oleh Gb.2 (International School of Photonics, 2004)
Absorbsi Optik
Pemanasan Volume Teriluminasi
Ekspansi Termal
Gelombang Tekanan
Deteksi Akustik
Gambar 2 Evolusi sinyal fotoakustik
51
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
ISSN: 1412-0917
GELOMBANG AKUSTIK DAN TERMAL Setiap cahaya yang diserap oleh zat padat sebagian atau seluruhnya akan diubah menjadi panas melalui proses deeksitasi nonradiatif dalam zat padat. Sel pada Gb.3 (Rosencwaig dan Gersho, 1976) berbentuk silinder dengan diameter D dan panjang L. Diasumsikan bahwa panjang L lebih kecil dibandingkan dengan panjang gelombang sinyal akustik dan mikrofon akan mendeteksi tekanan rata-rata yang dihasilkan oleh sel. Selain itu, diasumsikan pula bahwa gas dan bahan backing tidak menyerap cahaya. Beberapa parameter yang terkait dengan persamaan aliran panas yaitu: ki merupakan konduktivitas termal bahan i (cal cm-1 sec-1 °C-1); i adalah kerapatan bahan i (g cm-3); Ci merupakan kapasitas panas bahan (cal g-1 °C-1), i k i i Ci adalah difusivitas termal bahan i (cm-2 sec-1); ai 2 i 1 2 adalah koefisien difusi termal bahan i (cm-1); i 1 i adalah panjang difusi termal bahan i (cm); i = s, g dan b untuk zat padat, gas dan bahan backing; ω merupakan frekuensi modulasi berkas cahaya datang (rad sec-1).
Sampel
Backing
Gas
Window Cahaya datang x
-(ls +lb)
- ls
0
lg
lg + lw
Gambar 3 Skema sel fotoakustik beserta koordinat yang digunakan
Ditinjau proses akustik dan termal satu dimensi dalam sel fotoakustik seperti yang ditunjukkan pada Gb.3. Penyerapan berkas cahaya termodulasi dianggap memberikan suatu sumber energi termal sesaat. Hal ini mensyaratkan penggunaan persamaan tergandeng untuk gerak gelombang akustik dan difusi termal dalam sampel dan gas, dengan sejumlah syarat batas yang benar. Difusi termal ditinjau pada jendela dan backing. Persamaan linear tergandeng untuk tekanan dan temperatur yang bergantung waktu dalam suatu fluida yang viskositasnya diabaikan ditulis sebagai 2 p
52
0 2 p B t 2
0
2 t 2
,
(1)
ISSN: 1412-0917
2
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005 p H x, t T0 , t t
(2)
dengan 0 dan T0 adalah kerapatan dan temperatur ambien, B adalah modulus bulk isotermal, adalah koefisien ekspansi termal, adalah konduktivitas termal, adalah difusivitas termal (= 0 C p , dengan Cp kapasitas panas pada tekanan tetap), dan H(x,t) (hanya ada pada sampel) mewakili sumber energi termal yang dihasilkan dari absorbsi optik. Jika intensitas cahaya monokromatik (dengan panjang gelombang ) yang datang I0 (dalam W/cm2) dimodulasi oleh chopper secara sinusoidal dengan frekuensi sudut , maka intensitas termodulasi dapat dinyatakan sebagai
I 12 I 0 1 cost . Rapat panas pada suatu titik x dalam sampel sebagai hasil penyerapan pada titik tersebut diberikan oleh
H x, t I e x , sehingga diperoleh H x, t 12 I 0 e x 1 cost ,
dengan adalah koefisien serapan optik dari sampel (dalam cm 1) pada panjang gelombang . Berdasarkan Gb.3 dapat diketahui bahwa x bernilai negatif karena berada pada rentang x = 0 sampai x = -ls dengan cahaya datang pada x = 0. Pada Gb. 3 tampak pula bahwa kolom udara menempati x = 0 sampai x = lg dan backing berada pada x = -ls sampai x = -(ls + lb). Untuk p dan yang memiliki kebergantungan waktu berbentuk e jt (variabel fisis diwakili oleh bagian real), maka pers. (1) dan (2) menjadi 2 p x 2
2 x
2
0 2
j
B
p 0 2 ,
T Ae x 1 e jt j 0 p,
(3)
(4)
dengan A I 0 2 . Gelombang bidang yang merupakan solusi persamaan tergandeng yang sesuai (untuk temperatur dan perpindahan) dengan zat padat isotropik memiliki bentuk yang sama dengan gelombang bidang solusi pers. (3) dan (4). Hal ini merupakan dasar pendekatan untuk penerapan pada sampel zat cair dan zat padat isotropik. Solusi pers. (3) dan (4) untuk sampel, dengan mengabaikan kebergantungan terhadap waktu, memiliki bentuk
p PA e jkx PA e x ,
(5)
53
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
ISSN: 1412-0917
T A e jkx T A e x .
(6)
Substitusi bentuk ini pada pers. (3) dan (4) menghasilkan empat persamaan faktor amplitudo, yaitu:
2 0 2 k PA 0 2T A 0 B jT0 2 j PA 0 k T A
(7)
2 0 2 PA 0 2T A 0 B jT0 2 j PA A T A
(8)
Kombinasi pers. (7) menghasilkan persamaan kuadratik dalam k2 yang memiliki penyelesaian 1 j 0 1 j 0 2 k2 2 B 4 B 2
1
2 j 3 0 T0 1 B 2
2
(9)
Nilai khas konstanta untuk gelombang akustik tak tergandeng dalam sampel adalah
k 2 k s2
0 2
(10)
B
sedangkan untuk “gelombang” termal k 2 s2
j
(11)
(gelombang termal ini melemah dengan cepat). Pada frekuensi yang biasa digunakan dalam SFA (10-1000 Hz), koreksi untuk ungkapan pendekatan ini diabaikan. Rasio amplitudo ( PA TA ) memiliki nilai yang nyata untuk setiap solusi k 2 . Untuk mode akustik, k 2 ks2 pada pers. (10) disubstitusikan pada pers. (7) menghasilkan d as PA TA a
54
, T0
(12)
ISSN: 1412-0917
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
dan untuk mode termal k 2 s2 pada pers. (11) disubstitusikan pada pers. (7) menghasilkan d ts PA T A t
(13)
j 0 .
Pers. (8) memiliki penyelesain T A
A
s2 2
,
(14)
PA d T A ,
(15)
dengan d
B 0 2 B 2 0 2
.
(16)
Selanjutnya, solusi harmonik dalam sampel adalah
s As e jks x Bs e jks x Cs e s x Ds e s x TA e x ,
(17)
ps d as As e jks x Bs e jks x d ts Cs e s x Ds e s x d T A e x .
(18)
Solusi pers. (3) dan (4) untuk gas (dengan mengabaikan bentuk S) diperoleh dengan cara yang sama. Gas dianggap mematuhi hukum gas ideal, sehingga B P0 (tekanan ambien), dan T0 1 . Karena pada gas tidak ada sumber maka pers. (3) dan (4) untuk gas dapat diungkapkan sebagai berikut.
2 p x 2 2 x
2
g 2 B
j
p g 2 ,
T j 0 p. g g
Gas dianggap mematuhi hukum gas ideal sehingga B P0 dan 1 T0 , dengan C p Cv , maka persamaan di atas menjadi
2 p x 2 2 x 2
g 2
j
B
p g 2 ,
j p. g g
(19)
(20)
55
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
ISSN: 1412-0917
Penyelesaian untuk pers. (19) dan (20) berbentuk
p PA e jkx ,
(21)
TA e jkx
(22)
Substitusi pers. (21) dan (22) ke dalam pers. (19) dan (20) menghasilkan dua persamaan untuk faktor amplitudo 2 2 k 2 g P g T 0; A A p0 T0 2 j k T A j PA 0. g g
(23)
Kedua persamaan tersebut digabungkan sehingga diperoleh persamaan dalam variabel k berikut:
j 0 2 k4 g p0
2 j 3 0 j 3 0 k 0 p0 g g
(24)
Solusi persamaan di atas berbentuk 1 j 0 2 k2 2 g p0
1 j 0 2 4 g p0
1
2
p g T0 j 3 0 0 g p0 g g T0
2 .
(25)
Dengan demikian solusi k 2 untuk mode akustik dan mode termal adalah k g2
g2
2 g P0 j
(26)
(27)
g
Rasio antara tekanan dan temperatur untuk masing-masing mode diperoleh dengan substitusi pers. (26) dan (27) pada pers. (23) menghasilkan
d sg
g g
d tg
,
(28)
j g g T0
Solusi harmonik untuk gas ialah
56
(29)
ISSN: 1412-0917
g Ag e
jkg x
pg d ag Ag e
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
Bg e jkg x
jkg x
Cg e
Bg e
jkg x
g x
Dg e
d C e tg
gx
g x
g
,
(30)
Dg e
gx
.
(31)
Distribusi tempereatur pada jendela dan backing dinyatakan dalam bentuk
Ce x Dex ,
(32)
yang berlaku baik untuk jendela maupun untuk backing, dengan menggunakan subskrip w dan b masing-masing untuk , , C dan D . Keduabelas koefisien As , Bs , dalam pers. (24), (25) dan (30)–(32) dapat ditentukan dengan menerapkan delapan syarat batas kontinuitas temperatur dan fluks panas, tiga syarat batas kecepatan, dan syarat batas kontinuitas tekanan. Syarat batas kontinuitas temperatur dan fluks panas meliputi
w 0 pada x l w lg , maka 0 Cw e
w g Cw e
w lw lg
D e w lw lg w
pada x lg , maka
w lg
Dw e
w lg
g w g g x x
w
(33)
Cw w k w e
w lg
Ag e
jkglg
Bg e
jkglg
Cg e
glg
Dg e
glg
pada x l g , maka
k w w Dw e
jk g Ag e
jkg lg
g s
pada x 0 , maka
jk g Bg e
jkg lg
w lg
g Cg e
g lg
g Dg e
g lg
(35) .
Ag Bg Cg Dg As Bs Cs Ds TA
g
g
x
s
s x
(36)
pada x 0 , maka
k g jk g Ag jk g Bg g Cg g Dg
k s jk s As jk s Bs s Cs s Ds TA .
s b
(34)
(37)
pada x ls , maka
As e jksls Bs e jks ls Cs e sls Ds e sls TA e ls Cb ebls Db e bls
(38)
57
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
s
s b b x x
ISSN: 1412-0917
pada x ls , maka
k s jk s As e jksls jk s Bs e jksls s Cs e sls s Ds e sls TA e ls
k b b Cb e
bls
b Db e
bls
.
(39)
b 0 pada x ls l b , maka Cb eb ls lb Db e b ls lb 0.
(40)
Syarat batas untuk kecepatan dengan v j 0 dp dx ,
vg 0 pada x lg , maka j 0
d ag jk g Ag e j kglg jk g Bg e j kglg d tg g C g e glg g Dg e glg 0.
vg vs
(41)
pada x 0 , maka
d ag jk g Ag jk g Bg d tg g Cg g D g
(42)
d as jk s As jk s Bs d ts s Cs s Ds d TA .
vs 0 pada x ls , maka
d a s jk s As e j ksls jk s Bs e j ksls d ts s Cs e sls s Ds e sls d TA e ls 0.
(43)
Syarat batas kontinuitas tekanan
p g ps
pada x 0 , maka
d ag Ag Bg d tg C g D g d as As Bs d ts C s Ds d TA .
(44)
Pers. (33) dan (40) membolehkan eliminasi dua koefisien Dw dan Cb. Substitusi Dw dari pers. (33) ke pers. (34) dan (35) menghasilkan
Cw e
w lg
1 e
w Cw k w e
2 w l w
wlg
A e
1 e
g
wl w
jkg lg
Bg e
jk
jkg lg
g Ag e
Cg e
jkg lg
g lg
jk g Bg e
Dg e jkg lg
g lg
g Cg e
(45) g lg
g Dg e
g lg
(46) Selanjutnya substitusi Cb dari pers. (40) ke pers. (38) dan (39) menghasilkan
As e jksls Bs e jks ls Cs e sls Ds e sls TA e ls Db e bls 1 e 2 blb
58
(47)
ISSN: 1412-0917
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
k s jk s As e jksls jk s Bs e jksls s Cs e sls s Ds e sls TA e ls
k b b Db e
bls
1 e
2 b l b
(48)
sehingga pada akhirnya diperoleh sepuluh persamaan linear tak homogen dengan sepuluh koefisien, yaitu pers. (45), (46), (36), (37), (47), (48) (41), (42), (43), dan (44). Hal ini menunjukkan bahwa pers. (24) dan (25) akan memberikan elemen-elemen tak homogen dari seperangkat persamaan-persamaan tersebut. Persamaan-persamaan itu dapat diselesaikan secara aljabar (diperkirakan hasilnya berbelit-belit) maupun secara numerik. Penyelesaian persamaan-persamaan di atas diharapkan dapat menunjukkan ketergantungan terhadap ruang (spasial) yang akan memperkenalkan beberapa parameter panjang yang bermanfaat untuk analisis lebih lanjut. Panjang difusi termal l merupakan jarak tempuh gelombang termal yang teredam dengan faktor e i . Jika ai Re i , maka l l-i 2 i i 2 dengan indeks i menunjukkan gas, sampel, dan lain-lain. Panjang serapan optis adalah jarak tempuh intensitas cahaya yang berkurang dengan faktor e i dan berlaku kaitan i (Sutiadi dkk, 2004). Seluruh perangkat persamaan tergandeng yang diperikan oleh pers. (24), (25) dan (30)(32) dapat digunakan untuk memprediksikan mengenai sinyal fotoakustik. Gambaran yang jelas mengenai sinyal FA dapat diperoleh jika persamaan-persamaan tersebut diselesaikan dengan memanfaatkan syarat batas yang diperikan oleh pers. (33)-(43). Serapan cahaya yang termodulasi menyebabkan pemanasan secara periodik dalam sampel yang menghasilkan perambatan variasi tekanan yang disebut dengan gelombang akustik. Superposisi gelombang-gelombang tersebut pada permukaan sampel membangkitkan gerak permukaan yang kemudian dipergunakan sebagai syarat batas untuk gelombang akustik dalam gas. Aliran panas secara periodik pada gas secara bersamaan menyebabkan ekspansi dan kontraksi dalam lapisan batas yang tipis pada sampel, yang oleh Rosencwaig dan Gersho (RG) telah dikarakterisasi sebagai piston bergetar. Jika piston termal ini dilapiskan ke atas permukaan mesin yang bergerak, maka piston komposit yang terbentuk akan menghasilkan sinyal yang terdeteksi dalam gas.
KESIMPULAN Distribusi temperatur dalam sistem sangat dipengaruhi oleh suku yang tergandeng pada tekanan yang ditunjukkan pada ruas kanan persamaan linear tergandeng untuk tekanan dan temperatur yang bergantung waktu dengan viskositas yang diabaikan. Persamaanpersamaan umum berkaitan dengan akustik dan temperatur dapat diperoleh dari persamaan tergandeng.Dengan memasukkan syarat batas diperoleh persamaan-persamaan yang melibatkan ukuran sampel dan ukuran lain yang berkaiatan dengan posis sampel dalam sel, serta konstanta-konstanta yang berkaitan dengan persamaan akustik dan temperatur. Gambaran yang jelas mengenai sinyal FA dapat diperoleh jika konstanta-konstanta yang terlibat pada persamaan-persamaan umum tersebut terlebih dahulu ditentukan, misalnya dengan cara numerik. Dengan kata lain persamaan-persamaan yang diperoleh berkaitan dengan proses akustik dan termal ini diperlukan untuk analisis lebih lanjut mengenai sinyal fotoakustik
59
Jurnal Pengajaran MIPA, Vol. 6 No. 1 Juni 2005
ISSN: 1412-0917
DAFTAR PUSTAKA International School of Photonics, 2004. Photoacoustics Effect [Online]. Tersedia: http://www.photonics.cusat.edu/Research_Photoacoustics.html [24 September 2004]. Rosencwaig, A. dan Gersho, A., 1976. Theory of Photoacoustic Effect with Solid, J. Appl. Physics, Vol. 47, No. 1, 64 – 69. Rosencwaig, A., 1980. Photoacoustic and Photoacoustic Spectroscopy, John Wiley and Sons, New York. Sutiadi, A. Sinaga, P., dan Setiawan, A., 2004. Kontribusi Efek Vibrasi Mekanik pada Bentuk Formulasi Efek Fotoakustik untuk Menentukan Pengaruh Parameter Sel Fooakustik terhadap Sinyal Fotoakustik, Laporan Hibah Penelitian Program SP-4 UPI, Program Studi Fisika FPMIPA UPI, Bandung Svanberg, S., 1991. Atomic and Molecular Spectroscopy, Springer-Verlag, Berlin, Germany
60