The Comprehensive Approach to Preparation and Investigation of the Eu3+ Doped Hydroxyapatite/poly(L-lactide) Nanocomposites: Promising Materials for Theranostics Application

In response to the need for new materials for theranostics application, the structural and spectroscopic properties of composites designed for medical applications, received in the melt mixing process, were evaluated. A composite based on medical grade poly(L-lactide) (PLLA) and calcium hydroxyapatite (HAp) doped with Eu3+ ions was obtained by using a twin screw extruder. Pure calcium Hap, as well as the one doped with Eu3+ ions, was prepared using the precipitation method and then used as a filler. XRPD (X-ray Powder Diffraction) and IR (Infrared) spectroscopy were applied to investigate the structural properties of the obtained materials. DSC (Differential Scanning Calorimetry) was used to assess the Eu3+ ion content on phase transitions in PLLA. The tensile properties were also investigated. The excitation, emission spectra as well as decay time were measured to determine the spectroscopic properties. The simplified Judd–Ofelt (J-O) theory was applied and a detailed analysis in connection with the observed structural and spectroscopic measurements was made and described.


Introduction
Regenerative medicine comes in many forms, promising to develop new biomedical treatments for people who suffer due to the burden of trauma, congenital defects and degenerative diseases. Recently, the greatest effort has been put in seeking materials that could help and accelerate the regenerative process by stimulating the body's own repair mechanisms to functionally heal previously irreparable tissues or organs [1]. Nowadays, the materials with great biological properties, broadly used in medicine, are calcium phosphates and poly(L-lactide) (PLLA) [2].
The biocompatible and biodegradable poly(L-lactide) is one of the most popular polymers produced from renewable raw materials [3]. It can enhance the adhesion and elasticity of composite materials, therefore, it is broadly used in medical and biomedical applications, i.e., for vascular stent production [4], for drug delivery systems [5], for tissue engineering [6], and as surgical sutures, implants and screws. Moreover, it is used in such fields as water purification by oil adsorption [7,8], food europium nitrate, and then the europium nitrate hydrates were re-crystallized three times in order to eliminate of the HNO 3 excess. Afterwards, the stoichiometric amount of calcium nitrate and europium nitrate was dissolved in deionized water and mixed together. Then the stoichiometric amount of (NH 4 ) 2 HPO 4 was added to the previous mixture, leading to the fast precipitation of the intermediate product. The pH value of the suspension was adjusted to 10 by ammonia. The reaction mixture was heated and stirred for 3 h. Subsequently, the obtained products were washed several times with de-ionized water and dried at 70 • C for 24 h. As the final products, three types of europium-doped hydroxyapatites (1 mol% Eu 3+ :HAp, 3 mol% Eu 3+ :HAp, 5 mol% Eu 3+ :HAp) were obtained. Calcium hydroxyapatite (HAp) was also prepared as a reference sample.

Preparation of the PLLA/Eu 3+ -Doped Hydroxyapatite Composites
PLLA/Eu 3+ :HAp composites having 10 wt.% of the fillers were prepared using the Thermo Scientific Process 11 (Waltham, MA, USA) co-rotating twin-screw micro-extruder (D = 20 mm, L/D = 40) with screw rotation speed of 200 min −1 and barrel temperature profile of 200-180 • C (from hopper to die) in a nitrogen atmosphere. The extruder screw geometry was presented in our previous work [44]. The composites were extruded in a one-step process. PLLA-and Eu 3+ -doped HAp were dried in 80 • C under vacuum for 4 h before compounding. After extrusion, the composites were cooled down in the air and pelletized.
The composites were formed into foils through the casting extrusion technique using the Ultra Micro Cast Film extruder (Labtech Engineering, Sweden/Thailand) having flat die with a width of 75 mm, a conical screw with a diameter ranging from 18 to 8 mm (from hopper to die), L/d = 24, temperature of extrusion of 200 • C and screw speed of 100 rpm. Using these parameters, foils with a thickness of~100 µm were obtained. The obtained materials are summarized in Table 1.

Characterization
Powder X-ray diffraction patterns were measured in the 2θ range of 2-60 • by using a Rigaku Ultima IV (Tokyo, Japan) X-ray diffractometer equipped with Ni-filtered Cu Kα 1 radiation (Kα 1 = 1.54060 Å, U = 40 kV, I = 30 mA). The experimental XRD patterns were compared with the standards obtained from the Inorganic Crystal Structure Database (ICSD) and analyzed. The degree of crystallinity from XRD patterns was calculated using deconvolution made in Origin 8.0.
High-resolution transmission electron microscopy (HRTEM) images were done by a Philips CM-20 SuperTwin microscope (Eindhoven, The Netherlands), operating at 200 kV. Sample was prepared by dispersing a small amount of specimen in methanol and putting a droplet of the suspension on a copper microscope gird covered with carbon.
The surface morphology and element mapping of the nanocomposites were observed with a scanning electron microscope equipped with energy dispersive spectroscopy FEI Nova NanoSEM 230 Nanomaterials 2019, 9, 1146 4 of 16 (Hillsboro, OR, USA) with an EDS spectrometer (EDAX Genesis XM4) at an acceleration voltage of 18 kV and spot 3.0. Before observation, a layer of graphite was sprayed uniformly over the samples. IR spectra were acquired using a Thermo Scientific Nicolet iS10 FT-IR Spectrometer (Waltham, MA, USA) equipped with an Automated Beamsplitter exchange system (iS50 ABX containing a DLaTGS KBr detector), a built-in all-reflective diamond ATR module (iS50 ATR), Thermo Scientific Polaris™ and a HeNe laser as an IR radiation source. IR spectra were recorded at 295 K temperature in the 4000-500 cm −1 range in KBr pellets with a spectral resolution of 2 cm −1 .
The emission, excitation spectra, and luminescence kinetics were recorded using an FLS980 fluorescence spectrometer (Edinburgh Instruments, Kirkton Campus, UK). A 450 W xenon lamp was used as an excitation source. The radiation from the lamp was filtered by a 300 mm monochromator equipped with holographic grating (1800 grooves per mm, blaze 250 nm) for the emission and excitation spectra. A microsecond flashlamp (µF2) was used for the measurements of luminescence kinetics, whereas a Hamamatsu R928P photomultiplier was used as a detector. All the emission and excitation spectra were corrected according to the apparatus characteristics and the excitation source intensity. The luminescence kinetics was recorded at 616 nm according to electric dipole transition ( 5 D 0 → 7 F 2 ). The powders were placed in quarc tube and the composites were placed directly in the holder in the spectrometer.
The temperatures of phase transitions were evaluated using the DSC1 STARe Differential Scanning Calorimeter System from Mettler-Toledo (Giessen, Germany). The research was conducted in a nitrogen environment for a temperature range of 25 ÷ 200 • C. The heating/cooling rate was 5 K/min and the gas flow rate remained at 20 mL/min (two thermal cycles).
Samples for mechanical tests were cut out from the obtained foils, from the middle region having 100 µm, using a 5A type cutter (PN-EN ISO 527-2). Tensile properties tests were conducted on a universal Instron 5966 (Norwood, MA, USA) machine with the speed of 1 mm/min (for Young's modulus) and 10 mm/min (for tensile strength and strain at break). Before measurements, samples were kept in 80 • C for 2 h to reduce internal stresses.

Morphology
The morphology of the received hydroxyapatite nanopowders doped with 3 mol% of Eu 3+ was examined using the TEM and Selected Area Electron Diffraction (SAED) techniques. The received nanoapatite is nanocrystalline, loosely aggregated, and its shape is irregular (see Figure 1). The crystal phase purity of the obtained hydroxyapatite was additionally proved by a SAED analysis. The particles size distribution based on SEM images of pure HAp and 3 mol% Eu 3+ :HAp particles, as well as representative SEM images of powders were presented on Figure S1. DLaTGS KBr detector), a built-in all-reflective diamond ATR module (iS50 ATR), Thermo Scientific Polaris™ and a HeNe laser as an IR radiation source. IR spectra were recorded at 295 K temperature in the 4000-500 cm −1 range in KBr pellets with a spectral resolution of 2 cm −1 . The emission, excitation spectra, and luminescence kinetics were recorded using an FLS980 fluorescence spectrometer (Edinburgh Instruments, Kirkton Campus, UK). A 450 W xenon lamp was used as an excitation source. The radiation from the lamp was filtered by a 300 mm monochromator equipped with holographic grating (1800 grooves per mm, blaze 250 nm) for the emission and excitation spectra. A microsecond flashlamp (μF2) was used for the measurements of luminescence kinetics, whereas a Hamamatsu R928P photomultiplier was used as a detector. All the emission and excitation spectra were corrected according to the apparatus characteristics and the excitation source intensity. The luminescence kinetics was recorded at 616 nm according to electric dipole transition ( 5 D0 → 7 F2). The powders were placed in quarc tube and the composites were placed directly in the holder in the spectrometer.
The temperatures of phase transitions were evaluated using the DSC1 STARe Differential Scanning Calorimeter System from Mettler-Toledo (Giessen, Germany). The research was conducted in a nitrogen environment for a temperature range of 25 ÷ 200 °C. The heating/cooling rate was 5 K/min and the gas flow rate remained at 20 mL/min (two thermal cycles).
Samples for mechanical tests were cut out from the obtained foils, from the middle region having ~100 μm, using a 5A type cutter (PN-EN ISO 527-2). Tensile properties tests were conducted on a universal Instron 5966 (Norwood, MA, USA) machine with the speed of 1 mm/min (for Young's modulus) and 10 mm/min (for tensile strength and strain at break). Before measurements, samples were kept in 80 °C for 2 h to reduce internal stresses.

Morphology
The morphology of the received hydroxyapatite nanopowders doped with 3 mol% of Eu 3+ was examined using the TEM and Selected Area Electron Diffraction (SAED) techniques. The received nanoapatite is nanocrystalline, loosely aggregated, and its shape is irregular (see Figure 1). The crystal phase purity of the obtained hydroxyapatite was additionally proved by a SAED analysis. The particles size distribution based on SEM images of pure HAp and 3 mol% Eu 3+ :HAp particles, as well as representative SEM images of powders were presented on Figure S1. The SEM pictures of the PLLA/HAp composites are shown in Figures 2, S2 and S3. All the materials consist of 90 wt.% of PLLA and 10 wt.% of hydroxyapatite. The main difference between the materials is the Eu 3+ content in the hydroxyapatite. The pictures indicate a random distribution of fillers in all the systems. The hydroxyapatite particle size in all the composites ranges from a few to over one hundred nanometers. The Eu 3+ addition does not affect the distribution of HAp in the composites. The EDAX picture ( Figure S4) also shows the distribution of basic elements (Ca, P) in the PLLA/5 mol% Eu 3+ :HAp. Additionally, the figure presents the distribution of Eu 3+ in the system. to over one hundred nanometers. The Eu 3+ addition does not affect the distribution of HAp in the composites. The EDAX picture ( Figure S4) also shows the distribution of basic elements (Ca, P) in the PLLA/5 mol% Eu 3+ :HAp. Additionally, the figure presents the distribution of Eu 3+ in the system.

Structural Analysis of Eu 3+ -Doped Composites
The crystal phase purity of Ca10(PO4)6(OH)2 nanocrystals doped with x Eu 3+ ions (where x = 1, 3, 5 mol%) was checked with the powder XRD technique and was compared with the reference standard of the hexagonal Ca10(PO4)6(OH)2 lattice ascribed to the P63/m space group [45] (ICSD-180315) (see Figure S5). The diffraction patterns of the Eu 3+ :Ca10(PO4)6(OH)2 nanoparticles embedded into the poly(L-lactide) composite having the function of optically active ions concentration obtained with the extrusion method are presented in Figure 3. The neat poly(L-lactide) crystallized in the orthorhombic α′-polylactide with P212121 space group. In this case, the degree of crystallinity was approximately 11% ( Figure S6a). The deconvolution of XRD curves of the hydroxyapatite-doped PLLA, which revealed amorphous PLLA in all the investigated cases, was performed (see Figure S6be). Interestingly, mesophase with an amount in the range of 1.9% to 5.9% was also found in PLLA doped with HAp. The highest mesophase content was found in the systems having the highest Eu 3+ content (3 and 5 mol%).
Infrared spectra were recorded for the entire materials to get a deeper insight into the structure. The IR spectra of HAp consist of typical active vibrational bands of phosphate and hydroxyl groups (see Figure S7).
The IR spectra of the PLLA/HAp composites comprise vibrational bands related to the hydroxyapatite and poly(L-lactide) (see Figure 4). The typical poly(L-lactide) active stretching and banding vibrations of the -CH3 and -CH2 groups are observed at 2848.3 and 2915.4 cm −1 , respectively. The most intense peak at 1749.1 cm −1 is connected with -C = O stretching vibration. The symmetric and asymmetric stretching of the C-C(=O)-O group was observed at 1452.1 and 1150 cm −1 . Peaks at this range can be connected with ester groups existing in polylactide and lactide molecules. The 1380.7 cm −1 and 1361.0 cm −1 peaks were assigned to the scissor vibration δs(-CH3) group and the bending vibration of the δ1(-CH3) group. The vibration of the ester group (-C-O-) derived from poly(L-lactide) molecules can also be clearly observed at 1180.7 cm −1 . Peaks lying at higher energy corresponded to the vibration of groups belonging to HAp. Regarding the infrared spectra of pure hydroxyapatite, it was possible to determine all of the functional groups of hydroxyapatite. In the case of polylactide composites, the position of phosphate bands is shifted slightly. The triply degenerate ν3 antisymmetric stretching of the phosphate groups is observed at 1039.4 cm −1 (ν3) and 1082.8 cm −1 (ν3). The ratio intensity of these peaks in reference to nanopowders ( Figure 4) is different. In the composite, the peak around 1080 cm −1 is more intense, because the ester group vibration (-C-O-) of poly(L-lactide) and the phosphate group vibration (v3(PO4 3− )) of hydroxyapatite are overlapped. The symmetric stretching (ν1) of the PO4 3− groups was assigned to the peak at 955.5 cm −1 .

Structural Analysis of Eu 3+ -Doped Composites
The crystal phase purity of Ca 10 (PO 4 ) 6 (OH) 2 nanocrystals doped with x Eu 3+ ions (where x = 1, 3, 5 mol%) was checked with the powder XRD technique and was compared with the reference standard of the hexagonal Ca 10 (PO 4 ) 6 (OH) 2 lattice ascribed to the P6 3 /m space group [45] (ICSD-180315) (see Figure S5). The diffraction patterns of the Eu 3+ :Ca 10 (PO 4 ) 6 (OH) 2 nanoparticles embedded into the poly(L-lactide) composite having the function of optically active ions concentration obtained with the extrusion method are presented in Figure 3. The neat poly(L-lactide) crystallized in the orthorhombic α -polylactide with P2 1 2 1 2 1 space group. In this case, the degree of crystallinity was approximately 11% ( Figure S6a). The deconvolution of XRD curves of the hydroxyapatite-doped PLLA, which revealed amorphous PLLA in all the investigated cases, was performed (see Figure S6b-e). Interestingly, mesophase with an amount in the range of 1.9% to 5.9% was also found in PLLA doped with HAp. The highest mesophase content was found in the systems having the highest Eu 3+ content (3 and 5 mol%).
Infrared spectra were recorded for the entire materials to get a deeper insight into the structure. The IR spectra of HAp consist of typical active vibrational bands of phosphate and hydroxyl groups (see Figure S7).
The IR spectra of the PLLA/HAp composites comprise vibrational bands related to the hydroxyapatite and poly(L-lactide) (see Figure 4). The typical poly(L-lactide) active stretching and banding vibrations of the -CH 3 and -CH 2 groups are observed at 2848.3 and 2915.4 cm −1 , respectively. The most intense peak at 1749.1 cm −1 is connected with -C = O stretching vibration. The symmetric and asymmetric stretching of the C-C(=O)-O group was observed at 1452.1 and 1150 cm −1 . Peaks at this range can be connected with ester groups existing in polylactide and lactide molecules. The 1380.7 cm −1 and 1361.0 cm −1 peaks were assigned to the scissor vibration δ s (-CH 3 ) group and the bending vibration of the δ 1 (-CH 3 ) group. The vibration of the ester group (-C-O-) derived from poly(L-lactide) molecules can also be clearly observed at 1180.7 cm −1 . Peaks lying at higher energy corresponded to the vibration of groups belonging to HAp. Regarding the infrared spectra of pure hydroxyapatite, it was possible to determine all of the functional groups of hydroxyapatite. In the case of polylactide composites, the position of phosphate bands is shifted slightly. The triply degenerate ν 3 antisymmetric stretching of the phosphate groups is observed at 1039.4 cm −1 (ν 3 ) and 1082.8 cm −1 (ν 3 ). The ratio intensity of these peaks in reference to nanopowders (Figure 4) is different. In the composite, the peak around 1080 cm −1 is more intense, because the ester group vibration (-C-O-) of poly(L-lactide) and the phosphate group vibration (v 3 (PO 4 3− )) of hydroxyapatite are overlapped. The symmetric stretching (ν 1 ) of the PO 4 3− groups was assigned to the peak at 955.5 cm −1 . The spectra of all obtained composites show great similarity and are in agreement with the literature data [46][47][48].
Nanomaterials 2019, 9, x FOR PEER REVIEW 6 of 16 The spectra of all obtained composites show great similarity and are in agreement with the literature data [46][47][48].   The spectra of all obtained composites show great similarity and are in agreement with the literature data [46][47][48].

Thermal Properties
More information about the supramolecular structure of PLLA in the composites was obtained from DSC measurements. Upon the first heating, the DSC curves of PLLA and the composites are presented in Figure 5a, characteristic phase transitions are visible, such as glass transition with T g at around 61 • C, cold crystallization at the mean temperature range of 80 • C ÷ 120 • C, and melting with preceding small exothermic effect corresponding to alpha'-alpha reorganization [49,50]. The thermal parameters estimated form the first heating DSC curves are collected in Table S1. around 61°C, cold crystallization at the mean temperature range of 80 °C ÷ 120 °C, and melting with preceding small exothermic effect corresponding to alpha'-alpha reorganization [49,50]. The thermal parameters estimated form the first heating DSC curves are collected in Table S1.
The influence of HAp particles and Eu 3+ ions doping HAp on melt crystallization of PLLA can be analyzed based on the cooling DSC curves presented in Figure 5b and Table S2. For comparison of the cooling DSC curves, it is clearly visible that melt crystallization exotherms of PLLA in the analyzed samples have different temperature ranges. Neat PLLA and PLLA with unmodified HAp crystallize upon cooling at lower temperature ranges than PLLA with HAp modified with Eu 3+ . The Eu 3+ ions have an influence on crystallization kinetics of PLLA. In the presence of Eu 3+ , PLLA crystallizes faster. Moreover, it is known that at lower temperatures, it crystallizes into alpha and alpha' crystals whereas at higher temperatures, PLLA crystallizes in the alpha form [49]. It can be expected that the presence of Eu 3+ has an influence on the crystalline form of PLLA that forms on melt crystallization. The verification of this assumption can be found in the second heating DSC curves (Figure 5c, Table S3).  The influence of HAp particles and Eu 3+ ions doping HAp on melt crystallization of PLLA can be analyzed based on the cooling DSC curves presented in Figure 5b and Table S2. For comparison of the cooling DSC curves, it is clearly visible that melt crystallization exotherms of PLLA in the analyzed samples have different temperature ranges. Neat PLLA and PLLA with unmodified HAp crystallize upon cooling at lower temperature ranges than PLLA with HAp modified with Eu 3+ . The Eu 3+ ions have an influence on crystallization kinetics of PLLA. In the presence of Eu 3+ , PLLA crystallizes faster. Moreover, it is known that at lower temperatures, it crystallizes into alpha and alpha' crystals whereas at higher temperatures, PLLA crystallizes in the alpha form [49]. It can be expected that the presence of Eu 3+ has an influence on the crystalline form of PLLA that forms on melt crystallization. The verification of this assumption can be found in the second heating DSC curves (Figure 5c, Table S3).

Spectroscopic Properties
The excitation emission spectra of all the obtained materials were measured at room temperature, with a recording emission wavelength at 616 nm, which responds to the maximum of the most intense electric dipole transition ( 5 D 0 → 7 F 2 ). The spectra were corrected to the intensity of the excitation source and normalized to the most intense transition. All spectra contain characteristic sharp lines associated with the intraconfigurational 4f-4f transitions of the Eu 3+ ions, as well as an intense and broad band related to the ligand-to-metal charge transfer (CT) O 2− → Eu 3+ transition located in the UV region. The excitation emission spectra presented in Figure S8 belong to Eu 3+ :HAp nanopowders and in Figure 6 belong to PLLA/Eu 3+ :HAp composites. The 4f orbitals of lanthanide ions are well isolated by the external 5s, 5p and 5d shells and well protected against the influence of the crystal field. Due to this feature, the barycenters of the f-f lanthanide electron transitions are weakly affected by the ligand field, and the position of these peaks remains almost independent of the host lattice structure [51,52]. The narrow peaks observed at 299.2 nm (33 422 cm −1 ) were ascribed to the 7 F 0 → 5 F (4,3,2,1) , 3 P 0 transitions, at 319.7 nm (31 279 cm −1 ) to 7 F 0 → 5 H (6,5,4,7,3) , at 363.4 nm (27 518 cm −1 ) to 7 F 0 → 5 D 4 , 5 L 8 at 377.0 nm (26 525 cm −1 ) to 5 L 8 , 7 F 0 → G 2 , 5 L 7 , 5 G 3 , at 394.4 nm (25 355 cm −1 ) to 7 F 0 → 5 L 6 , at 416.2 nm (24 027 cm −1 ) to 7 F 0 → 5 D 3 at 465.9 nm (21 464 cm −1 ) to 7 F 0 → 5 D 2 , as well as at 527.0 nm (18 975 cm −1 ) to 7 F 0 → 5 D 1 . The allowed CT transition is strongly affected by electron-lattice coupling, and the peak position depends on the surrounding symmetry of the ion. In the case of PLLA/HAp composites, the CT maximum is located at 253.5 nm (39 448 cm −1 ). It is well known that the substitution of divalent Ca 2+ by trivalent Eu 3+ cations requires the charge compensation mechanism. The two mechanisms are well known for the apatite host lattice and have been described previously in [51,53,54].

Spectroscopic Properties
The excitation emission spectra of all the obtained materials were measured at room temperature, with a recording emission wavelength at 616 nm, which responds to the maximum of the most intense electric dipole transition ( 5 D0 → 7 F2). The spectra were corrected to the intensity of the excitation source and normalized to the most intense transition. All spectra contain characteristic sharp lines associated with the intraconfigurational 4f-4f transitions of the Eu 3+ ions, as well as an intense and broad band related to the ligand-to-metal charge transfer (CT) O 2− → Eu 3+ transition located in the UV region. The excitation emission spectra presented in Figure S8 belong to Eu 3+ :HAp nanopowders and in Figure 6 belong to PLLA/Eu 3+ :HAp composites. The 4f orbitals of lanthanide ions are well isolated by the external 5s, 5p and 5d shells and well protected against the influence of the crystal field. Due to this feature, the barycenters of the f-f lanthanide electron transitions are weakly affected by the ligand field, and the position of these peaks remains almost independent of the host lattice structure [51,52]. The narrow peaks observed at 299.2 nm (33 422 cm −1 ) were ascribed to the 7 F0 → 5 F(4,3,2,1), 3 P0 transitions, at 319.7 nm (31 279 cm −1 ) to 7 F0 → 5 H(6,5,4,7,3), at 363.4 nm (27 518 cm −1 ) to 7 F0 → 5 D4, 5 L8 at 377.0 nm (26 525 cm −1 ) to 5 L8, 7 F0 → G2, 5 L7, 5 G3, at 394.4 nm (25 355 cm −1 ) to 7 F0 → 5 L6, at 416.2 nm (24 027 cm −1 ) to 7 F0 → 5 D3 at 465.9 nm (21 464 cm −1 ) to 7 F0 → 5 D2, as well as at 527.0 nm (18 975 cm −1 ) to 7 F0 → 5 D1. The allowed CT transition is strongly affected by electron-lattice coupling, and the peak position depends on the surrounding symmetry of the ion. In the case of PLLA/HAp composites, the CT maximum is located at 253.5 nm (39 448 cm −1 ). It is well known that the substitution of divalent Ca 2+ by trivalent Eu 3+ cations requires the charge compensation mechanism. The two mechanisms are well known for the apatite host lattice and have been described previously in [51,53,54]. The room temperature emission spectra of the entire nanopowders and nanocomposites were measured in the spectral range of 500-750 nm under 394.5 nm excitation as a function of optically The room temperature emission spectra of the entire nanopowders and nanocomposites were measured in the spectral range of 500-750 nm under 394.5 nm excitation as a function of optically active ions concentration. The emission spectra of Eu 3+ :HAp nanopowders are collected in Figure S9 and the emission spectra of PLLA/Eu 3+ :HAp composites are collected in Figure 7. The spectra were normalized to the 5 D 0 → 7 F 1 transition. As can be seen, the emission spectra of Eu 3+ ions are composed of five bands related to the 5 D 0 → 7 F 0,1,2,3,4 transitions occurring respectively at ca. 573.5 nm (17 437 cm −1 ), 590.0 nm (16 949 cm −1 ), 616.2 nm (16 228 cm −1 ), 653.2 nm (15 309 cm −1 ) and 700.9 nm (14 267 cm −1 ). The most important transitions in the study of the structural and spectroscopic properties of Eu 3+ ions are the 5 D 0 → 7 F 0,1,2 transitions. Nanomaterials 2019, 9, x FOR PEER REVIEW 9 of 16 active ions concentration. The emission spectra of Eu 3+ :HAp nanopowders are collected in Figure S9 and the emission spectra of PLLA/Eu 3+ :HAp composites are collected in Figure 7. The spectra were normalized to the 5 D0 → 7 F1 transition. As can be seen, the emission spectra of Eu 3+ ions are composed of five bands related to the 5   The ratio of the integral intensities of the 5 D0 → 7 F2 electric dipole transition to the 5 D0 → 7 F1 magnetic dipole transition is needed to evaluate the asymmetry of the coordination polyhedron of europium(III) ions and the variations in the local point symmetry. This is possible due to the fact that the intensity of the 5 D0 → 7 F2 transition is very sensitive to even small changes in the local environment of Eu 3+ ions in the crystal field, whereas the intensity of the 5 D0 → 7 F1 transition is nearly independent of those influences. If an Eu 3+ ion is located in a centrosymmetric site, the only permitted transition is the magnetic one. In the opposite case, the electric dipole transition is dominant. The ratio of the relative emission intensities (R) is defined by the equation: The higher the ratio between these transitions is, the less centrosymmetric the local environment around Eu 3+ ions becomes. The impact of the Eu 3+ ions concentration in powders and composites on the R values are presented in Table 2. It is difficult to observe a straight tendency in the case of both materials. With the influence of the Eu 3+ ions concentration, the R factor increases, so the Eu 3+ ions The ratio of the integral intensities of the 5 D 0 → 7 F 2 electric dipole transition to the 5 D 0 → 7 F 1 magnetic dipole transition is needed to evaluate the asymmetry of the coordination polyhedron of europium(III) ions and the variations in the local point symmetry. This is possible due to the fact that the intensity of the 5 D 0 → 7 F 2 transition is very sensitive to even small changes in the local environment of Eu 3+ ions in the crystal field, whereas the intensity of the 5 D 0 → 7 F 1 transition is nearly independent of those influences. If an Eu 3+ ion is located in a centrosymmetric site, the only permitted transition is the magnetic one. In the opposite case, the electric dipole transition is dominant. The ratio of the relative emission intensities (R) is defined by the equation: The higher the ratio between these transitions is, the less centrosymmetric the local environment around Eu 3+ ions becomes. The impact of the Eu 3+ ions concentration in powders and composites on the R values are presented in Table 2. It is difficult to observe a straight tendency in the case of both materials. With the influence of the Eu 3+ ions concentration, the R factor increases, so the Eu 3+ ions local environment became more distorted, but in the case of 5 mol% Eu 3+ :HAp, it decreased in both powders and in composites. Comparing both materials, the R value was higher in composites, which indicates that in these materials, the surroundings of the Eu 3+ ions is more distorted. Table 2. Decay rates of radiative (A rad ), non-radiative (A nrad ) and total (A tot ) processes of 5 D 0 → 7 F J transitions, luminescence lifetimes (τ), intensity parameters (Ω 2 , Ω 4 ), quantum efficiency (η) and asymmetry ratio (R) of the powders and composites. The simplified Judd-Ofelt approach was implemented to provide the intensity parameters Ω 2 and Ω 4 , as well as to have a deeper insight into the structure. The results of this approach are shown in Table 2. The value of the Ω 2 parameter indicates some changes in the distortion of europium coordination polyhedra caused by such factors as ions concentration, annealing temperature, etc., and could be related to an increase of the Eu 3+ − O 2− bond covalency. It is worth noting that the value of the Ω 2 parameter has the same tendency as the R factor in powders and composites. The value of the Ω 4 factor supplies some information about changes in the electron density around Eu 3+ cations. This value cannot be straightly correlated with changes in the Eu 3+ ions symmetry, but it can add some information about the electron density variations of the O 2− anions surrounding that influences the CT band position.
To determine the detailed spectroscopic properties of nanopowders and nanocomposites, the luminescence kinetics was analyzed. The luminescence life times were recorded at room temperature. The materials were excited by a 394.5 nm line and monitored at 618 nm corresponding to the 5 D 0 → 7 F 2 transition. The luminescence kinetic curves of the Eu 3+ :HAp nanopowders are shown in Figure S10, and for the Eu 3+ :Ca 10 (PO 4 ) 6 (OH) 2 nanoparticles embedded into poly(L-lactide) composites, they are shown in Figure 8. local environment became more distorted, but in the case of 5 mol% Eu 3+ :HAp, it decreased in both powders and in composites. Comparing both materials, the R value was higher in composites, which indicates that in these materials, the surroundings of the Eu 3+ ions is more distorted. The simplified Judd-Ofelt approach was implemented to provide the intensity parameters Ω2 and Ω4, as well as to have a deeper insight into the structure. The results of this approach are shown in Table 2. The value of the Ω2 parameter indicates some changes in the distortion of europium coordination polyhedra caused by such factors as ions concentration, annealing temperature, etc., and could be related to an increase of the Eu 3+ − O 2− bond covalency. It is worth noting that the value of the Ω2 parameter has the same tendency as the R factor in powders and composites. The value of the Ω4 factor supplies some information about changes in the electron density around Eu 3+ cations. This value cannot be straightly correlated with changes in the Eu 3+ ions symmetry, but it can add some information about the electron density variations of the O 2− anions surrounding that influences the CT band position. Table 2. Decay rates of radiative (Arad), non-radiative (Anrad) and total (Atot) processes of 5 D0 → 7 FJ transitions, luminescence lifetimes (τ), intensity parameters (Ω2, Ω4), quantum efficiency (η) and asymmetry ratio (R) of the powders and composites.  To determine the detailed spectroscopic properties of nanopowders and nanocomposites, the luminescence kinetics was analyzed. The luminescence life times were recorded at room temperature. The materials were excited by a 394.5 nm line and monitored at 618 nm corresponding to the 5 D0 → 7 F2 transition. The luminescence kinetic curves of the Eu 3+ :HAp nanopowders are shown in Figure  S10, and for the Eu 3+ :Ca10(PO4)6(OH)2 nanoparticles embedded into poly(L-lactide) composites, they are shown in Figure 8.

Tensile Properties
In Figure 9, we presented the tensile results for PLLA-based composites with 10 wt.% of HAp with different Eu 3+ ions concentration. The parameters are very important factors in designing 3D scaffolds for bone tissue engineering. From among three basic methods (polymerization in situ, solution-casting and melt mixing), we chose melt mixing-extrusion in a twin-screw extruder. This method ensures good filler distribution in the polymer. The preparation method is the most important in the case of nanocomposites preparation. Nanomaterials 2019, 9, x FOR PEER REVIEW 11 of 16

Tensile Properties
In Figure 9, we presented the tensile results for PLLA-based composites with 10 wt.% of HAp with different Eu 3+ ions concentration. The parameters are very important factors in designing 3D scaffolds for bone tissue engineering. From among three basic methods (polymerization in situ, solution-casting and melt mixing), we chose melt mixing-extrusion in a twin-screw extruder. This method ensures good filler distribution in the polymer. The preparation method is the most important in the case of nanocomposites preparation.

Thermal Properties
The strongest differences in thermal properties of PLLA and the composites concern cold crystallization. In the presence of HAp particles, the onset of cold crystallization (Tcc onset ), as well as that of peak temperature (Tcc), is about 3 °C lower compared to neat PLLA. This result indicates the well-known nucleating activity of HAp particles on PLLA cold crystallization [50]. It is worth underlining that the stronger effectiveness of nucleation activity towards cold crystallization is exhibited by HAp particles doped with Eu 3+ ions. The lowest Tcc onset was registered for the PLLA/3 mol% Eu 3+ :HAp composite, the cold crystallization onset is about 9.4 °C lower than for neat PLLA.
Moreover, the enthalpy of cold crystallization (ΔHcc) is lower in the case of the PLLA/Eu 3+ :HAp composite than for neat PLLA and PLLA with the unmodified HAp particles composite. The lower values of cold crystallization enthalpy prove the higher crystallinity of PLLA in the presence of Eu 3+ :HAp. The presence of Eu 3+ ions and its molar content have a strong influence on the crystallinity degree of PLLA in composites.
The lack of exothermic effects of the alpha'-alpha transition in the second heating DSC curves in the case of composites with Eu 3+ confirms that, upon subsequent cooling, PLLA crystallized in the alpha form. By contrast, neat PLLA and PLLA with unmodified HAp crystallized upon cooling as a

Thermal Properties
The strongest differences in thermal properties of PLLA and the composites concern cold crystallization. In the presence of HAp particles, the onset of cold crystallization (T cc onset ), as well as that of peak temperature (T cc ), is about 3 • C lower compared to neat PLLA. This result indicates the well-known nucleating activity of HAp particles on PLLA cold crystallization [50]. It is worth underlining that the stronger effectiveness of nucleation activity towards cold crystallization is exhibited by HAp particles doped with Eu 3+ ions. The lowest T cc onset was registered for the PLLA/3 mol% Eu 3+ :HAp composite, the cold crystallization onset is about 9.4 • C lower than for neat PLLA. Moreover, the enthalpy of cold crystallization (∆H cc ) is lower in the case of the PLLA/Eu 3+ :HAp composite than for neat PLLA and PLLA with the unmodified HAp particles composite. The lower values of cold crystallization enthalpy prove the higher crystallinity of PLLA in the presence of Eu 3+ :HAp. The presence of Eu 3+ ions and its molar content have a strong influence on the crystallinity degree of PLLA in composites.
The lack of exothermic effects of the alpha'-alpha transition in the second heating DSC curves in the case of composites with Eu 3+ confirms that, upon subsequent cooling, PLLA crystallized in the alpha form. By contrast, neat PLLA and PLLA with unmodified HAp crystallized upon cooling as a mixture of alpha and alpha' crystals, which is designated by the presence of an exothermic peak at 165 • C corresponding to the alpha'-alpha reorganization.
The PLLA crystallizes faster and in alpha and alpha' form in the presence of Eu 3+ ions embedded in the apatite. The HAp doped with Eu 3+ ions shown high nucleation activity towards cold crystallization as well as high crystallinity degree of PLLA in composites.

Spectroscopic Properties
The red emission was observed from the Eu 3+ :HAp powders and the PLLA/Eu 3+ :HAp composites. The analysis of the 5 D 0 → 7 F 0 transition can provide data about the number of crystallographic sites substituted by Eu 3+ ions into the host structure. Additionally, the existence of this transition approves that the Eu 3+ ions are located at a low-symmetry environment, and it is observed only if the Eu 3+ ions occupy sites with a local symmetry of C n , C nv or C s . The 5 D 0 → 7 F 0 transition is split into three components located at 573.5 nm (17 437 cm −1 ) and 577.2 nm (17 325 cm −1 ) and 578.6 nm (17 283 cm −1 ) in the case of Eu 3+ :HAp nanopowders. The Eu 3+ ions occupy three different types of crystallographic sites, which is well known in the literature on apatite host matrixes [51][52][53][54]. The 5 D 0 → 7 F 1 and 5 D 0 → 7 F 2 transitions consist of many overlapped Stark components. The most intense emission was observed for the hypersensitive 5 D 0 → 7 F 2 electronic transition. Its intensity is highly influenced by the local symmetry of the Eu 3+ ions as well as the type of ligands rather than the intensities of the other electronic transitions. This transition is notably used to calculate the asymmetry of the Eu 3+ site. In the case of europium(III) emission in the PLLA/1 mol% Eu 3+ :HAp composite, it is visible that the 5 D 0 → 7 F 0 transition is split into three components, but in the case of higher concentration, it is possible to distinguish only a single line. Surprisingly, the major difference is visible in the intensity of the 5 D 0 → 7 F 4 transition which is much more intense in the case of the composite in the entire range of optically active ions concentration.
The decay profiles are non-exponential. Since the physical meaning of multi-exponential fitting is complicated to explain, the lifetime values were calculated from the effective emission decay times using the following equation: where I(t) is the luminescence intensity at time t corrected for the background, and the integrals are calculated over the range of 0 < t < t max , where t max >> τ m . The value of decay time is shorter in the case of the Eu 3+ ions-doped hydroxyapatite embedded into poly(L-lactide). This observation could be connected with different response of the Eu 3+ ions present on the surface in nanopowders than in composites. The decay times shortened with an increase of Eu 3+ ions concentration in HAp embedded into PLLA.

Tensile Properties
The tensile strength of the PLLA/HAp composite is~40.1 MPa, which is comparable to all the three PLLA/Eu 3+ :HAp composites (40.4-44.3 MPa with a standard deviation 2.4-2.9 MPa). In this case, no influence of Eu 3+ ions on the tensile strength of the composites was observed. However, tensile strength of all the composites is lower compared to neat PLLA. The effect of mechanical properties reduction was described earlier [55] for PLLA doped with cellulose nanocrystals or carbon nanotubes but also for hydroxyapatite [16]. The reason for this effect is polymer-filler adhesion/interaction. In the research, neat calcium hydroxyapatite without any surface modification has been used. In this case, adhesion is rather poor and is not improved by the Eu 3+ ions incorporation. Another reason is the filler content. It has been shown before that, in the extruded composites, the higher the content of hydroxyapatite, the lower the tensile strength [16]. Young's modulus all the values for the composites are in the range of 3.1-3.2 GPa with a standard deviation of 0.1-0.3 GPa. When comparing the strain at break parameter, the reduction for PLLA/HAp (2.4%) as compared to neat PLLA (6.7%) can be observed along with further reduction of PLLA/Eu 3+ :HAp (1.4-1.7%) composites as compared to the PLLA/HAp system. The last effect is very slight (0.7-1.0%); however, the measurement error ranges from 0.1% to 0.2%. The tensile strength, Young's modulus and strain at break parameters value are the highest in the case of neat PLLA and are comparable to each other for all PLLA/HAp composites.

Conclusions
In the present research, pure as well as Eu 3+ -doped hydroxyapatite were synthesized by using the precipitation method and were confirmed by the XRD patterns analysis. The poly(L-lactide)/europium(III)-doped hydroxyapatite composite foils were successfully fabricated using the twin-screw co-rotating micro-extrusion technique in the weight ratio 9:1. The obtained composites were found to have similar properties to neat PLLA. Differential Scanning Calorimetry revealed some differences in thermal properties of PLLA and the composites. In the presence of HAp particles, the cold crystallization is lower in comparison with neat PLLA. This result indicated a well-known nucleating activity of HAp particles towards the cold crystallization. Moreover, stronger effectiveness of nucleation activity towards cold crystallization was shown in the case of HAp particles doped with Eu 3+ ions. The HAp nanopowders and composites doped with Eu 3+ ions are characterized by red emission under UV radiation with the 5 D 0 → 7 F 2 transition as the most intense. The luminescence kinetics shows a non-exponential curve indicating that more than one emitting center is present in the investigated materials. The value of luminescence decays is lower in the case of the Eu 3+ ions-doped hydroxyapatite embedded into poly(L-lactide).
The PLLA/Eu 3+ :HAp composites were obtained as prospective candidates to theranostic applications (therapy and diagnostics) due to support of bone healing by hydroxyapatite and bio-imaging possibility of Eu 3+ ions.