Malignant cancers that lead to fatal outcomes for patients may remain dormant for very long periods of time. Although individual mechanisms such as cellular dormancy, angiogenic dormancy and immunosurveillance have been proposed, a comprehensive understanding of cancer dormancy and the “switch” from a dormant to a proliferative state still needs to be strengthened from both a basic and clinical point of view. Computational modeling enables one to explore a variety of scenarios for possible but realistic microscopic dormancy mechanisms and their predicted outcomes. The aim of this paper is to devise such a predictive computational model of dormancy with an emergent “switch” behavior. Specifically, we generalize a previous cellular automaton (CA) model for proliferative growth of solid tumor that now incorporates a variety of cell-level tumor-host interactions and different mechanisms for tumor dormancy, for example the effects of the immune system. Our new CA rules induce a natural “competition” between the tumor and tumor suppression factors in the microenvironment. This competition either results in a “stalemate” for a period of time in which the tumor either eventually wins (spontaneously emerges) or is eradicated; or it leads to a situation in which the tumor is eradicated before such a “stalemate” could ever develop. We also predict that if the number of actively dividing cells within the proliferative rim of the tumor reaches a critical, yet low level, the dormant tumor has a high probability to resume rapid growth. Our findings may shed light on the fundamental understanding of cancer dormancy.

Numerous recent investigations have been devoted to the determination of the equilibrium phase behavior and packing characteristics of hard nonspherical particles, including ellipsoids, superballs, and polyhedra, to name but just a few shapes. 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, verifying the result of a previous qualitative study [J. Chem. Phys. 2011, 135, 151101]. The freezing- and melting-point packing fractions for this transition are respectively ϕ_F = 0.496 ± 0.006 and ϕ_M = 0.591 ± 0.005.

At higher densities, we find that the CT phase undergoes another first-order phase transition to one associated with the densest-known crystal, with coexistence densities in the range ϕ ∈ [0.780 ± 0.002, 0.802 ± 0.003]. 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. Specifically, we systematically study the structural characteristics of the MRJ packings, including the centroidal pair correlation function, structure factor and orientational pair correlation function. 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 (≈ 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.