Next Article in Journal
Polyaniline-Supported Nickel Oxide Flower for Efficient Nitrite Electrochemical Detection in Water
Next Article in Special Issue
Molecular Weight Segregation and Thermal Conductivity of Polydisperse Wax–Graphene Nanocomposites
Previous Article in Journal
Cd(II)/Mn(II)/Co(II)/Ni(II)/Zn(II) Coordination Polymers Built from Dicarboxylic Acid/Tetracarboxylic Acid Ligands: Their Structural Diversity and Fluorescence Properties
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Modelling across Multiple Scales to Design Biopolymer Membranes for Sustainable Gas Separations: 1—Atomistic Approach

by
Kseniya Papchenko
1,
Eleonora Ricci
2,† and
Maria Grazia De Angelis
1,3,*
1
Institute for Materials and Processes, School of Engineering, University of Edinburgh, Sanderson Building, Robert Stevenson Road, Edinburgh EH9 3FB, UK
2
Department of Civil, Chemical Environmental and Materials Engineering, DICAM, University of Bologna, Via Terracini 28, 40131 Bologna, Italy
3
National Interuniversity Consortium of Materials Science and Technology INSTM, Via G. Giusti, 58100 Firenze, Italy
*
Author to whom correspondence should be addressed.
Current address: National Centre for Scientific Research “Demokritos”, Agia Paraskevi, 15341 Athens, Greece.
Polymers 2023, 15(7), 1805; https://doi.org/10.3390/polym15071805
Submission received: 27 February 2023 / Revised: 28 March 2023 / Accepted: 4 April 2023 / Published: 6 April 2023

Abstract

:
In this work, we assessed the CO2 and CH4 sorption and transport in copolymers of 3-hydroxybutyrate and 3-hydroxyvalerate (PHBV), which showed good CO2 capture potential in our previous papers, thanks to their good solubility–selectivity, and are potential biodegradable alternatives to standard membrane-separation materials. Experimental tests were carried out on a commercial material containing 8% of 3-hydroxyvalerate (HV), while molecular modelling was used to screen the performance of the copolymers across the entire composition range by simulating structures with 0%, 8%, 60%, and 100% HV, with the aim to provide a guide for the selection of the membrane material. The polymers were simulated using molecular dynamics (MD) models and validated against experimental density, solubility parameters, and X-ray diffraction. The CO2/CH4 solubility–selectivity predicted by the Widom insertion method is in good agreement with experimental data, while the diffusivity–selectivity obtained via mean square displacement is somewhat overestimated. Overall, simulations indicate promising behaviour for the homopolymer containing 100% of HV. In part 2 of this series of papers, we will investigate the same biomaterials using a macroscopic model for polymers and compare the accuracy and performance of the two approaches.

1. Introduction

Membrane-based separation of gaseous mixtures is recognised as a low-energy and low-carbon alternative to traditional separation techniques [1,2,3,4]. The use of membranes can reduce the overall footprint of chemical processes, but most polymeric materials used for this purpose cannot be considered sustainable since they are fossil-based. Most industrial applications require a full replacement of the membrane module every three to five years [5,6], and recent works focused on finding renewable alternatives for membrane materials, such as bio-based and biodegradable polymers [2,7,8].
Polymers of the poly(hydroxyalkanoate)s (PHAs) family are attractive, as their origin is fully renewable, and biodegradability can be achieved in many environments [9,10]. Indeed, PHAs are linear polyesters synthesised naturally by several microorganisms as an energy reserve, while also well-established chemical routes for controlled PHA synthesis were developed [9,11,12,13]. Poly(3-hydroxybutyrate) (PHB), poly(3-hydroxyvalerate) (PHV), and their copolymers are currently the most studied and produced bio-polyesters, as the application range of such materials appears to be quite broad [14,15,16]. At the same time, the experimental characterisation of fluid transport properties, to this day, was reported only for a low number of gases in a small subset of PHAs [17,18,19].
Despite the attractiveness of such materials, several issues might delay their commercialisation, such as the difficulty in finding an optimal formulation which has good overall properties. On the other hand, a systematic experimental study is expensive and time-consuming, as there are many possible candidates. Computational approaches that can provide structure-property relationships faster need to be considered to reduce the time-to-market of biopolymers like the ones inspected here.
In this work, we focus on transport properties, specifically gas solubility, diffusivity, and permeability, for which different modelling techniques have been proposed in the literature and are readily-available to use [20]. In particular, molecular dynamics (MD) simulations have been proven efficient in predicting a range of polymer properties [21,22,23]. Recently, different groups applied MD to study new possible applications for PHAs [24] and to successfully investigate the thermal and mechanical properties of such materials [25,26,27].
In this work, an analogous approach was used to explore the gas transport in PHAs, with the aim of using them as sustainable membranes for gas separation. MD simulations were conducted to predict the CO2/CH4 separation performance of PHBV copolymers with different compositions, which is of interest for CO2 capture, biomethane purification and negative-emission processes. The computational screening of different structures to identify the most promising ones allows for the reduction of tests and guides future experimental efforts. The polymer models were validated using data coming from the literature [28,29,30,31,32,33], while for the gas-polymer mixtures, specific gas-sorption tests were conducted in this work.

2. Materials and Methods

2.1. Experimental

Materials

Bacterial poly(3-hydroxybutyrate-co-3-hydroxyvalerate) (PHBV), 8 mol% hydroxyvalerate (HV) units, was purchased from Sigma-Aldrich (St. Louis, MO, USA) in the form of film and used as received. The chemical structure of the biopolymer used in this study is shown in Figure 1. The thickness was measured to be 19 ± 1 µm. The crystallinity was not provided by the producers and was determined to be equal to 42 ± 1% through DSC analysis, as described in Section S1 of the Supporting Information (SI) [19]. Gases used for sorption tests, namely methane (CH4) and carbon dioxide (CO2), were purchased from Fluido Tecnica (Campi Bisenzio, FI, Italy) with purities ≥ 99.5%. The density of the semi-crystalline membrane was assumed equal to 1.214 g·cm−3, as reported by Mitomo et al. for an 8% HV copolymer [28].

2.2. Gas Solubility and Diffusivity Measurements

Solubility and diffusivity of CH4 and CO2 in PHBV were determined at 35 °C at pressures up to 10 bar in a pressure-decay apparatus [34]. The measurements were repeated at least twice for each gas. The full equipment set-up is well-represented and described elsewhere [35,36]. In this technique, a known amount of gas is fed into the sample chamber and the mass uptake is evaluated by measuring the pressure decrease over time. At equilibrium, the pressure becomes constant, as does the amount of gas absorbed by the membrane; after such values are recorded, a new pressure is imposed, and the procedure is repeated in a similar manner.
Sorption isotherms are constructed by connecting the equilibrium concentration values of the penetrant molecule in the polymer at each experimental pressure. The solubility coefficient S is the ratio between the concentration of the gas in the sample, c , and the pressure, p , as follows:
S = c p
Dense homogeneous matrices at temperatures above the glass transition, T g , such as the ones inspected here, generally follow Henry’s law, which describes a linear correlation between the gas uptake and the equilibrium pressure, so that S is constant with pressure. In each sorption step, the penetrant diffusivity in the film, D , can be evaluated from the sorption kinetics by considering Fickian diffusion and the variation of interfacial concentration during the experiment, as reported in the literature [35,37,38] and shown in Section S2 of the SI.
The performance of the membrane in gas separation is usually evaluated by the permeability, P , proportional to the molar flux of the gas across the material, and the selectivity, α , which is equal to the ratio between gas permeabilities. Under the assumption of the validity of the solution diffusion model [39], the permeability coefficient can be split into the product of the diffusion coefficient, D i , and the solubility coefficient, S i , as follows:
P i = D i S i
Under the same assumptions, the ideal selectivity between gas i and gas j, α i j , becomes the ratio of pure-gas permeabilities, and can be split into the product of the ideal diffusivity–selectivity, α i j D , and the ideal solubility–selectivity, α i j S :
α i j = P i P j = D i S i D j S j = α i j D α i j S

2.3. MD Simulations

Eight different PHBV copolymers were considered for density and solubility parameter estimation with 0, 8, 16, 24, 40, 60, 80, and 100 mol% of HV units, which in the following will be referred to as PHBV0, PHBV8, …, and PHBV100.
Monodispersed, amorphous polymeric melts at different compositions were generated using the Amorphous Builder plugin of the Materials and Processes Simulations (MAPS) software [40,41]. All systems were generated at 600 K. First, an energy minimization and short equilibration in the NVT (isothermal) and NPT (isobaric–isothermal) ensemble at 600 K was performed. Such equilibration allows to eliminate the overlaps that could have been introduced during the initial configuration generation and obtain a homogeneous density throughout the system. Then, a decreasing temperature ramp of 50 K/ns in the NPT ensemble was applied in order to reach the target value of 298 K, allowing the system to change the density according to the change in temperature imposed. Slower cooling rates were tested as well, and no appreciable difference in results was observed. Each system consisted of 5 chains of 150 monomers each, leading to molecular weight in the range between ~13,000 and 15,000 g/mol, number of atoms between ~9000 and 11,000, and box size between ~45 and 50 Å. For every system, three independent initial configurations were generated and simulated, allowing the extraction of the average value and standard deviation for each property.
Simulations were performed in full atomistic detail during all equilibration and production runs using the LAMMPS package (3 March 2020 stable release) [42,43]. The polymer consistent force field (PCFF) [44] was used for polymers, methane, and polymer-CH4 systems, while the COMPASS force field [45] was used to describe CO2 and its interactions with the polymer system.
All systems were simulated with periodic boundary conditions, and the cut-off for Lennard-Jones potential and Coulombic interactions was set to 12 Å as the maximum value above which the box energy no longer changes. In order to account for the long-range van der Waals interactions, tail corrections were included, while long-range electrostatics were computed with a particle–particle particle–mesh (pppm) method (relative error in forces calculations set to 10−6) [46]. Nosé–Hoover thermostat and barostat with dumping parameters of 100 fs and 500 fs were used for temperature and pressure control, respectively [47]. A timestep of 1 fs was adopted for runs in the NVT and NPT ensembles, while a shorter timestep of 0.5 fs was used for NVE (constant-energy) runs.
After the initial configurations were generated and cooled to 298 K, a short equilibration in the NVT ensemble was performed, and production NPT runs of 10 ns were performed at 298 K and 1 bar in order to extract average values of density and cohesive energy and compare them with data found in the literature. The systems were then heated to 308 K at 50 K/ns in the NPT ensemble, and NPT runs of 10 ns at 308 K and 1 bar were performed.
As the transport-properties analysis is computationally more expensive with respect to other properties retrieved in this work, we investigated such properties in a lower number of copolymers. The last 5 ns of the trajectories on 4 copolymers, namely PHBV0, PHBV8, PHBV60, and PHBV100 systems, were used to perform the Widom test particle insertions [48] in order to calculate the excess chemical potential of CO2 and CH4 at infinite dilution, as described elsewhere [23,49], and from that the solubility. Although several methods exist to simulate the solubility in polymers, for which we address the reader to competent and extended reviews [20,50], the Widom insertion method is deemed optimal for the low-concentration sorption of small molecules, as is the case inspected here.
A total of 106 insertions per gas per pure polymer system was sufficient for the convergence of the value of the solubility coefficient, calculated as follows:
1 S i = ρ R T M i lim x i 0 e x p μ i e x R T
where ρ is the density of the pure polymer system, M i is the molar mass of CO2 or CH4, and μ i e x is the excess chemical potential of CO2 or CH4 at infinite dilution.
The solubility coefficient obtained from the Widom analysis and averaged over three configurations for each system was used to calculate the corresponding number of gas molecules to insert into each system at atmospheric pressure. The calculation was done under the assumption of a quasi-constant solubility coefficient at temperatures higher than the glass transition (Tg) in PHBV, as confirmed by previous experimental work [19]. After the energy minimization, a short NVT run, a 10 ns NPT run, and a second 1 ns NVT run were performed for equilibration. Afterwards, 50 to 100 ns production runs in the NVE ensemble were performed in order to extract the average mean square displacements (MSDs) of gas molecules, to calculate the self-diffusion coefficients through Einstein’s relation, as follows [47,51]:
D = lim t 1 6 t r t r 0 2
where r ( t ) and r ( 0 ) are the positions of the centre of mass of the gas molecule at time t and at the initial time t = 0 , respectively. Self-diffusivities are a good approximation of binary diffusivities in the case of an infinitely dilute system or when the sorption isotherm is linear [52]. The NVE ensemble is considered the most suitable to study the dynamic properties of the system, as no thermostat or barostat is used in such integration and thus no external influence is imposed on the system [53]. A multiple time origin scheme was considered in the calculation, and the Fickian regime of MSDs was reached in all cases. A representative example of diffusion coefficient calculation from MD simulations is shown in Section S3 of the SI. The NVE ensemble was also used to calculate the radial distribution functions, g ( r ) , to evaluate the intermolecular interactions between gas molecules and the four polymeric structures, as described elsewhere [47]. In order to visualize the structural difference between the polymeric matrices, short NVE simulations were performed after production NPT runs at 308K and 1 bar, and short NVT equilibration runs. The static structure factor, S ( q ) , was then calculated from the radial distribution functions of all atom types pairs and the atomic-form factors, f i ( q ) [54], as described elsewhere [23]. The peaks at a specific magnitude of the wave vector, q , can be compared to the peaks at a specific scattering angle in X-ray scattering patterns.
The complete schematic of the workflow used in this work is shown in Figure 2.
Gas solubility and diffusivity values obtained from simulations were compared to the experimental values from this and previous work [19]. It is important to point out that, while in the case of pure polymer properties such as density and solubility parameters amorphous state data were available for validation, the experimental gas transport tests were carried out on real semi-crystalline samples. Full atomistic simulations do not allow observation of any crystallization in typical simulation times and are representative only of the free amorphous domains. The crystalline phase, however, is impenetrable and does not contribute directly neither to gas sorption or to gas diffusion, so some rules of thumb can be used.
In the case of solubility, the experiments measure the total semi-crystalline value, S i s c , whereas the simulations yield the solubility of a hypothetical free-amorphous phase only, S i a m [50]. Assuming, as discussed above, that the crystalline domains have zero solubility, a simple rule can be used to compare experimental and simulated values:
S i s c = S i a m 1 X c
where X c is the degree of crystallinity.
The scaling of the diffusivity coefficient on the crystalline content of the sample is somehow less straightforward because the crystals are impermeable and lengthen the diffusive path, as the gas molecules have to circumvent crystals to go from one amorphous area to the other. Such an effect depends on the number of crystals present but also on their shape and distribution. Presumably, the effect of crystals on diffusion is larger than what would be predicted by an additive rule such as the one represented by Equation (6), as we foresee overestimates as high as one order of magnitude. It is expected, however, that crystallinity effects are independent of the gas type and that diffusivity–selectivity simulations could be more reliable. Therefore we will only put in direct comparison the simulated and experimental values of the diffusivity–selectivity, leaving aside the gases' individual diffusion coefficients.

3. Results and Discussion

3.1. Validation of the Molecular Polymer Model

When a specific force field is used to describe polymer behaviour, it is essential to confirm its capability to correctly reproduce volumetric and energetic parameters before predictions of the sorption and transport properties of different gases can be undertaken [20,21,22,23]. Figure 3 shows the snapshots of the equilibrated structures for two homopolymers as an example.
It must be noticed that the polymer molecular models represent only the amorphous phase, as the current limits on the fully atomistic simulations considered here do not allow for observation of the polymer crystallization. For the pure polymer, however, the data used for model validation refer indeed to amorphous structures. In particular, PHBV copolymers were studied in the amorphous state by Mitomo et al. [28], who measured densities using the gradient column technique on samples obtained by melting at 185 °C for 2 min and immediately quenching in liquid nitrogen to avoid crystallization. As Figure 4a shows, the experimental density values, as well as their dependence on the amount of HV units present in the copolymer, are well-reproduced by the MD model, with deviations < 1% for all systems, except for PHBV24, where the deviation reaches 1.6%. The deviations were measured considering the linear decrease in the amorphous density, with different slopes at compositions lower and higher than 30 mol% of HV units, as described by Mitomo et al. [28]. In particular, one can see that the density decreases with increasing fraction of HV units in the polymer; this is due to the bulkier side group of the HV monomer, which inhibits the tight packing of the polymer chains.
The second quantity used to validate the molecular polymer model is an energetic one, namely Hildebrand’s solubility parameter, δ ; a thermodynamic property that defines the miscibility and compatibility of the polymer with fluids. Such value is equal to the square root of the intermolecular cohesive energy of the system per unit volume, being higher for substances with strong internal attractive interactions such as hydrogen bonds. For low molecular weight liquids, this value can be estimated from the vaporization energy, while for polymers, it should be determined either by trial and error using experimental solubility data of the polymer in various solvents, or with group contribution methods, which are normally affected by errors of ±10%. In molecular simulations, as in the case of the present work, the intermolecular energy E c o h may be evaluated as the difference between the sum of the intramolecular energy of every single polymer chain and the total potential energy of the simulation box [22]. The solubility parameter is then calculated as E c o h / V .
For pure PHB (PHBV0), the average literature values for δ vary between 19.45 [29] and 21.0 [30] MPa0.5, depending on the method used, while for PHBV100, values range between 18.6 and 19.4 MPa0.5 according to different group contribution methods [29]. The latter methods provide values for amorphous polymers, so they do not need correction before being compared to the ones obtained from atomistic simulations.
Figure 4b shows that the MD-simulated values fall within the uncertainty range estimated for the literature values (±10% [29]). The simulations also catch the qualitative trend of the solubility parameter decreasing with an increasing amount of HV units, which is due to a combination of the lower intrinsic polarity of the polymers richer in HV, as confirmed by decreasing cohesive energy (Figure 4c), and lower density.
In Figure S4 in the SI, we also reported values of the fractional free volume (FFV) for a probe of negligible size and the corresponding accessible surface area; both quantities agree with the trend of density and increase monotonically with the HV content. This is possibly due to the longer alkyl side group of the HV monomer, with respect to HB, which disrupts the packing of the polymer chains even in the amorphous state. The solubility parameter decreases with HV fraction due to the higher amount of non-polar groups brought by the HV monomer, that have a lower value of δ .
It is interesting to analyze the simulated models also in terms of their static structure factor S(q). This can provide further validation of the generated structures, as the features at low values of q reflect intermolecular correlations in the bulk of the polymeric phase that can be compared directly to the peaks observed at low scattering angles in wide-angle X-ray diffraction patterns.
We analysed such features in the same polymers chosen for the study of transport properties, thus in PHBB0, PHBV8, PHBV60, and PHBV100. Figure 5 reports the structure factor weighted on the wave vector q . Interestingly, all four structures present quite a pronounced peak in a narrow range q = 1.63 ± 0.05 Å 1 , which in terms of scattering angle 2 θ translates to 23.8 ± 0.8 ° . This agrees well with the diffraction spectra of PHB and PHBV copolymers reported in the literature [31,32,33], where the amorphous peaks are found in the range of scattering angles 2 θ 22 24 ° . It can be evinced from Figure 5 that the structure of two homopolymers is very similar, with a narrower main peak with respect to the two copolymers. Such a result can be attributed to the fact that homopolymers have a more ordered and regular structure as opposed to random copolymers.
Given the satisfactory comparison with PHBV copolymer density, solubility parameter values, and structure features from the literature, the force field and simulation protocol adopted in this work were validated, and the same model was used to simulate the gas solubility and diffusivity.

3.2. Gas Separation Performance: Experimental Data and MD Simulations

In this work, we obtained direct data on CO2 and CH4 sorption and diffusion, through sorption tests on the PHBV8 copolymer, which can also be used to estimate permeability according to the solution diffusion model. Additional data come from a previous work where a PHBV25 copolymer was characterized [19], specifically with CO2 solubility, diffusivity and permeability values from sorption and permeation measurements and CH4 diffusivity and permeability values from permeation tests.
The measured values of CO2 and CH4 solubility coefficients, normalized for the crystallinity content, as well as the selectivities, are reported in Table S1 and later in the text.
The experimental solubility coefficient in PHBV8, appropriately scaled with the crystallinity fraction with Equation (6), is higher for CO2 than for CH4, which is due to the different condensability and the favourable energetic interactions of CO2 with the ester groups of the polymer. A similar behaviour was observed in PHBV25 in a previous work [19]. The resulting CO2/CH4 solubility–selectivity, evaluated on a molar basis, is equal to 5.0, while in the case of PHBV25, it was 13.9 [19].
Conversely, the diffusion coefficient of CO2 measured in PHBV8 is smaller than that of CH4, despite the fact that such a gas molecule has a larger kinetic diameter. Such behaviour results in a rather low value of diffusivity–selectivity for PHBV8 (0.6). The same value measured in PHBV25 was higher (1.9) but still rather limited considering the different molecular sizes of the two gases. We believe that the rather limited diffusivity–selectivity of CO2 in these matrices is due to the strong interactions that such penetrant has with the polyesters, which hinders its mobility, ultimately making the solubility the main driver of separation in PHBV matrices.
The above considerations can be co”robo’ated by computing the radial distribution functions, g ( r ) , between polymeric atoms and the centre of mass of gas molecules, here represented by the carbon atom in CO2 and CH4, indicated as C-CO2 and C-CH4, respectively. Complete results are reported in Section S4 and Figure S3 of the SI. Figure 6a shows that the peak relative to the interaction between C-CO2 and the oxygen atom on the carbonyl group of the copolymer is at a shorter distance with respect to the one relative to C-CH4. This does not happen systematically with other atom types, as it can be seen in the SI; therefore, it cannot be attributed solely to the smaller size of the CO2 molecule with respect to CH4, but rather to a specific interaction between CO2 molecules and ester groups of PHBVs. On the other hand, Figure 6b shows that both penetrant molecules are at the same distance with respect to carbons of the methyl and ethyl groups on the side chains, but the correlation is stronger in the case of CH4, as shown by the peak heights, suggesting a preferential positioning of the CH4 molecules towards the alkyl groups. However, the difference between CO2-related and CH4-related peaks is not as marked as in Figure 6a, which may justify the overall CO2-selective behaviour of such matrices. Another thing to notice is that, unlike the pure polymer molecular model, which shows very regular trends of density and solubility parameters with HV fraction, the gas-polymer models have RDF trends not monotonous with the content of HV. This is possibly due to the fact that FFV and δ have opposite trends with HV fraction, thus resulting in a more complex dependence. Indeed, as we will see it in the following, the behaviour of gas-polymer diffusivity, solubility, and selectivity will be a less regular function of the copolymer composition.
In Table S1 and Figure 7, we reported the simulated solubility coefficients of CO2 and CH4 in PHBV0, PHBV8, PHBV60, and PHBV100, estimated through the Widom particle insertion method. It can be observed that the predicted value of gas solubility in the copolymer containing 8% of HV is 40% lower than the experimental counterpart in the case of CO2 and 67% lower in the case of CH4. However, such a discrepancy is in line with the current accuracy of the MD prediction of gas solubility in polymers. According to the specialist literature, a difference between measured and simulated values of S within a factor of 3 can be regarded as acceptable when approaches like the Widom particle insertion or the grand-canonical Monte Carlo (GCMC) are used [55,56,57,58,59,60].
The simulated solubility data at different compositions are characterised by a certain scattering, which is common for MD prediction of solubility in these structures, as discussed above, although interpolation of the values with a linear trend would indicate that CO2 solubility increases with the fraction of HV units (Figure 7a). The experimental trend observed on PHBV8 and PHBV25, and on PHBV with 8, 12, and 33 mol% HV from the literature [61] indicates that CO2 solubility increases with HV content, although a plateau is observed above 30% of HV [61], and polymer grades with HV fractions higher than 33% were not explored. Previous works attributed the higher solubility to the higher flexibility of copolymers containing more HV [61], which is consistent with our observed increase of FFV in materials richer in this monomer. On the other hand, the predicted CH4 solubility is constant or slightly decreases with an increasing amount of HV, as confirmed by the two experimental values obtained by our group on PHBV8 and PHBV25, but in this case, there are no independent measurements from other sources for comparison. As mentioned above, increasing the HV content has compounding effects on the binary gas-polymer systems, which show a less regular behaviour than the pure polymer alone.
The simulated value of CO2/CH4 solubility–selectivity for the PHBV8 copolymer is in good agreement with the experimental value, confirming that simulated selectivity factors are more reliable than individual transport coefficients (Figure 7c). The higher values of α S are recorded for PHBV8 and PHBV100, which can be explained by the fact that those polymers show the highest interactions between CO2 and polymer carbonyls, according to the RDF spectra of Figure 6a. Overall, CO2 solubility and CO2/CH4 solubility–selectivity increase respectively by a factor of 2 and 3.2 going from PHBV0 to PHBV100 [62,63].
Diffusivity values were calculated from the mean square displacements of CO2 and CH4 molecules extracted from NVE trajectories. As the direct comparison between simulated values of diffusivity in the amorphous phase and measured ones in the semi-crystalline phase is not meaningful, because the first ones are significantly higher, we compared only the diffusivity–selectivity values in Figure 8 [56,58,59,64]. The model over-predicts this value by a factor of 5 in the case of PHBV8. A clear trend with composition is less visible in this case, although, as for all the other performance indicators, the value of diffusion selectivity for PHBV100 is the highest.
Finally, we computed the values of permselectivity, estimated as the product between solubility and diffusivity contributions, and compared them with experimental data for PHBV8 and PHBV25, as shown in Table S1 and Figure 9. Additionally, the result obtained with a linear interpolation of the simulated permselectivity data is reported with a dashed line in Figure 9.
Based on such estimations, PHBV100 is the most promising candidate in terms of separation capabilities for the CO2/CH4 gas pair, as the permselectivity is higher than 50, an order of magnitude higher than the same value in PHBV0.
It must be pointed out that these results can be regarded as an order of magnitude estimates in light of the approximations regarding the semi-crystalline nature of these materials discussed previously. However, such predictions are useful to focus on the most promising candidates for a specific application.

4. Conclusions

In this work, MD simulations were applied to study the effect of composition on the separation properties of PHBV copolymers beyond current experimental knowledge. PHBVs were identified in previous works as promising sustainable alternatives to traditional membrane materials with an interesting sorption-driven CO2-selective nature. In this work, we focused on the CO2/CH4 separation performance of PHBV copolymers with increasing content of HV. The aim is that of providing hints about the optimal biopolymer formulation with target properties for specific gas separations, accelerating the introduction of sustainable materials to the market by reducing the experimental effort required.
The polymer amorphous molecular model and force field were validated using pure-polymer density and solubility parameter values from the literature for the amorphous state, as well as X-ray diffraction data. The density and solubility parameters of the copolymers decrease with the HV fraction, while the FFV increases due to the longer side chain of this monomer, which disrupts the chain packing. Random copolymers have a less ordered structure than pure homopolymers. The polymer-gas systems were simulated to extract the sorption and diffusion behaviour as a function of HV content and compared to experimental sorption and diffusion data on commercial semi-crystalline copolymers containing 8% and 25% of HV. This comparison is complicated by the semi-crystalline nature of the polymers inspected experimentally.
It is known that PHAs have a limited diffusivity–selectivity but a good solubility–selectivity for CO2, due to the interactions of this molecule with the ester groups of the polymer, which were investigated and explained in this work via the simulated RDFs. CO2 molecules show a preferential positioning around carbonyl polymer groups, indicating strong interactions which may enhance sorption but slow down diffusion with respect to CH4. The simulations can represent satisfactorily the CO2/CH4 solubility–selectivity value, which is particularly interesting in the case of these copolymers, whose CO2 capture performance relies mainly on the solubility contribution. The simulation of the diffusion selectivity of the amorphous systems overpredicts the experimental values obtained on semi-crystalline PHBV8. In general, gas-polymer systems show less regular trends with HV content than the pure polymer ones, because such monomer has opposite effects on the polymer FFV and solubility parameter δ .
However, the highest value of permselectivity is that of the pure homopolymer PHBV 100. This trend is consistent with previous experimental works that indicate CO2 solubility increases with HV content.
Although still limited in its accuracy, molecular simulation represents a tool to explore the performance of materials not available experimentally and to provide fundamental insight into molecular-scale behaviour. In part 2 of this series of papers, we will explore the use of a macroscopic modelling tool based on an equation of state and compare its performance with the one of the molecular simulations inspected in this work.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/polym15071805/s1, Section S1. Differential Scanning Calorimetry (DSC); Section S2. Method for the calculation of diffusivity from a sorption test; Section S3. Method for the calculation of diffusivity from Molecular Dynamics simulations; Figure S1. Representative data for MSD calculation from MD simulations. The arrow indicates the data portion used for calculations; Section S4. Radial Distribution Functions; Figure S2. Atom types present in polymeric structures; Figure S3. Radial distribution functions for intermolecolar interactions between different atom types and C-CO2 (a,c,e,g,k,n) and C-CH4 (b,d,f,h,m,o) in PHBB0, PHBV8, PHBV60, and PHBV100; Section S5. Free Fractional Volume and Accessible Surface Area; Figure S4. Free Fractional Volume (a) and Accessible Surface Area (b) for PHBV copolymers as function of HV molar percentage, calculated using penetrant sphere of 0 Å; Table S1. Solubility coefficient, Diffusivity, and selectivities for the CO2/CH4 gas pair in PHBV copolymers [17,19,37,40,47,51,52,65,66].

Author Contributions

Conceptualization, M.G.D.A.; methodology, E.R. and K.P.; software, E.R. and K.P.; validation, K.P.; investigation, K.P.; resources, K.P. and E.R.; data curation, K.P.; writing—original draft preparation, K.P.; writing—review and editing, M.G.D.A. and E.R.; supervision, M.G.D.A. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Acknowledgments

We are thankful to Benedetta Masia for performing gas sorption tests. We gratefully acknowledge the CINECA award under the ISCRA initiative (Call 2019/1) for the availability of high-performance computing resources and support. We gratefully acknowledge Micaela Degli Esposti for the DSC equipment. We gratefully acknowledge Niki Vergadou for the development of an initial version of the Widom test particle-insertion code.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bernardo, P.; Drioli, E.; Golemme, G. Membrane Gas Separation: A Review/State of the Art. Ind. Eng. Chem. Res. 2009, 48, 4638–4663. [Google Scholar] [CrossRef]
  2. Russo, F.; Galiano, F.; Iulianelli, A.; Basile, A.; Figoli, A. Biopolymers for Sustainable Membranes in CO2 Separation: A Review. Fuel Process. Technol. 2021, 213, 106643. [Google Scholar] [CrossRef]
  3. Galizia, M.; Chi, W.S.; Smith, Z.P.; Merkel, T.C.; Baker, R.W.; Freeman, B.D. 50th Anniversary Perspective: Polymers and Mixed Matrix Membranes for Gas and Vapor Separation: A Review and Prospective Opportunities. Macromolecules 2017, 50, 7809–7843. [Google Scholar] [CrossRef]
  4. Sidhikku Kandath Valappil, R.; Ghasem, N.; Al-Marzouqi, M. Current and Future Trends in Polymer Membrane-Based Gas Separation Technology: A Comprehensive Review. J. Ind. Eng. Chem. 2021, 98, 103–129. [Google Scholar] [CrossRef]
  5. Baker, R.W.; Low, B.T. Gas Separation Membrane Materials: A Perspective. Macromolecules 2014, 47, 6999–7013. [Google Scholar] [CrossRef]
  6. Owen, G.; Bandi, M.; Howell, J.A.; Churchouse, S.J. Economic Assessment of Membrane Processes for Water and Waste Water Treatment. J. Membr. Sci. 1995, 102, 77–91. [Google Scholar] [CrossRef]
  7. Galiano, F.; Briceño, K.; Marino, T.; Molino, A.; Christensen, K.V.; Figoli, A. Advances in Biopolymer-Based Membrane Preparation and Applications. J. Membr. Sci. 2018, 564, 562–586. [Google Scholar] [CrossRef]
  8. Dong, X.; Lu, D.; Harris, T.A.L.; Escobar, I.C. Polymers and Solvents Used in Membrane Fabrication: A Review Focusing on Sustainable Membrane Development. Membranes 2021, 11, 309. [Google Scholar] [CrossRef]
  9. Doi, Y.; Kitamura, S.; Abe, H. Microbial Synthesis and Characterization of Poly(3-hydroxybutyrate-co-3-hydroxyhexanoate). Macromolecules 1995, 28, 4822–4828. [Google Scholar] [CrossRef]
  10. Sudesh, K.; Abe, H.; Doi, Y. Synthesis, Structure and Properties of Polyhydroxyalkanoates: Biological Polyesters. Prog. Polym. Sci. 2000, 25, 1503–1555. [Google Scholar] [CrossRef]
  11. Winnacker, M. Polyhydroxyalkanoates: Recent Advances in Their Synthesis and Applications. Eur. J. Lipid Sci. Technol. 2019, 121, 1900101. [Google Scholar] [CrossRef] [Green Version]
  12. Costa, S.S.; Miranda, A.L.; de Morais, M.G.; Costa, J.A.V.; Druzian, J.I. Microalgae as Source of Polyhydroxyalkanoates (PHAs)—A Review. Int. J. Biol. Macromol. 2019, 131, 536–547. [Google Scholar] [CrossRef] [PubMed]
  13. Hori, Y.; Yamaguchi, A.; Hagiwara, T. Chemical Synthesis of High Molecular Weight Poly(3-hydroxybutyrate-co-4-hydroxybutyrate). Polymer 1995, 36, 4703–4705. [Google Scholar] [CrossRef]
  14. Boey, J.Y.; Mohamad, L.; Khok, Y.S.; Tay, G.S.; Baidurah, S. A Review of the Applications and Biodegradation of Polyhydroxyalkanoates and Poly(lactic acid) and Its Composites. Polymers 2021, 13, 1544. [Google Scholar] [CrossRef] [PubMed]
  15. Vert, M. Aliphatic Polyesters: Great Degradable Polymers That Cannot Do Everything. Biomacromolecules 2005, 6, 538–546. [Google Scholar] [CrossRef] [PubMed]
  16. Degli Esposti, M.; Chiellini, F.; Bondioli, F.; Morselli, D.; Fabbri, P. Highly Porous PHB-Based Bioactive Scaffolds for Bone Tissue Engineering by in Situ Synthesis of Hydroxyapatite. Mater. Sci. Eng. C 2019, 100, 286–296. [Google Scholar] [CrossRef] [PubMed]
  17. Follain, N.; Chappey, C.; Dargent, E.; Chivrac, F.; Crétois, R.; Marais, S. Structure and Barrier Properties of Biodegradable Polyhydroxyalkanoate Films. J. Phys. Chem. C 2014, 118, 6165–6177. [Google Scholar] [CrossRef]
  18. Siracusa, V.; Ingrao, C.; Karpova, S.G.; Olkhov, A.A.; Iordanskii, A.L. Gas Transport and Characterization of Poly(3 hydroxybutyrate) Films. Eur. Polym. J. 2017, 91, 149–161. [Google Scholar] [CrossRef]
  19. Papchenko, K.; Degli Esposti, M.; Minelli, M.; Fabbri, P.; Morselli, D.; De Angelis, M.G. New Sustainable Routes for Gas Separation Membranes: The Properties of Poly(hydroxybutyrate-co-hydroxyvalerate) Cast from Green Solvents. J. Membr. Sci. 2022, 660, 120847. [Google Scholar] [CrossRef]
  20. Ricci, E.; Minelli, M.; De Angelis, M.G. Modelling Sorption and Transport of Gases in Polymeric Membranes across Different Scales: A Review. Membranes 2022, 12, 857. [Google Scholar] [CrossRef]
  21. Theodorou, D.N.; Suter, U.W. Detailed Molecular Structure of a Vinyl Polymer Glass. Macromolecules 1985, 18, 1467–1478. [Google Scholar] [CrossRef]
  22. Antoniadis, S.J.; Samara, C.T.; Theodorou, D.N. Molecular Dynamics of Atactic Polypropylene Melts. Macromolecules 1998, 31, 7944–7952. [Google Scholar] [CrossRef]
  23. Ricci, E.; Vergadou, N.; Vogiatzis, G.G.; De Angelis, M.G.; Theodorou, D.N. Molecular Simulations and Mechanistic Analysis of the Effect of CO2 Sorption on Thermodynamics, Structure, and Local Dynamics of Molten Atactic Polystyrene. Macromolecules 2020, 53, 3669–3689. [Google Scholar] [CrossRef]
  24. González Torres, M.; Villarreal-Ramírez, E.; Moyaho-Bernal, M.D.L.A.; Álvarez, M.; González-Valdez, J.; Gutiérrez Uribe, J.A.; Leyva Gómez, G.; Cortez, J.R.C. Insights into the Application of Polyhydroxyalkanoates Derivatives from the Combination of Experimental and Simulation Approaches. J. Mol. Struct. 2019, 1175, 536–541. [Google Scholar] [CrossRef]
  25. Bejagam, K.K.; Iverson, C.N.; Marrone, B.L.; Pilania, G. Molecular Dynamics Simulations for Glass Transition Temperature Predictions of Polyhydroxyalkanoate Biopolymers. Phys. Chem. Chem. Phys. 2020, 22, 17880–17889. [Google Scholar] [CrossRef] [PubMed]
  26. Bejagam, K.K.; Iverson, C.N.; Marrone, B.L.; Pilania, G. Composition and Configuration Dependence of Glass-Transition Temperature in Binary Copolymers and Blends of Polyhydroxyalkanoate Biopolymers. Macromolecules 2021, 54, 5618–5628. [Google Scholar] [CrossRef]
  27. Bejagam, K.K.; Gupta, N.S.; Lee, K.-S.; Iverson, C.N.; Marrone, B.L.; Pilania, G. Predicting the Mechanical Response of Polyhydroxyalkanoate Biopolymers Using Molecular Dynamics Simulations. Polymers 2022, 14, 345. [Google Scholar] [CrossRef]
  28. Mitomo, H.; Morishita, N. Structural Changes of Poly(3-hydroxybutyrate-co-3-hydroxyvaIerate) Fractionated with Acetone-Water Solution. Polymer 1995, 36, 2573–2578. [Google Scholar] [CrossRef]
  29. Terada, M.; Marchessault, R.H. Determination of Solubility Parameters for Poly(3-hydroxyalkanoates). Int. J. Biol. Macromol. 1999, 25, 207–215. [Google Scholar] [CrossRef] [PubMed]
  30. Jacquel, N.; Lo, C.-W.; Wu, H.-S.; Wei, Y.-H.; Wang, S.S. Solubility of Polyhydroxyalkanoates by Experiment and Thermodynamic Correlations. AIChE J. 2007, 53, 2704–2714. [Google Scholar] [CrossRef]
  31. Capitán, M.J.; Rueda, D.R.; Ezquerra, T.A. Inhibition of the Crystallization in Nanofilms of Poly(3-hydroxybutyrate). Macromolecules 2004, 37, 5653–5659. [Google Scholar] [CrossRef]
  32. Hurrell, B.L.; Hurrell, B.L. A Wide-Angle X-Ray Scattering Study of the Ageing of Poly(hydroxybutyrate). J. Mater. Sci. 1998, 33, 1709–1713. [Google Scholar] [CrossRef]
  33. Righetti, M.C.; Aliotta, L.; Mallegni, N.; Gazzano, M.; Passaglia, E.; Cinelli, P.; Lazzeri, A. Constrained Amorphous Interphase and Mechanical Properties of Poly(3-hydroxybutyrate-co-3-hydroxyvalerate). Front. Chem. 2019, 7, 790. [Google Scholar] [CrossRef]
  34. Minelli, M.; Cocchi, G.; Ansaloni, L.; Baschetti, M.G.; De Angelis, M.G.; Doghieri, F. Vapor and Liquid Sorption in Matrimid Polyimide: Experimental Characterization and Modeling. Ind. Eng. Chem. Res. 2013, 52, 8936–8945. [Google Scholar] [CrossRef]
  35. Ferrari, M.C.; Galizia, M.; De Angelis, M.G.; Sarti, G.C. Gas and Vapor Transport in Mixed Matrix Membranes Based on Amorphous Teflon AF1600 and AF2400 and Fumed Silica. Ind. Eng. Chem. Res. 2010, 49, 11920–11935. [Google Scholar] [CrossRef]
  36. Galizia, M.; De Angelis, M.G.; Finkelshtein, E.; Yampolskii, Y.P.; Sarti, G.C. Sorption and Transport of Hydrocarbons and Alcohols in Addition-Type Poly(trimethyl silyl norbornene). I: Experimental Data. J. Membr. Sci. 2011, 385–386, 141–153. [Google Scholar] [CrossRef]
  37. Crank, J. The Mathematics of Diffusion; Oxford University Press: Oxford, UK, 1979; ISBN 978-0-19-853411-2. [Google Scholar]
  38. Minelli, M.; Baschetti, M.G.; Doghieri, F.; Ankerfors, M.; Lindström, T.; Siró, I.; Plackett, D. Investigation of Mass Transport Properties of Microfibrillated Cellulose (MFC) Films. J. Membr. Sci. 2010, 358, 67–75. [Google Scholar] [CrossRef]
  39. Wijmans, J.G.; Baker, R.W. The Solution-Diffusion Model: A Review. J. Membr. Sci. 1995, 107, 1–21. [Google Scholar] [CrossRef]
  40. MAPS Platform. Available online: https://www.scienomics.com (accessed on 15 March 2023).
  41. Ramos, J.; Peristeras, L.D.; Theodorou, D.N. Monte Carlo Simulation of Short Chain Branched Polyolefins in the Molten State. Macromolecules 2007, 40, 9640–9650. [Google Scholar] [CrossRef]
  42. Thompson, A.P.; Aktulga, H.M.; Berger, R.; Bolintineanu, D.S.; Brown, W.M.; Crozier, P.S.; in ’t Veld, P.J.; Kohlmeyer, A.; Moore, S.G.; Nguyen, T.D.; et al. LAMMPS—A Flexible Simulation Tool for Particle-Based Materials Modeling at the Atomic, Meso, and Continuum Scales. Comput. Phys. Commun. 2022, 271, 108171. [Google Scholar] [CrossRef]
  43. LAMMPS Molecular Dynamics Simulator. Available online: https://www.lammps.org (accessed on 15 March 2023).
  44. Sun, H.; Mumby, S.J.; Maple, J.R.; Hagler, A.T. An Ab Initio CFF93 All-Atom Force Field for Polycarbonates. J. Am. Chem. Soc. 1994, 116, 2978–2987. [Google Scholar] [CrossRef]
  45. Yang, J.; Ren, Y.; Tian, A.; Sun, H. COMPASS Force Field for 14 Inorganic Molecules, He, Ne, Ar, Kr, Xe, H2, O2, N2, NO, CO, CO2, NO2, CS2, and SO2, in Liquid Phases. J. Phys. Chem. B 2000, 104, 4951–4957. [Google Scholar] [CrossRef]
  46. Hockney, R.W.; Eastwood, J.W. Computer Simulation Using Particles, 1st ed.; CRC Press: Boca Raton, FL, USA, 1988. [Google Scholar] [CrossRef]
  47. Frenkel, D.; Smit, B. Understanding Molecular Simulation: From Algorithms to Applications; Elsevier: Amsterdam, The Netherlands, 2001. [Google Scholar]
  48. Widom, B. Some Topics in the Theory of Fluids. J. Chem. Phys. 1963, 39, 2808–2812. [Google Scholar] [CrossRef]
  49. Shing, K.S.; Chung, S.T. Computer Simulation Methods for the Calculation of Solubility in Supercritical Extraction Systems. J. Phys. Chem. 1987, 91, 1674–1681. [Google Scholar] [CrossRef]
  50. Atiq, O.; Ricci, E.; Baschetti, M.G.; De Angelis, M.G. Modelling Solubility in Semi-Crystalline Polymers: A Critical Comparative Review. Fluid Phase Equilibria 2022, 556, 113412. [Google Scholar] [CrossRef]
  51. Allen, M.P.; Tildesley, D.J.; Tildesley, D.J. Computer Simulation of Liquids; Oxford Science Publications: Oxford, UK, 1987; Reprinted in Clarendon Press: Oxford, UK, 2017; ISBN 978-0-19-855375-5. [Google Scholar]
  52. Freeman, B.; Yampolskii, Y.; Pinnau, I.; Yampolskii, Y. Materials Science of Membranes for Gas and Vapor Separation; Yampolskii, Y., Pinnau, I., Freeman, B., Eds.; John Wiley & Sons, Ltd.: Chichester, UK, 2006; ISBN 978-0-470-02903-9. [Google Scholar]
  53. Vogiatzis, G.G.; Theodorou, D.N. Local Segmental Dynamics and Stresses in Polystyrene–C60 Mixtures. Macromolecules 2014, 47, 387–404. [Google Scholar] [CrossRef] [Green Version]
  54. Atomic Form Factors. Available online: http://lampx.tugraz.at/~hadley/ss1/crystaldiffraction/atomicformfactors/formfactors.php (accessed on 15 March 2023).
  55. Voyiatzis, E.; Stroeks, A. Atomistic Modeling of Hydrogen and Oxygen Solubility in Semicrystalline PA-6 and HDPE Materials. J. Phys. Chem. B 2022, 126, 6102–6111. [Google Scholar] [CrossRef]
  56. Hofman, D.; Ulbrich, J.; Fritsch, D.; Paul, D. Molecular Modelling Simulation of Gas Transport in Amorphous Polyimide and Poly(amide imide) Membrane Materials. Polymer 1996, 37, 4773–4785. [Google Scholar] [CrossRef]
  57. Memari, P.; Lachet, V.; Rousseau, B. Molecular Simulations of the Solubility of Gases in Polyethylene below Its Melting Temperature. Polymer 2010, 51, 4978–4984. [Google Scholar] [CrossRef]
  58. Sun, D.; Zhou, J. Molecular Simulation of Oxygen Sorption and Diffusion in the Poly(lactic acid). Chin. J. Chem. Eng. 2013, 21, 301–309. [Google Scholar] [CrossRef]
  59. Yang, Y.; Narayanan Nair, A.K.; Sun, S. Sorption and Diffusion of Methane and Carbon Dioxide in Amorphous Poly(alkyl acrylates): A Molecular Simulation Study. J. Phys. Chem. B 2020, 124, 1301–1310. [Google Scholar] [CrossRef] [PubMed]
  60. Goubko, M.; Miloserdov, O.; Yampolskii, Y.; Alentiev, A.; Ryzhikh, V. A Novel Model to Predict Infinite Dilution Solubility Coefficients in Glassy Polymers. J. Polym. Sci. Part B Polym. Phys. 2017, 55, 228–244. [Google Scholar] [CrossRef]
  61. Cravo, C.; Duarte, A.R.C.; Duarte, C.M.M. Solubility of Carbon Dioxide in a Natural Biodegradable Polymer: Determination of Diffusion Coefficients. J. Supercrit. Fluids 2007, 40, 194–199. [Google Scholar] [CrossRef]
  62. Williams, L.L.; Rubin, J.B.; Edwards, H.W. Calculation of Hansen Solubility Parameter Values for a Range of Pressure and Temperature Conditions, Including the Supercritical Fluid Region. Ind. Eng. Chem. Res. 2004, 43, 4967–4972. [Google Scholar] [CrossRef]
  63. Panayiotou, C.; Hatzimanikatis, V. The Solubility Parameters of Carbon Dioxide and Ionic Liquids: Are They an Enigma? Fluid Phase Equilibria 2021, 527, 112828. [Google Scholar] [CrossRef]
  64. Tocci, E.; Bellacchio, E.; Russo, N.; Drioli, E. Diffusion of Gases in PEEKs Membranes: Molecular Dynamics Simulations. J. Membr. Sci. 2002, 206, 389–398. [Google Scholar] [CrossRef]
  65. Barham, P.J.; Keller, A.; Otun, E.L.; Holmes, P.A. Crystallization and morphology of a bacterial thermoplastic: Poly-3-hydroxybutyrate. J. Mater. Sci. 1984, 19, 2781–2794. [Google Scholar] [CrossRef]
  66. Kamiya, N.; Sakurai, M.; Inoue, Y.; Chujo, R. Isomorphic behavior of random copolymers: Thermodynamic analysis of cocrystallization of poly(3-hydroxybutyrate-co-3-hydroxyvalerate). Macromolecules 1991, 24, 3888–3892. [Google Scholar] [CrossRef]
Figure 1. Chemical structure of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) PHBV used in this work.
Figure 1. Chemical structure of poly(3-hydroxybutyrate-co-3-hydroxyvalerate) PHBV used in this work.
Polymers 15 01805 g001
Figure 2. Simulation workflow.
Figure 2. Simulation workflow.
Polymers 15 01805 g002
Figure 3. Snapshots of simulation boxes for PHBV0 (a) and PHBV100 (b) in expanded (left) and closed (right) form after equilibration.
Figure 3. Snapshots of simulation boxes for PHBV0 (a) and PHBV100 (b) in expanded (left) and closed (right) form after equilibration.
Polymers 15 01805 g003
Figure 4. Density (a), Hildebrand solubility parameters (b), and Cohesive Energy (c) for PHBV copolymers as a function of HV molar percentage, as obtained from simulations in this work and experiments from the literature. Experimental density for amorphous PHBV from ref. [28]. Solubility parameter for PHBV0 and PHBV100 from [29,30].
Figure 4. Density (a), Hildebrand solubility parameters (b), and Cohesive Energy (c) for PHBV copolymers as a function of HV molar percentage, as obtained from simulations in this work and experiments from the literature. Experimental density for amorphous PHBV from ref. [28]. Solubility parameter for PHBV0 and PHBV100 from [29,30].
Polymers 15 01805 g004
Figure 5. q—weighted static structure factor for several PHBV copolymers.
Figure 5. q—weighted static structure factor for several PHBV copolymers.
Polymers 15 01805 g005
Figure 6. Radial distribution functions for intermolecular interactions between C−CO2 (straight line) or C−CH4 (dotted line) and Carbonyl O (a) or side chain C (b) in PHBB0, PHBV8, PHBV60, and PHBV100.
Figure 6. Radial distribution functions for intermolecular interactions between C−CO2 (straight line) or C−CH4 (dotted line) and Carbonyl O (a) or side chain C (b) in PHBB0, PHBV8, PHBV60, and PHBV100.
Polymers 15 01805 g006
Figure 7. Solubility coefficient values for CO2 (a) and CH4 (b) and CO2/CH4 solubility−selectivity (c) in PHBV copolymers as a function of the molar concentration of HV units. Green squares: MD−simulated values. Yellow dots: experimental solubility data. For CO2 in PHBV25, data from [19]. For CH4 in PHBV25, data obtained as S = P / D in [19]. Data are expressed per kilograms of amorphous polymeric phase (Equation (6)). Dashed lines are linear interpolations of the simulated values.
Figure 7. Solubility coefficient values for CO2 (a) and CH4 (b) and CO2/CH4 solubility−selectivity (c) in PHBV copolymers as a function of the molar concentration of HV units. Green squares: MD−simulated values. Yellow dots: experimental solubility data. For CO2 in PHBV25, data from [19]. For CH4 in PHBV25, data obtained as S = P / D in [19]. Data are expressed per kilograms of amorphous polymeric phase (Equation (6)). Dashed lines are linear interpolations of the simulated values.
Polymers 15 01805 g007
Figure 8. CO2/CH4 Diffusivity−selectivity values in PHBV copolymers as a function of the molar concentration of HV units. Blue squares: MD−simulated values on amorphous PHBV boxes. Yellow dots: experimental data obtained on semi−crystalline PHBV samples (diffusivity values are the timelag of permeation for the PHBV25 sample in [19]).
Figure 8. CO2/CH4 Diffusivity−selectivity values in PHBV copolymers as a function of the molar concentration of HV units. Blue squares: MD−simulated values on amorphous PHBV boxes. Yellow dots: experimental data obtained on semi−crystalline PHBV samples (diffusivity values are the timelag of permeation for the PHBV25 sample in [19]).
Polymers 15 01805 g008
Figure 9. Permselectivity for the CO2/CH4 couple in PHBV copolymers as a function of the molar concentration of HV units. Violet squares: MD−simulated values on amorphous PHBV boxes. Yellow dots: experimental data obtained on semi−crystalline PHBV samples (experimental permeability values are estimated as P = D S for the PHBV8 sample, directly from permeation tests for the PHBV25 sample [19]). The dashed line is obtained through the linear interpolation of the simulated values reported.
Figure 9. Permselectivity for the CO2/CH4 couple in PHBV copolymers as a function of the molar concentration of HV units. Violet squares: MD−simulated values on amorphous PHBV boxes. Yellow dots: experimental data obtained on semi−crystalline PHBV samples (experimental permeability values are estimated as P = D S for the PHBV8 sample, directly from permeation tests for the PHBV25 sample [19]). The dashed line is obtained through the linear interpolation of the simulated values reported.
Polymers 15 01805 g009
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Papchenko, K.; Ricci, E.; De Angelis, M.G. Modelling across Multiple Scales to Design Biopolymer Membranes for Sustainable Gas Separations: 1—Atomistic Approach. Polymers 2023, 15, 1805. https://doi.org/10.3390/polym15071805

AMA Style

Papchenko K, Ricci E, De Angelis MG. Modelling across Multiple Scales to Design Biopolymer Membranes for Sustainable Gas Separations: 1—Atomistic Approach. Polymers. 2023; 15(7):1805. https://doi.org/10.3390/polym15071805

Chicago/Turabian Style

Papchenko, Kseniya, Eleonora Ricci, and Maria Grazia De Angelis. 2023. "Modelling across Multiple Scales to Design Biopolymer Membranes for Sustainable Gas Separations: 1—Atomistic Approach" Polymers 15, no. 7: 1805. https://doi.org/10.3390/polym15071805

APA Style

Papchenko, K., Ricci, E., & De Angelis, M. G. (2023). Modelling across Multiple Scales to Design Biopolymer Membranes for Sustainable Gas Separations: 1—Atomistic Approach. Polymers, 15(7), 1805. https://doi.org/10.3390/polym15071805

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

Article Metrics

Back to TopTop