Next Article in Journal
Polyhydroxylated Nanosized Graphite as Multifunctional Building Block for Polyurethanes
Next Article in Special Issue
Synthesis of Silver Nanocomposites for Stereolithography: In Situ Formation of Nanoparticles
Previous Article in Journal
Synthesis of Rigid Polyurethane Foams Incorporating Polyols from Chemical Recycling of Post-Industrial Waste Polyurethane Foams
Previous Article in Special Issue
Alumina Manufactured by Fused Filament Fabrication: A Comprehensive Study of Mechanical Properties and Porosity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modeling the Enhanced Efficacy and Curing Depth of Photo-Thermal Dual Polymerization in Metal (Fe) Polymer Composites for 3D Printing

1
Medical Photon Inc., New Taipei City 242, Taiwan
2
Department of Electrical and Engineering, National Taiwan University, Taipei 100, Taiwan
3
CNRS, Université de Haute-Alsace, F-68100 Mulhouse, France
4
Department of Nuclear Medicine and PET Center, China Medical University Hospital, Taichung 400, Taiwan
5
Department of Healthcare Administration, Asia University, Taichung City 413, Taiwan
6
Department of Biomedical Imaging and Radiological Science, China Medical University, Taichung 400, Taiwan
*
Authors to whom correspondence should be addressed.
Submission received: 7 February 2022 / Revised: 4 March 2022 / Accepted: 8 March 2022 / Published: 14 March 2022
(This article belongs to the Special Issue Polymers and Polymer-Based Composites for Additive Manufacturing)

Abstract

:
This article presents, for the first time, the efficacy and curing depth analysis of photo-thermal dual polymerization in metal (Fe) polymer composites for 3D printing of a three-component (A/B/M) system based on the proposed mechanism of our group, in which the co initiators A and B are Irgacure-369 and charge–transfer complexes (CTC), respectively, and the monomer M is filled by Fe. Our formulas show the depth of curing (Zc) is an increasing function of the light intensity, but a decreasing function of the Fe and photoinitiator concentrations. Zc is enhanced by the additive [B], which produces extra thermal radical for polymerization under high temperature. The heat (or temperature) increase in the system has two components: (i) due to the light absorption of Fe filler and (ii) heat released from the exothermic photopolymerization of the monomer. The heat is transported to the additive (or co-initiator) [B] to produce extra radicals and enhance the monomer conversion function (CF). The Fe filler leads to a temperature increase but also limits the light penetration, leading to lower CF and Zc, which could be overcome by the additive initiator [B] in thick polymers. Optimal Fe for maximal CF and Zc are explored theoretically. Measured data are analyzed based on our derived formulas.

Graphical Abstract

1. Introduction

To overcome the poor mechanical properties of pure polymers, metal polymer composites (MPCs) that combine the functionalities of metals and the advantages of polymers have been proposed [1,2,3,4,5,6,7]. However, MPCs still suffer some problems such as extremely poor metal–polymer compatibility, high metal filler content, and functional singularity [8,9,10]. Light-emitting diodes (LED) photoinitiated polymerization offers many advantages such as its solvent-free formulation, mild conditions (ambient temperature, without monomer purification or stabilizers removal etc.), fast reaction rates, and alignment with the green chemistry [11,12,13,14,15]. Copper complex photoredox catalyst photocatalysts have been reported by our group for free radical/cationic hybrid photopolymerization [16,17]. Nevertheless, very few polymeric composites containing metal fillers are prepared by the photopolymerization process due to the poor light penetration depth, which also limits the curing depth.
We have recently reported rimethylolpropane triacrylate (TMPTA) as the model monomer, iron powder as the filler, and Irgacure 369 (Irg 369) as the photoinitiator with the aid of two charge–transfer complexes (CTCs) as dual thermal/photochemical initiators using LED at 405 nm, in which the lack of light penetration was indirectly overcome with the help of CTC in association with heat release during the radical photopolymerization [18,19]. In the dual thermal/photo polymerization process, we have successfully achieved 10 times more depth of curing than that of photopolymerization alone. Our experiment [19] showed that the depth of curing (Zc) is an increasing function of the light intensity but a decreasing function of the Fe and photoinitiator concentration. Furthermore, Zc is enhanced by the additive (CTC), which produces extra thermal radical for polymerization under a high temperature.
In supporting our recent measured data, Bonardi et al. [17] and Ma et al. [18,19], this article will present the kinetics, the conversion profiles features, and the depth of curing of a three-component (A/B/M) system based on the proposed mechanism of our group [19]. We will demonstrate the conversion function (CF). Analytic formulas for the CF and depth of curing and heat (temperature) profiles will be derived rigorously. Our measured data, Ma et al. [19], will be analyzed based on derived formulas and the numerical curves.

2. Methods and Modeling Systems

Photochemical Kinetics

In association with the proposed scheme of our group published in Ma et al. [19], we propose, in Figure 1, a three-component system (A/B/M) defined by the co-initiators [A] and [B] with monomer (M), in which the initiator [A] is excited to its first-excited state A* and a triplet excited state T having a quantum yield (q). The triplet state T interacts with [A] to produce radical (R), which couples with M for photopolymerization and releases heat (H). Additional heat (H) is produced by the UV light absorption of the filler, iron mixed in the monomer. Heat is also produced by the light absorption of the iron filled in the monomer. The heat is transported to the additive (or co-initiator) [B] to produce extra radical R’. Therefore, the polymerization can be induced by the photo-radical (R) and enhanced by the thermal-radical (R’). A specific measured system related to Figure 1 was reported by Ma et al. [19] for a three-component system of PI/B/M, where PI is the photoinitiator Irgacure 369, B is a co-initiator of the charge–transfer complex (CTC), and the M is the monomer rimethylolpropane triacrylate (TMPTA), which is filled with filler iron (Fe).
A second example was proposed by our group published in Bonardi et al. [17], in which A is an absorbing phosphine dye of an infrared light, B is an iodonium salt as the co-initiator, and BlocBuilder MA is the thermal initiator. The light heating in the system of Ma et al. is due to the light absorption of Fe filler in the monomer TMPTA, whereas it is due to the direct light absorption of the thermal initiator in the second example. In both systems, the heat (or temperature increase) in the monomer has two components: the light heater and the heat released from the exothermic photopolymerization of the monomer.
The kinetic equations for our previous systems [20,21,22] are revised for the new system of Figure 1 (A/B/H/M) as follows. We will use the short-hand notations for the concentrations: [A] for photoinitiator Irgacure 369, [B] for the additive CTC, F for the filler iron (in wt%), and R and R’ for the photo and thermal radical, respectively.
d A d t = b I A + R E G
d B d t = k t B
d T d t = b I A K + k A + k M + P [ O 2 ) T  
d R d t = k A T K M R P [ O 2 ] R V R 2
d R d t = k t B K M R
In Equations (3) and (4), we have also included the oxygen inhibition effects (OIH), the P(z)[O2] term with P(z) = P0 exp (−Dz), and D is the oxygen diffusion constant. We note that the OIH reduces the active radical (R), leading to a delayed rising profile of the conversion efficacy for systems in air [23]. P’(z) has similar z dependence as that of P(z), i.e., near the sample surface (with small z) has higher OIH than large depth.
In Equations (2) and (5), k’(t) is the time-dependent decomposition rate of B to produce thermal radical R’, given by an Arhennius formula k’(t) = f exp[−Ed/(rH(t)]; f is the frequency factor for thermal initiator decomposition, Ed is the activation energy, and r is a gas constant. H is the heat (or temperature) in the monomer due to two components: (i) the aFI(z,t) term, due to light absorption of iron filler given by the absorption constant (a), light intensity (I), and iron concentration (F); and (ii) the heat released from the exothermic photopolymerization of the monomer (the KMR term) and (iii) the heat transport to the co-initiator, B. Therefore, H is given by a generalized heat diffusion equation
d H z , t d t = a F I z , t + K R M b H + c 2 H
where b’ is the heat transport constant, and the last term is the thermal conduction term in depth(z).
The UV light intensity has the depth (z) dependence due to the absorption of PI, [A], and iron filler (F) given by
  d I z , t d z = b A + a F I z , t
The monomer (M) conversion rate equation has three components: triplet state (T) direct coupling.
Under the above quasi-steady-state solutions and for unimolecular case for R, we obtain the simplified equations (with the thermal conduction term ignored) as follows.
d A d t = b I A R E G  
d B d t = k t B
d H d t = a F I + k A T b H
The conversion efficacy (CF) defined by CF = 1 − M(t)/M0, with M(t) given by the solution of
d M d t = k M + k A / 1 + O I H   T k B

3. Methods and Results

A full numerical simulation is required for the solutions of Equations (9) to (12), which will be presented elsewhere. We will focus on comprehensive analysis for special features and the key factors of dual radical enhancement for efficient conversion related to the measured data of Ma et al. [19], based on the analytic formulas to be derived as follows.

3.1. Analytic Results

In general, Equations (11) and (12) require numerical simulation. For analytic formulas, we further approximate g = (1 − D)/(k[A]), with D(t) = k”M(t)/(k[A]) being the second-order correction term, such that T = bI(1 − D)/(k[A]); Equations (11) and (12) are reduced to
d H d t = a F + b I z , t b H b D t I z , t
d M d t = 1 / 1 + O I H + k / k M / A 1 D   b I z , t k t B
Analytic formulas need further assumption of RGE taken as a mean reduction factor (f’), such that bI[A] − RGE = f’bI[A] having a value of f’ = 0 to 1.0. The first-order solution of Equation (9) is given by [A] = A0 exp(−dt), with d = bI(z), for f’ = 1, and I(z,t) assumed to be time-independent, that is, the increase of I(z,t) due to the depletion of [A] is neglected. For the perfect recycle case, f = 0, and [A] = A0 is a constant. Using an average of [A(z,t)] over z and t (defined as A’), we obtain I(z) = I0 exp[−Gz], with G = aF + bA’.
H z , t = H 0     e x p b t + d / b 1 e x p b t Q t
where H0 is the initial value of H(t = 0), and d’ = (aF + b)I(z)= (aF + b)I0 exp[−Gz]. We have also added a correction factor Q(t) = k”M0/(kA0) [1 − exp(−D’t)]D’ = qt (for small t), with D’ = 1 + k”. We note that H(t) has a transient state value H0 (1 − b’t) + (d’ − q)t and a steady state value (aF + b)I(z)/b’. Moreover, for z > 0, d’ has an optimal value at F = F*, and it is a decreasing function of F, for F > F*. This optimal value is due to the competing term (aF) in d’ = (aF + b)I0 exp[−(aF + bA’)]. Therefore, a higher F > F* and/or A’ leads to less light penetration due to the absorption factor G = aF + bA’ in I(z) = I0 exp[−Gz], and also leads to lower heat released (temperature) in the monomer, H(t).
Using Equation (15) to solve for Equation (10), with approximated k’(t) = f exp[−Ed/(rH(t)] = f[1 − Ed/(rH(t)] = f[1 − (1 − Q’t)Ed/(rH0)], with Q’ = (d’ − q)/H0 − b’, we obtain
B z , t = B 0     e x p P t  
where P(t) = f(P’ + 0.5P”t)t, with P’ = 1 − Ed/(rH0), P” = Ed/(rH0)Q’.
Using the first-order solution of Equation (14), with D = OIH = k’(t) = 0, and M(t) = M0 exp(−k”dt) and [A] = A0 exp(−dt), with d = bI, we find the second-order solution (including k’[B]) of Equation (14) given by
M(t) = M0 − [dt + (k”/k)(M0/A0)V(t)] − k’ B0V’(t)
where V(t) = [1 − exp(−st)]/s, with s = (1 + k”)bI; V’(t) = [1 − exp(−fP’t)]/(fP’), with P’ = 1 − Ed/(rH0).
Two special cases are analyzed for Equation (17) as follows.
(a) For steady-state: V(t) = 1/s, V’(t) = 1/(fP’), we obtain the conversion efficacy (CF) defined by CF = 1 − M(t)/M0,
CF(z,t) = bIt/M0 + k”/[(1 + k”)kA0] + k’ B0/(fP’M0)
We note CF is proportional to the light dose (It) and that this feature is only for the case of unimolecular termination of the radical R, which is proportional to (bI). For the case of bimolecular [20], R is proportional to the power of (bI) 0.5, and the steady-state value of CF is proportional to (bI)−0.5.
(b) For transient state, V(t) = V’(t) = t, we obtain the CF
CF(z,t) = b(It)[1/M0 + k”/(kA0)] + k’tB0/M0.
We note that Equation (17) has three terms proportional to the light dose E = I(z)t, and initial concentrations 1/A0, 1/M0, and B0 and the absorption constants, b. However, CF (for z > 0) is a decreasing function of Fe, since I(z) = I0 exp[−(aF + bA’)z].

3.2. Depth of Curing (DoC)

As shown by Equation (7) the light intensity is an exponentially decreasing function of the depth(z) according to a Beer–Lambert law (BLL). However, for the situation of time-dependent intensity I(z,t) due to the depletion of the PI concentration, A(t) = A0 exp(−dt), a generalized BLL was developed by Lin [20], such that the solution of Equation (8) becomes I(z,t) = I0 exp[−Gz], in which G is an average of [A(z,t)] over z and t (defined as A’), given by G = aF + bA’.
A curing depth (ZC) is defined by when the conversion efficacy is higher than a critical value, CF > C*. Using I(z,t) = I0 exp[−Gz], and Equation (18) for CF, we obtain
Z C = ln S I 0 / C * S   / a F + b A
where S = bt[1/M0 + k”/(kA0)]; and the enhancement factor, S’ = k’t(B0H0)/M0.
We note that ZC is proportional to the light dose I0t, but it is a decreasing function of the effective absorption constant (aF + bA’), with A’ being the z and time-averaged PI concentration [A], given by A’. Equation (19) shows that Zc is enhanced by the additive [B] via S’ = k’t(B0H0)/M0, which is proportional to k’B0H0, i.e., the limited light penetration due to Fe is overcome by the additive CTC, or [B]. Equation (19) shows that Zc is proportional to ln(SI0) under the unimolecular case, with neglected VR2 term in Equation (4). For the bimolecular dominant case (with VR2 >> KRM), Equation (19) is revised to [21]
Z C = 2 ln S I 0 0.5 / C * S   / a F + b A
We note that more detailed discussion and numerical results of Zc and curing time may be found in our previous work [20,22].
As shown by Equation (15), when Fe is too high, the DoC is too low to create effective polymerization, and it is limited to thin films. For lower Fe of 10% to 30%, the DoC is larger than the sample thickness, and sharp edges are produced. Moreover, high Fe leads to strong light absorption, and the light dose is smaller than the polymerization threshold at a certain sample depth.

3.3. General Features and New Findings

As shown by Equations (9) to (18), the following are the significant features of the [A]/[B]/M system, which was also proposed by Ma et al. [19]. Our modeling has explored the following general features, in which some new findings are not yet explored in our previous experiment of Ma et al. [23], but might be justified in our future works.
(i) As shown by Equation (6), the heat (or temperature) increase in the system has two components: (a) the aFI(z,t) term, due to light absorption of iron filler given by the absorption constant (a), light intensity (I), and iron concentration (F); and (b) the heat released (the KMR term) from the exothermic photopolymerization of the monomer.
(ii) The solution of heat H(t) given by Equation (14) has two terms, which have opposite trends in their time profiles; the first term is a decaying function, whereas the second term is an increasing profile. We note that H(t) has a transient state value H0 (1 − b’t) + d’t, and a steady state value (aF + b)I/b’. Therefore, for F > F*, higher F and/or A’ leads to less light penetration due to the term G = aF + bA’ in I(z) = I0 exp[−Gz], and lower heat released (temperature) in the monomer. This feature of optimalization is valid only for z > 0. It does not exist for surface (with z = 0). Numerical simulation will be shown later for this optimal feature.
(iii) As shown by Equations (8) and (12), (14), the monomer conversion has three components: (a) from the direct coupling of T and [A], (b) the coupling of photon-radical (R) with M, and (c) coupling of thermal-radical (R’) and M, in which R’ is an enhancement radical produced from the co-initiator, [B], under a high temperature given by the heat function, H(t). We note that the Fe filler leads to temperature increase but also limits the light penetration, leading to lower CF and DoC, which could be overcome by the additive initiator [B] for thick samples.
(iv) As shown by Equation (17), the CF has terms proportional to the light dose E = It, and initial concentrations A0, H0 and the last term of Equation (17), k’B0/M0, which is the enhancement factor in the presence of co-initiator B. The above features are for the case with unimolecular coupling of KRM term in Equation (4). For the situation of a strong bimolecular termination, we need to include the coupling of VR2 in Equation (4), and leading a CF that is proportional to the (bI)0.5 (for transient state) and inverse proportional to (bI)−0.5 at steady-state, i.e., higher intensity leads to a lower steady-state CF than lower intensity. The enhancement factor is given by S’ = k’t(B0H0)/M0, shown in Equations (17) and (18). The CF enhancement in the presence of co-initiator B was also reported experimentally by Bonardi et al. [17].
(v) As show by Equation (19), the depth of curing (DoC), Zc, is an increasing function of the light intensity. However, there is an optimal Zc, based on the parameter of (aF + b), i.e., there is an optimal value of the iron concentration (F) that can be found mathematically by taking the derivative of Zc over d’, and let it equal 0, which can be also found numerically [20]. This optimal value is due to the competing term in d’ = (aF + b)I0 exp[−(aF + A’)]. The enhanced Zc is given by S’ = k’t(B0H0)/M0. Equation (18) also defines the minimal enhancement factor (S*) to reach a sample depth of Z’ given by the condition of when (aF + b)Z’ > ln (SI0/(C* − S’), or S’ > S*, with S* = C* − SI0 exp[−(aF + b)Z’], or k’t(B0H0)/M0 > C* − SI0 exp[−(aF + b)Z’].
(vi) Our modeling also shows that there is an optimal condition for the initiator PI concentration given by d(CF)/dF = 0, with CF given by Equation (17), which also requires numerical stimulation [20].

3.4. Theoretical Predictions, Numerical Data, and Future Directions

We will now present numerical data based on our derived formulas, Equation (14) for H(t) and Equation (18) for DoC as follows. We also carefully choose the ranges of each key parameters, such that their roles could be well-demonstrated.
Based on Equation (14), we show, in Figure 2, the volume (at z = 1.0 cm) temperature profiles in the monomer in the absence of PI (with b = 0), with heat produced by light absorption of Fe only (or b = 0). It shows that larger Fe mixed in the monomer leads to lower temperature profiles due to the light penetration depth decreases as shown by I(z,t) = I0 exp(−Gz), with G = aF + bA’, which is a strong decreasing function of G. The saturation is due to the heat transport term, b’H.
Figure 3 shows the temperature profiles (at z = 1.0 cm) in the presence of co-initiators, [A] and [B], with aFI0 = (0, 0.1, 0.2) (1/s), for a fixed bI0 = 10 (1/s), shown by curves (1, 2, 3), which have much higher temperature than the case of pure Fe absorption (with b = 0), curve-4; Curves 1 to 3 also have much faster increasing profiles due to the strong heat transport term of b’H (with b’ = 0.035), compared to the weak coupling curve-4 (with b’ = 0.002). The drop of Curves 1 to 3 is due to Q(t) = qt in Equation (14), with q = 1.5 for Figure 3.
Figure 4 shows the optimal feature for a maximal heat produced by the light absorption of Fe, based on Equation (14). The optimal feature exists only for z > 0, due to the light intensity given by I(z,t) = I0 exp(−Gz), with G = aF + bA’, which is a strong decreasing function of (aF + bA’), for z > 0. These features are show in Figure 4, for z = (0, 0.4, 0.5, 0.6) cm.
Based on our formula, Equation (18), we also show, in Figure 5, the DoC (Zc) as the function of various parameters. Figure 5A shows that Zc is an increasing function of S’, the enhancement factor given by S’ = k’t(B0H0)/M0, for a given curing threshold C* = 0.4, and fixed SI0 = 0.7, aF + bA’ = 1.0. Figure 5B shows that Zc is a decreasing function of the factor (aF + bA’), Figure 5C shows that Zc is an increasing function of SI0, which is proportional to the light dose tI0, and the coupling constant, b, as shown by Equation (19).
Based on our modeling predictions described above, we propose to conduct the following new measurements, which are not yet explored in [24].
(i) CF profiles for various initiator concentrations, A0 and B0, to justify our predicted enhancement factor given by S, in Equations (14) and (18).
(ii) CF profiles for various light intensity (I0) to justify the steady-state CF scaling law of (bI0)m, with m = 1.0 and −0.5;
(iii) curing depth (Zc) for various enhanced factor given by S’ = k’t(B0H0)/M0, as shown by Figure 5.
(iv) the optimal feature shown by Figure 4, for various aFI0, for a maximal heat produced by the light absorption of Fe.
Besides the enhancement effects discussed in the present article using an additive of CTC, there are many other strategies that have been reported by our group, including the use of multi-wavelength light and strategies of reduction the oxygen inhibition [24,25,26,27,28]. A critical review of photopolymerization for 3D bioprinting can be found in [24].

3.5. Analysis of Measured Data

The measured data shown in the images (figures) from Ref. [19] will be analyzed as follows, where a PI Irgacure 369 is exposed to a LED at 405 nm (0.11 W cm−2) for various iron (Fe) concentrations.
Image 1A from Ref. [19] shows the measured conversion profiles of TMPDA for various Fe concentrations. These profiles can be analyzed by our CF shown in Equation (18), which shows a typical saturation profile defined by Q’(t) = [1 − exp(−b’t)]/b’ in Equation (17).
Image 1B from Ref. [19] shows the measured DoC for various Fe concentrations. These profiles can be analyzed by our Zc shown in Equation (19), which shows Zc is an increasing function of the light intensity I0 but a decreasing function of the Fe (F) and PI concentration (A’), shown by the term (aF + bA’). The DoC is enhanced by the additive [B] via the S function, which is proportional to (bt)[(1/M0 − A0) + k’B0/(b’M0)]. Our measured data had yielded (10 times) thick composites under the same dose using a LED irradiation at 405 nm.
Image 1C from Ref. [19] shows the temperature profiles in the absence of PI. They are related to the Equations (12) and (14) with b = 0, but b’ > 0. This feature is also predicted by our theory as shown in Figure 2.
Image 1A from Ref. [19] shows the temperature profiles for various Fe concentrations. These profiles can be analyzed by Equation (14), which shows that H(t) has a transient state value H0 (1 − b’t) + d’t and a steady state value d’/b’, with d’ = (aF + b)I(z), with I(z) = I0 exp[−Gz]. Therefore, higher F and/or A’ leads to less light penetration due to the absorption factor G = aF + bA’ and also leads to a lower heat released (temperature) in the monomer, as shown by the measure curves of Image 1D of Ref. [19]. These features are also predicted by our theory, as shown in Figure 3. The CF enhancement in the presence of co-initiator B was also reported by Bonardi [17].
We found that wt10% and wt30% Fe filler resin can achieve a perfect three-dimensional structure, as shown in Image 5A and 5B from Ref. [19] with great spatial resolution for the photoproducts in the presence of the iron fillers. However, when the content of iron keeps up to 50% wt%, the pattern was not flat, and some area under irradiation was not fabricated (shown by image 5C of Ref. [19]). See more detailed of measurement method in Ref. [19]. We note that possibly metal plus polymer composites are discussed in [28].

4. Conclusions

We have demonstrated theoretically and justified by measured CF profiles of a three-component (A/B/M) system based on the proposed mechanism of our group, Ma et al. [19], in which the co initiators are A, Irgacure 369, and B, charge–transfer complexes (CTC), and the monomer M is TMPDA. Our formulas show the depth of curing (Zc) is an increasing function of the light intensity but a decreasing function of the Fe and photoinitiator concentration. Zc is enhanced by the additive [B], which produces an extra thermal radical for polymerization under high temperature. The heat (or temperature) increase in the system has two components: due to light absorption of iron filler, and heat released from the exothermic photopolymerization of the monomer. The heat is transported to the additive (or co-initiator) [B] to produce extra radical R’ and enhance the monomer conversion function (CF).

Author Contributions

Conceptualization, J.-T.L. and J.L.; data curation, C.-H.K., K.-H.L. and D.-C.C.; formal analysis, J.-T.L.; funding acquisition, D.-C.C.; investigation, J.L., C.-H.K., K.-H.L. and D.-C.C.; methodology, J.-T.L. and Y.-Z.L.; project administration, J.L. and D.-C.C.; resources, J.L. and D.-C.C.; software, Y.-Z.L.; supervision, J.-T.L., J.L. and D.-C.C.; validation, J.-T.L., Y.-Z.L., C.-H.K. and D.-C.C.; writing—original draft, J.-T.L.; writing—review and editing, J.-T.L. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by China Medical University under grant number CMU110-ASIA-11. And The APC was funded by CMU110-ASIA-11.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

J.-T.L. is the CEO of Medical Photon Inc.

References

  1. Patel, D.; Sakhaei, A.H.; Layani, M.; Zhang, B.; Ge, Q.; Magdassi, S. Highly Stretchable and UV Curable Elastomers for Digital Light Processing Based 3D Printing. Adv. Mater. 2017, 29, 1606000. [Google Scholar] [CrossRef]
  2. Narupai, B.; Nelson, A. 100th Anniversary of Macromolecular Science Viewpoint: Macromolecular Materials for Additive Manufacturing. ACS Macro Lett. 2020, 9, 627–638. [Google Scholar] [CrossRef]
  3. Xu, Z.; Gao, C. In situ Polymerization Approach to Graphene-Reinforced Nylon-6 Composites. Macromolecules 2010, 43, 6716–6723. [Google Scholar] [CrossRef]
  4. Zhang, Y.; Josien, L.; Salomon, J.-P.; Simon-Masseron, A.; Lalevée, J. Photopolymerization of Zeolite/Polymer-Based Composites: Toward 3D and 4D Printing Applications. ACS Appl. Polym. Mater. 2020, 3, 400–409. [Google Scholar] [CrossRef]
  5. Zhang, Y.; Xu, Y.; Simon-Masseron, A.; Lalevée, J. Radical photoinitiation with LEDs and applications in the 3D printing of composites. Chem. Soc. Rev. 2021, 50, 3824–3841. [Google Scholar] [CrossRef] [PubMed]
  6. Liu, Y.; Lin, Y.; Jiao, T.; Lu, G.; Liu, J. Photocurable modification of inorganic fillers and their application in photopolymers for 3D printing. Polym. Chem. 2019, 10, 6350–6359. [Google Scholar] [CrossRef]
  7. Xin, Y.; Peng, H.; Xu, J.; Zhang, J. Ultrauniform Embedded Liquid Metal in Sulfur Polymers for Recyclable, Conductive, and Self-Healable Materials. Adv. Funct. Mater. 2019, 29, 1808989. [Google Scholar] [CrossRef]
  8. Lei, Z.; Wu, P. A supramolecular biomimetic skin combining a wide spectrum of mechanical properties and multiple sensory capabilities. Nat. Commun. 2018, 9, 1–7. [Google Scholar] [CrossRef] [PubMed]
  9. Dadashi-Silab, S.; Doran, S.; Yagci, Y. Photoinduced Electron Transfer Reactions for Macromolecular Syntheses. Chem. Rev. 2016, 116, 10212–10275. [Google Scholar] [CrossRef]
  10. Corrigan, N.; Yeow, J.; Judzewitsch, P.; Xu, J.; Boyer, C. Seeing the Light: Advancing Materials Chemistry through Photopolymerization. Angew. Chem. Int. Ed. Engl. 2019, 58, 5170–5189. [Google Scholar] [CrossRef] [PubMed]
  11. Garra, P.; Caron, A.; Mousawi, A.; Graff, B.; Morlet-Savary, F.; Dietlin, C.; Yagci, Y.; Fouassier, J.-P.; Lalevée, J. Photochemical, Thermal Free Radical, and Cationic Polymerizations Promoted by Charge Transfer Complexes: Simple Strategy for the Fabrication of Thick Composites. Macromolecules 2018, 51, 7872–7880. [Google Scholar] [CrossRef]
  12. Garra, P.; Dietin, F.; Morlet-Savary, B.; Dumur, F.; Gigmes, D.; Fouassier, J.P.; Lalevée, J. Photopolymerization processes of thick films and in shadow areas: A review for the access to composites. Polym. Chem. 2017, 8, 7088–7101. [Google Scholar] [CrossRef]
  13. Xu, Y.; Ding, Z.; Zhu, H.; Zhang, Y.; Knopf, S.; Xiao, P.; Lalevée, J. Preparation of Iron Filler-Based Photocomposites and Application in 3D Printing. Macromol. Mater. Eng. 2021, 306, 2000720. [Google Scholar] [CrossRef]
  14. Garra, P.; Graff, B.; Morlet-Savary, F.; Dietlin, C.; Becht, J.M.; Fouassier, J.P.; Lalevee, J. Charge Transfer Complexesas Pan-Scaled Photoinitiating Systems: From 50 μm 3D Printed Polymers at 405 nm to Extremely Deep Photopolymerization (31 cm). Macromolecules 2018, 51, 57–70. [Google Scholar] [CrossRef]
  15. Mokbel, H.; Anderson, D.; Plenderleith, R.; Dietlin, C.; Morlet-Savary, F.; Dumur, F.; Gigmes, D.; Fouassier, J.-P.; Lalevée, J. Copper photoredox catalyst “G1”: A new high performance photoinitiator for near-UV and visible LEDs. J. Polym. Chem. 2017, 8, 5580–5592. [Google Scholar] [CrossRef]
  16. Mokbel, H.; Anderson, D.; Plenderleith, R.; Dietlin, C.; Morlet-Savary, F.; Dumur, F.; Gigmes, D.; Fouassier, J.; Lalevée, J. Simultaneous initiation of radical and cationic polymerization reactions using the “G1” copper complex as photoredox catalyst: Applications of free radical/cationic hybrid photopolymerization in the composites and 3D printing fields. Prog. Org. Coat. 2019, 132, 50–61. [Google Scholar] [CrossRef]
  17. Bonardi, A.H.; Bonardi, F.; Noirbent, G.; Dumur, F.; Gigmes, G.; Lalevee, J. Free-Radical Polymerization upon Near-Infrared Light Irradiation, Merging Photochemical and Photothermal Initiating Methods. J. Polym. Sci. 2020, 58, 300–308. [Google Scholar] [CrossRef]
  18. Ma, Q.; Liu, S.; Le Dot, M.; Mokbel, H.; Zhang, Y.; Graff, B.; Lalevée, J. Imidazole based dual photo/thermal initiators for highly efficient radical polymerization under air with a metal-free approach. Polym. Chem. 2021, 12, 6386–6391. [Google Scholar] [CrossRef]
  19. Ma, Q.; Zhang, Y.; Launay, V.; Le Dot, M.; Liu, S.; Lalevée, J. How to overcome the light penetration issue in photopolymerization? An example for the preparation of high content iron-containing opaque composites and application in 3D printing. Eur. Polym. J. 2022, 165, 111011. [Google Scholar] [CrossRef]
  20. Lin, J.-T.; Cheng, D.-C. Modeling the efficacy profiles of UV-light activated corneal collagen crosslinking. PLoS ONE 2017, 12, e0175002. [Google Scholar] [CrossRef] [PubMed]
  21. Lin, J.T. Efficacy S-formula and kinetics of oxygen-mediated (type-II) and non-oxygen-mediated (type-I) corneal cross-linking. Ophthalmol. Res. 2018, 8, 1–11. [Google Scholar] [CrossRef]
  22. Lin, J.-T.; Liu, H.-W.; Chen, K.-T.; Cheng, D.-C. Modeling the Kinetics, Curing Depth, and Efficacy of Radical-Mediated Photopolymerization: The Role of Oxygen Inhibition, Viscosity, and Dynamic Light Intensity. Front. Chem. 2019, 7, 760. [Google Scholar] [CrossRef] [PubMed]
  23. Lin, J.-T.; Lalevee, J.; Cheng, D.-C. A Critical Review for Synergic Kinetics and Strategies for Enhanced Photopolymerizations for 3D-Printing and Additive Manufacturing. Polymers 2021, 13, 2325. [Google Scholar] [CrossRef] [PubMed]
  24. Lin, J.; Cheng, D.; Chen, K.; Chiu, Y.; Liu, H. Enhancing UV Photopolymerization by a Red-light Preirradiation: Kinetics and Modeling Strategies for Reduced Oxygen Inhibition. J. Appl. Polym. Sci. 2020, 58, 683–691. [Google Scholar] [CrossRef]
  25. Lin, J.-T.; Chen, K.-T.; Chang, P.-K.; Cheng, D.-C. Enhancing blue-light-initiated photopolymerization in a three-component system: Kinetic and modeling of conversion strategies. J. Polym. Res. 2021, 28, 1–9. [Google Scholar] [CrossRef]
  26. Lin, J.-T.; Cheng, D.-C.; Chen, K.-T.; Liu, H.-W. Dual-Wavelength (UV and Blue) Controlled Photopolymerization Confinement for 3D-Printing: Modeling and Analysis of Measurements. Polymers 2019, 11, 1819. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Lin, J.T.; Liu, H.W.; Chen, K.T.; Cheng, D.C. 3-wavelength (UV, blue, red) controlled photopolymerization: Improved conversion and confinement in 3D-printing. IEEE Access 2020, 8, 49353–49362. [Google Scholar] [CrossRef]
  28. Brostow, W.; Hagg Lobland, H.E. Materials: Introduction and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2016. [Google Scholar]
Figure 1. The schematics of a dual function system, (A/B/H). The heat (H) can be produced by the light absorption of the iron filled in the monomer (M) and the heat released by the photopolymerization, in which R and R’ lead to a photo-thermal dual polymerization (see text for more details).
Figure 1. The schematics of a dual function system, (A/B/H). The heat (H) can be produced by the light absorption of the iron filled in the monomer (M) and the heat released by the photopolymerization, in which R and R’ lead to a photo-thermal dual polymerization (see text for more details).
Polymers 14 01158 g001
Figure 2. The calculated temperature (at z = 1.0 cm) profiles of monomer based on Equation (15) for various weight concentration of Fe, with aFI0 = (0.3, 1.0, 1.6) (1/s) and b’ = 0.002 (1/s), in the absence of PI (with b = 0), with heat produced by light absorption of Fe only, with initial temperature H0 = 20 °C. The arrow indicates the direction of increased values of aFI0.
Figure 2. The calculated temperature (at z = 1.0 cm) profiles of monomer based on Equation (15) for various weight concentration of Fe, with aFI0 = (0.3, 1.0, 1.6) (1/s) and b’ = 0.002 (1/s), in the absence of PI (with b = 0), with heat produced by light absorption of Fe only, with initial temperature H0 = 20 °C. The arrow indicates the direction of increased values of aFI0.
Polymers 14 01158 g002
Figure 3. Same as Figure 2, but in the presence of co-initiators, [A] and [B], with aFI0 = (0, 0.5, 1.0, 1.5) (1/s), for Curves 1, 2, 3 (with bI0 = 10) and Curve-4 in the absence of PI (with b = 0). The arrow indicates the direction of increased values of aFI0.
Figure 3. Same as Figure 2, but in the presence of co-initiators, [A] and [B], with aFI0 = (0, 0.5, 1.0, 1.5) (1/s), for Curves 1, 2, 3 (with bI0 = 10) and Curve-4 in the absence of PI (with b = 0). The arrow indicates the direction of increased values of aFI0.
Polymers 14 01158 g003
Figure 4. The calculated temperature of monomer versus G = aF + bA’ based on Equation (15) at various sample depth z = (0, 0.4, 0.5, 0.6) cm, for curves (1, 2, 3, 4), at t = 200 s.
Figure 4. The calculated temperature of monomer versus G = aF + bA’ based on Equation (15) at various sample depth z = (0, 0.4, 0.5, 0.6) cm, for curves (1, 2, 3, 4), at t = 200 s.
Polymers 14 01158 g004
Figure 5. The calculated DoC (Zc) based on Equation (19) for various parameters: (A) for Zc vs. S’ with fixed SI0 = 0.7, aF + bA’ = 1.0 and C* = 0.4; (B) for Zc vs. aF + bA’, with fixed SI0 = 0.5 and S’ = 0.5; (C) for Zc vs. SI0 = 0.5, with fixed and S’ = 0.2 and aF + bA’ = 1.0.
Figure 5. The calculated DoC (Zc) based on Equation (19) for various parameters: (A) for Zc vs. S’ with fixed SI0 = 0.7, aF + bA’ = 1.0 and C* = 0.4; (B) for Zc vs. aF + bA’, with fixed SI0 = 0.5 and S’ = 0.5; (C) for Zc vs. SI0 = 0.5, with fixed and S’ = 0.2 and aF + bA’ = 1.0.
Polymers 14 01158 g005
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lin, J.-T.; Lee, Y.-Z.; Lalevee, J.; Kao, C.-H.; Lin, K.-H.; Cheng, D.-C. Modeling the Enhanced Efficacy and Curing Depth of Photo-Thermal Dual Polymerization in Metal (Fe) Polymer Composites for 3D Printing. Polymers 2022, 14, 1158. https://0-doi-org.brum.beds.ac.uk/10.3390/polym14061158

AMA Style

Lin J-T, Lee Y-Z, Lalevee J, Kao C-H, Lin K-H, Cheng D-C. Modeling the Enhanced Efficacy and Curing Depth of Photo-Thermal Dual Polymerization in Metal (Fe) Polymer Composites for 3D Printing. Polymers. 2022; 14(6):1158. https://0-doi-org.brum.beds.ac.uk/10.3390/polym14061158

Chicago/Turabian Style

Lin, Jui-Teng, Yi-Ze Lee, Jacques Lalevee, Chia-Hung Kao, Kuan-Han Lin, and Da-Chuan Cheng. 2022. "Modeling the Enhanced Efficacy and Curing Depth of Photo-Thermal Dual Polymerization in Metal (Fe) Polymer Composites for 3D Printing" Polymers 14, no. 6: 1158. https://0-doi-org.brum.beds.ac.uk/10.3390/polym14061158

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop