Next Article in Journal
Mechanical Properties of PolyJet 3D-Printed Composites Inspired by Space-Filling Peano Curves
Previous Article in Journal
The Preparation and Characterization of Polylactic Acid Composites with Chitin-Based Intumescent Flame Retardants
Previous Article in Special Issue
Removal of Oily Contaminants from Water by Using the Hydrophobic Ag Nanoparticles Incorporated Dopamine Modified Cellulose Foam
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics

by
José Antonio González-Mijangos
1,
Enrique Lima
2,
Roberto Guerra-González
3,
Fernando Iguazú Ramírez-Zavaleta
1 and
José Luis Rivera
1,*
1
Facultad de Ciencias Físico-Matemáticas, Universidad Michoacana de San Nicolás de Hidalgo, Morelia 58000, Mexico
2
Laboratorio de Fisicoquímica y Reactividad de Superficies (LaFReS), Instituto de Investigaciones en Materiales, Universidad Nacional Autónoma de Mexico, Circuito Exterior S/N, CU, Del. Coyoacán, Ciudad de Mexico 04510, Mexico
3
Facultad de Ingeniería Química, Universidad Michoacana de San Nicolás de Hidalgo, Morelia 58000, Mexico
*
Author to whom correspondence should be addressed.
Polymers 2021, 13(20), 3515; https://doi.org/10.3390/polym13203515
Submission received: 15 September 2021 / Revised: 4 October 2021 / Accepted: 7 October 2021 / Published: 13 October 2021
(This article belongs to the Special Issue Advances of Poly- and Nano-Based Membranes in Separation Processes)

Abstract

:
The mechanical stability of nanothin free-standing films made of melted polyethylene chains was predicted via molecular dynamics simulations in the range of 373.15–673.15 K. The predicted critical thickness, t c , increased with the square of the temperature, T, with additional chains needed as T increased. From T = 373.15 K up to the thermal limit of stability for polyethylene, t c values were in the range of nanothin thicknesses (3.42–5.63 nm), which approximately corresponds to 44–55 chains per 100 nm2. The density at the center of the layer and the interfacial properties studied (density profiles, interfacial thickness, and radius of gyration) showed independence from the film thickness at the same T. The polyethylene layer at its t c showed a lower melting T (<373.15 K) than bulk polyethylene.

Graphical Abstract

1. Introduction

Polymeric membranes are a widespread separation technology used to solve many scientific and industrial problems. Commercial membranes used for gas separation are very thin, in the range of 100 to 500 nm, and the permeability to selectively separate gases increases as the thickness of the layer decreases [1,2,3]. Polyethylene membranes are permeated by gases with low critical temperatures independently of the pressure differential, while gases with high critical temperatures permeate the membranes in processes that are highly dependent on the pressure differential, which can indicate an affinity between small molecules and the polyethylene layer [4]. The solution–diffusion model is the most accepted transport mechanism of gases through a membrane [5,6], and it has three barriers: (1) adsorption of permeable gases at the polymer surface (commonly exposed to the permeable gas at a high pressure), (2) diffusion of the adsorbed gases in the polymer layer, and (3) desorption of permeated gases at the opposed interface (commonly exposed to the permeated bulk gas phase at a lower pressure). Therefore, if the thickness of the layer is decreased, the role of the second barrier is minimized.
Previous studies have shown that very thin layers of polymers [7] and atomistic systems [8] are mechanically unstable and break down into droplets below a certain critical thickness, t c . Usually, the polymeric membranes are supported on porous substrates to enhance the mechanical stability of the system [9]. It is theorized that, in unsupported very-thin layers, attractive surface forces (disjoining pressures) increase the amplitude of capillary waves, causing the thin films to shatter, provided they are sufficiently thin. De Vries postulated that if the films are thinner than the t c , the two surfaces will come into contact due to the corrugations of the capillary waves, causing the film to break [10], and this theory have been utilized in several studies of layer stability [11,12,13]. Due to computational limitations, molecular simulations of thin layers are commonly carried out in simulation cells with very small interfacial areas, Ai, in the scale of a few nm2, where the full extent of the capillary waves cannot manifest. Simulations on Ai as large as 600 nm × 600 nm [14] did not develop capillary waves or other thermal fluctuations large enough to break the films; therefore, other phenomena should arise at this peculiar state. Additionally, other phenomena occur as the thickness of the layer is reduced. Two important properties of the polyethylene layer shift to a lower T: the glass transition temperature, Tg [15,16,17,18], and the melting temperature, Tm [19,20], which indicate an effect of the layer thickness on its cohesivity.
In this paper, we investigated nanothin films (<100 nm) of melted polyethylene chains in a vacuum, without impurities or large thermal fluctuations, through molecular dynamics simulations. These simulations were used to predict the limits of mechanical stability over a wide range of T values, below the bulk Tm, and up to the limits of thermal stability. We examined interfacial and structural properties (layer and interface thickness, bulk liquid densities, and radius of gyration, Rg, distributions) for thin layers at their specific t c values and compared their interfacial behaviors with those of wider layers.

2. Methodology

The liquid–vacuum interface of melted polyethylene chains was studied directly through the simulation of a thin layer of liquid in a vacuum using the molecular dynamics method. The system was expected to behave similarly to ionic liquid films with very low vapor pressure and thermal degradation at high T values [21,22]. We simulated layers at T values of 373.15, 473.15, 573.15, and 673.15 K, below the bulk Tm of linear polyethylene (403.65 K for a molecular weight of 100.5 kg/mol) [20], beyond the highest Tg reported (253 K) [23,24], up to T values near the limit of thermal stability (losses of 5% of the weight at 695.55 K, and the maximum loss begins at 747.45 K) [25,26]. The simulation cell consisted of a parallelepiped with A i = 145   Å × 145   Å , which was large enough to contain a stretched chain in each lateral direction, had a variable length in the inhomogeneous direction (depending on the number of chains simulated), and contained between 93 and 1568 C200 polyethylene chains. Chains of 200 carbon units were chosen as a tradeoff between computation time and chain size; larger chains required a larger Ai to find the tc. The chosen Ai is not large enough to avoid size effects; based on the results for atomistic systems at a T close to the triple point, a simulation cell ~20 times larger the lateral length used in this work is needed to minimize the size effects [8,14]; it is estimated that our calculations for tc will have underpredictions of ~30% at temperatures close to the shifted Tm.
The initial systems consisted of ordered, stretched, and oriented chains that were parallel to the interfacial surface, with a chain–chain separation of 5 Å in all directions. The initial position of the sites in the configuration were obtained using an in-house code with randomly distributed velocities. The chains were brought to the liquid–vacuum equilibrium slowly over a period of 1 ns, and the mechanical stability was studied using additional simulations in the NVT ensemble (constant number of molecules, rigid volume of simulation cell, and constant average T) ensemble, with a timestep of 1 fs. The thermostat used was that of Nosé [27], implemented in the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [28].
The chains interacted through the TraPPE potential [29], which is a united atom potential that considers the CH3 and CH2 functional groups as Lennard–Jones sites of interaction. The non-bonded Lennard–Jones sites interacted through their corresponding potential,
U L J = 4 π ε i j [ ( σ i j r i j ) 12 ( σ i j r i j ) 6 ]
where rij represents the reduced separation between Lennard–Jones sites i and j. σij and εij are parameters dependent on the interaction site; for unlike sites, standard arithmetic ( σ i j ) and geometric (εij) combining rules were used. We used a long cutoff radius, r c = 7.5   σ C H 3 C H 3 , to consider all potential interactions and avoid the use of long-range corrections at the end of the simulations, or approximated corrections during the simulation, which masks the true dynamics of the system [8,14,30]. The use of a short r c with long-range corrections can lead to correct predictions of some liquid–vapor and liquid–vacuum equilibrium properties, but the dynamics of the system are dependent on the r c employed [14,21]. The TraPPE potential for linear alkanes uses harmonic potentials for bond and angle bending and a cosine function on the dihedral angles [31]. The TraPPE potential has been employed to study the structural conformations of single polyethylene chains (C1000) [32], semicrystalline polyethylene (C112) [33], and melted polyethylene (C30–C150) [34] and polyethylene blends (C320) [35]. The use of the TraPPE potential to study the crystallization process of polyethylene (C192) has produced results that are in better agreement with experiments compared to when other potentials are used [36].
Liquid interfaces are dynamic, and under vacuum or surrounded by vapor phases, they breathe; the thickness of the layer and the thickness of each interface expands and contracts in a characteristic periodical behavior. Therefore, to properly characterize the distributions of their structural properties we need to demarcate the different regions of the layer: bulk liquid-like and interfacial regions. A commonly used delimiting point of the interfacial regions is the location in the inhomogeneous direction of the Gibbs dividing surface, which is commonly used in binary and multicomponent phase equilibria (including a solvent) to delimit homogeneous fluid phases (vapor or liquid). In this paper, the Gibbs dividing surface delimits the liquid-like regions from the homogeneous vacuum. For each interface (the free-standing layer has 2 liquid–vacuum interfaces), the Gibbs dividing surface is the point in the interface that divides each interfacial region into two sub-regions (one rich and another poor in solvent), and the criteria for the location of the point is that the sub-region poor in solvent will contain an amount of solvent equal to the amount of solvent “lost” (compared to the bulk phase) in the sub-region rich in solvent. For one-component systems, the only component is considered as the solvent of the system.
To obtain more realistic profiles, we allowed the system to move in the inhomogeneous direction and calculated density profiles every 100 steps. In the end, we averaged the profiles by correcting the positions of the profiles according to the position in the inhomogeneous direction of the center of the layer, which was calculated as the midpoint between the two positions of the Gibbs dividing surface,   z G d s . This allowed us to obtain more realistic profiles of the calculated properties. To obtain the positions of z G d s , each of the density profiles (calculated using bins of 0.1 Å) were adjusted to the hyperbolic tangent expression commonly used in liquid–vacuum phase equilibrium studies [21,22]:
ρ ( z ) = 1 2 ρ l   [ 1 t a n h ( z z G d s d ) ]
where ρ l is the average bulk density of the liquid phase. d is a measure of the thickness of the interface and describes the length in the inhomogeneous direction, where the density changes from the bulk liquid-like phase to zero in the vacuum. The “10–90” interfacial thickness, t i = 2.1790 d [37,38], represents the interfacial region in the density profile where the density changes from 10% ( z 10 % ) to 90% ( z 90 % ) of the average bulk ρ l .
Molecular dynamics simulations of thin fluid layers of melted chains of Lennard–Jones sites representing the behavior of linear polyethylene (C200) were performed in an inhomogeneous arrangement to predict the tc at which the thin layers were at the limit of mechanical stability at a given T (Figure 1); narrower layers than those at the tc shrank, and nano-sized holes formed in the layer (Figure 1). The simulation cell was under periodic boundary conditions and represented a thin layer of infinite Ai.
The search for the minimum number of chains that form stable layers was by trial and error. We started with a stable layer with 200 chains, in which we removed a few chains from the interface and equilibrated the system for 10 ns. If the layer was stable, we continued and eliminated a few chains more, and allowed the system to again reach equilibrium (for additional simulations of 10 ns); otherwise, we used the last stable layer. We repeated this procedure, eliminating a progressively lower number of chains (close to the tc, we eliminated only 1 chain at each step) until we encountered the thinnest, most stable liquid layer that would not retract itself to form nano-sized holes. Once we found the number of chains required at tc, we tested the equilibrium and performed additional simulations (100 ns) on the stable thin layer to verify its stability at longer times, because some of the layers became unstable at periods longer than 70 ns; when this happened, we simulated the system with 1 more chain and verified again its stability for 100 ns, until we found the stable system. In a previous publication, we predicted the tc of atomistic Lennard–Jones free-standing thin layers [8], but those systems did not require long periods to become unstable; they only needed a few ns to become either stable or unstable, and remained stable or unstable for longer periods (100 ns). We focused on investigating systems that maintain mechanical stability and have the smallest possible liquid thickness, which corresponds to the tc of mechanical stability.
The equilibration of the system was studied through time profiles of the studied properties, including the total energy of the layer; the system is considered equilibrated because the studied properties reached a plateau, with values oscillating around a constant value. Some properties such as total energy equilibrated faster, but the time of equilibration of the system is established when all properties reached a plateau. To avoid local minima states, additional tests of equilibration were performed through perturbations of Ai; we increased Ai to 5% of its original value in a period of 1 ns and returned to the original Ai value in an additional period of 1 ns. After the perturbation finished, the properties returned to their original values. We tried several initial configurations within our computational restrictions, but most of them failed to produce a stable layer (Figure 2a); the initial configuration that always produced stable layers contained non uniform distributions of chains (Figure 2b). We expected that the tests and perturbations performed, along with the long periods of simulation (100 ns), will produce layers unbiased by the initial configuration.

3. Results

Among the properties available for investigation in these very thin layers using molecular dynamics simulations, we studied how the density changes along the two interfaces generated in the liquid–vacuum system, the possibility of adsorption at the interfaces, the width of the homogeneous zone corresponding to the bulk liquid, any distortion in the structure of the chains due to vicinity of the interfaces, and predicted the limits of mechanical stability.

3.1. Density Profiles

Figure 3a shows the density profiles of two thin liquid layers of melted polyethylene chains (C200) with thicknesses corresponding to their tc values at 373.15 K and 673.15 K. At 373.15 K, the critical layer consisted of 96 chains, while at 673.15 K, the critical layer needed 116 chains to form a mechanically stable free-standing layer at the same Ai. For comparison, we also show a thin layer at 673.15 K with a thickness wider than the corresponding thin layer at its t C . The total density profiles did not exhibit adsorption at the interfaces, probably due to the long length of the chains. Similar density profiles have been reported in the literature for large linear alkanes [39,40,41]. We observed that the system at the lower T did not develop a central zone characterized by a flat density, which corresponds to the homogeneous bulk liquid. Meanwhile, the system at the higher T had a sufficiently wide thickness to develop a thin bulk liquid-like. When compared with wider thin layers, such as the one shown in Figure 3a, at the same T (673.15 K), we found that the average ρ l at the center of the thin layer (−1 Å, 1 Å) at its t c reached the same value as the bulk ρ l of wider layers; this observation was also made for thin layers at the lower T of 373.15 K.
End CH3 sites are freer to move and associate than backbone CH2 sites. We wonder if CH3 sites preferred to accumulate in a specific region of the layers, and how such an accumulation affects the properties of the layers. Density profiles of only the end CH3 sites of the polyethylene chains are shown in Figure 3b. Density profiles for backbone CH2 sites should be very similar to the total profiles, because they represent 99.5% of the total sites. At the higher T of 673.15 K, the CH3 profiles showed the same trend as the density profiles of the whole system, the only difference being that the CH3 profiles reached positions farther away from the bulk liquid (magnified in Figure 3b due to the scale used in the y-axis). At the lower T of 373.15 K, the density profile showed two peaks of adsorption at the interfaces, which likely indicates that, at this T, the system was all interfacial, the system could be considered as two united interfaces, and the adsorption at the interfaces was the result of low mobility of the CH3 sites at this lower T. The peaks showed the system’s maxima at regions outside the bulk liquid-like regions, where some CH2 sites are present, and probably indicate that CH3 sites are phobic to the nature of the bulk liquid-like region. The two peaks of density reached values that represent only 1% of the total bulk liquid-like density, even though they are 0.5% of the total number of sites; therefore, its effect on other properties should also be very limited. The valley between the interfaces likely indicates a barrier for the CH3 sites exchange in this region; crossing from one interfacial peak to the other may require more thermal energy than is available at this T.

3.2. Bulk Liquid Density

By using systems with 288 C200 chains and the same A i , which form wider thin layers than those at their tc (more than double the number of chains), we estimated the average bulk ρ l and its standard deviations, which are shown in Figure 4 as a function of T. Experimental results for the ρ l of linear polyethylene at 0.1 MPa are also shown in Figure 4 [42]. Ideally, we should have compared our results with measurements in a vacuum, but due to the lack of experiments in vacuum or ultravacuum conditions, we compared our results to the available experimental data at the lowest reported pressure (0.1 MPa), which we expected would be close to the measurements at vacuum conditions. The results obtained from the simulations successfully quantitatively reproduced the available experimental results and their trends at T in the range of 473.15–673.15 K.
Between 394 and 414.1 K, the experimental data had a discontinuity, which has been attributed to the transition to molten states; such a transition was not observed in the simulation result at 373.15 K, which followed the trend observed in this work at higher T values. The simulated results were obtained after 100 ns of equilibration of the layers, but we observed that, at a T in the range of 473.15–673.15 K, the layers equilibrated in shorter periods (10 ns), while the system at 373.15 K required the full 100 ns to equilibrate, which indicates that the system at this T was likely near or at the region of metastability. It is also well known that the thickness of polymer layers affects the Tg [15,16,17,18] and Tm [19,20], so the observed trend in this work may indicate that the Tm at nanothin thicknesses was lower than the T reported in the experimental work at bulk liquid conditions (403.65 K for a molecular weight of 100.5 kg/mol) [19,20]. Previously reported simulations using the TraPPE potential have quantitively agreed with the experimental Tm (using C50 chains) [43] and Tg (using C192 chains) [36] bulk values.
The system is not expected to crystallize under simulation, because bulk simulations studying the specific heat of crystalizing melts have shown that crystallization starts at T values below 350 K, with a maximum in the specific heat at 325 K when the cooling rate is of 0.05 K/ns and starts to crystalize below 325 K for a heating rate of 0.5 K/ns [36]. To our knowledge, there are no reported simulations at lower heating rates (probably due to computational restrictions), which will increase the crystallization temperatures.
We also investigated the effect of a smaller ( r c · ( 2.5   σ C H 2 C H 2 ), which is commonly employed in simulations as a tradeoff due to the computational restrictions of many large systems; the system becomes less cohesive due to a large number of unaccounted interactions, and the ρ l reduces (thickness becomes larger) by up to 22% of the value predicted using the full potential.

3.3. Radius of Gyration at the Interfaces

The interfaces of the thin layers at their t c were characterized by histograms (calculated using 200 bins) of the Rg. Guided by the results of the density profiles at these critical conditions, the layers could be regarded as the union of two interfaces, and all chains were considered in the histograms. The histograms for the components, lateral and normal to the interfacial surface, and the total value of Rg for the critical thin layers at 373.15 and 673.15 K are plotted in Figure 5a,b, respectively. The histograms of Rg could be characterized as right-skewed distributions with long tails. The distributions of the lateral components of Rg overlapped each other at 673.15 K, while at 373.15 K, there was a small mismatch close to the top of the distribution, probably indicating that the system was not fully equilibrated even after 100 ns of simulation. A comparison of the lateral distributions of Rg at both T values revealed an overlapping of the distributions. The comparison between the distributions of Rg in the normal direction at both T values showed a more physically intuitive behavior; the mode of the two distributions were the same, but a shorter peak and a longer tail were observed at the larger T (673.15 K), which resulted in a larger average value at 673.15 K. The mode in the normal direction were lower by approximately 1.7 Å than those observed in the tangential direction. The distributions of the total Rg showed small differences; the mode and mean values grew by around 1 Å when the T changed from 373.15 to 673.15 K. Similar changes in the total Rg have been reported in atomistic simulations of bulk linear polyethylene [44]. The distributions obtained in this work can be interpreted as a deformation of the chain coils in the normal direction to the interfacial surface. Such a deformation increased as the T values were lowered. The deformation of the chains in the normal direction as a function of the film thickness has been reported previously from experiments of thin film polymers expanding freely on the surface of water [45], and also observed in the free surface of supported polymers [46]. The lateral and normal snapshots of a chain coil with Rg values corresponding to the modes of their corresponding distributions (373.15 K) showed that chains are contained in well-defined interfaces, and they cover only part of the layer in the normal direction, even though they had space within the layer to expand in that direction; therefore, the observed deformations in the normal direction are probably due to the interaction with the closest interface.
To compare the values of Rg obtained at the t c with those at bulk liquid phases, we computed the Rg profiles for nanothin layers containing 288 C200 chains, which represented layers with more than double the number of chains at the t c at the same T. The profiles for the components of Rg and its total value at 673.15 K and at its t c are shown in Figure 6a,c, respectively. The profiles were averaged over the total period of properties estimation (10 ns). Compared to the density profiles, the Rg profiles were not uniform, because the Rg calculation produced only one point per chain at its center of mass, and the density profiles produced 200 points (corresponding to the CH2 and CH3 sites). Only bins in the profiles with more than 1% of the average number of chains in the bulk ρ l were plotted. The position in the normal direction to the surface interface, at which the Rg value was assigned, corresponded to the center of mass of the chain. The lateral components remained constant along the axis in the normal direction, while the average magnitude of the normal component was constantly reduced in the interfacial region. The reduction in the average normal Rg reached around 4.4 Å at the outermost layer of the interface plotted. Molecular simulations at the free surface of supported polymer films have revealed a similar behavior: deformations of the chain coil in the normal direction at interfaces with very low sensitivity to T and insensitivity to the bulk ρ l [46]. The profile for the total Rg showed the same behavior as the component in the lateral directions, with an average value in the center of the layer of 18.67 ± 0.14 Å.
We compared the Rg profiles to those obtained using a smaller r c of 2.5   σ C H 3 C H 3 , and the components and total value of Rg are shown in Figure 6b,c, respectively. The profiles behaved very similarly to those using a longer rc, and they were wider due to the lower cohesivity and density of the system (Figure 4). The total Rg showed a slightly higher average value in the center of the layer of 18.77 ± 0.15 Å, and the maximum deformation of the chains in the normal direction was approximately only 2.1 Å. The insensitivity of the total Rg to the bulk ρ l of both systems at different rc values indicates that interactions within each chain (intrachain) were probably insensitive to rc (the same rc was used for inter- and intra-chain interactions). Interactions between different chains are sensitive to rc and dictated the density of the layer. The smaller variation of the deformation in the normal direction using a smaller rc was likely the result of unaccounted bulk and interfacial chain–chain interactions.
The effect of the layer thickness on interfacial Rg distributions was evaluated through a comparison of the distributions of Rg at two different thicknesses at the same T. In Figure 7, the distributions of Rg for the total value and its components are shown for thicknesses at its tc (116 chains) and for a thicker layer with 288 chains, both at a T of 673.15 K. For the thicker layer, we accounted for the chains located not only at the interface defined as t i , but also deeper in the bulk liquid, where the normal components of Rg started to deform: (− , −28 Å) for the left interface and (28 Å, ) for the right interface. The distributions of the total value and the components in the lateral directions did not change, but the distribution in the normal direction increased in magnitude close to the region of the mode and decreased the number of chains with higher values. This can be interpreted as a slight increment in the number of molecules that were highly deformed in this region, which was the result of the increment of the thickness of the layer. This effect would likely disappear at a larger thickness once the cohesivity of the liquid layer reached a saturation value.

3.4. Critical Thickness

A measurement of the thickness of thin layers at their limit of mechanical stability was obtained by adjusting the density profiles, averaged every 100 fs using Equation (2). From these adjustments, we obtained 100,000 values (10 ns of properties evaluation) for two parameters that characterized the components of the system (interfaces and “bulk” liquid zone): t G d s and t i . t G d s is the separation in the normal direction between the z G d s values of the interfaces. t G d s and t i increased as T increased and are shown in Figure 8a,b, respectively. The increase in t G d s was not only due to the change of density due to increasing T, but additional chains were also needed to mechanically stabilize the thin layer. The data obtained for both properties showed normal distributions. The standard deviations for t i are plotted in Figure 8b as the errors in the measurements (for t G d s ; the deviations are smaller than the symbols); the standard deviations also increased with T. The sets of data for both properties adjusted well to empiric quadratic expressions on T. At the largest T studied, where the interfacial thickness was at a minimum, the characteristic values were 25.76 Å ( t G d s ) and 8.23 Å ( t i ). If comparisons were made with experimental measurements, it would likely be difficult to determine experimentally where the interfaces ended due to the breathing of the interfaces. Depending on the sensitivity of the experimental technique, different values could be reported. Based on this probable experimental uncertainty, we proposed a criterion for the prediction of t C in a range of separations between the two interfaces, at surfaces at each interface located between 10% ( t G d s + t i ) and 90% of the bulk ρ l ( t G d s t i ). The experimental detection of layers with less than 10% of the total ρ l would likely be very difficult to achieve due to the dynamics of the interfaces (breathing). At 373.15 K, t c would be in the range of 17.67–34.20 Å. Both thicknesses, t G d s t i and t G d s + t i , are represented in Figure 3a by the vertical limits for the thin layer at its t c at 673.15 K. The proposed criterion was based on the importance of long-range dispersion forces (van der Waals) at the interfaces, which create disjoining pressures between the interfaces and were more influential as the thickness of the layers decreased [11,47,48,49,50]. The range of values for t c at 373.15 K represented between 4.5 and 8.7 times the value of the Lennard–Jones parameter σ C H 3 C H 3 . The minimum in the Lennard–Jones potential is at 2 1 / 6 σ   ~   1.12 σ , and the separation between contiguous sites is around this value. Therefore, at 373.15 K, each interface of the bulk layer could be viewed as two monolayers, with the outermost monolayer moving freely into the vacuum, up to distances corresponding to two monolayers ( ~   2.24 σ C H 3 C H 3 ). At 673.15 K, t c was in the range of 6.2–14.3 times σ C H 3 C H 3 and could be viewed as interfaces consisting of three monolayers, with the outermost monolayer moving into the vacuum up to distances corresponding to four monolayers.

4. Conclusions

Molecular dynamics simulations under vacuum conditions predicted t c values for free-standing melted and linear polyethylene C200 chains as small as 1.77–3.42 nm at 373.15 K (depending on where the thickness was measured) and as high as 2.45–5.63 nm at T values near the limit of thermal stability (673.15 K). The required thickness to maintain mechanical stability increased with the square of T. The t c not only increased due to the reduction in the layer density as T increased, but additional chains were also needed to ensure mechanical stability; approximately 44 chains per 100 nm2 at 373.15 K and 55 chains at 673.15 K. At the t c , the system was mostly formed by two united interfaces with interfacial thicknesses between 0.83 nm at 373.15 K and 1.59 nm near the limit of thermal stability.
The layer thickness did not affect interfacial behavior, because the interfaces at the t c developed the same average density profiles and interfacial thicknesses as wider layers. There was not a proper bulk liquid at the t c as there was in wider layers, but at the center of the layer (−1 Å, 1 Å), the density reached values equal to those found in the average bulk liquid of wider layers, which reproduced experimental values previously reported. The liquid density of the layer at 373.15 K followed the trend observed in this work at higher values of T, indicating that the melting point of these nanothin layers was affected by the thickness of the layer and should be below 373.15 K.
The distributions of the Rg showed the independence of the interfacial dynamics of the chains with the layer thickness, and they also demonstrated a deformation in the normal direction of the chain coils, probably due to the high attraction between the two conjoined interfaces. The mode of the right-skewed distributions of the Rg was independent of T; as T increased, the right tail of the distribution also increased, producing larger average values.
The studied systems did not show capillary waves large enough to be the factor that make the layers to break. The studied structural and interfacial properties did not show an enhancement (or reduction) as the thickness of the layer was reduced, not even when they reached its critical thickness; therefore, in future work, we will investigate the role of other properties as the interfacial forces on the mechanical stability at the critical thickness.

Author Contributions

Conceptualization, J.A.G.-M. and J.L.R.; Data curation, J.A.G.-M., E.L. and J.L.R.; Formal analysis, J.A.G.-M., E.L. and J.L.R.; Funding acquisition, R.G.-G. and J.L.R.; Investigation, J.A.G.-M., R.G.-G. and J.L.R.; Methodology, F.I.R.-Z. and J.L.R.; Project administration, J.L.R.; Software, R.G.-G. and J.L.R.; Supervision, F.I.R.-Z. and J.L.R.; Validation, J.A.G.-M. and J.L.R.; Visualization, J.A.G.-M., R.G.-G. and J.L.R.; Writing—original draft, E.L., F.I.R.-Z. and J.L.R.; Writing—review & editing, J.L.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by CONACYT—México (Grant INFR-2016-01-268652) and Universidad Michoacana de San Nicolás de Hidalgo.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Acknowledgments

We thank CONACYT (México) for an infrastructure fellowship and a graduate scholarship. We also thank Universidad Michoacana de San Nicolás de Hidalgo for research funds.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Robeson, L.M. Polymer membranes for gas separation. Curr. Opin. Solid State Mater. Sci. 1999, 4, 549–552. [Google Scholar] [CrossRef]
  2. Vinogradov, N.E.; Kagramanov, G.G. The development of polymer membranes and modules for air separation. J. Phys. Conf. Ser. 2016, 751, 12038. [Google Scholar] [CrossRef]
  3. Selyanchyn, R.; Fujikawa, S. Membrane thinning for efficient CO2 capture. Sci. Technol. Adv. Mater. 2017, 18, 816–827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Stern, S.A.; Mullhaupt, J.T.; Gareis, P.J. The effect of pressure on the permeation of gases and vapors through polyethylene. Usefulness of the corresponding states principle. AIChE J. 1969, 15, 64–73. [Google Scholar] [CrossRef]
  5. Wijmans, J.G.; Baker, R.W. The solution-diffusion model: A review. J. Memb. Sci. 1995, 107, 1–21. [Google Scholar] [CrossRef]
  6. Alqaheem, Y.; Alomair, A.; Vinoba, M.; Pérez, A. Polymeric Gas-Separation Membranes for Petroleum Refining. Int. J. Polym. Sci. 2017, 2017, 4250927. [Google Scholar] [CrossRef]
  7. Angarska, J.K.; Dimitrova, B.S.; Danov, K.D.; Kralchevsky, P.A.; Ananthapadmanabhan, K.P.; Lips, A. Detection of the Hydrophobic Surface Force in Foam Films by Measurements of the Critical Thickness of the Film Rupture. Langmuir 2004, 20, 1799–1806. [Google Scholar] [CrossRef]
  8. Rivera, J.L.; Douglas, J.F. Influence of film thickness on the stability of free-standing Lennard-Jones fluid films. J. Chem. Phys. 2019, 150, 144705. [Google Scholar] [CrossRef]
  9. Agboola, O.; Popoola, P.; Sadiku, R.; Sanni, S.E.; Babatunde, D.E.; Ayoola, A.; Abatan, O.G. Fabrication and Potential Applications of Nanoporous Membranes for Separation Processes BT—Environmental Nanotechnology. In Nanoscience in Food and Agriculture 2; Dasgupta, N., Ranjan, S., Lichtfouse, E., Mishra, B.N., Eds.; Springer: Cham, Switzerland, 2021; pp. 119–171. ISBN 978-3-030-73010-9. [Google Scholar]
  10. De Vries, A.J. Foam stability: Part IV. Kinetics and activation energy of film rupture. Recl. Trav. Chim. Pays-Bas 1958, 77, 383–399. [Google Scholar] [CrossRef]
  11. Vrij, A. Possible mechanism for the spontaneous rupture of thin, free liquid films. Discuss. Faraday Soc. 1966, 42, 23–33. [Google Scholar] [CrossRef]
  12. Ivanov, I.B.; Dimitrov, D.S. Hydrodynamics of thin liquid films. Colloid Polym. Sci. 1974, 252, 982–990. [Google Scholar] [CrossRef]
  13. Ivanov, I.B. Effect of surface mobility on the dynamic behavior of thin liquid films. Pure Appl. Chem. 1980, 52, 1241. [Google Scholar] [CrossRef]
  14. Rivera, J.L.; Douglas, J.F. Reducing uncertainty in simulation estimates of the surface tension through a two-scale finite-size analysis: Thicker is better. RSC Adv. 2019, 9, 35803–35812. [Google Scholar] [CrossRef] [Green Version]
  15. Kim, J.H.; Jang, J.; Zin, W.-C. Thickness Dependence of the Glass Transition Temperature in Thin Polymer Films. Langmuir 2001, 17, 2703–2710. [Google Scholar] [CrossRef]
  16. Fryer, D.S.; Nealey, P.F.; de Pablo, J.J. Thermal Probe Measurements of the Glass Transition Temperature for Ultrathin Polymer Films as a Function of Thickness. Macromolecules 2000, 33, 6439–6447. [Google Scholar] [CrossRef]
  17. Bhattacharya, M.; Sanyal, M.K.; Geue, T.; Pietsch, U. Glass transition in ultrathin polymer films: A thermal expansion study. Phys. Rev. E 2005, 71, 41801. [Google Scholar] [CrossRef] [Green Version]
  18. Inoue, R.; Kanaya, T.; Miyazaki, T.; Nishida, K.; Tsukushi, I.; Shibata, K. Glass transition and thermal expansivity of polystyrene thin films. Mater. Sci. Eng. A 2006, 442, 367–370. [Google Scholar] [CrossRef]
  19. Wang, Y.; Rafailovich, M.; Sokolov, J.; Gersappe, D.; Araki, T.; Zou, Y.; Kilcoyne, A.D.L.; Ade, H.; Marom, G.; Lustiger, A. Substrate Effect on the Melting Temperature of Thin Polyethylene Films. Phys. Rev. Lett. 2006, 96, 28303. [Google Scholar] [CrossRef]
  20. Mohammadi, H.; Vincent, M.; Marand, H. Investigating the equilibrium melting temperature of linear polyethylene using the non-linear Hoffman-Weeks approach. Polymer 2018, 146, 344–360. [Google Scholar] [CrossRef]
  21. Rivera, J.L.; Molina-Rodríguez, L.; Ramos-Estrada, M.; Navarro-Santos, P.; Lima, E. Interfacial properties of the ionic liquid [bmim][triflate] over a wide range of temperatures. RSC Adv. 2018, 8, 10115–10123. [Google Scholar] [CrossRef] [Green Version]
  22. Arroyo-Valdez, J.A.; Viramontes-Gamboa, G.; Guerra-Gonzalez, R.; Ramos-Estrada, M.; Lima, E.; Rivera, J.L. Cation folding and the thermal stability limit of the ionic liquid [BMIM+][BF4−] under total vacuum. RSC Adv. 2021, 11, 12951–12960. [Google Scholar] [CrossRef]
  23. McKechnie, D.; Cree, J.; Wadkin-Snaith, D.; Johnston, K. Glass transition temperature of a polymer thin film: Statistical and fitting uncertainties. Polymer 2020, 195, 122433. [Google Scholar] [CrossRef]
  24. Reiter, G.; Kindl, P. Positron lifetime investigations on linear polyethylene compared to branched polyethylene. Phys. Status Solidi 1990, 118, 161–168. [Google Scholar] [CrossRef]
  25. Sustaita-Rodríguez, J.M.; Medellín-Rodríguez, F.J.; Olvera-Mendez, D.C.; Gimenez, A.J.; Luna-Barcenas, G. Thermal Stability and Early Degradation Mechanisms of High-Density Polyethylene, Polyamide 6 (Nylon 6), and Polyethylene Terephthalate. Polym. Eng. Sci. 2019, 59, 2016–2023. [Google Scholar] [CrossRef]
  26. Zong, R.; Wang, Z.; Liu, N.; Hu, Y.; Liao, G. Thermal degradation kinetics of polyethylene and silane-crosslinked polyethylene. J. Appl. Polym. Sci. 2005, 98, 1172–1179. [Google Scholar] [CrossRef]
  27. Noseé, S. A unified formulation of the constant temperature molecular dynamics methods. J. Chem. Phys. 1984, 81, 511. [Google Scholar] [CrossRef] [Green Version]
  28. Plimpton, S. Fast Parallel Algorithms for Short-Range Molecular-Dynamics. J. Comput. Phys. 1995, 117, 1–19. [Google Scholar] [CrossRef] [Green Version]
  29. Martin, M.G.; Siepmann, J.I. Transferable Potentials for Phase Equilibria. 1. United-Atom Description of n-Alkanes. J. Phys. Chem. B 1998, 102, 2569–2577. [Google Scholar] [CrossRef]
  30. Trokhymchuk, A.; Alejandre, J. Computer simulations of liquid/vapor interface in Lennard-Jones fluids: Some questions and answers. J. Chem. Phys. 1999, 111, 8510–8523. [Google Scholar] [CrossRef]
  31. Jorgensen, W.L.; Madura, J.D.; Swenson, C.J. Optimized intermolecular potential functions for liquid hydrocarbons. J. Am. Chem. Soc. 1984, 106, 6638–6646. [Google Scholar] [CrossRef]
  32. Hagita, K.; Fujiwara, S.; Iwaoka, N. Structure formation of a quenched single polyethylene chain with different force fields in united atom molecular dynamics simulations. AIP Adv. 2018, 8, 115108. [Google Scholar] [CrossRef] [Green Version]
  33. Kumar, V.; Locker, C.R.; in’t Veld, P.J.; Rutledge, G.C. Effect of Short Chain Branching on the Interlamellar Structure of Semicrystalline Polyethylene. Macromolecules 2017, 50, 1206–1214. [Google Scholar] [CrossRef] [Green Version]
  34. Takahashi, K.Z.; Nishimura, R.; Yasuoka, K.; Masubuchi, Y. Molecular Dynamics Simulations for Resolving Scaling Laws of Polyethylene Melts. Polymers 2017, 9, 24. [Google Scholar] [CrossRef] [Green Version]
  35. Moyassari, A.; Gkourmpis, T.; Hedenqvist, M.S.; Gedde, U.W. Molecular dynamics simulation of linear polyethylene blends: Effect of molar mass bimodality on topological characteristics and mechanical behavior. Polymer 2019, 161, 139–150. [Google Scholar] [CrossRef]
  36. Ramos, J.; Vega, J.F.; Martínez-Salazar, J. Molecular Dynamics Simulations for the Description of Experimental Molecular Conformation, Melt Dynamics, and Phase Transitions in Polyethylene. Macromolecules 2015, 48, 5016–5027. [Google Scholar] [CrossRef]
  37. Juárez-Guerra, F.M.; Rivera, J.L.; Zúñiga-Moreno, A.; Galicia-Luna, L.A.; Rico, J.L.; Lara, J. Molecular modeling of thiophene in the vapor-liquid equilibrium. Sep. Sci. Technol. 2006, 41, 261–281. [Google Scholar] [CrossRef]
  38. Rivera, J.L.J.L.; Nicanor-Guzman, H.; Guerra-Gonzalez, R. The Intramolecular Pressure and the Extension of the Critical Point’s Influence Zone on the Order Parameter. Adv. Condens. Matter Phys. 2015, 2015, 258601. [Google Scholar] [CrossRef] [Green Version]
  39. Wen, B.; Sun, C.; Bai, B.; Gatapova, E.Y.; Kabov, O.A. Ionic hydration-induced evolution of decane–water interfacial tension. Phys. Chem. Chem. Phys. 2017, 19, 14606–14614. [Google Scholar] [CrossRef] [Green Version]
  40. Choudhary, N.; Narayanan Nair, A.K.; Che Ruslan, M.F.A.; Sun, S. Bulk and interfacial properties of decane in the presence of carbon dioxide, methane, and their mixture. Sci. Rep. 2019, 9, 19784. [Google Scholar] [CrossRef] [Green Version]
  41. Katiyar, P.; Singh, J.K. The effect of ionisation of silica nanoparticles on their binding to nonionic surfactants in oil–water system: An atomistic molecular dynamic study. Mol. Phys. 2018, 116, 2022–2031. [Google Scholar] [CrossRef]
  42. Sato, Y.; Hashiguchi, H.; Inohara, K.; Takishima, S.; Masuoka, H. PVT properties of polyethylene copolymer melts. Fluid Phase Equilib. 2007, 257, 124–130. [Google Scholar] [CrossRef]
  43. Zhang, W.; Larson, R.G. A metastable nematic precursor accelerates polyethylene oligomer crystallization as determined by atomistic simulations and self-consistent field theory. J. Chem. Phys. 2019, 150, 244903. [Google Scholar] [CrossRef]
  44. Lu, H.; Zhou, Z.; Hao, T.; Ye, X.; Ne, Y. Temperature Dependence of Structural Properties and Chain Configurational Study: A Molecular Dynamics Simulation of Polyethylene Chains. Macromol. Theory Simul. 2015, 24, 335–343. [Google Scholar] [CrossRef]
  45. Kim, J.-H.; Jang, K.-L.; Ahn, K.; Yoon, T.; Lee, T.-I.; Kim, T.-S. Thermal expansion behavior of thin films expanding freely on water surface. Sci. Rep. 2019, 9, 7071. [Google Scholar] [CrossRef] [Green Version]
  46. Hanakata, P.Z.; Douglas, J.F.; Starr, F.W. Local variation of fragility and glass transition temperature of ultra-thin supported polymer films. J. Chem. Phys. 2012, 137, 244901. [Google Scholar] [CrossRef] [Green Version]
  47. Sheludko, A. Thin liquid films. Adv. Colloid Interface Sci. 1967, 1, 391–464. [Google Scholar] [CrossRef]
  48. Ruckenstein, E.; Jain, R.K. Spontaneous rupture of thin liquid films. J. Chem. Soc. Faraday Trans. 2 Mol. Chem. Phys. 1974, 70, 132–147. [Google Scholar] [CrossRef]
  49. Morariu, M.D.; Schäffer, E.; Steiner, U. Capillary instabilities by fluctuation induced forces. Eur. Phys. J. E 2003, 12, 375–381. [Google Scholar] [CrossRef]
  50. Reiter, G.; Sharma, A.; Casoli, A.; David, M.-O.; Khanna, R.; Auroy, P. Thin Film Instability Induced by Long-Range Forces. Langmuir 1999, 15, 2551–2558. [Google Scholar] [CrossRef]
Figure 1. Normal views of the representation of a thin layer (interfaces grow infinitely on the paper plane) of linear polyethylene (C200), made of Lennard–Jones sites in equilibrium at T = 473 K. The simulation cell used A i = 145   Å × 145   Å . Thin layers with 102 or more chains (left) were stable up to 0.1 μs, while layers with 101 or fewer chains (right) contracted and developed stable nano-sized holes.
Figure 1. Normal views of the representation of a thin layer (interfaces grow infinitely on the paper plane) of linear polyethylene (C200), made of Lennard–Jones sites in equilibrium at T = 473 K. The simulation cell used A i = 145   Å × 145   Å . Thin layers with 102 or more chains (left) were stable up to 0.1 μs, while layers with 101 or fewer chains (right) contracted and developed stable nano-sized holes.
Polymers 13 03515 g001
Figure 2. Normal view to the interfacial surface of the representation of uniform (a) and non-uniform (b) initial configurations of a thin layer of linear polyethylene (C200) made of Lennard–Jones sites. The simulation cell used A i = 145   Å × 145   Å with periodic boundary conditions (continuous black line).
Figure 2. Normal view to the interfacial surface of the representation of uniform (a) and non-uniform (b) initial configurations of a thin layer of linear polyethylene (C200) made of Lennard–Jones sites. The simulation cell used A i = 145   Å × 145   Å with periodic boundary conditions (continuous black line).
Polymers 13 03515 g002
Figure 3. (a) Density profiles for the whole system and (b) density profiles for only the end CH3 groups as a function of the inhomogeneous position within the simulation cell at their corresponding t c , for systems of polyethylene chains (C200) at T = 373.15 K (black) and T = 673.15 K (red). Each point represents the average density over 1 ps of simulation. Blue points correspond to a thin layer at 673.15 K with a larger number of chains than the layer at its t c . Vertical dashed and dotted red lines represent z 90 % and z 10 % , respectively, for the layer at 673.15 K and at its t c .
Figure 3. (a) Density profiles for the whole system and (b) density profiles for only the end CH3 groups as a function of the inhomogeneous position within the simulation cell at their corresponding t c , for systems of polyethylene chains (C200) at T = 373.15 K (black) and T = 673.15 K (red). Each point represents the average density over 1 ps of simulation. Blue points correspond to a thin layer at 673.15 K with a larger number of chains than the layer at its t c . Vertical dashed and dotted red lines represent z 90 % and z 10 % , respectively, for the layer at 673.15 K and at its t c .
Polymers 13 03515 g003
Figure 4. Bulk ρ l of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium as a function of T. Black circles represent the results of this work using simulations. The blue circle represents a system simulated at a smaller r c ( 2.5   σ C H 3 C H 3 ). Error bars represent the standard deviation of the data. Red circles represent experimental results of Sato et al. for linear polyethylene at 0.1 MPa [Adapted from Ref. [42] with permission from Elsevier under License Number 5163721342541 (7 October 2021)].
Figure 4. Bulk ρ l of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium as a function of T. Black circles represent the results of this work using simulations. The blue circle represents a system simulated at a smaller r c ( 2.5   σ C H 3 C H 3 ). Error bars represent the standard deviation of the data. Red circles represent experimental results of Sato et al. for linear polyethylene at 0.1 MPa [Adapted from Ref. [42] with permission from Elsevier under License Number 5163721342541 (7 October 2021)].
Polymers 13 03515 g004
Figure 5. Histograms of the Rg of thin liquid layers of melted polyethylene chains (C200) at their t c in a liquid–vacuum equilibrium at (a) 373.15 K and (b) 673.15 K. Black and red points represent the components in the lateral directions to the interfacial surface, blue points represent components normal to the interfacial surface, and green points represent the total values of Rg. Vertical lines represent the mode of the distributions. Snapshots in (a) represent normal and lateral views of a chain with Rg values equal to the mode values of its correspondent distribution (highlighted with red spheres). The rest of the chains are represented with gray lines.
Figure 5. Histograms of the Rg of thin liquid layers of melted polyethylene chains (C200) at their t c in a liquid–vacuum equilibrium at (a) 373.15 K and (b) 673.15 K. Black and red points represent the components in the lateral directions to the interfacial surface, blue points represent components normal to the interfacial surface, and green points represent the total values of Rg. Vertical lines represent the mode of the distributions. Snapshots in (a) represent normal and lateral views of a chain with Rg values equal to the mode values of its correspondent distribution (highlighted with red spheres). The rest of the chains are represented with gray lines.
Polymers 13 03515 g005
Figure 6. Rg profiles as a function of the inhomogeneous position within the simulation cell at thickness layers wider than the corresponding tc for systems of 288 polyethylene chains (C200) at T = 673.15 K, simulated using a r c of (a) 7.5   σ C H 3 C H 3 and (b) 2.5   σ C H 3 C H 3 . Black and red points represent the lateral Rg, while blue points represent the normal Rg. (c) Total Rg profiles for systems simulated using r c = 7.5   σ C H 3 C H 3 (orange) and r c = 2.5   σ C H 3 C H 3 (green). Vertical dashed and dotted lines in (a) and (c) represent z G d s and z 90 % , respectively.
Figure 6. Rg profiles as a function of the inhomogeneous position within the simulation cell at thickness layers wider than the corresponding tc for systems of 288 polyethylene chains (C200) at T = 673.15 K, simulated using a r c of (a) 7.5   σ C H 3 C H 3 and (b) 2.5   σ C H 3 C H 3 . Black and red points represent the lateral Rg, while blue points represent the normal Rg. (c) Total Rg profiles for systems simulated using r c = 7.5   σ C H 3 C H 3 (orange) and r c = 2.5   σ C H 3 C H 3 (green). Vertical dashed and dotted lines in (a) and (c) represent z G d s and z 90 % , respectively.
Polymers 13 03515 g006
Figure 7. Histograms of Rg of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium at their tc at 673.15 K. Black and red points represent the components in the lateral directions to the interfacial surface, blue points represent the normal direction, and green points represent the total values. Continuous lines represent the histograms for the interfacial chains of a wider thin layer than the layer at its tc containing 288 chains; the color code in continuous lines is the same as the lines of points. Vertical lines represent the mode of the distributions.
Figure 7. Histograms of Rg of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium at their tc at 673.15 K. Black and red points represent the components in the lateral directions to the interfacial surface, blue points represent the normal direction, and green points represent the total values. Continuous lines represent the histograms for the interfacial chains of a wider thin layer than the layer at its tc containing 288 chains; the color code in continuous lines is the same as the lines of points. Vertical lines represent the mode of the distributions.
Polymers 13 03515 g007
Figure 8. (a) t G d s (black circles) of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium as a function of T. Standard deviations are smaller than the symbols. The black line represents the adjustment to a quadratic function on T. (b) t i (black circles) as a function of T. The solid line represents the adjustment to a quadratic function on T. The bars represent standard deviations of the measurements. Red and blue lines in (a) represent t G d s t i and t G d s + t i , respectively.
Figure 8. (a) t G d s (black circles) of thin liquid layers of melted polyethylene chains (C200) in a liquid–vacuum equilibrium as a function of T. Standard deviations are smaller than the symbols. The black line represents the adjustment to a quadratic function on T. (b) t i (black circles) as a function of T. The solid line represents the adjustment to a quadratic function on T. The bars represent standard deviations of the measurements. Red and blue lines in (a) represent t G d s t i and t G d s + t i , respectively.
Polymers 13 03515 g008
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

González-Mijangos, J.A.; Lima, E.; Guerra-González, R.; Ramírez-Zavaleta, F.I.; Rivera, J.L. Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics. Polymers 2021, 13, 3515. https://doi.org/10.3390/polym13203515

AMA Style

González-Mijangos JA, Lima E, Guerra-González R, Ramírez-Zavaleta FI, Rivera JL. Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics. Polymers. 2021; 13(20):3515. https://doi.org/10.3390/polym13203515

Chicago/Turabian Style

González-Mijangos, José Antonio, Enrique Lima, Roberto Guerra-González, Fernando Iguazú Ramírez-Zavaleta, and José Luis Rivera. 2021. "Critical Thickness of Free-Standing Nanothin Films Made of Melted Polyethylene Chains via Molecular Dynamics" Polymers 13, no. 20: 3515. https://doi.org/10.3390/polym13203515

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