Next Article in Journal
High-Entropy Perovskites Pr1−xSrx(Cr,Mn,Fe,Co,Ni)O3−δ (x = 0–0.5): Synthesis and Oxygen Permeation Properties
Next Article in Special Issue
Polystyrene Sulfonate Particles as Building Blocks for Nanofiltration Membranes
Previous Article in Journal
Charged Lipids Influence Phase Separation in Cell-Sized Liposomes Containing Cholesterol or Ergosterol
Previous Article in Special Issue
Impacts of Calcium Addition on Humic Acid Fouling and the Related Mechanism in Ultrafiltration Process for Water Treatment
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Water Molecular Dynamics in the Porous Structures of Ultrafiltration/Nanofiltration Asymmetric Cellulose Acetate–Silica Membranes

by
João Cunha
1,2,
Miguel P. da Silva
1,3,
Maria J. Beira
1,2,*,
Marta C. Corvo
4,
Pedro L. Almeida
4,5,
Pedro J. Sebastião
1,2,
João L. Figueirinhas
1,2 and
Maria Norberta de Pinho
1,3
1
Center of Physics and Engineering of Advanced Materials (CeFEMA), Laboratory for Physics of Materials and Emerging Technologies (LaPMET), Instituto Superior Técnico (IST), Universidade de Lisboa (ULisboa), Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal
2
Department of Physics (DF), Instituto Superior Técnico (IST), Universidade de Lisboa (ULisboa), Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal
3
Department of Chemical Engineering (DEQ), Instituto Superior Técnico (IST), Universidade de Lisboa (ULisboa), Av. Rovisco Pais 1, 1049-001 Lisboa, Portugal
4
Centro de Investigação em Materiais (CENIMAT), Faculdade de Ciências e Tecnologia, Universidade Nova de Lisboa, Campus da Caparica, 2829-516 Caparica, Portugal
5
Department of Physics, ISEL, R. Conselheiro Emídio Navarro 1, 1959-007 Lisboa, Portugal
*
Author to whom correspondence should be addressed.
Membranes 2022, 12(11), 1122; https://doi.org/10.3390/membranes12111122
Submission received: 10 October 2022 / Revised: 26 October 2022 / Accepted: 1 November 2022 / Published: 9 November 2022
(This article belongs to the Special Issue Advanced Membrane Technologies for Wastewater Treatment and Recycling)

Abstract

:
This study presents the characterization of water dynamics in cellulose acetate–silica asymmetric membranes with very different pore structures that are associated with a wide range of selective transport properties of ultrafiltration (UF) and nanofiltration (NF). By combining 1 H NMR spectroscopy, diffusometry and relaxometry and considering that the spin–lattice relaxation rate of the studied systems is mainly determined by translational diffusion, individual rotations and rotations mediated by translational displacements, it was possible to assess the influence of the porous matrix’s confinement on the degree of water ordering and dynamics and to correlate this with UF/NF permeation characteristics. In fact, the less permeable membranes, CA/SiO 2 -22, characterized by smaller pores induce significant orientational order to the water molecules close to/interacting with the membrane matrix’s interface. Conversely, the model fitting analysis of the relaxometry results obtained for the more permeable sets of membranes, CA/SiO 2 -30 and CA/SiO 2 -34, did not evidence surface-induced orientational order, which might be explained by the reduced surface-to-volume ratio of the pores and consequent loss of sensitivity to the signal of surface-bound water. Comparing the findings with those of previous studies, it is clear that the fraction of more confined water molecules in the CA/SiO 2 -22-G20, CA/SiO 2 -30-G20 and CA/SiO 2 -34-G20 membranes of 0.83, 0.24 and 0.35, respectively, is in agreement with the obtained diffusion coefficients as well as with the pore sizes and hydraulic permeabilities of 3.5, 38 and 81 kg h 1 m 2 bar 1 , respectively, reported in the literature. It was also possible to conclude that the post-treatment of the membranes with Triton X-100 surfactants produced no significant structural changes but increased the hydrophobic character of the surface, leading to higher diffusion coefficients, especially for systems associated with average smaller pore dimensions. Altogether, these findings evidence the potential of combining complementary NMR techniques to indirectly study hydrated asymmetric porous media, assess the influence of drying post-treatments on hybrid CA/SiO 2 membrane’ surface characteristics and discriminate between ultra- and nano-filtration membrane systems.

1. Introduction

It is established that the structure and dynamics properties of pore-confined molecules are greatly affected by the morphology of porous media [1,2,3]. In membrane separation, the state of water within a membrane’s three-dimensional porous network plays a role in elucidating the mechanisms of its selective mass transfer task. Concertedly, the separation performance of a membrane can be gauged by the interplay of factors such as the pore size, electrical charge, hydrophilic/hydrophobic characteristics of the membrane polymeric or hybrid matrix and the solutes [4,5]. Therefore, the membranes’ porous structure and the state of water within its porous matrix are crucial to understanding the mechanisms of membrane selective transport.
The determination of the accurate morphological features of porous media still represents a challenge as many properties depend not only on the void size’s distribution but also on their connectivity and liquid–surface interactions [6]. Although there is a vast amount of scientific literature focused on microscopic and spectroscopic characterisation for elucidating the mechanisms of membrane selective transport in the active layer structures of integrally skinned cellulose acetate (CA) or cellulose esters membranes [4,7,8,9,10,11,12], this subject is more complex in the study of hybrid CA and silica, CA/SiO2, asymmetric membranes constituting the system of this work [13,14]. Previous studies by de Pinho et al. [15,16] on the characterisation of the water order and dynamics in asymmetric CA/SiO2 hybrid membranes, covering a wide range of ultrafiltration (UF) and nanofiltration (NF) permeation properties, pointed to an essential indication that Nuclear Magnetic Resonance (NMR) relaxometry observables, which are strongly dependent on water–surface interactions due to confinement, can be reliably correlated with the membranes’ asymmetric porous structures and selective permeation performance.
Nuclear Magnetic Resonance (NMR) relaxometry is a widely used experimental technique that enables the study of a large variety of chemical compounds, such as liquid crystals, polymers, ionic liquids and complex food systems, just to name a few [17,18,19,20]. The 1H NMR longitudinal relaxation rate dispersion (R 1 in the function of the 1 H Larmor frequency) is sensitive to molecular motions occurring at timescales ranging from milli- to picoseconds and from slower collective motions in liquid crystalline phases to fast molecular rotations. 1 H NMR relaxometry is especially sensitive to the existence of some degree of confinement, enabling an indirect study of a confining matrix by introducing a well-known liquid, usually water, into its structure. Relaxation-inducing interactions of the probing liquid with the surrounding surfaces, often referred to as rotations mediated by translational displacements, enables the characterization of a given matrix in terms of the effective mean square displacement of the liquid molecules confined in the porous system as well as the degree of order induced by these interactions [21,22,23,24].
In the present work, the main objective is to probe the water molecular dynamics within the porous structure of asymmetric CA/SiO2 hybrid membranes over a wide range of UF and NF permeation properties by 1H NMR relaxometry as a means to assess the effect of the drying post-treatments on the membranes’ asymmetric structure modification.

2. Experimental Section

2.1. Membrane Preparation and Characterization

A series of flat asymmetric CA/SiO2 hybrid membranes were made in a laboratory by coupling the wet phase inversion [25] with sol–gel techniques [26]. The synthesis methodology is described by de Pinho et al. [13]. Membranes were made from casting solutions containing 16.4 wt.% cellulose acetate (CA) polymer (≈30,000 average molecular weight), supplied by Sigma-Aldrich (Steinheim, Germany), a SiO2 content equal to 5 wt.%, and three different solvent system ratios of formamide (enhancing pore-forming agent) and acetone. The acid catalysed hydrolysis of the SiO2 alkoxide sol–gel precursor was promoted in situ by adding deionised water, tetraethyl orthosilicate (TEOS), supplied by Sigma-Aldrich (Steinheim, Germany), and nitric acid to the polymer casting solution. All chemicals were of reagent grade and 65% nitric acid was of technical grade. Membrane films were cast with the aid of a 250 μm calibrated doctor blade, followed by evaporation for 30 s before coagulation in an ice-cold deionised water bath. Table 1 shows the membranes’ casting solution compositions and film-casting conditions used in the preparation of three membranes with distinct UF porous structures, labelled as CA/SiO2-22, CA/SiO2-30 and CA/SiO2-34. In these membrane labels, the second field is represented by numbers 22, 30 and 34, which correspond to the formamide contents of 21.3%, 29% and 32.9% (wt.%), respectively, in the casting solutions.
Following preparation, the asymmetric CA/SiO2 hybrid membranes were conditioned in surfactant mixtures by a procedure adapted from Vos et al. [27]. This treatment was carried out using aqueous solutions of non-ionic surface-active agents composed of glycerol, supplied by PanReac (Darmstadt, Germany), and/or triton X-100, supplied by VWR (Briare, France). In that regard, membrane films were immersed for 15 min in one of the following solutions: (a) an aqueous solution of glycerol 20 vol.% (G20) or (b) an aqueous solution of triton X-100 4 vol.% and glycerol 20 vol.% (GT). All chemicals used in the treatments were of reagent grade and the conductivity of the deionised water was lower than 10 μS cm 1 . For NMR sample preparation, to access the water behavior within the membranes’ porous matrices, the membrane films were immersed in deionised water for 48 h. Excess surface water was gently removed before enclosing a roll of hydrated membrane film in a sealed 5 mm outer diameter NMR tube. The membranes are identified throughout this work by a three-field code: the first code refers to the membrane hybrid matrix (CA/SiO2), followed by a second field relative to the formamide content (in wt.%) in the casting solutions (of 22, 30 and 34), and the third corresponds to the drying membrane post-treatment of G20 or GT.
The membranes were characterised in terms of pure water hydraulic permeability (L p ) and a molecular weight cut-off (MWCO) referring to the molecular weight of the solute that is 95 % retained by the membrane. Details on the characterisation of the membranes studied are described in da Silva et al. [15].

2.2. Methods

1 H NMR Spectroscopy: The series of spectra obtained from the high resolution 1 H NMR relaxometry experiments performed at 7T was analyzed in order to extract the number of Lorentzian components and their respective longitudinal relaxation rates and signal contribution.
1 H NMR Diffusometry: At controlled temperatures and using a probe head with field gradient coils (Bruker Diff 30, Billerica, MA, USA) and a Bruker 7T superconductor connected to a Bruker Avance III NMR console, it was possible to measure the self-diffusion coefficient, D, of water molecules entrapped in the membrane matrix. The applied Pulsed Gradient Stimulated Echo (PGSE) sequence produces an attenuation of the signal intensity for increasing magnetic field gradient strengths, expressed by Equation (1):
I = I 0 e x p γ 1 H 2 g 2 D δ 3 Δ δ 1 3 ,
where γ 1 H is the proton gyromagnetic ratio, g is the gradient strength, δ is the length of the gradient pulses and Δ is the delay between pulsed gradients. Expression (1) does not take into account that water molecules are confined, which means that the obtained diffusion coefficients can be viewed as having apparent values with an order of magnitude that is well-estimated. More exact estimations of the diffusion coefficients would require the development of robust models that take into account the experimental conditions, namely magnetic field gradient pulse durations, which, as far as the authors know, were not yet achieved. In the case of the studied systems, except for pure water, multi-exponential decays were observed, which lead to the addition of the corresponding number of components to Equation (1).
1 H NMR Relaxometry: The longitudinal relaxation rate, R 1 , was measured across a broad frequency range at controlled temperatures. For 1 H Larmor frequencies ranging between 10 kHz and 9 MHz, the measurements were made using a home-developed Fast Field Cycling (FFC) relaxometer [28]. For the remaining frequencies, the conventional inversion recovery technique was applied using the Bruker Avance II console paired with a variable field iron-core Bruker BE-30 electromagnet (10–100 MHz) or with a Bruker Widebore 7T superconductor magnet for the measurements at 300 MHz.

3. Results and Discussion

3.1. Membrane Characterization

Table 2 shows the hydraulic permeability, L p , and molecular weight cut-off, MWCO, of the asymmetric CA/SiO 2 hybrid membranes.
As it can be observed by looking at the hydraulic permeabilities previously obtained by de Pinho et al. [15] for the membrane systems studied in the present work, the CA/SiO 2 -30 and CA/SiO 2 -34 membranes present marked ultrafiltration characteristics, whereas the CA/SiO 2 -22 membrane has a hydraulic permeability tjhat is one order of magnitude lower, thus standing within the border between nano- and ultrafiltration.

3.2. 1H NMR Spectroscopy

Generally, the results from relaxometry experiments are obtained by integrating over the entire 1 H NMR spectrum and fitting the varying amplitudes, proportional to the magnetization along the fixed external magnetic field, to Equation (2). In Figure 1, the model fitting results following spectral integration are exemplified.
M z = M + ( M 0 M ) e τ R 1
In the case of the present work, the high resolution spectrum, obtained at a 7T external magnetic field, was divided into the minimum number of Lorentzian components for which it was possible to determine the longitudinal relaxation rate and the fraction of the population corresponding to each contribution. The obtained results are presented in Appendix A.1.
For the majority of the studied systems, two components are observed. In these cases, the fraction of more confined water molecules, q, relates with the shortest relaxation time, T 1 = R 1 1 , which is highlighted in red in Figure A1Figure A3. For samples CA/SiO 2 -22 G20 and CA/SiO 2 -34 G20, four and three contributions were, respectively, detected. In the case of sample CA/SiO 2 -22 G20, one of the contributions was immediately disregarded in view of the extremely small T 1 (0.038 s), which would make it undetectable at lower frequencies. The other contribution having the shortest longitudinal relaxation time represented a very small percentage of the signal (3%), and it was, therefore, also not considered. For the CA/SiO 2 -34 G20 system, we simply considered the contribution with the shortest relaxation time to be that of water molecules in a more confined environment. The list of more confined water population fraction is presented in Table 3.
As it can be immediately concluded from the results presented in Table 3, the CA/SiO 2 -22 G20 membrane is dramatically different from CA/SiO 2 -30 G20 and CA/SiO 2 -34 G20 systems in terms of more confined population fraction or, in other words, the surface-to-volume ratio is much larger for CA/SiO 2 -22 G20. This result is consistent with the smaller pores observed for the CA/SiO 2 -22 membranes and the consequent lower hydraulic permeability of this system (see Table 2). The post-treatment with triton X-100 (GT) appears to have uniformized the confined population ratio for the three membrane compositions.

3.3. 1H NMR Diffusometry

Figure 2 shows the model fitting analyses made for each of the studied hydrated membranes, and Table 4 presents the obtained diffusion coefficients. The model fitting to the diffusometry and relaxometry data was performed using the open access online platform at fitteia.org (accessed in 1 September 2022), fitteia®, which applies the non-linear least squares minimization method with a global minimum target provided by the powerful MINUIT numerical routine from the CERN library [29,30].
As it can be observed, all hydrated membranes present at least two diffusion coefficients that can be associated with the water molecules experiencing different degrees of confinement. For the CA/SiO 2 -30 (G20 and GT) and CA/SiO 2 -34 (G20 and GT) systems, three diffusion components were observed. In all these cases, the third residual component can only be observed in the logarithmic scale. In the case of CA/SiO 2 -22 (G20 and GT), the slowest component is probably not observable due to its smaller value, which may fall out of the measurable range for this technique.
From Table 4, it is possible to conclude that the CA/SiO 2 -22 systems present much smaller diffusion coefficients than the CA/SiO 2 -30 and CA/SiO 2 -34 systems, which is expected in view of their smaller pores. Membranes CA/SiO 2 -30 and CA/SiO 2 -34 seem to be harder to distinguish in terms of the diffusion coefficient, possibly because their higher permeability increases the relative amount of less confined water. The slower and intermediate diffusion coefficients, D s l o w and D i n t , respectively, seem to be smaller for CA/SiO 2 -30 systems, which is consistent with the smaller pore sizes observed for these membranes [15]. However, the faster diffusion component is larger for membranes CA/SiO 2 -30 than for membranes CA/SiO 2 -34, which might be a consequence of a the pore size distribution in membranes CA/SiO 2 -30 varying across a broader range of characteristic lengths. The fact that previous SEM studies have shown a wide distribution of pore sizes in these membranes makes it difficult to compare the 1 H NMR diffusometry results obtained for membranes CA/SiO 2 -30 and CA/SiO 2 -34 [15]. Nevertheless, the CA/SiO 2 -22 systems are markedly less permeable and lead to much smaller diffusion coefficients, rendering the comparison between this and the CA/SiO 2 -30 and CA/SiO 2 -34 systems meaningful.

3.4. 1H NMR Relaxometry

3.4.1. Raw Data and Theoretical Models

In Figure 3, the 1 H NMR relaxometry profiles obtained for the membranes studied in the present work are obtained. In order to enable a comparison between the profiles of membranes that were subject to a different drying process, the results previously obtained by de Pinho et al. [16] for membranes dried using the solvent exchange procedure were also added to the figure.
As it can be immediately concluded from the observation of the longitudinal relaxation profiles displayed in Figure 3, for the systems studied in the present work (black circles—G20—and blue squares—GT), the CA/SiO 2 -30 and CA/SiO 2 -34 membranes present rather similar relaxometry profiles, while CA/SiO 2 -22 membranes presented significantly different results.
It is also possible to see that post-treatment with triton X-100 leads to very small differences for CA/SiO 2 -30 and CA/SiO 2 -34 membranes, while it produces a significant R 1 decrease across the lower frequency range in the case of membranes CA/SiO 2 -22.
Furthermore, comparing the results obtained in the present work with those related to the solvent-exchange—(SE) dried membranes (red diamonds)—it is possible to observe a significant difference for the CA/SiO 2 -34 membranes while the CA/SiO 2 -22 systems seem almost insensitive to the drying process. This result may be explained by the fact that membranes that are more permeable, such as CA/SiO 2 -34, are bound to be more impacted by the drying process than systems with smaller pores. In fact, the hydraulic permeabilities found for these systems, reported in previous studies, also show that the permeability of the CA/SiO 2 -22 membrane is almost unaffected by post-treatment drying processes, while permeabilities obtained for the CA/SiO 2 -30 and CA/SiO 2 -34 systems vary over a wider range of values, especially when comparing the SE drying process to the G20 and GT post-treatments [15].
The curves presented in Figure 3 representing the longitudinal relaxation rate, R 1 , obtained at different magnetic fields (or 1 H Larmor frequencies) and called NMR dispersion (NMRD) curves encode information on the molecular dynamics of the systems under analysis. In the present work, it was considered that the water entrapped in the membranes’ pores may relax as a result of rotational and translational diffusions and rotations mediated by translational displacements, which are motivated by the interactions of water molecules with the porous matrix. Furthermore, assuming that these mechanisms are effective at different time scales and, thus, independent of one another, the total relaxation rate may be written as the sum of the individual rates (Equation (3)):
R 1 = R 1 R o t + ( 1 q ) R 1 S D + q R 1 R M T D ,
where q is the fraction of water molecules interacting with the pore walls, which was determined with the analysis of the spectral components of the signals obtained at a 1 H Larmor frequency of 300 MHz:
  • Rotational diffusion (Rot):
    The model by Bloembergen, Purcell and Pound, better known as the BPP model, was applied in order to describe rotations of water molecules in the membranes [31,32]. The contribution of this mechanism to the NMR dispersion curves of water 1 H spins is given by Equation (4).
    R 1 R o t = A R o t τ R o t 1 + ω 2 τ R o t 2 + 4 τ R o t 1 + 4 ω 2 τ R o t 2 .
    The prefactor A R o t depends on the effective intramolecular distance between 1 H nuclear spins, r e f f ( 1.58 Å in the case of the water molecule), via Expression (5), which can easily be calculated for the water molecule:
    A R o t = 3 10 μ 0 4 π 2 γ I 4 2 1 r e f f 6 ,
    with μ 0 denoting the vacuum magnetic permeability ( 4 π × 10 7 H/m), γ I denoting the magnetic ratio of the nucleus with spin I and = h / ( 2 π ) denoting the reduced Planck constant ( 1.0545718 × 10 34 m 2 Kg/s). Given that A r o t can be estimated and fixed, the only parameter in Equation (4) that needs to be determined via the model-fitting analysis is the rotational correlation time, τ R o t .
  • Translational Diffusion (SD):
    Self-diffusion of water molecules may be accounted for using the Torrey model [33,34]. Torrey assumed that molecules have equal probabilities of jumping in any direction from an initial state into another, reaching a random jump-like solution. The associated longitudinal relaxation rate frequency dependence is described by Equation (6).
    R 1 S D = 3 2 μ 0 4 π 2 γ I 4 2 I ( I + 1 ) j ( 1 ) ( ω , τ D , d , r , n ) + j ( 2 ) ( 2 ω , τ D , d , r , n ) .
    Parameter n is the 1 H spin density, and d is the average intermolecular interspin distance. τ D , the translational diffusion correlation time, < r 2 > , the mean square jump distance, and the diffusion coefficient, D, are related by the following equation.
    < r 2 > = 6 τ D D .
    The functions j ( i ) ( ω , τ D , d , r , and n ) are the spectral density functions described in references [33,34].
  • Rotations mediated by translational displacements (RMTD):
    The water motion in the confined system gives rise to a relaxation mechanism associated with rotations mediated by translational displacements. This model describes the movement of water molecules near the pores’ walls and, therefore, is related to the interaction of those molecules with the membranes’ surfaces. The contribution of this model to the longitudinal relaxation rate is given by [35,36] the following:
    R 1 R M T D = A R M T D ν p G ( ν , ν m a x , ν m i n ) = A R M T D f ν m a x ν f ν m i n ν ν p + 4 f ν m a x 2 ν f ν m a x 2 ν 2 ν p ,
    where
    f ( x ) = 1 π a r c t a n ( 2 x + 1 ) + a r c t a n ( 2 x 1 ) a r c t a n 2 x x + 1
    This contribution exhibits one high cut-off frequency, ν m a x , and one low cut-off frequency, ν m a x , which are, respectively, associated with the largest and smallest possible translational relaxation modes and, therefore, to the smallest and largest possible average displacements, respectively: ν m a x 1 = l m i n 2 π / 2 D and ν m i n 1 = l m a x 2 π / 2 D , where D is the diffusion coefficient and l is the average displacement. Exponent p can vary between 0.5 and 1, where p = 0.5 corresponds to a situation where there is an isotropic distribution of coupled rotations and self-diffusion motions along the pore/channel’s surfaces, while for p = 1, there is a preferential orientation of the rotations/translations relaxation modes along the constraining surfaces. The parameter A R M T D is inversely proportional to the square root of the diffusion coefficient and to the range of wave numbers related to the motional modes induced by the surface, Δ q . This parameter is proportional to the square of the fraction of molecules interacting with the surface and to the square of the order parameter, representing the long time limit residual correlation of restricted tumbling.

3.4.2. Model Fitting

In Figure 4 and Figure 5, the model fitting results produced by the Master module of the online platform fitteia® [37] are presented. Figure 4 show the results obtained for pure water and the CA/SiO 2 -22 G20 and GT-hydrated membranes. The model fitting analysis of CA/SiO 2 -30 and CA/SiO 2 -34 systems is presented in Figure 5. The model fitting parameters resulting from the NMRD curves analysis are summarized in Table 5.
The model proposed in this work and the combination of 1 H NMR relaxometry and diffusometry experimental techniques allowed for a consistent analysis of all the studied hydrated membranes, as it can be concluded by the good quality of the fits.
In Figure 4, it is possible to observe the striking NMRD profile difference when comparing pure water with confined water. Water molecules entrapped in the matrix have the additional RMTD relaxation pathway, which significantly increases the longitudinal relaxation rate. Moreover, confined water presented diffusion coefficients that are up to three orders of magnitude smaller than that measured for free water (see inserted image in Figure 4a and Table 5).
The parameter q was fixed to the value obtained from the analysis of the spectral components. D S D and D R M T D were set equal to the value of D f a s t and D i n t obtained from the diffusometry analysis and presented in Table 4, respectively. D S D corresponds to a less confined fraction of water that does not interact directly with the matrix, while D R M T D corresponds to a more confined fraction of water that relaxes as a result of interactions with the surface.
Despite the apparent similarities between the relaxometry profiles obtained for the G20 and GT versions of membranes CA/SiO 2 -30 and CA/SiO 2 -34, the model-fitting analysis evidences a decrease in the self-diffusion relaxation rate contribution for the membranes that were post-treated with triton X-100 (compare the dashed red line of the sub figures with that of the respective inserted image in Figure 5). This contribution decrease is more significant for the CA/SiO 2 -22 membranes, as observed in Figure 4. This observation is consistent with triton X-100 increasing the hydrophobicity of the cellulose acetate matrix, making the water less bound to it and leading to higher diffusion coefficients (see D S D and D R M T D in Table 5 or, respectively, D f a s t and D i n t in Table 4). The fact that this increase is more significant for the CA/SiO 2 -22 hydrated membranes may be explained by the fact that smaller pore sizes relate to a larger ratio of water/surface interactions.Furthermore, the increased hydrophobicity suggested from this model-fitting analysis might explain the uniformization of the bound water fraction, q, found for the GT porous membranes (see Table 3).
This analysis enabled the estimation of the characteristic pore size given by the parameter l m a x , that, on average, induces more effective 1 H NMR relaxation through rotations mediated by translational displacements. As it can be observed, the additional treatment with triton X-100 does not significantly affect this dimension, except in the case of CA/SiO 2 -34 systems. Combining the previously described increased hydrophobicity with the fact that this membrane is the most permeable, it is possible that the signal from more bound water molecules is masked by the signal of unbound water, leading to an apparently larger characteristic dimension.
Regarding the fact that A R M T D is inversely proportional to the square root of the diffusion coefficient, the values obtained for this parameter seem to be consistent for all the samples and further support the increased hydrophobicity conferred upon treating the matrix with triton X-100 (GT).
Parameter p shows that there is an isotropic distribution of coupled rotations and self-diffusion motions along the matrix’ pores for all systems, except for the CA/SiO 2 -22 pair, where some degree of anisotropy is detected. The fact that CA/SiO 2 -22 membranes have smaller pores is expected to increase the degree of confinement, thus evidencing water-ordering induced by the surface.

4. Conclusions

In this study, 1 H NMR spectroscopy, diffusometry and relaxometry were successfully combined in order to consistently analyze three pairs of hydrated ultrafiltration/nanofiltration asymmetric cellulose acetate–silica membranes. Each CA/SiO 2 -22, CA/SiO 2 -30 and CA/SiO 2 -34 pair of membranes was composed of one membrane in which the post-treatment involved an aqueous solution of glycerol with 4 vol.% of triton X-100 (GT) and another where triton was not involved in the post- treatment (G20).
The results seem to be consistent with the post-treatment with triton X-100 rendering the matrix surfaces more hydrophobic and increasing the self-diffusion coefficients obtained for water molecules in different confinement environments. This impact is more significant when the characteristic pore sizes are smaller given the increased probability of water/matrix interactions.
Comparing the results obtained in the present work with those related to membranes dried using the solvent-exchange procedure, presented in previous studies, it becomes clear that the drying process has a much less pronounced impact on the cellulose acetate–silica matrix when the pores are characterized by smaller dimensions.
The surface-bound water population variation observed between the CA/SiO 2 -22 G20 and the two analogous membrane systems is in line not only with the diffusion coefficients obtained in the present work but also with the hydraulic permeabilities reported in previous studies.
On the whole, this work evidences the advantage of combining complementary experimental techniques with a relatively simple relaxation model to study and differentiate between ultrafiltration/nanofiltration porous media and track their sensitivity to different post-treatments/drying processes.

Author Contributions

Conceptualization, all authors; software, P.J.S. and J.L.F.; formal analysis, J.C., J.L.F., M.J.B. and P.J.S.; investigation, J.C., M.P.d.S., J.L.F., M.J.B., M.C.C. and P.L.A.; resources, M.C.C., P.L.A., M.P.d.S. and M.N.d.P.; writing—original draft preparation, J.C., J.L.F., M.N.d.P., M.P.d.S., M.J.B. and P.J.S.; writing—review and editing, all authors; supervision, P.J.S., J.L.F. and M.N.d.P. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Portuguese Fundação para a Ciência e a Tecnologia (FCT) grant number PTDC/CTM-BIO/6178/2014, UID/CTM/04540/2019, UIDB/04540/2020 and UIDB/50025/2020-2023; and MCC researcher’s contract (2021.03255.CEECIND). M. J. Beira was funded by FCT grant number PD/BD/142858/2018.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The NMR spectrometers at IST and FCT NOVA are part of the National NMR network (PTNMR), supported by FCT (ROTEIRO/0031/2013-PINFRA/22161/2016). The authors would like to thank European COST Action EURELAX CA15029 (2016–2020) and acknowledge the article processing charge’s full waiver granted by the Membranes journal.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

Abbreviations

The following abbreviations are used in this manuscript:
NMRNuclear magnetic resonance;
1HProton, hydrogen-1;
PGSEPulse Gradient Stimulated Echo;
RotRotations;
SDSelf-diffusion;
RMTDRotations Mediated by Translational Displacements;
CACellulose acetate;
NFNanofiltration;
UFUltrafiltration;
CA/SiO2Cellulose acetate/silica;
G20Surfactant conditioning with an aq. sol. of glycerol 20 vol.%;
GTSurfactant conditioning with an aq. sol. of glycerol 20 vol.% and triton x-100 4 vol.%.

Appendix A

Appendix A.1. 1 H NMR Spectra Obtained from the Relaxometry Experiments Performed at 7T

Figure A1. 1H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -22 G20 (a) and GT (b).
Figure A1. 1H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -22 G20 (a) and GT (b).
Membranes 12 01122 g0a1
Figure A2. 1 H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -30 G20 (a) and GT (b).
Figure A2. 1 H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -30 G20 (a) and GT (b).
Membranes 12 01122 g0a2
Figure A3. 1 H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -34 G20 (a) and GT (b).
Figure A3. 1 H NMR spectroscopy analysis of the inversion recovery experimental spectra for membranes CA/SiO 2 -34 G20 (a) and GT (b).
Membranes 12 01122 g0a3

References

  1. Hansen, E.W.; Fonnum, G.; Weng, E. Pore Morphology of Porous Polymer Particles Probed by NMR Relaxometry and NMR Cryoporometry. J. Phys. Chem. B 2005, 109, 24295–24303. [Google Scholar] [CrossRef] [PubMed]
  2. Chakraborty, S.; Kumar, H.; Dasgupta, C.; Maiti, P.K. Confined Water: Structure, Dynamics, and Thermodynamics. Accounts Chem. Res. 2017, 50, 2139–2146. [Google Scholar] [CrossRef] [PubMed]
  3. Mietner, J.B.; Brieler, F.J.; Lee, Y.J.; Fröba, M. Properties of Water Confined in Periodic Mesoporous Organosilicas: Nanoimprinting the Local Structure. Angew. Chem. Int. Ed. 2017, 56, 12348–12351. [Google Scholar] [CrossRef] [PubMed]
  4. Stamatialis, D.F.; Dias, C.R.; de Pinho, M.N. Structure and Permeation Properties of Cellulose Esters Asymmetric Membranes. Biomacromolecules 2000, 1, 564–570. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, F.; Wei, M.; Zhang, X.; Song, Y.; Zhou, W.; Wang, Y. How Pore Hydrophilicity Influences Water Permeability? Research 2019, 2019, 2581241. [Google Scholar] [CrossRef] [Green Version]
  6. Silletta, E.V.; Velasco, M.I.; Gomez, C.G.; Strumia, M.C.; Stapf, S.; Mattea, C.; Monti, G.A.; Acosta, R.H. Enhanced Surface Interaction of Water Confined in Hierarchical Porous Polymers Induced by Hydrogen Bonding. Langmuir 2016, 32, 7427–7434. [Google Scholar] [CrossRef] [Green Version]
  7. Murphy, D.; de Pinho, M.N. An ATR-FTIR study of water in cellulose acetate membranes prepared by phase inversion. J. Membr. Sci. 1995, 106, 245–257. [Google Scholar] [CrossRef]
  8. Dias, C.R.; Rosa, M.J.; de Pinho, M.N. Structure of water in asymmetric cellulose ester membranes—And ATR-FTIR study. J. Membr. Sci. 1998, 138, 259–267. [Google Scholar] [CrossRef]
  9. Stamatialis, D.F.; Dias, C.R.; Norberta De Pinho, M. Atomic force microscopy of dense and asymmetric cellulose-based membranes. J. Membr. Sci. 1999, 160, 235–242. [Google Scholar] [CrossRef]
  10. Ong, R.; Chung, T.; Helmer, B.; De Wit, J. Characteristics of water and salt transport, free volume and their relationship with the functional groups of novel cellulose esters. Polymer 2013, 54, 4560–4569. [Google Scholar] [CrossRef]
  11. del Gaudio, I.; Hunter-Sellars, E.; Parkin, I.P.; Williams, D.; Da Ros, S.; Curran, K. Water sorption and diffusion in cellulose acetate: The effect of plasticisers. Carbohydr. Polym. 2021, 267, 118185. [Google Scholar] [CrossRef] [PubMed]
  12. Ioan, S.; Necula, A.M.; Stoica, I.; Olaru, N.; Olaru, L.; Ioanid, G.E. Surface Properties of Cellulose Acetate. High Perform. Polym. 2010, 22, 598–608. [Google Scholar] [CrossRef]
  13. Mendes, G.; Faria, M.; Carvalho, A.; Gonçalves, M.C.; de Pinho, M.N. Structure of water in hybrid cellulose acetate-silica ultrafiltration membranes and permeation properties. Carbohydr. Polym. 2018, 189, 342–351. [Google Scholar] [CrossRef] [PubMed]
  14. Faria, M.; Moreira, C.; Eusébio, T.; Brogueira, P.; de Pinho, M.N. Hybrid flat sheet cellulose acetate/silicon dioxide ultrafiltration membranes for uremic blood purification. Cellulose 2020, 27, 3847–3869. [Google Scholar] [CrossRef]
  15. da Silva, M.P.; Beira, M.J.; Nogueira, I.D.; Sebastião, P.J.; Figueirinhas, J.L.; de Pinho, M.N. Tailoring the Selective Permeation Properties of Asymmetric Cellulose Acetate/Silica Hybrid Membranes and Characterisation of Water Dynamics in Hydrated Membranes by Deuterium Nuclear Magnetic Resonance. Membranes 2022, 12, 559. [Google Scholar] [CrossRef]
  16. Beira, M.J.; Silva, M.P.; Condesso, M.; Cosme, P.; Almeida, P.L.; Corvo, M.C.; Sebastião, P.J.; Figueirinhas, J.L.; de Pinho, M.N. Molecular order and dynamics of water in hybrid cellulose acetate–silica asymmetric membranes. Mol. Phys. 2019, 117, 975–982. [Google Scholar] [CrossRef]
  17. de Almeida Martins, J.P.; Chávez, F.V.; Sebastião, P.J. NMR molecular dynamics study of chromonic liquid crystals Edicol Sunset Yellow doped with salts. Magn. Reson. Chem. 2014, 52, 540–545. [Google Scholar] [CrossRef]
  18. Overbeck, V.; Appelhagen, A.; Rößler, R.; Niemann, T.; Ludwig, R. Rotational correlation times, diffusion coefficients and quadrupolar peaks of the protic ionic liquid ethylammonium nitrate by means of 1H fast field cycling NMR relaxometry. J. Mol. Liq. 2021, 322, 114983. [Google Scholar] [CrossRef]
  19. Calucci, L.; Pizzanelli, S.; Mandoli, A.; Birczyński, A.; Lalowicz, Z.T.; De Monte, C.; Ricci, L.; Bronco, S. Unravelling Main- and Side-Chain Motions in Polymers with NMR Spectroscopy and Relaxometry: The Case of Polyvinyl Butyral. Polymers 2021, 13, 2686. [Google Scholar] [CrossRef]
  20. Ates, E.G.; Domenici, V.; Florek-Wojciechowska, M.; Gradišek, A.; Kruk, D.; Maltar-Strmečki, N.; Oztop, M.; Ozvural, E.B.; Rollet, A.L. Field-dependent NMR relaxometry for Food Science: Applications and perspectives. Trends Food Sci. Technol. 2021, 110, 513–524. [Google Scholar] [CrossRef]
  21. Ordikhani Seyedlar, A.; Stapf, S.; Mattea, C. Nuclear magnetic relaxation and diffusion study of the ionic liquids 1-ethyl- and 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl)imide confined in porous glass. Magn. Reson. Chem. 2019, 57, 818–828. [Google Scholar] [CrossRef] [PubMed]
  22. Zavada, T.; Kimmich, R. The anomalous adsorbate dynamics at surfaces in porous media studied by nuclear magnetic resonance methods. The orientational structure factor and Lévy walks. J. Chem. Phys. 1998, 109, 6929–6939. [Google Scholar] [CrossRef]
  23. Carvalho, A.; Sebastiao, P.J.; Fonseca, I.; Matos, J.; Goncalves, M.C. Silica and silica organically modified nanoparticles: Water dynamics in complex systems. Microporous Mesoporous Mater. 2015, 217, 102–108. [Google Scholar] [CrossRef]
  24. Kumar, A.; Cruz, C.; Figueirinhas, J.L.; Sebastião, P.J.; Trindade, A.C.; Fernandes, S.N.; Godinho, M.H.; Fossum, J.O. Water Dynamics in Composite Aqueous Suspensions of Cellulose Nanocrystals and a Clay Mineral Studied through Magnetic Resonance Relaxometry. J. Phys. Chem. B 2021, 125, 12787–12796. [Google Scholar] [CrossRef]
  25. Kunst, B.; Sourirajan, S. An approach to the development of cellulose acetate ultrafiltration membranes. J. Appl. Polym. Sci. 1974, 18, 3423–3434. [Google Scholar] [CrossRef]
  26. Brinker, C.J.; Scherer, G.W. Sol-Gel Science; Academic Press: San Diego, CA, USA, 1990. [Google Scholar] [CrossRef]
  27. Vos, K.D.; Burris, F.O. Drying cellulose acetate reverse osmosis membranes. Ind. Eng. Chem. Prod. Res. Dev. 1969, 8, 84–89. [Google Scholar] [CrossRef]
  28. Sousa, D.; Domingos Marques, G.; Cascais, J.M.; Sebastião, P.J. Desktop fast-field cycling nuclear magnetic resonance relaxometer. Solid State NMR 2010, 38, 36–43. [Google Scholar] [CrossRef]
  29. Sebastião, P.J. The art of model fitting to experimental results. Eur. J. Phys. 2014, 35, 015017. [Google Scholar] [CrossRef]
  30. James, F. MINUIT Function Minimization and Error Analysis. In MINUIT Function Minimization and Error Analysis. Reference Manual Version 94.1; CERN: Geneva, Switzerland, 1994. [Google Scholar]
  31. Bloembergen, N.; Purcell, E.M.; Pound, R.V. Relaxation effect in nuclear magnetic resonance absorption. Phys. Rev. 1948, 73, 679–712. [Google Scholar] [CrossRef] [Green Version]
  32. Abragam, A. The Principles of Nuclear Magnetism; Clarendon Press: Oxford, UK, 1961. [Google Scholar]
  33. Torrey, H.C. Nuclear Spin Relaxation By Translational Diffusion. Phys. Rev. 1953, 92, 962–969. [Google Scholar] [CrossRef]
  34. Harmon, J.F.; Muller, B.H. Nuclear Spin Relaxation by Translational Diffusion in Liquid Ethane. Phys. Rev. 1969, 182, 400–410. [Google Scholar] [CrossRef]
  35. Sebastião, P.J.; Sousa, D.; Ribeiro, A.C.; Vilfan, M.; Lahajnar, G.; Seliger, J.; Žumer, S. Field-cycling NMR relaxometry of a liquid crystal above TNI in mesoscopic confinement. Phys. Rev. E 2005, 72, 061702. [Google Scholar] [CrossRef] [PubMed]
  36. Vilfan, M.; Apih, T.; Sebastião, P.J.; Lahajnar, G.; Žumer, S. Liquid crystal 8CB in random porous glass: NMR relaxometry study of molecular diffusion and director fluctuations. Phys. Rev. E 2007, 76, 051708. [Google Scholar] [CrossRef] [PubMed]
  37. Sebastiao, P.J.O.; Beira, M.J.; Cordeiro, R.M.O.; Kumar, A.; Fernandes, J.F.; Ferraz, A.M.P.; Gonçalves, L.N. The art of fitting ordinary differential equations models to experimental results. Eur. J. Phys. 2022, 43, 035807. [Google Scholar] [CrossRef]
Figure 1. Sample curve showing the magnetization recovery for each inversion recovery delay, τ .
Figure 1. Sample curve showing the magnetization recovery for each inversion recovery delay, τ .
Membranes 12 01122 g001
Figure 2. Diffusometry model fitting results obtained for the G20 and GT versions of membranes CA/SiO 2 -22—(a,b); CA/SiO 2 -30—(c,d); and CA/SiO 2 -34—(e,f)—membranes at 22 °C. The dashed-red line represents the fast diffusion contribution, the dot-dashed-blue line represents the intermediate diffusion contribution and the dot-dot-dashed-violet line represents the slow diffusion contribution.
Figure 2. Diffusometry model fitting results obtained for the G20 and GT versions of membranes CA/SiO 2 -22—(a,b); CA/SiO 2 -30—(c,d); and CA/SiO 2 -34—(e,f)—membranes at 22 °C. The dashed-red line represents the fast diffusion contribution, the dot-dashed-blue line represents the intermediate diffusion contribution and the dot-dot-dashed-violet line represents the slow diffusion contribution.
Membranes 12 01122 g002
Figure 3. NMRD profiles obtained for the CA/SiO 2 -22—(a); CA/SiO 2 -30—(b); and CA/SiO 2 -34—(c)—membranes at 22 °C. (*) Data extracted from previous works by de Pinho et al. [16] related to membranes dried using the solvent exchange procedure.
Figure 3. NMRD profiles obtained for the CA/SiO 2 -22—(a); CA/SiO 2 -30—(b); and CA/SiO 2 -34—(c)—membranes at 22 °C. (*) Data extracted from previous works by de Pinho et al. [16] related to membranes dried using the solvent exchange procedure.
Membranes 12 01122 g003
Figure 4. NMRD profiles and model fitting results obtained for the pure water—(a); CA/SiO 2 -22 G20—(b); and CA/SiO 2 -22 GT—(c)—membranes. The dot-dashed-brown lines represents the RMTD contribution, the dashed-red-line represents the self-diffusion contributions and the dot-dot-dashed-blue line represents the rotations/reorientations contribution to the longitudinal relaxation rate profiles.
Figure 4. NMRD profiles and model fitting results obtained for the pure water—(a); CA/SiO 2 -22 G20—(b); and CA/SiO 2 -22 GT—(c)—membranes. The dot-dashed-brown lines represents the RMTD contribution, the dashed-red-line represents the self-diffusion contributions and the dot-dot-dashed-blue line represents the rotations/reorientations contribution to the longitudinal relaxation rate profiles.
Membranes 12 01122 g004
Figure 5. NMRD profiles and model fitting results obtained for systems CA/SiO 2 -30 G20 and GT—(a); and CA/SiO 2 -34 G20 and GT—(b). The dot-dashed-brown lines represents the RMTD contribution, the dashed-red-line represents the self-diffusion contributions and the dot-dot-dashed-blue line represents the rotations/reorientations contribution to the longitudinal relaxation rate profiles.
Figure 5. NMRD profiles and model fitting results obtained for systems CA/SiO 2 -30 G20 and GT—(a); and CA/SiO 2 -34 G20 and GT—(b). The dot-dashed-brown lines represents the RMTD contribution, the dashed-red-line represents the self-diffusion contributions and the dot-dot-dashed-blue line represents the rotations/reorientations contribution to the longitudinal relaxation rate profiles.
Membranes 12 01122 g005
Table 1. Asymmetric CA/SiO2 hybrid membranes film casting solutions and casting conditions.
Table 1. Asymmetric CA/SiO2 hybrid membranes film casting solutions and casting conditions.
Casting Solution Composition (wt.%)
Membrane CA/SiO2-22 CA/SiO2-30 CA/SiO2-34
CA16.416.416.4
Formamide21.329.032.9
Acetone58.851.147.2
TEOS (SiO2 precursor)3.03.03.0
H2O0.50.50.5
HNO34 drops (pH ≈ 2)4 drops (pH ≈ 2)4 drops (pH ≈ 2)
Casting Conditions
Temperature of casting solution (°C)20–25
Temperature of casting atmosphere (°C)20–25
Relative humidity of casting atmosphere (%)40–50
Solvent evaporation time (min)0.5
Gelation mediumIce-cold deionised water (2 h)
Table 2. Characteristics of the asymmetric CA/SiO 2 hybrid membranes [15].
Table 2. Characteristics of the asymmetric CA/SiO 2 hybrid membranes [15].
MembraneHydraulic Permeability,Molecular Weight Cut-Off,
L p (kg h 1 m 2 bar 1 )MWCO (kDa)
CA/SiO 2 -22G203.5 ± 0.24
GT2.2 ± 0.23
CA/SiO 2 -30G2038 ± 214
GT40 ± 329
CA/SiO 2 -34G2081 ± 435
GT62 ± 421
Table 3. Considered fraction of more confined water molecules to apply in the 1 H NMR relaxometry analysis.
Table 3. Considered fraction of more confined water molecules to apply in the 1 H NMR relaxometry analysis.
CA/SiO 2 -22CA/SiO 2 -30CA/SiO 2 -34
G20GTG20GTG20GT
More confined population ratio, q0.830.400.240.400.350.50
Table 4. Diffusion coefficients obtained from the PGSE 1 H NMR experiments performed at 25 °C. The model fitting was performed considering by an uncertainty equal to 5% of the signal intensity for each point.
Table 4. Diffusion coefficients obtained from the PGSE 1 H NMR experiments performed at 25 °C. The model fitting was performed considering by an uncertainty equal to 5% of the signal intensity for each point.
Membrane D f a s t ( 10 10   m 2 /s) D i n t ( 10 11   m 2 /s) D s l o w ( 10 11   m 2 /s)
CA/SiO 2 -22G200.090.11
GT1.400.97
CA/SiO 2 -30G206.509.500.38
GT9.509.601.10
CA/SiO 2 -34G205.5010.02.20
GT6.9016.03.10
Table 5. Parameters obtained from the NMRD model fitting analysis made on the studied membranes. The model fitting was performed considering an uncertainty of the relaxation rate equal to 10% of its value. The 1 H spin density, n, needed for the Torrey model for translational self-diffusion was fixed to the calculated value of 6.69 × 10 28 1 H nuclear spins per cubic meter. Parameters D S D and D R M T D were fixed to the D f a s t and D i n t diffusion coefficients presented in Table 4, respectively. In the case of pure water, the diffusion coefficient was fixed to that presented in Figure 4. The fraction q, representing the more confined water, was also not a free parameters, and its value was set to that presented in Table 3 for each studied hydrated membrane. Parameter A r o t was also calculated and fixed as explained in the rotations model section.
Table 5. Parameters obtained from the NMRD model fitting analysis made on the studied membranes. The model fitting was performed considering an uncertainty of the relaxation rate equal to 10% of its value. The 1 H spin density, n, needed for the Torrey model for translational self-diffusion was fixed to the calculated value of 6.69 × 10 28 1 H nuclear spins per cubic meter. Parameters D S D and D R M T D were fixed to the D f a s t and D i n t diffusion coefficients presented in Table 4, respectively. In the case of pure water, the diffusion coefficient was fixed to that presented in Figure 4. The fraction q, representing the more confined water, was also not a free parameters, and its value was set to that presented in Table 3 for each studied hydrated membrane. Parameter A r o t was also calculated and fixed as explained in the rotations model section.
Parameters CA/SiO 2 -22 CA/SiO 2 -30 CA/SiO 2 -34 H 2 O
G20 GT G20 GT G20 GT
A R o t ( 10 10 s 2 )1.081.081.081.08
τ R o t ( 10 12 s)135665
D S D ( 10 10 m 2 /s)0.091.406.509.505.506.9021
r(Å)3.03.03.03.0
d(Å)2.72.72.72.7
q0.830.400.240.400.350.50
D R M T D ( 10 11 m 2 /s)0.110.979.509.6010.016.0
A R M T D ( 10 3 s ( 1 + p ) )27185.83.73.82.9
l m a x ( 10 9 m)5.35.313121016
p0.560.510.500.50
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Cunha, J.; da Silva, M.P.; Beira, M.J.; Corvo, M.C.; Almeida, P.L.; Sebastião, P.J.; Figueirinhas, J.L.; de Pinho, M.N. Water Molecular Dynamics in the Porous Structures of Ultrafiltration/Nanofiltration Asymmetric Cellulose Acetate–Silica Membranes. Membranes 2022, 12, 1122. https://doi.org/10.3390/membranes12111122

AMA Style

Cunha J, da Silva MP, Beira MJ, Corvo MC, Almeida PL, Sebastião PJ, Figueirinhas JL, de Pinho MN. Water Molecular Dynamics in the Porous Structures of Ultrafiltration/Nanofiltration Asymmetric Cellulose Acetate–Silica Membranes. Membranes. 2022; 12(11):1122. https://doi.org/10.3390/membranes12111122

Chicago/Turabian Style

Cunha, João, Miguel P. da Silva, Maria J. Beira, Marta C. Corvo, Pedro L. Almeida, Pedro J. Sebastião, João L. Figueirinhas, and Maria Norberta de Pinho. 2022. "Water Molecular Dynamics in the Porous Structures of Ultrafiltration/Nanofiltration Asymmetric Cellulose Acetate–Silica Membranes" Membranes 12, no. 11: 1122. https://doi.org/10.3390/membranes12111122

APA Style

Cunha, J., da Silva, M. P., Beira, M. J., Corvo, M. C., Almeida, P. L., Sebastião, P. J., Figueirinhas, J. L., & de Pinho, M. N. (2022). Water Molecular Dynamics in the Porous Structures of Ultrafiltration/Nanofiltration Asymmetric Cellulose Acetate–Silica Membranes. Membranes, 12(11), 1122. https://doi.org/10.3390/membranes12111122

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