For coaxial-annular through silicon vias (CA-TSV) structure with superior performances, characteristic impedance, power, time constant and analytical models of parasitic parameters are proposed and effects of structural parameters on electrical properties are studied. S21 parameter was verified by software HFSS. It shows that increasing inner diameter of CA-TSV or reducing outer diameter reduces characteristic impedance, while reducing inner diameter of CA-TSV or increasing outer diameter reduces its power consumption effectively. Increasing inner diameter of CA-TSV or outer diameter reduces time constant of RC equivalent circuit, whereas increasing inner diameter of CA-TSV or reducing outer diameter reduces time constant of RL equivalent circuit. Increasing inner diameter of CA-TSV or outer diameter reduces resistance effectively and capacitance can be increased significantly. It provides reference for electrical properties of three-dimensional integrated circuits based on TSV interconnects.

An expression of thermal resistance matrix is given.Transient temperature characteristics of three-dimensional chip-muhiprocessors(3D-CMP) are studied.Effect of heat capacity,thermal resistance and power consumption on temperature is analyzed. It shows that steady temperature of 3D-CMP is limited effectively by reducing thermal resistance and power consumption.Heat capacity influences rise time of temperature.It does not affect steady temperature.

With through silicon via(TSV) area scale factor r,an analytical thermal model for top layer of three-dimensional integrated circuits(3D IC) taking TSV into account was proposed.It is shown that temperature is lower after considering TSVs under same working conditions;the greater the scale factor r,the lower the temperature is;For more layers and smaller r,temperature increases sharply with decrease of r;The best range of TSV area ratio factor r is 0.5% to 1% for an 8-layer 3D IC.

With consideration of via effect and heat fringing effect, a thermoelectric simulation method is proposed which modifies node heat flow due to temperature distribution. Based on thermoelectric duality, thermal resistance models including inner/inter-layer and vias are presented. Take advantage of feedback relationship between heat and electric, the node network heat flow model is modified with temperature distribution. Multilevel interconnects temperature distribution with polymer and silicon oxide as insulator dielectric are analyzed. Compared with results of finite element, the relative standards deviation of the proposed method can be reduced by 71.2% and 12. 9% respectively than those of available models. With consideration of via effect and heat fringing effect, we calculate peak temperature rise in different technology nodes. It shows that interconnect temperature distribution is overestimated in traditional models.

Binding energy of substrate,coupling energy of nearest neighbor particles and strain field are introduced to diffusion barrier.Cluster growth process at different temperatures is investigated with Monte Carlo simulation.It shows that as 400 K≤T ≤ 480 K the average width of ramified clusters is independent on temperature T and is almost equal to the diameter of single particle.As 500 K≤T ≤ 680 K,however,width of the ramified cluster increases gradually to that of 4 particles with increasing temperature.As T increases further,clusters with large number of particles disappear,due to strong activity of particles.Average number of particles in each cluster is smaller than 2.Evolution of cluster morphology during growth process and number of cluster at temperatures are studied as well.

Based on a brief description of Corrected Smoothed Particle Hydrodynamics (CSPH) method, the discretization scheme of conservation equations in axisymmetric coordinate systems is developed and a method to determine the number of particles in the computational domain is proposed. Numerical examples for hypervelocity impact are carried out and the results clearly show that CSPH is not only a suitable method for the simulation of large deformation problems in impact dynamics, but also a new method with the advantages that can not be replaced by finite element/ finite difference method.