2.2.1. Tham sô toơng quát
Các phađn tử C6H4F2, C6H4Cl2, C6H4Br2, C6H4I2, đeău có ba dáng lieđn kêt giữa các nguyeđn tử nhóm Halogen với vòng benzene là 1-4 (lieđn kêt đôi xứng), 1-3, 1-2 (Hình 2.3). Veă maịt lịch sử, đađy là những hợp chât đã được các nhà hóa hĩc tređn thê giới nghieđn cứu toơng hợp, phađn tích tính chât hóa lý moơt cách đaăy đụ từ 1991 trở veă trước. Các tráng thái vaơt chât cụa phađn tử và nhieơt đoơ chuyeơn pha được theơ hieơn tređn Bạng 2.1.
Bạng 2.1: Nhieơt đoơ chuyeơn các tráng thái raĩn – lỏng – khí cụa các chât [35]
Chât lieơu
Chuyeơn pha C6H4F2 C6H4Cl2 C6H4Br2 C6H4I2 Pha raĩn -> Pha lỏng 249,45 K 326,1 K 360 K 402 K Pha lỏng -> Pha khí 361 K 447,3 K 491,7 K 558,2 K
Do keđnh phađn tử yeđu caău lieđn kêt đôi xứng với đieơn cực D và S neđn chư có kieơu câu trúc 1,4 là thỏa mãn yeđu caău đaịt ra. Chính vì lý do đó, chúng ta tiêp túc đi sađu nghieđn cứu veă câu trúc này đeơ lây các tham sô cho chương trình mođ phỏng transistor phađn tử.
2.2.2. Phaăn meăm mođ phỏng hóa hĩc CAChe
Ngày nay, nhờ áp dúng tiên boơ cụa tin hĩc vào hóa hĩc lượng tử, các phađn tử vòng benzene lieđn kêt 1-4 với các nguyeđn tô thuoơc nhóm Halogen được tính toán các tham sô chính xác hơn, đaăy đụ hơn, mođ phỏng câu trúc 3D moơt cách trực quan hơn. Hieơn tái, tređn thê giới có nhieău phaăn meăm mođ phỏng hóa hĩc. Trong đeă tài này, tác giạ chĩn CAChe bởi đađy là phaăn meăm được Phòng thí nghieơm Cođng ngheơ nanođ cụa Đái hĩc Quôc gia Thành phô Hoă Chí Minh hieơn
đang sử dúng đeơ giạng dáy và nghieđn cứu. Phaăn meăm mođ phỏng hóa hĩc CAChe được xađy dựng từ 1989 bởi Oxford Molecular Limited. Từ naím 2000, Fujitsu Limited tiêp túc phát trieơn cho đên ngày nay.
Phaăn meăm mođ phỏng hóa hĩc CAChe cho phép thiêt kê các phađn tử và tôi ưu hóa câu trúc, kieơm tra lieđn kêt electron, xác định kích thước phađn tử, các mức naíng lượng và maơt đoơ mức trong dại. Đađy là những tham sô caăn thiêt đeơ giại phương trình dòng IDS, từ đó vẽ neđn đoă thị đaịc trưng dòng - thê cụa transistor phađn tử.
2.2.3. Các tham sô cụa phađn tử làm keđnh dăn xác định trong CAChe
Mở moơt cửa soơ ChemicalSample trong CAChe và dùng các cođng cú đeơ thiêt kê moơt vòng benzene roăi gaĩn hai nguyeđn tử cụa moơt trong sô các nguyeđn tô thuoơc nhóm Halogen vào vị trí 1,4. Dùng menu: Experiment chĩn New đeơ vào hoơp thối mới. Trong hoơp thối này, tiêp túc chĩn Property: Optimized geometry roăi nhâp chuoơt vào nút Start đeơ tôi ưu hóa câu trúc phađn tử và được kêt quạ như Hình 2.4. Các phađn tử đeău có tính đôi xứng.
Hình 2.4: Câu trúc phađn tử sau khi tôi ưu hóa a) C6H4F2, b) C6H4Cl2, c) C6H4Br2, d) C6H4I2
Trở lái menu View chĩn Show Electrons ta kieơm tra được đieău kieơn ba caịp đieơn tử lieđn kêt gép đođi moêi phía từ các nguyeđn tử thuoơc nhóm Halogen như Hình 2.5.
Từ menu: Experiment chĩn New và cháy chương trình tính các múc trong hoơp thối Property: UV-visible transitions, HOMO & LUMO, All molecular orbitals, van der Waals surface… Các kêt quạ được đieău khieơn tự đoơng nhớ vào file dữ lieơu trong cùng thư múc với file *.csf (phađn tử mà ta đã thiêt kê).
Rât nhieău tham sô Vaơt lý và Hóa hĩc cụa phađn tử được cung câp sau khi cháy chương trình trong hoơp thối Property. Đeơ tính dòng IDS cụa transistor phađn tử caăn xác định được vùng hóa trị, vùng câm thođng qua chư sô mức naíng lượng thâp nhât trong vùng dăn - Lowest Unoccupied Molecular Orbital (LUMO) và mức naíng lượng cao nhât trong vùng hóa trị - Highest Occupied Molecular Orbital (HOMO). Vùng dăn được xác định qua LUMO, LUMO+1, … LUMO+5 và các mức LUMO tiêp theo. Ngoài ra chúng ta caăn xem xét maơt đoơ mức trong moêi LUMO và sô mức trong đó. Tuy nhieđn, giới hán tređn cụa vùng dăn chính là mức 0 eV. Các mức naíng lượng cao hơn mức khođng ta coi như các đieơn tử (nêu có) trong đó là chuyeơn đoơng tự do đôi với phađn tử. Đaịc trưng IDS cụa transistor phađn tử khođng xét tới các mức naíng lượng > 0eV. Khi xem xét các tham sô mức
naíng lượng do CAChe cung câp ta thây có khạ naíng choăng chât các mức cụa các dại trong vùng dăn. Tuy nhieđn, sai sô này khođng ạnh hưởng leđn dáng đường đaịc trưng mà chư gađy neđn sai sô giá trị dòng bão hòa mà thođi.
Sau khi cháy chương trình với các múc caăn thiêt ta gĩi bạng sô lieơu từ menu “Analyze”. Bạng sô lieơu “Chemical Properties Spreadsheet” cung câp nhieău thođng tin veă phađn tử (Hình 2.7). Chĩn Molecular Orbitals tređn thanh cođng cú ở dưới màn hình ta có trị sô các mức naíng lượng tương ứng với sô đieơn tử đã có tređn đó. Hai mức naíng lượng lieđn tiêp đã có hai đieơn tử và chưa có đieơn tử chính là HOMO và LUMO. Chúng ta cũng có thođng tin này baỉng hình ạnh và sô lieơu nêu vào menu “Analyze” chĩn Show Surfaces. Trong hoơp thối Displayed Surfaces chĩn sô quỹ đáo tương ứng với mức caăn xem xét. Ở phađn tử C6H4Cl2, nêu chĩn MO22 ta được kêt quạ như Hình 2.8.
Hình 2.6: Hoơp thối đeơ cháy chương trình tính các mức naíng lượng trong phaăn meăm mođ phỏng hóa hĩc CAChe
Chọn “Atoms” tređn thanh cođng cú cụa bảng số liệu ta có thođng tin veă “LUMO Density”, đó chính là maơt đoơ mức do các nguyeđn tử trong phađn tử gađy ra.
Hình 2.8: Mức LUMO tređn ạnh phađn tử ở tráng thái tương ứng
Hình 2.7: Sô lieơu veă Molecular Orbitals cung câp cho thođng tin các mức naíng lượng. Bạng sô lieơu Atoms cung câp thođng tin veă maơt đoơ mức
Chọn “Sleter Basis Functions” tređn thanh cođng cụ của bảng số liệu ta có thođng tin về sô mức naíng lượng trong dại (Hình 2.9). Để đđo kích thước phađn tử
ngay trong giao diện của CAChe ta cần chạy chương trình ở Property: van der
Hình 2.9: Xác định sô mức naíng lượng trong dại
Waals surface và lấy kết quả từ menu “Analyze” với mục “Show Surfaces…”. Từ đađy, thực hiện thao tác đo ba kích thước dài, rộng, cao của phađn tử.
Kêt quạ đo cho thây có sự lượng tử hóa kích thước phađn tử (Hình 2.11). Chính vì vaơy, khi tính dòng IDS cụa transistor phađn tử caăn phại xét ạnh hưởng cụa hieơu ứng này leđn hĩ đáp trưng dòng - thê.
Baỉng cách làm tương tự như vaơy, ta có bạng sô lieơu cho cạ bôn phađn tử vòng benzene lieđn kêt 1-4 với bôn nguyeđn tô thuoơc nhóm Halogen (Fluorine, Chlorine, Bromine, Iodion). Nguyeđn tô thứ naím thuoơc nhóm Halogen là Astatine có tính phóng xá neđn khođng xem xét trong đeă tài này. Nêu coi kích thước phađn tử là giới hán beă maịt đẳng trị theo lực van der Waals thì cạ bôn phađn tử đeău có dáng hình đĩa. Maịt caĩt theo chieău roơng như moơt bođng hoa sáu cánh.
Bạng 2.2: Các tham sô phađn tử caăn thiêt cho vieơc tính dòng IDS xác định được qua phaăn meăm mođ phỏng hóa hĩc CAChe [4]
Phađn tưû
Tham sô C6H4I2 C6H4Br2 C6H4F2 C6H4Cl2 LUMO+4 2,361eV 2,216eV 2,113eV 2,292eV LUMO+3 -0,065eV -0,182eV 2,060eV 1,330eV LUMO+2 -0,140eV -0,202eV 2,042eV 1,294eV LUMO+1 -0,484eV -0,295eV -0,295eV -0,052eV LUMO -0,653eV -0,318eV -0,332eV -0,243eV HOMO -8,947eV -9,869eV -9,870eV -9,235eV Sô mức naíng
lượng / vùng dăn 49 49 24 24
Đoơ cao H 0,3530 nm 0,3803 nm 0,3839 nm 0,3366 nm Đoơ dài L 1,0321 nm 1,0025 nm 0,7957 nm 0,9340 nm Đoơ roơng W 0,6637 nm 0,6714 nm 0,6782 nm 0,6696 nm
2.2.4. Lieđn kêt cụa phađn tử với nguyeđn tử Vàng ở đieơn cực D và S
Các phađn tử C6H4F2, C6H4Cl2, C6H4Br2, C6H4I2 đeău có ba caịp đieơn tử lieđn kêt ghép đođi moêi phía. Ta có theơ dùng phaăn meăm hóa hĩc CAChe đeơ mođ phỏng sự lieđn kêt này và xem xét ạnh hưởng cụa các nguyeđn tử Vàng leđn các mức naíng lượng trong phađn tử làm keđnh dăn. Kêt quạ tređn Hình 2.12 cho thây, sau khi lieđn kêt với Vàng và tôi ưu hóa, câu trúc phađn tử làm keđnh dăn khođng bị biên dáng. Chính vì lý do đó, nêu ta kieơm tra beă maịt đẳng trị lực van der Waals thì kích thước ba chieău cụa phađn tử văn khođng thay đoơi. Sô lieơu tređn Bạng 2.3 cho kêt quạ là khođng tìm thây ạnh hưởng cụa các nguyeđn tử Vàng leđn mức naíng lượng cụa phađn tử làm keđnh dăn. Tređn cơ sở đó, vieơc laơp trình tính dòng IDS cụa transistor có theơ sử dúng các kêt quạ ở Bạng 2.2.
Hình 2.12: Phađn tử C6H4F2 lieđn kêt với Vàng ở ba hướng quan sát khác nhau Bạng 2.3: Kêt quạ tính toán cho thây các tham sô naíng lượng cụa phađn tử khođng
Kêt quạ thiêt kê và tính toán các tham sô kích thước, mức naíng lượng cụa các phađn tử C6H4Br2, C6H4F2, C6H4Cl2, C6H4I2 theơ hieơn qua hình ạnh và bieơu bạng được trình bày ở Phú lúc B.
2.3. SỬ DÚNG GUI TRONG MATLAB ĐEƠ MOĐ PHỎNG TRANSISTOR PHAĐN TỬ PHAĐN TỬ
2.3.1. Táo Slider đeơ nhaơp sô lieơu
Yêu tô quan trĩng nhât đeơ đánh giá khạ naíng ứng dúng cụa moơt transistor là xác định hĩ đaịc trưng dòng - thê theo phương trình
( DS)U const DS G U f I = = (2.2)
Ở transistor phađn tử, IDS còn phú thuoơc vào nhieơt đoơ, kích thước và các tham sô naíng lượng cụa keđnh dăn. Vì lý do đó, vieơc chĩn GUI trong MATLAB đeơ thiêt kê giao dieơn mođ phỏng là rât phù hợp. Các tham sô caăn khạo sát được xử lý thođng qua các Slider neđn vieơc thay đoơi thực hieơn rât nhanh chóng và chính xác. Giao dieơn mođ phỏng cho ta moơt cách nhìn trực quan hơn đeơ hieơu được những vân đeă thuoơc veă bạn chât, câu táo cụa transistor trong khi cođng ngheơ chê táo ra nó caăn phại moơt thời gian nữa. Có bôn Slider được táo ra trong giao dieơn mođ phỏng transistor phađn tử ở đeă tài này.
- Slider “Đoơ dài keđnh dăn L” cho phép thay đoơi kích thước theo đúng như sự lượng tử hóa tham sô này khi thực hieơn đo trong phaăn meăm mođ phỏng hóa hĩc CAChe.
- Slider “Đoơ roơng keđnh dăn W” đeơ thay đoơi kích thước roơng theo sáu mức như kêt quạ đo trong CAChe.
- Slider “Thê coơng Vg” đeơ chĩn đieơn áp đieău khieơn transistor nhaỉm vẽ neđn nhieău đường đoă thị táo neđn moơt hĩ đường đaịc trưng dòng - thê. Các MOSFET hieơn nay thường dùng đieơn áp đieău khieơn VG vào khoạng 0 đên 0,5V. Đeơ xem xét tính tương đương cụa nguoăn tín hieơu vào giữa MOSFET với transiator phađn tử, chương trình mođ phỏng chĩn giá trị VG thay đoơi được trong khoạng neđu tređn.
- Slider “Nhieơt đoơ T” đeơ khạo sát ạnh hưởng cụa nhieơt đoơ leđn đaịc trưng dòng – thê là tham sô quan trĩng nhaỉm đánh giá khạ naíng ứng dúng vào các mođi trường trong thực tê. Biêt được hieơu ứng nhieơt cụa transistor, người sử dúng deê dàng thiêt kê các mách đieơn lieđn quan nhaỉm tự đoơng oơn định các tham sô kỹ thuaơt cụa thiêt bị đieơn tử.
2.3.2. Thuaơt toán cụa chương trình tính dòng IDS
Đường đaịc trưng I-V cụa transistor trường trong bán dăn khôi được mođ tạ baỉng các phương trình khác nhau ở phaăn có đoơ dôc lớn và phaăn gaăn bão hòa. Các phương trình này có theơ được suy ra baỉng phương pháp khớp phím hàm từ đoă thị thực nghieơm.
Trong transistor phađn tử xét ở đeă tài này, đường đaịc trưng đươc vẽ leđn từ vieơc nôi các đieơm cụa kêt quạ giại moơt phương trình tính dòng IDS ứng với các tham sô câu trúc, nhieơt đoơ, VG, VD xác định (Hình 2.13).
- Tham sô VD xác định từ sô đieơm caăn tính tređn moơt đường đoă thị và khoạng cháy thođng qua dòng leơnh:
IV=63; %Tính 63 đieơm giá trị dòng IDS cho moơt đường đoă thị
VV=linspace (0,1.0,IV)%Có 63 giá trị VD cách đeău nhau trong khoạng 0V đên 1V.
Đeơ đường đoă thị được lieđn túc và giạm gâp khúc ta caăn khai báo taíng sô đieơm tính toán I-V. Khi sô đieơm tính toán taíng thì thời gian xử lý đeơ vẽ ra đường
đoă thị sẽ dài hơn, đoăng nghĩa với vieơc hieơn thị kêt quạ bị chaơm lái. Tuy nhieđn, đieău này còn phú thuoơc tôc đoơ xử lý cụa máy tính đang dùng.
a. Tính đieơn dung các tú đieơn
Từ tham sô kích thước dài L, kích thước roơng W lây từ các Slider tương ứng, kêt hợp với khai báo các haỉng sô vaơt lý ban đaău, đieơn dung tú CG được tính theo bieơu thức (2.1)
t WL CG 4ε0
=
Trong đó t là đoơ dày cụa lớp đieơn mođi SiO2 đóng vai trò cách ly keđnh dăn với cực đieău khieơn G. Đôi với transistor nanođ CS = CD = 0,05CG [34, tr.24]. Dòng leơnh đeơ tính đieơn dung cụa các tú đieơn trong chương trình:
Cg=epsr*eps0*W*L/t; %Tính đieơn dung tú CG
Cs=0.05*Cg; Cd=0.05Cg; %Tính đieơn dung tú CS và CD
b. Tính DO và NO
Sô tráng thái cụa keđnh phađn tử dựa vào bieơu thức (1.49) có theơ viêt dưới dáng:
D(E) = mcWL/πħ2 υ(E - EC) = D0 υ(E - EC)
Khi đó, D0 trong chương trình được tính bởi dòng leơnh:
D0=m*q*W*L/(pi*hbar*hbar);%Tính maơt đoơ tráng thái ban đaău
Khi tính được sô mức naíng lượng trong vùng dăn, maơt đoơ mức và giá trị LUMO từ phaăn meăm mođ phỏng CAChe thì dòng leơnh tính maơt đoơ tráng thái D sẽ baỉng tích cụa D0 với moơt ma traơn hàng có sô phaăn tử baỉng sô mức naíng lượng trong dại. Với dE đã biêt và giá trị VG cú theơ nào đó thì sô đieơn tử ở tráng thái ban đaău trong keđnh (khi VD = 0) được tính baỉng phép lây toơng (sum trong
MATLAB) cụa tích từng phaăn tử cụa ma traơn hàng D với hàm Fermi f0 thođng qua dòng leđnh:
N0=2*dE*sum(D.*f0);
Trong dòng leơnh này, ta đã tính cho cạ 2 đieơn tử có spin ngược chieău nhau tređn moêi mức naíng lượng.
c. Tính UL
Khi các tham sô do kích thước phađn tử được xác định ở moơt mức lượng tử nào đó thì các tú đieơn CS; CD; CG và CE là hoàn toàn xác định. Thê Laplace UL
cụa keđnh phađn tử là hàm cụa các biên VG và VD. UL được tính theo bieơu thức (1.34). ( ) ( D) E D G E G L qV C C qV C C U = − + − Trong đó, CE = CS + CG + CD.
Naíng lượng đieơn thê cụa keđnh phađn tử ở moơt thời đieơm là toơng cụa thê Laplace với sô háng tư leơ dăn đên sự thay đoơi sô đieơn tử trong keđnh như bieơu thức (1.37). N U N C q U U E L+ ∆ = ∆ = 0 2
Với q2/ CE = U0 là naíng lượng đieơn thê cụa keđnh phađn tử thay đoơi khi náp theđm moơt đieơn tử. Khi đieơn tử từ cực S boơ sung vào keđnh hay đieơn tử từ trong keđnh thoát ra cực D khođng baỉng nhau (do VG) thì giá trị naíng lượng đieơn thê cụa keđnh sẽ biên đoơng. Vòng laịp N → U → Nnew→ Unew như trình bày ở Múc 1.2.6.7 được thực hieơn qua đốn mã cụa chương trình dùng while trong MATLAB:
while dU>1e-6
f2=1./(1+exp((E+UL+U+ep-mu2)./kT)); N(iV)=dE*sum(D.*((f1.*g1/g)+(f2.*g2/g))); Unew=U0*(N(iV)-N0); dU=abs(U-Unew); U=U+0.1*(Unew-U); end
e. Tính IDS ứng với các tham sô có giá trị xác định truyeăn từ Slider sang
Khi đieơn thê cụa keđnh đụ đeơ đưa naíng lượng leđn theđm moơt mức mới thì dòng đieơn I được tính theo dòng leơnh:
I(iV)=dE*I0*(sum(D.*(f1-f2)))*g1*g2/g;
Ở đađy, phép lây tích phađn như bieơu thức (1.47) được thay thê baỉng phép lây toơng do dE cùng ma traơn D nhađn với hieơu sô hai hàm Fermi f1, f2 đã biêt. Các ký hieơu toán hĩc trong phaăn lý thuyêt transistor phađn tử trình bày ở Chương 1 phại thay đoơi khi viêt chương trình trong MATLAB là đeơ tránh sự trùng laịp với các hàm đã được MATLAB định nghĩa maịc nhieđn hoaịc khođng hieơu moơt sô ký tự Latin được dùng trong các phương trình đó.