1
|
Auclair P. Mean-field approach to random Apollonian packing. Phys Rev E 2023; 107:034129. [PMID: 37073012 DOI: 10.1103/physreve.107.034129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2022] [Accepted: 03/02/2023] [Indexed: 04/20/2023]
Abstract
We revisit the scaling properties of growing spheres randomly seeded in d=2,3, and 4 dimensions using a mean-field approach. We model the insertion probability without assuming a priori a functional form for the radius distribution. The functional form of the insertion probability shows an unprecedented agreement with numerical simulations in d=2,3, and 4 dimensions. We infer from the insertion probability the scaling behavior of the random Apollonian packing and its fractal dimensions. The validity of our model is assessed with sets of 256 simulations each containing 20×10^{6} spheres in two, three, and four dimensions.
Collapse
Affiliation(s)
- Pierre Auclair
- Cosmology, Universe and Relativity at Louvain (CURL), Institute of Mathematics and Physics, University of Louvain, 2 Chemin du Cyclotron, 1348 Louvain-la-Neuve, Belgium
| |
Collapse
|
2
|
Ballow A, Linton P, Priour DJ. Percolation through voids around toroidal inclusions. Phys Rev E 2023; 107:014902. [PMID: 36797924 DOI: 10.1103/physreve.107.014902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2022] [Accepted: 01/04/2023] [Indexed: 01/28/2023]
Abstract
In the case of media comprised of impermeable particles, fluid flows through voids around impenetrable grains. For sufficiently low concentrations of the latter, spaces around grains join to allow transport on macroscopic scales, whereas greater impenetrable inclusion densities disrupt void networks and block macroscopic fluid flow. A critical grain concentration ρ_{c} marks the percolation transition or phase boundary separating these two regimes. With a dynamical infiltration technique in which virtual tracer particles explore void spaces, we calculate critical grain concentrations for randomly placed interpenetrating impermeable toroidal inclusions; the latter consist of surfaces of revolution with circular and square cross sections. In this manner, we study continuum percolation transitions involving nonconvex grains. As the radius of revolution increases relative to the length scale of the torus cross section, the tori develop a central hole, a topological transition accompanied by a cusp in the critical porosity fraction for percolation. With a further increase in the radius of revolution, as constituent grains become more ringlike in appearance, we find that the critical porosity fraction converges to that of high-aspect-ratio cylindrical counterparts only for randomly oriented grains.
Collapse
Affiliation(s)
- A Ballow
- Department of Physics and Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| | - P Linton
- Department of Physics and Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| | - D J Priour
- Department of Physics and Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| |
Collapse
|
3
|
He B, Ye A, Chi S, Mi P, Ran Y, Zhang L, Zou X, Pu B, Zhao Q, Zou Z, Wang D, Zhang W, Zhao J, Avdeev M, Shi S. CAVD, towards better characterization of void space for ionic transport analysis. Sci Data 2020; 7:153. [PMID: 32444597 PMCID: PMC7244509 DOI: 10.1038/s41597-020-0491-x] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/09/2020] [Accepted: 04/15/2020] [Indexed: 12/03/2022] Open
Abstract
Geometric crystal structure analysis using three-dimensional Voronoi tessellation provides intuitive insights into the ionic transport behavior of metal-ion electrode materials or solid electrolytes by mapping the void space in a framework onto a network. The existing tools typically consider only the local voids by mapping them with Voronoi polyhedra vertices and then define the mobile ions pathways using the Voronoi edges connecting these vertices. We show that in some structures mobile ions are located on Voronoi polyhedra faces and thus cannot be located by a standard approach. To address this deficiency, we extend the method to include Voronoi faces in the constructed network. This method has been implemented in the CAVD python package. Its effectiveness is demonstrated by 99% recovery rate for the lattice sites of mobile ions in 6,955 Li-, Na-, Mg- and Al-containing ionic compounds extracted from the Inorganic Crystal Structure Database. In addition, various quantitative descriptors of the network can be used to identify and rank the materials and further used in materials databases for machine learning.
Collapse
Affiliation(s)
- Bing He
- School of Computer Engineering and Science, Shanghai University, Shanghai, 200444, China
| | - Anjiang Ye
- School of Computer Engineering and Science, Shanghai University, Shanghai, 200444, China
| | - Shuting Chi
- School of Computer Engineering and Science, Shanghai University, Shanghai, 200444, China
| | - Penghui Mi
- School of Computer Engineering and Science, Shanghai University, Shanghai, 200444, China
| | - Yunbing Ran
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China
| | - Liwen Zhang
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China
| | - Xinxin Zou
- School of Computer Engineering and Science, Shanghai University, Shanghai, 200444, China
| | - Bowei Pu
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China
| | - Qian Zhao
- Materials Genome Institute, Shanghai University, Shanghai, 200444, China
| | - Zheyi Zou
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China
| | - Da Wang
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China
| | - Wenqing Zhang
- Department of Physics and Shenzhen Institute for Quantum Science & Technology, Southern University of Science and Technology, Shenzhen, 518055, China
| | - Jingtai Zhao
- School of Materials Science and Engineering, Guilin University of Electronic Technology, Guilin, 541004, China
| | - Maxim Avdeev
- Australian Nuclear Science and Technology Organisation, Locked Bag 2001, Kirrawee DC, NSW, 2232, Australia
- School of Chemistry, The University of Sydney, Sydney, 2006, Australia
| | - Siqi Shi
- State Key Laboratory of Advanced Special Steel, School of Materials Science and Engineering, Shanghai University, Shanghai, 200444, China.
- Materials Genome Institute, Shanghai University, Shanghai, 200444, China.
| |
Collapse
|
4
|
Ramirez LS, Centres PM, Ramirez-Pastor AJ. Percolation phase transition by removal of k^{2}-mers from fully occupied lattices. Phys Rev E 2019; 100:032105. [PMID: 31640014 DOI: 10.1103/physreve.100.032105] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2019] [Indexed: 11/07/2022]
Abstract
Numerical simulations and finite-size scaling analysis have been carried out to study the problem of inverse site percolation by the removal of k×k square tiles (k^{2}-mers) from square lattices. The process starts with an initial configuration, where all lattice sites are occupied and, obviously, the opposite sides of the lattice are connected by occupied sites. Then the system is diluted by removing k^{2}-mers of occupied sites from the lattice following a random sequential adsorption mechanism. The process finishes when the jamming state is reached and no more objects can be removed due to the absence of occupied sites clusters of appropriate size and shape. The central idea of this paper is based on finding the maximum concentration of occupied sites, p_{c,k}, for which the connectivity disappears. This particular value of the concentration is called the inverse percolation threshold and determines a well-defined geometrical phase transition in the system. The results obtained for p_{c,k} show that the inverse percolation threshold is a decreasing function of k in the range 1≤k≤4. For k≥5, all jammed configurations are percolating states, and consequently, there is no nonpercolating phase. In other words, the lattice remains connected even when the highest allowed concentration of removed sites is reached. The jamming exponent ν_{j} was measured, being ν_{j}=1 regardless of the size k considered. In addition, the accurate determination of the critical exponents ν, β, and γ reveals that the percolation phase transition involved in the system, which occurs for k varying between one and four, has the same universality class as the standard percolation problem.
Collapse
Affiliation(s)
- L S Ramirez
- Departamento de Física, Instituto de Física Aplicada, Universidad Nacional de San Luis-CONICET, Ejército de Los Andes 950, D5700HHW, San Luis, Argentina
| | - P M Centres
- Departamento de Física, Instituto de Física Aplicada, Universidad Nacional de San Luis-CONICET, Ejército de Los Andes 950, D5700HHW, San Luis, Argentina
| | - A J Ramirez-Pastor
- Departamento de Física, Instituto de Física Aplicada, Universidad Nacional de San Luis-CONICET, Ejército de Los Andes 950, D5700HHW, San Luis, Argentina
| |
Collapse
|
5
|
Gniewek P, Hallatschek O. Fluid flow through packings of elastic shells. Phys Rev E 2019; 99:023103. [PMID: 30934257 PMCID: PMC6542697 DOI: 10.1103/physreve.99.023103] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/13/2018] [Indexed: 11/07/2022]
Abstract
Fluid transport in porous materials is commonly studied in geological samples (soil, sediments, etc.) or idealized systems, but the fluid flow through compacted granular materials, consisting of substantially strained granules, remains relatively unexplored. As a step toward filling this gap, we study a model of liquid transport in packings of deformable elastic shells using finite-element and lattice-Boltzmann methods. We find that the fluid flow abruptly vanishes as the porosity of the material falls below a critical value, and the flow obstruction exhibits features of a percolation transition. We further show that the fluid flow can be captured by a simplified permeability model in which the complex porous material is replaced by a collection of disordered capillaries, which are distributed and shaped by the percolation transition. To that end, we numerically explore the divergence of hydraulic tortuosity τ_{H} and the decrease of a hydraulic radius R_{h} as the percolation threshold is approached. We interpret our results in terms of scaling predictions derived from the percolation theory applied to random packings of spheres.
Collapse
Affiliation(s)
- Pawel Gniewek
- Biophysics Graduate Group, University of California, Berkeley, USA
| | - Oskar Hallatschek
- Departments of Physics and Integrative Biology, University of California, Berkeley, USA
| |
Collapse
|
6
|
Priour DJ, McGuigan NJ. Percolation through Voids around Randomly Oriented Polyhedra and Axially Symmetric Grains. PHYSICAL REVIEW LETTERS 2018; 121:225701. [PMID: 30547614 DOI: 10.1103/physrevlett.121.225701] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/30/2018] [Revised: 07/10/2018] [Indexed: 06/09/2023]
Abstract
Porous materials made up of impermeable grains constrain fluid flow to voids around the impenetrable inclusions. A percolation transition marks the boundary between densities of grains permitting bulk transport and concentrations blocking traversal on macroscopic scales. With dynamical infiltration of void spaces using virtual tracer particles, we treat inclusion geometries exactly. We calculate the critical number density per volume ρ_{c} for a variety of axially symmetric shapes and faceted solids with the former including cylinders, ellipsoids, cones, and tablet shaped grains from highly oblate (platelike) to highly prolate (needlelike) extremes. For the latter, results suggest a common asymptotic value identical to the counterpart for aligned cylindrical grains. We find percolation thresholds for each of the five platonic solids (i.e., tetrahedra, cubes, octahedra, dodecahedra, and icosahedra) as well as truncated icosahedra. For each polyhedron type, we consider aligned and randomly oriented grains, finding distinct percolation thresholds for the former versus the latter only for cubes. The anomalous diffusion exponents we find differ from those of the universality class for discrete models on 3D lattices.
Collapse
Affiliation(s)
- D J Priour
- Department of Physics & Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| | - N J McGuigan
- Department of Physics & Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| |
Collapse
|
7
|
Waterflooding of Surfactant and Polymer Solutions in a Porous Media Micromodel. COLLOIDS AND INTERFACES 2018. [DOI: 10.3390/colloids2020023] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
8
|
Three-Dimensional Lattice Boltzmann Simulations of Single-Phase Permeability in Random Fractal Porous Media with Rough Pore–Solid Interface. Transp Porous Media 2018. [DOI: 10.1007/s11242-017-0938-5] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/17/2022]
|
9
|
Ongari D, Boyd PG, Barthel S, Witman M, Haranczyk M, Smit B. Accurate Characterization of the Pore Volume in Microporous Crystalline Materials. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2017; 33:14529-14538. [PMID: 28636815 PMCID: PMC5745516 DOI: 10.1021/acs.langmuir.7b01682] [Citation(s) in RCA: 76] [Impact Index Per Article: 10.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/19/2017] [Revised: 06/20/2017] [Indexed: 05/19/2023]
Abstract
Pore volume is one of the main properties for the characterization of microporous crystals. It is experimentally measurable, and it can also be obtained from the refined unit cell by a number of computational techniques. In this work, we assess the accuracy and the discrepancies between the different computational methods which are commonly used for this purpose, i.e, geometric, helium, and probe center pore volumes, by studying a database of more than 5000 frameworks. We developed a new technique to fully characterize the internal void of a microporous material and to compute the probe-accessible and -occupiable pore volume. We show that, unlike the other definitions of pore volume, the occupiable pore volume can be directly related to the experimentally measured pore volumes from nitrogen isotherms.
Collapse
Affiliation(s)
- Daniele Ongari
- Laboratory
of Molecular Simulation, Institut des Sciences et Ingeénierie
Chimiques, Ecole Polytechnique Fédérale
de Lausanne (EPFL), Rue de l’Industrie 17, CH-1951 Sion, Valais, Switzerland
| | - Peter G. Boyd
- Laboratory
of Molecular Simulation, Institut des Sciences et Ingeénierie
Chimiques, Ecole Polytechnique Fédérale
de Lausanne (EPFL), Rue de l’Industrie 17, CH-1951 Sion, Valais, Switzerland
| | - Senja Barthel
- Laboratory
of Molecular Simulation, Institut des Sciences et Ingeénierie
Chimiques, Ecole Polytechnique Fédérale
de Lausanne (EPFL), Rue de l’Industrie 17, CH-1951 Sion, Valais, Switzerland
| | - Matthew Witman
- Department
of Chemical and Biomolecular Engineering, University of California, Berkeley, California 94760, United States
| | - Maciej Haranczyk
- IMDEA
Materials Institute, C/Eric Kandel 2, 28906 Getafe, Madrid, Spain
- Computational
Research Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
| | - Berend Smit
- Laboratory
of Molecular Simulation, Institut des Sciences et Ingeénierie
Chimiques, Ecole Polytechnique Fédérale
de Lausanne (EPFL), Rue de l’Industrie 17, CH-1951 Sion, Valais, Switzerland
- E-mail:
| |
Collapse
|
10
|
Modeling Oil Recovery for Mixed Macro- and Micro-Pore Carbonate Grainstones. Sci Rep 2017; 7:9780. [PMID: 28851934 PMCID: PMC5575129 DOI: 10.1038/s41598-017-09507-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2017] [Accepted: 07/27/2017] [Indexed: 11/08/2022] Open
Abstract
In general, modeling oil-recovery is a challenging problem involving detailed fluid flow calculations with required structural details that challenge current experimental resolution. Recent laboratory experiments on mixed micro- and macro-pore suggest that there is a systematic relationship between remaining oil saturation (ROS) and the fraction of micro-pores. Working with experimental measurements of the pores obtained from X-ray tomography and mercury intrusion capillary pressure porosimetry, we define a digital rock model exemplifying the key structural elements of these carbonate grainstones. We then test two fluid-flow models: invasion percolation model and effective medium model. Although invasion percolation identifies the important impact of macro-pore percolation on permeability, it does not describe the dependence of ROS on micro-pore percentage. We thus modified the effective medium model by introducing a single-parameter descriptor, reff. Oil from pores r ≥ reff is fully removed, while for the remaining pores with r < reff, their contribution is scaled by (r/reff)2. Applying this straightforward physics to pore size distributions for the mixed-pore grainstones reproduces the experimental ROS dependence.
Collapse
|
11
|
|
12
|
Jin Y, Charbonneau P. Dimensional study of the dynamical arrest in a random Lorentz gas. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:042313. [PMID: 25974497 DOI: 10.1103/physreve.91.042313] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/01/2014] [Indexed: 06/04/2023]
Abstract
The random Lorentz gas (RLG) is a minimal model for transport in heterogeneous media. Upon increasing the obstacle density, it exhibits a growing subdiffusive transport regime and then a dynamical arrest. Here, we study the dimensional dependence of the dynamical arrest, which can be mapped onto the void percolation transition for Poisson-distributed point obstacles. We numerically determine the arrest in dimensions d=2-6. Comparison of the results with standard mode-coupling theory reveals that the dynamical theory prediction grows increasingly worse with d. In an effort to clarify the origin of this discrepancy, we relate the dynamical arrest in the RLG to the dynamic glass transition of the infinite-range Mari-Kurchan-model glass former. Through a mixed static and dynamical analysis, we then extract an improved dimensional scaling form as well as a geometrical upper bound for the arrest. The results suggest that understanding the asymptotic behavior of the random Lorentz gas may be key to surmounting fundamental difficulties with the mode-coupling theory of glasses.
Collapse
Affiliation(s)
- Yuliang Jin
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
- Dipartimento di Fisica, Sapienza Università di Roma and INFN, Sezione di Roma I, IPFC-CNR, Piazzale Aldo Moro 2, I-00185 Roma, Italy
- LPT, École Normale Supérieure, UMR 8549 CNRS, 24 Rue Lhomond, 75005 Paris, France
| | - Patrick Charbonneau
- Department of Chemistry, Duke University, Durham, North Carolina 27708, USA
- Department of Physics, Duke University, Durham, North Carolina 27708, USA
| |
Collapse
|
13
|
Charbonneau P, Jin Y, Parisi G, Zamponi F. Hopping and the Stokes-Einstein relation breakdown in simple glass formers. Proc Natl Acad Sci U S A 2014; 111:15025-30. [PMID: 25288722 PMCID: PMC4210276 DOI: 10.1073/pnas.1417182111] [Citation(s) in RCA: 65] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
One of the most actively debated issues in the study of the glass transition is whether a mean-field description is a reasonable starting point for understanding experimental glass formers. Although the mean-field theory of the glass transition--like that of other statistical systems--is exact when the spatial dimension d → ∞, the evolution of systems properties with d may not be smooth. Finite-dimensional effects could dramatically change what happens in physical dimensions,d = 2, 3. For standard phase transitions finite-dimensional effects are typically captured by renormalization group methods, but for glasses the corrections are much more subtle and only partially understood. Here, we investigate hopping between localized cages formed by neighboring particles in a model that allows to cleanly isolate that effect. By bringing together results from replica theory, cavity reconstruction, void percolation, and molecular dynamics, we obtain insights into how hopping induces a breakdown of the Stokes-Einstein relation and modifies the mean-field scenario in experimental systems. Although hopping is found to supersede the dynamical glass transition, it nonetheless leaves a sizable part of the critical regime untouched. By providing a constructive framework for identifying and quantifying the role of hopping, we thus take an important step toward describing dynamic facilitation in the framework of the mean-field theory of glasses.
Collapse
Affiliation(s)
| | - Yuliang Jin
- Departments of Chemistry and Dipartimento di Fisica, Sapienza Università di Roma, I-00185 Rome, Italy;
| | - Giorgio Parisi
- Dipartimento di Fisica, Sapienza Università di Roma, I-00185 Rome, Italy; Istituto Nazionale di Fisica Nucleare, Sezione di Roma I, Istituto per i Processi Chimico Fisici, Consiglio Nazionale delle Ricerche, I-00185 Rome, Italy; and
| | - Francesco Zamponi
- Laboratoire de Physique Théorique, Ecole Normale Supérieure, UMR 8549 CNRS, 75005 Paris, France
| |
Collapse
|
14
|
Hebert V, Garing C, Luquot L, Pezard PA, Gouze P. Multi-scale X-ray tomography analysis of carbonate porosity. ACTA ACUST UNITED AC 2014. [DOI: 10.1144/sp406.12] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2022]
Abstract
AbstractThe porosity (from values of <1% up to more than 50%, with a mean value of 36.5%) and permeability (from very low values to 20.4 D, with a mean value around 800 mD) of 38 plugs, with a diameter 40 mm, have been measured. These plugs were taken approximately every 1.5 m along a 100 m-long core, sampling a series of reefal carbonate platform lithofacies. A non-unique relationship between porosity and permeability is obtained. In addition, the porosity was measured from processing X-ray micro-tomography images with resolutions ranging from 0.42 to 190.0 µm and sampling volumes ranging from less than 1 µm3 to few cm3. Depending on the structure of the pore network, the computed porosity is, in some cases, controlled by the X-ray image resolution and the sampled volume. While high resolution is required to image micro-porous material, large samples with lower resolution images are necessary for identifying distributed vugs and millimetre-scale structures. This study shows that multi-resolution X-ray micro-tomography is an effective tool for characterizing the multi-scale pore structure of carbonate rocks and understanding how it may control key petrophysical parameters such as porosity and permeability.
Collapse
Affiliation(s)
- V. Hebert
- Géosciences Montpellier, UMR 5243 CNRS/INSU, Université de Montpellier 2, 34095 Montpellier, France
| | - C. Garing
- Géosciences Montpellier, UMR 5243 CNRS/INSU, Université de Montpellier 2, 34095 Montpellier, France
| | - L. Luquot
- Géosciences Montpellier, UMR 5243 CNRS/INSU, Université de Montpellier 2, 34095 Montpellier, France
| | - P. A. Pezard
- Géosciences Montpellier, UMR 5243 CNRS/INSU, Université de Montpellier 2, 34095 Montpellier, France
| | - P. Gouze
- Géosciences Montpellier, UMR 5243 CNRS/INSU, Université de Montpellier 2, 34095 Montpellier, France
| |
Collapse
|
15
|
Priour DJ. Percolation through voids around overlapping spheres: a dynamically based finite-size scaling analysis. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:012148. [PMID: 24580213 DOI: 10.1103/physreve.89.012148] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/01/2012] [Indexed: 06/03/2023]
Abstract
The percolation threshold for flow or conduction through voids surrounding randomly placed spheres is calculated. With large-scale Monte Carlo simulations, we give a rigorous continuum treatment to the geometry of the impenetrable spheres and the spaces between them. To properly exploit finite-size scaling, we examine multiple systems of differing sizes, with suitable averaging over disorder, and extrapolate to the thermodynamic limit. An order parameter based on the statistical sampling of stochastically driven dynamical excursions and amenable to finite-size scaling analysis is defined, calculated for various system sizes, and used to determine the critical volume fraction ϕc=0.0317±0.0004 and the correlation length exponent ν=0.92±0.05.
Collapse
Affiliation(s)
- D J Priour
- Department of Science, Kansas City Kansas Community College, Kansas City, Kansas 66112, USA and Department of Physics and Astronomy, Youngstown State University, Youngstown, Ohio 44555, USA
| |
Collapse
|
16
|
Characterization and comparison of pore landscapes in crystalline porous materials. J Mol Graph Model 2013; 44:208-19. [PMID: 23876827 DOI: 10.1016/j.jmgm.2013.05.007] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2013] [Revised: 05/03/2013] [Accepted: 05/11/2013] [Indexed: 11/23/2022]
Abstract
Crystalline porous materials have many applications, including catalysis and separations. Identifying suitable materials for a given application can be achieved by screening material databases. Such a screening requires automated high-throughput analysis tools that characterize and represent pore landscapes with descriptors, which can be compared using similarity measures in order to select, group and classify materials. Here, we discuss algorithms for the calculation of two types of pore landscape descriptors: pore size distributions and stochastic rays. These descriptors provide histogram representations that encode the geometrical properties of pore landscapes. Their calculation involves the Voronoi decomposition as a technique to map and characterize accessible void space inside porous materials. Moreover, we demonstrate pore landscape comparisons for materials from the International Zeolite Association (IZA) database of zeolite frameworks, and illustrate how the choice of pore descriptor and similarity measure affects the perspective of material similarity exhibiting a particular emphasis and sensitivity to certain aspects of structures.
Collapse
|
17
|
Höfling F, Franosch T. Anomalous transport in the crowded world of biological cells. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2013; 76:046602. [PMID: 23481518 DOI: 10.1088/0034-4885/76/4/046602] [Citation(s) in RCA: 584] [Impact Index Per Article: 53.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
A ubiquitous observation in cell biology is that the diffusive motion of macromolecules and organelles is anomalous, and a description simply based on the conventional diffusion equation with diffusion constants measured in dilute solution fails. This is commonly attributed to macromolecular crowding in the interior of cells and in cellular membranes, summarizing their densely packed and heterogeneous structures. The most familiar phenomenon is a sublinear, power-law increase of the mean-square displacement (MSD) as a function of the lag time, but there are other manifestations like strongly reduced and time-dependent diffusion coefficients, persistent correlations in time, non-Gaussian distributions of spatial displacements, heterogeneous diffusion and a fraction of immobile particles. After a general introduction to the statistical description of slow, anomalous transport, we summarize some widely used theoretical models: Gaussian models like fractional Brownian motion and Langevin equations for visco-elastic media, the continuous-time random walk model, and the Lorentz model describing obstructed transport in a heterogeneous environment. Particular emphasis is put on the spatio-temporal properties of the transport in terms of two-point correlation functions, dynamic scaling behaviour, and how the models are distinguished by their propagators even if the MSDs are identical. Then, we review the theory underlying commonly applied experimental techniques in the presence of anomalous transport like single-particle tracking, fluorescence correlation spectroscopy (FCS) and fluorescence recovery after photobleaching (FRAP). We report on the large body of recent experimental evidence for anomalous transport in crowded biological media: in cyto- and nucleoplasm as well as in cellular membranes, complemented by in vitro experiments where a variety of model systems mimic physiological crowding conditions. Finally, computer simulations are discussed which play an important role in testing the theoretical models and corroborating the experimental findings. The review is completed by a synthesis of the theoretical and experimental progress identifying open questions for future investigation.
Collapse
Affiliation(s)
- Felix Höfling
- Max-Planck-Institut für Intelligente Systeme, Heisenbergstraße 3, 70569 Stuttgart, and Institut für Theoretische Physik IV, Universität Stuttgart, Pfaffenwaldring 57, 70569 Stuttgart, Germany
| | | |
Collapse
|
18
|
Pinheiro M, Martin RL, Rycroft CH, Haranczyk M. High accuracy geometric analysis of crystalline porous materials. CrystEngComm 2013. [DOI: 10.1039/c3ce41057a] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
19
|
Maiti M, Lakshminarayanan A, Sastry S. Characterization of void space in polydisperse sphere packings: Applications to hard-sphere packings and to protein structure analysis. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2013; 36:5. [PMID: 23355091 DOI: 10.1140/epje/i2013-13005-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/20/2012] [Revised: 10/27/2012] [Accepted: 12/19/2012] [Indexed: 06/01/2023]
Abstract
The implementation of a method for the exact evaluation of the volume and surface area of cavities and free volumes in polydisperse sphere packings is described. The generalization of an algorithm for Voronoi tessellation by Tanemura et al. is presented, employing the radical plane construction, as a part of the method. We employ this method to calculate the equation of state for monodisperse and polydisperse hard-sphere fluids, crystals, and for the metastable amorphous branch up to random close packing or jamming densities. We compute the distribution of free volumes, and compare with previous results employing a heuristic definition of free volume. We show the efficacy of the method for analyzing protein structure, by computing various quantities such as the distribution of sizes of buried cavities and pockets, the scaling of solvent accessible area to the corresponding occupied volume, the composition of residues lining cavities, etc.
Collapse
Affiliation(s)
- Moumita Maiti
- Theoretical Sciences Unit, JNCASR, 560065, Jakkur, Banglore, India
| | | | | |
Collapse
|
20
|
Cho HW, Kwon G, Sung BJ, Yethiraj A. Effect of polydispersity on diffusion in random obstacle matrices. PHYSICAL REVIEW LETTERS 2012; 109:155901. [PMID: 23102336 DOI: 10.1103/physrevlett.109.155901] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/01/2012] [Revised: 08/13/2012] [Indexed: 06/01/2023]
Abstract
The dynamics of tracers in disordered matrices is of interest in a number of diverse areas of physics such as the biophysics of crowding in cells and cell membranes, and the diffusion of fluids in porous media. To a good approximation the matrices can be modeled as a collection of spatially frozen particles. In this Letter, we consider the effect of polydispersity (in size) of the matrix particles on the dynamics of tracers. We study a two dimensional system of hard disks diffusing in a sea of hard disk obstacles, for different values of the polydispersity of the matrix. We find that for a given average size and area fraction, the diffusion of tracers is very sensitive to the polydispersity. We calculate the pore percolation threshold using Apollonius diagrams. The diffusion constant, D, follows a scaling relation D~(φ(c)-φ(m))(μ-β) for all values of the polydispersity, where φ(m) is the area fraction and φ(c) is the value of φ(m) at the percolation threshold.
Collapse
Affiliation(s)
- Hyun Woo Cho
- Department of Chemistry, Sogang University, Seoul, Republic of Korea
| | | | | | | |
Collapse
|
21
|
Novak IL, Gao F, Kraikivski P, Slepchenko BM. Diffusion amid random overlapping obstacles: similarities, invariants, approximations. J Chem Phys 2011; 134:154104. [PMID: 21513372 PMCID: PMC3094463 DOI: 10.1063/1.3578684] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2010] [Accepted: 03/26/2011] [Indexed: 11/14/2022] Open
Abstract
Efficient and accurate numerical techniques are used to examine similarities of effective diffusion in a void between random overlapping obstacles: essential invariance of effective diffusion coefficients (D(eff)) with respect to obstacle shapes and applicability of a two-parameter power law over nearly entire range of excluded volume fractions (φ), except for a small vicinity of a percolation threshold. It is shown that while neither of the properties is exact, deviations from them are remarkably small. This allows for quick estimation of void percolation thresholds and approximate reconstruction of D(eff) (φ) for obstacles of any given shape. In 3D, the similarities of effective diffusion yield a simple multiplication "rule" that provides a fast means of estimating D(eff) for a mixture of overlapping obstacles of different shapes with comparable sizes.
Collapse
Affiliation(s)
- Igor L Novak
- Richard D. Berlin Center for Cell Analysis and Modeling, Department of Cell Biology, University of Connecticut Health Center, Farmington, Connecticut 06030, USA.
| | | | | | | |
Collapse
|
22
|
Höfling F, Munk T, Frey E, Franosch T. Critical dynamics of ballistic and Brownian particles in a heterogeneous environment. J Chem Phys 2008; 128:164517. [PMID: 18447469 DOI: 10.1063/1.2901170] [Citation(s) in RCA: 46] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The dynamic properties of a classical tracer particle in a random, disordered medium are investigated close to the localization transition. For Lorentz models obeying Newtonian and diffusive motion at the microscale, we have performed large-scale computer simulations, demonstrating that universality holds at long times in the immediate vicinity of the transition. The scaling function describing the crossover from anomalous transport to diffusive motion is found to vary extremely slowly and spans at least five decades in time. To extract the scaling function, one has to allow for the leading universal corrections to scaling. Our findings suggest that apparent power laws with varying exponents generically occur and dominate experimentally accessible time windows as soon as the heterogeneities cover a decade in length scale. We extract the divergent length scales, quantify the spatial heterogeneities in terms of the non-Gaussian parameter, and corroborate our results by a thorough finite-size analysis.
Collapse
Affiliation(s)
- Felix Höfling
- Arnold Sommerfeld Center for Theoretical Physics (ASC) and Center for NanoScience, Fakultät für Physik, Ludwig-Maximilians-Universität München, Theresienstrasse 37, 80333 München, Germany.
| | | | | | | |
Collapse
|
23
|
Sung BJ, Yethiraj A. The effect of matrix structure on the diffusion of fluids in porous media. J Chem Phys 2008; 128:054702. [PMID: 18266457 DOI: 10.1063/1.2823735] [Citation(s) in RCA: 38] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The effect of matrix structure on the transport properties of adsorbed fluids is studied using computer simulations and percolation theory. The model system consists of a fluid of hard spheres diffusing in a matrix of hard spheres fixed in space. Three different arrangements of the fixed spheres, random, templated, and polymeric, are investigated. For a given matrix volume fraction the diffusion coefficient of the fluid, D, is sensitive to the manner in which the matrix is constructed, with large differences between the three types of matrices. The matrix is mapped onto an effective lattice composed of vertices and bonds using a Voronoi tessellation method where the connectivity of bonds is determined using a geometric criterion, i.e., a bond is connected if a fluid particle can pass directly between the two pores the bond connects, and disconnected otherwise. The percolation threshold is then determined from the connectivity of the bonds. D displays universal scaling behavior in the reduced volume fraction, i.e., D approximately (1-phi(m)phi(c))(gamma), where phi(m) is the matrix volume fraction and phi(c) is the matrix volume fraction at the percolation threshold. We find that gamma approximately 2.2, independent of matrix type, which is different from the result gamma approximately 1.53 for diffusion in lattice models, but similar to that for conduction in Swiss cheese models. Lattice simulations with biased hopping probabilities are consistent with the continuous-space simulations, and this shows that the universal behavior of diffusion is sensitive to details of local dynamics.
Collapse
Affiliation(s)
- Bong June Sung
- Department of Chemistry, Sogang University, Seoul 121-742, Republic of Korea
| | | |
Collapse
|
24
|
Yi YB. Void percolation and conduction of overlapping ellipsoids. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 74:031112. [PMID: 17025599 DOI: 10.1103/physreve.74.031112] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/07/2006] [Revised: 07/17/2006] [Indexed: 05/12/2023]
Abstract
The void percolation and conduction problems for equisized overlapping ellipsoids of revolution are investigated using the discretization method. The method is validated by comparing the estimated percolation threshold of spheres with the precise result found in literature. The technique is then extended to determine the threshold of void percolation as a function of the geometric aspect ratio of ellipsoidal particles. The finite element method is also applied to evaluate the equivalent conductivity of the void phase in the system. The results confirm that there are no universalities for void percolation threshold and conductivity in particulate systems, and these properties are clearly dependent on the geometrical shape of particles. As a consequence, void percolation and conduction associated with ellipsoidal particles of large aspect ratio should be treated differently from spheres.
Collapse
Affiliation(s)
- Y B Yi
- Department of Engineering, University of Denver, Denver, Colorodo 80208, USA
| |
Collapse
|
25
|
Sung BJ, Yethiraj A. Lateral diffusion and percolation in membranes. PHYSICAL REVIEW LETTERS 2006; 96:228103. [PMID: 16803348 DOI: 10.1103/physrevlett.96.228103] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/16/2006] [Indexed: 05/10/2023]
Abstract
An algorithm based on Voronoi tessellation and percolation theory is presented to study the diffusion of model membrane components (solutes) in the plasma membrane. The membrane is modeled as a two-dimensional space with integral membrane proteins as static obstacles. The Voronoi diagram consists of vertices, which are equidistant from three matrix obstacles, joined by edges. An edge between two vertices is said to be connected if solute particles can pass directly between the two regions. The percolation threshold, pc, determined using this passage criterion is pc approximately equal to 0.53. This is smaller than if the connectivity of edges were assigned randomly, in which case the percolation threshold pr=2/3, where p is the fraction of connected edges. Molecular dynamics simulations show that diffusion is determined by percolation of clusters of edges.
Collapse
Affiliation(s)
- Bong June Sung
- Theoretical Chemistry Institute and Department of Chemistry, University of Wisconsin, Madison, Wisconsin 53706, USA
| | | |
Collapse
|
26
|
Abstract
We consider growing spheres seeded by random injection in time and space. Growth stops when two spheres meet leading eventually to a jammed state. We study the statistics of growth limited by packing theoretically in d dimensions and via simulation in d=2, 3, and 4. We show how a broad class of such models exhibit distributions of sphere radii with a universal exponent. We construct a scaling theory that relates the fractal structure of these models to the decay of their pore space, a theory that we confirm via numerical simulations. The scaling theory also predicts an upper bound for the universal exponent and is in exact agreement with numerical results for d=4.
Collapse
|
27
|
Abstract
The average packing density inside proteins is as high as in crystalline solids. Does this mean proteins are well-packed? We go beyond average densities, and look at the full distribution functions of free volumes inside proteins. Using a new and rigorous Delaunay triangulation method for parsing space into empty and filled regions, we introduce formal definitions of interior and surface packing densities. Although proteins look like organic crystals by the criterion of average density, they look more like liquids and glasses by the criterion of their free volume distributions. The distributions are broad, and the scalings of volume-to-surface, volume-to-cluster-radius, and numbers of void versus volume show that the interiors of proteins are more like randomly packed spheres near their percolation threshold than like jigsaw puzzles. We find that larger proteins are packed more loosely than smaller proteins. And we find that the enthalpies of folding (per amino acid) are independent of the packing density of a protein, indicating that van der Waals interactions are not a dominant component of the folding forces.
Collapse
Affiliation(s)
- J Liang
- Department of Bioengineering, University of Illinois at Chicago, Chicago, Illinois 60607-7052, USA.
| | | |
Collapse
|
28
|
Vocka R, Dubois MA. Pore network as a model of porous media: comparison between nonhierarchical and hierarchical organizations of pores. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 2000; 62:5216-5224. [PMID: 11089083 DOI: 10.1103/physreve.62.5216] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/16/1998] [Revised: 06/02/2000] [Indexed: 05/23/2023]
Abstract
A pore network is used as a model of a porous medium with hierarchical and single level organization of the pores. A microscopic diffusion or adsorption process is introduced that permits the study of both temporal evolution and effective properties of transport. The relationships between the model parameters and the measurable quantities are derived. The transport in media with a hierarchical organization of pores is shown to be qualitatively different from that in media lacking this organization. The question of experimental distinction between the two types of pore organization is also studied. We show that mercury intrusion porosimetry (MIP) measurements furnish valuable information about the pore organization in a sample, but that the exact deconvolution of the real pore size distribution from the experimental data is not straightforward. Our work provides indications for a correct interpretation of MIP results. Qualitative comparison with experiments show the pertinence of hierarchical models.
Collapse
Affiliation(s)
- R Vocka
- Service de Physique de l'Etat Condense, SPEC/DRECAM, CEN Saclay, F-91191 Gif-sur-Yvette, France
| | | |
Collapse
|
29
|
Rintoul MD. Precise determination of the void percolation threshold for two distributions of overlapping spheres. PHYSICAL REVIEW. E, STATISTICAL PHYSICS, PLASMAS, FLUIDS, AND RELATED INTERDISCIPLINARY TOPICS 2000; 62:68-72. [PMID: 11088435 DOI: 10.1103/physreve.62.68] [Citation(s) in RCA: 55] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2000] [Indexed: 11/06/2022]
Abstract
The void percolation threshold is calculated for a distribution of overlapping spheres with equal radii, and for a binary-sized distribution of overlapping spheres, where half of the spheres have radii twice as large as the other half. Using systems much larger than previous work, we determine a much more precise value for the percolation thresholds and correlation length exponent. The value of the percolation threshold for the monodisperse case is shown to be 0. 0301+/-0.0003, whereas the value for the bidisperse case is shown to be p(c)=0.0287+/-0.0005. The fact that these are significantly different is in contrast with previous, less precise works that speculated that the threshold might be universal with respect to sphere size distribution.
Collapse
Affiliation(s)
- MD Rintoul
- MS 1111, Sandia National Laboratories, P.O. Box 5800, Albuquerque, New Mexico 87185-1111, USA
| |
Collapse
|
30
|
Sahimi M, Tsotsis TT. Transient Diffusion and Conduction in Heterogeneous Media: Beyond the Classical Effective-Medium Approximation. Ind Eng Chem Res 1997. [DOI: 10.1021/ie960602k] [Citation(s) in RCA: 33] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Muhammad Sahimi
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089-1211
| | - Theodore T. Tsotsis
- Department of Chemical Engineering, University of Southern California, Los Angeles, California 90089-1211
| |
Collapse
|