Wettability of soft PLGA surfaces predicted by experimentally augmented atomistic models

A challenging topic in surface engineering is predicting the wetting properties of soft interfaces with different liquids. However, a robust computational protocol suitable for predicting wettability with molecular precision is still lacking. In this article, we propose a workflow based on molecular dynamics simulations to predict the wettability of polymer surfaces and test it against the experimental contact angle of several polar and nonpolar liquids, namely water, formamide, toluene, and hexane. The specific case study addressed here focuses on a poly(lactic-co-glycolic acid) (PLGA) flat surface, but the proposed experimental-modeling protocol may have broader fields of application. The structural properties of PLGA slabs have been modeled on the surface roughness determined with microscopy measurements, while the computed surface tensions and contact angles were validated against standardized characterization tests, reaching a discrepancy of less than 3% in the case of water. Overall, this work represents the initial step toward an integrated multiscale framework for predicting the wettability of more complex soft interfaces, which will eventually take into account the effect of surface topology at higher scales and synergically be employed with experimental characterization techniques. Controlling the wettability of surfaces has important implications for energy (e.g., self-cleaning solar panels), mechanical (e.g., enhanced heat transfer), chemical (e.g., fluids separation), and biomedical (e.g., implants biocompatibility) industries. Wetting properties arise from a combination of chemical and physical features of surfaces, which are inherently intertwined and multiscale. Therefore, tailoring wettability to target functionalities is a time-intensive process, especially if relying on a trial-and-error approach only. This becomes even more challenging with soft materials, since their surface configuration depends on the solid-liquid interactions at the molecular level and could not be defined a priori. The improved accuracy of atomistic models allows detailing how the effective properties of materials arise from their nanoscale features. In this article, we propose and validate a new molecular dynamics protocol for assessing the wettability of soft interfaces with polar and nonpolar liquids. The prediction capabilities of simulations are augmented by a close comparison with microscopy and contact angle experiments. Since smooth copolymer surfaces are considered, here the effort mainly focuses on the effect of chemical features on wettability. In perspective, the proposed atomistic in silico approach could be coupled with computational models at higher scales to include the effect of surface microstructures, eventually easing the development of multi-scale surfaces with tunable wettability.


Introduction
The study of soft interfaces is an attractive topic for scientists and engineers due to the broad range of industrial applications dependent on surface properties, such as energy-saving films, 1 antibacterial coatings, 2 biodegradable scaffolding for tissue regeneration, 3,4 rheological and mechanical control, 5,6 phase-separation membranes and interfaces, [7][8][9] and polymer sensors. 10 From an experimental standpoint, different techniques have been used to characterize soft surfaces, giving an insight into their chemical, physical, structural, and topological properties. Spectroscopy analysis is commonly used to highlight the molecular composition, 11,12 electronic microscopy, 13 and scanning probe microscopy 14 to analyze the topology, whereas contact angle (CA) measurement 15 to quantify the surface free energy (SFE, using Girifalco or Owens-Wendt models 16,17 ). Moreover, recent developments in contact-mechanics-based protocols have been proposed as alternative means to investigate SFE at the microscale via the Johnson-Kendall-Roberts (JKR) or Derjaguin-Muller-Toporov (DMT) models. 18 However, the principal limitation of experimental characterization alone is the lengthy trial-and-error approach, which involves looping over the cycle of materials design, production, and analysis, thus extending the time of prototyping and, consequently, commercialization. 19 In this scenario, computational experiments can help predict the role played by the chemical and topological characteristics on the effective properties of soft interfaces across multiple scales, therefore providing a rational design of their composition and structure according to the required target performance. 20, 21 Given the time and space scales involved in wetting phenomena, 22 the adopted computational methods must include a combination of models from the continuum to the atomistic scale.
Since the pioneering study of Young, 23 several continuum models have been proposed to describe the dependence of the static contact angle and dynamical wetting process with surface morphology. Generally, the CA of a liquid on a surface can be correlated to the surface tension components solid/ vapor ( γ sg ), solid/liquid ( γ sl ), and liquid/vapor ( γ lg ), through the Young equation: 24 where θ Y is the Young CA, which depends only on the chemistry of the liquid and solid phases. Despite its simplicity, this relation is valid only for ideal flat surfaces. The CA obtained experimentally corresponds to the apparent CA, which depends also on the surface roughness and morphology. According to the Wenzel and Cassie-Baxter models of static contact angle on rough surfaces, [24][25][26] the Young or local contact angle can be included in continuum models as a boundary condition on the wetted wall. This value depends on the surface chemistry of the material, and it can be provided by experiments performed on a perfectly flat surface or-as highlighted in the work of Johansson et al. 27 29,30 to evaluate the CA of water on graphene and, consequently, to calibrate their interaction parameters. Similarly, Zhang et al. 31 employed the FEP to parametrize the interaction of water with fresh and aged molybdenum disulfide (MoS 2 ) surfaces, whereas Rajan et al. 32 used the FEP and the Green-Kubo relation to evaluate the work of adhesion and friction coefficient of water and other solvents on a hexagonal boron nitride surface. In our previous work, 33 we highlighted that a synergistic use of MD simulations and experiments allows discriminating how surface chemistry and topology 1 γ sg = γ sl + γ lg cos θ Y , of a hydrophobic coating contribute to increasing the contact angle of water. The CA assessment by MD can be carried out also via the droplet method, where a liquid drop is positioned above the tested surface and then left to relax until equilibrium conditions are reached. Using this method, Yaghoubi et al. and  studied the effect of roughness on the wettability of graphite surfaces, whereas Sresht et al. assessed MoS 2 surfaces. 37 In such simulations, the solid substrate is typically modeled as a rigid crystal, thus neglecting deformation effects at the solid-liquid interface. Similarly, Liu et al., 38 Zhu et al., 39 and Etha et al. 40 applied the droplet method to soft surfaces, keeping the substrate rigid to preserve the planar configuration of the interface. However, in the case of soft surfaces such as for polymers, the interfacial deformation of the solid substrate and its effect on wettability cannot be neglected to accurately reproduce the actual experimental conditions. Hence, the characteristic deformable boundary of soft interfaces requires new (and validated) 3D simulations and post-processing protocols for the CA measurement by MD, especially for the definition of the mean plane where the contact angle is evaluated.
In this work, we propose and validate a new MD protocol, based on the sessile droplet method, for assessing the wettability of smooth, soft interfaces with polar and nonpolar liquids. In detail, we: 1. Validate the force field by computing the density and glass-transition temperature of a soft material (a polymer), the density and surface tension of the probe liquids, and comparing the computed values with the experimental data available. 2. Create a MD sample whose surface matches the experimental topology of the soft material. 3. Simulate the 3D sessile droplet experiment by evaluating the apparent contact angle from the MD trajectory.
The proposed protocol is employed to study the wettability of poly(lactic-co-glycolic acid) (PLGA) surfaces with water, formamide, toluene, and hexane ( Figure 1). PLGA is a copolymer formed by repeating units of lactic acid (LA) and glycolic acid (GA). PLGA is a relevant case study because its wetting properties determine its biocompatibility and biodegradability in many biomedical applications. [41][42][43] In perspective, this validated protocol can be easily transferred to study flat and rough surfaces of both soft and hard materials, with a promising synergistic potential with respect to experimental characterization techniques.

Results and discussion
The characterization protocol discussed in this work employs MD simulations to predict the wetting properties of soft interfaces. The details related to the model building process are reported in the "Materials and methods" section. The wetting properties of PLGA surfaces are assessed after carefully evaluating the bulk properties of the polymer slab and the probe liquids. First, we evaluate the density WEttabiLity OF sOFt PLGa sUrFacEs PrEdictEd by ExPEriMEntaLLy aUGMEntEd atOMistic MOdELs and glass-transition temperature (T g ) of the simulated PLGA bulk. The density of PLGA computed at 300 K during the equilibrated MD trajectory is equal to 1.275 ± 0.025 g/cm 3 , in excellent agreement with the experimental value (1.30 g/cm 3 ) 44 . The T g was evaluated by performing a scan in temperature and monitoring the variation of density. 45 In Figure 2a, we report the density of PLGA against temperature, identifying the two distinct branches of the plastic and glassy states. These branches are interpolated with two linear regressions, intercepting at a temperature T g = 328 K that accurately matches the experimental value of 320-330 K. 44 Similarly, we evaluate the density of the probe liquids at 300 K and 1 atm, finding an excellent agreement with the experimental evidence from the literature (see Table SI in the Supplementary information). For the specific case of water in contact with the polymer substrate, we investigate the coupling of the two phases by considering a liquid film onto the polymer slab and evaluating the density profile along the normal direction to the polymer surface. 46 Results in Figure 2b show an overlap between the two phases in the contact region due to the soft (polymer/water) interface, while the nominal densities of PLGA and pure water are eventually recovered after the interface. In Figure 5a, we report the calculated surface tension for the four liquids, finding that MD results are in good agreement with both experimental and other numerical results available in the literature [47][48][49][50][51] (see Table SI in the Supplementary information for details). This evidence supports the reliability of the force field used to describe the PLGA and probe solvents, thus providing confidence for predicting the wetting properties of surfaces. Second, we characterize the topography of the PLGA surface in terms of the root-mean-square (rms) roughness of the MD sample, comparing this value with that obtained by atomic force microscopy (AFM) measurements. The smallest viable experimental scans are of 100 nm × 100 nm, from which representative zoomed sections were extracted with a lateral dimension comparable to that of the MD samples (see Figure 3). The rms roughness of the numerical and experimental samples is 3.2 Å and 2.2 Å, respectively. According to the Wenzel model for hydrophilic surfaces 25,52,53 (i.e., A proj the roughness factor, being A eff the effective area of the polymer surface and A proj the projected area on the acquisition window), we found the experimental roughness factor in the range of 1.0570-1.1078 depending on the scan size (see Table SII in the Supplementary information) and 1.5250 in the MD sample. Considering the nonperiodicity of the PLGA surface features (as analyzed from 2D fast Fourier transform of AFM scans), the influence of surface topology on the measured CA is smaller than 5% for the conducted experimental tests, thus within the statistical uncertainty of the measure. Consequently, the CA obtained from experiments and simulations can be considered as the Young CA with good approximation ( θ ≈ θ Y ). The only notable exception is the case of formamide on PLGA, where the low value of θ leads to larger deviations from θ Y (see Table  SII in the Supplementary information).
Third, the CA of the four liquids on the flat PLGA surface are computed and compared against the experimental measurements. In the case of water droplets, the CA of the SPC-E and the Tip4P-05 models was found to be 61.7 • ± 3.8 • and 64.6 • ± 5.0 • , respectively, in excellent agreement with the experimental measure of 65 • ± 4 • (see Figure S1 in the Supplementary information) and previous works from the literature. 54 The hydrophilic nature of the PLGA sample is also highlighted considering the surface free energy of 43 ± 3 mJ/m 2 evaluated through the Owen-Wendt method 16 (see Note 1 in the Supplementary information). The evolution of a water droplet onto the PLGA surface is shown in Figure   4a, where equilibrium conditions are reached after some 500 ps of simulation. The CA is measured during the last 2 ns of trajectory, where the CA oscillates around a stable mean value (see Figure 4b, Tip4P-05 model). The error on the CA is computed by propagating the error on each snapshot due to different choices of δ (see the "Materials and methods" section). As reported in Figure S2 of the Supplementary information, the average CA obtained with the SPC-E water model and the standard deviation of each snapshot is similar to the Tip4P-05 model, proving the reliability of our numerical protocol regardless of the chosen water model. In the case of formamide, we observe a mismatch between the computed (27.3 • ± 5.5 • ) and experimental (52 • ± 4 • , see Figure S3 in the Supplementary information) CA. Considering the Wenzel model, however, the measured apparent CA (27.3 • ) can be used to extract the Young CA (54.3 • ), which results in excellent agreement with the experimental value of θ Y (see Table SII in the Supplementary information). Furthermore, given the good solubility of formamide in water, we can assume that atmospheric moisture in the testing environment could be absorbed in the formamide droplet during the experimental measurements. However, the perfectly anhydrous conditions in MD simulations can partially also explain the mismatch between numerical and experimental values of CA for formamide. We tested this hypothesis by performing additional simulations with droplets of formamide containing 10% wt and 20% wt of water, using the same computational procedure (see Figure S4-S6 in the Supplementary information). In the first case, we did not observe a significant variation in the evaluated CA, whereas in the second one, the CA reached its equilibrium at a higher CA value, namely 32 • ± 9 • .
In the case of nonpolar hydrocarbon liquids, the simulations reproduce the high affinity with the PLGA surface observed experimentally with good accuracy. For toluene, MD shows an initial metastable state of the droplet with a CA of 11 • ± 2 • for approximately 500 ps before a complete spreading is eventually achieved (see Figure S7 in the Supplementary information). During experiments, we observe a b c a stable contact angle of 23 • ± 5 • for approximately 15 s before a complete spreading (see Figure S8 in the Supplementary information). Indeed, it is worth noting that toluene is a good solvent of PLGA and that the size of the MD sample is orders of magnitude smaller than the experimental one. Hence a shorter time scale is expected for the interaction between the two phases. In the case of hexane, we observe a complete spreading of the droplet onto the PLGA surface (0  Figure 5b. The main advantage of the numerical method proposed and validated here is to simplify and accelerate the initial stages of material surface design and prototyping, offering an integrated framework that allows predicting the wettability of a wide range of soft interfaces with good accuracy. However, MD simulations can suffer from artefacts due to finite-size effects. 55 The size of the samples investigated in this work was capped at approximately 350,000 atoms to match the computational resources available. These constraints limit the size of the surface investigated, the maximum roughness that can be realistically modeled on a periodic surface, and of course, the size of the droplet, whose diameter should not exceed half the periodic box vectors. An elegant way of increasing the size of the droplets while keeping the computational cost at a manageable level is by adopting a pseudo-2D model, where the droplet is replaced by a periodic cylinder. 27,28 Another significant limitation is the time scale accessible through MD simulations, here on the order of 10 ns per sample. With the parameters used in our MD simulations, most samples reach thermodynamic equilibrium within a few nanoseconds unless other phenomena with more long time scales are involved. [56][57][58] Although it cannot be ruled out that the observed droplets were in a metastable state lasting longer than the simulated time, we have verified that the contact angle reaches a consistent equilibrium value for droplets of increasing size. In this work, we have used droplets containing approximately 13,000 molecules to limit finite-size effects and reduce the numerical noise of simulations.
A way to overcome the limitations of atomistic models is to decrease the resolution of the system and its internal degrees of freedom, a process known as coarse-graining. A recently developed coarse-grained model, which reproduces the morphology and density of organic semiconductors 59 and polymers, 60 and the structural and dynamic properties of water 61 in quantitative agreement with atomistic MD simulations, will be used in a follow-up study to predict the wetting behavior of soft interfaces, addressing the known limitations of atomistic models but without compromising on the quality of the predicted properties.
This considered, the evaluation of water droplet sizedependent behavior on hydrophilic surfaces from literature data 62 reported that a remarkable effect of surface asperities, within the microstate pinning behavior, is such that a narrowing of contact angle evaluated error from MD simulations occurs for a nanoscale roughened surface. This is in accordance with authors' preliminary evaluations for specific selection of the simulated sessile droplet volume (while roughness itself determines an increase in the CA mean value). As, however, accordance with experimental data are stated throughout this work, it is worth highlighting that our modeling activities reached the pinning state equilibrium. With this goal in mind, this work represents the first step toward the creation of a multiscale model able to predict and optimize the wettability of microstructured soft surfaces, for instance, by providing fast and accurate values of the local or Young contact angle, which in turn can be used as a boundary condition in multiphase continuum models (e.g., finiteelements or finite-volume simulations).

Conclusions
In this work, we reproduced sessile droplet experiments with molecular dynamics simulations, proposing a methodological procedure for predicting the wettability of soft interfaces. In detail, we studied PLGA surfaces wetted by polar and nonpolar solvents (water, formamide, toluene, and hexane) and validated the obtained results with experiments.
We initially evaluated the surface topography of the MD and experimental samples through the solvent accessible surface and AFM, respectively. Both numerical and experimental surfaces present a rms roughness smaller than 10 Å, while the roughness factor (R f ) is between 1.0570 and 1.1078 in the experimental samples and 1.5250 in the MD sample. According to the Wenzel theory, the local or Young contact angle can be assessed with droplet experiments under these conditions since negligible effects from surface topology are expected. For the numerical model, the surface topology may affect the measurement of liquids that present a low apparent contact angle, such as formamide. The glass-transition temperature of bulk PLGA and the surface tension of the probe liquids computed from MD simulations were found to be in excellent agreement with the experimental evidence, thus validating the force field chosen to describe the materials in this study. Finally, we performed the sessile droplet experiment, leaving a droplet of solvent to relax onto the PLGA surface and then measuring the apparent contact angle by post-processing the MD trajectory. The results obtained agree with the experimental measures, therefore supporting the proposed protocol. In the case of water droplets, both the SPC-E and the Tip4P-05 models show an excellent agreement with experiments. With the formamide model, we observe a reduction of approximately 30% with respect to the experimental CA: we attributed this mismatch to the effects of the roughness factor (R f ) according to the Wenzel theory and to the perfectly anhydrous conditions in the MD simulations of formamide droplets, differently from the testing environment. The MD simulations predict a strong affinity of nonpolar solvents with the polymer surface, which was confirmed by the experimental measurements. Toluene spreads over the PLGA surface after a transient metastable state, whereas the hexane droplet spreads immediately.
This work is a relevant case study of convergence and synergy between digital and experimental measurements toward a new class of hybrid characterization standards (e.g., coupled MODA 63 and CHADA 64 workflows). In perspective, the proposed methodology will speed up the development of soft material interfaces with controllable wetting properties. The protocol presented in this work has been shown to accurately predict the CA of polar and nonpolar solvents on a soft polymeric surface. This information can be used as a boundary condition in continuum simulations able to assess the effect on wettability produced by modifications of the surface topology at the nanoand microscales, thus connecting the property of materials at the device level with their underlying molecular structure.

Model building
The experimental length of PLGA chains typically comprises between 240 and 390 repeating units, a size that is computationally inaccessible for practical applications based on atomistic simulations. Literature works on PLGA have reported MD simulations with a polymer length between 10 and 90 monomers. [65][66][67][68] In this work, we used PLGA chains composed of 48 monomers, corresponding to 12 repetitions of 3 LA monomers and 1 of GA monomer, thus keeping the 75:25 LA to GA ratio of the experimental samples.
The structure of a single PLGA chain was built with Avogadro, 69,70 whereas the condensed PLGA phase was generated starting from a randomly sparse configuration of chains in a simulation box with periodic boundary conditions in the three orthogonal directions (see Figure 1a-b). An initial guess composed of 90 polymer chains in a cubic box with side of 900 Å was created with Moltemplate, 71 whereas the MD simulations were carried out with LAMMPS. 72,73 The OPLS-AA force field 74 was used to describe the interaction potentials of PLGA. The wettability of PLGA was evaluated with four liquids: water (both Tip4P-05 75 and SPC-E 76 models), formamide, toluene, and hexane, which were modeled with the OPLS-AA force field (Figure 1c-f). A cutoff of 1.3 nm was used for nonbonded interactions, with cross-term Lennard-Jones parameters computed with geometric mixing rules. Long-range interactions were computed with the PPPM method with 10 −4 accuracy.
The equilibration of PLGA samples began with a controlled multistep compression at 1 K, in which the volume of the simulation box was progressively reduced by 15% in each step for 49 steps. The electrostatic interactions were switched off during this process to emulate the screening induced by a virtual solvent. Then, the bulk sample was relaxed in the isobaric-isothermal ensemble (NPT) at 300 K and 1 atm for 2 ns, with nominal electrostatic interactions. Surfaces of PLGA were created by removing the periodic boundary in the z-direction, turning the bulk into a slab geometry. The system was then equilibrated in the canonical ensemble (NVT) at 300 K to allow surface reconstruction for approximately 4 ns. A correction term for systems with slab geometry was added to the long-range electrostatic interactions. 77 The equilibrated surface was then replicated three times in the x and y directions, yielding a square surface with lateral dimension * of approximately 220 Å. This supercell was used for CA simulations, as its larger size helped to avoid finite-size effects, which could affect the wetting properties of the soft interface. 34,35 The topological characterization of PLGA surfaces is based on the solvent accessible surface (SAS) of the slab, computed with the program Jmol 78 on the final configuration of the NVT simulation. The resulting SAS was exported as a triangular mesh and edited with Meshlab 79 to remove the boundary walls. Then, the surface roughness of PLGA was computed with Gwyddion, 80 taking the meshed SAS as input, and compared with the experimentally determined one.
The density and surface tension of the four liquids were used to gauge the accuracy of the force field, by comparing the computed values with the experimental ones. The density was measured on samples made of 1000 molecules in a periodic box reproducing bulk conditions. Each sample was equilibrated in the NPT ensemble at 300 K and 1 atm for 2 ns with coupling constants of 0.1 and 1 ps for temperature and pressure, respectively. The surface tension was calculated on a slab geometry according to the Kirkwood-Buff relation: 81 where P z is the liquid pressure in the direction orthogonal to the surface, P x and P y are the pressures in the transverse directions, and L z is the length of the simulation box in the orthogonal direction to the surface, whereas brackets indicate an NVT ensemble average. A production run of 5 ns was used for data acquisition after the slab was equilibrated for 5 ns.
The apparent contact angle of the different liquids on the PLGA slab was computed with numerical experiments following the droplet method. The initial guess of each droplet was a cube of equilibrated liquid with box edges of approximately 73 Å (see Figure 6a). The droplet was equilibrated on top of the PLGA slab until the potential energy and the evaluated CA reached a stable value with time. The total number of atoms of the slab/solvent systems was approximately 350'000.
The contact angle was calculated from the MD trajectory using a post-processing code developed for this purpose. We initially define a mean plane representative of the interface, comparing the positions of liquid and polymer surface atoms. The droplet axis was identified from the average position of all the droplet atoms, and the droplet height h from the interface plane. Then, we consider two contiguous liquid slabs ( s 1 and s 2 ) of height δ at a distance δ from the interface plane. We evaluate the radial density profile of the liquid in s 1 and s 2 to estimate the droplet radius r 1 (at z 1 , that is the mean z coordinate of s 1 ) and r 2 (at z 2 , that is the mean z coordinate of s 2 ) in the x − y plane. We measure the radius at the liquid-air interface, defined as the point where the density falls in a range below 50% of the bulk value (see the code available at Reference 82 for details). Finally, the CA is evaluated as θ = arctan(�z/�r) , where r = r 1 − r 2 and z = z 2 − z 1 (see Figure 6b). This procedure is subject to the choice of δ for the evaluation of r 1 and r 2 . Then, * The lateral size of the supercell is approximately two orders of magnitude bigger than the experimental surface roughness.
we replicate the procedure with different values of δ , expressed as a percentage of the drop height (within 7% and 10%), and average the CA final value from these results.

Sample preparation
Typical commercial proportions of LA:GA in PLGA samples range from 50:50 to 85:15, 50:50 and 75:25 being the most common ones. 44 The presence of CH 3 side groups makes LA more hydrophobic than GA. Therefore, a copolymer with higher LA content absorbs less water and then degrades more slowly, the 50:50 ratio being the one showing the fastest degradation rate. 41,83 In this work, we focus on 75:25 copolymers (i.e., 66'000-107'000 g/mol molecular weight) to achieve a good balance between durability and hydrophilicity. The PLGA solution was prepared in a mixture of chloroform (Fisher Chemical, USA) with a concentration of 5% in weight. The solution rested for a homogenization time of 12 h. Then, 2 ml of solution was poured into a dry petri dish with a diameter of 12 cm, and finally dried until the solvent was fully evaporated and the flat slab of PLGA obtained.

Experimental characterization
The topological characterization of the PLGA samples was performed in terms of the surface root mean square roughness, evaluated using a Bruker atomic force microscope (AFM dimension icon, Bruker) operating in standard tapping mode at a scan rate of 0.996 Hz. The AFM is equipped with a TESP-SS Bruker silicon probe featuring a nominal cantilever elastic constant of 42 N/m and a tip with a nominal radius of 2 nm. The selection of the AFM operating mode and tip radius value is based and cross-validated considering the relative errors (from the true value of rms roughness theoretically approachable using an infinitely sharp tip) arising from eventual convolution effects on surface topography, as elucidated by Sedin et al. 84 Indeed, for small lateral scans (≤500 nm), the image root mean square roughness decreases as tip size increases: in the case of the investigations performed, it has been evaluated an error for a tip radius of 2 nm lesser than 3%, eventually falling within the experimental error associated with surface morphology repeatability. Each measurement was performed on a scan area of 100 nm × 100 nm or 500 nm × 500 nm, where two sensor images were recorded simultaneously (corresponding to forward and backward scans). Images were then analyzed and processed with the software Gwyddion. 80 The solid-liquid interaction between PLGA and the different fluids was validated by sessile droplet contact angle measurements, employing normative-compliant equipment (UNI EN 828, UNI 9752, ASTM D-5725-99). CA testing was conducted at room temperature (20 • C) and relative humidity between 40% and 50 percent. A minimum of six droplets (with a volume of 3 µ L each) of water, formamide (polar solvents), toluene, or hexane (nonpolar hydrocarbons) were deposited on the sample surface to measure the apparent contact angle and then estimate the overall surface free energy (SFE) of the PLGA samples. After 60 s of stabilization, images of the liquid drops were captured (LAS-EZ software, Leica Camera AG) and subsequently elaborated (ImageJ software) using a spherical approximation of the droplet shape. This method allows assessing the angle formed between the two tangents of the fitted sphere with the straight line corresponding to the sample surface, thus obtaining the CA value ( θ). WEttabiLity OF sOFt PLGa sUrFacEs PrEdictEd by ExPEriMEntaLLy aUGMEntEd atOMistic MOdELs carried out at the Inter-Departmental Laboratory of Electron Microscopy (LIME), Università degli studi Roma Tre (http:// www. lime. uniro ma3. it). The Grant of Excellence Departments, MIUR (ARTICOLO 1, COMMI 314 -337 LEGGE 232/2016), to the Department of Engineering, Università degli studi Roma Tre, is also gratefully acknowledged.

Funding
Open access funding provided by Politecnico di Torino within the CRUI-CARE Agreement.

Data availability
All data generated or analyzed during this study are included in this published article as supporting information files. All simulation files and post-processing codes are available at the Zenodo archive associated with this work (https:// doi. org/ 10. 5281/ zenodo. 66294 27).

Conflict of interest
On behalf of all authors, the corresponding authors state that there is no conflict of interest.

Open access
This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.