DFT-modeling of Ammonia Molecules Protonation on a p-type Silicon Surface

Abstract

У роботі проводилося DFT-моделювання з метою встановлення умов, необхідних для протонування молекул аміаку на поверхні кремнію. Моделювання показало, що для того, щоб протонування було енергетично вигідним, необхідно виконання трьох умов: закріплення молекул NH3 хоча б на двох ОНгрупах, участь одної молекули води та наявність віддаленого атома бору. В цьому випадку енергія протонування (різниця між енергіями модельного кластера в непроторованому і протонованому станах) Eprot = 0,01еВ. Протонування відбувається при одночасному переході двох протонів: від поверхневої силанової групи до молекули Н2О та від Н2О до NH3. Такий перехід веде до утворення позитивно зарядженого іона NH4+ та негативно зарядженого pb-центра. Після цього відбувається перехід електрона від pb-центра до віддаленого атому бору, тобто його пасивація. Саме пасивація акцепторної домішки іонами NH4+ може призводити до зниження концентрації вільних дірок та провідності кремнієвих структур р-типу. Запропонована модель протонування молекул NH3 також може пояснити процеси поверхнево стимульованої лазерній десорбції-іонізації (SALDI). Вільні дірки, утворені під дією лазерного випромінювання, можуть рекомбінувати з електроном, локалізованим на атомі бору, тобто знімати його негативний заряд. Після цього енергія десорбції NH4+ знижується від 3,75 до 1,05 еВ, що є близьким до експериментально отриманих значень енергій лазерної десорбції іонів аміносполук. Моделювання також показало, що відсутність віддаленого атома бору суттєво знижує ефективність протонування, Eprot = – 0,50 еВ, оскільки негативний заряд зосереджується не на віддаленому атомі бору, а на pb-центрі, розміщеному поруч з NH4+. Значне кулонівське притягання між іоном NH4+ та зарядженим pb-центром веде до суттєвого зміщення атома кремнію, що є енергетично невигідним. Закріплення молекул NH3 лише на одній ОН-групі та відсутність проміжної молекули води також роблять протонування малоймовірним: значення Eprot в цих випадках становлять – 0,37 і – 0,08 еВ, відповідно. Це можна пояснити тим, що молекули Н2О і ОН-групи створюють енергетично вигідні водневі зв’язки з іоном NH4+ та екранують його електричне поле.

In this work, DFT modeling was carried out in order to establish the conditions necessary for the protonation of ammonia molecules on the silicon surface. Simulations have shown that for an energetically favorable protonation, three conditions are necessary: the fixing of NH3 molecules on at least two OH-groups, the participation of one water molecule, and the presence of a distant boron atom. In this case, the protonation energy (decrease in the cluster energy after protonation) is Eprot = 0.01 eV. Protonation occurs with the simultaneous transition of two protons: one from the surface silane group to the H2O molecule and the other from H2O to NH3. Such a transition leads to the formation of a positively charged NH4+ ion and a negatively charged pb-center (Si atom with a dangling bond). After that, the electron is transferred from the pb-center to the distant boron atom, and passivates it. Such passivation of acceptors by NH4+ ions can lead to a decrease in the concentration of free holes (and a decrease in conductivity) in p-type silicon. The proposed model of NH3 molecules protonation can also explain the processes of surface-assisted laser desorption/ionization (SALDI) of amino compounds. Free holes formed by laser radiation can recombine with electrons localized on passivated boron atoms, i.e., remove their negative charge. After this, the desorption energy of NH4+ decreases from 3.67 to 1.05 eV, which is close to the experimental values for amino ions. The simulation also showed that the absence of a distant boron atom significantly reduces the protonation efficiency, Eprot = – 0.50 eV, since the negative charge is localized not at the distant boron atom, but at the pb-center located near NH4+. Significant Coulomb attraction between the NH4+ ion and the charged pbcenter leads to a substantial displacement of the silicon atom, which is energetically unfavorable. The fixation of NH3 molecules on only one OH-group or the absence of water molecules also reduces the protonation probability: the value of Eprot in these cases is – 0.37 and – 0.08 eV, respectively. This can be explained by the fact that H2O molecules and OH-groups create energetically favorable hydrogen bonds with the NH4+ ion and shield its electric field.