Circular Image

Z. Li

info

Please Note

9 records found

A combined nanoindentation test and hybrid Potts-phase field simulation study

Journal article (2026) - Junwei Chen, Zezhan Li, Bin Yang, Xiao Hu, Wenyu Li, Zichuan Li, Xuyang Yan, Guoqi Zhang, Jiajie Fan, More authors...
The reliability of through-glass via (TGV) interconnects is critical for advanced semiconductor packaging. This work investigates microstructural and mechanical evolution in electroplated TGV–Cu subjected to long-term aging at 250 °C. TGV samples were fabricated via laser-induced etching and double-sided copper electroplating, then aged for up to 1008 h. Nanoindentation revealed region-dependent reductions in hardness (from 2.0–2.5 GPa to below 0.5 GPa) and modulus (from 110–130 GPa to 40–90 GPa), with surface-near regions most affected. The glass substrate maintained stable mechanical properties until microcracks formed after 1008 h. EBSD quantification showed grain-size enlargement from 0.46 µm to 1.86 µm and a concurrent decrease in dislocation density. Molecular dynamics simulations of 3, 4, 5 nm grains corroborated the inverse relationship between grain size and micro-mechanical properties. A hybrid Potts-phase field model further linked grain coarsening to stress relaxation and elastic-energy minimization, revealing that as grains grow, the overall von Mises stress in the structure decreases; high-modulus grains retain relatively higher local stresses, while low-modulus, low-stress grains exhibit faster growth rates. Electrical I–V measurements confirmed stable ohmic behavior, despite a drop in insulation resistance. These integrated experimental and computational insights provide theoretical guidance for optimizing TGV interposer design and ensuring long-term operational reliability in heterogeneous integration technologies. (Figure presented.) ...
Journal article (2026) - Jiajie Fan, Xuyang Yan, Chuantong Chen, Leiming Du, Zichuan Li, Junran Zhang, Willem Van Driel, Guoqi Zhang
Sintered Cu nanoparticles (Cu NPs) are promising interconnection materials for high-temperature power electronics, yet how their authentic three-dimensional pore architecture governs microscale deformation remains unclear. Here, synchrotron nano-computed tomography (nano-CT) was combined with in-situ micropillar compression, explicit dynamic elastoplastic finite element analysis, and TEM/TKD characterization to interrogate sintered Cu NPs. The nano-CT voxel size was 45 nm, and the reconstructed volume corresponded to a cylinder 16 µm in diameter and 10 µm in height. The average sectional porosity was 12.44%, with a systematic discrepancy between two-dimensional and three-dimensional porosity quantification. During loading, the porosity decreased to 9.55% while the pore aspect ratio increased from 1.82–2.35. Finite element analysis further showed pronounced pore-adjacent stress/strain localization at the elastic–plastic transition, with local stress and equivalent plastic strain reaching 650 MPa and 1.7 × 10−2, compared with 250 MPa and 1.1 × 10−3 in adjacent regions. The GND density increased by 95.9% at a compressive strain of 26%, linking pore-induced strain gradients to dislocation accumulation. These results quantitatively connect authentic three-dimensional pore architecture, local deformation localization, and dislocation-mediated strengthening in sintered Cu NPs. Highlights Synchrotron nano-CT (45 nm voxel size) reconstructed a 16 × 10 µm cylindrical volume of sintered Cu NPs and resolved the authentic 3D pore network. Sectional porosity was 12.44%, and 2D/3D quantification showed a systematic discrepancy, with porosity decreasing to 9.55% and pore aspect ratio increasing from 1.82 to 2.35 during compression. Pore-adjacent localization was quantified at the elastic–plastic transition, with local stress/PEEQ reaching 650 MPa and 1.7 × 10−2 versus 250 MPa and 1.1 × 10−3 in adjacent regions. A 95.9% increase in GND density at 26% compressive strain links pore-induced strain gradients to dislocation accumulation and strain-gradient-driven strengthening. ...
Journal article (2026) - Zichuan Li, Jiajie Fan, Guoqi Zhang
This study is motivated by a conceptual inconsistency in the physical interpretation of eight-chain hyperelastic theory, which arises from the combined effect of two distinct issues: the use of the marginal projection distribution pz(|rz|) as a surrogate for the full probability density of end-to-end distance p(r̄), and the subsequent reliance on a root mean square (RMS) approximation step in the micro–macro averaging of chain stretch. We first revisit this probabilistic mismatch by reformulating the probability density function of freely-jointed chains (FJCs) in terms of the squared end-to-end vector r2, thereby restoring consistency on chain-level statistics. Building on this formulation, the micro–macro mapping averaging of chain conformational free energy is constructed directly in terms of r2, leading to a one-step mean-field approximation that avoids RMS averaging. The modified probability transformation is examined by Monte Carlo sampling at the microscopic level. To account for interchain interactions, q-mean statistical description of micro tube confinement was incorporated, leading to the appearance of the general invariant Iq1q2q3q. The resulting continuum constitutive model is assessed against multiaxial experimental data for several polymer networks, including vulcanized natural rubber, Entec Enflex S4035A thermoplastic elastomer, Tetra-PEG, and isoprene rubber vulcanizate. Comparisons with three existing hyperelastic strain energy formulations, the extended eight-chain, extended tube models, and the four-parameter ”comprehensive” model, demonstrate comparable phenomenological accuracy of the current model while providing a clearer and more consistent micro–macro physical interpretation of model parameters. A parametric study further illustrates how the dimensionless parameters n and q govern the shape of the macroscopic stress–strain responses. The present formulation provides a consistent theoretical basis within the scope of hyperelasticity and admits potential extensions toward more complex irreversible phenomena. ...
Journal article (2024) - Zichuan Li, Luutzen Wymenga, Jiajie Fan, Guoqi Zhang
This study demonstrates a breakdown analysis of the dynamics of a liquid crystal elastomer (LCE) including quality check, geometric measurement, thermal characterization, and comparison of heat- and light-induced contractions. A blue light-responsive acrylate side chain LCE with 1% azobenzene dye was characterized. From a classical viewpoint, photo-thermal contraction is considered a dominating effect, while direct photo-mechanical deformation can be neglected due to a low dye percentage. However, the findings of this research suggest that a low percentage of azobenzene dye does not necessarily lead to heat-dominating dynamics of LCE. This phenomenon has not yet been quantitatively studied before. The approach reported in this Letter can potentially be used to extract the data to improve the dynamics models of light-driven LCEs. ...
Conference paper (2024) - Xu Liu, Dong Hu, Zichuan Li, Xuejun Fan, Guoqi Zhang, Jiajie Fan
The nano-copper particles are widely used in the sintering processes of packaging wide bandgap semiconductors. Despite the significant success in the industry, the mechanism bridging the sintering process to the mechanical properties of sintered nano-copper is not yet well-modeled. In this paper, the impacts of different sintering temperatures and initial porosities caused by different stacking patterns on the uniaxial tensile performance of the sintered layer were studied via a molecular dynamics approach. Two stacking patterns, simple cubic and face-centered cubic, were simulated, respectively. Evolution of their structure at temperatures of 300, 400, 500, and 600 K were simulated as the sintering process. Afterward, the sintered structures were subjected to uniaxial tensile with rates of 0.01 and 0.04 Å/ps at different temperatures to compare the mechanical properties. The results show that the sintering rate and density of the sintered structure increase with a higher temperature. However, the tensile strength of the sintered structure is less relevant to the difference in stacking pattern. This study proves that porosity has a greater effect on sintering quality. ...
Conference paper (2024) - Dong Hu, Chieh Wang, Zichuan Li, Nikhil Gupta, René H. Poelma, Ziliang Shi, Jiajie Fan, Guoqi Zhang
Driven by the increasing demand for high-power systems, ceramic substrates have received more attention for handling higher power density. Warpage in active metal brazed (AMB) ceramic substrate becomes a critical issue as it can deteriorate the reliability performance. This study comprises three phases, including investigation of the cause of the warpage, validation of the proposed model, and optimization for effective warpage management. At first, the coefficient of thermal expansion (CTE) and yield strength of the copper (Cu) layer in AMB were characterized and adopted in a two-dimensional (2D) finite element model. The evolution of simulated strain and moments revealed the cause of the warpage during the manufacturing processes. Furthermore, the 2D model was extended to a three-dimensional (3D) model. The finite element method (FEM) and experiments were conducted on different heat treatment conditions for 3D model validation. The validated 3D model was applied to carry out a design of experiments (DoEs) for design optimization to reduce the warpage. Consequently, the factor analysis in DoEs was demonstrated by different pattern designs using subtractive milling techniques. ...
Conference paper (2023) - Dong Hu, Zichuan Li, Jiajie Fan, Guoqi Zhang
As a critical part of speeding up industrial electrification, power electronics, and its packaging technology are undergoing rapid development. Cu nanoparticle sintering technology has therefore received extensive attention for its excellent performance in the die-attachment layer, where the mechanical properties are essential to be known for design for reliability. Both sintering and subsequent nanoindentation were studied by simulation. The effect of porosity on the nanoindentation response was investigated by setting different initial packing densities. In addition, the impact of indenter size and indentation positions on the nanoindentation response were discussed. The nanoindentation behaviors were studied by extracting loading-displacement (P-h) curves and calculating elastic modulus and hardness. In addition, the microstructure evolution was adopted using atomic configuration to study the nanoindentation mechanism. This work offers valuable insights into the Cu sinter paste preparation phase for sintering technology. ...
Journal article (2022) - Jue Wang, Weiyi Yuan, Zichuan Li, Yuri Trofimov, Sergey Lishik, Jiajie Fan
Some atrial contractile assist devices applied on the heart surface can be regarded as a laminated Liquid crystal elastomer (LCE) plate under steady temperature loads and a contact mechanical force. An exact solution for the deformation of the laminated LCE plate under combined thermal and mechanical loads is derived by solving the three-dimensional (3D) equilibrium equations including heat conduction and thermoelastic theory. The validity of mathematical formula and computer programming is proved by convergence and comparison examples with finite element method (FEM). In order to simplify the complex calculation of exact solution, a back propagation neural network (BPNN) is further trained with a database containing 9504 sets of thermo-mechanical load conditions and their corresponding deformation which is solved by the exact solutions. Then the deformations of LCE plate subject to combined thermo-mechanical load can be predicted by this BP neural network instead of complex numerical calculation. Moreover, it is also applied to inverse the contact mechanical force at the bottom surface of LCE plate with a given deformation and temperature conditions. The results show that: (1) The results from the exact theoretical solution are in consistence with that from FEM but have a higher computational efficiency and stability; (2) The deformation of the laminated plate is more sensitive to the layered thickness of LCE than the variation of the temperature; (3) 3-D elasticity solutions of a laminated LCE plate under the combined thermos-mechanical load can be effectively predicted by a trained BP neural network. ...
Journal article (2022) - Jue Wang, Weiyi Yuan, Zichuan Li, Yingcan Zhu, Thebano Santos, Jiajie Fan
This paper presents a convenient and efficient method to predict the mechanical solutions of a laminated Liquid Crystal Elastomers (LCEs) system subjected to combined thermo-mechanical load, based on a back propagation (BP) neural network which is trained by machine learning from a database established by analytical solutions. Firstly, the general solutions of temperature, displacement, and stress of any single layer in the LCEs system are obtained by solving the two-dimensional (2D) governing equations of both heat conduction and thermoelasticity. Then, the unknown coefficients in above general solutions are determined by a transfer-matrix method based on the continuity condition at the interface of adjacent layers and the combined thermo-mechanical loads condition at the surface of the LCEs system. The formula derivation and calculator program are verified through convergence studies and comparisons with FEM results. Finally, a database with displacements of LCEs system in a temperature field subjected to 561 sets of mechanical loads is established based on the presented analytical model. The BP neural network based on above database is further applied to establish the relationship between deformation and mechanical load to predict the elastic deformation of the LCEs system in a temperature field subjected to a mechanical load. Moreover, the BP network can also inverse the coefficients of mechanical load which induces the specific deformation in a temperature field. The numerical examples show that: (1) The deformation of a laminated LCEs system due to thermal load is limited within the range of human temperature changes from 36 °C to 40 °C. (2) The thickness of the LCE is a sensitive parameter on the deformation at the bottom surface of the system. (3) The accuracy of predicted displacements induced by the thermo-mechanical load and the inversed mechanical load based on deformation of the LCEs system in a temperature field using BP neural network reaches 99.6% and 98.5% respectively. ...