• 沒有找到結果。

Synthesis, properties and pyrolysis of siloxane- and imide-modified epoxy resin cured with siloxane-containing dianhydride

N/A
N/A
Protected

Academic year: 2021

Share "Synthesis, properties and pyrolysis of siloxane- and imide-modified epoxy resin cured with siloxane-containing dianhydride"

Copied!
6
0
0

加載中.... (立即查看全文)

全文

(1)

DOI: 10.1002/pi.1863

Synthesis, properties and pyrolysis of

siloxane- and imide-modified epoxy resin

cured with siloxane-containing dianhydride

Hsun-Tien Li, Hsun-Ren Chuang, Ming-Wei Wang and Mu-Shih Lin

Department of Applied Chemistry, National Chiao-Tung University, Hsinchu, 30050, Taiwan

Abstract: Hydrosilylation of nadic anhydride with tetramethyl disiloxane yielded 5,5 -(1,1,3,3-tetramethyl disiloxane-1,3-diyl)-bis-norborane-2,3-dicarboxylic anhydride (I), which further reacted with 4-aminophenol to give N,N-bis(4-hydroxyphenyl)-5,5-bis-(1,1,3,3-tetramethyl di1,3-diyl)-bis-norborane-2,3-dicarboximide (II). Epoxidation of II with excess epichlorohydrin formed a siloxane-and imide-modified epoxy oligomer (ie diglycidyl ether of N,N-bis(4-hydroxyphenyl)-5,5 -bis(1,1,3,3-tetramethyl disiloxane-1,3-diyl)-bis-norborane-2,3-dicarboximide) (III). Equivalent ratios of III/I of 1/1 and 1/0.8 were prepared and cured to produce crosslinked materials. Thermal mechanical and dynamic mechanical properties were investigated by TMA and DMA, respectively. It was noted that each of these two materials showed a glass transition temperature (Tg) higher than 160◦C with moderate moduli. The thermal degradation kinetics was studied with dynamic thermogravimetric analysis (TGA) and the estimated apparent activation energies were 111.4 kJ mol−1(in N2), 117.1 kJ mol−1(in air) for III/I= 1/0.8, and 149.2 kJ mol−1(in N2), 147.6 kJ mol−1(in air) for III/I= 1/1. The white flaky residue of the TGA char was confirmed to be silicon dioxide, which formed a barrier at the surface of the polymer matrix and, in part, accounted for the unique heat resistance of this material.

 2005 Society of Chemical Industry

Keywords: siloxane- and imide-modified epoxy oligomer; crosslinking; dynamic mechanical properties; thermal degradation kinetics; apparent activation energy; SiO2char residual

INTRODUCTION

In the development of electronic devices, electronic packaging has a tendency towards lower cost, finer pitch, higher electrical performance and better reliability. Flip chip technology prevails among chip packaging candidates to meet these trends.1 – 9 The

widespread use of epoxy resins in the flip chip package interests people in the development of higher performance for encapsulants. Epoxy resin exhibits many desirable properties, such as high strength and modulus, excellent chemical and solvent resistance, good thermal and electrical properties, outstanding adhesion, and easy processing. However, cured epoxy resins are generally hard and relatively brittle because of their rigid crosslinked networks.

Modification with elastomers such as carboxyl-terminated butadiene-acrylonitrile copolymer (CTBN) and silicone rubber has been considered as a successful method to improve the toughness of epoxy resin after many efforts, while there are inevitable drawbacks in mechanical properties and thermal sta-bility. Moreover, if epoxy is directly blended with sili-cone elastomer, phase separation occurs because of the incompatibility of the two components, and this wors-ens the mechanical properties.10 – 17 Consequently,

we are interested in developing modified epoxy to enhance the flexibility and toughness while sustain-ing the thermo-mechanical properties with improvsustain-ing dielectric constants. In order to improve the mechan-ical properties, modification of silicone-epoxy with a rigid segment, such as an imide group, is essential owing to the superior thermo-mechanical properties obtained: high crosslinking ability, high glass transi-tion temperature (Tg), high thermal stability, high char

yield, and excellent heat resistance and mechanical modulus. In this study, we synthesized and character-ized a novel siloxane- and imide-modified epoxy resin and cured it with siloxane-containing dianhydride for possible electronic applications. The thermal mechan-ical properties and thermal degradation kinetics of the cured epoxy were also investigated. The main pyroly-sis product of the materials was characterized as SiO2

in this study.

EXPERIMENTAL Materials

5-Norbornene-2,3-dicarboxylic anhydride (nadic anhydride), and tetramethyl disiloxane were purchased from Merck Co. Platinum divinyl tetramethyl

Correspondence to: Mu-Shih Lin, Department of Applied Chemistry, National Chiao-Tung University, Hsinchu, 30050, Taiwan

E-mail: [email protected]

Contract/grant sponsor: National Science Council of Taiwan; contract/grant number: NSC 93-2216-E-009-020 (Received 17 November 2004; revised version received 9 April 2005; accepted 22 April 2005)

(2)

disiloxane complex was purchased from UCT Co. Epichlorohydrin was purchased from TCI. All chem-icals were purified before use. Electronic grade of diglycidyl ether of bisphenol-F (830LVP) was obtained from DIC Co, and 2-ethyl-4-methyl-1-ethylcyanoimidazole (2E4MZCN) was purchased from Shikoku Chemicals Co and used as the catalyst for curing the epoxy and anhydride.

Synthesis of 5,5-(1,1,3,3-tetramethyl disiloxane-1,3-dilyl)-bis-norborane-2,3-dicarboxylic

anhydride (I)

The synthesis of siloxane-containing dianhydride (I) was performed according to the method reported in the literature (Scheme 1).18 The chemical

struc-ture of the white product was confirmed with proton and 13C nuclear magnetic resonance spec-troscopy (Unity 300 MHz NMR, Varian Inc., Wal-nut Creek, CA, USA). 1H-NMR (CDCl3, ppm) δ:

0.03 – 0.05 (m, 12H), 0.65 (t, 2H), 1.55 – 1.66 (m, 8H), 2.73 – 2.78 (m, 2H), 2.84 (m, 2H), 3.39 – 3.43 (m, 4H). 13C-NMR (CDCl3, ppm) δ: −1.20–0.94, 25.75, 26.74, 40.32, 41.00– 41.03, 41.64, 49.46, 52.62, 171.98 – 172.25. Synthesis of N,N-bis(4-hydroxylphenyl)-5,5 -bis(1,1,3,3-tetramethyl disiloxane-1,3-diyl)-bis-norborane-2,3-dicarboximide (II)

The siloxane-containing dianhydride (I) (30 g, 0.065 mol) was dissolved in 90 ml of dimethyl-formamide (DMF), and 4-aminophenol (14.9 g, 0.136 mol) in 40 ml of DMF was gradually added (Scheme 1). The solution was stirred for 8 h at ice-bath conditions, followed by imidization using a Dean-Stark instrument. Compound II in solid powder was obtained after vacuum drying, and it was recrystallized from isopropanol (38 g, yield= 91 %, mp = 123◦C). An Avatar 360 FT-IR (Nicolet Instrument Corpora-tion, Madison, WI, USA) was used for IR analysis.

1H-NMR (CDCl 3, ppm) δ: 0.01 –0.02 (m, 12H), 0.61 (m, 2H), 1.54 – 1.62 (m, 8H), 2.74 (m, 2H), 2.78 (m, 2H), 3.12 – 3.17 (m, 4H), 6.70 – 6.73 (d, 4H), 6.90 – 6.94(d, 4H), 7.42 (s, 2H). 13C-NMR (CDCl3, ppm) δ:−0.88, −0.84, 25.64, 26.54, 39.70, 40.57, 41.33, 48.54, 51.09, 116.16, 123.45, 127.81, 156.62, 178.57, 178.60, 178.82. IR: imide 1789 and 1720 cm−1, OH 3100 – 3500 cm−1(broad band).

Synthesis of diglycidyl ether of N,N -bis(4-hydroxylphenyl)-5,5-bis(1,1,3,3-tetramethyl disiloxane-1,3-diyl)-bis-norborane-2,3-dicar-boximide oligomer (III)

The synthesis of siloxane- and imide-containing epoxy (III) was carried out according to a method described previously.19 A three-necked flask equipped with

an addition funnel was filled with a solution of II (12.88 g, 0.02 mol) and epichlorohydrin (37 g, 0.4 mol) (Scheme 1). A solution of NaOH (1.6 g NaOH, 0.04 mol in 3.2 ml of H2O) was added

drop-wise to the mixture, which was stirred for 6 h at 65◦C. The product was extracted three times with toluene. The solvent and residual reactants were removed by vacuum drying. The product was then purified by column chromatography using n-hexane/ethyl acetate (v/v= 1/1) as eluent to obtain an oligomer of siloxane-and imide-containing epoxy (III) with an epoxy equiv-alent weight (EEW) of 477 g eq−1 by titration with HBr, and a calculated average repeating unit of 0.34 (III in Scheme 1, where n= 0.34; yield = 92.6 %; mp= 100–103◦C). Elemental analysis for III: (calcu-lated) C= 64.84 %, N = 3.78 %, H = 6.53 %, O = 17.27 %; (found) C= 64.00 %, N = 3.78 %, H = 6.73 %, O= 17.16 %. IR: absorption of epoxide at 910 cm−1. 1H-NMR (CDCl 3, ppm) δ: 0.01 –0.16, 0.63, 1.54 – 1.64, 2.69 – 2.89, 3.20 – 3.21, 3.29 – 3.30, 3.92 – 3.94, 4.18 – 4.22, 6.93 – 6.97, 7.09 – 7.12. 13 C-NMR (CDCl, ppm) δ: −0.85, 25.62, 26.54, 39.76, Si O CH3 CH3 Si CH3 CH3 H H O O O Si O Si CH3 CH 3 CH3 CH3 O O O O O O + 2 Pt toluene, 80 °C I Si O Si CH3 CH3 CH3 CH3 N N O O O O CHCH2O OCH2CHCH2O H 2C O OH Si O Si CH3 CH3 CH3 CH3 N N O O O O OCH2CH CH2 O n III Si O Si CH3 CH3 CH3 CH3 O O O O O O H 2N OH + 2 −H2O Benzene/DMF Si O Si CH3 CH3 CH3 CH3 N N O O O O HO OH I II Si O Si CH3 CH3 CH3 CH3 N N O O O O HO OH H 2C O CH2Cl + 20 NaOH / H2O methyl ethyl ketone, 80 °C

II

(3)

40.62, 41.33, 44.50, 48.50, 49.90, 51.07, 68.95, 115.17, 125.00, 127.72, 158.30, 177.35 – 177.57. Preparation of epoxy curing systems

The siloxane- and imide-modified epoxy oligomer (III) and the siloxane-containing dianhydride (I) in molar ratios of 1/1 and 1/0.8, together with 1 phr of 2-ethyl-4-methyl-1-ethylcyanoimidazole, were dissolved in acetone and then vacuum dried at 50◦C to obtain well-mixed blends. For comparison, control systems of similar blend with bisphenol-F epoxy in place of III were prepared exactly by the same formulation and procedure.

Characterization

Thermal mechanical properties

The epoxy curing systems were poured onto alu-minium plates and cured according to the profile 130◦C for 30 min, 150◦C for 180 min and 200◦C for 30 min. The thermal mechanical behaviour was inves-tigated by TMA (TA 2940, DuPont) and the dynamic mechanical properties were analysed using a DuPont DMA2980 with a dual cantilever head on cured film with dimensions of 60 mm× 8 mm × 0.1 mm at a fre-quency of 1 Hz and heating rate of 3◦C min−1 from 30◦C to 250◦C.

Thermal degradation kinetics

The thermal degradation kinetics of the cured epoxy were investigated using thermogravimetric analysis (TGA) (DuPont TA 2950). The system was operated in the dynamic mode over the temperature range 50 – 800◦C, at various heating rates: 5, 10, 15 and 20◦C min−1. All the experiments were carried out under dry nitrogen and under air atmospheres

RESULTS AND DISCUSSION Thermal mechanical properties

The thermal mechanical properties of these four curing materials are shown in Fig 1 and summarized in Table 1. It should be noted that the Tg values of

all these systems were higher than 160◦C. Dynamic mechanical properties

Figures 2 and 3 compare DMA data for the four crosslinked materials. It is noted that the Tgvalues of

the III/I systems showed a consistent tendency with the data obtained from TMA (Fig 1). In addition, the III/I

systems showed higher damping peaks (tan δ), higher area of tan δ peaks and lower loss modulus than those of bisphenol-F epoxy/I systems. Among these systems, the one with III/I= 1/0.8 had the lowest storage modulus E but showed the highest damping peak (tan δmax), and also the highest area of tan δ. All of the

DMA data confirmed the stress – release, enhancement of flexibility and toughness of the III/I systems. Although III/I systems had lower Tg values than

the corresponding bisphenol-F/I systems (Figs 1, 2 and 3), all these systems exhibited excellent flexibility and toughness, because of the proper balance of soft siloxane and rigid imide segments in the polymer matrices. In addition, it was found that both the tan δ and Ecurves of the III/I systems had a large main transition peak with a small second transition, implying

1900 900 −100 Dimension change (µ m) 0 50 100 150 200 250 300 Temperature (°C) A: III/I=1/1 B: III/I=1/0.8 C: 830LVP/I=1/0.8 D: 830LVP/I=1/1 A B C D

Figure 1. TMA thermograms of 830LVP/I and III/I systems.

50 100 150 200 250 300 −0.2 0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 C D A B tan d Temperature (°C) A: III/I=1/1 B: III/I=1/0.8 D: 830LVP/I=1/0.8 C: 830LVP/I=1/1

Figure 2. Damping curves of 830LVP/I and III/I systems.

Table 1. Comparisons of thermal mechanical properties of III/I and 830LVP/I systems

System Tga tan δ (◦C) Peak heightb Area of tan δ Eglass(MPa) Erubber(Mpa)

III/I= 1/1 173.2 185.6 0.970 36.90 1903.7 6.604

III/I= 1/0.8 171.1 187.6 0.996 37.76 1122.0 6.604

830LVP/I= 1/1 204.1 227.7 0.163 14.37 2012.5 92.17

830LVP/I= 1/0.8 201.9 237.2 0.128 14.91 1436.9 95.92

aT

gvalues obtained from the TMA thermograms.

(4)

50 100 150 200 250 300 −200 0 200 400 600 800 1000 1200 1400 1600 1800 2000 2200 2400 D C B A Storage modulus (MP α ) Temperature (°C) A: III/I=1/1 B: III/I=1/0.8 D: 830LVP/I=1/0.8 C: 830LVP/I=1/1

Figure 3. Storage moduli of 830LVP/I and III/I systems.

compatibility of the system with possibly some phase separation. This small second transition may also possibly have been due to the segmental motion of the flexible siloxane. Nevertheless, it was noted that the siloxane moiety was well dispersed in the crosslinked network and it was hard to see phase separation in the cured systems. This inference was confirmed by the homogeneity of carbon-mapping scanning electron microscopy (SEM), oxygen-mapping SEM and silicon-mapping SEM of the cured systems (Fig 4). All the mapping SEM results indicated homogeneity and compatibility of the cured matrices. Thermal degradation behaviour

The TGA thermograms of cured epoxy provided valuable information regarding their thermal stability and thermal degradation behaviour. TGA data of these systems (Table 2 and Fig 5) indicated that 5 % weight loss in nitrogen atmosphere for III/I=

Figure 4. Typical mapping SEM images of III/I= 1/1: (upper left)

SEM photo; (upper right) carbon-mapping; (lower left) oxygen-mapping; (lower right) silicon-mapping.

100 80 60 40 20 0 Weight (%) 60 260 460 660 Temperature (°C) A: III/I=1/1 10 °C min−1 in N 2 C: III/I=1/0.8 10 °C min−1 in N 2

B: III/I=1/1 10 °C min−1 in air

D: III/I=1/0.8 10 °C min−1 in air

C D A

B

Figure 5. TGA thermograms of III/I system.

Table 2. Dynamic TGA and activation energy Eaaby the Kissinger method

Atmosphere Heating rate (◦C min−1) 5 % weight loss T (◦C) Peak 1 T (◦C) Yield (%) Ea1 (kJ mol−1) Peak 2 T (◦C) Yield (%) Char yield (%) Ea2 (kJ mol−1) III/I= 1/0.8 N2 20 365.9 408.4 84.17 111.4 642.0 36.09 14.35 75.1 15 350.3 396.3 82.99 610.3 29.29 12.44 10 344.6 384.3 83.12 593.3 34.29 18.81 5 320.7 366.9 81.65 543.3 33.68 23.22 Air 20 347.0 394.9 79.35 117.1 593.8 31.02 16.49 72.9 15 357.5 400.6 81.06 634.2 33.00 18.30 10 346.3 382.1 85.33 579.1 34.35 18.46 5 321.8 365.1 81.94 557.2 28.85 18.26 III/I= 1/1 N2 20 356.6 401.3 78.45 149.2 576.1 35.46 19.13 78.0 15 350.3 399.7 78.85 602.1 28.71 24.94 10 342.1 386.1 80.12 553.4 35.10 24.71 5 319.7 373.2 74.67 532.7 27.80 24.62 Air 20 356.9 397.5 81.14 147.6 596.1 31.11 18.01 82.7 15 350.2 386.0 83.53 586.2 28.47 18.04 10 317.1 382.8 74.60 558.3 29.11 16.59 5 310.2 365.6 76.64 517.1 32.81 16.76 aEa= activation energy

(5)

5800 5600 5400 5200 5000 4800 4600 4400 4200 4000 3800 3600 3400 3200 3000 2800 2600 2400 2200 2000 1800 1600 1400 1200 1000 800 600 400 200 0 0 Counts C 1 2 3 4 5 6 7 8 9 10 11 12 13 14 Energy (keV) S1 Au Au Au

Figure 6. ESCA spectrum showing the TGA residual of the III/I system, identified as silicon dioxide; the residue itself is shown in the SEM photo.

1/0.8 occurred in the range 320.7 –365.7◦C, and in the range 319.7 – 356.4◦C for III/I= 1/1. The 5 % weight loss in air atmosphere for III/I= 1/0.8 occurred 321.7 – 357.5◦C, and at 310.2 – 356.9◦C for III/I= 1/1. Furthermore, flaky residuals with char yields in the range 12 – 24.9 % were found as the pyrolysis products. The main composition of the residual was characterized to be silicon dioxide, as shown in the SEM and electron spectroscopy for chemical analysis (ESCA) spectrum in Fig 6. All of the experimental results indicated that the transformation of these silicon-containing polymers into stable silicon dioxide could form a flaky layer on the polymer surface, which is believed to inhibit the propagation of decomposition, and thus to improve the heat resistance of the cured materials.

Thermal degradation kinetics

Dynamic TGA has frequently been used to study the thermal degradation kinetics of polymers. We applied Kissinger’s method20 – 22 to determine the activation

energy (Ea) of solid – gas reactions using the following

equation: ln  β Tmax2  = ln  AR Ea  + ln(n(1 − αmax)n−1)Ea RTmax

where β is the heating rate, Tmax is the temperature

of the inflection point of the thermograms that

corresponds to the maximum reaction rate, A is the pre-exponential factor, αmax is the maximum

conversion, n is the reaction order, and R is the molar gas constant (8.314 J K−1mol−1). The straight lines were obtained according to the least-squares method for both systems (Fig 7). The dynamic method proposed herein gave apparent activation energies of degradation of 111.4 kJ mol−1(in N2), 117.1 kJ mol−1

(in air) for III/I= 1/0.8 system, and 149.2 kJ mol−1 (in N2), 147.6 kJ mol−1(in air) for III/I= 1/1 system.

The higher the siloxane and imide content in the crosslinked material, the higher the activation energy. Moreover, there was a unique phenomenon that the apparent activation energies obtained in air were very close to those obtained in nitrogen atmosphere. It can be inferred that the formation of thermally stable silicon dioxide creates a barrier on the material surface and prevents further decomposition in the inner polymer matrix, thus explaining this experimental finding.

CONCLUSION

A novel siloxane- and imide-modified epoxy oligomer (III) was successfully synthesized, and characterized by elemental analysis, IR,1H-NMR, and13C-NMR.

This novel epoxy oligomer was then cured with a siloxane-containing dianhydride (I). Unlike with most silicone-containing compounds, there was no apparent phase separation in the cured epoxy. The TMA and DMA data revealed that the flexible siloxane group

(6)

(a) R2= 0.9899 −12 −11 −10 −9 1.48 1.49 1.52 1.53 1.55 1/T × 10−3 (K−1) ln (R /T 2) 1/T × 10−3 (K−1) (b) R2 = 0.9599 −12 −11 −10 −9 1.49 1.52 1.52 1.57 ln (R/ T 2)

Figure 7. Typical dynamic TGA of III/I systems: (a) III/I= 1/1 in N2;

(b) III/I= 1/1 in air; R2corresponds to the correlation of the points.

and the rigid imide group improved the toughness as well as the mechanical properties. The thermal degradation kinetic data, the SEM micrographs and the ESCA results showed that during pyrolysis the siloxane groups had been transformed into flaky silicon dioxide as char residual, which further improved its heat resistance by inhibiting the propagation of thermal decomposition in the polymer matrix. This evidence accounts for the fact that this siloxane- and imide-containing epoxy material exhibited good thermal stability with similar thermal degradation behaviour in both N2and air atmospheres.

ACKNOWLEDGEMENT

The authors thank the National Science Council of Taiwan for financial support under contract number NSC 93-2216-E-009-020.

REFERENCES

1 Chen KC, Nemoto N, Huang SC, Kitamura K, Tsuji T,

Fukui T, Li HT and Lee TM, EPTC International Confer-ence, p 141 (2003).

2 Zhang ZQ, Shi SH and Wong CP, IEEE Trans Components Packaging Technol 24:59 (2001).

3 Luo SJ and Wong CP, International Symposium on Advanced Packaging Materials, p 183 (2000).

4 Yin MJ and Paik KW, Electronic Components and Technology Conference, p 899 (2000).

5 Li H, Wang L, Jacob K and Wong CP, J Polym Sci, Part A: Polym Chem 40:1796 (2002).

6 Jang MS and Crivello JJ, J Polym Sci, Part A: Polym Chem

41:3056 (2003).

7 Cheng KC, Yu SY and Chiu WY, J Appl Polym Sci 83:274 (2002).

8 Chiang CL and Ma CC, Eur Polym J 38:2219 (2002). 9 Lin RH, J Polym Sci Part, A: Polym Chem 38:2934 (2000). 10 Abad MJ, Barral L, Cano J, Lopez J, Nogueira P, Ramirez C

and Torres A, Eur Polym J 37:1613 (2001).

11 Lee JY, Shim MJ and Kim SW, J Appl Polym Sci 81:479 (2001).

12 Kumar AA, Alagar M and Rao RMVGK, Polymer 43:693 (2002).

13 Ochi M and Shimaoka S, Polymer 40:1305 (1999). 14 Yee AF and Pearson RA, J Mater Sci 21:2475 (1986). 15 Bucknall CB and Partridge IK, Polym Eng Sci 26:54 (1986). 16 Bascom WD, Cottington RL, Jones RL and Peyser P, J Appl

Polym Sci 19:2425 (1975).

17 Sultan JN and McGarry FJ, Polym Eng Sci 13:29 (1973). 18 Eddy VJ, Hallgren JE and Robert E, J Polym Sci, Part A: Polym

Chem 28:2417 (1990).

19 Wang WJ, Perng LH, Hsiue GH and Chang FC, Polymer

41:6113 (2000).

20 Kissinger HE, Anal Chem 29:1703 (1957).

21 Nair CPR, Krishnan K and Ninan KN, Thermochim Acta 39:61 (2000).

22 Yang LF, Yao KD and Koh W, J Appl Polym Sci 73:1501 (1998).

數據

Figure 2. Damping curves of 830LVP/I and III/I systems.
Table 2. Dynamic TGA and activation energy E a a by the Kissinger method
Figure 6. ESCA spectrum showing the TGA residual of the III/I system, identified as silicon dioxide; the residue itself is shown in the SEM photo.
Figure 7. Typical dynamic TGA of III/I systems: (a) III/I = 1/1 in N 2 ;

參考文獻

相關文件

substance) is matter that has distinct properties and a composition that does not vary from sample

To stimulate creativity, smart learning, critical thinking and logical reasoning in students, drama and arts play a pivotal role in the..

 Promote project learning, mathematical modeling, and problem-based learning to strengthen the ability to integrate and apply knowledge and skills, and make. calculated

Wang, Solving pseudomonotone variational inequalities and pseudocon- vex optimization problems using the projection neural network, IEEE Transactions on Neural Networks 17

Then, we tested the influence of θ for the rate of convergence of Algorithm 4.1, by using this algorithm with α = 15 and four different θ to solve a test ex- ample generated as

Particularly, combining the numerical results of the two papers, we may obtain such a conclusion that the merit function method based on ϕ p has a better a global convergence and

Define instead the imaginary.. potential, magnetic field, lattice…) Dirac-BdG Hamiltonian:. with small, and matrix

We investigate some properties related to the generalized Newton method for the Fischer-Burmeister (FB) function over second-order cones, which allows us to reformulate the