Proton and gamma radiation simulation modeling method and system for tsv-based microsystem interconnect modules

Abstract

The present application relates to the technical field of three-dimensional integrated microsystem, and particularly relates to a proton and gamma radiation simulation modeling method and system of a microsystem interconnection module based on TSV. The method comprises: obtaining geometric parameters and material characteristics of a through silicon via, a redistribution layer and a micro bump, establishing a three-dimensional finite element model and analyzing its scattering characteristics, and constructing a correlation model of material characteristics and scattering parameters; extracting circuit parameters in an ADS compilation environment, establishing an RLC equivalent circuit and its correlation with scattering parameters; through a proton and gamma irradiation experiment, obtaining experimental data of scattering parameters under different doses, and then establishing a correlation model between radiation dose and scattering parameters and RLC parameters, to realize accurate prediction of the electrical performance of a three-dimensional integrated microsystem module under different radiation conditions, and to provide a scientific basis for reliability design of a three-dimensional integrated microsystem module in a radiation environment such as aerospace.

Description

Technical Field

[0001] This invention relates to the field of three-dimensional integrated microsystems technology, and in particular to a method and system for simulating and modeling proton and gamma radiation of microsystem interconnect modules based on TSV. Background Technology

[0002] 3D integration technology, as a key means to improve circuit integration and performance, significantly reduces signal transmission delay and optimizes system power consumption and package size by stacking and interconnecting multiple chips. The introduction of through-silicon via (TSV) technology, enabling vertical chip interconnection, is the core of 3D integration technology. However, in specialized applications such as aerospace, 3D integrated microsystems face the challenge of harsh radiation environments.

[0003] High-energy particles in the space radiation environment, such as protons and gamma rays, can significantly impact the signal integrity of three-dimensional microsystems. Proton irradiation leads to lattice defects and radiation-induced traps, increasing carrier recombination and trapping, thus affecting the switching characteristics and signal transmission speed of devices. Gamma irradiation induces electron-hole pairs in materials, forming radiation traps and defects, altering the impedance and capacitance characteristics of circuits in microsystems, and increasing signal attenuation. The combined effect of proton and gamma irradiation exacerbates radiation damage in materials, leading to more severe electrical performance degradation and affecting the long-term stability of the system.

[0004] Currently, simulation studies of three-dimensional interconnect structures mainly focus on structural design and electrical performance optimization. Finite element analysis software (such as HFSS) and circuit simulation tools are widely used for the electrical performance analysis of microsystems. However, in reliability studies under radiation environments, existing technologies still have the following shortcomings:

[0005] 1. There is a lack of accurate radiation damage modeling methods for three-dimensional interconnect structures (especially TSV structures). Existing models struggle to accurately describe radiation-induced changes in material parameters and their impact on signal transmission characteristics.

[0006] 2. The simulation results were not adequately considered due to defects in the actual sample manufacturing process. Defects of varying degrees that occur during the manufacturing process of actual samples lead to subtle changes in the structure of the test samples. The randomness of the defect distribution makes accurate modeling impossible in the software, resulting in discrepancies between the simulation results and actual measurements.

[0007] 3. The equivalent circuit model is oversimplified and cannot accurately reflect the electrical characteristics of complex three-dimensional interconnect structures. Existing equivalent circuit models may have ignored some key electrical components during the simplification process, leading to discrepancies between simulation results and measured data.

[0008] Therefore, it is necessary to further develop accurate radiation damage modeling methods, consider the impact of manufacturing defects, and optimize equivalent circuit models, so as to provide more accurate theoretical and technical support for the reliability study of three-dimensional integrated microsystems under radiation environment. Summary of the Invention

[0009] The purpose of this invention is to overcome the problem of low reliability of three-dimensional integrated microsystem modules under proton irradiation, thereby providing a proton irradiation simulation modeling method and system for microsystem interconnection modules based on TSV.

[0010] To solve the above-mentioned technical problems, the present invention provides a method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV, comprising:

[0011] Step 1: Obtain the geometric parameters and material properties of the through-silicon vias, redistribution layers, and microbumps of the 3D microsystem;

[0012] Step 2: Based on geometric parameters and material properties, establish a finite element model of a three-dimensional microsystem including through-silicon vias, redistribution layers, and microbumps using finite element analysis software;

[0013] Step 3: By adjusting the material parameters of the finite element model of the three-dimensional microsystem, the scattering parameters under different material property conditions are obtained, and a correlation model between the material property parameters and the scattering parameters is constructed.

[0014] Step 4: In the ADS compilation environment, extract the circuit parameters of through silicon vias, redistribution layers, and microbumps from the finite element model, and construct the RLC equivalent circuit; based on the RLC equivalent circuit, establish a correlation model between the RLC equivalent circuit and the scattering parameters.

[0015] Step 5: Conduct proton irradiation experiments with a preset radiation flux, followed by gamma radiation experiments with different radiation doses to test scattering parameters and obtain experimental data on scattering parameters; based on the experimental data on scattering parameters, construct a correlation model between gamma radiation dose and scattering parameters.

[0016] Step 6: Construct a correlation model between gamma radiation dose and RLC equivalent circuit, and based on the correlation model between gamma radiation dose and RLC equivalent circuit, construct correlation models between gamma radiation dose and resistance, inductance and capacitance in RLC equivalent circuit respectively.

[0017] As an improvement to the above method, in step 1, the geometric parameters of the through-silicon via include: radius, length, and insulating layer parameters; the geometric parameters of the redistribution layer include: width and height; the geometric parameters of the microbumps include: size; and the material properties include: electrical conductivity, magnetic permeability, and dielectric constant.

[0018] As an improvement to the above method, in step 2, the finite element analysis software used is HFSS simulation software, specifically including: using HFSS simulation software, based on geometric parameters and material properties, preset boundary conditions and port excitation, to establish a finite element model of a three-dimensional microsystem containing through silicon vias, redistribution layers and microbumps, wherein the boundary conditions include the electrical conductivity, dielectric constant and magnetic permeability of the material; the frequency range of the port excitation is set to 0.1-26.5 GHz.

[0019] As an improvement to the above method, step 3 specifically includes: adjusting the material property parameters of the through-silicon vias, redistribution layers, and microbumps in HFSS simulation software; extracting scattering parameters under different material property conditions; and constructing a correlation model between the material property parameters and scattering parameters at a fixed frequency point using a polynomial fitting method; wherein,

[0020] The model for conductivity S1 is:

[0021] S1=aμ 2 +bμ+c

[0022] The model for the dielectric constant S2 is:

[0023] S2=aε 2 +bε+c

[0024] Where a is the first coefficient, b is the second coefficient, c is the constant term, μ is the relative permeability of the material, and ε is the relative permittivity of the material.

[0025] As an improvement to the above method, step 4, which involves extracting the equivalent circuit parameters of the redistribution layer, includes:

[0026] Extract the total resistance R of the redistribution layer. RDL :

[0027]

[0028] Among them, R RDL,dc R is the DC resistance of the redistribution layer. RDL,ac For the AC resistance of the redistribution layer; where,

[0029]

[0030] Where, ρ RDL For the material resistivity of the redistribution layer, l RDL w is the length of the redistribution layer. RDL t is the conductor width of the rewiring layer. RDL σ represents the conductor thickness of the rewiring layer. RDL The conductivity of the material used in the redistribution layer, where,

[0031]

[0032] Where μ0 is the free permeability, μ r,RDL S represents the relative permeability of the redistribution layer material. RDL This represents the distance between adjacent redistribution layers.

[0033] Extract the capacitor C of the insulating layer below the redistribution layer. RDLtoSub :

[0034]

[0035] Where ε0 is the vacuum permittivity, ε RDLox1 is the dielectric constant of the insulating layer below the redistribution layer, and lox is the length of the insulating layer;

[0036] Extracting the redistribution layer in silicon equivalent capacitance C RDLinSub :

[0037]

[0038] ε r,eff h is the effective dielectric constant. eff For effective height, among which,

[0039]

[0040] Where h is the measurement height of the redistribution layer.

[0041] As an improvement to the above method, step 4, which involves extracting the equivalent circuit parameters of the through-silicon vias and microbumps, includes:

[0042] Extracting the DC resistance R of a through-silicon via TSV,dc :

[0043]

[0044] Extracting the AC resistance R of a through-silicon via TSV,ac :

[0045]

[0046] Where, ρ TSV h is the resistivity of the material in a through-silicon via. TSV r is the height of the through-silicon via. TSV Let δ be the radius of the through-silicon via. skindepth,TSV Indicates the depth of current in a through-silicon via;

[0047] Extracting the DC resistance R of the micro-bump bump,dc :

[0048]

[0049] Extracting the AC resistance R of the micro-bump BUMPac :

[0050]

[0051] Where, ρ BUMP h is the resistivity of the material with micro-bumps. Bump r is the height of the micro-bump. Bump Let δ be the radius of the micro-bump. skindepth,BUMP This indicates the depth of the current in the microbump;

[0052] Extracting the inductance L of the micro-bump Bump :

[0053]

[0054] Among them, h TSV r is the height of the through-silicon via. TSV p is the radius of the through-silicon via. TSV The distance between two adjacent silicon via pillars is denoted by f, and f is the signal frequency when the circuit is in use.

[0055]

[0056] Among them, C TSV-ax ε represents the total oxide capacitance of the through-silicon via. TSVox Let r be the dielectric constant of the oxide layer material of the through-silicon via. dep t is the outer radius of the maximum depletion layer of the oxide layer. ox C represents the structural thickness of the oxide layer. ox For silicon dioxide dielectric capacitors with oxide layers, C dep This is the depletion region capacitance when the depletion layer thickness is at its maximum in a silicon substrate with an oxide layer.

[0057] As an improvement to the above method, in step 5, the preset radiation fluence is E12 P / cm. 2 Different radiation doses were included: 100 krad, 300 krad, 500 krad, and 900 krad; scattering parameters were tested using a network analyzer; using a polynomial fitting method, based on the experimental data of the scattering parameters, the relationship between gamma radiation dose and scattering parameters was established, linear and nonlinear influence coefficients were determined, and a correlation model between gamma radiation dose and scattering parameters was constructed.

[0058] S11=dD 3 +eD 2 +fD+g;

[0059] Where S11 is the reflection coefficient, D is the radiation dose, d is the cubic coefficient characterizing the higher-order nonlinear effect of the radiation dose on the reflection coefficient S11, e is the quadratic coefficient reflecting the nonlinear effect of the radiation dose on the reflection coefficient S11, f is the linear coefficient representing the degree of linear effect of the radiation dose on the reflection coefficient S11, and g represents the initial reflection coefficient value under no-irradiation conditions.

[0060] As an improvement to the above method, step 6 is followed by a verification step: taking the equivalent circuit parameters under different gamma radiation dose conditions as input, performing circuit simulation in the ADS simulation environment, and obtaining the simulation results of the scattering parameters; performing error analysis between the simulation results of the scattering parameters and the experimental data of the scattering parameters to verify the accuracy of the correlation model between gamma radiation dose and scattering parameters.

[0061] As an improvement to the above method, the error analysis includes the following steps:

[0062] The degree of matching (R) between the correlation model of gamma radiation dose and scattering parameters and the experimental data of scattering parameters was obtained. 2 :

[0063]

[0064] Where SSE represents the predicted value of scattering parameter S21 fitted by the correlation model between gamma radiation dose and scattering parameters. The measured value y of the scattering parameter S21 obtained from the experiment i The squared difference between them, SST represents the measured value of the scattering parameter S21 obtained in the experiment. i The measured mean of the scattering parameter S21 obtained from the experiment The difference of squares between them; where,

[0065]

[0066] To achieve another objective of the present invention, the present invention also provides a proton and gamma radiation simulation modeling system for a microsystem interconnect module based on TSV, comprising:

[0067] The finite element model building module is used to obtain the geometric parameters and material properties of through-silicon vias, redistribution layers, and microbumps in a three-dimensional microsystem. Based on the geometric parameters and material properties, a finite element model of the three-dimensional microsystem containing through-silicon vias, redistribution layers, and microbumps is built using finite element analysis software.

[0068] The module for building a correlation model between material property parameters and scattering parameters is used to adjust the material parameters of the finite element model of a three-dimensional microsystem, obtain the scattering parameters under different material property conditions, and build a correlation model between material property parameters and scattering parameters.

[0069] The module for constructing the correlation model between RLC equivalent circuit and scattering parameters is used in the ADS compilation environment to extract the circuit parameters of through silicon vias, redistribution layers and microbumps in the finite element model and construct the RLC equivalent circuit; based on the RLC equivalent circuit, a correlation model between the RLC equivalent circuit and scattering parameters is established.

[0070] A module for constructing a correlation model between gamma radiation dose and scattering parameters is used to conduct proton irradiation experiments with preset radiation fluence, and also to conduct gamma radiation experiments with different radiation doses to test scattering parameters and obtain experimental data on scattering parameters; based on the experimental data on scattering parameters, a correlation model between gamma radiation dose and scattering parameters is constructed; and,

[0071] The RLC equivalent circuit correlation model construction module is used to construct a correlation model between gamma radiation dose and RLC equivalent circuit, and based on the correlation model between gamma radiation dose and RLC equivalent circuit, construct correlation models between gamma radiation dose and resistance, inductance and capacitance in RLC equivalent circuit respectively.

[0072] Compared to existing technologies, the advantages of this invention lie in its method and system for proton irradiation simulation modeling of TSV-based microsystem interconnect modules, which solves the reliability problem of three-dimensional integrated microsystem modules under proton irradiation and gamma radiation. Addressing the complex effects of three-dimensional integrated microsystem modules in radiation environments, this invention establishes a correlation model between radiation dose, frequency, equivalent circuit parameters, and scattering parameters through finite element analysis and multivariate model construction, combined with univariate analysis calculations and model diagnostic verification. This overcomes the shortcomings of existing technologies in simulating and modeling the radiation effects of three-dimensional integrated microsystem modules, enabling accurate prediction of the electrical performance of three-dimensional integrated microsystem modules under different radiation conditions, and providing a scientific basis for the reliability design of three-dimensional integrated microsystem modules in radiation environments such as aerospace. Attached Figure Description

[0073] Figure 1 A flowchart illustrating the proton irradiation simulation modeling method for TSV-based microsystem interconnect modules;

[0074] Figure 2 A 3D view of TSV;

[0075] Figure 3 First schematic diagram showing the relationship between scattering parameter S21 and frequency as measured by experiments and simulations respectively;

[0076] Figure 4 The second schematic diagram shows the relationship between the scattering parameter S21 and the frequency as measured in experiments and simulations, respectively. Detailed Implementation

[0077] The technical solutions provided by the present invention will be further illustrated below with reference to the embodiments.

[0078] Example 1

[0079] like Figure 1 As shown, the proton irradiation simulation modeling method for a TSV-based microsystem interconnect module provided in this embodiment includes the following steps:

[0080] 1. Sample parameter acquisition

[0081] Geometric parameters and material properties of the 3D microsystem modules are extracted, including: the radius, length, filler material, and insulating layer parameters of the TSV; the width, height, and material properties of the redistribution layer; the size and material properties of the microbumps; and the conductivity and dielectric constant of the substrate material. High-precision parameter extraction is ensured using specialized measurement equipment.

[0082] 2. Geometric Model Construction

[0083] In the HFSS simulation environment, a three-dimensional microsystem finite element model containing structures such as TSV, RDL, and BUMP is established based on the extracted sample parameters. Appropriate boundary conditions and excitation sources are set to ensure that the model can realistically simulate electrical behavior and response under radiation conditions.

[0084] 3. Correlation model between material parameters and scattering parameters

[0085] In HFSS, scattering parameter curves under different material conditions are obtained by adjusting the material parameters in the finite element model. Key frequency points are selected, and the relationship between material parameters and scattering parameters is established using polynomial fitting.

[0086] 4. Construction of equivalent circuit model

[0087] In the ADS environment, based on the geometric model analysis results, the equivalent circuit parameters of RDL, TSV, and BUMP are extracted. An RLC equivalent circuit model is constructed, and the parasitic resistance and parasitic capacitance and their calculation methods are clarified.

[0088] 5. Correlation model between circuit characteristics and scattering parameters

[0089] Based on the equivalent circuit model extraction method, and combined with the correlation model between material parameters and scattering parameters, the relationship between S-parameters and RLC equivalent circuit parameters is established.

[0090] 6. Combined Proton and Gamma Irradiation Experiment

[0091] Total dose radiation experiments were conducted on the samples, with different dose points set, and the scattering parameter curves under the radiation environment were recorded. A network analyzer was used to test and ensure the accuracy of the data.

[0092] 7. Modeling the relationship between radiation dose and scattering parameters

[0093] Using polynomial fitting, a relationship between radiation dose and scattering parameters was established based on experimental data.

[0094] 8. Correlation Model between Radiation Dose and Equivalent Circuit Parameters

[0095] Based on the correlation model between scattering parameters and RLC parameters, a model of the influence of radiation dose on equivalent circuit parameters is further established.

[0096] 9. Model Validation and Optimization

[0097] The equivalent circuit parameters under different radiation dose conditions were used as input to perform circuit simulation in the ADS simulation environment. The simulated scattering parameters were obtained and compared with experimental data to perform error analysis and optimize the accuracy of the model.

[0098] 10. Comprehensive analysis of multiple frequencies and multiple doses

[0099] Comprehensive simulations and analyses were performed under multiple frequency and dose conditions in HFSS and ADS environments. The applicability of the model under different electrical conditions was verified, providing comprehensive data support for the reliability design of three-dimensional microsystem modules.

[0100] The following provides a more detailed description of each step.

[0101] (1) Obtaining sample parameters

[0102] Obtain the geometric parameters and material properties of the three-dimensional microsystem modules, and comprehensively extract the dimensions and material properties of key structures such as TSVs (Through Silicon Vias), RDLs (Rewiring Layers), and BUMPs (Microbumps), including the radius, length, filling material, and insulating layer of the TSVs, the width, height, and material properties of the rewiring layers, as well as the electrical conductivity and dielectric constant of the materials. Figure 2 A three-dimensional view of a through-silicon via (TSV) is shown.

[0103] (2) Constructing a geometric model

[0104] Using HFSS software, based on these material parameters and geometric parameters and other boundary conditions, a finite element model of a microsystem module containing structures such as silicon vias, microbumps, and redistribution was established.

[0105] Define boundary conditions, including the material's electrical conductivity, dielectric constant, and magnetic permeability.

[0106] Configure port excitation with a frequency range of 0.1-26.5GHz.

[0107] (3) Construct a correlation model between material parameters and scattering parameters

[0108] By adjusting the material parameters (copper conductivity, silicon dielectric constant, silicon dioxide dielectric constant) in HFSS, scattering parameters under different material conditions are extracted.

[0109] At a fixed frequency point, a polynomial fitting method is used to establish the relationship between scattering parameters and material parameters.

[0110] Conductivity: S1=aμ 2 +bμ+c

[0111] Dielectric constant: S² = aε 2 +bε+c

[0112] Where a is the first coefficient, b is the second coefficient, c is the constant term, μ is the relative permeability of the material, and ε is the relative permittivity of the material.

[0113] (4) Constructing an equivalent circuit model

[0114] Based on the ADS (ARM Developer Suite) compilation environment and geometric model, the circuit parameters of redistribution layers, through-silicon vias, and microbumps are extracted through equivalent analysis of module parameters. The equivalent methods for key parameters such as parasitic resistance, capacitance, and inductance are clarified, and an RLC equivalent circuit model is established.

[0115] Redeployment layer parameter extraction:

[0116]

[0117] Where R RDL R is the total resistance of the redistribution layer. RDL,dc R is the DC resistance of the redistribution layer. RDL,ac For the AC resistance of the redistribution layer, ρ RDL For the material resistivity of the redistribution layer, w RDL t is the conductor width of the rewiring layer. RDL σ represents the conductor thickness of the rewiring layer. RDL The conductivity of the material used for the redistribution layer.

[0118]

[0119] Where μ0 is the free permeability, μ r,RDL S represents the relative permeability of the redistribution layer material. RDL This represents the distance between adjacent rewiring layers.

[0120] Capacitor C of the insulating layer below the redistribution layer RDLtoSub And the equivalent capacitance C in silicon of the redistribution layer RDLinSub The calculation formula is as follows. Where ε RDLox1 t is the dielectric constant of the insulating layer below the redistribution layer. RDLoxε represents the height of the insulating layer below the redistribution layer. r,eff h is the effective dielectric constant. eff For effective height, l RDL This is the length of the redistribution layer.

[0121]

[0122] Extraction of parasitic parameters for through-silicon vias and microbumps:

[0123]

[0124]

[0125] Among them, R TSV,dc R is the DC resistance of the through-silicon via. TSV,ac R is the AC resistance of the through-silicon via. bump,dc ρ is the DC resistance of the micro-bump. TSV h is the resistivity of the through-silicon via material. TSV r is the height of the through-silicon via. TSV Where δ is the radius of the through-silicon via. skin-depth-TSV h represents the depth of current in a through-silicon via. Bump r is the height of the micro-bump. Bump Let δ be the radius of the micro-bump. skin-depth-Bump This indicates the depth of the current in the microbump.

[0126]

[0127] Among them, L Bump For inductors with micro-scratches, h TSV r is the height of the through-silicon via. TSV p is the radius of the through-silicon via. TSV This represents the spacing between two adjacent through-silicon via pillars.

[0128]

[0129] Among them, C TSV-ax ε represents the total oxide capacitance of the through-silicon via. TSVox Let r be the dielectric constant of the oxide layer material of the through-silicon via. dep t is the outer radius of the maximum depletion layer of the oxide layer. ox C represents the structural thickness of the oxide layer. ox For silicon dioxide dielectric capacitors with oxide layers, C dep This is the depletion region capacitance when the depletion layer thickness is at its maximum in a silicon substrate with an oxide layer.

[0130] (5) Establish a correlation model between circuit characteristics and scattering parameters

[0131] Based on the equivalent circuit extraction method and combined with the correlation model between material parameters and scattering parameters, the relationship between S-parameters and equivalent circuit parameters (R, L, C) is established, and the key coefficients are identified.

[0132] (6) Combined Proton and Gamma Irradiation Experiment

[0133] First, the radiation fluence is E12 P / cm 2 Proton irradiation experiments were conducted, followed by gamma radiation experiments with different radiation doses (100 krad, 300 krad, 500 krad, 900 krad).

[0134] The electrical characteristics of the experimental chip were tested, and the scattering parameters were measured using a network analyzer (VNA) to ensure the accuracy of the data.

[0135] (7) Constructing a model relating proton radiation dose to scattering parameters:

[0136] Using a polynomial fitting method, the relationship between radiation dose and scattering parameters was established based on experimental data, and the linear and nonlinear influence coefficients were determined.

[0137] S11=dD 3 +eD 2 +fD+g;

[0138] Where S11 is the reflection coefficient, D is the radiation dose, d is the cubic coefficient characterizing the higher-order nonlinear effect of the radiation dose on the reflection coefficient S11, e is the quadratic coefficient reflecting the nonlinear effect of the radiation dose on the reflection coefficient S11, f is the linear coefficient representing the degree of linear effect of the radiation dose on the reflection coefficient S11, and g represents the initial reflection coefficient value under no-irradiation conditions.

[0139] (8) Construct a correlation model between radiation dose and equivalent circuit parameters

[0140] Based on the correlation model between scattering parameters and RLC parameters, a relationship model between radiation dose and resistance, inductance, and capacitance in the equivalent circuit is further established.

[0141] (9) Model Validation

[0142] Equivalent circuit parameters under different radiation dose conditions were used as input to perform circuit simulation in the ADS simulation environment to obtain scattering parameters. Error analysis was then performed between the simulation results and experimental data to verify the accuracy of the model. R... 2 This is used to measure the degree of fit between the model and the experimental data. The following is R... 2 Calculation method:

[0143]

[0144] Where SSE represents the squared difference between the model's predicted value and the actual value, the formula is:

[0145]

[0146] SST represents the squared difference between a data point and the data mean, and the formula is:

[0147]

[0148] in, The predicted value of scattering parameter S21, y, is fitted to the correlation model between gamma radiation dose and scattering parameters. i The measured value of scattering parameter S21 obtained from the experiment. The mean value of the scattering parameter S21 obtained from the experiment.

[0149] (10) Multi-frequency and multi-dose integrated analysis

[0150] Comprehensive simulations and analyses were performed under multiple frequency and dose conditions in HFSS and ADS environments. The applicability of the model under different electrical conditions was verified, providing comprehensive data support for the reliability design of three-dimensional microsystem modules.

[0151] Table 1 shows the experimentally determined relationship between electrical resistance and radiation dose:

[0152]

[0153] Table 1. Relationship between electrical resistance and radiation dose as determined experimentally.

[0154] Table 2 shows the experimentally determined relationship between inductance and RDL resistance and radiation dose:

[0155]

[0156] Table 2. Relationship between experimentally determined inductance and RDL resistance and radiation dose.

[0157] The data in Tables 1 and 2 can be used to determine the coefficients. Specifically, after the model is established, the coefficients are determined during the practical operation.

[0158] Figure 3 and Figure 4 The relationship between the fitting results and actual measurements is shown, demonstrating that the TSV-based microsystem interconnect module proton irradiation simulation modeling method has a high degree of fitting and good simulation effect, thus facilitating prediction.

[0159] Example 2

[0160] This embodiment provides a system for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV, including:

[0161] The finite element model building module is used to obtain the geometric parameters and material properties of through-silicon vias, redistribution layers, and microbumps in a three-dimensional microsystem. Based on the geometric parameters and material properties, a finite element model of the three-dimensional microsystem containing through-silicon vias, redistribution layers, and microbumps is built using finite element analysis software.

[0162] The module for building a correlation model between material property parameters and scattering parameters is used to adjust the material parameters of the finite element model of a three-dimensional microsystem, obtain the scattering parameters under different material property conditions, and build a correlation model between material property parameters and scattering parameters.

[0163] The module for constructing the correlation model between RLC equivalent circuit and scattering parameters is used in the ADS compilation environment to extract the circuit parameters of through silicon vias, redistribution layers and microbumps in the finite element model and construct the RLC equivalent circuit; based on the RLC equivalent circuit, a correlation model between the RLC equivalent circuit and scattering parameters is established.

[0164] A module for constructing a correlation model between gamma radiation dose and scattering parameters is used to conduct proton irradiation experiments with preset radiation fluence, and also to conduct gamma radiation experiments with different radiation doses to test scattering parameters and obtain experimental data on scattering parameters; based on the experimental data on scattering parameters, a correlation model between gamma radiation dose and scattering parameters is constructed; and,

[0165] The RLC equivalent circuit correlation model construction module is used to construct a correlation model between gamma radiation dose and RLC equivalent circuit, and based on the correlation model between gamma radiation dose and RLC equivalent circuit, construct correlation models between gamma radiation dose and resistance, inductance and capacitance in RLC equivalent circuit respectively.

[0166] As can be seen from the above detailed description of the present invention, the present invention can describe multilayer three-dimensional structures under different conductivity and dielectric constants, achieving high-precision prediction of scattering parameters and significantly improving the accuracy of modeling and simulation. By combining HFSS and ADS simulation tools, a correlation model between scattering parameters and conductivity, and between dielectric constant and equivalent circuit parameters, is systematically constructed.

[0167] This invention quantifies the influence of conductivity and dielectric constant on scattering parameters, enabling the description of electrical characteristic variations under different gamma radiation doses. This method optimizes the reliability design of multilayer three-dimensional structures in radiation environments, providing a scientific basis for module reliability assessment and optimization, and significantly reducing design errors and experimental costs.

[0168] Finally, it should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to the embodiments, those skilled in the art should understand that modifications or equivalent substitutions to the technical solutions of the present invention do not depart from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV, comprising: Step 1: Obtain the geometric parameters and material properties of the through-silicon vias, redistribution layers, and microbumps of the 3D microsystem; Step 2: Based on geometric parameters and material properties, establish a finite element model of a three-dimensional microsystem including through-silicon vias, redistribution layers, and microbumps using finite element analysis software; Step 3: By adjusting the material parameters of the finite element model of the three-dimensional microsystem, the scattering parameters under different material property conditions are obtained, and a correlation model between the material property parameters and the scattering parameters is constructed. Step 4: In the ADS compilation environment, extract the circuit parameters of through silicon vias, redistribution layers, and microbumps from the finite element model, and construct the RLC equivalent circuit; based on the RLC equivalent circuit, establish a correlation model between the RLC equivalent circuit and the scattering parameters. Step 5: Conduct proton irradiation experiments with a preset radiation flux, followed by gamma radiation experiments with different radiation doses to test scattering parameters and obtain experimental data on scattering parameters; based on the experimental data on scattering parameters, construct a correlation model between gamma radiation dose and scattering parameters. Step 6: Construct a correlation model between gamma radiation dose and RLC equivalent circuit, and based on the correlation model between gamma radiation dose and RLC equivalent circuit, construct correlation models between gamma radiation dose and resistance, inductance and capacitance in RLC equivalent circuit respectively.

2. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 1, characterized in that, In step 1, the geometric parameters of the through-silicon via include: radius, length, and insulating layer parameters; the geometric parameters of the redistribution layer include: width and height; the geometric parameters of the microbumps include: size; and the material properties include: electrical conductivity, magnetic permeability, and dielectric constant.

3. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 2, characterized in that, In step 2, the finite element analysis software used is HFSS simulation software, specifically including: using HFSS simulation software, based on geometric parameters and material properties, preset boundary conditions and port excitation, to establish a finite element model of a three-dimensional microsystem containing through silicon vias, redistribution layers and microbumps, wherein the boundary conditions include the electrical conductivity, dielectric constant and magnetic permeability of the material; the frequency range of the port excitation is set to 0.1-26.5 GHz.

4. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 3, characterized in that, Step 3 specifically includes: adjusting the material property parameters of the through-silicon vias, redistribution layers, and microbumps in the HFSS simulation software; extracting scattering parameters under different material property conditions; and constructing a correlation model between the material property parameters and scattering parameters at a fixed frequency point using a polynomial fitting method; wherein, The model for conductivity S1 is: S1 = aμ 2 + bμ + c The model for the dielectric constant S2 is: S2 = aε 2 + bε + c Where a is the first coefficient, b is the second coefficient, c is the constant term, μ is the relative permeability of the material, and ε is the relative permittivity of the material.

5. The method for simulating and modeling proton and gamma radiation of a TSV-based microsystem interconnect module according to claim 4, characterized in that, Step 4, which involves extracting the equivalent circuit parameters of the redistribution layer, includes: extracting the total resistance R of the rewiring layer of the rewiring layer RDL : wherein R RDL,dc is the DC resistance of the rewiring layer, R RDL,ac is the AC resistance of the rewiring layer; wherein, wherein p RDL is the resistivity of the material of the re-wiring layer, l RDL is the length of the re-wiring layer, w RDL is the wire width of the re-wiring layer, t RDL is the wire thickness of the re-wiring layer, and σ RDL is the conductivity of the material used for the re-wiring layer, wherein Where μ0 is the free permeability, μ r,RDL S represents the relative permeability of the redistribution layer material. RDL This represents the distance between adjacent redistribution layers. Extract the capacitor C of the insulating layer below the redistribution layer. RDLtoSub : Where ε0 is the vacuum permittivity, ε RDLox1 is the dielectric constant of the insulating layer below the redistribution layer, and lox is the length of the insulating layer; Extracting the redistribution layer in silicon equivalent capacitance C RDLinSub : ε r,eff h is the effective dielectric constant. eff For effective height, among which, Where h is the measurement height of the redistribution layer.

6. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 4, characterized in that, Step 4, which involves extracting the equivalent circuit parameters of through-silicon vias and microbumps, includes the following steps: Extracting the DC resistance R of a through-silicon via TSV,dc : Extracting the AC resistance R of a through-silicon via TSV,ac : Where, ρ TSV h is the resistivity of the material in a through-silicon via. TSV r is the height of the through-silicon via. TSV Let δ be the radius of the through-silicon via. skindepth,TSV Indicates the depth of current in a through-silicon via; Extracting the DC resistance R of the micro-bump bump,dc : Extracting the AC resistance R of the micro-bumps BUMPac : Where, ρ BUMP h is the resistivity of the material with micro-bumps. Bump r is the height of the micro-bump. Bump Let δ be the radius of the micro-bump. skindepth,BUMP This indicates the depth of the current in the microbump; Extracting the inductance L of the micro-bump Bump : Among them, h TSV r is the height of the through-silicon via. TSV p is the radius of the through-silicon via. TSV The distance between two adjacent silicon via pillars is denoted by f, and f is the signal frequency when the circuit is in use. Among them, C TSV-ax ε represents the total oxide capacitance of the through-silicon via. TSVox Let r be the dielectric constant of the oxide layer material of the through-silicon via. dep t is the outer radius of the maximum depletion layer of the oxide layer. ox C represents the structural thickness of the oxide layer. ox For silicon dioxide dielectric capacitors with oxide layers, C dep This is the depletion region capacitance when the depletion layer thickness is at its maximum in a silicon substrate with an oxide layer.

7. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 1, characterized in that, In step 5, the preset radiation fluence is E12P / cm. 2 Different radiation doses were included: 100 krad, 300 krad, 500 krad, and 900 krad; scattering parameters were tested using a network analyzer; using a polynomial fitting method, based on the experimental data of the scattering parameters, the relationship between gamma radiation dose and scattering parameters was established, linear and nonlinear influence coefficients were determined, and a correlation model between gamma radiation dose and scattering parameters was constructed. S11=dD 3 +eD 2 +fD+g; Where S11 is the reflection coefficient, D is the radiation dose, d is the cubic coefficient characterizing the higher-order nonlinear effect of the radiation dose on the reflection coefficient S11, e is the quadratic coefficient reflecting the nonlinear effect of the radiation dose on the reflection coefficient S11, f is the linear coefficient representing the degree of linear effect of the radiation dose on the reflection coefficient S11, and g represents the initial reflection coefficient value under no-irradiation conditions.

8. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 1, characterized in that, Step 6 is followed by a verification step: using the equivalent circuit parameters under different gamma radiation dose conditions as input, performing circuit simulation in the ADS simulation environment to obtain the simulation results of the scattering parameters; performing error analysis between the simulation results of the scattering parameters and the experimental data of the scattering parameters to verify the accuracy of the correlation model between gamma radiation dose and scattering parameters.

9. The method for simulating and modeling proton and gamma radiation of a microsystem interconnect module based on TSV according to claim 8, characterized in that, The error analysis includes the following steps: The degree of matching (R) between the correlation model of gamma radiation dose and scattering parameters and the experimental data of scattering parameters was obtained. 2 : Where SSE represents the predicted value of scattering parameter S21 fitted by the correlation model between gamma radiation dose and scattering parameters. The measured value y of the scattering parameter S21 obtained from the experiment i The squared difference between them, SST represents the measured value of the scattering parameter S21 obtained in the experiment. i The measured mean of the scattering parameter S21 obtained from the experiment The difference of squares between them; where, 10. A proton and gamma radiation simulation modeling system for a TSV-based microsystem interconnect module, comprising: The finite element model building module is used to obtain the geometric parameters and material properties of through-silicon vias, redistribution layers, and microbumps in a three-dimensional microsystem. Based on geometric parameters and material properties, a finite element model of a three-dimensional microsystem including through silicon vias, redistribution layers, and microbumps was established using finite element analysis software. The module for building a correlation model between material property parameters and scattering parameters is used to adjust the material parameters of the finite element model of a three-dimensional microsystem, obtain the scattering parameters under different material property conditions, and build a correlation model between material property parameters and scattering parameters. The module for constructing the correlation model between RLC equivalent circuit and scattering parameters is used in the ADS compilation environment to extract the circuit parameters of through silicon vias, redistribution layers and microbumps in the finite element model and construct the RLC equivalent circuit; based on the RLC equivalent circuit, a correlation model between the RLC equivalent circuit and scattering parameters is established. A module for constructing a correlation model between gamma radiation dose and scattering parameters is used to conduct proton irradiation experiments with preset radiation fluence, and also to conduct gamma radiation experiments with different radiation doses to test scattering parameters and obtain experimental data on scattering parameters; based on the experimental data on scattering parameters, a correlation model between gamma radiation dose and scattering parameters is constructed; and, The RLC equivalent circuit correlation model construction module is used to construct a correlation model between gamma radiation dose and RLC equivalent circuit, and based on the correlation model between gamma radiation dose and RLC equivalent circuit, construct correlation models between gamma radiation dose and resistance, inductance and capacitance in RLC equivalent circuit respectively.

Images