Download this article
Download this article For screen
For printing
Recent Issues

Volume 19
Issue 5, 747–835
Issue 4, 541–746
Issue 3, 303–540
Issue 2, 157–302
Issue 1, 1–156

Volume 18, 5 issues

Volume 17, 5 issues

Volume 16, 5 issues

Volume 15, 5 issues

Volume 14, 5 issues

Volume 13, 5 issues

Volume 12, 5 issues

Volume 11, 5 issues

Volume 10, 5 issues

Volume 9, 5 issues

Volume 8, 8 issues

Volume 7, 10 issues

Volume 6, 9 issues

Volume 5, 6 issues

Volume 4, 10 issues

Volume 3, 10 issues

Volume 2, 10 issues

Volume 1, 8 issues

The Journal
About the journal
Ethics and policies
Peer-review process
 
Submission guidelines
Submission form
Editorial board
 
Subscriptions
 
ISSN 1559-3959 (online)
ISSN 1559-3959 (print)
 
Author index
To appear
 
Other MSP journals
Phase field simulation of the inclusion instability and splitting processes in interconnects due to interface diffusion induced by electromigration

Linyong Zhou, Peizhen Huang and Jiaming Zhang

Vol. 18 (2023), No. 1, 39–58
Abstract

Based on the bulk free energy density and the degenerate mobility constructed by the quartic double-well potential function, a phase field model is established to simulate the evolution of inclusions in interconnects due to interface diffusion in an electric field. The corresponding phase field governing equations are derived and the reliability of the program is proved by the agreement between the numerical simulation results and the theoretical analysis of the inclusion. The evolution of elliptical inclusions under different electric fields χ, different aspect ratios β, and different electric conductivity ratios λ is calculated using the mesh adaptation finite element method. The results show that the drift velocity of the circular inclusion is proportional to the electric field and inversely proportional to the electric conductivity ratio. There exist critical values of the electric field χc , the aspect ratio βc and the electric conductivity ratio λc . When λ λc , χ χc or β βc , the elliptical inclusions will split into several small inclusions. When λ > λc , χ < χc or β < βc , the elliptical inclusions will drift along the direction of the electric field as a relatively stable shape. The smaller the electric conductivity ratio of the inclusions λ, the greater the electric field χ or the aspect ratio β, the easier it is for the elliptical inclusions to split. Moreover, the time required for splitting increases with increasing the electric conductivity ratio or the aspect ratio, and decreases with increasing the electric field. In addition, the interconnect line with two inclusions is more complex, and the inclusions will split into more small ones under the high electric field strength, and the phenomenon of multiple merging and multiple splitting will occur.

Keywords
phase field method, electromigration, interface diffusion, finite element method, inclusion evolution
Milestones
Received: 20 November 2021
Revised: 5 June 2022
Accepted: 10 August 2022
Published: 24 March 2023
Authors
Linyong Zhou
The State Key Laboratory of Mechanics and Control of Mechanical Structures
Nanjing University of Aeronautics and Astronautics
Nanjing
China
Peizhen Huang
The State Key Laboratory of Mechanics and Control of Mechanical Structures
Nanjing University of Aeronautics and Astronautics
Nanjing
China
Jiaming Zhang
The State Key Laboratory of Mechanics and Control of Mechanical Structures
Nanjing University of Aeronautics and Astronautics
Nanjing
China