1
|
Sun Y, Wang C, Yang J, Shi W, Pang Q, Wang Y, Li J, Hu B, Xia C. Evident structural anisotropies arising from near-zero particle asphericity in granular spherocylinder packings. Phys Rev E 2024; 110:014903. [PMID: 39161035 DOI: 10.1103/physreve.110.014903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2024] [Accepted: 07/02/2024] [Indexed: 08/21/2024]
Abstract
With magnetic resonance imaging experiments, we study packings of granular spherocylinders with merely 2% asphericity. Evident structural anisotropies across all length scales are identified. Most interestingly, the global nematic order decreases with increasing packing fraction, while the local contact anisotropy shows an opposing trend. We attribute this counterintuitive phenomenon to a competition between gravity-driven ordering aided by frictional contacts and a geometric frustration effect at the marginally jammed state. It is also surprising to notice that such slight particle asphericity can trigger non-negligible correlations between contact-level and mesoscale structures, manifested in drastically different nonaffine structural rearrangements upon compaction from that of granular spheres. These observations can help improve statistical mechanical models for the orientational order transformation of nonspherical granular particle packings, which involves complex interplays between particle shape, frictional contacts, and external force field.
Collapse
Affiliation(s)
| | | | | | | | | | - Yujie Wang
- School of Physics and Astronomy, Shanghai Jiao Tong University, Shanghai 200240, China
- State Key Laboratory of Geohazard Prevention and Geoenvironment Protection, Chengdu University of Technology, Chengdu 610059, China
- Department of Physics, College of Mathematics and Physics, Chengdu University of Technology, Chengdu 610059, China
| | | | | | | |
Collapse
|
2
|
Wang W, Chen Z, Gao Y, Chen C, Jiao Y, Zhang S. Spheroid models to elaborate the broken symmetry and equivalent volume of molecules in crystalline phase. Phys Rev E 2024; 109:064603. [PMID: 39020901 DOI: 10.1103/physreve.109.064603] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2023] [Accepted: 05/19/2024] [Indexed: 07/20/2024]
Abstract
Dense packing of particles has provided powerful models to elaborate the important structural features of matter in various systems such as liquid, glassy, and crystalline phases. The simplest sphere packing models can represent and capture salient properties of the building blocks for covalent, metallic, and ionic crystals; it, however, becomes insufficient to reflect the broken symmetry of the commonly anisotropic molecules in molecular crystals. Here, we develop spheroid models with a minimal degree of anisotropy, which serve as a simple geometrical representation for a rich spectrum of molecules-including both isotropic and anisotropic, convex and concave ones-in crystalline phases. Our models are determined via an inverse packing approach: Given a molecular crystal, an optimal spheroid model is constructed using a contact diagram, which depicts the packing relationship between neighboring molecules within the crystal. The spheroid models are capable of accurately capturing the broken symmetry and characterizing the equivalent volume of molecules in the crystalline phases. Moreover, our model retrieves such molecular information from low-quality x-ray diffraction data with poorly resolved structures, and by using soft spheroids, it can also describe the packing behavior in cocrystals.
Collapse
|
3
|
Schaller FM, Punzmann H, Schröder-Turk GE, Saadatfar M. Mixing properties of bi-disperse ellipsoid assemblies: mean-field behaviour in a granular matter experiment. SOFT MATTER 2023; 19:951-958. [PMID: 36633168 DOI: 10.1039/d2sm00922f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
The structure and spatial statistical properties of amorphous ellipsoid assemblies have profound scientific and industrial significance in many systems, from cell assays to granular materials. This paper uses a fundamental theoretical relationship for mixture distributions to explain the observations of an extensive X-ray computed tomography study of granular ellipsoidal packings. We study a size-bi-disperse mixture of two types of ellipsoids of revolutions that have the same aspect ratio of α ≈ 0.57 and differ in size, by about 10% in linear dimension, and compare these to mono-disperse systems of ellipsoids with the same aspect ratio. Jammed configurations with a range of packing densities are achieved by employing different tapping protocols. We numerically interrogate the final packing configurations by analyses of the local packing fraction distributions calculated from the Voronoi diagrams. Our main finding is that the bi-disperse ellipsoidal packings studied here can be interpreted as a mixture of two uncorrelated mono-disperse packings, insensitive to the compaction protocol. Our results are consolidated by showing that the local packing fraction shows no correlation beyond their first shell of neighbours in the binary mixtures. We propose a model of uncorrelated binary mixture distribution that describes the observed experimental data with high accuracy. This analysis framework will enable future studies to test whether the observed mean-field behaviour is specific to the particular granular system or the specific parameter values studied here or if it is observed more broadly in other bi-disperse non-spherical particle systems.
Collapse
Affiliation(s)
- F M Schaller
- Friedrich-Alexander Universität Erlangen-Nürnberg, Institut für Theoretische Physik, Staudtstr. 7B, 91058 Erlangen, Germany.
- Karlsruhe Institute of Technology (KIT), Institut für Stochastik, 76131 Karlsruhe, Germany
| | - H Punzmann
- The Australian National University, Research School of Physics, Canberra ACT 2601, Australia
| | - G E Schröder-Turk
- The Australian National University, Research School of Physics, Canberra ACT 2601, Australia
- Murdoch University, College of Science, Technology, Engineering and Mathematics, 90 South St, Murdoch WA 6150, Australia
| | - M Saadatfar
- The Australian National University, Research School of Physics, Canberra ACT 2601, Australia
- The University of Sydney, School of Civil Engineering, NSW 2006, Australia.
| |
Collapse
|
4
|
Phua A, Smith J, Davies CH, Cook PS, Delaney GW. Understanding the structure and dynamics of local powder packing density variations in metal additive manufacturing using set Voronoi analysis. POWDER TECHNOL 2023. [DOI: 10.1016/j.powtec.2023.118272] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
|
5
|
Iniyatova G, Yermukhambetova A, Boribayeva A, Golman B. Approximate Packing of Binary Mixtures of Cylindrical Particles. MICROMACHINES 2022; 14:36. [PMID: 36677097 PMCID: PMC9862688 DOI: 10.3390/mi14010036] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2022] [Revised: 12/14/2022] [Accepted: 12/20/2022] [Indexed: 06/17/2023]
Abstract
Particle packing plays an essential role in industry and chemical engineering. In this work, the discrete element method is used to generate the cylindrical particles and densify the binary cylindrical particle mixtures under the poured packing conditions. The influences of the aspect ratio and volume fraction of particles on the packing structure are measured by planar packing fraction. The Voronoi tessellation is used to quantify the porous structure of packing. The cumulative distribution functions of local packing fractions and the probability distributions of the reduced free volume of Voronoi cells are calculated to describe the local packing characteristics of binary mixtures with different volume fractions. As a result, it is observed that particles with larger aspect ratios in the binary mixture tend to orient randomly, and the particles with smaller aspect ratios have a preferentially horizontal orientation. Results also show that the less dense packings are obtained for mixtures with particles of higher aspect ratios and mixtures with a larger fraction of elongated cylindrical particles.
Collapse
|
6
|
Zhang C, Zhao S, Zhao J, Zhou X. Three-dimensional Voronoi analysis of realistic grain packing: An XCT assisted set Voronoi tessellation framework. POWDER TECHNOL 2021. [DOI: 10.1016/j.powtec.2020.10.054] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
7
|
Voggenreiter M, Roller J, Geiger J, Ebner L, Zumbusch A, Meijer JM. Preparation and Tracking of Oblate Core-Shell Polymethyl-Methacrylate Ellipsoids. LANGMUIR : THE ACS JOURNAL OF SURFACES AND COLLOIDS 2020; 36:13087-13095. [PMID: 33085481 DOI: 10.1021/acs.langmuir.0c02597] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
Although single-particle level studies on prolate ellipsoidal colloids are relatively abundant, similar studies on oblate ellipsoids are rare because suitable model systems are scarcely available. Here, we present the preparation of monodisperse hard core-shell oblate ellipsoids that can be imaged and tracked in 3D with confocal laser scanning microscopy. Using a thermomechanical squeezing method, we transform spherical core-shell polymethyl-methacrylate (PMMA) particles into oblate ellipsoids. We show how the shape polydispersity as well as the aspect ratio of the obtained oblate ellipsoids can be controlled. In addition, we discuss how the core-shell geometry limits the range of aspect ratios because of the different viscoelastic properties of the cross-linked PMMA core and linear PMMA shell. We further demonstrate imaging of the core-shell oblate dispersions on a single-particle level in real space and time and the tracking of position and orientation using our recently developed tracking algorithm for anisotropic core-shell colloids. Our results thus provide the tools for the future investigation of the behavior of oblate ellipsoids, especially in dense suspensions.
Collapse
Affiliation(s)
- Markus Voggenreiter
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany
| | - Jörg Roller
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany
| | - John Geiger
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany
| | - Lukas Ebner
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany
| | - Andreas Zumbusch
- Department of Chemistry, University of Konstanz, Universitätsstraße 10, 78457 Konstanz, Germany
| | - Janne-Mieke Meijer
- Department of Applied Physics, Eindhoven University of Technology, P.O. Box 513, 5600 MB Eindhoven, The Netherlands
| |
Collapse
|
8
|
Kim Y, Lee S, Lim J, Weon BM. X-ray nanotomography of dry colloidal packings. Sci Rep 2020; 10:17222. [PMID: 33057174 PMCID: PMC7560702 DOI: 10.1038/s41598-020-74315-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Accepted: 09/25/2020] [Indexed: 11/18/2022] Open
Abstract
Random packings are crucial in understanding arrangement and geometry of particles. Random packings of dry small particles may be subject to adhesion or friction, as expected theoretically and numerically. We explore experimentally random packings of dry colloids with X-ray nanotomography that directly provides three-dimensional structural and geometric information of dry colloidal packings. We find that dry colloidal packings, as characterized by contact number and packing density, are quite consistent with adhesive loose packings that significantly deviate from random loose packings for hard spheres. This study may offer direct evidence for adhesive loose packings comprising dry small particles, as proven by X-ray nanotomography.
Collapse
Affiliation(s)
- Yeseul Kim
- Soft Matter Physics Laboratory, SKKU Advanced Institute of Nanotechnology (SAINT), School of Advanced Materials Science and Engineering, Sungkyunkwan University, Suwon, Gyeonggi, 16419, South Korea
| | - Sangsul Lee
- Industrial Technology Convergence Center, Pohang Light Source, Pohang, Gyeongbuk, 37673, South Korea
| | - Jun Lim
- Industrial Technology Convergence Center, Pohang Light Source, Pohang, Gyeongbuk, 37673, South Korea
| | - Byung Mook Weon
- Soft Matter Physics Laboratory, SKKU Advanced Institute of Nanotechnology (SAINT), School of Advanced Materials Science and Engineering, Sungkyunkwan University, Suwon, Gyeonggi, 16419, South Korea. .,Research Center for Advanced Materials Technology, Sungkyunkwan University, Suwon, Gyeonggi, 16419, South Korea. .,Department of Biomedical Engineering, Johns Hopkins University, Baltimore, MD, 21218, USA.
| |
Collapse
|
9
|
Liu L, Li S. Uniform shape elongation effects on the random packings of uniaxially variable superellipsoids. POWDER TECHNOL 2020. [DOI: 10.1016/j.powtec.2020.08.017] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
10
|
Yuan Y, Deng W, Li S. Structural universality in disordered packings with size and shape polydispersity. SOFT MATTER 2020; 16:4528-4539. [PMID: 32356543 DOI: 10.1039/d0sm00110d] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We numerically investigate disordered jammed packings with both size and shape polydispersity, using frictionless superellipsoidal particles. We implement the set Voronoi tessellation technique to evaluate the local specific volume, i.e., the ratio of cell volume over particle volume, for each individual particle. We focus on the average structural properties for different types of particles binned by their sizes and shapes. We generalize the basic observation that the larger particles are locally packed more densely than the smaller ones in a polydisperse-sized packing into systems with coupled particle shape dispersity. For this purpose, we define the normalized free volume vf to measure the local compactness of a particle and study its dependency on the normalized particle size A. The definition of vf relies on the calibrated monodisperse specific volume for a certain particle shape. For packings with shape dispersity, we apply the previously introduced concept of equivalent diameter for a non-spherical particle to define A properly. We consider three systems: (A) linear superposition states of mixed-shape packings, (B) merely polydisperse-sized packings, and (C) packings with coupled size and shape polydispersity. For (A), the packing is simply considered as a mixture of different subsystems corresponding to monodisperse packings for different shape components, leading to A = 1, and vf = 1 by definition. We propose a concise model to estimate the shape-dependent factor αc, which defines the equivalent diameter for a certain particle. For (B), vf collapses as a function of A, independent of specific particle shape and size polydispersity. Such structural universality is further validated by a mean-field approximation. For (C), we find that the master curve vf(A) is preserved when particles possess similar αc in a packing. Otherwise, the dispersity of αc among different particles causes the deviation from vf(A). These findings show that a polydisperse packing can be estimated as the combination of various building blocks, i.e., bin components, with a universal relation vf(A).
Collapse
Affiliation(s)
- Ye Yuan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China.
| | - Wei Deng
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China.
| | - Shuixiang Li
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China.
| |
Collapse
|
11
|
Zhao S, Zhao J, Guo N. Universality of internal structure characteristics in granular media under shear. Phys Rev E 2020; 101:012906. [PMID: 32069573 DOI: 10.1103/physreve.101.012906] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2019] [Indexed: 11/07/2022]
Abstract
We examine the signatures of internal structure emerged from quasistatic shear responses of granular materials based on three-dimensional discrete element simulations. Granular assemblies consisting of spheres or nonspherical particles of different polydispersity are sheared from different initial densities and under different loading conditions (drained or undrained) steadily to reach the critical state (a state featured by constant stress and constant volume). The radial distribution function used to measure the packing structure is found to remain almost unchanged during the shearing process, regardless of the initial states or loading conditions of an assembly. Its specific form, however, varies with polydispersities in both grain size and grain shape. Set Voronoi tessellation is employed to examine the characteristics of local volume and anisotropy, and deformation. The local inverse solid fraction and anisotropy index are found following inverse Weibull and log-normal distributions, respectively, which are unique at the critical states. With further normalization, an invariant distribution for local volume and anisotropy is observed whose function can be determined by the polydispersities in both particle size and grain shape but bears no relevance to initial densities or loading conditions (or paths). An invariant Gaussian distribution is found for the local deformation for spherical packings, but no invariant distribution can be found for nonspherical packings where the distribution of normalized local volumetric strain is dependent on initial states.
Collapse
Affiliation(s)
- Shiwei Zhao
- Department of Civil and Environmental Engineering, The Hong Kong University of Science and Technology, Clearwater Bay, Kowloon, Hong Kong
| | - Jidong Zhao
- Department of Civil and Environmental Engineering, The Hong Kong University of Science and Technology, Clearwater Bay, Kowloon, Hong Kong
| | - Ning Guo
- College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
| |
Collapse
|
12
|
Yuan Y, VanderWerf K, Shattuck MD, O'Hern CS. Jammed packings of 3D superellipsoids with tunable packing fraction, coordination number, and ordering. SOFT MATTER 2019; 15:9751-9761. [PMID: 31742301 PMCID: PMC6902436 DOI: 10.1039/c9sm01932d] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/09/2023]
Abstract
We carry out numerical studies of static packings of frictionless superellipsoidal particles in three spatial dimensions. We consider more than 200 different particle shapes by varying the three shape parameters that define superellipsoids. We characterize the structural and mechanical properties of both disordered and ordered packings using two packing-generation protocols. We perform athermal quasi-static compression simulations starting from either random, dilute configurations (Protocol 1) or thermalized, dense configurations (Protocol 2), which allows us to tune the orientational order of the packings. In general, we find that superellipsoid packings are hypostatic, with coordination number zJ < ziso, where ziso = 2df and df = 5 or 6 depending on whether the particles are axi-symmetric or not. Over the full range of orientational order, we find that the number of quartic modes of the dynamical matrix for the packings always matches the number of missing contacts relative to the isostatic value. This result suggests that there are no mechanically redundant contacts for ordered, yet hypostatic packings of superellipsoidal particles. Additionally, we find that the packing fraction at jamming onset for disordered packings of superellipsoidal depends on at least two particle shape parameters, e.g. the asphericity A and reduced aspect ratio β of the particles.
Collapse
Affiliation(s)
- Ye Yuan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China. and Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, USA
| | - Kyle VanderWerf
- Department of Physics, Yale University, New Haven, Connecticut 06520, USA
| | - Mark D Shattuck
- Benjamin Levich Institute and Physics Department, The City College of New York, New York, New York 10031, USA
| | - Corey S O'Hern
- Department of Mechanical Engineering and Materials Science, Yale University, New Haven, Connecticut 06520, USA and Department of Physics, Yale University, New Haven, Connecticut 06520, USA and Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA and Graduate Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
13
|
Lovrić J, Kaliman S, Barfuss W, Schröder-Turk GE, Smith AS. Geometric effects in random assemblies of ellipses. SOFT MATTER 2019; 15:8566-8577. [PMID: 31637393 DOI: 10.1039/c9sm01067j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/10/2023]
Abstract
Assemblies of anisotropic particles commonly appear in studies of active many-body systems. However, in two dimensions, the geometric ramifications of the finite density of such objects are not entirely understood. To fully characterize these effects, we perform an in-depth study of random assemblies generated by a slow compression of frictionless elliptical particles. The obtained configurations are then analysed using the Set Voronoi tessellation, which takes the particle shape into account. Not only do we analyse most scalar and vectorial morphological measures, which are commonly discussed in the literature or which have recently been addressed in experiments, but we also systematically explore the correlations between them. While in a limited range of parameters similarities with findings in 3D assemblies could be identified, important differences are found when a broad range of aspect ratios and packing fractions are considered. The data discussed in this study should thus provide a unique reference set such that geometric effects and differences from random assemblies could be clearly identified in more complex systems, including ones with soft and active particles that are typically found in biological systems.
Collapse
Affiliation(s)
- Jakov Lovrić
- Division of Physical Chemistry, Ruđer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
| | | | | | | | | |
Collapse
|
14
|
Shiraishi K, Mizuno H, Ikeda A. Vibrational properties of two-dimensional dimer packings near the jamming transition. Phys Rev E 2019; 100:012606. [PMID: 31499851 DOI: 10.1103/physreve.100.012606] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/10/2019] [Indexed: 11/07/2022]
Abstract
Jammed particulate systems composed of various shapes of particles undergo the jamming transition as they are compressed or decompressed. To date, sphere packings have been extensively studied in many previous works, where isostaticity at the transition and scaling laws with the pressure of various quantities, including the contact number and the vibrational density of states, have been established. Additionally, much attention has been paid to nonspherical packings, and particularly recent work has made progress in understanding ellipsoidal packings. In this work, we study the dimer packings in two dimensions, which have been much less understood than systems of spheres and ellipsoids. We first study the contact number of dimers near the jamming transition. It turns out that packings of dimers have "rotational rattlers," each of which still has a free rotational motion. After correcting this effect, we show that dimers become isostatic at the jamming, and the excess contact number obeys the same critical law and finite-size scaling law as those of spheres. We next study the vibrational properties of dimers near the transition. We find that the vibrational density of states of dimers exhibits two characteristic plateaus that are separated by a peak. The high-frequency plateau is dominated by the translational degree of freedom, while the low-frequency plateau is dominated by the rotational degree of freedom. We establish the critical scaling laws of the characteristic frequencies of the plateaus and the peak near the transition. In addition, we present detailed characterizations of the real space displacement fields of vibrational modes in the translational and rotational plateaus.
Collapse
Affiliation(s)
- Kumpei Shiraishi
- Graduate School of Arts and Sciences, University of Tokyo, Komaba, Tokyo 153-8902, Japan
| | - Hideyuki Mizuno
- Graduate School of Arts and Sciences, University of Tokyo, Komaba, Tokyo 153-8902, Japan
| | - Atsushi Ikeda
- Graduate School of Arts and Sciences, University of Tokyo, Komaba, Tokyo 153-8902, Japan.,Research Center for Complex Systems Biology, Universal Biology Institute, University of Tokyo, Komaba, Tokyo 153-8902, Japan
| |
Collapse
|
15
|
Liu L, Yu Z, Jin W, Yuan Y, Li S. Uniform and decoupled shape effects on the maximally dense random packings of hard superellipsoids. POWDER TECHNOL 2018. [DOI: 10.1016/j.powtec.2018.06.033] [Citation(s) in RCA: 22] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
16
|
Liu L, Yuan Y, Deng W, Li S. Evolutions of packing properties of perfect cylinders under densification and crystallization. J Chem Phys 2018; 149:104503. [PMID: 30218998 DOI: 10.1063/1.5049562] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Cylindrical particles are ubiquitous in nature and industry, and a cylinder is a representative shape of rod-like particles. However, the disordered packing results of cylinders in previous studies are quite inconsistent with each other. In this work, we obtain the MRJ (maximally random jammed) packings and the MDRPs (maximally dense random packings) of perfect cylinders with the aspect ratio (height/diameter) 0.2 ≤ w ≤ 6.0 using the ASC (adaptive shrinking cell) algorithm and the IMC (inverse Monte Carlo) method, respectively. The optimal aspect ratio corresponding to the maximal packing density is w = 0.9 in the MRJ state, while the value is w = 1.2 in the MDRP state. Then we investigate the evolutions of packing properties of perfect cylinders under densification and crystallization. We compare the different final packing states generated via the two methods with different compression rates and order constraints. In the densification procedure, we generate jammed and random packings of cylinders with various compression rates via the ASC and IMC method, respectively. When decreasing the compression rate, we find that the packing density increases but the optimal w remains the same in both methods. In the crystallization procedure, the order constraint in the IMC method is gradually released which means the degree of order in the packings is allowed to increase, and we find that the optimal w shifts from 1.2 to 0.9 while the packing density increases as well. Meanwhile, the random packings evolve into the jammed packings in the crystallization procedure which reflects the competition mechanism between the randomness and jamming. These results also indicate that the optimal w is solely related to the degree of order in the cylinder packings but not determined by the protocol or packing density. Furthermore, a uniform shape elongation effect on the random-packing densities of various shaped particles is found via a new proposed definition of the scaled aspect ratio. Finally, a rough linear relationship between the mean and standard deviation of the reduced Voronoi cell volumes is obtained only for the random packings. Our findings should lead to a better understanding toward the jammed and random packings and are helpful in guiding the granular material design.
Collapse
Affiliation(s)
- Lufeng Liu
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China
| | - Ye Yuan
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China
| | - Wei Deng
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China
| | - Shuixiang Li
- Department of Mechanics and Engineering Science, College of Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
17
|
Wagner JW, Dannenhoffer-Lafage T, Jin J, Voth GA. Extending the range and physical accuracy of coarse-grained models: Order parameter dependent interactions. J Chem Phys 2018; 147:044113. [PMID: 28764380 DOI: 10.1063/1.4995946] [Citation(s) in RCA: 47] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Order parameters (i.e., collective variables) are often used to describe the behavior of systems as they capture different features of the free energy surface. Yet, most coarse-grained (CG) models only employ two- or three-body non-bonded interactions between the CG particles. In situations where these interactions are insufficient for the CG model to reproduce the structural distributions of the underlying fine-grained (FG) model, additional interactions must be included. In this paper, we introduce an approach to expand the basis sets available in the multiscale coarse-graining (MS-CG) methodology by including order parameters. Then, we investigate the ability of an additive local order parameter (e.g., density) and an additive global order parameter (i.e., distance from a hard wall) to improve the description of CG models in interfacial systems. Specifically, we study methanol liquid-vapor coexistence, acetonitrile liquid-vapor coexistence, and acetonitrile liquid confined by hard-wall plates, all using single site CG models. We find that the use of order parameters dramatically improves the reproduction of structural properties of interfacial CG systems relative to the FG reference as compared with pairwise CG interactions alone.
Collapse
Affiliation(s)
- Jacob W Wagner
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Thomas Dannenhoffer-Lafage
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Jaehyeok Jin
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| | - Gregory A Voth
- Department of Chemistry, James Franck Institute, and Institute for Biophysical Dynamics, The University of Chicago, Chicago, Illinois 60637, USA
| |
Collapse
|
18
|
|
19
|
VanderWerf K, Jin W, Shattuck MD, O'Hern CS. Hypostatic jammed packings of frictionless nonspherical particles. Phys Rev E 2018; 97:012909. [PMID: 29448406 PMCID: PMC6295208 DOI: 10.1103/physreve.97.012909] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/25/2017] [Indexed: 11/07/2022]
Abstract
We perform computational studies of static packings of a variety of nonspherical particles including circulo-lines, circulo-polygons, ellipses, asymmetric dimers, dumbbells, and others to determine which shapes form packings with fewer contacts than degrees of freedom (hypostatic packings) and which have equal numbers of contacts and degrees of freedom (isostatic packings), and to understand why hypostatic packings of nonspherical particles can be mechanically stable despite having fewer contacts than that predicted from naive constraint counting. To generate highly accurate force- and torque-balanced packings of circulo-lines and cir-polygons, we developed an interparticle potential that gives continuous forces and torques as a function of the particle coordinates. We show that the packing fraction and coordination number at jamming onset obey a masterlike form for all of the nonspherical particle packings we studied when plotted versus the particle asphericity A, which is proportional to the ratio of the squared perimeter to the area of the particle. Further, the eigenvalue spectra of the dynamical matrix for packings of different particle shapes collapse when plotted at the same A. For hypostatic packings of nonspherical particles, we verify that the number of "quartic" modes along which the potential energy increases as the fourth power of the perturbation amplitude matches the number of missing contacts relative to the isostatic value. We show that the fourth derivatives of the total potential energy in the directions of the quartic modes remain nonzero as the pressure of the packings is decreased to zero. In addition, we calculate the principal curvatures of the inequality constraints for each contact in circulo-line packings and identify specific types of contacts with inequality constraints that possess convex curvature. These contacts can constrain multiple degrees of freedom and allow hypostatic packings of nonspherical particles to be mechanically stable.
Collapse
Affiliation(s)
- Kyle VanderWerf
- Department of Physics, Yale University, New Haven, Connecticut 06520, USA
| | - Weiwei Jin
- Department of Mechanics and Engineering Science, Peking University, Beijing 100871, China
- Department of Mechanical Engineering & Materials Science, Yale University, New Haven, Connecticut 06520, USA
| | - Mark D Shattuck
- Benjamin Levich Institute and Physics Department, The City College of New York, New York, New York 10031, USA
| | - Corey S O'Hern
- Department of Physics, Yale University, New Haven, Connecticut 06520, USA
- Department of Mechanical Engineering & Materials Science, Yale University, New Haven, Connecticut 06520, USA
- Department of Applied Physics, Yale University, New Haven, Connecticut 06520, USA
- Graduate Program in Computational Biology and Bioinformatics, Yale University, New Haven, Connecticut 06520, USA
| |
Collapse
|
20
|
Schella A, Weis S, Schröter M. Charging changes contact composition in binary sphere packings. Phys Rev E 2017; 95:062903. [PMID: 28709334 DOI: 10.1103/physreve.95.062903] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2017] [Indexed: 11/07/2022]
Abstract
Equal volume mixtures of small and large polytetrafluorethylene spheres are shaken in an atmosphere of controlled humidity which allows one to also control their tribocharging. We find that the contact numbers are charge dependent: As the charge density of the beads increases, the number of same-type contacts decreases and the number of opposite-type contacts increases. This change is not caused by a global segregation of the sample. Hence, tribocharging can be a way to tune the local composition of a granular material.
Collapse
Affiliation(s)
- André Schella
- Max Planck Institute for Dynamics and Self-Organization Göttingen, 37077 Göttingen, Germany
| | - Simon Weis
- Institute of Theoretical Physics I, University Erlangen-Nürnberg, Staudtstraße 7, 91058 Erlangen, Germany
| | - Matthias Schröter
- Max Planck Institute for Dynamics and Self-Organization Göttingen, 37077 Göttingen, Germany.,Institute for Multiscale Simulation, Friedrich-Alexander University, 91052 Erlangen, Germany
| |
Collapse
|
21
|
|
22
|
Wang D, Zheng H, Behringer RP. A Granular System of Ellipses under Linear Shear. EPJ WEB OF CONFERENCES 2017. [DOI: 10.1051/epjconf/201714006003] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Weis S, Schröter M. Analyzing X-ray tomographies of granular packings. THE REVIEW OF SCIENTIFIC INSTRUMENTS 2017; 88:051809. [PMID: 28571396 DOI: 10.1063/1.4983051] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Starting from three-dimensional volume data of a granular packing, as, e.g., obtained by X-ray Computed Tomography, we discuss methods to first detect the individual particles in the sample and then analyze their properties. This analysis includes the pair correlation function, the volume and shape of the Voronoi cells, and the number and type of contacts formed between individual particles. We mainly focus on packings of monodisperse spheres, but we will also comment on other monoschematic particles such as ellipsoids and tetrahedra. This paper is accompanied by a package of free software containing all programs (including source code) and an example three-dimensional dataset which allows the reader to reproduce and modify all examples given.
Collapse
Affiliation(s)
- Simon Weis
- Institut für Theoretische Physik I, Friedrich-Alexander-Universität, 91058 Erlangen, Germany
| | - Matthias Schröter
- Institute for Multiscale Simulation, Friedrich-Alexander-Universität, 91052 Erlangen, Germany
| |
Collapse
|
24
|
Jin W, Jiao Y, Liu L, Yuan Y, Li S. Dense crystalline packings of ellipsoids. Phys Rev E 2017; 95:033003. [PMID: 28415357 DOI: 10.1103/physreve.95.033003] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Indexed: 11/07/2022]
Abstract
An ellipsoid, the simplest nonspherical shape, has been extensively used as a model for elongated building blocks for a wide spectrum of molecular, colloidal, and granular systems. Yet the densest packing of congruent hard ellipsoids, which is intimately related to the high-density phase of many condensed matter systems, is still an open problem. We discover an unusual family of dense crystalline packings of self-dual ellipsoids (ratios of the semiaxes α:sqrt[α]:1), containing 24 particles with a quasi-square-triangular (SQ-TR) tiling arrangement in the fundamental cell. The associated packing density ϕ exceeds that of the densest known SM2 crystal [ A. Donev et al., Phys. Rev. Lett. 92, 255506 (2004)10.1103/PhysRevLett.92.255506] for aspect ratios α in (1.365, 1.5625), attaining a maximal ϕ≈0.75806... at α=93/64. We show that the SQ-TR phase derived from these dense packings is thermodynamically stable at high densities over the aforementioned α range and report a phase diagram for self-dual ellipsoids. The discovery of the SQ-TR crystal suggests organizing principles for nonspherical particles and self-assembly of colloidal systems.
Collapse
Affiliation(s)
- Weiwei Jin
- Department of Mechanics and Engineering Science, Peking University, Beijing 100871, China
| | - Yang Jiao
- Materials Science and Engineering, Arizona State University, Tempe, Arizona 85287, USA
| | - Lufeng Liu
- Department of Mechanics and Engineering Science, Peking University, Beijing 100871, China
| | - Ye Yuan
- Department of Mechanics and Engineering Science, Peking University, Beijing 100871, China
| | - Shuixiang Li
- Department of Mechanics and Engineering Science, Peking University, Beijing 100871, China
| |
Collapse
|
25
|
|
26
|
Schaller FM, Punzmann H, Schröder-Turk GE, Saadatfar M. Towards minimal models for realistic granular materials: Tomographic analysis of bidispersed assemblies of ellipsoids. EPJ WEB OF CONFERENCES 2017. [DOI: 10.1051/epjconf/201714006030] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
Delaney GW. Structure Generation and Analysis from the Loosest to the Densest Random Packings. EPJ WEB OF CONFERENCES 2017. [DOI: 10.1051/epjconf/201714006029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
28
|
Malmir H, Sahimi M, Tabar MRR. Packing of nonoverlapping cubic particles: Computational algorithms and microstructural characteristics. Phys Rev E 2016; 94:062901. [PMID: 28085418 DOI: 10.1103/physreve.94.062901] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2016] [Indexed: 06/06/2023]
Abstract
Packing of cubic particles arises in a variety of problems, ranging from biological materials to colloids and the fabrication of new types of porous materials with controlled morphology. The properties of such packings may also be relevant to problems involving suspensions of cubic zeolites, precipitation of salt crystals during CO_{2} sequestration in rock, and intrusion of fresh water in aquifers by saline water. Not much is known, however, about the structure and statistical descriptors of such packings. We present a detailed simulation and microstructural characterization of packings of nonoverlapping monodisperse cubic particles, following up on our preliminary results [H. Malmir et al., Sci. Rep. 6, 35024 (2016)2045-232210.1038/srep35024]. A modification of the random sequential addition (RSA) algorithm has been developed to generate such packings, and a variety of microstructural descriptors, including the radial distribution function, the face-normal correlation function, two-point probability and cluster functions, the lineal-path function, the pore-size distribution function, and surface-surface and surface-void correlation functions, have been computed, along with the specific surface and mean chord length of the packings. The results indicate the existence of both spatial and orientational long-range order as the the packing density increases. The maximum packing fraction achievable with the RSA method is about 0.57, which represents the limit for a structure similar to liquid crystals.
Collapse
Affiliation(s)
- Hessam Malmir
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089-1211, USA
| | - Muhammad Sahimi
- Mork Family Department of Chemical Engineering and Materials Science, University of Southern California, Los Angeles, California 90089-1211, USA
| | - M Reza Rahimi Tabar
- Department of Physics, Sharif University of Technology, Tehran 11365-9161, Iran
| |
Collapse
|
29
|
|
30
|
Baranau V, Zhao SC, Scheel M, Tallarek U, Schröter M. Upper bound on the Edwards entropy in frictional monodisperse hard-sphere packings. SOFT MATTER 2016; 12:3991-4006. [PMID: 27020114 DOI: 10.1039/c6sm00567e] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/05/2023]
Abstract
We extend the Widom particle insertion method [B. Widom, J. Chem. Phys., 1963, 39, 2808-2812] to determine an upper bound sub on the Edwards entropy in frictional hard-sphere packings. sub corresponds to the logarithm of the number of mechanically stable configurations for a given volume fraction and boundary conditions. To accomplish this, we extend the method for estimating the particle insertion probability through the pore-size distribution in frictionless packings [V. Baranau, et al., Soft Matter, 2013, 9, 3361-3372] to the case of frictional particles. We use computer-generated and experimentally obtained three-dimensional sphere packings with volume fractions φ in the range 0.551-0.65. We find that sub has a maximum in the vicinity of the Random Loose Packing Limit φRLP = 0.55 and decreases then monotonically with increasing φ to reach a minimum at φ = 0.65. Further on, sub does not distinguish between real mechanical stability and packings in close proximity to mechanical stable configurations. The probability to find a given number of contacts for a particle inserted in a large enough pore does not depend on φ, but it decreases strongly with the contact number.
Collapse
Affiliation(s)
- Vasili Baranau
- Department of Chemistry, Philipps-Universität Marburg, Hans-Meerwein-Strasse 4, D-35032 Marburg, Germany.
| | | | | | | | | |
Collapse
|
31
|
Wang C, Dong K, Yu A. Structural characterization of the packings of granular regular polygons. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 92:062203. [PMID: 26764678 DOI: 10.1103/physreve.92.062203] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/03/2015] [Indexed: 06/05/2023]
Abstract
By using a recently developed method for discrete modeling of nonspherical particles, we simulate the random packings of granular regular polygons with three to 11 edges under gravity. The effects of shape and friction on the packing structures are investigated by various structural parameters, including packing fraction, the radial distribution function, coordination number, Voronoi tessellation, and bond-orientational order. We find that packing fraction is generally higher for geometrically nonfrustrated regular polygons, and can be increased by the increase of edge number and decrease of friction. The changes of packing fraction are linked with those of the microstructures, such as the variations of the translational and orientational orders and local configurations. In particular, the free areas of Voronoi tessellations (which are related to local packing fractions) can be described by log-normal distributions for all polygons. The quantitative analyses establish a clearer picture for the packings of regular polygons.
Collapse
Affiliation(s)
- Chuncheng Wang
- Laboratory for Simulation and Modeling of Particulate Systems, Department of Chemical Engineering, Monash University, Clayton, Victoria 3800, Australia
| | - Kejun Dong
- Institute for Infrastructure Engineering, Western Sydney University, Penrith NSW 2751, Australia
| | - Aibing Yu
- Laboratory for Simulation and Modeling of Particulate Systems, Department of Chemical Engineering, Monash University, Clayton, Victoria 3800, Australia
| |
Collapse
|
32
|
Ferrage E, Hubert F, Tertre E, Delville A, Michot LJ, Levitz P. Modeling the arrangement of particles in natural swelling-clay porous media using three-dimensional packing of elliptic disks. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:062210. [PMID: 26172708 DOI: 10.1103/physreve.91.062210] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/15/2015] [Indexed: 06/04/2023]
Abstract
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.
Collapse
Affiliation(s)
- Eric Ferrage
- Université de Poitiers, CNRS, UMR 7285, Equipe HydrASA Laboratoire, Institut de Chimie des Milieux et Matériaux de Poitiers, F-86022 Poitiers, France
| | - Fabien Hubert
- Université de Poitiers, CNRS, UMR 7285, Equipe HydrASA Laboratoire, Institut de Chimie des Milieux et Matériaux de Poitiers, F-86022 Poitiers, France
| | - Emmanuel Tertre
- Université de Poitiers, CNRS, UMR 7285, Equipe HydrASA Laboratoire, Institut de Chimie des Milieux et Matériaux de Poitiers, F-86022 Poitiers, France
| | - Alfred Delville
- Université d'Orléans, CNRS, UMR 7374, Interfaces, Confinement, Matériaux et Nanostructures, F-45071 Orléans, France
| | - Laurent J Michot
- Université Pierre et Marie Curie, CNRS, UMR 8234, Physicochimie des Electrolytes et Nanosystèmes Interfaciaux, F- 72522 Paris, France
| | - Pierre Levitz
- Université Pierre et Marie Curie, CNRS, UMR 8234, Physicochimie des Electrolytes et Nanosystèmes Interfaciaux, F- 72522 Paris, France
| |
Collapse
|
33
|
Hanifpour M, Francois N, Robins V, Kingston A, Allaei SMV, Saadatfar M. Structural and mechanical features of the order-disorder transition in experimental hard-sphere packings. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2015; 91:062202. [PMID: 26172700 DOI: 10.1103/physreve.91.062202] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/28/2015] [Indexed: 06/04/2023]
Abstract
Here we present an experimental and numerical investigation on the grain-scale geometrical and mechanical properties of partially crystallized structures made of macroscopic frictional grains. Crystallization is inevitable in arrangements of monosized hard spheres with packing densities exceeding Bernal's limiting density ϕ(Bernal)≈0.64. We study packings of monosized hard spheres whose density spans over a wide range (0.59<ϕ<0.72). These experiments harness x-ray computed tomography, three-dimensional image analysis, and numerical simulations to access precisely the geometry and the 3D structure of internal forces within the sphere packings. We show that clear geometrical transitions coincide with modifications of the mechanical backbone of the packing both at the grain and global scale. Notably, two transitions are identified at ϕ(Bernal)≈0.64 and ϕ(c)≈0.68. These results provide insights on how geometrical and mechanical features at the grain scale conspire to yield partially crystallized structures that are mechanically stable.
Collapse
Affiliation(s)
- M Hanifpour
- Department of Physics, University of Tehran, Tehran 14395-547, Iran
- School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran
| | - N Francois
- Department of Applied Mathematics, Research School of Physics and Engineering, Australian National University, Canberra, Australia
| | - V Robins
- Department of Applied Mathematics, Research School of Physics and Engineering, Australian National University, Canberra, Australia
| | - A Kingston
- Department of Applied Mathematics, Research School of Physics and Engineering, Australian National University, Canberra, Australia
| | - S M Vaez Allaei
- Department of Physics, University of Tehran, Tehran 14395-547, Iran
- School of Physics, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran
| | - M Saadatfar
- Department of Applied Mathematics, Research School of Physics and Engineering, Australian National University, Canberra, Australia
| |
Collapse
|