Hasil Pencarian  ::  Simpan CSV :: Kembali

Hasil Pencarian

Ditemukan 85835 dokumen yang sesuai dengan query
cover
Efta Yudiarsah
"Kerapatan keadaan vibrasi lokal dan fungsi korelasi pada beberapa struktur ikatan di dalam silikon karbon amorf (a-Si1-xCx:H) telah hitting dengan menggunakan model potensial Born dan metode "cluster" kisi Bethe. Mode-mode vibrasi lokal yang disebabkan oleh unit monohidrit (SiH dan CH), dihidrit (SiH2 dan CH2), dan trihidrit (SiH3 dan CH3) telah diidentifrkasi. Studi dilakukan pada dua jenis distribusi atom-atom penyusun bahan amorf yaitu "random sequence" dan "chemically ordered". Skema interpolasi digunakan untuk mendiskripsikan probabilitas ikatan dengan rasio percabangan pohon Cayley. Mode-mode vibrasi yang disebabkan keberadaan atom H di dalam pita fonon bulk silikon karbon sangat dipengaruhi oleh jenis atom tetangga terdekat kedua dari atom H. Mode-mode vibrasi yang disebabkan keberadaan atom H di luar pita fonon bulk silikon karbon hampir tidak dipengaruhi oleh jenis atom tetangga terdekat kedua dari atom H Hasil perhitungan yang diperoleh sesuai dengan data-data eksperimen.

Vibration Local Modes of Amorphous Silicon Carbon (a-Si1-xCx:H): Theoretical Study by Using Cluster Bethe Lattice Methods. Local vibrational density of states and correlation function for various bonding structures in amorphous silicon carbon (a-Si1-xCx:H) have been calculated by employing Born potential models and cluster Bethe lattice method . The local modes induced by the monohydride (SiH and CH), dihydride (SiH2 and CH2), and trihydride (SiH3 and CH3) units are identified The study has been made for two types of distributions of constituent atoms of the alloy: a random sequence and a chemically ordered The interpolation scheme can relate the bonds probabilities to the branching ratios of a Cayley tree. The vibration modes induced by H atoms inside the bulk phonon region of an amorphous silicon carbon (a-Si1-xCx) are influenced very much by the presence of the kind of atoms lying on the next-nearest-neighboring sites of the H atom. On the other hand, the modes induced by H atoms lying outside the bulk phonon region of an amorphous silicon carbon (a-Si1-xCx) are nearly undisturbed by the presence of the kind of atoms lying on the next-nearest neighboring sites of the H atom. The calculated results are in a good agreement with the available experimental data."
2000
T10328
UI - Tesis Membership  Universitas Indonesia Library
cover
Ayu Citraningtyas
"ABSTRAK
Material semikonduktor Cu2ZnSnS4 (CZTS) dikenal sebagai semikonduktor tipe-p dengan energi celah pita ideal dan koefisien penyerapan tinggi untuk lapisan penyerap pada aplikasi sel surya. Fabrikasi menggunakan metode Successive Ionic Layer Adsorption and Reaction (SILAR) menjadi kombinasi yang tepat untuk menghasilkan sel surya berbasis lapisan tipis yang terjangkau dan rendah toksisitas. Siklus pencelupan menjadi salah satu faktor penting proses yang dapat mempengaruhi struktur dan sifat optis lapisan tipis CZTS yang terbentuk pada permukaan substrat. Dengan menggunakan variabel 25, 30, 35, dan 40 siklus, serta perlakuan anil tanpa dan dengan suasana sulfur, penelitian ini melakukan investigasi pengaruhnya terhadap struktur dan sifat optis berupa nilai energi celah pita. Hasil XRD menunjukkan penurunan nilai kristalinitas dengan kenaikan jumlah siklus pencelupan. Topografi permukaan lapisan tipis CZTS hasil SEM menunjukkan adanya retak dan penggumpalan partikel pada permukaan sampel yang diduga sebagai fasa kedua berdasarkan analisis hasil EDX. Nilai energi celah pita pada sampel hasil anil tanpa suasana sulfur dan sampel hasil anil dalam suasana sulfur pun mengalami penurunan seiring dengan peningkatan jumlah siklus pencelupan.

ABSTRACT
Semiconducting Cu2ZnSnS4 (CZTS) material is known as p-type semiconductor which has ideal direct band gap and high absorption coefficient for absorber layer in thin-film solar cells application. Fabricated by Successive Ionic Layer Adsorption and Reaction (SILAR) method, this could be a promising technique to produce low cost and low toxicity thin-film solar cells. Immersion cycle is one of the important factors in SILAR method that may effect on structure and optical properties of CZTS thin film. By using following variables: 25, 30, 35, and 40 immersion cycles, and annealing treatment in non-sulfur condition and annealing treatment in sulfur condition as well, this investigation focused on their effects to structure and optical properties. The XRD results give decreased crystallinity with the increasing of immersion cycles. Surface topography of CZTS thin film, as the results of SEM examination, indicate the presence of cracks and coalescence particles on the surface of samples, suspected as second phases according to the results of EDX examination. Besides, as the immersion cycles are going up, it leads to decreasing on band gap energy on both annealed samples."
Fakultas Teknik Universitas Indonesia, 2014
S53881
UI - Skripsi Membership  Universitas Indonesia Library
cover
Pulung Karo Karo
"Telah dibuat lapisan tipis CdS dengan metode koevaporasi CdS dan S dan telah diuji dengan XRD, UV-VIS Spectrophotometer, dan pengukuran hambatan. Lapisan yang terbentuk dipengaruhi beberapa parameter seperti, temperatur substrat, jarak antara sumber dan substrat, dan jarak antara sumber CdS dan sumber S. Hasil XRD menunjukkan bahwa lapisan yang terbentuk adalah CdS dengan preferred orientation pada bidang (0002}, dengan struktur hexagonal. Perlakuan panas yang diberikan pada temperatur 200°C di ruang vakum akan menurunkan besar hambatan, dan menaikkan respon terhadap cahaya."
Depok: Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Indonesia, 1998
T-Pdf
UI - Tesis Membership  Universitas Indonesia Library
cover
Engelin Shintadewi Julian
"Jika dua bahan semikonduktor yang tidak sama digabungkan, sambungannya disebut heterojunction. Karena kedua bahan semikonduktor tersebut mempunyai bandgap yang berbeda, akan terjadi diskontinuitas energi pada persambungannya. Sambungan yang terbentuk dapat berupa sambungan abrupt atau graded. Sambungan abrupt terjadi jika perubahan dari satu bahan ke bahan lainnya terjadi secara tiba-tiba, sedangkan sambungan graded terjadi jika perubahan dari bahan satu ke bahan lainnya terjadi secara perlahan-lahan.
Pada penelitian ini dibuat rancangan HBT (Heterojunction Bipolar Transistor) Si/SiGe sambungan abrupt yang dapat memberikan frekuensi transit (fT) dan penguatan arus yang optimum. Untuk memperoleh frekuensi transit tersebut digunakan model struktur divais dengan NE=I01 S cm-3, NB=5.1019cm"3, dan WE=-50 nm. Lebar basis W8 dibuat bervariasi antara 10 - 50 nm dengan kenaikan 10 nm dan fraksi mol Ge dibuat bervariasi antara 0.15 sarnpai 0.25 dengan kenaikan 0.025. Dalam perancangan digunakan model HBT yang memasukkan mekanisme drift, difusi dan emisi termionik. Model tersebut juga digunakan untuk menentukan densitas arus basis, arus kolektor dan penguatan arus untuk lebar basis dan fraksi mol Ge yang bervariasi.
Hasil perhitungan menunjukkan bahwa pada HBT yang diteliti komponen arus emisi termionik sekitar 2.8 sampai 4.2 kali lebih besar dibanding dengan anus drift-difusi. Pengurangan lebar basis dari 50 nm menjadi 10 nm pada HBT dengan N E.--1018 cm-3, NB=5.109 cm 3, WE=50 nm, dan x .25 dapat meninakatkan frekuensi transit dari 52 GHz menjadi 178 GHz, sedangkan besarnya penguatan arus untuk WB=1O nm adalah 312 kali.

When two different semiconductor materials are used to form a junction, the junction is called a heterojunction. Since the two materials used to form a heterojunction will have different bandgaps, the energy band will have a discontinuity at the junction interface. The type of the junction could be an abrupt or graded junction. In abrupt junction, the semiconductor changes abruptly from one material to the other material. on the other.hand in graded junction, the composition of materials are graded.
In this work, we designed abrupt junction Si/SiGe HBT (Heterojunction Bipolar Transistor) with optimum transit frequency VT) and current gain. To achieve that transit frquency, the divals structure model with NE=1418 cm 3, Na 5.1019cm 3, and WE-50 nm is used. Base width WS was changed between 10 - 50 run with 10 nm increment and the Ge mole fraksi was between 0.15 - 0,25 with 0.025 increment. The model of HBT which include the drift, difussion and thermionic emission is used. The HBT model is also used to calculate the base and collector current density and the current gain for variable base width and Ge mole fraction.
The result show that the thermionic emission current is 2.8 until 4.2 times higher than the drift-difussion current. If base width is decreased from 50 nm to 10 nm of an HBT with NE=1018 cm-3, N8=5.1019 cm-3, WE=50 nm, and x=0.25 , the transit frequency is increased from 52 GHz to 178 GHz. The current gain for WB= 10 nm is 312 times larger.
"
Depok: Fakultas Teknik Universitas Indonesia, 1999
T-Pdf
UI - Tesis Membership  Universitas Indonesia Library
cover
Seno Aji
"Sifat ferromagnetik yang muncul pada DMS berbasis ZnO dengan Tc di atas suhu kamar memiliki prospek yang baik untuk aplikasi spintronik. Oleh karena itu, diperlukan model yang dapat memprediksi dengan akurat sifat-sifat fisis material tersebut. Model yang digunakan pada penelitian ini mengasumsikan bahwa atom-atom impuritas (atom Mn) yang diberikan pada sistem ZnO membentuk momen magnetik lokal. Sifat magnetik muncul karena adanya magnetic exchange interaction antara spin elektron bebas dan momen magnetik lokal atom-atom impuritas tersebut.
Hamiltonian model yang digunakan terdiri atas suku kinetik yang berasal dari pendekatan k.p 8 band dan suku interaksi magnetik antara spin elektron bebas dan momen magnetik lokal. Model ini diselesaikan dengan metode Dynamical Mean Field Theory (DMFT). Perhitungan dilakukan dengan memberi variasi terhadap kopling Hund JH sebesar 1 eV, 2 eV, dan 3 eV, variasi konsentrasi ion dopan x sebesar 2 %, 5 %, 7 %, dan 10 % untuk masing-masing JH serta variasi temperatur T untuk masing-masing konsentrasi. Hasil perhitungan menunjukkan munculnya impurity band pada kurva densitas keadaan (DOS) yang disebabkan oleh interaksi magnetik antara spin atom-atom impuritas dan spin elektron bebas yang menempati pita valensi. Nilai Tc tertinggi yang diperoleh menggunakan model ini yaitu sekitar 442 K dengan konsentrasi hole sebesar 1:71020 cm3.

Ferromagnetic properties that appear in ZnO-based DMS with Tc above room temperature are promising for spintronics applications. To make accurate predictions on other important physical properties, a good model for these materials is needed. The model used in this study assumes that impurity atoms (Mn atoms) doped in the ZnO system form local magnetic moments. The magnetic properties in this model arise due to magnetic exchange interaction between spins of electrons and local magnetic moments of the impurity atoms.
The model Hamiltonian consists of a kinetic term derived from k.p 8-band approximation and the magnetic exchange interaction term arising from the magnetic interactions between spins of electrons and local magnetic moments. The model is solved using Dynamical Mean Field Theory (DMFT) method. In this calculations, we use the following parameter values: Hund's coupling JH about 1 eV, 2 eV, and 3 eV, concentration of ion dopant x about 2%, 5%, 7 %, and 10 % for each JH, with variation of temperature for each concentration. The impurity bands appear in our density of states (DOS) calculations due to the magnetic interaction between spins of the Mn atoms and spins of the electrons occupying the valence bands. The highest Tc obtained using this model is about 442 K with the hole concentration of 1:71020 cm3.
"
Depok: Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Indonesia, 2013
S47599
UI - Skripsi Membership  Universitas Indonesia Library
cover
Altifani Rizky Hayyu
"ABSTRACT
Current electronics technology still relies on semiconductors for several purposes. The special property that distinguishes a semiconductor from another material is band gap energy inside of which the Fermi energy is located. Most conventional semiconductors such as Ge, Si, or GaAs may be classified into none or weakly correlated systems in which electron electron e e interactions do not play any significant role, where the band gaps insensitive to temperature change. Meanwhile, in semiconductors containing transition metal elements having d orbitals in their valence and or conduction band s , which we consider as strongly correlated semiconductors, e e interactions may play a more significant role. We hypothesize that such kind of semiconductors would have band structures, including their band gaps, being rather sensitive to temperature change due to e e interactions. We propose to explore this hypothesis theoretically by modeling the semiconductor band structure through Tight Binding Approximation with the GW method numerically in the Matsubara frequency domain. At the end, we use Pad approximant to obtain the self energy defined in the real frequency domain. Using this self energy we can calculate and analyze the Density Of States DOS at various temperatures, by giving the certain range of bare interaction line V. Our research results the correlation effects become stronger as we decrease the temperature. In short ranged interaction confirms that the semiconductor band gap increases and chemical potential shifts to a higher energy. Whereas with long ranged interaction shows that the range of two bands semiconductor becomes farther apart as compared to the bare DOS.

ABSTRACT
Teknologi elektronik saat ini masih mengandalkan material semikonduktor untuk berbagai kebutuhan. Properti yang membedakan semikonduktor dari material lain adalah celah pita energi di mana energi Fermi berada di dalamnya. Umumnya semikonduktor seperti Ge, Si, atau GaAs dapat diklasifikasikan ke dalam sistem non terkorelasi / terkorelasi lemah di mana interaksi elektron-elektron e-e tidak memiliki peran penting, sehingga pita energi nya mungkin tidak sensitif terhadap perubahan suhu. Sedangkan semikonduktor yang mengandung unsur logam transisi dan orbital d dalam pita valensi dan / atau pita konduksi dianggap sebagai semikonduktor terkorelasi kuat, dimana interaksi e-e memainkan peran lebih signifikan. Kami berhipotesis bahwa jenis semikonduktor ini akan memiliki lebar celah pita yang agak sensitif terhadap perubahan suhu. Kami mengeksplorasi hipotesis ini secara teoritis dengan pemodelan struktur pita semikonduktor melalui pendekatan Tight-Binding, kemudian menerapkan interaksi e-e dalam metode GW secara numerik dalam domain frekuensi Matsubara. Setelah itu kami menggunakan Pad approximant untuk memperoleh self-energy yang didefinisikan dalam domain frekuensi riil. Dengan menggunakan self-energy ini kita dapat menghitung dan menganalisis Density Of States DOS seiring dengan bertambahnya suhu, dengan memberikan variasi rentang dari bare interaction line V . Penelitian kami membuktikan bahwa efek-efek korelasi menjadi lebih kuat seiring menurunnya suhu. Untuk variasi interaksi rentang pendek, pita energi melebar dan potensial kimia semakin bergeser ke tingkat energi lebih tinggi. Sedangkan pada variasi interaksi rentang panjang, jarak antara dua pita energi semikonduktor semakin melebar jika dibandingkan dengan DOS murni."
2017
S67535
UI - Skripsi Membership  Universitas Indonesia Library
cover
Asep Budiyanto
"ABSTRAK
Study Magnetoimpedance material nano soft magnetic Fe-(Al,Ga)-(P,C,B,Si),
yang merupakan material amorf berbasiskan Fe dalam bentuk pita dengan lebar 4
mm dan tebal 40 µm yang dipreparasi dengan metode melt-spinning dan diberikan
perlakuan panas melalui teknik joule heating dilakukan dengan tujuan untuk
mengetahui hubungan magnetoimpedance dari material ini dengan struktur mikronya
baik sebelum ataupun setelah mengalami perlakuan panas serta bagaimana
pengaruhnya terhadap sifat-sifat magnetik yang sudah ada.
Karakterisasi struktur terhadap sampel dari material ini menunjukkan bahwa
perlakuan annil melalui teknik joule heating sampai 1.6 A pada sampel menyebabkan
terjadinya relaksasi struktural yang membuat penurunan anisotropi dan berakibat pada
peningkatan sifat kemagnetan lunak dari sampel.
Pengaruh medan magnet luar terhadap sampel pada frekuensi rendah
menyebabkan terjadinya pengurangan nilai impedansi dari sampel."
2007
T21283
UI - Tesis Membership  Universitas Indonesia Library
cover
Janice Aphrodite Tombeg
"Untuk memenuhi keperluan teknologi termoelektrik, para ilmuwan mengembangkan perangkat termoelektrik yang dapat mengubah energi panas menjadi energi listrik. Salah satu material yang digunakan dalam penelitian ini adalah semikonduktor Half Heusler. Beberapa penelitian sebelumnya menunjukkan bahwa material ini adalah kandidat yang baik untuk perangkap termoelektrik. Material ini memiliki stoikiometri XYZ, dimana X dan Y merepresentasikan metal transisi dan Z merepresentasikan unsur golongan utama. Penelitian ini dilaksanakan dengan metode komputasi menggunakan Quantum ESPRESSO untuk membentuk grafik struktur pita dan Density of states dari material ini dan BoltzTraP untuk menghitung properti termoelektrik dari material 12 jenis material Half Heusler. Keduabelas material ini memiliki celah pita yang menunjukkan material ini adalah semikonduktor, yang merupakan kandidat yang baik untuk menghasilkan konduktivitas elektrik dan koefisien Seebeck yang tinggi, serta nilai figure of merit (ZT) yang mengkarakterisasi efisiensi perangkat termoelektrik. Selain itu, penelitian kami menunjukkan bahwa semikonduktor Half Heusler lebih baik digunakan untuk alat termoelektrik tipe p karena nilai properti termoelektrik pada area tipe p lebih tinggi dari area tipe n.

For providing the need of thermoelectric technology, scientists develop thermoelectric device that can convert heat energy to electrical energy. One of the material that can be used in this research is Half Heusler semiconductor. Recent studies shown that this material is a good candidate to be a material for thermoelectric device. This material has a stoichiometry of XYZ, where X and Y represent transition metal and Z is a main group element. This research use a computational method using Quantum ESPRESSO to generate the graph of Band structure and Density of states of the material and BoltzTraP to calculate the thermoelectric properties of 12 kind of half-Heusler materials. All of these 12 material have band gap that indicate that they are asemiconductors, which are good candidates for high value good electrical conductivity, Seebeck coefficient and of figure of merit (ZT) that characterize the efficiency of a thermoelectric device. In addition, our result indicates that Half Heusler semiconductor is better to be used for p-type thermoelectric devices because the thermoelectric properties in p-type region is higher than in n-type region."
Depok: Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Indonesia, 2021
S-pdf
UI - Skripsi Membership  Universitas Indonesia Library
cover
Imam Agung Setiawan
"Telah berhasil dilakukan sintesis semikonduktor Cu2SnS3 (CTS) sebagai absorber semikonduktor yang dilakukan dengan metode fisik menggunakan mesin high energy planetary ball mill dan bola-bola ZrO2 (zirconia) dengan variasi waktu milling dan komposisi unsur. Hasil milling kemudian dilakukan annealing pada temperatur 450oC selama 30 menit. Kehadiran fasa CTS meningkat pada setiap peningkatan waktu milling, dengan waktu milling terlama selama 240 menit menghasilkan pola difraksi yang paling mendekati referensi dan paling tinggi dari semua sampel. Penambahan unsur Sn pada sampel Sn-rich berpengaruh terhadap kehalusan hasil uji difraksi, akan tetapi mengurangi fasa CTS yang terbentuk dan meningkatkan pembentukan senyawa pengotor. Nilai energi celah pita yang lebih tinggi pada sampel Sn-rich bukan dihasilkan oleh absorbansi CTS, melainkan oleh fasa-fasa pengotor yang diduga memiliki kemampuan sebagai absorber sebagai hasil dari proses fisik dan penambahan unsur Sn pada sampel uji. Didapatkan hasil uji X-ray Diffraction yang tinggi serta nilai energi celah pita sebesar 1,1-1,2 eV yang sesuai dengan literatur. Didapatkan hasil uji scanning electron microscope yang menandai terjadinya penyatuan butir-butir CTS dan voids sebagai hasil dari proses fisik dan annealing.

The semiconductor syntesis of Cu2SnS3 has been succesfully carried as a semiconductor absorber by physical method using high energy planetary ball mill machine and ZrO2 (zirconia) balls with milling time and powder composition variations. The results from this process then annealed in an electric furnace in 450oC for 30 minutes. The presence of CTS phase increased simultaneously with the increase of milling time, with the longest milling time 240 minutes produced the closest and the highest x-ray diffraction pattern compared to the literatures of all samples. The addition of Sn element in Sn-rich samples affected the refinement of x-ray diffraction patterns, however reduced the CTS phases formed and shown increased impurity phases. The higher value of bandgap energies from Sn-rich samples may be not produced by the absorbance of CTS, but of the impurities as the other results from physical process and the addition of Sn element to the test samples. X-ray diffraction test confirm the high peak of CTS with 1,1-1,2 eV bandgap energy. Scanning electron microscope showed fusion of CTS grains and voids as the result of physical and annealing process."
Depok: Fakultas Teknik Universitas Indonesia, 2019
S-Pdf
UI - Skripsi Membership  Universitas Indonesia Library
cover
Dita Oktri
"Perilaku magnetik BaFe12O19 (BaM) berbentuk hexagonal dalam berbagai ukuran telah diinvestigasi dengan menggunakan perangkat lunak simulasi mikromagnetik untuk umum, OOMMF berbasis persamaan LLG. BaM berbentuk hexagonal dengan ukuran 50 nm hingga 500 nm disimulasi dalam keadaan tanpa medan luar pada temperatur 0 K dan konstanta redaman α = 0.1. Struktur domain yang dibentuk oleh sistem pada ukuran sampai dengan 420 nm berupa struktur singledomain. Saat ukuran BaM 430 - 500 nm, struktur domain yang ditunjukkan bukan lagi single-domain melainkan two-domain. Menariknya, Bloch wall muncul sebagai dinding yang memisahkan dua domain dengan orientasi magnetik antiparallel pada BaM dengan struktur two-domain tersebut. Simulasi BaM berbentuk hexagonal dengan ukuran 50 – 100 nm dalam keadaan diberi medan quasi-static ke arah hard- dan easy axis juga dilakukan guna memahami pengaruh ukuran terhadap kurva histeresis.
Hasilnya, kurva histeresis memperlihatkan koersivitas nol saat medan eksternal diberikan ke arah hard-axisnya dan koersivitas yang besar (>375 mT) saat medan eksternal diberikan ke arah easyaxisnya yangmana merupakan tipikal kurva histeresis material hard magnet dengan uniaxial axis. Selain itu, medan nukleasi dan medan saturasi menunjukkan adanya pergeseran seiring dengan berubahnya ukuran BaM. Perilaku magnetik bergantung pada ukuran (size-dependent) yang ditunjukkan pada hasil penelitian ini dapat dijelaskan sebagai kompetisi antara interaksi magnetostatik, exchange, dan anisotropi yang direpresentasikan oleh densitas energi-energinya.

Magnetic behavior of hexagonal-shaped BaFe12O19 (BaM) with various size has been investigated by using public micromagnetic simulation software, OOMMF based on LLG equation. Hexagonal-shaped BaM with size ranging from 50 to 500 nm has been simulated under no applied field at 0 K with damping constant α = 0.1. BaM with size up to 420 nm formed single-domain structure. As the size increased at 430 - 500 nm, domain structure showed by BaM was no longer single-domain but two-domain. Interestingly, Bloch-type domain wall appeared to separate two domains with antiparallel orientation in this two-domain state. Micromagnetic simulation of hexagonal-shaped BaM with size ranging from 50 to 500 nm under quasi-static field parallel and perpendicular with its hard axis was also carried out in order to understand the influence of size to hysteresis curve of BaM.
The results showed hysteresis curve with zero coercivity when the external field was applied to its hard axis and large coercivity (>375 mT) as the external field was applied to its easy axis which is typical of hard magnet material with uniaxial axis. Moreover, nucleation and saturation field were shifted as the BaM’s size varied. Size-dependent magnetic behavior showed in these results could be explained as the competition between magnetostatic, exchange, and anisotropy which are represented by its energy densities.
"
Depok: Fakultas Matematika dan Ilmu Pengetahuan Alam Universitas Indonesia, 2013
S46634
UI - Skripsi Membership  Universitas Indonesia Library
<<   1 2 3 4 5 6 7 8 9 10   >>