Influence of γ -Irradiation on the Electronic Structure and the Chemical and Mechanical Properties of Poly(hydroxybutyrate-valerate)/Poly(caprolactone) Blends: Insights from Experimental Data and Computational Approaches

: In this study, we investigated the influence of γ -irradiation (0, 50, and 100 kGy) doses on the chemical and mechanical properties of biodegradable poly(hydroxybutyrate-valerate)/poly (caprolactone) (PHBV/PCL) polymer blends rich in low-molar-mass PCL, which were prepared using a co-rotating twin-screw extruder. In parallel, the density functional theory (DFT) and the time-dependent DFT (TD-DFT) methods were used together with a model containing four monomer units to provide an insight into the electronic structure, chemical bonds, and spectroscopic (such as Nuclear Magnetic Resonance (NMR) and Ultraviolet-visible (UV-vis)) properties of PHBV and PCL blend phases, which are critical for predicting and designing new materials with desired properties. We found that an increase in γ -irradiation doses caused splitting instead of crosslinks in the polymer chains, which led to evident deformation and an increase in tensile strength at break of 2.0 to 5.7 MPa for the PHBV/PCL blend. Further, this led to a decrease in crystallinity and proved the occurrence of a more favorable interaction between the blend phases.

Various synthesis strategies (e.g., chemoenzymatic methods, enzymatic ring opening polymerization, radical ring opening polymerization, coordinative ring opening polymerization, anionic ring opening polymerization, photo-initiated ring opening polymerization, and so on) have been widely employed in the preparation of a diverse mixture of biodegradable polymer blends with different compositions, phases, and morphologies [7,[9][10][11][12][13][14][15][16][17].Among biodegradable polymeric blends, more recently, poly(hydroxybutyrate-valerate) (PHBV) and poly(caprolactone) (PCL) have garnered extensive interest concerning a wide range of technological applications, thereby providing opportunities to develop emergent technologies in the future [18][19][20][21][22].In particular, PHBV is a copolymer of hydroxybutyrate and valerate monomers and has several interesting properties that make it useful for a variety of applications.For example, it has a high melting point (above 160 • C) and is crystalline in nature, giving it good mechanical properties such as high tensile strength and modulus [23][24][25].One of the main applications of PHBV is in the production of biodegradable plastics [26].These plastics based on PHBV can be used in a variety of applications where traditional petroleum-based plastics would normally be used, but can be broken down by microorganisms over time, reducing the environmental impact of their use [27][28][29].PHBV can also be used to make films, fibers, and other products [30].Other potential applications of PHBV, more recently explored in the literature, include its use in drug delivery systems, tissue engineering, and as a material for 3D printing [31][32][33][34][35][36][37].However, it is well-known that the large-scale production of these biopolymers occurs mainly through a bacterial fermentation process, which remains somewhat costly [38][39][40][41].A rational strategy to reduce production costs and modify the physicochemical properties and degradation rates of such blends is of high value [1,2,42].
On the other hand, PCL is a biodegradable petroleum-based polymer, linear polyester with a low melting temperature (around 65 • C) and good mechanical properties [43][44][45].However, it is well-known that several microbial enzymes, including esterase, cutinase, and lipase, can efficiently break down PCL into valuable chemicals [46][47][48][49].In addition, the hydrolytic degradation rate of PCL is significantly slow, most likely due to its wellknown hydrophobic nature [50,51].It is also biocompatible, making it suitable for use in medical and pharmaceutical applications [52][53][54][55].When blended, PHBV and PCL can form materials with a range of mechanical and thermal properties, depending on the blend ratio.It has been found that an increase in PHBV within the PCL matrix leads to a marked decrease in its miscibility [56].In terms of mechanical properties, PHBV/PCL blends have a combination of the properties of the individual polymers.For example, PHBV has a high melting point and good mechanical properties such as high tensile strength and modulus, while PCL has a lower melting point and good flexibility [57,58].As a result, PHBV/PCL blends can have a balance of high strength and good flexibility, as well as biodegradability, depending on the ratio of the two polymers in the blend.
For instance, Liu et al. have studied the biodegradation assay of PHBV/PCL blends using Pseudomonas mendocina [59].After 96 h of degradation, the PHBV blends undergo extensive decomposition, leaving behind several channel-like structures, while the PCL portion remains apparently intact.Notably, it is well-known that the PCL polymer degrades over months to years depending on its molecular weight, degree of crystallinity, morphology, porosity, and the surrounding environment [22,[60][61][62].
In some cases, biodegradable polymers may be used in applications where they will be exposed to γ-radiation, such as in medical devices or food packaging [63,64].It has been observed that using ionized radiation during the processing of these systems improves thermal stability, miscibility, and other relevant properties [18,65,66].Furthermore, in the polymers that undergo crosslinking when subjected to ionizing radiation, chemical bonds are formed between adjacent molecules, which lead to permanent linking [67].In these cases, it is important to understand how the polymer will be affected by the radiation to ensure that it will still function as intended.The specific effects of γ-radiation on the chemical and mechanical properties of particular biodegradable polymer blends will depend on the structure and composition of each polymer, as well as the specific conditions under which they are irradiated [68][69][70][71].However, the effects of γ-radiation doses on the chemical and mechanical properties of biodegradable PHBV/PCL polymeric blends remain to be elucidated.
Herein, we have tried to understand further the effects of γ-radiation (0, 50, and 100 kGy) doses on the chemical and mechanical properties of PHBV/PCL (25:75) systems rich in low molar mass PCL and prepared by means of a co-rotating twin-screw extruder.This 25:75 ratio for the PHBV/PCL blend is due to its optimal suitability for biomedical applications.These blends were characterized using scanning electron microscopy (SEM), differential scanning calorimetry (DSC), dynamic mechanical thermal analysis (DMTA), and mechanical tensile tests, together with computational studies based on density functional theory (DFT) and time-dependent DFT (TD-DFT).By combining theoretical simulations with experimental techniques, it is possible to better understand the interaction between polymer blends and to predict their stability, reactivity, and spectroscopic behavior in different applications.Additionally, to elucidate the mechanical response of these systems, we evaluated the elastic modulus of these samples by changing the doses of γ-irradiation.Notably, the experimental results revealed the general degradation mechanism of these PHBV/PCL systems.Our theoretical results contribute to a greater physical and chemical understanding of this PHBV/PCL (25:75) blend.

Characterization
The morphological analysis of these samples was conducted using an SEM microscope (PHILIPS XL-30 FEG, Hillsboro, OR, USA) with a secondary electron detector.To obtain SEM micrographs, all the samples were cryogenically fractured from ambient conditions (~20 • C and 960.7 mmHg) and coated with gold.An AT QS100 device was used to conduct DSC analysis at temperatures ranging from ambient temperature (~20 • C and 960.7 mmHg) to 200 • C with a heating rate of 20 • C/min.Analysis of the DSC curves for the raw materials of the PHBV and PCL indicated multiple melting peaks, which were consistent with the previous studies [2,21].However, owing to the high crystallinity of these polymers and the sensitivity limitation of their DSC analysis, the determination of glass transition temperatures for the PHBV and PCL raw samples was not possible.However, in this study, to define the glass transition, we used the temperature corresponding to the peak of the tan δ (T d ) curve.This selection was made because the constant glass transition temperature value was defined by this variable, which can be obtained by the relationship with the elasticity modulus (E ′′ /E ′ ) [73].DMTA analysis was performed using a DMA 2980 Dynamic Mechanical Analyser (TA Instruments, Inc., New Castle, DE, USA) in the three-point bending mode, according to the ASTM D-4065-95 standard [74].This analysis was performed at a temperature range of −80 to 40 • C, a frequency of 1 Hz, and a heating rate of 2 • C/min.Mechanical tensile tests were performed using an Instron 5500R universal testing machine with a load cell of 50 kN; an extensometer was employed to determine the elastic modulus at 0.8% strain.The distance between grips was 105 mm (10), and a rate of 100 mm/min was maintained.The tensile tests were performed with the specimen dimensions and analysis procedure specified in the ASTM D-638 standard, and the results were obtained from the average of the values corresponding to seven samples [72].

Computational Details
The theoretical models selected for this study contain polymer chains with four monomer units prepared for both pure PHBV and PCL (see Figure 1A,B).Next, based on experimental stoichiometry for PHBV/PCL blend (1:3), we have evaluated for our model all connection possibilities between these polymers, such as PHBV-PCL-PCL-PCL (de-noted as blend BA), PCL-PHBV-PCL-PCL (denoted as blend BB), PCL-PCL-PHBV-PCL (denoted as blend BC), and PCL-PCL-PCL-PHBV (denoted as blend BD), as shown in Figure 1C-F.Although the blend structures may look symmetrical (in the case between BA and BD and between BB and BC), they are not equivalent, as is shown by the functional group connecting the fragments and their terminations (e.g., the PCL termination in BA and BD fragments are hydroxide and carbonyl groups, respectively).All of the DFT and TDDFT calculations were performed using the standard B3LYP/6-31+G(d,p) level as implemented in the Gaussian 09 package [75].It is widely recognized that the B3LYP level can introduce errors of up to 5 kcal/mol in thermochemical energy calculations [76][77][78][79].Full geometry optimizations and their frequencies confirmed that all models correspond to a minimum of the potential energy surface (PES).Also, 1 H and 13 C NMR calculations were performed for these structures using the gauge-including-atomicorbital (GIAO) method, where tetramethylsilane (TMS) was used as reference [80,81].In order to provide a better description of the electronic properties of the systems under investigation, single-point calculations at the TDDFT/B3LYP/6-31+G(d,p) level of theory were performed.
Nonconvalent interaction (NCI), electron localization function (ELF), electrostatic surface potential (ESP), NMR, and total density of states (TDOS) analyses were performed with the Multiwfn software utilizing the TDDFT checkpoint and output files [82,83].The NCI and ESP maps were printed on Visual Molecular Dynamics (VMD) with isosurface values of 0.65 and 0.03, respectively [84].As for the electronic properties derived from the energies of the molecular orbitals, i.e., hardness η, softness S, Mulliken electronegativity x, and electrophilicity ω, the values were estimated by using the equations below [85].

Results and Discussion
By examining the surface features of the irradiated sample using SEM micrographs, we can gain insights into how the irradiation has affected the material at the microscale.Figure 2A-C shows the SEM micrographs of PHBV/PCL polymeric blends with and without irradiation.The SEM analysis revealed that the PHBV phase in the blend exhibited distributed surface micropores.Presumably, these micropores originated during the evaporation of waste liquids (water and waste solvent), which were eliminated in processing.In the SEM micrographs of the irradiated samples, the presence of microcracks in the PHBV was observed because of the radiation effect on the splitting of the chains that weakened the polymer.Contrarily, without irradiation, the PHBV/PCL blend did not exhibit any microcracks on its structure.As a result, the extent of chain scission and the resulting weakening of the polymer will depend on the specific type of radiation used, the dose of radiation, as well as the properties of the polymer (or blend) itself.Similar results has been reported by Tubio et al. [86].These authors also highlight the microstructure evolution in a PHBV/PCL blend with compositional variation.
Figure 2A,B shows the DSC curves for pure PHBV and PCL polymers as well as their PHBV/PCL (25:75) blend with and without irradiation, respectively.DSC results were used to determine crystalline melting temperature (Tm), enthalpy of fusion (∆Hm), and percentage of crystallinity (% C), which are summarized in Table 1.The percentages of PHBV and PCL crystallinity were calculated from the theoretical fusion heat data of the hypothetical 100% crystalline PHBV (∆Hm 100% = 146 J/g) and the hypothetical 100% crystalline PCL (∆Hm 100% = 136 J/g) [87,88].In addition, these results suggested a small decrease in the crystalline melting temperature of the γ-irradiated PHBV/PCL blend.Therefore, the influence of irradiation on the PHBV structure resulted in the disintegration of chains instead of crosslinks, which caused an increase in the crystallinity of the pure and blend forms.In addition, the splitting occurred more frequently in large chains, which narrowed the molar mass distribution curve and shifted it to the left.This favored the order and regularity of the chains over a long range [69,89].Contrarily, there was a significant decrease in the crystallinity of the PCL oligomer owing to the low molar mass of PCL.It was observed that under irradiation, there was a larger decrease in molar mass, as well as a clear increase in the number of very small chains, and there may have been a more favorable interaction between the blend phases.DMTA curves for the PHBV/PCL blend (with and without irradiation) are shown in Figure 2F-H.These curves indicate that the irradiation did not affect the mobility of chains in the amorphous phase.As is well known, the glass transition temperature is influenced not only by factors such as a decrease in molecular weight and polymer compatibility but also by the mobility of chains within the amorphous phase [72, [89][90][91][92].In this context, electron paramagnetic resonance (EPR) emerges as a well-established and optimal technique for studying the amorphous phase of polymers [93][94][95].By determining the probe-radical correlation time, EPR offers a nuanced exploration of chain dynamics, providing a more holistic perspective on the molecular behavior within the amorphous regions.Therefore, the use of this methodology will be the target of our future study.According to the results presented in Table 1, in PHBV, rigidity and fragility were observed, which were characterized by the low deformation of the material, its sudden rupture, and the high elastic modulus.These properties were associated with the high crystallinity of PHBV [90].Weakening of the PHBV was prolonged for some time until complete crystallization (~3 days).This phenomenon is known as the 'aging of PHBV'.It has not yet been clearly explained and has been attributed to the occurrence of secondary crystallization in the polymeric matrix [91,92].For the PHBV/PCL blend, after irradiation, there was a slight variation in the values of the polymers' elastic moduli.The tensile strength at break significantly increased for the pure PHBV polymer.For the PHBV/PCL blend, the increase in the irradiation dose from 50 to 100 kGy caused small changes in the tensile strength; in contrast, for the pure PHBV polymer, it caused an increase in the tensile strength.All the materials suffered a significant reduction in deformation at the break with irradiation doses of 50 and 100 kGy.When compared with pure PHBV, the PHBV/PCL blend without irradiation exhibited an increase in the deformation at the break.As can clearly be seen in Figure 1, the pure PHBV polymer has a nonlinear structure due to its methyl and ethyl ramification points, while PCL presents a linear chain with alternate (zigzag) carbonyl groups and exact 180 • dihedrals.This structural nonlinearity leads to an increase in the dipole moment of 2.44 Debye to 3.92 Debye for PCL and PHBV, respectively.In the case of the blend models, the exchange of a PCL monomer unit for a PHBV fragment induces a strong structural polarization, which is shown with the angle disparity between the PHBV and PCL fragments and the dihedrals diverging from the PCL by 180 • around the PHBV/PCL frontier regions, likely resulting in a significant decrease in crystallinity.In this same framework, the BA blend also presents higher stability in terms of total energy, being more stable than other conformations studied exhibiting relative values of +0.32 kcal/mol compared to BB, +0.43 kcal/mol compared to BC, and +0.90 kcal/mol compared to BD.
Figure 3 depicts the ESP surface maps of pure and blend polymers, where the red and blue regions represent the negative and positive charge extremes, respectively.Looking at the PHBV surface, the negative charges are strongly centered on carbonyl oxygens, while the positive charges are localized on the hydrogens in the opposite extremity of the carbonyl group.On the other hand, the PCL polymer has negative charges which are also localized on the carbonyl oxygens and the positive charges on the hydrogens.As expected, the exchange of a PHBV unit in the PCL chain introduces extra charge concentrations around the carbonyl oxygens and hydrogens.As the PHBV fragment is dislocated between the extremities, the blend polarity is shifted due to the present groups in each extremity, which polarize the structure.Hence, the ESP surfaces are overall consistent with the dipolar moments and the total energy shift between the blends, as seen above.
Expanding the discussion regarding the pure and blend polymer structures, 13 C NMR and 1 H spectra for each system are found in Figure 4A,B.Overall, the 1 C NMR spectra are divided into two regions, being assigned to higher chemical shifts (160-195 ppm) within the poorly shielded carbonyl group regions and lower chemical shifts (10-80 ppm) within the strongly shielded fragments' midparts and ramifications, while this difference is more seamless for 1 H spectra. Therefore, due to the ramification groups in the PHBV structure, several extra signals are detected in both 1 H and 13 C NMR spectra in comparison to the PCL polymer.In particular, the peaks around 10, 20, 30, and 35 ppm in the 13 C spectrum are associated with the strong shielding of outer ethyl, methyl, and inner ethyl C nuclei, showing equivalent methyl nuclei and an ethyl C that is less shielded than the other fragments in the alcohol termination.The intense mode at 40 ppm is related to the shielding of C nuclei adjacent to the carbonyl groups, while the two peaks around 47 and 57 ppm belongs to nuclei between the carbonyl groups and the ethyl and methyl ramifications in the alcohol and carbonyl terminations, respectively.In the region between 65 and 90 ppm are five peaks related to the main chain C nuclei, where the modes are, respectively, associated to the C nucleus bonded to the methyl group in the carbonyl termination (69 ppm), followed by the slightly less shielded nuclei bonded to other methyl ramifications (72 ppm), the nucleus bonded to the ethyl and the alcohol termination (74 ppm), the slightly poorly shielded C bonded to the ethyl groups (76 ppm), and the nucleus bonded to the ethyl ramification next to the carbonyl termination (77 pm).In the poorly shielded region, the mode at 166 ppm is related to the carbonyl group C between the ramifications in the carbonyl termination, followed by the peak at 168 ppm from other carbonyl C nuclei, and the last mode around 195 ppm associated with the carbonyl termination nucleus that is very poorly shielded in comparison to the alcohol termination C at 74 ppm.The same tendency is seen for the 1 H NMR spectrum, in which H nuclei in ramifications groups are more shielded and show chemical shifts <3 ppm, while the chain H atoms are less shielded and have shifts >3 ppm.On the other hand, the PCL polymer shows a less shielded structure than PHBV, which reflects a system more sensitive to the chemical environment due to its linear nature.As evident in the NMR spectra of PCL, due to its linear chain with no ramifications and symmetrical nuclei it can be easily distinguished from PHBV.Instead, the PCL structure has a mode around 25 ppm related to the carbonyl termination, followed by the region of 26-40 ppm with main chain C modes, which are less shielded and show higher chemical shifts as the nuclei are closer to the carbonyl groups.The peak at 50 ppm belongs to a nucleus adjacent to the carbonyl termination, while the C adjacent to the alcohol termination is found at 66 pm, and the modes at 69 ppm are associated with nuclei between C-C-O in the frontier between fragments.Like PHBV, the carbonyl nuclei are less shielded and related to three modes around 168, 169, and 170 ppm for each of the frontier regions between fragments, where the nuclei are slightly poorly shielded as they get closer to the alcohol termination, and there is the last peak around 195 ppm from the carbonyl termination.PCL shows a cleaner 1 H spectrum due to its symmetrical nature, having three pairs of H nuclei in the same axis of the frontier O nuclei, and two pairs in the opposite direction of the axis.Hence, the first mode at 0.3 ppm belongs to the alcohol H, while the peaks around 1.1, 1.7, and 2.2 are associated with nuclei in the middle of the three pairs, the opposite H pairs, and the pairs next to the carbonyl group, respectively, with the carbonyl termination lightly shifted.The peaks around 3.6 and 3.9 ppm are related to the pair next to the alcohol termination and the pairs adjacent to the frontier O nuclei, respectively.Lastly, the mode at 10.3 belongs to the poorly shielded nucleus from the carbonyl termination.
As expected, the blends show a hybrid behavior between PHBV and PCL, where the PHBV fragments have heavily shielded ramifications and decrease next to the carbonyl groups, while the PCL fragments are more shielded in the middle regions than the extremities in the fragments' frontier zones.These modes, discussed for each part of the PHBV and PCL polymers, are found in the blends' pattern, with slight shifts depending on where the PHBV-PCL frontier is present.In general, the exchange of a PCL fragment for PHBV in the polymer chain has a very low impact on the structure in terms of shielding, which, except for the BC blend and its mostly equivalent nuclei, shows a seamless difference with a chemical shift around 2 ppm for carbonyl C nuclei in the PHBV fragment.However, this small difference makes the PHBV frontier regions slightly more sensitive to the chemical environment than PCL zones.In this framework, the same tendencies observed in the 1 H NMR spectra of both polymers and blends can also be identified in the experimental work of Xia et al. [93], in which nuclei in the carbonyl groups show higher chemical shifts and are significantly less shielded than other regions.In addition, as clearly shown in the higher regions of the NMR spectra, nuclei in the termination parts are the least shielded for all structures, thus making the terminations the most sensitive regions with respect to the chemical environment.
A UV-vis spectrum for each polymer is found in Figure 5A.All systems have two bands, one with a higher intensity around 215 nm and the other with a lower intensity near 288 nm.The higher peak of PCL is slightly blue shifted to around 211 nm due to transitions within 195 nm.Hence, the blends have the same blue-shifted band between the PCL and PHBV positions but with lower intensities, except for the BC blend, which shows higher intensity than PCL.The excitations involved in the UV-vis peaks are characterized by electronic transitions of the following nature: HOMO-3 → LUMO+2 (5.79 eV) for PHBV, HOMO-1 → LUMO+2 (6.35 eV) for PCL, HOMO-5 → LUMO+1 (5.77 eV), HOMO-4 → LUMO+1 (5.79 eV), HOMO-2 → LUMO+5 (5.78 eV), and HOMO-1 → LUMO+1 (5.75 eV) for the BA, BB, BC, and BD blends, respectively.In the same perspective, all excitations at 288 nm are due to HOMO → LUMO (4.31 eV) electronic transitions.Our results are aligned with the experimental literature, as seen in the work of Abdelrazek et al. [94] regarding PCL/PMMA blends, where they reported a characteristic absorption band centered around 280 nm related to carbonyl groups for pure PCL.Thus, the band described by the authors is likely caused by HOMO → LUMO transitions, while the HOMO and LUMO are indeed localized within carbonyl groups in the termination regions as shown in the MOs maps in the Supplementary Information.In addition, Nanda et al. (2011) in their work with PHBV/PLA blends also reported an absorption band around 350 nm for pure PHBV polymer, likely being related to HOMO → LUMO transitions as well [95].Figure 5B-G shows the TDOS of the polymers.As is shown in the TDOS spectra, the HOMO and LUMO have relatively low densities of states in comparison to the valence and conduction bands, which show two bands each with very high contributions, from −14 to −8 eV for the valence band, and from 0 to 5 eV for the conduction band.In general, both HOMO and LUMO are localized within the carbonyl terminations in all polymers, which HOMOs are mostly composed of p atomic orbitals from carbonyl O (~67%), as well as with a small contribution of p and s atomic orbitals from neighbor C (~18%) and H (~11%) nuclei, respectively.On the other hand, LUMOs show a major contribution of p atomic orbitals from carbonyl C (~66%) and O (~28%) nuclei, and lesser contribution of s atomic orbitals from H (~5%), thus in agreement with the MOs maps shown in the Supplementary Information.
Although the PHBV has a very slight band gap energy difference to PCL (0.04 eV), and the blends have a tendency of intermediate band gap between the pure polymers, the energy difference is still significantly small and within the error range of DFT approximations, in which an indication of which system has a lower or higher band gap would be inaccurate.Overall, information regarding the polymers and blend properties, such as the band gap, is very scarce, due to being very dependent on variables like molar weight.However, Abdelrazek et al. (2015) studied the spectroscopic properties of PCL polymer and derivative blends, where, in their work, a direct and indirect band gap of 4.71 eV and 4.26 eV were estimated using the Devis and Mott model, while the Tauc method reached 3.83 eV [94].As mentioned beforehand, this high shift in the band gap is caused by the huge difference in molar weight between a PCL polymer and our model, in which the polymerization and molar weight leads to a heavy increase in intermediate states between the valence and conduction bands and the band gap decrease.In addition, despite PCL and PHBV having similar band gap values in our models, in reality the band gaps should show a significant discrepancy between these polymers because of the difference between the linear PCL and the ramified PHBV structures, having a strong influence in the polymer packing and density.
As a complement to this discussion, Table 2 lists some electronic properties derived from the HOMO and LUMO energies.As shown in the table, all values are mostly identical and with very slight shifts, such as the BA blend with properties very slightly inclined to the PHBV polymer while the BB, BC, and BD blends are more PCL alike, although the DFT error must be taken in consideration.Nonetheless, all systems show an average band gap of 6.2 eV, hardness of 3.1 eV, softness of 0.32 eV, and Muliken Electronegativity of −4.3 eV, while the electrophilicity is more sensitive to the energy values and ranges from 34.93 eV for the BC blend to 39.30 eV for BA. Figure 6A-F shows the NCI maps and their respective reduced density gradient (RDG) maps, where the blue, green, and red regions represent strong interactions (i.e., hydrogen and halogen bonds), van der Waals (vdW) interactions, and strong repulsions (i.e., steric effects), respectively.NCI analysis helps in understanding the nature and strength of intermolecular interactions in complex systems [96].In this framework, the PHBV structure (Figure 6A) has strong steric effects between the carbonyl oxygens, and weak vdW interactions between the same oxygens and the adjacent hydrogen atoms from the methyl and ethyl ramifications, contributing to a non-linear structure.On the other hand, those same interactions are of significantly lower intensity for the PCL polymer (Figure 6B), showing a structure with strong repulsions between the carbonyl oxygens in a smaller area and vdW interactions between the chain hydrogens, complementing the linear geometry.Hence, both patterns are also observed in the blends' structures, where the PCL monomer units show smaller surfaces of steric effects and vdW interactions, while the PHBV fractions have wider surfaces due to interactions with the methyl and ethyl ramifications that strongly contribute to the structural distortion and break of linearity in the PHBV/PCL boundary, as seen clearly in the BB blend (Figure 6D).ELF surface maps shown in Figure 7 were made from half planes cut from the Xand Y-axis with the intersection at the Z-axis origin.ELF surface maps indicate electron localization, where red values are associated with higher localization.Hence, the ELF analysis provides insights into the bonding and delocalization of electrons in a molecular system [97].Although the electron localization density outside of the intersecting plane is without the ELF surface range due to the nonlinear nature of the PHBV polymer (Figure 7A), there is very high electron localization around the hydrogen nuclei due to lone pairs and high localization within the oxygen nuclei valence shells and the C-C and C-O bond midpoints, while an electron-gas-like pair probability from external shells surrounds the polymer chain.In contrast, the electron localization is more clearly shown in the PCL polymer (Figure 7B) because of its linear nature, in which very high localization is seen around the terminal hydrogen nuclei due to lone pairs and high localization around the polymer chain bonding midpoints and the oxygen nuclei external shells, with wider areas on double-bonding oxygens.In the case of PCL, the localized electrons around the chain hydrogen nuclei are not as visible as in the PHBV polymer because of their angle in relation to the surface.The same patterns are observed within the blends' electron localization surfaces (Figure 7C   Both NCI and ELF analyses, in general, contribute to our understanding of molecular structure and interactions of blends, that is, shedding light on the fundamental physical properties and behavior of such compounds [96][97][98].

Conclusions
In this study, we tried to evaluate the physical and mechanical properties of PHBV by simultaneous or isolated modifications based on low-molar-mass PCL at γ-irradiation doses of 50 and 100 kGy.In addition, we have used computational simulation to study the interactions between PHBV and PCL at the molecular level, which in turn may provide relevant chemical insights into their spectroscopic and physicochemical behavior.It was concluded that the processing conducted for all materials, i.e., extrusion, drying, and injection, was adequate for providing uniform specimens.The 50 and 100 kGy irradiation doses degraded (via splitting of macromolecules) the pure polymer and the blend, as indicated by the mechanical properties (more evident from the deformation owing to the tensile strain at break).The radiation caused splitting instead of crosslinks in the polymer chains.In short, studying the effects of γ-radiation on biodegradable polymers can also help to improve our understanding of the degradation mechanisms of these PHBV/PCL materials, which can be useful in developing strategies for improving their stability and performance in the future.
-F) as well due to their PCL and PHBV monomer units, in which highly localized lone pairs are found on hydrogen nuclei and high localized electrons are around the oxygen nuclei valence shells and the C-C and C-O bonds midpoints.Nanomanufacturing 2024, 4, FOR PEER REVIEW 7

Table 2 .
Electronic properties derived from the HOMO and LUMO energies.