Disks vs. spheres: Contrasting properties of random packings

Collections of random packings of rigid disks and spheres have been generated by computer using a previously described concurrent algorithm. Particles begin as infinitesimal moving points, grow in size at a uniform rate, undergo energy-onconserving collisions, and eventually jam up. Periodic boundary conditions apply, and various numbers of particles have been considered (N2000 for disks,N8000 for spheres). The irregular disk packings thus formed are clearly polycrystalline with mean grain size dependent upon particle growth rate. By contrast, the sphere packings show a homogeneously amorphous texture substantially devoid of crystalline grains. This distinction strongly influences the respective results for packing pair correlation functions and for the distributions of particles by contact number. Rapidly grown disk packings display occasional vacancies within the crystalline grains; no comparable voids of such distinctive size have been found in the random sphere packings. Rattler particles free to move locally but imprisoned by jammed neighbors occur in both the disk and sphere packings.This paper is dedicated to Jerry Percus on the occasion of his 65th birthday.     

直径任意分布球填充的数值模拟

提出球填充数值算法的新分类方法.改进原有的松弛算法,使其能够模拟直径任意分布的球填充问题,采用可变循环周期使不同球数情形下的填充率基本保持不变.算例数据表明,该算法的填充率和配位数均高于原算法.由于采用背景网格搜索和双向链表组数据结构,使得邻接球搜索效率有相当大的提高,算法的时间复杂度为O(N)(N为球数).在一台AMD Athlon 3200+PC上,对于10000个等径球的随机密排列,只需217s,填充率即可达到0.64.     

Dense packing in the monodisperse hard-sphere system: A numerical study

We report a numerical study of the close packing of monodisperse hard spheres. The close packings of hard spheres are produced by the Lubachesky-Stillinger (LS) compression algorithm and span the range from the disordered states to the ordered states. We provide quantitative evidence for the claim that the density and structural order of the arrested close packing can be determined by the compression rate, i.e., with slower rates producing denser and more ordered structures. Through deeply analyzing the structure of the resulting arrested close packings, a transition region has been identified in the plane of density and reciprocal compression rate, in between what have been historically thought of as amorphous and crystalline packings. We also find clear system size dependences in studying the structural properties of the packings from the disordered ones to the ordered ones. These detailed investigations, on the structure of the arrested close packings, may provide a link between the glassy states and the crystalline states in the hard spheres.     

Quantum corrections to thermodynamics of polar hard sphere fluids

The quantum corrections to the thermodynamic properties of polar hard sphere fluids and fluid mixtures are estimated taking into account the influence of dipole and quadrupole moments. Expressions are given for the second virial coefficient, free energy and pressure and results are given for different values ofμ* andϑ*. The first order quantum correction arises due to the translational contribution only. The quantum effect increases with density,μ* andϑ*. Numerical results are also estimated for binary mixtures of (i) hard spheres and dipole hard spheres and (ii) hard spheres and quadrupole hard spheres. The ‘excess’ free energy for dipole hard sphere binary mixture is also reported. It is found that the ‘excess’ quantum effect depends on the concentration and the particle diameter ratio and increases with increase ofμ* andϑ*.     

Shear stress correlations in hard and soft sphere fluids

The shear stress autocorrelation function has been studied by molecular dynamics simulation using the q^?n potential for very large n. The results are analysed and interpreted here by comparing them with the shear stress response function for hard spheres. It is shown that the hard sphere response function has a singular contribution, and that this is reproduced accurately by the simulations for large n. A simple model for the stress autocorrelation function at finite n is proposed, based on the required hard sphere limiting form.     

Numerical Simulation of Random Close Packings in Particle Deformation from Spheres to Cubes

Variation of packing density in particle deforming from spheres to cubes is studied. A new model is presented to describe particle deformation between different particle shapes. Deformation is simulated by relative motion of component spheres in the sphere assembly model of a particle. Random close packings of particles in deformation form spheres to cubes are simulated with an improved relaxation algorithm. Packings in both 2D and 3D cases are simulated. With the simulations, we find that the packing density increases while the particle sphericity decreases in the deformation. Spheres and cubes give the minimum (0.6404) and maximum (0.7755) of packing density in the deformation respectively. In each deforming step, packings starting from a random configuration and from the final packing of last deforming step are both simulated. The packing density in the latter case is larger than the former in two dimensions, but is smaller in three dimensions. The deformation model can be applied to other particle shapes as well.     

Square-well diatomics Exact low density results

Exact semi-analytical expressions are obtained for the zero density site-site distribution function and the second virial coefficient for homonuclear square-well diatomics. The diatomic molecules considered here are composed of two fused square-well spheres with hard core diameter σ, well diameter λσ, and dimer bond length L, such that 0<Lˇ-σ and 1ˇ-λˇ-1+L/σ. Very accurate, although approximate, analytical expressions are also given for these functions.     

Statistical verification of crystallization in hard sphere packings under densification

The performance of various structure characteristics in the task of indicating structural peculiarities in packings of hard spheres is investigated. Various characteristics based on Voronoi polyhedra, spherical harmonics, and Delaunay simplices are considered together with the pair correlation function and the mean number of r-close triples. They are applied to a set of hard sphere packings of density φ from 0.62 to 0.72. It turns out that all used structure characteristics are able to indicate changes of order from non-crystalline to crystalline packings. However, not all of them are sensitive enough to indicate different stages of structure transformation under densification. The characteristics based on Delaunay simplices turn out to be the most sensitive for this purpose. For the models considered three principal structure classes are found: packings of densities lower than the known critical value 0.64 showing a non-crystalline behavior; packings with considerable crystalline regions for φ up to 0.66–0.67; rather complete crystals although with numerous defects for φ above 0.67.     

Numerical Simulation of Random Close Packing with Tetrahedra

The densest packing of tetrahedra is still an unsolved problem. Numerical simulations of random close packing of tetrahedra are carried out with a sphere assembly model and improved relaxation algorithm. The packing density and average contact number obtained for random close packing of regular tetrahedra is 0.6817 and 7.21 respectively, while the values of spheres are 0.6435 and 5.95. The simulation demonstrates that tetrahedra can be randomly packed denser than spheres. Random close packings of tetrahedra with a range of height are simulated as well. We find that the regular tetrahedron might be the optimal shape which gives the highest packing density of tetrahedra.     

Dense and nearly jammed random packings of freely jointed chains of tangent hard spheres

Dense packings of freely jointed chains of tangent hard spheres are produced by a novel Monte Carlo method. Within statistical uncertainty, chains reach a maximally random jammed (MRJ) state at the same volume fraction as packings of single hard spheres. A structural analysis shows that as the MRJ state is approached (i) the radial distribution function for chains remains distinct from but approaches that of single hard sphere packings quite closely, (ii) chains undergo progressive collapse, and (iii) a small but increasing fraction of sites possess highly ordered first coordination shells.     

Monte Carlo simulation of Hamaker nanospheres coated with dipolar particles

Parallel tempering Monte Carlo simulation is carried out in systems of N attractive Hamaker spheres dressed with n dipolar particles, able to move on the surface of the spheres. Different cluster configurations emerge for given values of the control parameters. Energy per sphere, pair distribution functions of spheres and dipoles as function of temperature, density, external electric field, and/or the angular orientation of dipoles are used to analyse the state of aggregation of the system. As a consequence of the non-central interaction, the model predicts complex structures like self-assembly of spheres by a double crown of dipoles. This interesting result could be of help in understanding some recent experiments in colloidal science and biology.     

Comparison of two representations of a random cut of identical sphere packing

We compare two-dimensional froths obtained by radical tessellation of random planar cuts made through disordered assemblies of monosize spheres at different packing fractions C from C=0 to C=0.64 with two-dimensional stereological cuts obtained through the three-dimensional froths made with the same packing. We have built numerically the packings using different algorithms. The study of both topological and metric properties shows significant differences between the two representations. Received 26 May 1999 and Received in final form 13 November 1999     

Phyllotactic description of hard sphere packing in cylindrical channels

We develop a simple analytical theory that relates dense sphere packings in a cylinder to corresponding disk packings on its surface. It applies for ratios R=D/d (where d and D are the diameters of the hard spheres and the bounding cylinder, respectively) up to R=1+1/sin(π/5). Within this range the densest packings are such that all spheres are in contact with the cylindrical boundary. The detailed results elucidate extensive numerical simulations by ourselves and others by identifying the nature of various competing phases.     

Regular exact models for a nonstatic gas sphere in general relativity

Some new exact models for an expanding or a contracting gaseous sphere (i.e., the density is to vanish at the outer boundary together with the pressurep) are given. The physical properties of the models are investigated, and it is found that both the pressure and the density are positive inside the outer boundary of the sphere, and their respective gradients are negative. The density is increasing for contracting spheres, and it is decreasing for expanding spheres. It is also shown that this is the case for the pressure at any moment for the layers close to the boundary of the spheres. For these layers it is further shown that the adiabatic speed of sound is less than the speed of light, and the trace of the energy-momentum tensor is positive. The rate of change of the circumference as measured by an observer riding on the boundary of the sphere is increasing for expanding spheres and it is decreasing for collapsing spheres. We also find that the physical radius is an increasing function of comoving radial coordinate. The mass function is further shown to be positive.     

Polytetrahedral nature of the dense disordered packings of hard spheres

We study the structure of numerically simulated hard sphere packings at different densities by investigating local tetrahedral configurations of the spheres. Clusters of tetrahedra adjacent by faces present relatively dense aggregates of spheres atypical for crystals. The number of spheres participating in such polytetrahedral configurations increases with densification of the packing, and at the Bernal's limiting density (the packing fraction around 0.64) all spheres of the packing become involved in such tetrahedra. Thus the polytetrahedral packing cannot provide further increase in the density, and alternative structural change (formation of crystalline nuclei) begins henceforth.     

Apollonian packings as physical fractals

The Apollonian packings (APs) of spheres are fractals that result from a space-filling procedure. We discuss the finite size effects for finite intervals s?∈?[s _min,?s _max] between the largest and the smallest sizes of the filling spheres. We derive a simple analytical generalization of the scale-free laws, which allows a quantitative study of such physical fractals. To test our result, a new efficient space-filling algorithm has been developed which generates random APs of spheres with a finite range of diameters: the correct asymptotic limit s _min/s _max?→?0 and the known APs' fractal dimensions are recovered and an excellent agreement with the generalized analytic laws is proved within the overall range of sizes.     

Gravity driven instabilities in miscible non-Newtonian fluid displacements in porous media

Gravity driven instabilities in model porous packings of 1 mm diameter spheres are studied by comparing the broadening of the displacement front between fluids of slightly different densities in stable and unstable configurations. Water, water–glycerol and water–polymer solutions are used to vary independently viscosity and molecular diffusion and study the influence of shear-thinning properties. Both injected and displaced solutions are identical but for a different concentration of NaNO₃ salt used as an ionic tracer and to introduce the density contrast. Dispersivity in stable configuration increases with polymer concentration – as already reported for double porosity packings of porous grains. Gravity-induced instabilities are shown to develop below a same threshold Péclet number Pe for water and water–glycerol solutions of different viscosities and result in considerable increases of the dispersivity. Measured threshold Pe values decrease markedly on the contrary with polymer concentration. The quantitative analysis demonstrates that the development of the instabilities is controlled by viscosity through a characteristic gravity number G (ratio between hydrostatic and viscous pressure gradients). A single threshold value of G accounts for results obtained on Newtonian and non-Newtonian solutions.     

Three-sphere low-Reynolds-number swimmer with a cargo container

A recently introduced model for an autonomous swimmer at low Reynolds number that is comprised of three spheres connected by two arms is considered when one of the spheres has a large radius. The Stokes hydrodynamic flow associated with the swimming strokes and net motion of this system can be studied analytically using the Stokes Green's function of a point force in front of a sphere of arbitrary radius R provided by Oseen. The swimming velocity is calculated, and shown to scale as 1/R ³ with the radius of the sphere.     

Radiation Emission of Fast Electrons in Collisions with “Ion‐Sphere” in Dense Plasmas

Radiative emission of fast electrons in collision with an “ion‐sphere” electron distribution in dense plasmas is under consideration. The electron structure of the ion sphere is calculated ab initio using self‐consistent solution of both bound and free electron distribution inside the sphere. Two radiation channels are included: emission of the colliding electron itself in static potential (conventional or static Bremsstrahlung) and the emission of “ion sphere” medium due to its polarization by the colliding electron (polarization Bremsstrahlung). The last one is calculated in the frame of local plasma density approximation. Interference between conventional and polarization Bremsstrahlung is taken into account. It is shown that spectral cross section of the process has characteristic features depending on plasma density and ionization stage of plasma ions. (© 2013 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)