Effect of Dimensionality on the Continuum Percolation of Overlapping Hyperspheres and Hypercubes: II. Simulation Results and Analyses

In the first paper of this series [S. Torquato, J. Chem. Phys. {\bf 136}, 054106 (2012)], analytical results concerning the continuum percolation of overlapping hyperparticles in $d$-dimensional Euclidean space $\mathbb{R}^d$ were obtained, including lower bounds on the percolation threshold. In the present investigation, we provide additional analytical results for certain cluster statistics, such as the concentration of $k$-mers and related quantities, and obtain an upper bound on the percolation threshold $\eta_c$. We utilize the tightest lower bound obtained in the first paper to formulate an efficient simulation method, called the {\it rescaled-particle} algorithm, to estimate continuum percolation properties across many space dimensions with heretofore unattained accuracy. This simulation procedure is applied to compute the threshold $\eta_c$ and associated mean number of overlaps per particle ${\cal N}_c$ for both overlapping hyperspheres and oriented hypercubes for $3 \le d \le 11$. These simulations results are compared to corresponding upper and lower bounds on these percolation properties. We find that the bounds converge to one another as the space dimension increases, but the lower bound provides an excellent estimate of $\eta_c$ and ${\cal N}_c$, even for relatively low dimensions. We confirm a prediction of the first paper in this series that low-dimensional percolation properties encode high-dimensional information. We also show that the concentration of monomers dominate over concentration values for higher-order clusters (dimers, trimers, etc.) as the space dimension becomes large. Finally, we provide accurate analytical estimates of the pair connectedness function and blocking function at their contact values for any $d$ as a function of density.Comment: 24 pages, 10 figure

Effect of Dimensionality on the Percolation Thresholds of Various dd-Dimensional Lattices

We show analytically that the $[0,1]$, $[1,1]$ and $[2,1]$ Pad{\'e} approximants of the mean cluster number $S(p)$ for site and bond percolation on general $d$-dimensional lattices are upper bounds on this quantity in any Euclidean dimension $d$, where $p$ is the occupation probability. These results lead to certain lower bounds on the percolation threshold $p_c$ that become progressively tighter as $d$ increases and asymptotically exact as $d$ becomes large. These lower-bound estimates depend on the structure of the $d$-dimensional lattice and whether site or bond percolation is being considered. We obtain explicit bounds on $p_c$ for both site and bond percolation on five different lattices: $d$-dimensional generalizations of the simple-cubic, body-centered-cubic and face-centered-cubic Bravais lattices as well as the $d$-dimensional generalizations of the diamond and kagom{\'e} (or pyrochlore) non-Bravais lattices. These analytical estimates are used to assess available simulation results across dimensions (up through $d=13$ in some cases). It is noteworthy that the tightest lower bound provides reasonable estimates of $p_c$ in relatively low dimensions and becomes increasingly accurate as $d$ grows. We also derive high-dimensional asymptotic expansions for $p_c$ for the ten percolation problems and compare them to the Bethe-lattice approximation. Finally, we remark on the radius of convergence of the series expansion of $S$ in powers of $p$ as the dimension grows.Comment: 37 pages, 5 figure

Diversity of Dynamics and Morphologies of Invasive Solid Tumors

Complex tumor-host interactions can significantly affect the growth dynamics and morphologies of progressing neoplasms. The growth of a confined solid tumor induces mechanical pressure and deformation of the surrounding microenvironment, which in turn influences tumor growth. In this paper, we generalize a recently developed cellular automaton model for invasive tumor growth in heterogeneous microenvironments [Y. Jiao and S. Torquato, PLoS Comput. Biol. 7, e1002314 (2011)] by incorporating the effects of pressure. Specifically, we explicitly model the pressure exerted on the growing tumor due to the deformation of the microenvironment and its effect on the local tumor-host interface instability. Both noninvasive-proliferative growth and invasive growth with individual cells that detach themselves from the primary tumor and migrate into the surrounding microenvironment are investigated. We find that while noninvasive tumors growing in "soft" homogeneous microenvironments develop almost isotropic shapes, both high pressure and host heterogeneity can strongly enhance malignant behavior, leading to finger-like protrusions of the tumor surface. Moreover, we show that individual invasive cells of an invasive tumor degrade the local extracellular matrix at the tumor-host interface, which diminishes the fingering growth of the primary tumor. The implications of our results for cancer diagnosis, prognosis and therapy are discussed.Comment: 21 pages, 5 figures, invited article for the special issue "Physics of Cancer" in AIP Advances, in pres

Accretion onto a Kiselev black hole

We consider accretion onto a Kiselev black hole. We obtain the fundamental equations for accretion without the back-reaction. We determine the general analytic expressions for the critical points and the mass accretion rate and find the physical conditions the critical points should fulfill. The case of polytropic gas are discussed in detail. It turns out that the quintessence parameter plays an important role in the accretion process.Comment: 7 page

Equilibrium Phase Behavior and Maximally Random Jammed State of Truncated Tetrahedra

Systems of hard nonspherical particles exhibit a variety of stable phases with different degrees of translational and orientational order, including isotropic liquid, solid crystal, rotator and a variety of liquid crystal phases. In this paper, we employ a Monte Carlo implementation of the adaptive-shrinking-cell (ASC) numerical scheme and free-energy calculations to ascertain with high precision the equilibrium phase behavior of systems of congruent Archimedean truncated tetrahedra over the entire range of possible densities up to the maximal nearly space-filling density. In particular, we find that the system undergoes two first-order phase transitions as the density increases: first a liquid-solid transition and then a solid-solid transition. The isotropic liquid phase coexists with the Conway-Torquato (CT) crystal phase at intermediate densities. At higher densities, we find that the CT phase undergoes another first-order phase transition to one associated with the densest-known crystal. We find no evidence for stable rotator (or plastic) or nematic phases. We also generate the maximally random jammed (MRJ) packings of truncated tetrahedra, which may be regarded to be the glassy end state of a rapid compression of the liquid. We find that such MRJ packings are hyperuniform with an average packing fraction of 0.770, which is considerably larger than the corresponding value for identical spheres (about 0.64). We conclude with some simple observations concerning what types of phase transitions might be expected in general hard-particle systems based on the particle shape and which would be good glass formers