• Open Access
  • Editors Choice
B: Soft Matter, Fluid Interfaces, Colloids, Polymers, and Glassy Materials

Atomistic Insights into Structure and Properties of ε-Caprolactone Oligomers
Click to copy article linkArticle link copied!

  • Mai Ahmed
    Mai Ahmed
    Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, Germany
    Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
    More by Mai Ahmed
  • Deniz Yilmaz
    Deniz Yilmaz
    Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, Germany
    Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
    More by Deniz Yilmaz
  • Purushottam Poudel
    Purushottam Poudel
    HIPOLE Jena (Helmholtz Institute for Polymers in Energy Applications Jena), Lessingstrasse 12-14, 07743 Jena, Germany
  • Felix H. Schacher
    Felix H. Schacher
    Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, Germany
    Jena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
    HIPOLE Jena (Helmholtz Institute for Polymers in Energy Applications Jena), Lessingstrasse 12-14, 07743 Jena, Germany
  • Eva Perlt*
    Eva Perlt
    Otto Schott Institute of Materials Research (OSIM), Faculty of Physics and Astronomy, Friedrich Schiller University Jena, Löbdergraben 32, 07743 Jena, Germany
    *Email: [email protected]
    More by Eva Perlt
Open PDFSupporting Information (1)

The Journal of Physical Chemistry B

Cite this: J. Phys. Chem. B 2026, 130, 5, 1675–1683
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcb.5c06385
Published January 27, 2026

Copyright © 2026 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Abstract

Click to copy section linkSection link copied!

The design of functional and sustainable materials requires a detailed understanding of the material properties and degradation mechanisms. In particular, the design of fully biodegradable polymers could allow a quick and controlled decomposition of materials before they accumulate in the environment and break down to micro- and nanoplastics. An important degradation pathway proceeds via the hydrolysis of polyesters. To obtain the best performing material candidates, a multiscale-level understanding that takes into account electronic structure combined with multiple configurations at the macroscopic scale is necessary. In this contribution, we present the extension of the multiscale Quantum Cluster Equilibrium method to oligomer materials. We showcase the first application of this methodology to oligomer systems, in particular oligo(ε-Caprolactone). The ε-Caprolactone oligomers were synthesized and characterized comprehensively by means of NMR, SEC, DSC, and TGA. Experimentally, two melting temperatures were observed, which were predicted by theoretical calculations and are in convincing agreement.

This publication is licensed under

CC-BY 4.0 .
  • cc licence
  • by licence
Copyright © 2026 The Authors. Published by American Chemical Society

Special Issue

Published as part of The Journal of Physical Chemistry B special issue “Physical Chemistry of Microplastics and Nanoplastics”.

Introduction

Click to copy section linkSection link copied!

Polymers have been widely incorporated into modern material science, with applications ranging from packaging (1,2) and textiles (3) to biomedical devices, (4,5) electronics, (6,7) and sustainable technologies. (8) The abundance of these stable and long-lasting materials combined with poor waste management has caused the problem of micro- and nanoplastics, which accumulate mainly in marine environment. (9) Promising strategies to prevent the formation of micro- and nanoplastics include the prevention of the degradation into micro- and nanoparticles (9) or the further degradation of microplastics by microorganisms. (10) Alternatively, the design of fully biodegradable polymers would allow the implementation of proper waste treatment and reduce the amount of plastic waste that accumulates in the environment before being broken down into micro- and nanoparticles, which are impossible to eliminate.
Polyesters can undergo hydrolysis as an important step in degradation, which makes them a prospective class of compounds for biodegradable polymers. Their accessible functionalization allows one to tailor their molecular structure to achieve specific mechanical, physical, and chemical properties and furthermore improve hydrophilic properties, thereby enhancing degradability. (11−14) Due to the ester linkages in these compounds, they are prone to hydrolytic degradation, i.e., the cleavage of ester bonds by water molecules, decreasing the polymer’s molecular weight. Functionalization by incorporation of hydroxyl (−OH) or carboxyl (−COOH) in the polymer backbone enhances hydrophilicity and water absorption, making the polymer more accessible for further break down. (15−18) Tailoring these properties, for example of polycaprolactone, could make these polymers suitable compounds for applications such as drug delivery systems. (19) To obtain a mechanistic understanding of these properties and processes and finally design functional degradable polymers, a theoretical model that includes quantum chemical data is essential.
Polymer systems inherently exhibit a multiscale character, and there are relevant processes happening across all scales from the atomistic monomer level up to macroscopic engineering scale. In a perspective, Schmid introduced the different theoretical techniques which address different types of polymers and, more importantly, their different properties and phenomena, which occur on different time and length scales. (20) The most commonly used polymer theories include quantum chemical and atomistic classical simulations at the monomer and oligomer scale, generic statistical mechanics models including lattice models and off-lattice models at the polymer scale, ultracoarse-grained particle-based models at the “blob” scale of interacting polymers, field theories at the mesoscopic scale, and finally continuous field theories and transport equations at the engineering scale. However, the different time and length scales strongly overlap in polymer systems, requiring multiscale methodologies. To overcome this limitation, various flavors of static and dynamic coarse graining techniques are available. Finally, machine learning-based approaches are gaining more and more attention on different scales as well. For example, machine-learned force fields, coarse-grained potentials or the direct prediction of polymer properties based on machine-learned structure–property relationships or using other predictors were already presented, see (20) and references therein.
Still, even with the wealth of modern theoretical tools, there remain open challenges. Two of these are strong inhomogeneities, which require a certain variety or flexibility in coarse-grained potentials and are therefore hard to cover. Similarly, defects present a type of very dilute inhomogeneity. This low concentration, or rareness, poses another challenge to the simulation, namely, the need for even larger system sizes.
Another example of a scale-bridging approach is the Quantum Cluster Equilibrium (QCE) theory, (21−23) which has been successfully applied to a number of pure associated liquids (21,24) as well as mixtures thereof. (25) It has already been applied to describe the freezing of water, so it is generally applicable to solid materials. (26) It combines a highly accurate electronic structure treatment at the cluster level with an enhanced statistical thermodynamic weighting of these clusters to obtain the macroscopic partition function and, thereby, virtually all thermodynamic potentials. Additionally, the populations of the individual clusters in combination with their quantum chemical characterization allow to compute liquid phase spectra, or even acidic constants or the ionic product of water. (27,28) If transferred to polymer materials, this type of approach might be able to address the issues raised above: inhomogeneities may be covered by including different representative oligomer fragments. The same holds for defects, which may be introduced at the cluster level and then may be present in only small amounts, as was shown at the example of dissociated water clusters in liquid water to determine the ionic product. (27)
In this combined experimental and theoretical study, we extend the multiscale QCE theory to polymers, as they may be found in micro- and nanoplastics. We targeted oligo-(ε-Caprolactone) (OεCL) with a degree of polarization (DP) of 10 for synthesis and characterization. In parallel, a selection of decamer clusters of the OεCL was characterized using density functional theory and subsequently processed with the QCE approach. We will demonstrate the general feasibility of QCE for polymer systems and present first estimates of a phase transitions in convincing agreement with experimental observations.

Experimental Methods

Click to copy section linkSection link copied!

Materials

Tin(II) 2-ethylhexanoate (∼92.5–100%) (SnOct2) and benzyl alcohol (BnOH) were acquired from Sigma-Aldrich. Hydrochloric acid (HCl) (∼37%) and ε-caprolactone were obtained from ThermoFischer Scientific. Calcium hydride (CaH2) was purchased from J&K Scientific. All chemicals, excluding SnOct2 and ε-Caprolactone were used as received if not stated otherwise.
ε-Caprolactone was dried over CaH2 at room temperature (RT) for 24 h, distilled under reduced pressure and SnOct2 was dried via azeotropic distillation of toluene (3 times) before oligomer synthesis. The solvents (methanol (MeOH), tetrahydrofuran (THF) (not stabilized, HiPerSolv CHROMANORM for HPLC) and toluene) were purchased from VWR Chemicals. Toluene was dried with molecular sieves before use. The deuterated solvent, deuterated chloroform (CDCl3), was obtained from Deutero. PTFE syringe filters were acquired from Carl Roth under Membrane Solutions.

Instrumentation

Proton Nuclear Magnetic Resonance (1H NMR) Spectroscopy measurements were performed using a 300 MHz Bruker AC spectrometer with deuterated chloroform (CDCl3) as the solvent. The recorded spectra were adjusted in Mestrenova, referenced to the residual solvent peak, and underwent automatic baseline and phase corrections.
Size Exclusion Chromatography (SEC) measurements were done on a Schimadzu 10er Series system with a CBM-20A controller, a DGU-14A degasser, a LC-10AD vp pump, a SIL-10AD vp autosampler, a CTO-10A vp oven, RID-10A, SPD-10AD VP, ETA- 2010, SLD 7100 detectors, and a PSS DV guard/Linear M (5 μm particle size) column. THF was utilized as the eluent with an elution rate of 1 mL min–1 at a temperature of 30 °C and the calibration was performed by polystyrene (PS) standards.
Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TGA) measurements were carried out on a NETZSCH DSC 204F1 Phoenix differential scanning calorimeter and on a NETZSCH TG 209F1 Libra thermobalance. TGA measurements were done under a nitrogen atmosphere (N2) in a temperature range from 25 to 600 °C with a heating rate of 20 K min–1. DSC measurements were performed under a nitrogen atmosphere (N2) (20 mL min–1) in a temperature range from −110 to 150 °C with heating and cooling rates of 10 K min–1. Second heating and cooling scans were used for evaluation.
Small angle X-ray scattering (SAXS) measurements were conducted on an Anton Paar SAXSpoint 5.0 SAXS system (Anton Paar, Graz, Austria) equipped with a Primux 100 microfocus X-ray source (Cu Kα radiation; λ = 1.54Å), ASTIX 2D multilayer X-ray optics, and a 2D EIGER2 R 1 M hybrid photon-counting detector shielded with a Mylar film (Dectris, Baden, Switzerland). The oligomer sample was put in a powder cell containing two polyimide windows and packed densely. The sample-source distance was adjusted to 1621 nm for SAXS measurements and 64.2 nm for wide-angle X-ray scattering (WAXS) measurements. The oligomer was heated up to 50 °C prior to its measurements to erase thermal history and cooled down to 25 °C at which the measurements were carried out. The data were processed with corrections for sample transmission, background scattering, and detector sensitivity. The analyses of the obtained data from SAXS and WAXS measurements were carried out in OriginPro 2022 through the peak analysis tool. The baseline of the WAXS pattern was manually fitted, and the peaks were manually selected. The determination of the degree of crystallinity (Xc) involved the utilization of calculation of the areas under the peaks and under the entire curve. The determination of the crystallite sizes (Dhkl) included the fitting of the most prominent diffraction peaks with a basic Gaussian function for the calculation of the full width at half-maximum.

Syntheses of the εCL Oligomer (OεCL)

The synthesis of the oligomer of εCL with a desired degree of polymerization of 10 followed a procedure adapted from the methods given by Ren et al. (29) and El Habnouni et al. (30) The synthesis revolved around the subsequent addition of εCL (5 mmol, 570 mg, 554 μL), anhydrous toluene (2 mL), SnOct2 (0.05 mmol, 20.2 mg, 16 μL) and BnOH (0.5 mmol, 54 mg, 52 μL) to a dry Schlenk flask under constant argon flow. The flask was then sealed and stirred in an oil bath at 120 °C for 1 h. The synthesis was then terminated with an excess of 1N HCl and the oligomer was recovered through its precipitation in cold methanol and was dried under vacuum. 1H NMR (300 MHz, CDCl3, δ ppm): 7.40–7.30 ppm (BnOH, aromatic, 5H), 5.12 ppm (BnOH, ArCH2O, 2H), 4.12–4.02 ppm (−CO(CH2)4CH2O −, 2H), 3.70–3.63 ppm (−CH2CH2OH, 2H), 2.40–2.20 ppm (−COCH2(CH2)3CH2O −, 2H), 1.76–1.20 ppm (−COCH2CH2CH2CH2CH2O −, 6H).

Computational Details

Click to copy section linkSection link copied!

For a computational investigation of oligo(ε-Caprolactone) with DP = 10, the corresponding decamer units were generated and subsequently subject to a Quantum Cluster Equilibrium calculation.
Within the QCE method, a macroscopic phase is represented by a set of finite cluster structures which are assumed to be in thermodynamic equilibrium. (21−23) Compared with routine thermochemistry calculations, it includes the following corrections and improvements. The cluster equilibrium results in a temperature-dependent distribution of different cluster motifs and therefore allows one to include different configurations, which are balanced regarding enthalpic and entropic effects. The model contains two correction parameters. The mean-field parameter amf scales the mean-field interaction potential, which covers the effective intercluster interaction (the intracluster interaction energy is treated by the quantum chemical calculation). Furthermore, the volume available for translation is corrected for the clusters’ volumes and scaled by the exclusion volume parameter bxv. Both parameters are optimized to reproduce reference data, which can be the density at one or more temperatures or the phase transition temperature.
To set up the cluster set, one cluster was designed to be representative of a crystalline domain by arranging the decamer chain in four parallel subunits, as shown in Figure 1 (left).

Figure 1

Figure 1. Example structures, as used in the QCE study. Crystalline motif (OεCL-c, left) as well as one exemplary amorphous structure (OεCL-a3, right).

The remaining clusters in the set represent amorphous structures in random orientations and were generated using openbabel. (31) The conformer search tool using the genetic algorithm was applied with the default settings. The RMSD was chosen as the score property, and 25 conformers were written. All clusters were characterized computationally using the Turbomole V. 7.8.1 program suite. (32,33) The efficient composite functional PBEh-3c was chosen. (34) It is based on the hybrid PBE0 functional (35) and used with the def2-mSVP basis set, the D3(BJ) dispersion correction (36) with Becke–Johnson damping (37) and the inclusion of three-body terms as well as geometrical counterpoise correction. The SCF convergence criterion was chosen to be 10–7.
A geometry optimization was performed with a convergence criterion of 1 × 10–4 for the norm of the Cartesian gradient. The derivatives of quadrature weights were included in the gradient evaluations using the option weight derivatives. Harmonic vibrational analysis was performed using Turbomole’s aoforce module to confirm that the geometry is a true minimum on the potential energy surface. Furthermore, vibrational frequencies are processed in the vibrational partition function.
Apart from the crystalline structure, seven out of the 25 amorphous structures turned out to be true minimum structures, which was considered sufficient for subsequent QCE calculations. Figure 1 (right) shows one of the amorphous structures. All snapshots of the clusters are shown in Figures S1–S8 in the Supporting Information.
Subsequently, QCE calculations were performed using Peacemaker V. 4.0.0. (23) The cluster set consisted of seven amorphous clusters, labeled OεCL-a1 to OεCL-a7, and the crystalline motif OεCL–c. Calculations were performed at standard pressure and over the temperature interval from 200 to 500 K with increments of 0.5 K. QCE calculations require electronic energies to be given as relative energies with respect to a reference cluster, which commonly is the isolated monomer in studies of associated liquids. In the present case, all clusters have the same size; therefore, the least stable structure OεCL-a6 was chosen as the reference cluster. Thus, the interaction energies were calculated as the energy difference with respect to this reference cluster in kJ mol–1; see Table 1.
Table 1. Cluster Labels and Interaction Energies in kJ mol–1 with Respect to the Reference Cluster OεCL-a6
cluster labelΔE [kJ mol–1]
OεCL-c–106.01
OεCL-a1–18.29
OεCL-a2–83.45
OεCL-a3–60.23
OεCL-a4–27.30
OεCL-a5–67.43
OεCL-a60.00
OεCL-a7–61.54
The model parameters were sampled in the range from 1.0 to 50.0 Jm3 mol–2 for amf and from 0.5 to 5.0 for bxv. Both intervals were chosen to be significantly larger than those for associated liquids in previous applications to allow more flexibility for this unknown system and a much larger reference species. Both parameters were optimized in four subsequent cycles, each refining the optimization interval, using the keyword grid iterations in Peacemaker. As reference data for the optimization, only the density of 1.145 g cm–3 at 298.15 K was chosen.

Results and Discussion

Click to copy section linkSection link copied!

Oligomer Synthesis and Characterization

The synthesis of the εCL oligomer with an intended degree of polymerization (DP) of 10 was centered on a coordination–insertion ring-opening mechanism of εCL. The mentioned oligomer formation, initiated by BnOH and catalyzed by SnOct2, is shown in Scheme 1.

Scheme 1

Scheme 1. Oligomer Formation of εCL with BnOH/SnOct2 as the Initiator/Catalyst
The characterization of the oligomerization product of εCL with the mentioned system (OεCL1) was performed through proton NMR spectroscopy and SEC measurements. The 1H NMR spectrum of the acquired product is shown in Figure S9 in the Supporting Information, and its SEC chromatogram is given in the top panel of Figure 2.

Figure 2

Figure 2. Characterization of the OεCL1 oligomer. Top, SEC chromatogram; middle, DSC heating and cooling scans; bottom, TGA curve.

The 1H NMR spectrum of the product confirmed the successful oligomerization of εCL, as evidenced by the desired signals of BnOH protons and those of the εCL backbone. The obtained spectrum was used to determine the DP and the Mn¯ of the oligomer. The number of εCL repeating units was acquired by comparing the integration ratios between the methylene group protons of εCL (4.12–4.02 ppm) and the aromatic protons of BnOH (7.40–7.30 ppm).
The SEC revealed a relatively narrow and monomodal molecular weight distribution. The Mn¯ value of the oligomer, calculated from the number of εCL repeating units obtained from the 1H NMR spectrum and its corresponding molecular weight, was also determined using SEC measurements.
The Mn¯ value obtained from the proton NMR spectrum was in accordance with the feed in terms of the degree of polymerization, whereas the one acquired from the SEC measurement was noted to be higher. The relatively low dispersity of OεCL1 and the higher Mn¯ value of OεCL1 obtained from SEC were indicative of good control over the reaction. Table 2 contains the repeating units of εCL by means of the composition of the oligomer, along with their Mn¯ values and dispersity.
Table 2. εCL Amount in Feed Respecting BnOH, Composition, Values, and Dispersity of the Oligomer
oligomerεCL/BnOH (mmol/mmol)compositionaa,b (kg mol–1)c (kg mol–1)dispersity Đ
OεCL15/0.5O(εCL12)1.52.71.24
a

Determined by 1H NMR spectroscopy (300 MHz, CDCl3).

b

Determined by .

c

Determined by SEC measurements (THF, PS calibration).

The εCL oligomer OεCL1, with a DP of 12, determined via 1H NMR spectroscopy, additionally underwent thermal and morphological characterization. The thermal characterization was performed through DSC and TGA measurements. The DSC heating and cooling scans of OεCL1 are depicted in Figure 2 (middle panel).
PεCL is a semicrystalline polyester having a melting temperature (Tm) ranging from 50 to 70 °C and a glass transition temperature (Tg) of around −60 °C. (38) The aforedepicted heating scan of OεCL1 possessed split melting peaks, one at around 42 °C and the other at around 46 °C. This splitting could be attributed to the melting of the crystalline phases of PεCL, similarly observed by Ju et al. (39) for a PεCL copolymer with pendant pyridyl disulfide groups, and to reorganizations as mentioned by Fernández-Tena et al. (38) for PεCL with low Mn¯ values. The relatively low Tm values of OεCL1 in comparison to those of PεCL, on the other hand, were explained by the low DP and thus the low molecular weight of the oligomer. The heating scan additionally showed a Tg of around −65 °C for OεCL1, relatively similar to the one reported for PεCL. (38)
The cooling scan of OεCL1, on the other hand, displayed a crystallization temperature (Tc) of around 16 °C. This value was underlined to be relatively close to the ones of PεCL samples with Mn¯ ranging from 4 kg mol–1 to 24 kg mol–1 given by Klonos et al. (40) Therein, the Tc values of the PεCLs were between 20 and 33 °C, respectively, with consideration of the mentioned Mn¯. It was thus plausible for OεCL1 to possess a Tc of around 16 °C, given its lower Mn¯ of 1.5 kg mol–1 by 1H NMR spectroscopy and of 2.7 kg mol–1 by SEC chromatography.
The acquired εCL oligomer, OεCL1, underwent TGA measurements in addition to DSC measurements, as previously mentioned. The TGA curve of OεCL1, which possesses crucial information on its stability and degradation behavior, is shown in Figure 2 (bottom panel).
The depicted TGA curve of OεCL1 is indicative of its stability up to around 205 °C. The oligomer then experienced a main degradation until around 343 °C, with almost complete weight loss. A similar stability behavior was reported by Persenaire et al. (41) for a PεCL of Mn¯ of 1.8 kg mol–1, with an onset of degradation at 230 °C under helium atmosphere, extending to 600 °C. The observed degradation of OεCL1 might be attributed to a possible combination of chain cleavage by pyrolysis and unzipping depolymerization.
The morphological characterization of OεCL1 was performed through SAXS and WAXS measurements. The 1D integrated SAXS pattern of the oligomer and the corresponding WAXS data are shown in Figure 3.

Figure 3

Figure 3. 1D integrated (a) SAXS patterns of OεCL1 corresponding to q and (b) WAXS patterns of OεCL1 corresponding to 2θ.

The SAXS pattern of the oligomer shows two prominent reflexes at q1 ≈ 0.57 nm–1 and q2 ≈ 1.17 nm–1. The reflex at a q value between 0.4 and 0.6 nm–1 was attributed to long periods of crystalline lamellar structures, similarly explaining the reflex at a q value of 0.04 Å–1 in the SAXS patterns of PεCL and its Pluronics in the works of Tenorio-Alfonso et al. (42) The near-integer ratio (q2/q1 ≈ 2) is presumably attributed to the second-order harmonics of lamellar stacking, suggesting a periodic arrangement of crystalline and amorphous layers rather than microphase separation. The first reflex corresponding to the long periods (q1) was later utilized to determine the d-spacing of the domains using the following equation (43)
q=2πd
(1)
The morphology of OεCL1 was further investigated through WAXS measurements. The WAXS pattern of the oligomer (Figure 3 b) exhibited the expected diffraction peaks, consistent with those reported for PεCL by Nanaki et al. (44) Here, the most prominent diffraction from (110) and (200) planes appeared at angles of 2θ ≈ 22.4° and 2θ ≈ 24.8°. The two crystallographic planes were utilized to determine the crystallite sizes (Dhkl) in the direction perpendicular to the (hkl) plane through the Scherrer equation (45)
Dhkl=KλBhklcosθ
(2)
where K is the shape factor (typically 0.9), λ is the X-ray wavelength (0.154 nm), Bhkl is the full width at half-maximum (FWHM) of the diffraction peak in radians after instrumental correction, and θ is the Bragg angle.
The WAXS pattern of the oligomer was also utilized to determine its degree of crystallinity (Xc) by relying on the diffractions from all of the planes through the equation given below (46)
Xc=AcAc+Aa
(3)
where Ac represents the integrated area of the crystalline peaks and Aa represents the integrated area of the amorphous region in the WAXS pattern.
The d-spacing of the domains (d), the crystallite sizes (Dhkl), and the degree of crystallinity (Xc) of the acquired oligomer, OεCL1, are given in Table 3.
Table 3. Crystalline Domain Spacing or d-Spacing (d), Crystallite Sizes (Dhkl), Degree of Crystallinity (Xc) of OεCL1
 Dhkl (nm) 
d (nm)D100D200Xc (%)
10.987.95 ± 0.837.20 ± 0.7668.90
Crystallite sizes in the range of a few nanometers were recalled to be typical for small oligomers and poorly ordered polymers, in accordance with the values obtained for OεCL1. The aforementioned similar Dhkl values obtained for the two planes (110) and (200), on the other hand, indicated nearly isotropic crystallite dimensions. The combined results of the SAXS and WAXS measurements supported a lamellar morphology for OεCL1. The WAXS-derived crystallite sizes (7.20–7.95 nm) were smaller than the SAXS-derived long period (10.98 nm), indicating the presence of intervening amorphous layers of approximately 3–4 nm. Such lamellar structures were known to be consistent with the behavior of short-chain PεCL oligomers, which have the ability to form ordered crystalline domains separated by amorphous regions. (47)

QCE Calculations

This first application of QCE to oligomer systems was completed successfully and smoothly, demonstrating that QCE is capable of being extended to polymer systems. The optimized parameters were determined to be amf = 6.230 Jm3 mol–2 and bxv = 0.890. An excellent agreement with the experimental reference was achieved, with a calculated density of 1.145 g cm–3 and an error of 0.2989 × 10–13. It is noted that including more reference data for the density would certainly improve predictions for thermodynamic quantities, especially their temperature-dependent behavior. However, the effect of both parameters on the populations is expected to be only minor, and the one reference density chosen in this study is considered sufficient for the following analysis.
The QCE parameters are found to be somewhat different from those of associated molecular liquids, which have been extensively studied using QCE in the past. The amf is found to be significantly larger for the OεCL compared to values for liquids, such as water, which are typically smaller than one. The reason for this is probably the large cluster size in this study. The mean-field interaction energy depends on the density and is calculated in the following way for each cluster
Eint=amfiV
(4)
where i denotes the number of monomers in the respective cluster and V is the phase volume. However, in our example, the number of repeat units in each cluster is ten, and at the same time, the number of monomers i entered in the QCE input is one (because a smaller number of monomers generally leads to a more stable calculation). This discrepancy has no consequences for the results, and only a larger amf is observed to compensate for the mean-field interaction.
Also, the parameter bxv differs from earlier studies, in which it was typically slightly larger than one. This is likely due to the close packing of crystalline motifs combined with the interpenetration of fragments of the amorphous structural motifs.
The populations of clusters with a significant population above 1% over the investigated temperature range are shown in Figure 4. Interestingly, only two clusters are populated. At lower temperatures, the crystalline cluster of OεCL-c is the only structure with a population of 1; see the blue curve in Figure 4. At ∼250 K, the population of this species decreases, while the cluster OεCL-a3, represented by the orange curve, becomes increasingly populated. Experimentally, at 25 °C, a degree of crystallinity of 68.90% was observed from the WAXS data; see Table 3. At this temperature, the QCE population of the crystalline motif is 83.44%, which is an excellent agreement, especially given the simplicity of the model, the limited cluster set, and the small number of reference data. Both curves intersect at 326.75 K (54.6 °C), and both clusters have a population of 50%, as indicated by the gray lines in Figure 4. Table 3 reports crystallite sizes of 7.95 (83) and 7.20 (76) nm. We emphasize that the cluster OεCL-c is characterized in vacuum, i.e., in the absence of periodic boundary conditions. However, constructing a hexagonal unit cell that can accommodate the structure results in cell dimensions of a ≈ 10Å, b ≈ 11.5Å, and c ≈ 23Å, which would indicate that crystallites may be represented as a 7 × 7 supercell of the crystalline motif.

Figure 4

Figure 4. Temperature-dependent cluster populations showing only significantly populated clusters (>1%). The gray line indicates the interception of both curves.

The crystalline structure is the most stable one, as shown in Table 1, and consequently dominates the entire phase at low temperatures, in agreement with experimental observations. The great potential of the QCE approach for such systems becomes apparent at larger temperatures: by taking entropic effects into account, the enthalpic preference for the crystalline motif is counterbalanced, which leads to an increasing population of the amorphous OεCL-a3 cluster, which is less stable than OεCL-c by as much as ∼45 kJ mol–1. Importantly, this change in the population of cluster species is an indicator of a phase transition. However, it is interesting to note that out of seven amorphous cluster structures, a cluster that falls within the medium range of energies is exclusively populated. Even enthalpically much more favorable clusters like OεCL-a2 do not contribute at all. This again underlines the importance of vibrational effects. Furthermore, for this system consisting of rather long chains, there is likely a substantial influence of the overall shape of the cluster (coiled versus stretched), which is incorporated into the QCE method via the rotational partition function.
The calculated phase transition temperature is in astonishing agreement with the experimentally observed melting temperature (46 °C), demonstrating the general applicability of the QCE method for oligomer systems again. The deviation in the phase transition temperature of only 8 K is relatively small and can be attributed to the simplicity of the QCE method, particularly the finite cluster set.
The phase volume for 1 mole of the OεCL decamers is shown in Figure 5 for the temperature range of 285–335 K. This isobar exhibits a sharp increase in volume at 314 K (40.85 °C). This phase transition is in agreement with the experimentally observed first melting temperature. This observation is rather interesting, as the computation predicts two distinct phase transitions, which is in agreement with the two values for Tm observed in the heat flow measurements.

Figure 5

Figure 5. Phase volume of 1 mol OεCL. At 314 K, a sharp phase transition in the volume is observed.

The lower-temperature phase transition computationally results in a discontinuous increase in volume, while the computational parameters and the resulting cluster populations exhibit no such behavior. The reason for this can be explained by looking at how the volume is determined within QCE. It is calculated from the volume polynomial, which is a third-order polynomial derived from the van der Waals-like equation and thus can have up to three roots. Out of these, the volume that yields the lowest Gibbs free energy is chosen. In the present case, at 314 K, the volume solution that corresponds to the lower Gibbs Free Energy changes and results in an expansion of the phase. While the temperature is in perfect agreement with the experimental measurements, the magnitude of the expansion is overestimated. This, however, can be explained by the nature of the QCE theory, which is a fluid phase method, combined with the very limited cluster set and the lack of environmental effects at the cluster level.
In contrast, at the higher-temperature phase transition, the cluster configuration changes from the ordered crystalline type to the amorphous type. No further change in volume is observed here, but the change of the dominating cluster motifs clearly shows the breaking of the short-range structure.

Conclusions

Click to copy section linkSection link copied!

In a joint experimental and theoretical approach to well-defined OεCL decamers, we demonstrated the general applicability of the QCE method in describing oligomer systems.
The successful synthesis of the εCL oligomer with a degree of polymerization of 12 was confirmed by a combination of chromatography and spectroscopy. Thermal characterization revealed the presence of split melting peaks and a glass transition temperature similar to that of poly-(ε-caprolactone). Finally, the oligomer OεCL1 remains stable until 205 °C before it degrades in a two-step process, in agreement with previous work.
Computationally, the QCE method was successfully applied to oligomer systems for the first time. The observed population of the crystalline cluster is in surprisingly good, although not quantitative, agreement with X-ray scattering data. If adequate cluster structures are chosen, the phase volume and a transition from the crystalline phase at lower temperatures to the amorphous phase at larger temperatures can be predicted. In agreement with experimental findings, we observed two phase transitions. The first one is computationally characterized by a discontinuous change in the phase volume, while the second, at a slightly larger temperature, is observed as a change in the dominating structural motifs from crystalline to amorphous. Both temperatures agree quantitatively with the experimental melting peaks.
While these results are yet of limited value due to the limited cluster set, they demonstrate the potential of QCE when applied to oligomers and potentially also to polymer materials. For example, if more clusters for the crystalline motifs are considered and environmental effects are taken into account for all clusters, a mechanistic understanding of the nature of the phase transition can be obtained, and a more quantitative agreement with X-ray measurements can be expected. Finally, this first demonstration of the applicability of the QCE to polymer systems has important implications for future studies. By exploiting further QCE features, such as mixtures in binary QCE (25) or the interpretation of low-populated clusters (27) regarding low-concentration defects, further questions, for example, regarding the role of additives or the properties of copolymers, can be addressed.

Supporting Information

Click to copy section linkSection link copied!

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcb.5c06385.

  • Snapshots of all cluster structures considered in the QCE calculation (Figures S1–S8) and the 1H NMR spectrum of OεCL1 (Figure S9) (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

Click to copy section linkSection link copied!

  • Corresponding Author
  • Authors
    • Mai Ahmed - Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, GermanyJena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, Germany
    • Deniz Yilmaz - Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, GermanyJena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, GermanyOrcidhttps://orcid.org/0009-0000-6291-4076
    • Purushottam Poudel - HIPOLE Jena (Helmholtz Institute for Polymers in Energy Applications Jena), Lessingstrasse 12-14, 07743 Jena, GermanyOrcidhttps://orcid.org/0000-0002-2452-2634
    • Felix H. Schacher - Institute of Organic Chemistry and Macromolecular Chemistry (IOMC), Friedrich Schiller University Jena, Humboldtstraße 10, 07743 Jena, GermanyJena Center for Soft Matter (JCSM), Friedrich Schiller University Jena, Philosophenweg 7, 07743 Jena, GermanyHIPOLE Jena (Helmholtz Institute for Polymers in Energy Applications Jena), Lessingstrasse 12-14, 07743 Jena, GermanyOrcidhttps://orcid.org/0000-0003-4685-6608
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

Click to copy section linkSection link copied!

The SAXS setup was funded by the Free State of Thuringia and REACT-EU (2021 FGI 0036─SAXS/WAXS). E.P., F.H.S., and P.P. were further supported by funding from the Carl Zeiss Foundation (“Durchbrüche”─Intelligente Substrate). Funding provided to M.A. within the ProChance program of FSU Jena is gratefully acknowledged.

References

Click to copy section linkSection link copied!

This article references 47 other publications.

  1. 1
    Gupta, V.; Biswas, D.; Roy, S. Materials 2022, 15 (17), 5899  DOI: 10.3390/ma15175899
  2. 2
    Tajeddin, B.; Arabkhedri, M. Polymer Science and Innovative Applications; AlMaadeed, M. A. A.; Ponnamma, D.; Carignano, M. A., Eds.; Elsevier, 2020; pp 525543.
  3. 3
    Hu, J.; Chen, S. J. Mater. Chem. 2010, 20, 33463355,  DOI: 10.1039/b922872a
  4. 4
    Moghadasi, K.; Mohd Isa, M. S.; Ariffin, M. A.; Mohd jamil, M. Z.; Raja, S.; Wu, B.; Yamani, M.; Bin Muhamad, M. R.; Yusof, F.; Jamaludin, M. F.; bin Ab Karim, M. S.; binti Abdul Razak, B.; bin Yusoff, N. J. Mater. Res. Technol. 2022, 17, 10541121
  5. 5
    Chen, W.-H.; Chen, Q.-W.; Chen, Q.; Cui, C.; Duan, S.; Kang, Y.; Liu, Y.; Liu, Y.; Muhammad, W.; Shao, S.; Tang, C.; Wang, J.; Wang, L.; Xiong, M.-H.; Yin, L.; Zhang, K.; Zhang, Z.; Zhen, X.; Feng, J.; Chen, X. Sci. China: Chem. 2022, 65, 10101075,  DOI: 10.1007/s11426-022-1243-5
  6. 6
    Ding, L.; Yu, Z.-D.; Wang, X.-Y.; Yao, Z.-F.; Lu, Y.; Yang, C.-Y.; Wang, J.-Y.; Pei, J. Chem. Rev. 2023, 123, 74217497,  DOI: 10.1021/acs.chemrev.2c00696
  7. 7
    Li, L.; Han, L.; Hu, H.; Zhang, R. Mater. Adv. 2023, 4, 726746,  DOI: 10.1039/D2MA00940D
  8. 8
    Campbell, I. R.; Lin, M.-Y.; Iyer, H.; Parker, M.; Fredricks, J. L.; Liao, K.; Jimenez, A. M.; Grandgeorge, P.; Roumeli, E. Annu. Rev. Mater. Res. 2023, 53, 81104,  DOI: 10.1146/annurev-matsci-080921-083655
  9. 9
    Jayakumar, A.; Radoor, S.; Siengchin, S.; Shin, G. H.; Kim, J. T. Sci. Total Environ. 2023, 878, 163156  DOI: 10.1016/j.scitotenv.2023.163156
  10. 10
    Sánchez, C. Biotechnol. Adv. 2020, 40, 107501  DOI: 10.1016/j.biotechadv.2019.107501
  11. 11
    EL-Ghoul, Y.; Alminderej, F. M.; Alsubaie, F. M.; Alrasheed, R.; Almousa, N. H. Polymers 2021, 13 (24), 4327  DOI: 10.3390/polym13244327
  12. 12
    Delfi, M.; Ghomi, M.; Zarrabi, A.; Mohammadinejad, R.; Taraghdari, Z. B.; Ashrafizadeh, M.; Zare, E. N.; Agarwal, T.; Padil, V. V. T.; Mokhtari, B.; Rossi, F.; Perale, G.; Sillanpaa, M.; Borzacchiello, A.; Kumar Maiti, T.; Makvandi, P. Prosthesis 2020, 2, 117139,  DOI: 10.3390/prosthesis2020012
  13. 13
    Nashchekina, Y.; Chabina, A.; Moskalyuk, O.; Voronkina, I.; Evstigneeva, P.; Vaganov, G.; Nashchekin, A.; Yudin, V.; Mikhailova, N. Polymers 2022, 14, 4654  DOI: 10.3390/polym14214654
  14. 14
    Shen, J.; Yuan, W.; Badv, M.; Moshaverinia, A.; Weiss, P. S. ACS Mater. Au 2023, 3, 540547,  DOI: 10.1021/acsmaterialsau.3c00027
  15. 15
    Leroux, A.; Ngoc Nguyen, T.; Rangel, A.; Cacciapuoti, I.; Duprez, D.; Castner, D. G.; Migonney, V. Biointerphases 2020, 15, 061006  DOI: 10.1116/6.0000429
  16. 16
    Feng, S.; Yue, Y.; Chen, J.; Zhou, J.; Li, Y.; Zhang, Q. Environ. Sci.: Processes Impacts 2020, 22, 23322344,  DOI: 10.1039/D0EM00340A
  17. 17
    Barber, V. J.; Borden, M. A.; Alty, J. W.; Tran, L. D.; Koerner, H.; Baldwin, L. A.; Alexanian, E. J.; Leibfarth, F. A. Macromolecules 2023, 56, 36793686,  DOI: 10.1021/acs.macromol.3c00125
  18. 18
    Castilla-Cortázar, I.; Más-Estellés, J.; Meseguer-Dueñas, J.; Escobar Ivirico, J.; Marí, B.; Vidaurre, A. Polym. Degrad. Stab. 2012, 97, 12411248,  DOI: 10.1016/j.polymdegradstab.2012.05.038
  19. 19
    Bhadran, A.; Shah, T.; Babanyinah, G. K.; Polara, H.; Taslimy, S.; Biewer, M. C.; Stefan, M. C. Pharmaceutics 2023, 15, 1977  DOI: 10.3390/pharmaceutics15071977
  20. 20
    Schmid, F. ACS Polym. Au 2023, 3, 2858,  DOI: 10.1021/acspolymersau.2c00049
  21. 21
    Kirchner, B. J. Chem. Phys. 2005, 123, 204116  DOI: 10.1063/1.2126977
  22. 22
    Kirchner, B. P. Phys. Rep. 2007, 440, 1111,  DOI: 10.1016/j.physrep.2006.11.005
  23. 23
    von Domaros, M.; Perlt, E.; Ingenmey, J.; Marchelli, G.; Kirchner, B. SoftwareX 2018, 7, 356359,  DOI: 10.1016/j.softx.2018.11.002
  24. 24
    Perlt, E.; Friedrich, J.; von Domaros, M.; Kirchner, B. ChemPhysChem 2011, 12, 34743482,  DOI: 10.1002/cphc.201100592
  25. 25
    Brüssel, M.; Perlt, E.; Lehmann, S. B. C.; von Domaros, M.; Kirchner, B. J. Chem. Phys. 2011, 135, 194113  DOI: 10.1063/1.3662071
  26. 26
    Ludwig, R.; Weinhold, F. J. Chem. Phys. 1999, 110, 508515,  DOI: 10.1063/1.478136
  27. 27
    Perlt, E.; von Domaros, M.; Kirchner, B.; Ludwig, R.; Weinhold, F. Sci. Rep. 2017, 7, 10244
  28. 28
    Blasius, J.; Ingenmey, J.; Perlt, E.; von Domaros, M.; Hollóczki, O.; Kirchner, B. Angew. Chem., Int. Ed. 2019, 58, 32123216,  DOI: 10.1002/anie.201811839
  29. 29
    Ren, J. M.; Fu, Q.; Blencowe, A.; Qiao, G. G. ACS Macro Lett. 2012, 1, 681686,  DOI: 10.1021/mz300169m
  30. 30
    El Habnouni, S.; Blanquer, S.; Darcos, V.; Coudane, J. J. Polym. Sci., Part A: Polym. Chem. 2009, 47, 61046115,  DOI: 10.1002/pola.23652
  31. 31
    O’Boyle, N. M.; Banck, M.; James, C. A.; Morley, C.; Vandermeersch, T.; Hutchison, G. R. J. Cheminf. 2011, 3, 33  DOI: 10.1186/1758-2946-3-33
  32. 32
    Balasubramani, S. G.; Chen, G. P.; Coriani, S.; Diedenhofen, M.; Frank, M. S.; Franzke, Y. J.; Furche, F.; Grotjahn, R.; Harding, M. E.; Hättig, C.; Hellweg, A.; Helmich-Paris, B.; Holzer, C.; Huniar, U.; Kaupp, M.; Marefat Khah, A.; Karbalaei Khani, S.; Müller, T.; Mack, F.; Nguyen, B. D.; Parker, S. M.; Perlt, E.; Rappoport, D.; Reiter, K.; Roy, S.; Rückert, M.; Schmitz, G.; Sierka, M.; Tapavicza, E.; Tew, D. P.; Van Wüllen, C.; Voora, V. K.; Weigend, F.; Wodyński, A.; Yu, J. M. J. Chem. Phys. 2020, 152, 184107  DOI: 10.1063/5.0004635
  33. 33
    Franzke, Y. J.; Holzer, C.; Andersen, J. H.; Begušić, T.; Bruder, F.; Coriani, S.; Della Sala, F.; Fabiano, E.; Fedotov, D. A.; Fürst, S.; Gillhuber, S.; Grotjahn, R.; Kaupp, M.; Kehry, M.; Krstić, M.; Mack, F.; Majumdar, S.; Nguyen, B. D.; Parker, S. M.; Pauly, F.; Pausch, A.; Perlt, E.; Phun, G. S.; Rajabi, A.; Rappoport, D.; Samal, B.; Schrader, T.; Sharma, M.; Tapavicza, E.; Treß, R. S.; Voora, V.; Wodyński, A.; Yu, J. M.; Zerulla, B.; Furche, F.; Hättig, C.; Sierka, M.; Tew, D. P.; Weigend, F. J. Chem. Theory Comput. 2023, 19, 68596890,  DOI: 10.1021/acs.jctc.3c00347
  34. 34
    Grimme, S.; Brandenburg, J. G.; Bannwarth, C.; Hansen, A. J. Chem. Phys. 2015, 143, 054107  DOI: 10.1063/1.4927476
  35. 35
    Adamo, C.; Barone, V. J. Chem. Phys. 1999, 110, 61586170,  DOI: 10.1063/1.478522
  36. 36
    Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. J. Chem. Phys. 2010, 132, 154104  DOI: 10.1063/1.3382344
  37. 37
    Grimme, S.; Ehrlich, S.; Goerigk, L. J. Comput. Chem. 2011, 32, 14561465,  DOI: 10.1002/jcc.21759
  38. 38
    Fernández-Tena, A.; Pérez-Camargo, R. A.; Coulembier, O.; Sangroniz, L.; Aranburu, N.; Guerrica-Echevarria, G.; Liu, G.; Wang, D.; Cavallo, D.; Müller, A. J. Effect of Molecular Weight on the Crystallization and Melt Memory of Poly(ε-caprolactone) (PCL). Macromolecules 2023, 56, 46024620,  DOI: 10.1021/acs.macromol.3c00234
  39. 39
    Ju, Y.; Zhang, M.; Zhao, H. Polym. Chem. 2017, 8, 54155426,  DOI: 10.1039/C7PY01012E
  40. 40
    Klonos, P. A.; Bikiaris, N. D.; Christodoulou, E.; Zamboulis, A.; Papageorgiou, G. Z.; Kyritsis, A. Polymer 2022, 242, 124603  DOI: 10.1016/j.polymer.2022.124603
  41. 41
    Persenaire, O.; Alexandre, M.; Degée, P.; Dubois, P. Biomacromolecules 2001, 2, 288294,  DOI: 10.1021/bm0056310
  42. 42
    Tenorio-Alfonso, A.; Ramos, E. V.; Martínez, I.; Ambrosi, M.; Raudino, M. Assessment of the structures contribution (crystalline and mesophases) and mechanical properties of polycaprolactone/pluronic blends. J. Mech. Behav. Biomed. Mater. 2023, 139, 105668  DOI: 10.1016/j.jmbbm.2023.105668
  43. 43
    Singh, P. S. Membrane Characterization; Hilal, N.; Ismail, A. F.; Matsuura, T.; Oatley-Radcliffe, D., Eds.; Elsevier, 2017; pp 95111.
  44. 44
    Nanaki, S. G.; Papageorgiou, G. Z.; Bikiaris, D. N. J. Therm. Anal. Calorim. 2012, 108, 633645,  DOI: 10.1007/s10973-011-2155-8
  45. 45
    Holzwarth, U.; Gibson, N. Nat. Nanotechnol. 2011, 6, 534,  DOI: 10.1038/nnano.2011.145
  46. 46
    Tencé-Girault, S.; Lebreton, S.; Bunau, O.; Dang, P.; Bargain, F. Crystals 2019, 9, 271  DOI: 10.3390/cryst9050271
  47. 47
    Wang, Z.; He, Y.; Müller, A. J. Polymer 2023, 271, 125783  DOI: 10.1016/j.polymer.2023.125783

Cited By

Click to copy section linkSection link copied!

This article has not yet been cited by other publications.

The Journal of Physical Chemistry B

Cite this: J. Phys. Chem. B 2026, 130, 5, 1675–1683
Click to copy citationCitation copied!
https://doi.org/10.1021/acs.jpcb.5c06385
Published January 27, 2026

Copyright © 2026 The Authors. Published by American Chemical Society. This publication is licensed under

CC-BY 4.0 .

Article Views

933

Altmetric

-

Citations

-
Learn about these metrics

Article Views are the COUNTER-compliant sum of full text article downloads since November 2008 (both PDF and HTML) across all institutions and individuals. These metrics are regularly updated to reflect usage leading up to the last few days.

Citations are the number of other articles citing this article, calculated by Crossref and updated daily. Find more information about Crossref citation counts.

The Altmetric Attention Score is a quantitative measure of the attention that a research article has received online. Clicking on the donut icon will load a page at altmetric.com with additional details about the score and the social media presence for the given article. Find more information on the Altmetric Attention Score and how the score is calculated.

  • Abstract

    Figure 1

    Figure 1. Example structures, as used in the QCE study. Crystalline motif (OεCL-c, left) as well as one exemplary amorphous structure (OεCL-a3, right).

    Scheme 1

    Scheme 1. Oligomer Formation of εCL with BnOH/SnOct2 as the Initiator/Catalyst

    Figure 2

    Figure 2. Characterization of the OεCL1 oligomer. Top, SEC chromatogram; middle, DSC heating and cooling scans; bottom, TGA curve.

    Figure 3

    Figure 3. 1D integrated (a) SAXS patterns of OεCL1 corresponding to q and (b) WAXS patterns of OεCL1 corresponding to 2θ.

    Figure 4

    Figure 4. Temperature-dependent cluster populations showing only significantly populated clusters (>1%). The gray line indicates the interception of both curves.

    Figure 5

    Figure 5. Phase volume of 1 mol OεCL. At 314 K, a sharp phase transition in the volume is observed.

  • References


    This article references 47 other publications.

    1. 1
      Gupta, V.; Biswas, D.; Roy, S. Materials 2022, 15 (17), 5899  DOI: 10.3390/ma15175899
    2. 2
      Tajeddin, B.; Arabkhedri, M. Polymer Science and Innovative Applications; AlMaadeed, M. A. A.; Ponnamma, D.; Carignano, M. A., Eds.; Elsevier, 2020; pp 525543.
    3. 3
      Hu, J.; Chen, S. J. Mater. Chem. 2010, 20, 33463355,  DOI: 10.1039/b922872a
    4. 4
      Moghadasi, K.; Mohd Isa, M. S.; Ariffin, M. A.; Mohd jamil, M. Z.; Raja, S.; Wu, B.; Yamani, M.; Bin Muhamad, M. R.; Yusof, F.; Jamaludin, M. F.; bin Ab Karim, M. S.; binti Abdul Razak, B.; bin Yusoff, N. J. Mater. Res. Technol. 2022, 17, 10541121
    5. 5
      Chen, W.-H.; Chen, Q.-W.; Chen, Q.; Cui, C.; Duan, S.; Kang, Y.; Liu, Y.; Liu, Y.; Muhammad, W.; Shao, S.; Tang, C.; Wang, J.; Wang, L.; Xiong, M.-H.; Yin, L.; Zhang, K.; Zhang, Z.; Zhen, X.; Feng, J.; Chen, X. Sci. China: Chem. 2022, 65, 10101075,  DOI: 10.1007/s11426-022-1243-5
    6. 6
      Ding, L.; Yu, Z.-D.; Wang, X.-Y.; Yao, Z.-F.; Lu, Y.; Yang, C.-Y.; Wang, J.-Y.; Pei, J. Chem. Rev. 2023, 123, 74217497,  DOI: 10.1021/acs.chemrev.2c00696
    7. 7
      Li, L.; Han, L.; Hu, H.; Zhang, R. Mater. Adv. 2023, 4, 726746,  DOI: 10.1039/D2MA00940D
    8. 8
      Campbell, I. R.; Lin, M.-Y.; Iyer, H.; Parker, M.; Fredricks, J. L.; Liao, K.; Jimenez, A. M.; Grandgeorge, P.; Roumeli, E. Annu. Rev. Mater. Res. 2023, 53, 81104,  DOI: 10.1146/annurev-matsci-080921-083655
    9. 9
      Jayakumar, A.; Radoor, S.; Siengchin, S.; Shin, G. H.; Kim, J. T. Sci. Total Environ. 2023, 878, 163156  DOI: 10.1016/j.scitotenv.2023.163156
    10. 10
      Sánchez, C. Biotechnol. Adv. 2020, 40, 107501  DOI: 10.1016/j.biotechadv.2019.107501
    11. 11
      EL-Ghoul, Y.; Alminderej, F. M.; Alsubaie, F. M.; Alrasheed, R.; Almousa, N. H. Polymers 2021, 13 (24), 4327  DOI: 10.3390/polym13244327
    12. 12
      Delfi, M.; Ghomi, M.; Zarrabi, A.; Mohammadinejad, R.; Taraghdari, Z. B.; Ashrafizadeh, M.; Zare, E. N.; Agarwal, T.; Padil, V. V. T.; Mokhtari, B.; Rossi, F.; Perale, G.; Sillanpaa, M.; Borzacchiello, A.; Kumar Maiti, T.; Makvandi, P. Prosthesis 2020, 2, 117139,  DOI: 10.3390/prosthesis2020012
    13. 13
      Nashchekina, Y.; Chabina, A.; Moskalyuk, O.; Voronkina, I.; Evstigneeva, P.; Vaganov, G.; Nashchekin, A.; Yudin, V.; Mikhailova, N. Polymers 2022, 14, 4654  DOI: 10.3390/polym14214654
    14. 14
      Shen, J.; Yuan, W.; Badv, M.; Moshaverinia, A.; Weiss, P. S. ACS Mater. Au 2023, 3, 540547,  DOI: 10.1021/acsmaterialsau.3c00027
    15. 15
      Leroux, A.; Ngoc Nguyen, T.; Rangel, A.; Cacciapuoti, I.; Duprez, D.; Castner, D. G.; Migonney, V. Biointerphases 2020, 15, 061006  DOI: 10.1116/6.0000429
    16. 16
      Feng, S.; Yue, Y.; Chen, J.; Zhou, J.; Li, Y.; Zhang, Q. Environ. Sci.: Processes Impacts 2020, 22, 23322344,  DOI: 10.1039/D0EM00340A
    17. 17
      Barber, V. J.; Borden, M. A.; Alty, J. W.; Tran, L. D.; Koerner, H.; Baldwin, L. A.; Alexanian, E. J.; Leibfarth, F. A. Macromolecules 2023, 56, 36793686,  DOI: 10.1021/acs.macromol.3c00125
    18. 18
      Castilla-Cortázar, I.; Más-Estellés, J.; Meseguer-Dueñas, J.; Escobar Ivirico, J.; Marí, B.; Vidaurre, A. Polym. Degrad. Stab. 2012, 97, 12411248,  DOI: 10.1016/j.polymdegradstab.2012.05.038
    19. 19
      Bhadran, A.; Shah, T.; Babanyinah, G. K.; Polara, H.; Taslimy, S.; Biewer, M. C.; Stefan, M. C. Pharmaceutics 2023, 15, 1977  DOI: 10.3390/pharmaceutics15071977
    20. 20
      Schmid, F. ACS Polym. Au 2023, 3, 2858,  DOI: 10.1021/acspolymersau.2c00049
    21. 21
      Kirchner, B. J. Chem. Phys. 2005, 123, 204116  DOI: 10.1063/1.2126977
    22. 22
      Kirchner, B. P. Phys. Rep. 2007, 440, 1111,  DOI: 10.1016/j.physrep.2006.11.005
    23. 23
      von Domaros, M.; Perlt, E.; Ingenmey, J.; Marchelli, G.; Kirchner, B. SoftwareX 2018, 7, 356359,  DOI: 10.1016/j.softx.2018.11.002
    24. 24
      Perlt, E.; Friedrich, J.; von Domaros, M.; Kirchner, B. ChemPhysChem 2011, 12, 34743482,  DOI: 10.1002/cphc.201100592
    25. 25
      Brüssel, M.; Perlt, E.; Lehmann, S. B. C.; von Domaros, M.; Kirchner, B. J. Chem. Phys. 2011, 135, 194113  DOI: 10.1063/1.3662071
    26. 26
      Ludwig, R.; Weinhold, F. J. Chem. Phys. 1999, 110, 508515,  DOI: 10.1063/1.478136
    27. 27
      Perlt, E.; von Domaros, M.; Kirchner, B.; Ludwig, R.; Weinhold, F. Sci. Rep. 2017, 7, 10244
    28. 28
      Blasius, J.; Ingenmey, J.; Perlt, E.; von Domaros, M.; Hollóczki, O.; Kirchner, B. Angew. Chem., Int. Ed. 2019, 58, 32123216,  DOI: 10.1002/anie.201811839
    29. 29
      Ren, J. M.; Fu, Q.; Blencowe, A.; Qiao, G. G. ACS Macro Lett. 2012, 1, 681686,  DOI: 10.1021/mz300169m
    30. 30
      El Habnouni, S.; Blanquer, S.; Darcos, V.; Coudane, J. J. Polym. Sci., Part A: Polym. Chem. 2009, 47, 61046115,  DOI: 10.1002/pola.23652
    31. 31
      O’Boyle, N. M.; Banck, M.; James, C. A.; Morley, C.; Vandermeersch, T.; Hutchison, G. R. J. Cheminf. 2011, 3, 33  DOI: 10.1186/1758-2946-3-33
    32. 32
      Balasubramani, S. G.; Chen, G. P.; Coriani, S.; Diedenhofen, M.; Frank, M. S.; Franzke, Y. J.; Furche, F.; Grotjahn, R.; Harding, M. E.; Hättig, C.; Hellweg, A.; Helmich-Paris, B.; Holzer, C.; Huniar, U.; Kaupp, M.; Marefat Khah, A.; Karbalaei Khani, S.; Müller, T.; Mack, F.; Nguyen, B. D.; Parker, S. M.; Perlt, E.; Rappoport, D.; Reiter, K.; Roy, S.; Rückert, M.; Schmitz, G.; Sierka, M.; Tapavicza, E.; Tew, D. P.; Van Wüllen, C.; Voora, V. K.; Weigend, F.; Wodyński, A.; Yu, J. M. J. Chem. Phys. 2020, 152, 184107  DOI: 10.1063/5.0004635
    33. 33
      Franzke, Y. J.; Holzer, C.; Andersen, J. H.; Begušić, T.; Bruder, F.; Coriani, S.; Della Sala, F.; Fabiano, E.; Fedotov, D. A.; Fürst, S.; Gillhuber, S.; Grotjahn, R.; Kaupp, M.; Kehry, M.; Krstić, M.; Mack, F.; Majumdar, S.; Nguyen, B. D.; Parker, S. M.; Pauly, F.; Pausch, A.; Perlt, E.; Phun, G. S.; Rajabi, A.; Rappoport, D.; Samal, B.; Schrader, T.; Sharma, M.; Tapavicza, E.; Treß, R. S.; Voora, V.; Wodyński, A.; Yu, J. M.; Zerulla, B.; Furche, F.; Hättig, C.; Sierka, M.; Tew, D. P.; Weigend, F. J. Chem. Theory Comput. 2023, 19, 68596890,  DOI: 10.1021/acs.jctc.3c00347
    34. 34
      Grimme, S.; Brandenburg, J. G.; Bannwarth, C.; Hansen, A. J. Chem. Phys. 2015, 143, 054107  DOI: 10.1063/1.4927476
    35. 35
      Adamo, C.; Barone, V. J. Chem. Phys. 1999, 110, 61586170,  DOI: 10.1063/1.478522
    36. 36
      Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. J. Chem. Phys. 2010, 132, 154104  DOI: 10.1063/1.3382344
    37. 37
      Grimme, S.; Ehrlich, S.; Goerigk, L. J. Comput. Chem. 2011, 32, 14561465,  DOI: 10.1002/jcc.21759
    38. 38
      Fernández-Tena, A.; Pérez-Camargo, R. A.; Coulembier, O.; Sangroniz, L.; Aranburu, N.; Guerrica-Echevarria, G.; Liu, G.; Wang, D.; Cavallo, D.; Müller, A. J. Effect of Molecular Weight on the Crystallization and Melt Memory of Poly(ε-caprolactone) (PCL). Macromolecules 2023, 56, 46024620,  DOI: 10.1021/acs.macromol.3c00234
    39. 39
      Ju, Y.; Zhang, M.; Zhao, H. Polym. Chem. 2017, 8, 54155426,  DOI: 10.1039/C7PY01012E
    40. 40
      Klonos, P. A.; Bikiaris, N. D.; Christodoulou, E.; Zamboulis, A.; Papageorgiou, G. Z.; Kyritsis, A. Polymer 2022, 242, 124603  DOI: 10.1016/j.polymer.2022.124603
    41. 41
      Persenaire, O.; Alexandre, M.; Degée, P.; Dubois, P. Biomacromolecules 2001, 2, 288294,  DOI: 10.1021/bm0056310
    42. 42
      Tenorio-Alfonso, A.; Ramos, E. V.; Martínez, I.; Ambrosi, M.; Raudino, M. Assessment of the structures contribution (crystalline and mesophases) and mechanical properties of polycaprolactone/pluronic blends. J. Mech. Behav. Biomed. Mater. 2023, 139, 105668  DOI: 10.1016/j.jmbbm.2023.105668
    43. 43
      Singh, P. S. Membrane Characterization; Hilal, N.; Ismail, A. F.; Matsuura, T.; Oatley-Radcliffe, D., Eds.; Elsevier, 2017; pp 95111.
    44. 44
      Nanaki, S. G.; Papageorgiou, G. Z.; Bikiaris, D. N. J. Therm. Anal. Calorim. 2012, 108, 633645,  DOI: 10.1007/s10973-011-2155-8
    45. 45
      Holzwarth, U.; Gibson, N. Nat. Nanotechnol. 2011, 6, 534,  DOI: 10.1038/nnano.2011.145
    46. 46
      Tencé-Girault, S.; Lebreton, S.; Bunau, O.; Dang, P.; Bargain, F. Crystals 2019, 9, 271  DOI: 10.3390/cryst9050271
    47. 47
      Wang, Z.; He, Y.; Müller, A. J. Polymer 2023, 271, 125783  DOI: 10.1016/j.polymer.2023.125783
  • Supporting Information

    Supporting Information


    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jpcb.5c06385.

    • Snapshots of all cluster structures considered in the QCE calculation (Figures S1–S8) and the 1H NMR spectrum of OεCL1 (Figure S9) (PDF)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.