Correlating Substrate Reactivity at Electrified Interfaces with the Electrolyte Structure in Synthetically Relevant Organic Solvent/Water Mixtures.

Optimizing electrosynthetic reactions requires fine tuning of a vast chemical space, including charge transfer at electrocatalyst/electrode surfaces, engineering of mass transport limitations, and complex interactions of reactants and products with their environment. Hybrid electrolytes, in which supporting salt ions and substrates are dissolved in a binary mixture of organic solvent and water, represent a new piece of this complex puzzle as they offer a unique opportunity to harness water as the oxygen or proton source in electrosynthesis. In this work, we demonstrate that modulating water-organic solvent interactions drastically impacts the solvation properties of hybrid electrolytes. Combining various spectroscopies with synchrotron small-angle X-ray scattering (SAXS) and force field-based molecular dynamics (MD) simulations, we show that the size and composition of aqueous domains forming in hybrid electrolytes can be controlled. We demonstrate that water is more reactive for the hydrogen evolution reaction (HER) in aqueous domains than when strongly interacting with solvent molecules, which originates from a change in reaction kinetics rather than a thermodynamic effect. We exemplify novel opportunities arising from this new knowledge for optimizing electrosynthetic reactions in hybrid electrolytes. For reactions proceeding first via the activation of water, fine tuning of aqueous domains impacts the kinetics and potentially the selectivity of the reaction. Instead, for organic substrates reacting prior to water, aqueous domains have no impact on the reaction kinetics, while selectivity may be affected. We believe that such a fine comprehension of solvation properties of hybrid electrolytes can be transposed to numerous electrosynthetic reactions.

Crystalline restacking of 2D-materials from their nanosheets suspensions.

We report here the highly ordered restacking of the layered phosphatoantimonic dielectric materials H3(1-x)M3xSb3P2O14, (where M = Li, Na, K, Rb, Cs and 0 ≤ x ≤ 1), from their nanosheets dispersed in colloidal suspension, induced by a simple pH change using alkaline bases. H3Sb3P2O14 aqueous suspensions are some of the rare examples of colloidal suspensions based on 2D materials exhibiting a lamellar liquid crystalline phase. Because the lamellar period can reach several hundred nanometers, the suspensions show vivid structural colors and because these colors are sensitive to various chemicals, the suspensions can be used as sensors. The structures of the lamellar liquid crystalline phase and the restacked phase have been studied by X-ray scattering (small and wide angle), which has followed the dependence of the lamellar/restacked phase equilibrium on the cation exchange rate, x. The X-ray diffraction pattern of the restacked phase is almost identical to that of the M3Sb3P2O14 crystalline phase, showing that the restacking is highly accurate and avoids the turbostratic disorder of the nanosheets classically observed in nanosheet stacking of other 2D materials. Strikingly, the restacking process exhibits features highly reminiscent of a first-order phase transition, with the existence of a phase coexistence region where both ∼1 nm (interlayer spacing of the restacked phase) and ∼120 nm lamellar periods can be observed simultaneously. Furthermore, this first-order phase transition is well described theoretically by incorporating a Lennard-Jones-type lamellar interaction potential into an entropy-based statistical physics model of the lamellar phase of nanosheets. Our work shows that the precise cation exchange produced at room temperature by a classical neutralization reaction using alkaline bases leads to a crystal-like restacking of the exfoliated free Sb3P2O143- nanosheets from suspension, avoiding the turbostratic disorder typical of van der Waals 2D materials, which is detrimental to the controlled deposition of nanosheets into complex integrated electronic, spintronic, photonic or quantum structures.

Electrical surface properties of nanoporous alumina membranes: influence of nanochannels' curvature, roughness and composition studied via electrokinetic experiments.

Among classical nanoporous oxide membranes, anodic aluminum oxide (AAO) membranes, made of non-connected, parallel and ordered nanochannels, are very interesting nanoporous model systems widely used for multiple applications. Since most of these applications involve local phenomena at the nanochannel surface, the fine description of the electrical surface behavior in aqueous solution is thus of primordial interest. Here, we use an original experimental approach combining several electrokinetic techniques (tangential and transverse streaming potential as well as electrophoretic mobility experiments) to measure the ζ-potential and determine the surface isoelectric points (IEPs) of several AAOs having different characteristic sizes and compositions. Using such an approach, all the different surfaces available in AAOs can be probed: outer surfaces (top and bottom planes), pore wall surfaces (i.e., inner surfaces) and surfaces created by the grinding of the AAOs. We find clear IEP differences between the outer, pore wall and ground surfaces and discuss these in terms of nanochannel and surface morphology (curvature and roughness) and of modifications of the chemical environment of the surface hydroxyl groups. These results highlight the heterogeneities between the different surfaces of these AAO membranes and emphasize the necessity to combine complementary electrokinetic techniques to properly understand the material, an approach which can be extended to many nanoporous systems.

Rheo-SAXS characterization of lead-treated oils: Understanding the influence of lead driers on artistic oil paint's flow properties.

From the 15th century onwards, painters began to treat their oils with lead compounds before grinding them with pigments. Such a treatment induces the partial hydrolysis of the oil triglycerides and the formation of lead soaps, which significantly modify the rheological properties of the oil paint. Organization at the supramolecular scale is thus expected to explain these macroscopic changes. Synchrotron Rheo-SAXS (Small Angle X-ray Scattering) measurements were carried out on lead-treated oils, with different lead contents. We can now propose a full picture of the relationship between structure and rheological properties of historical saponified oils. At rest, lead soaps in oil are organized as lamellar phases with a characteristic period of 50 Å. Under shear, the loss of viscoelastic properties can be linked to the modification of this organization. Continuous shear resulted in a preferential and reversible orientation of the lamellar domains which increased with the concentration of lead soaps. The parallel orientation predominates over the entire shear range (0-1000 s-1). Conversely, oscillatory shear coiled the lamellae into cylinders that oriented themselves vertically in the rheometer cell. This is the first report of such a vertical cylindrical structure obtained under shear from lamellae.

Do Aqueous Suspensions of Smectite Clays Form a Smectic Liquid-Crystalline Phase?

Bottom-up strategies for the production of well-defined nanostructures often rely on the self-assembly of anisotropic colloidal particles (nanowires and nanosheets). These building blocks can be obtained by delamination in a solvent of low-dimensionality crystallites. To optimize particle availability, determination of the delamination mechanism and the different organization stages of anisotropic particles in dispersion is essential. We address this fundamental issue by exploiting a recently developed system of fluorohectorite smectite clay mineral that delaminates in water, leading to colloidal dispersions of single-layer, very large (≈20 µm) clay sheets at high dilution. We show that when the clay crystallites are dispersed in water, they swell to form periodic one-dimensional stacks of fluorohectorite sheets with very low volume fraction (<1%) and therefore huge (≈100 nm) periods. Using optical microscopy and synchrotron X-ray scattering, we establish that these colloidal stacks bear strong similarities, yet subtle differences, with a smectic liquid-crystalline phase. Despite the high dilution, the colloidal stacks of sheets, called colloidal accordions, are extremely robust mechanically and can persist for years. Moreover, when subjected to AC electric fields, they rotate as solid bodies, which demonstrates their outstanding internal cohesion. Furthermore, our theoretical model captures the dependence of the stacking period on the dispersion concentration and ionic strength and explains, invoking the Donnan effect, why the colloidal accordions are kinetically stable over years and impervious to shear and Brownian motion. Because our model is not system specific, we expect that similar colloidal accordions frequently appear as an intermediate state during the delamination process of two-dimensional crystals in polar solvents.

Innovative pressure environment combining hydrostatic pressure gradient and mechanical compression for structural investigations of nanoporous soft films.

The development of a new sample environment enabling X-ray scattering measurements at small and large angles under mechanical compression and hydraulic flow is presented. The cell, which is adapted for moderate pressures, includes beryllium windows, and allows applying simultaneously a compressive pressure up to 2.5 kbar in the perpendicular direction to the flow and either a hydrostatic pressure up to 300 bar or a pressure gradient of the same amplitude. The development of high-pressure devices for synchrotron experiments is relevant for many scientific fields in order to unveil details of a material's structure under relevant conditions of stresses. In particular, mechanical constraints coupled to hydrostatic pressure or flow, leading to complex stress tensor and mechanical response, and therefore unexpected deformations (swelling and pore deformation), are poorly addressed. Here, first the design of the environment is described, and then its performance with measurements carried out on a regenerated cellulose membrane is demonstrated.

Fine tuning the structural colours of photonic nanosheet suspensions by polymer doping.

Aqueous suspensions of nanosheets are readily obtained by exfoliating low-dimensional mineral compounds like H3Sb3P2O14. The nanosheets self-organize, at low concentration, into a periodic stack of membranes, i.e. a lamellar liquid-crystalline phase. Due to the dilution, this stack has a large period of a few hundred nanometres, it behaves as a 1-dimensional photonic material and displays structural colours. We experimentally investigated the dependence of the period on the nanosheet concentration. We theoretically showed that it cannot be explained by the usual DLVO interaction between uniform lamellae but that the particulate nature of nanosheet-laden membranes must be considered. Moreover, we observed that adding small amounts of 100 kDa poly(ethylene oxide) (PEO) decreases the period and allows tuning the colour throughout the visible range. PEO adsorbs on the nanosheets, inducing a strong reduction of the nanosheet charge. This is probably due to the Lewis-base character of the EO units of PEO that become protonated at the low pH of the system, an interpretation supported by theoretical modeling. Oddly enough, adding small amounts of 1 MDa PEO has the opposite effect of increasing the period, suggesting the presence of an additional intermembrane repulsion not yet identified. From an applied perspective, our work shows how the colours of these 1-dimensional photonic materials can easily be tuned not only by varying the nanosheet concentration (which might entail a phase transition) but also by adding PEO. From a theoretical perspective, our approach represents a necessary step towards establishing the phase diagram of aqueous suspensions of charged nanosheets.

Gold-clay nanocomposite colloids with liquid-crystalline and plasmonic properties.

Imparting liquid-crystal (LC) materials with the plasmonic properties of metal nanoparticles is actively pursued for applications. We achieved this goal by synthetizing gold nanoparticles onto clay nanosheets, leading to nematic nanocomposite suspensions. Optical observations and structural analysis show the growth of the gold nanoparticles without altering the LC properties of the nanosheets. These colloids display plasmonic structural colours and they can be aligned by an electric field, which is relevant for fundamental and materials chemistry of colloidal LC.

Organisation of clay nanoplatelets in a polyelectrolyte-based hydrogel.

We investigate the organisation of clay nanoplatelets within a hydrogel based on modified ionenes, cationic polyelectrolytes forming physically crosslinked hydrogels induced by hydrogen bonding and π-π stacking. Combination of small angle X-ray and neutron scattering (SAXS, SANS) reveals the structure of the polyelectrolyte network as well as the organisation of the clay additives. The clay-free hydrogel network features a characteristic mesh-size between 20 and 30 nm, depending on the polyelectrolyte concentration. Clay nanoplatelets inside the hydrogel organise in a regular face-to-face stacking manner, with a large repeat distance, following rather closely the hydrogel mesh-size. The presence of the nanoplatelets does not modify the hydrogel mesh size. Further, the clay-compensating counterions (Na+, Ca2+ or La3+) and the clay type (montmorillonite, beidellite) both have a significant influence on nanoplatelet organisation. The degree of nanoplatelet ordering in the hydrogel is very sensitive to the negative charge location on the clay platelet (different for each clay type). Increased nanoplatelet ordering leads to an improvement of the elastic properties of the hydrogel. On the contrary, the presence of dense clay aggregates (tactoids), induced by multi-valent clay counterions, destroys the hydrogel network as seen by the reduction of the elastic modulus of the hydrogel.

Green earth pigments dispersions: Water dynamics at the interfaces.

HYPOTHESIS: The objective is to elucidate the multiscale dynamics of water within natural mixtures of minerals, green earth pigments that are mainly composed of phyllosilicates containing large amount of iron. In particular, the interaction of water with the different kinds of surfaces has to be probed. One issue is to examine the influence of surface type, basal or edge, on the dispersion quality. EXPERIMENT: The study was carried out using 1H variable field NMR relaxometry on various green earth pigment dispersions and concentrations. To analyse the data, a new analytical model was developed for natural phyllosilicates containing large amount of paramagnetic centres. FINDING: The proposed theoretical framework is able to fit the experimental data for various samples using few parameters. It allows to determining water diffusion and residence times in complex phyllosilicate dispersions. Furthermore, it makes it possible to differentiate the contribution of the basal and edge surfaces and their respective surface area in interaction with water. Moreover, NMR relaxation profile reveals to be highly sensitive to the structural aspect of the phyllosilicates and to the accessibility of water to iron, hence allowing to discriminate clearly between two very similar phyllosilicates (glauconite and celadonite) that are difficult to distinguish by standard structural methods.

Destabilization of the Nematic Phase of Clay Nanosheet Suspensions by Polymer Adsorption.

Complex aqueous mixtures comprised of swelling clays and hydrosoluble polymers naturally occur in soils and play a major role in pedogenesis. They are also very often used for formulating oil-well drilling fluids, paints, and personal-care products. The suspensions of some natural clays, thanks to their large nanoparticle aspect ratio, spontaneously form nematic liquid-crystalline phases where the particles align parallel to each other, which affects their flow properties. We observed that adding small amounts of hydrosoluble polymers to these clay suspensions destabilizes the nematic phase with respect to the isotropic (disordered) phase. The polymers that we used (poly(ethylene oxide) and dextran) were too small to adopt particle-bridging conformations and small-angle X-ray scattering experiments showed that the structure of the nematic phase is not altered by polymer doping. However, the adsorption isotherm shows that the macromolecules adsorb onto the clay nanosheets, effectively coating them with a polymer layer. Our extension of Onsager's theory for polymer-coated platelets properly captures the experimental phase diagram and shows how the nematic phase destabilization can be due to the polymer adsorbing more on the platelet faces than at the rim. Because the flow properties of the nematic phase are very different from those of the isotropic phase, the presence or absence of the former phase is an important factor to be determined and considered to explain the rheological behavior of these complex systems.

Green Earth pigments aqueous dispersions: NMR relaxation rates dataset.

The data presented here are related to the research paper entitled "Green Earth pigments dispersions: water dynamics at the interfaces". The nuclear magnetic resonance (NMR) relaxometry data are provided for various aqueous Green Earth (GE) pigments dispersions with volume fraction spanning approximately from 0.1 to 0.5. For two of them (Cyprus GE and Bohemian GE), the NMR relaxation profiles from 10 kHz to 30 MHz (1H frequency) is given for several temperatures spanning from 293 to 318K. In addition, the X-ray diffraction pattern is provided for France GE (Kremer pigments) for the identification of the main mineral component. The nitrogen gas isotherms are provided for Cyprus GE and Bohemian GE.

Aggregation of Plate-like Colloids Induced by Charged Polymer Chains: Organization at the Nanometer Scale Tuned by Polymer Charge Density.

We study the aggregation of charged plate-like colloids, Na-montmorillonite clays, in the presence of ionenes, oppositely charged polymer chains. The choice of the charged polymer allows tuning its linear charge density to match/mismatch the average charge separation on the clay surfaces. We assess the nanoscale structure of the aggregates formed by small-angle X-ray and neutron scattering. The nanoscale features of the formed clay aggregates are dominated by the presence of a stacking peak, giving clear evidence for the formation of clay tactoids, that is, a face-to-face aggregation geometry of the clay platelets. The chain charge density of ionenes influences not only the stacking repeat distance within the clay tactoids but also the extent of stacking and abundance of the tactoids. We may distinguish two regimes as a function of clay and ionene polymer charge densities (ρc and ρp, respectively). The first regime applies to ρp > ρc and ρp ≈ ρc, that is, for highly and "matching" charged chains. Under these conditions, the intercalated chains lie in a flat conformation within the tactoids, irrespective of the ionic strength (within the range studied, i.e., up to 0.05 M NaBr). For weakly charged chains, ρp < ρc, undulation of the ionene chains within the tactoid is seen. The degree of undulation increases with ionic strength due to the decreasing persistence length of the ionene chains. The extent of stacking (5-10 platelets per tactoid) is a general feature of all the systems, and its origin remains unknown. The system corresponding to the closest match in charge separations on the clay surface and on the polymer chain (ρp ≈ ρc) features the highest abundance of tactoids. This coincides with the highest macroscopic density as deduced from simple visual inspection of sediment volumes. This leads to the open question regarding the link between the density at the nanoscale and the macroscopic density and sedimentation behavior of the aggregate.

One Step up the Ladder of Prebiotic Complexity: Formation of Nonrandom Linear Polypeptides from Binary Systems of Amino Acids on Silica.

Evidence for the formation of linear oligopeptides with nonrandom sequences from mixtures of amino acids coadsorbed on silica and submitted to a simple thermal activation is presented. The amino acid couples (glutamic acid+leucine) and (aspartic acid+valine) were deposited on a fumed silica and submitted to a single heating step at moderate temperature. The evolution of the systems was characterized by X-ray diffraction, infrared spectroscopy, thermosgravimetric analysis, HPLC, and electrospray ionization mass spectrometry (ESI-MS). Evidence for the formation of amide bonds was found in all systems studied. While the products of single amino acids activation on silica could be considered as evolutionary dead ends, (glutamic acid+leucine) and, at to some extent, (aspartic acid+valine) gave rise to the high yield formation of linear peptides up to the hexamers. Oligopeptides of such length have not been observed before in surface polymerization scenarios (unless the amino acids had been deposited by chemical vapor deposition, which is not realistic in a prebiotic environment). Furthermore, not all possible amino acid sequences were present in the activation products, which is indicative of polymerization selectivity. These results are promising for origins of life studies because they suggest the emergence of nonrandom biopolymers in a simple prebiotic scenario.

How do Nucleotides Adsorb Onto Clays?

Adsorption of prebiotic building blocks is proposed to have played a role in the emergence of life on Earth. The experimental and theoretical study of this phenomenon should be guided by our knowledge of the geochemistry of the habitable early Earth environments, which could have spanned a large range of settings. Adsorption being an interfacial phenomenon, experiments can be built around the minerals that probably exhibited the largest specific surface areas and were the most abundant, i.e., phyllosilicates. Our current work aims at understanding how nucleotides, the building blocks of RNA and DNA, might have interacted with phyllosilicates under various physico-chemical conditions. We carried out and refined batch adsorption studies to explore parameters such as temperature, pH, salinity, etc. We built a comprehensive, generalized model of the adsorption mechanisms of nucleotides onto phyllosilicate particles, mainly governed by phosphate reactivity. More recently, we used surface chemistry and geochemistry techniques, such as vibrational spectroscopy, low pressure gas adsorption, X-ray microscopy, and theoretical simulations, in order to acquire direct data on the adsorption configurations and localization of nucleotides on mineral surfaces. Although some of these techniques proved to be challenging, questioning our ability to easily detect biosignatures, they confirmed and complemented our pre-established model.

Soft X-ray Heterogeneous Radiolysis of Pyridine in the Presence of Hydrated Strontium-Hydroxyhectorite and its Monitoring by Near-Ambient Pressure Photoelectron Spectroscopy.

The heterogeneous radiolysis of organic molecules in clays is a matter of considerable interest in astrochemistry and environmental sciences. However, little is known about the effects of highly ionizing soft X-rays. By combining monochromatized synchrotron source irradiation with in situ Near Ambient Pressure X-ray Photoelectron Spectroscopy (in the mbar range), and using the synoptic view encompassing both the gas and condensed phases, we found the water and pyridine pressure conditions under which pyridine is decomposed in the presence of synthetic Sr2+-hydroxyhectorite. The formation of a pyridine/water/Sr2+ complex, detected from the Sr 3d and N 1s core-level binding energies, likely presents a favorable situation for the radiolytic breaking of the O-H bond of water molecules adsorbed in the clay and the subsequent decomposition of the molecule. However, decomposition stops when the pyridine pressure exceeds a critical value. This observation can be related to a change in the nature of the active radical species with the pyridine loading. This highlights the fact that the destruction of the molecule is not entirely determined by the properties of the host material, but also by the inserted organic species. The physical and chemical causes of the present observations are discussed.

Swelling Inhibition of Liquid Crystalline Colloidal Montmorillonite and Beidellite Clays by DNA.

Exploring the interaction of nucleic acids with clay minerals is important to understand such issues as the persistence in soils of biomolecules and the appearance of genetic polymers in prebiotic environments. Colloidal dispersions of double stranded DNA and clay nanosheets may also provide interesting model systems to study the statistical physics of mixtures of semi-flexible rods and plates. Here, we show that adding very small amounts of DNA to liquid-crystalline montmorillonite and beidellite smectite clay suspensions strongly widens the isotropic/nematic phase coexistence region. Moreover, a spectroscopic study shows that, upon DNA addition, the first DNA molecules adsorb onto the clay particles. Remarkably, synchrotron small-angle X-ray scattering experiments reveal that the average distance between the clay sheets, in the nematic phase at coexistence, decreases with increasing DNA concentration and that the inhibition of swelling by DNA becomes almost independent of clay concentration. We interpret this DNA-mediated attraction between clay nanosheets by bridging conformations of DNA strands (plates on a string structure). In addition to bridging, DNA chains can form "loops" between sections adsorbed on the same particle, giving rise to sheet repulsions due to protruding loops. This interpretation agrees with the observed inter-clay spacings being dependent only on the DNA concentration.

Effects of salinity on the adsorption of nucleotides onto phyllosilicates.

In the context of the origin of life, phyllosilicate surfaces might favor the adsorption, concentration and reactivity of otherwise diluted prebiotic molecules. The primitive oceanic seafloor was certainly rich in Fe-Mg-rich phyllosilicates. The salinity of the primitive seawater remains largely unknown. Values ranging from 1 to 15 times modern salinity have been proposed and the salt composition of the primitive ocean also remains elusive although it may have played a role in the interactions between nucleotides and mineral surfaces. Therefore we studied the adsorption of 5'-monophosphate deoxyguanosine (dGMP) as a model nucleotide onto a Fe-rich swelling clay, i.e. nontronite, and an Al-rich phyllosilicate, i.e. pyrophyllite, for comparison. Experiments were carried out at atmospheric pressure, 25 °C and natural pH, with a series of salts NaCl, MgCl2, CaCl2, MgSO4, NaH2PO4 and LaCl3 in order to evaluate the effect of cations and anions on dGMP adsorption. The present study shows that nucleotides are adsorbed on both phyllosilicates via a ligand exchange mechanism. The phosphate group of the nucleotide is adsorbed on the lateral metal hydroxyls of the broken edges of phyllosilicates. The presence of divalent cations or molecular anions, such as phosphate or sulfate, tends to inhibit this interaction on mineral surfaces. However, in the presence of divalent cations, cationic bridging on the basal surfaces of the swelling clay also occurs and could induce a higher retention capacity of the swelling clays compared to non-swelling phyllosilicates in primitive and modern natural environments.

Flocculation of Clay Colloids Induced by Model Polyelectrolytes: Effects of Relative Charge Density and Size.

Flocculation and its tuning are of utmost importance in the optimization of several industrial protocols in areas such as purification of waste water and civil engineering. Herein, we studied the polyelectrolyte-induced flocculation of clay colloids on a model system consisting of purified clay colloids of well-defined size fractions and ionene polyelectrolytes presenting regular and tunable chain charge density. To characterize ionene-induced clay flocculation, we turned to the combination of light absorbance (turbidity) and ζ-potential measurements, as well as adsorption isotherms. Our model system allowed us to identify the exact ratio of positive and negative charges in clay-ionene mixtures, the (c+/c-) ratio. For all samples studied, the onset of efficient flocculation occurred consistently at c+/c- ratios significantly below 1, which indicated the formation of highly ionene-deficient aggregates. At the same time, the ζ-potential measurements indicated an apparent zero charge on such aggregates. Thus, the ζ-potential values could not provide the stoichiometry inside the clay-ionene aggregates. The early onset of flocculation in clay-ionene mixtures is reminiscent of the behavior of multivalent salts and contrasts that of monovalent salts, for which a large excess amount of ions is necessary to achieve flocculation. Clear differences in the flocculation behavior are visible as a function of the ionene charge density, which governs the conformation of the ionene chains on the clay surface.

Modeling the arrangement of particles in natural swelling-clay porous media using three-dimensional packing of elliptic disks.

Swelling clay minerals play a key role in the control of water and pollutant migration in natural media such as soils. Moreover, swelling clay particles' orientational properties in porous media have significant implications for the directional dependence of fluid transfer. Herein we investigate the ability to mimic the organization of particles in natural swelling-clay porous media using a three-dimensional sequential particle deposition procedure [D. Coelho, J.-F. Thovert, and P. M. Adler, Phys. Rev. E 55, 1959 (1997)]. The algorithm considered is first used to simulate disk packings. Porosities of disk packings fall onto a single master curve when plotted against the orientational scalar order parameter value. This relation is used to validate the algorithm used in comparison with existing ones. The ellipticity degree of the particles is shown to have a negligible effect on the packing porosity for ratios ℓ(a)/ℓ(b) less than 1.5, whereas a significant increase in porosity is obtained for higher values. The effect of the distribution of the geometrical parameters (size, aspect ratio, and ellipticity degree) of particles on the final packing properties is also investigated. Finally, the algorithm is used to simulate particle packings for three size fractions of natural swelling-clay mineral powders. Calculated data regarding the distribution of the geometrical parameters and orientation of particles in porous media are successfully compared with experimental data obtained for the same samples. The results indicate that the obtained virtual porous media can be considered representative of natural samples and can be used to extract properties difficult to obtain experimentally, such as the anisotropic features of pore and solid phases in a system.