主权项 |
一种铝栅CMP协同计算模型的仿真方法,其特征在于,包括步骤:获取铝栅初始表面高度;依据所述铝栅初始表面高度,采用铝栅CMP压力分布精确计算模型获取研磨垫和铝栅间的压力分布p(x,y,t);向铝栅CMP协同计算模型中输入模型参数,其中,所述铝栅CMP协同计算模型由所述铝栅CMP压力分布精确计算模型和铝栅CMP化学反应动力学模型确定;所述模型参数至少包括研磨垫材料属性参数、CMP工艺参数以及化学动力学参数,所述CMP工艺参数至少包括外部压力p<sub>0</sub>;采用所述铝栅CMP协同计算模型并依据所述模型参数进行预定时间段的仿真,预测所述铝栅研磨后的表面高度;其中,所述铝栅CMP协同计算模型为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>MRR</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>M</mi><msub><mi>ρ</mi><mn>0</mn></msub></mfrac><mfrac><mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>60</mn></msub><msub><mi>k</mi><mn>2</mn></msub><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>/</mo><msub><mi>p</mi><mn>0</mn></msub><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><msub><mi>k</mi><mn>3</mn></msub><mo>[</mo><mi>Oxi</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>4</mn></msub><mo>[</mo><mi>CA</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>50</mn></msub><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>/</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>k</mi><mn>2</mn></msub><msub><mi>k</mi><mn>3</mn></msub><mo>[</mo><mi>Oxi</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>2</mn></msub><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mo>[</mo><mi>In</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>4</mn></msub><mo>[</mo><mi>CA</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>50</mn></msub><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>/</mo><msub><mi>p</mi><mn>0</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000679267490000011.GIF" wi="1602" he="185" /></maths>所述铝栅CMP压力分布精确计算模型为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><mi>w</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mn>1</mn><mo>-</mo><msup><mi>υ</mi><mn>2</mn></msup></mrow><mi>πE</mi></mfrac><msub><mo>∫</mo><mi>I</mi></msub><mi>p</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mfrac><mn>1</mn><msqrt><msup><mrow><mo>(</mo><mi>x</mi><mo>-</mo><mi>ξ</mi><mo>)</mo></mrow><mn>2</mn></msup><mo>+</mo><msup><mrow><mo>(</mo><mi>y</mi><mo>-</mo><mi>η</mi><mo>)</mo></mrow><mn>2</mn></msup></msqrt></mfrac><mi>dξdη</mi></mrow>]]></math><img file="FDA0000679267490000012.GIF" wi="1134" he="150" /></maths>g(x,y,t)=h(x,y,t)+w(x,y,t)‑c≥0,(x,y,t)∈I<sub>C</sub>g(x,y,t)=0,p(x,y,t)>0,(x,y,t)∈I<sub>C</sub>g(x,y,t)>0,p(x,y,t)=0,<img file="FDA0000679267490000014.GIF" wi="266" he="73" />F<sub>0</sub>(t)=∫<sub>I</sub>p(x,y,t)dxdy;所述铝栅CMP化学反应动力学模型为:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>MRR</mi><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>=</mo><mfrac><mi>M</mi><msub><mi>ρ</mi><mn>0</mn></msub></mfrac><mfrac><mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>6</mn></msub><msub><mi>k</mi><mn>2</mn></msub><mo>+</mo><msub><mi>k</mi><mn>2</mn></msub><msub><mi>k</mi><mn>3</mn></msub><mo>[</mo><mi>Oxi</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>4</mn></msub><mo>[</mo><mi>CA</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><mo>)</mo></mrow></mrow><mrow><msub><mi>k</mi><mn>2</mn></msub><msub><mi>k</mi><mn>3</mn></msub><mo>[</mo><mi>Oxi</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><mrow><mo>(</mo><msub><mi>k</mi><mn>2</mn></msub><mo>+</mo><msub><mi>k</mi><mn>1</mn></msub><mo>[</mo><mi>In</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>)</mo></mrow><mrow><mo>(</mo><msub><mi>k</mi><mn>4</mn></msub><mo>[</mo><mi>CA</mi><mo>]</mo><mrow><mo>(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo>,</mo><mi>t</mi><mo>)</mo></mrow><mo>+</mo><msub><mi>k</mi><mn>5</mn></msub><mo>)</mo></mrow></mrow></mfrac><mo>;</mo></mrow>]]></math><img file="FDA0000679267490000013.GIF" wi="1631" he="166" /></maths>其中,MRR(x,y,t)为铝栅研磨去除率,M为铝的原子质量,ρ<sub>0</sub>为铝的密度,[Oxi](x,y,t)为研磨液中氧化剂的浓度,[In](x,y,t)为研磨液中抑制剂的浓度,[CA](x,y,t)为研磨液中螯合剂的浓度,k<sub>i</sub>为化学反应速率常数,i=1,2,…,6,k<sub>50</sub>为研磨粒子去除氧化膜、螯合物的反应速率常数,k<sub>60</sub>为研磨粒子去除铝单质的反应速率常数,p<sub>0</sub>为外部压力,x为选定坐标系沿x轴方向的坐标值,y为选定坐标系沿y轴方向的坐标值,t为铝栅CMP的仿真时间,w(x,y,t)是研磨垫的表面形变,υ是研磨垫的泊松比,E是研磨垫的弹性模量,g(x,y,t)是形变后研磨垫和铝栅表面之间的间隔,h(x,y,t)是研磨垫和铝栅表面之间的初始间隔,c是研磨垫沿外加载荷方向的整体位移,I<sub>C</sub>是研磨垫和铝栅表面的接触表面积,F<sub>0</sub>(t)是t时刻的外加载荷,I是整个铝栅表面积,ξ、η均为积分变量。 |