CO{}_{2} packing polymorphism under confinement in cylindrical nanopores

Co packing polymorphism under confinement in cylindrical nanopores

Ilaria Gimondi Thomas Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, UK.    Matteo Salvalaglio Thomas Young Centre and Department of Chemical Engineering, University College London, London WC1E 7JE, UK.
July 31, 2019

We investigate the effect of cylindrical nano-confinement on the phase behaviour of a rigid model of carbon dioxide using both molecular dynamics and well tempered metadynamics. To this aim we study a simplified pore model across a parameter space comprising pore diameter, CO-pore wall potential and CO density. In order to systematically identify ordering events within the pore model we devise a generally applicable approach based on the analysis of the distribution of intermolecular orientations. Our simulations suggest that, while confinement in nano-pores inhibits the formation of known crystal structures, it induces a remarkable variety of ordered packings unrelated to their bulk counterparts, and favours the establishment of short range order in the fluid phase. We summarise our findings by proposing a qualitative phase diagram for this model.

polymorphism, carbon dioxide, confinement, metadynamics, phase diagram

I Introduction

Confinement is known to play a role in the phase behaviour of molecular solids, most notably affecting polymorph selectionpolymorph_selection_2004 (); polymorph_control_2008 (). For instance, a paradigmatic example of the dramatic effects of confiment on the spatial arrangement of molecules is provided by water, which, as proven both experimentally and computationally, displays a counter-intuitively complex phase diagram under confinementManiwa2002 (); Striolo2005 (); Kyakuno2011 (); Algara-Siller2015 (); Soper2015 (); Agrawal2016 (); Chen2016 (); Chen2016a (). Understanding polymorphism in confined volumes is relevant both to describe natural processesBiomineralization2010 (); Biomineralization2013 () as well as for driving rational materials and process developmentconfined_surfaces_2005 (); pharmaceuticals_nanoconfinement_2007 (). Despite its importance a systematic understanding of confinement effects is still lacking.

Following up a recent work, in which we have investigated the thermodynamics and mechanism of phase transition between CO forms I and III in bulk, we set out to study the effect of confinement on CO condensed phases. We tackle this problem by carrying out a systematic analysis of CO phase behaviour confined in weakly interacting cylindrical nano pores. Our work has a two-fold aim: on the one hand understanding phase behaviour of confined CO is relevant due to its prominent role within the carbon cycle, and the surging needs for mitigating its emissions in atmosphere by implementing capture and storage technologies based on adsorption in porous solids Benson2005 (); Orr2009 (); macdowell2010overview (); Krevor2015 (). On the other hand due to its modest structural complexity accompanied with a rich phase diagram, CO represents a convenient model system to perform extensive sampling of polymorphic transitions at finite temperatureGimondi2017 () and gain insight on general aspects of molecular phase transitions under confinement.

Both experiments and theory highlight remarkable effects of confinements on CO phase behavior. For instance, the density of confined CO can significantly exceed the fluid phase bulk density Papadopoulos2001 (); Steriotis2004 (); Melnichenko2010 (); Rother2012 (); Sanghi2012 (); Elola2016 () reaching values comparable to that of solid CO and displaying signs of orientational correlations and structural rearrangementsSteriotis2004 (); Papadopoulos2001 ().

Several examples of molecular modelling studies can be found investigating the transport of CO in a confined fluid phase SlitClayPore2005 (); CO2inClayH2010 (); Sanghi2012 (), and the modelling literature characterising structural features of CO under confinement is limited. An insightful work from Elola and Rodriguez has recently proposed a detailed analysis of the CO fluid structure within silica pores, highlighting how confinement induces a layered arrangement in the fluid phase and determines a significant slowdown of both translational and rotational diffusion of CO moleculesElola2016 ().

In our work we aim complementing the state of the art by understanding whether confinement promotes or inhibits the formation of ordered phases. To this aim we systematically analyse the intermolecular structural organisation of CO confined in a simplified pore geometry consisting of a cylinder of Lennard-Jones (LJ) particles. We carry out a systematic study as a function of the cylinder radius (spanning from 1 to 5 nm), the nominal density of CO molecules (5 - 15 molecules/nm), and the interaction potential with the pore walls. We investigate this parameter space by systematically carrying out unbiased molecular dynamics simulations starting from liquid-like initial conditions, and by enhancing the exploration of ordered molecular packings with metadynamics.

Our analysis shows that, while the formation of known solid CO arrangements is inhibited, confinement unveils a rather complex behaviour. Depending on the location in parameter space, confined CO can either approach a completely disordered fluid state, exhibit short-range orientation correlations or arrange into long-range ordered packings. By systematically analysing relative orientations distributions we have systematically detected ordering events from molecular trajectories and classified ordered configurations.

This paper is organized as follows: firstly, in section 2, we introduce the model system, methods, parameters space, and the analysis tools developed to characterise ordering transitions. In section 3 we present and discuss our results, and finally in section 4 we summarise our findings by proposing a qualitative phase diagram in the space defined the nominal density of the pore and its normalised radius.

Ii Methods

Figure 1: Simulations setup. (a) Representation of the LJ potential for the interaction between the pore wall and the carbon atom of CO we report three pore sizes (d = 1, 2, 5 nm), for which the zero on the axis represents the cylinder axis and the walls are located at 0.5, 1 and 2.5 nm, respectively; is set to 0.34 nm. Arrows highlight (b) Intermolecular LJ potential characterising CC, CO and OO non bonded interactions. The axis represents the interatomic distance; the black arrow indicates . (c) Representation of the cross-section of the initial configuration for simulations with = 1, 2, and 5 nm.

To investigate the effect of confinement on liquid carbon dioxide we have employed both standard molecular dynamics (MD)Tuckerman2010 (); Frenkel2002 (); Palmer2015 () and well-tempered metadynamics (WTMetaD)Barducci2008 (). For a detailed description of metadynamics we refer to Barducci et al.Barducci2008 (); Barducci2011 (), and Valsson et al.Valsson2016 (), and for a brief overview of its applications in crystallisation studies to Giberti et al.Giberti2015 (). In the following we report the setup of the molecular model investigated, the details of MD simulations and WTmetaD simulations setup, and the analysis approach implemented to detect and characterise the emergence of an ordered phase from an initial liquid state.

Co Potential.

In order to model CO we employ the Transferable Potentials for Phase Equilibria (TraPPE) force fieldPOTOFF1999 (); Potoff2001 () (Table 1) with the introduction of dummy atoms to maintain molecules rigid with a 180 angle without inducing instabilitySanghi2012 (). This choice is consistent with our previous work on carbon dioxide polymorphism at high pressure Gimondi2017 (), and allows for a direct comparison with results obtained in bulk.

m m
[nm] [nm] [kJ/mol]
12 16 0.280 0.305 0.224
q q
[kJ/mol] [e] [e] [Å] []
0.657 0.70 -0.35 1.160 180
Table 1: Parameters for the TraPPE force field

Pore Model.

To develop a simplified model to study confinement in a weakly interacting porous medium, we build a single-wall cylindrical nanopore of LJ particles. We consider different pore diameters (1, 1.3, 2, 5 nm) with a constant wall density of 33.104 atoms/nm and a height of 10 nm for all diameters except the case with diameter 5 nm, where the height has been limited to 5 nm as well for the sake of computational efficiency. The cylinder axis aligns with z. We also investigate the effect induced by variations in the wall potential induced by keeping constant while varying the , in the interval between 0.253 to 0.405 nm.

Simulation set up.

For our simulations we choose temperature conditions of 323 K and nominal carbon dioxide densities compatible with supercritical COKrevor2015 (); Melnichenko2010 (); Sanghi2012 (). Firstly, we run NVT molecular dynamics simulations of the described system, with Bussi-Donadio-ParrinelloBussi2007 () thermostat. We apply periodic boundary conditions (pbc) to ensure continuity along z, i.e. the cylinder axis. To avoid overlapping effects in the radial direction from periodic images, the simulation box is much larger (around 10 times) than the pore diameter. Lennard-Jones potential with Lorentz-Berthelot combination rules and long-range corrections is employed. Typical MD simulations are 20 ns-long, starting from disordered configurations equilibrated for 100 ps.

We use the same setup to carry out WTMetaD simulations aimed at exploring the space of accessible packings beyond the timescale limit of unbiased MD. As collective variable we employ an order parameter, hereafter indicated with , that expresses the degree of crystallinity of a system based on the local environment around each molecule. A complete description of the mathematical formulation of has been reported by Giberti et al.Giberti2015_2 (). Explorative WTMetaD simuations are carried out biasing two distinct formulations of the order parameter: , tuned to capture the molecular arrangement of CO form I, which we recently developed to investigate CO polymorphism at high pressureGimondi2017 (); and , based instead on the characteristic angles and local density of the most abundant ordered structure observed under confinement (configuration B in Table 4). For details on the parameters used to define both and refer to Table 2. Further details about the WTMetaD set up are reported in the Supporting Information.

[] [] [] n [-] r [Å]
70.47 108.86 14.32 5 4
46.18 133.82 5.21 3 4
Table 2: Tuning of the -order parameters. The table reports , its supplementary , the associated width of the Gaussian, , which for symmetry reasons is the same for both angles. The cut-off values for the number of neighbours and the coordination shell are presented as well.

Parameter space representation.

In order to represent our results throughout parameter space in a rational and efficient form we define a single adimensional radius, , which accounts for variations in both LJ potential () and pore size and quantifies the void space available within the poreKyakuno2011 (). The adimensional radius is defined by the following expression:


where represents the distance from the wall to the minimum of the C-wall potential (as in Figure 1(a), for = 0.34 nm) and is the position of the minimum of the C-C interaction. Values of for all simulations performed in this work are reported in Table 3. The second parameter is the nominal density of CO molecules within the pore , computed accounting for the volume of the cylinder of diameter .

Figure 2: (a) The distribution of relative orientations is built considering all the angles between an assigned molecule and its nearest neighbours, falling within a sphere of radius represented in transparent green. The relative orientation is computed as the elevation angle between the axis of two CO molecules. (b) Characteristic angle distribution for bulk phases, in particular melt (green) and phase I (blue); the analytical expression for the random distribution in melt is also reported in black.
Figure 3: Probability density profiles associated with the position of CO molecules in the radial direction of the pore, for the range of densities investigated at different values of . In particular, (a) refers to = 2.075, and (b) to = 0.484. Along the abscissa, i.e. the pore radius, zero corresponds to the cylinder axis, while the maximum value reported to the radius of the pore. For the case in (b), (e) further analyses the height of the adsorbed layer peak and the number of peaks as a function of . For the case in (a), instead, (c) reports snapshots of hollow to filled structures for increasing , from 5.41 to 8.15 to 12.73 molecules/nm. More in detail, (d) presents the probability density of the angle as a function of the radial distance, for = 2.075 and = 12.73 molecules/nm. The angle , reported as (cos), is evaluated between CO molecular axis and the vector normal to the pore walls; to enhance the differences on the plot the probability density is represented by its logarithm. The -axis convention is the same as before, with cylinder axis in zero.
d [nm]
1 1.3 2 5
[nm] 0.253 0.639 1.116 2.230 7.003
0.315 0.528 1.006 2.119 6.892
0.340 0.484 0.961 2.075 6.847
0.372 0.427 0.904 2.017 6.790
0.405 0.368 0.845 1.959 6.731
Table 3: Value of for the combinations of d and investigated.

Systematic detection of ordered arrangements.

In the definition of ordered parameters such as the identification of characteristic relative orientations is keySalvalaglio2012 (); Giberti2015 (); Gimondi2017 (). This is typically based on the analysis of the distribution of relative orientations that characterises a relevant structure, as extensively discussed for bulk CO in Ref. Gimondi2017 (). Indeed the relative orientation between molecular axes of neighbouring molecules represents a fingerprint of each arrangement and in principle allows to distinguish not only between liquid and solid, but also among different solid structures. An example of such distributions for bulk dry ice and liquid CO is reported in Figure 2. In the following we build on this observation to systematically detect ordering phenomena in within molecular trajectories.

As reported in Figure 2, crystal and liquid carbon dioxide display significantly different characteristic orientation distributions. In order to quantitatively capture such difference we apply the Bhattacharyya distanceBhattacharyya1943 (), a metric that quantifies the dissimilarity of two probability densities. The Bhattacharyya distance D between two angle distributions and is defined as:


where is the Bhattacharyya coefficient, which measures the overlap between and defined as:


In our analysis we take advantage of the fact that a completely random arrangement of linear molecules displays a relative orientation probability density characterised by the functional form . For a derivation of the reference distribution functional form see the Supporting Information. As shown in Figure  2 the distribution of angles in the bulk of liquid CO closely resembles the theoretical random distribution .

Given a molecular trajectory, we can therefore define a time dependent measure of the deviation from an ideal random packing of CO molecules as:


where is the instantaneous probability density of relative orientations.

The construction of is straightforward for carbon dioxide, due to its linear, symmetrical geometry. For each molecule we evaluate the polar angle between its molecular axis and the analogous vector in each one of its nearest neighbours, i.e. molecules within a cut-off radius of 4 Å for ordered phases and 5 Å for the liquid as shown in Figure 2(a).


Pores are built using an in-house FORTRAN code, while we use PackmolTruhlar2009 () to insert carbon dioxide molecules at random positions and orientations in order to avoid any creation of patterns. For MD and WTMetaD simulations we employ Gromacs 5.2.1Abraham2015 () and Plumed 2.2Tribello2014 (). Post-processing of the results is carried on with MATLAB (R2015a), Visual Molecular Dynamics (VMD)HUMP96 () and Plumed 2.2Tribello2014 ().

Iii Results and Discussion

Figure 4: Bhattacharyya analysis of an ordering phenomenon. (a) Temporal evolution of the angle distribution, , for = 0.484 - = 11.78 molecules/nm, where colour dark blue to light blue represents time increment, while the black distribution is the probability density for melt, employed as reference for evaluating D. (b) Bhattacharyya distance, D, as a function of time for an example of ordering (violet, - as (a)) and fluid (green, = 0.961 - = 5.41 molecules/nm conditions; the probability density of D is also reported in the form of a histogram.

In this section we present and discuss the results of our molecular dynamics investigation of the effect of confinement on CO ordering. Preliminary simulations were carried out to investigate the equilibrium distribution of a single CO molecule in the presence of a cylindrical LJ wall and are discussed in Supporting Information. In the following we begin by discussing the arrangement of CO molecules within the pore as a function of pore size and density, and we continue by reporting an analysis of the ordered arrangements identified within our parameter space.

Molecular arrangement within the pore.

We begin by analysing the arrangement of CO molecules in the pores at increasing . In the case of low densities we observe the formation of a single layer of adsorbed CO molecules in contact with the pore surface associated with a hollow structure. With increasing larger pores (d = 2 and 5 nm), can accommodate a second layer of CO, and then display bulk-like filling as shown in Figure 3(a,c). This finding is in agreement with the results obtained by Elola and RodriguezElola2016 () for CO confined in silica pores.

Furthermore, we analyse CO molecular orientation within the pore by evaluating the angle between the molecular axis and the normal to the pore surface. As shown in Figure 3 (d), we can clearly note the high probability regions in the r, plane, corresponding to individual layers. We note that in both layers is narrowly distributed around 90. This indicates that CO molecules arrange parallel to the pore surface, in agreement with observations reported in the literatureVishnyakov1999 (); Steriotis2004 (); Elola2016 (); Sanghi2012 ().

In smaller pores (d = 1 and 1.3 nm) a layered structure is not observed and the decrease in height of the dominant peak (peak 1 in Figure 3(b)), accompanied by the appearance of a secondary peak (peak 2 in Figure 3(b)) is instead associated with a reorganisation of the CO packing within the pore (see Figure 3(b,e)). For instance, we note that at d=1.3 nm, CO arranges as a chain along the z axis for densities equal or above 9.87 molecules/nm.

Detecting emerging order.

After discussing the position and alignment of CO molecules with respect to the pore wall, we delve into the analysis of intermolecular arrangement. In particular, we analyse the arrangement of the local environment around each CO molecule by monitoring the distribution of relative orientations. As mentioned in section 2, this information constitutes a structural fingerprint and can discriminates not only between ordered and disordered arrangements but also among different ordered packings.

As discussed in section 2, to detect the unfolding of ordering events we compute the Bhattacharyya distance, D with respect to an ideal random packing for each MD trajectory. When ordered arrangements nucleate within a pore, the D metric displays a sharp increase. An example of this behaviour in a typical trajectory is reported in Figure 4, where in (b) we compare the trend of D(t) for a trajectory undergoing an ordering transition (violet) and a trajectory that remains stable in the fluid state (green). In Figure 4(a) the evolution of the instantaneous intermolecular angle distribution is reported together with the corresponding to an ideal random packing.

Figure 5: Radial pair distribution function, g(r) (black) and angle distribution (blue) for illustrative cases of liquid (a-b), transition (c-d), and ordered structures (e-f). In particular, (a) refers to = 2.075 - = 5.41 molecules/nm, (b) to = 2.075 - = 12.73 molecules/nm, (c) to = 0.484 - = 5.41 molecules/nm, (d) to = 0.484 - = 8.15 molecules/nm, (e) to = 0.484 - = 11.78 molecules/nm, and (f) to = 0.484 - = 12.73 molecules/nm. Green areas on the phase diagram refer to fluid state, violet to ordered, while violet-to-green colour to transition structures

Mapping phase behaviour in the parameter space.

The systematic analysis of the deviation from an ideal random packing and the contextual detection of the appearance of ordered phases from a liquid configuration allows to clearly map the qualitative behaviour of confined CO within the parameter space. In particular we can identify areas where CO under confinement spontaneously evolves to ordered packings within the characteristic timescale of in MD simulations. To further analyse our finding, we complement the analysis of the angle distribution by reporting the C-C radial pair distribution function g(r). Inspecting both g(r) and relative orientation distribution we can identify three regions in parameter space: a disordered region (Figure 5 green area, examples (a-b)), an ordered region (Figure 5 violet area, examples (e-f)), and a transition region between them (Figure 5 green-to-violet area, examples (c-d)). We note that in the disordered region, at large and low density, the distribution of relative orientations approaches the ideal random distribution and the g(r) displays the short-range order typical of a liquid. On the contrary, in the ordered region both the angle distribution and the g(r), deviate from the typical liquid behaviour denoting long-range order through the presence of characteristic peaks in both orientation distribution and g(r). Finally, the transition region denotes a peculiar behaviour. In this region of parameter space, while g(r) maintains a close resemblance to that of a liquid phase, the distribution of relative orientations displays significant deviations from an ideal random packing suggesting the development of correlation in the relative orientation of CO molecules.

From the analysis of the transitions obtained in the structured region of parameter space we notice a finer grained level of complexity emerging. As displayed by the qualitative comparison of both the distribution of relative orientations and the g(r) reported in examples Figure 5(e) and (f) multiple ordered packings can be obtained.

Stable ordered packings.

We now move on to present a detailed analysis of the ordered structures encountered within the ordered region of the phase diagram. We shall note that, in order to confirm the finding emerging from unbiased simulations, we have applied WTMetaD simulations to further sample the configuration space of confined CO. This allows to verify that in the disordered or transition regions ordered arrangements cannot be found or, when found do not states with a finite lifetime and spontaneously revert to liquid in unbiased MD.
Nevertheless, in the region of the phase diagram where confinement-induced order emerges, we identify four different stable arrangements. We assign to such arrangements a label progressing from A to D, based on their structural features. As we can see in Table 4, progressing from A to D these structures present an increasing number of molecules in their repeating unit. Moreover, the same classification applies when analysing the geometry of the apparent ring that characterises the horizontal cross-section of the pore (see Table 4 ). The number of ring members is 3 in A, 4 in B, and 6 in both C and D, with D displaying one molecule in the middle of the pore oriented parallel to the pore longitudinal axis.

Label Snaphots Relative Orientation Distribution
Table 4: Molecular representation of the arrangement and characteristic angles distribution for the ordered structures identified. Such structures are named A to D according to the apparent number of ring members when observed from the top. The angles distributions are obtained at = 0.484 - = 9.87, at = 0.484 - = 11.78, at = 0.484 - = 12.73, at = 0.961 - = 12.73, for A, B, C, and D, respectively.
Figure 6: Qualitative phase diagram of CO confined in weakly attractive cylindrical nanopores. In green is represented the region of phase space characterised by disordered configurations, while we report in violet the region of phase space where ordered structures are found. Configurations yielding one or more amongst packing A, B, C, and D are appropriately marked. Qualitative phase boundaries have been reported in order to highlight how simulations yielding similar structures tend to cluster in parameter space, suggesting that the spontaneous ordering process observed from MD reflects the underlying thermodynamic stability of the packings.

As far as the packing of these four ordered structures is concerned, it can be seen that the characteristic angle distribution represents a clear fingerprint of their molecular arrangement (Table 4), while the coordination number appears to be a function of the total density inside the pore and cannot clearly discriminate between different arrangements. In particular, the orientation distribution of B presents well-defined peaks in = 46.2 and symmetrically in = 133.8, while C a very sharp single one in = 90. Structure A, instead, appears to be an intermediate between B and C both from the inspection of the coordinates and from its characteristic angle distribution. Indeed, while the most densely populated orientations closely resemble those of B, the peak at 90 is more pronounced than in B. Form D present instead wider peaks in = 57.9 and the additional peak = 122.1.

Interestingly, three out of four ordered structures, more precisely A, B, and D, form spontaneously in unbiased MD simulations; phase C, instead, is generated from metadynamics simulations performed in conditions that spontaneously yielded other ordered arrangements (violet area in the / phase diagram in Figures 5 and6); structure C is confirmed to be a metastable state for the system possessing a finite lifetime through 20 ns long unbiased MD simulations, initialised in configuration C, which do not display any sign of instability transition.
It should be noted that the exploration carried out with WTMetaD yields additional ordered arrangements, which, however, do not survive unbiased MD and are therefored considered to be unstable (see SI).

In addition metadynamics further confirms that for - conditions that induce spontaneous ordering, the disordered liquid state is strongly disfavoured with respect to the ordered phases. Indeed even in extended WTMetaD simulations, characterised by the deposition of significant bias, no melting event can be observed, while transitions between ordered arrangements are routinely detected. This finding is consistent irrespective of the set of CVs used to enhance the exploration and suggests that the energetic barrier associated with melting is higher than the barriers associated with solid-solid interconversion.

Furthermore, the conversion from hollow (A, B, C) to filled (D) ordered structures and vice versa is not sampled during WTMetaD simulations; this transition appears unlikely in the conditions investigated since, as discussed before, different values of can accommodate different numbers of CO layers.

On the other hand, transition between phase A, B and C is possible and consistently sampled. More precisely, it takes place though a gradual reorganisation of the layers that starts in a localised region and then spreads through the entire pore length. Interestingly, WTMetaD allows to sample mixed arrangements of different packings, (e.g. half ordered as B and half as C) which reveal to be stable in subsequent MD simulation, confirming that solid-solid interconversion is an activated event. Moreover, in these structures it is common to find boundary regions in which the order is not well-defined, effectively acting as defects within the pore.

Iv Conclusions

In this work we present a systematic investigation of the effects of confinement in cylindrical nanopores on CO, as a function of pore size, wall potential and CO density. To systematically detect the formation of ordered packings we have devised a general analysis approach based on the calculation of the time-dependent Bhattacharyya distance between a reference probability density and the instantaneous probability density of the relative orientation of nearest neighbours. Here we have adopted as a reference an ideal random packing, which closely approximates the distribution of relative orientations in the bulk of the liquid phase. This analysis approach is general and we anticipate its applicability to a much wider range of phase transitions in molecular solids.

We find that cylindrical confinement induces polymorph selection by inhibiting the nucleation of known solid forms of CO while at the same time inducing the organisation of CO into a series of distinct ordered arrangements, which do not resemble any of the known CO polymorphs. Using a combination of unbiased MD and WTMetaD we identify four ordered packing that are stable in the region of parameter space corresponding to small pore radii and large densities. These configurations have been labelled A, B, C and D following the progression in the number of CO molecules visible in a cross-section of a pore. Interestingly, none of these configurations corresponds to known bulk structures. It should be noted that, while A, B, and D emerge spontaneously from unbiased MD, C is found from the enhanced exploration of the configuration space achieved with WTmetaD.

In Figure 6 we report a qualitative phase diagram that summarizes the identified structures and the regions of parameter space where such structures have obtained. Markers point the location in parameter space where specific packings have been observed. Cases in which multiple configurations are indicated on the phase diagram represent points in which multiple arrangements have been observed either with unbiased MD or through WTMetaD. We observe that, in the ordered region of the phase diagram WTMetaD was quite efficient in promoting transitions between ordered packings, however it was unable to induce melting. Despite its qualitative character, this observation suggests that within the ordered region of the phase diagram solid-solid interconversion is energetically favoured with respect to liquid mediated transitions.

Finally we observe that, even in the disordered region confinement induces a short-range organisation of the liquid phase, in which the distribution of relative orientations departs from an ideal random packing, as shown in Figure 5. At large volumes and low densities the bulk behaviour is recovered as expected.

To conclude we note how such a simple model system such as CO confined in a simple nanometric pore model unveils remarkable structural complexity, showcasing the profound effect of confinement on the phase behaviour of molecular solids. This work paves the way for a systematic assessment of CO packing polymorphism in systems characterised by realistic pore geometries and surface chemistry.

Conflict of interest

There are no conflicts to declare.


The authors acknowledge EPSRC (Engineering and Physical Sciences Research Council) for PhD scholarship, and UCL Legion High Performance Computing Facility for access to Legion@UCL and associated support services, in the completion of this work.

Supplementary Material

See supplementary material for single molecule probability distribution in the model pore, derivation of the ideal random distribution , details of the setup and results of explorative metadynamics simulations, additional analysis of the radial density profiles. .



  • [1] Jeong-Myeong Ha, Johanna H. Wolf, Marc A. Hillmyer, and Michael D. Ward. Polymorph selectivity under nanoscopic confinement. J. Am. Chem. Soc., 126(11):3382–3383, 2004. PMID: 15025439.
  • [2] Antonio Llinàs and Jonathan M. Goodman. Polymorph control: past, present and future. Drug. Discov. Today, 13(5):198 – 210, 2008.
  • [3] Yutaka Maniwa, Hiromichi Kataura, Masatoshi Abe, Shinzo Suzuki, Yohji Achiba, Hiroshi Kira, and Kazuyuki Matsuda. Phase Transition in Confined Water Inside Carbon Nanotubes. J. Phys. Soc. Jpn., 71(12):2863–2866, dec 2002.
  • [4] A. Striolo, A. A. Chialvo, K. E. Gubbins, and P. T. Cummings. Water in carbon nanotubes: Adsorption isotherms and thermodynamic properties from molecular simulation. J. Chem. Phys., 122(23):234712, jun 2005.
  • [5] Haruka Kyakuno, Kazuyuki Matsuda, Hitomi Yahiro, Yu Inami, Tomoko Fukuoka, Yasumitsu Miyata, Kazuhiro Yanagi, Yutaka Maniwa, Hiromichi Kataura, Takeshi Saito, Motoo Yumura, and Sumio Iijima. Confined water inside single-walled carbon nanotubes: Global phase diagram and effect of finite length. J. Chem. Phys., 134(24):244501, jun 2011.
  • [6] G Algara-Siller, O Lehtinen, F C Wang, R R Nair, U Kaiser, H a Wu, a K Geim, and I V Grigorieva. Square ice in graphene nanocapillaries. Nature, 519(7544):443–5, 2015.
  • [7] Alan K. Soper. Physical chemistry: Square ice in a graphene sandwich. Nature, 519(7544):417–418, mar 2015.
  • [8] Kumar Varoon Agrawal, Steven Shimizu, Lee W. Drahushuk, Daniel Kilcoyne, and Michael S. Strano. Observation of extreme phase transition temperatures of water confined inside isolated carbon nanotubes. Nat. Nanotechnol., 12(3):267–273, nov 2016.
  • [9] Ji Chen, Georg Schusteritsch, Chris J. Pickard, Christoph G. Salzmann, and Angelos Michaelides. Two Dimensional Ice from First Principles: Structures and Phase Transitions. Phys. Rev. Lett., 116(2):025501, jan 2016.
  • [10] Ji Chen, Andrea Zen, Jan Gerit Brandenburg, Dario Alfè, and Angelos Michaelides. Evidence for stable square ice from quantum Monte Carlo. Phys. Rev. B, 94(22):220102, dec 2016.
  • [11] Christopher J. Stephens, Sophie F. Ladden, Fiona C. Meldrum, and Hugo K. Christenson. Amorphous calcium carbonate is stabilized in confinement. Adv. Funct. Mater., 20(13):2108–2115, 2010.
  • [12] Yun-Wei Wang, Hugo K. Christenson, and Fiona C. Meldrum. Confinement leads to control over calcium sulfate polymorph. Adv. Funct. Mater., 23(45):5615–5623, 2013.
  • [13] Alfred Y. Lee, In Sung Lee, Severine S. Dette, Jana Boerner, and Allan S. Myerson. Crystallization on confined engineered surfaces:? a method to control crystal size and generate different polymorphs. J. Am. Chem. Soc., 127(43):14982–14983, 2005. PMID: 16248610.
  • [14] M. Beiner, Rengarajan, S. Pankaj, D. Enke, and M. Steinhart. Manipulating the crystalline state of pharmaceuticals by nanoconfinement. Nano Lett., 7(5):1381–1385, 2007. PMID: 17439189.
  • [15] Sally M. Benson. Overview of Geologic Storage of CO2. In Carbon Dioxide Capture for Storage in Deep Geologic Formations, volume 2, pages 665–672. Elsevier, 2005.
  • [16] F M Orr. Onshore Geologic Storage of CO2. Science, 325(5948):1656–1658, sep 2009.
  • [17] Niall MacDowell, Nick Florin, Antoine Buchard, Jason Hallett, Amparo Galindo, George Jackson, Claire S Adjiman, Charlotte K Williams, Nilay Shah, and Paul Fennell. An overview of co 2 capture technologies. Energ. Environ. Sci., 3(11):1645–1669, 2010.
  • [18] Samuel Krevor, Martin J. Blunt, Sally M. Benson, Christopher H. Pentland, Catriona Reynolds, Ali Al-Menhali, and Ben Niu. Capillary trapping for geologic carbon dioxide storage ? From pore scale physics to field scale implications. Int. J. Greenh. Gas Con., 40:221–237, sep 2015.
  • [19] Ilaria Gimondi and Matteo Salvalaglio. Co2 packing polymorphism under pressure: mechanism and thermodynamics of the i-iii polymorphic transition. J. Chem. Phys., 2017.
  • [20] George K. Papadopoulos. Influence of orientational ordering transition on diffusion of carbon dioxide in carbon nanopores. J. Chem. Phys., 114(18):8139–8144, may 2001.
  • [21] Th.A Steriotis, K.L Stefanopoulos, N.K Kanellopoulos, A.Ch Mitropoulos, and A. Hoser. The structure of adsorbed CO2 in carbon nanopores: a neutron diffraction study. Colloid. Surface. A, 241(1-3):239–244, jul 2004.
  • [22] Yuri B. Melnichenko, H. Mayama, G. Cheng, and T. Blach. Monitoring Phase Behavior of Sub- and Supercritical CO 2 Confined in Porous Fractal Silica with 85% Porosity. Langmuir, 26(9):6374–6379, may 2010.
  • [23] Gernot Rother, Elizabeth G. Krukowski, Dirk Wallacher, Nico Grimm, Robert J. Bodnar, and David R. Cole. Pore Size Effects on the Sorption of Supercritical CO 2 in Mesoporous CPG-10 Silica. J. Phys. Chem. C, 116(1):917–922, jan 2012.
  • [24] T Sanghi and N R Aluru. Coarse-grained potential models for structural prediction of carbon dioxide (CO 2 ) in confined environments. J. Chem. Phys., 136(2):024102, jan 2012.
  • [25] M. Dolores Elola and Javier Rodriguez. Excess Sorption of Supercritical CO 2 within Cylindrical Silica Nanopores. J. Phys. Chem. C, 120(2):1262–1269, jan 2016.
  • [26] Xiaoning Yang and Cuijuan Zhang. Structure and diffusion behavior of dense carbon dioxide fluid in clay-like slit pores by molecular dynamics simulation. Chem. Phys. Lett., 407(4):427 – 432, 2005.
  • [27] Alexandru Botan, Benjamin Rotenberg, Virginie Marry, Pierre Turq, and Benoît Noetinger. Carbon dioxide in montmorillonite clay hydrates: Thermodynamics, structure, and transport from molecular simulation. J. Phys. Chem. C, 114(35):14962–14969, 2010.
  • [28] Mark E. Tuckerman. Statistical Mechanics: Theory and Molecular Simulation. OXFORD UNIVERSITY PRESS, 2010.
  • [29] Daan Frenkel and Berend Smit. Understanding Molecular Simulation - From Algorithms to Applications. ACADEMIC PRESS, 2nd edition, 2002.
  • [30] Jeremy C. Palmer and Pablo G. Debenedetti. Recent advances in molecular simulation: A chemical engineering perspective. AIChE J., 61(2):370–383, feb 2015.
  • [31] Alessandro Barducci, Giovanni Bussi, and Michele Parrinello. Well-tempered metadynamics: A smoothly converging and tunable free-energy method. Phys. Rev. Lett., 100(2):1–4, 2008.
  • [32] Alessandro Barducci, Massimiliano Bonomi, and Michele Parrinello. Metadynamics. Wires Comput. Mol. Sci., 1(5):826–843, sep 2011.
  • [33] Omar Valsson, Pratyush Tiwary, and Michele Parrinello. Enhancing Important Fluctuations: Rare Events and Metadynamics from a Conceptual Viewpoint. Annu. Rev. Phys. Chem., 67(1):159–184, may 2016.
  • [34] Federico Giberti, Matteo Salvalaglio, and Michele Parrinello. Metadynamics studies of crystal nucleation. IUCrJ, 2(2):256–266, mar 2015.
  • [35] J. J. Potoff, J. R. Errington, and A. Z. Panagiotopoulos. Molecular simulation of phase equilibria for mixtures of polar and non-polar components. Mol. Phys., 97(10):1073–1083, nov 1999.
  • [36] Jeffrey J Potoff and J Ilja Siepmann. Vapor–liquid equilibria of mixtures containing alkanes, carbon dioxide, and nitrogen. AIChE J., 47(7):1676–1682, jul 2001.
  • [37] Giovanni Bussi, Davide Donadio, and Michele Parrinello. Canonical sampling through velocity rescaling. J. Chem. Phys., 126(1):014101, 2007.
  • [38] Federico Giberti, Matteo Salvalaglio, Marco Mazzotti, and Michele Parrinello. Insight into the nucleation of urea crystals from the melt. Chem. Eng. Sci., 121:51–59, jan 2015.
  • [39] Matteo Salvalaglio, Thomas Vetter, Federico Giberti, Marco Mazzotti, and Michele Parrinello. Uncovering Molecular Details of Urea Crystal Growth in the Presence of Additives. J. Ame. Chem. Soc., 134(41):17221–17233, oct 2012.
  • [40] A Bhattacharyya. No Title. Bull. Calcutta Math. Soc., 35:99–109, 1943.
  • [41] Donald G. Truhlar. Valence Bond Theory for Chemical Dynamics. J. Comput. Chem., 28(1):73–86, 2009.
  • [42] M.J. Abraham, D. van der Spoel, E. Lindahl, B. Hess, and the GROMACS development Team. GROMACS User Manual version 5.1.1., 2015.
  • [43] Gareth A. Tribello, Massimiliano Bonomi, Davide Branduardi, Carlo Camilloni, and Giovanni Bussi. PLUMED 2: New feathers for an old bird. Comput. Phys. Comm., 185(2):604–613, feb 2014.
  • [44] William Humphrey, Andrew Dalke, and Klaus Schulten. VMD – Visual Molecular Dynamics. J. Mol. Graphics, 14:33–38, 1996.
  • [45] Aleksey Vishnyakov, Peter I. Ravikovitch, and Alexander V. Neimark. Molecular Level Models for CO 2 Sorption in Nanopores. Langmuir, 15(25):8736–8742, dec 1999.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description