Catalog footprint

What is connected

40works
24topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

40 published item(s)

preprint2022arXiv

Disordered Hyperuniform Quasi-1D Materials

Carbon nanotubes are quasi-one-dimensional systems that possess superior transport, mechanical, optical, and chemical properties. In this work, we generalize the notion of disorder hyperuniformity, a recently discovered exotic state of matter with hidden long-range order, to quasi-one-dimensional materials. As a proof of concept, we then apply the generalized framework to quantify the density fluctuations in amorphous carbon nanotubes containing randomly distributed Stone-Wales defects. We demonstrate that all of these amorphous nanotubes are hyperuniform, i.e., the infinite-wavelength density fluctuations of these systems are completely suppressed, regardless of the diameter, rolling axis, number of rolling sheets, and defect fraction of the nanotubes. We find that these amorphous nanotubes are energetically more stable than nanotubes with periodically distributed Stone-Wales defects. Moreover, certain semiconducting defect-free carbon nanotubes become metallic as sufficiently large amounts of defects are randomly introduced. This structural study of amorphous nanotubes strengthens our fundamental understanding of these systems, and suggests possible exotic physical properties, as endowed by their disordered hyperuniformity. Our findings also shed light on the effect of dimensionality reduction on the hyperuniformity property of materials.

preprint2022arXiv

MORE: Multi-Order RElation Mining for Dense Captioning in 3D Scenes

3D dense captioning is a recently-proposed novel task, where point clouds contain more geometric information than the 2D counterpart. However, it is also more challenging due to the higher complexity and wider variety of inter-object relations contained in point clouds. Existing methods only treat such relations as by-products of object feature learning in graphs without specifically encoding them, which leads to sub-optimal results. In this paper, aiming at improving 3D dense captioning via capturing and utilizing the complex relations in the 3D scene, we propose MORE, a Multi-Order RElation mining model, to support generating more descriptive and comprehensive captions. Technically, our MORE encodes object relations in a progressive manner since complex relations can be deduced from a limited number of basic ones. We first devise a novel Spatial Layout Graph Convolution (SLGC), which semantically encodes several first-order relations as edges of a graph constructed over 3D object proposals. Next, from the resulting graph, we further extract multiple triplets which encapsulate basic first-order relations as the basic unit, and construct several Object-centric Triplet Attention Graphs (OTAG) to infer multi-order relations for every target object. The updated node features from OTAG are aggregated and fed into the caption decoder to provide abundant relational cues, so that captions including diverse relations with context objects can be generated. Extensive experiments on the Scan2Cap dataset prove the effectiveness of our proposed MORE and its components, and we also outperform the current state-of-the-art method. Our code is available at https://github.com/SxJyJay/MORE.

preprint2022arXiv

Task-aware Similarity Learning for Event-triggered Time Series

Time series analysis has achieved great success in diverse applications such as network security, environmental monitoring, and medical informatics. Learning similarities among different time series is a crucial problem since it serves as the foundation for downstream analysis such as clustering and anomaly detection. It often remains unclear what kind of distance metric is suitable for similarity learning due to the complex temporal dynamics of the time series generated from event-triggered sensing, which is common in diverse applications, including automated driving, interactive healthcare, and smart home automation. The overarching goal of this paper is to develop an unsupervised learning framework that is capable of learning task-aware similarities among unlabeled event-triggered time series. From the machine learning vantage point, the proposed framework harnesses the power of both hierarchical multi-scale sequence autoencoders and Gaussian Mixture Model (GMM) to effectively learn the low-dimensional representations from the time series. Finally, the obtained similarity measure can be easily visualized for explaining. The proposed framework aspires to offer a stepping stone that gives rise to a systematic approach to model and learn similarities among a multitude of event-triggered time series. Through extensive qualitative and quantitative experiments, it is revealed that the proposed method outperforms state-of-the-art methods considerably.

preprint2022arXiv

vdW-DF-ahcx: a range-separated van der Waals density functional hybrid

Hybrid density functionals replace a fraction of an underlying generalized-gradient approximation (GGA) exchange description with a Fock-exchange component. Range-separated hybrids (RSHs) also effectively screen the Fock-exchange component and thus open the door for characterizations of metals and adsorption at metal surfaces. The RSHs are traditionally based on a robust GGA, such as PBE [PRL $\textbf{77}$, 3865 (1996)], for example, as implemented in the HSE design [JPC $\textbf{118}$, 8207 (2003)]. Here we define a RSH extension to the van der Waals density functional (vdW-DF) method [ROPP $\textbf{78}$, 066501 (2015)], launching vdW-DF-ahcx. We use an analytical-hole (AH) framework [JCP $\textbf{128}$, 194105 (2008)] to characterize the GGA-type exchange in the vdW-DF-cx version [PRB $\textbf{89}$, 075148 (2014)], isolate the short-ranged exchange component, and define the new RSH. We find that the performance vdW-DF-ahcx compares favorably to (dispersion-corrected) HSE for descriptions of bulk (broad molecular) properties. We also find that it provides accurate descriptions of noble-metal surface properties, including CO adsorption.

preprint2021arXiv

Optical Flow Estimation via Motion Feature Recovery

Optical flow estimation with occlusion or large displacement is a problematic challenge due to the lost of corresponding pixels between consecutive frames. In this paper, we discover that the lost information is related to a large quantity of motion features (more than 40%) computed from the popular discriminative cost-volume feature would completely vanish due to invalid sampling, leading to the low efficiency of optical flow learning. We call this phenomenon the Vanishing Cost Volume Problem. Inspired by the fact that local motion tends to be highly consistent within a short temporal window, we propose a novel iterative Motion Feature Recovery (MFR) method to address the vanishing cost volume via modeling motion consistency across multiple frames. In each MFR iteration, invalid entries from original motion features are first determined based on the current flow. Then, an efficient network is designed to adaptively learn the motion correlation to recover invalid features for lost-information restoration. The final optical flow is then decoded from the recovered motion features. Experimental results on Sintel and KITTI show that our method achieves state-of-the-art performances. In fact, MFR currently ranks second on Sintel public website.

preprint2021arXiv

Quantifying Microstructural Evolution via Time-Dependent Reduced-Dimension Metrics Based on Hierarchical $n$-Point Polytope Functions

We devise reduced-dimension metrics for effectively measuring the distance between two points (i.e., microstructures) in the microstructure space and quantifying the pathway associated with microstructural evolution, based on a recently introduced set of hierarchical $n$-point polytope functions $P_n$. The $P_n$ functions provide the probability of finding particular $n$-point configurations associated with regular $n$-polytopes in the material system, and a special sub-set of the standard $n$-point correlation functions $S_n$ that effectively decomposes the structural features in the systems into regular polyhedral basis with different symmetry. The $n$-th order metric $Ω_n$ is defined as the $\mathbb{L}_1$ norm associated with the $P_n$ functions of two distinct microstructures. By choosing a reference initial state (i.e., a microstructure associated with $t_0 = 0$), the $Ω_n(t)$ set quantifies the evolution of distinct polyhedral symmetries and can in principle capture emerging polyhedral symmetries that are not apparent in the initial state. To demonstrate their utility, we apply the $Ω_n$ metrics to a 2D binary system undergoing spinodal decomposition to extract the phase separation dynamics via the temporal scaling behavior of the corresponding $Ω_n(t)$, which reveals mechanisms governing the evolution. Moreover, we employ $Ω_n(t)$ to analyze pattern evolution during vapor-deposition of phase-separating alloy films with different surface contact angles, which exhibit rich evolution dynamics including both unstable and oscillating patterns. The $Ω_n$ metrics have potential applications in establishing quantitative processing-structure-property relationships, as well as real-time processing control and optimization of complex heterogeneous material systems.

preprint2020arXiv

Information Content of Hierarchical n-Point Polytope Functions for Quantifying and Reconstructing Disordered Systems

Disordered systems are ubiquitous in physical, biological and material sciences. Examples include liquid and glassy states of condensed matter, colloids, granular materials, porous media, composites, alloys, packings of cells in avian retina and tumor spheroids, to name but a few. A comprehensive understanding of such disordered systems requires, as the first step, systematic quantification, modeling and representation of the underlying complex configurations and microstructure, which is generally very challenging to achieve. Recently, we introduce a set of hierarchical statistical microstructural descriptors, i.e., the n-point polytope functions Pn, which are derived from the standard n-point correlation functions Sn, and successively include higher-order n-point statistics of the morphological features of interest in a concise, explainable, and expressive manner. Here we investigate the information content of the Pn functions via optimization-based realization rendering. This is achieved by successively incorporating higher order Pn functions up to n = 8 and quantitatively assessing the accuracy of the reconstructed systems via un-constrained statistical morphological descriptors (e.g., the lineal-path function). We examine a wide spectrum of representative random systems with distinct geometrical and topological features. We find that generally, successively incorporating higher order Pn functions, and thus, the higher-order morphological information encoded in these descriptors, leads to superior accuracy of the reconstructions. However, incorporating more Pn functions into the reconstruction also significantly increases the complexity and roughness of the associated energy landscape for the underlying stochastic optimization, making it difficult to convergence numerically.

preprint2020arXiv

Microstructure Design of Low-Melting-Point Alloy (LMPA)/ Polymer Composites for Dynamic Dry Adhesion Tuning in Soft Gripping

Tunable dry adhesion is a crucial mechanism in compliant manipulation. The gripping force, mainly originated from the van der Waals force between the adhesive composite and the object to be gripped, can be controlled by reversibly varying the physical properties (e.g., stiffness) of the composite via external stimuli. The maximal gripping force Fmax and its tunability depend on, among other factors, the stress distribution on the gripping interface and its fracture dynamics (during detaching), which in turn are determined by the composite microstructure. Here, we present a computational framework for the modeling and design of a class of binary smart composites containing a porous low-melting-point alloy (LMPA) phase and a polymer phase, in order to achieve desirable dynamically tunable dry adhesion. In particular, we employ spatial correlation functions to quantify, model and represent the complex bi-continuous microstructure of the composites, from which a wide spectrum of realistic virtual 3D composite microstructures can be generated using stochastic optimization. A recently developed volume-compensated lattice-particle (VCLP) method is then employed to model the dynamic interfacial fracture process to compute Fmax for different composite microstructures. We focus on the interface defect tuning (IDT) mechanism for dry adhesion tuning enabled by the composite, in which the thermal expansion of the LMPA phase due to Joule heating initializes small cracks on the adhesion interface, subsequently causing the detachment of the gripper from the object due to interfacial fracture. We find that for an optimal microstructure among the ones studied here, a 10-fold dynamic tuning of Fmax before and after the thermal expansion of the LMPA phase can be achieved. Our computational results can provide valuable guidance for experimental fabrication of the LMPA-polymer composites.

preprint2020arXiv

Optical nonreciprocal response and conversion in a Tavis-Cummings coupling optomechanical system

We propose a scheme to realize optical nonreciprocal response and conversion in a Tavis-Cummings coupling optomechanical system, where a single cavity mode interacts with the vibrational mode of a flexible membrane with an embedded ensemble of two-level quantum emitters. Due to the introduction of the Tavis-Cummings interaction, we find that the phases between the mechanical mode and the optical mode, as well as between the mechanical mode and the dopant mode, are correlated with each other, and further give the analytical relationship between them. By optimizing the system parameters, especially the relative phase between two paths, the optimal nonreciprocal response can be achieved. Under the frequency domain, we derive the transmission matrix of the system analytically based on the input-output relation and study the influence of the system parameters on the nonreciprocal response of the quantum input signal. Moreover, compared with the conventional optomechanical systems, the Tavis-Cummings coupling optomechanical system exhibits richer nonreciprocal conversion phenomena among the optical mode, mechanical mode, and dopant mode, which provide a new applicable way of achieving the phonon-photon transducer and the optomechanical circulator in future practice.

preprint2020arXiv

Stone-Wales Defects Preserve Hyperuniformity in Amorphous Two-Dimensional Materials

Crystalline two-dimensional (2D) materials such as graphene possess unique physical properties absent in their bulk form, enabling many novel device applications. Yet, little is known about their amorphous counterparts, which can be obtained by introducing the Stone-Wales (SW) topological defects via proton radiation. Here we provide strong numerical evidence that SW defects preserve hyperuniformity in hexagonal 2D materials, a recently discovered new state of matter characterized by vanishing normalized infinite-wavelength density fluctuations, which implies that all amorphous states of these materials are hyperuniform. Specifically, the static structure factor S(k) of these materials possesses the scaling S(k) ~ k^α for small wave number k, where 1<=α(p)<=2 is monotonically decreasing as the SW defect concentration p increases, indicating a transition from type-I to type-II hyperuniformity at p ~= 0.12 induced by the saturation of the SW defects. This hyperuniformity transition marks a structural transition from perturbed lattice structures to truly amorphous structures, and underlies the onset of strong correlation among the SW defects as well as a transition between distinct electronic transport mechanisms associated with different hyperuniformity classes.

preprint2020arXiv

Translate Reverberated Speech to Anechoic Ones: Speech Dereverberation with BERT

Single channel speech dereverberation is considered in this work. Inspired by the recent success of Bidirectional Encoder Representations from Transformers (BERT) model in the domain of Natural Language Processing (NLP), we investigate its applicability as backbone sequence model to enhance reverberated speech signal. We present a variation of the basic BERT model: a pre-sequence network, which extracts local spectral-temporal information and/or provides order information, before the backbone sequence model. In addition, we use pre-trained neural vocoder for implicit phase reconstruction. To evaluate our method, we used the data from the 3rd CHiME challenge, and compare our results with other methods. Experiments show that the proposed method outperforms traditional method WPE, and achieve comparable performance with state-of-the-art BLSTM-based sequence models.

preprint2018arXiv

Modeling Three-dimensional Invasive Solid Tumor Growth in Heterogeneous Microenvironment under Chemotherapy

A systematic understanding of the evolution and growth dynamics of invasive solid tumors in response to different chemotherapy strategies is crucial for the development of individually optimized oncotherapy. Here, we develop a hybrid three-dimensional (3D) computational model that integrates pharmacokinetic model, continuum diffusion-reaction model and discrete cell automaton model to investigate 3D invasive solid tumor growth in heterogeneous microenvironment under chemotherapy. Specifically, we consider the effects of heterogeneous environment on drug diffusion, tumor growth, invasion and the drug-tumor interaction on individual cell level. We employ the hybrid model to investigate the evolution and growth dynamics of avascular invasive solid tumors under different chemotherapy strategies. Our simulations reproduce the well-established observation that constant dosing is generally more effective in suppressing primary tumor growth than periodic dosing, due to the resulting continuous high drug concentration. In highly heterogeneous microenvironment, the malignancy of the tumor is significantly enhanced, leading to inefficiency of chemotherapies. The effects of geometrically-confined microenvironment and non-uniform drug dosing are also investigated. Our computational model, when supplemented with sufficient clinical data, could eventually lead to the development of efficient in silico tools for prognosis and treatment strategy optimization.

preprint2016arXiv

Algorithms for Automatic Ranking of Participants and Tasks in an Anonymized Contest

We introduce a new set of problems based on the Chain Editing problem. In our version of Chain Editing, we are given a set of anonymous participants and a set of undisclosed tasks that every participant attempts. For each participant-task pair, we know whether the participant has succeeded at the task or not. We assume that participants vary in their ability to solve tasks, and that tasks vary in their difficulty to be solved. In an ideal world, stronger participants should succeed at a superset of tasks that weaker participants succeed at. Similarly, easier tasks should be completed successfully by a superset of participants who succeed at harder tasks. In reality, it can happen that a stronger participant fails at a task that a weaker participants succeeds at. Our goal is to find a perfect nesting of the participant-task relations by flipping a minimum number of participant-task relations, implying such a "nearest perfect ordering" to be the one that is closest to the truth of participant strengths and task difficulties. Many variants of the problem are known to be NP-hard. We propose six natural $k$-near versions of the Chain Editing problem and classify their complexity. The input to a $k$-near Chain Editing problem includes an initial ordering of the participants (or tasks) that we are required to respect by moving each participant (or task) at most $k$ positions from the initial ordering. We obtain surprising results on the complexity of the six $k$-near problems: Five of the problems are polynomial-time solvable using dynamic programming, but one of them is NP-hard.

preprint2015arXiv

Heterogeneous Force Chains in Cellularized Biopolymer Network

Biopolymer Networks play an important role in coordinating and regulating collective cellular dynamics via a number of signaling pathways. Here, we investigate the mechanical response of a model biopolymer network due to the active contraction of embedded cells. Specifically, a graph (bond-node) model derived from confocal microscopy data is used to represent the network microstructure, and cell contraction is modeled by applying correlated displacements at specific nodes, representing the focal adhesion sites. A force-based stochastic relaxation method is employed to obtain force-balanced network under cell contraction. We find that the majority of the forces are carried by a small number of heterogeneous force chains emitted from the contracting cells. The force chains consist of fiber segments that either possess a high degree of alignment before cell contraction or are aligned due to the reorientation induced by cell contraction. Large fluctuations of the forces along different force chains are observed. Importantly, the decay of the forces along the force chains is significantly slower than the decay of radially averaged forces in the system. These results suggest that the fibrous nature of biopolymer network structure can support long-range force transmission and thus, long-range mechanical signaling between cells.

preprint2015arXiv

Schottky barrier formation and band bending revealed by first principles calculations

An atomistic insight into potential barrier formation and band bending at the interface between a metal and an n-type semiconductor is achieved by ab initio simulations and model analysis of a prototype Schottky diode, i.e., niobium doped rutile titania in contact with gold (Au/Nb:TiO$_2$). The local Schottky barrier height is found to vary between 0 and 1.26 eV depending on the position of the dopant. The band bending is caused by a dopant induced dipole field between the interface and the dopant site, whereas the pristine Au/TiO$_2$ interface does not show any band bending. These findings open the possibility for atomic scale optimization of the Schottky barrier and light harvesting in metal-semiconductor nanostructures.

preprint2014arXiv

A Cellular Automaton Model for Tumor Dormancy: Emergence of a Proliferative Switch

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.

preprint2014arXiv

Avian photoreceptor patterns represent a disordered hyperuniform solution to a multiscale packing problem

Optimal spatial sampling of light rigorously requires that identical photoreceptors be arranged in perfectly regular arrays in two dimensions. Examples of such perfect arrays in nature include the compound eyes of insects and the nearly crystalline photoreceptor patterns of some fish and reptiles. Birds are highly visual animals with five different cone photoreceptor subtypes, yet their photoreceptor patterns are not perfectly regular. By analyzing the chicken cone photoreceptor system consisting of five different cell types using a variety of sensitive microstructural descriptors, we find that the disordered photoreceptor patterns are ``hyperuniform'' (exhibiting vanishing infinite-wavelength density fluctuations), a property that had heretofore been identified in a unique subset of physical systems, but had never been observed in any living organism. Remarkably, the photoreceptor patterns of both the total population and the individual cell types are simultaneously hyperuniform. We term such patterns ``multi-hyperuniform'' because multiple distinct subsets of the overall point pattern are themselves hyperuniform. We have devised a unique multiscale cell packing model in two dimensions that suggests that photoreceptor types interact with both short- and long-ranged repulsive forces and that the resultant competition between the types gives rise to the aforementioned singular spatial features characterizing the system, including multi-hyperuniformity.

preprint2014arXiv

Dense periodic packings of tori

Dense packings of nonoverlapping bodies in three-dimensional Euclidean space are useful models of the structure of a variety of many-particle systems that arise in the physical and biological sciences. Here we investigate the packing behavior of congruent ring tori, which are multiply connected nonconvex bodies of genus 1, as well as horn and spindle tori. We analytically construct a family of dense periodic packings of unlinked tori guided by the organizing principles originally devised for simply connected solid bodies [Torquato and Jiao, PRE 86, 011102 (2012)]. We find that the horn tori as well as certain spindle and ring tori can achieve a packing density higher than the densest known packing of both sphere and ellipsoids. In addition, we study dense packings of cluster of pair-linked ring tori (i.e., Hopf links).

preprint2014arXiv

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.

preprint2014arXiv

Maximally dense packings of two-dimensional convex and concave noncircular particles

Dense packings of hard particles have important applications in many fields, including condensed matter physics, discrete geometry and cell biology. In this paper, we employ a stochastic search implementation of the Torquato-Jiao Adaptive-Shrinking-Cell optimization scheme [Nature 460, 876 (2009)] to find maximally dense particle packings in d-dimensional Euclidean space $\mathbb{R}^d$. While the original implementation was designed to study spheres and convex polyhedra in $d\ge3$, our implementation focuses on $d=2$ and extends the algorithm to include both concave polygons and certain convex or concave non-polygonal particle shapes. We verify the robustness of this packing protocol by successfully reproducing the known putative optimal packings of congruent copies of regular pentagons and octagons, then employ it to suggest dense packing arrangements of congruent copies of certain families of concave crosses, convex and concave curved triangles (incorporating shapes resembling the Mercedes-Benz logo), and "moon-like" shapes. Analytical constructions are determined subsequently to obtain the densest known packings of these particle shapes. For the examples considered, we find that the densest packings of both convex and concave particles with central symmetry are achieved by their corresponding optimal Bravais lattice packings; for particles lacking central symmetry, the densest packings obtained are non-lattice periodic packings, which are consistent with recently-proposed general organizing principles for hard particles. Moreover, we find that the densest known packings of certain curved triangles are periodic with a four-particle basis, and we find that the densest known periodic packings of certain moon-like shapes possess no inherent symmetries. Our work adds to the growing evidence that particle shape can be used as a tuning parameter to achieve a diversity of packing structures.

preprint2013arXiv

Effect of Dimensionality on the Percolation Threshold of Overlapping Nonspherical Hyperparticles

A set of lower bounds on the continuum percolation threshold $η_c$ of overlapping convex hyperparticles of general nonspherical (anisotropic) shape with a specified orientational probability distribution in $d$-dimensional Euclidean space have been derived [S. Torquato, J. Chem. Phys. {\bf 136}, 054106 (2012)]. The simplest of these lower bounds is given by $η_c \ge v/v_{ex}$, where $v_{ex}$ is the $d$-dimensional exclusion volume of a hyperparticle and $v$ is its $d$-dimensional volume. In order to study the effect of dimensionality on the threshold $η_c$ of overlapping nonspherical convex hyperparticles with random orientations here, we obtain a scaling relation for $η_c$ that is based on this lower bound and a conjecture that hyperspheres provide the highest threshold among all convex hyperparticle shapes for any $d$. This scaling relation exploits the principle that low-dimensional continuum percolation behavior encodes high-dimensional information. We derive a formula for the exclusion volume $v_{ex}$ of a hyperparticle in terms of its $d$-dimensional volume $v$, surface area $s$ and {\it radius of mean curvature} ${\bar R}$ (or, equivalently, {\it mean width}). These basic geometrical properties are computed for a wide variety of nonspherical hyperparticle shapes with random orientations across all dimensions, including, among other shapes, various polygons for $d=2$, Platonic solids, spherocylinders, parallepipeds and zero-volume plates for $d=3$ and their appropriate generalizations for $d \ge 4$. We then compute the lower bound and scaling relation for $η_c$ for this comprehensive set of continuum percolation models across dimensions. We show that the scaling relation provides accurate {\it upper-bound} estimates of the threshold $η_c$ across dimensions and becomes increasingly accurate as the space $d$ increases.

preprint2013arXiv

Effect of Dimensionality on the Percolation Thresholds of Various $d$-Dimensional Lattices

We show analytically that the $[0,1]$, $[1,1]$ and $[2,1]$ Pad{é} 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{é} (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.

preprint2012arXiv

A Superior Descriptor of Random Textures and Its Predictive Capacity

Two-phase random textures abound in a host of contexts, porous and composite media, ecological structures, biological media and astrophysical structures. Questions surrounding the spatial structure of such textures continue to pose many theoretical challenges. For example, can two-point correlation functions be identified that can be both manageably measured and yet reflect nontrivial higher-order structural information about the textures? We present a novel solution to this question by probing the information content of the widest class of different types of two-point functions examined to date using inverse "reconstruction" techniques. This enables us to show that a superior descriptor is the two-point cluster function $C_2({\bf r})$, which is sensitive to topological {\it connectedness} information. We demonstrate the utility of $C_2({\bf r})$ by accurately reconstructing textures drawn from materials science, cosmology and granular media, among other examples. Our work suggests an entirely new theoretical pathway to predict the bulk physical properties of random textures, and also has important ramifications for atomic and molecular systems.

preprint2012arXiv

Density of States for a Specified Correlation Function and the Energy Landscape

The degeneracy of two-phase disordered microstructures consistent with a specified correlation function is analyzed by mapping it to a ground-state degeneracy. We determine for the first time the associated density of states via a Monte Carlo algorithm. Our results are described in terms of the roughness of the energy landscape, defined on a hypercubic configuration space. The use of a Hamming distance in this space enables us to define a roughness metric, which is calculated from the correlation function alone and related quantitatively to the structural degeneracy. This relation is validated for a wide variety of disordered systems.

preprint2012arXiv

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.

preprint2012arXiv

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 $η_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 $η_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 $η_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.

preprint2012arXiv

Microstructural Degeneracy associated with a Two-Point Correlation Function and its Information Content

Two-point correlation functions provide crucial yet incomplete characterization of microstructures because different microstructures may have the same correlation function. In an earlier Letter [Phys. Rev. Lett. 108, 080601 (2012)], we addressed the degeneracy question: What is the number of microstructures compatible with a specified correlation function? We computed this degeneracy, i.e., configurational entropy, in the framework of reconstruction methods, which enabled us to map the problem to the determination of ground-state degeneracies. Here, we provide a more comprehensive presentation and additional results. Since the configuration space of a reconstruction problem is a hypercube on which a Hamming distance is defined, we can calculate analytically an energy profile corresponding to the average energy of all microstructures at a given Hamming distance from a ground state. The steepness of this profile is a measure of the roughness of the energy landscape, which can be used as a proxy for ground-state degeneracy. The relationship between roughness metric and ground-state degeneracy is calibrated using a Monte Carlo algorithm for determining the degeneracy of a variety of microstructures, including hard disks and Poisson point processes as well as those with known degeneracies (single disks of various sizes and a particular crystalline microstructure). We show that our results can be expressed in terms of the information content of the two-point correlation functions. From this perspective, the a priori condition for a reconstruction to be accurate is that the information content, expressed in bits, should be comparable to the number of pixels in the unknown microstructure. We provide a formula to calculate the information content of any two-point correlation function, which makes our results broadly applicable to any field in which correlation functions are employed.

preprint2012arXiv

Organizing Principles for Dense Packings of Nonspherical Hard Particles: Not All Shapes Are Created Equal

We have recently devised organizing principles to obtain maximally dense packings of the Platonic and Archimedean solids, and certain smoothly-shaped convex nonspherical particles [Torquato and Jiao, Phys. Rev. E 81, 041310 (2010)]. Here we generalize them in order to guide one to ascertain the densest packings of other convex nonspherical particles as well as concave shapes. Our generalized organizing principles are explicitly stated as four distinct propositions. We apply and test all of these organizing principles to the most comprehensive set of both convex and concave particle shapes to date, including Catalan solids, prisms, antiprisms, cylinders, dimers of spheres and various concave polyhedra. We demonstrate that all of the densest known packings associated with this wide spectrum of nonspherical particles are consistent with our propositions. Among other applications, our general organizing principles enable us to construct analytically the densest known packings of certain convex nonspherical particles, including spherocylinders, "lens-shaped" particles, square pyramids and rhombic pyramids. Moreover, we show how to apply these principles to infer the high-density equilibrium crystalline phases of hard convex and concave particles. We also discuss the unique packing attributes of maximally random jammed packings of nonspherical particles.

preprint2012arXiv

Quantitative Characterization of the Microstructure and Transport Properties of Biopolymer Networks

Biopolymer networks are of fundamental importance to many biological processes in normal and tumorous tissues. In this paper, we employ the panoply of theoretical and simulation techniques developed for characterizing heterogeneous materials to quantify the microstructure and effective diffusive transport properties (diffusion coefficient $D_e$ and mean survival time $τ$) of collagen type I networks at various collagen concentrations. In particular, we compute the pore-size probability density function $P(δ)$ for the networks and present a variety of analytical estimates of the effective diffusion coefficient $D_e$ for finite-sized diffusing particles. The Hashin-Strikman upper bound on the effective diffusion coefficient $D_e$ and the pore-size lower bound on the mean survival time $τ$ are used as benchmarks to test our analytical approximations and numerical results. Moreover, we generalize the efficient first-passage-time techniques for Brownian-motion simulations in suspensions of spheres to the case of fiber networks and compute the associated effective diffusion coefficient $D_e$ as well as the mean survival time $τ$, which is related to nuclear magnetic resonance (NMR) relaxation times. Specifically, the Torquato approximation provides the most accurate estimates of $D_e$ for all collagen concentrations among all of the analytical approximations we consider. We formulate a universal curve for $τ$ for the networks at different collagen concentrations. We apply rigorous cross-property relations to estimate the effective bulk modulus of collagen networks from a knowledge of the effective diffusion coefficient computed here.

preprint2012arXiv

Rigidity of spherical codes

A packing of spherical caps on the surface of a sphere (that is, a spherical code) is called rigid or jammed if it is isolated within the space of packings. In other words, aside from applying a global isometry, the packing cannot be deformed. In this paper, we systematically study the rigidity of spherical codes, particularly kissing configurations. One surprise is that the kissing configuration of the Coxeter-Todd lattice is not jammed, despite being locally jammed (each individual cap is held in place if its neighbors are fixed); in this respect, the Coxeter-Todd lattice is analogous to the face-centered cubic lattice in three dimensions. By contrast, we find that many other packings have jammed kissing configurations, including the Barnes-Wall lattice and all of the best kissing configurations known in four through twelve dimensions. Jamming seems to become much less common for large kissing configurations in higher dimensions, and in particular it fails for the best kissing configurations known in 25 through 31 dimensions. Motivated by this phenomenon, we find new kissing configurations in these dimensions, which improve on the records set in 1982 by the laminated lattices.

preprint2011arXiv

Analytical Construction of A Dense Packing of Truncated Tetrahedra

Dense polyhedron packings are useful models of a variety of condensed matter and biological systems and have intrigued scientists mathematicians for centuries. Recently, organizing principles for the types of structures associated with the densest polyhedron packings have been put forth. However, finding the maximally dense packings for specific shapes remains nontrivial and is a subject of intensive research. Here, we analytically construct the densest known packing of truncated tetrahedra with $ϕ=207/208=0.995~192...$, which is amazingly close to unity and strongly implies that this packing is maximally dense. This construction is based on a generalized organizing principle for polyhedra that lack central symmetry. Moreover, we find that the holes in this putative optimal packing are small regular tetrahedra, leading to a new tiling of space by regular tetrahedra and truncated tetrahedra. We also numerically study the equilibrium melting properties of what apparently is the densest packing of truncated tetrahedra as the system undergoes decompression. Our simulations reveal two different stable crystal phases, one at high densities and the other at intermediate densities, as well as a first-order liquid-crystal phase transition.

preprint2011arXiv

Hyperuniform long-range correlations are a signature of disordered jammed hard-particle packings

We show that quasi-long-range (QLR) pair correlations that decay asymptotically with scaling $r^{-(d+1)}$ in $d$-dimensional Euclidean space $\mathbb{R}^d$, trademarks of certain quantum systems and cosmological structures, are a universal signature of maximally random jammed (MRJ) hard-particle packings. We introduce a novel hyperuniformity descriptor in MRJ packings by studying local-volume-fraction fluctuations and show that infinite-wavelength fluctuations vanish even for packings with size- and shape-distributions. Special void statistics induce hyperuniformity and QLR pair correlations.

preprint2011arXiv

Hyperuniformity, quasi-long-range correlations, and void-space constraints in maximally random jammed particle packings. I. Polydisperse spheres

Hyperuniform many-particle distributions possess a local number variance that grows more slowly than the volume of an observation window, implying that the local density is effectively homogeneous beyond a few characteristic length scales. Previous work on maximally random strictly jammed sphere packings in three dimensions has shown that these systems are hyperuniform and possess unusual quasi-long-range pair correlations, resulting in anomalous logarithmic growth in the number variance. However, recent work on maximally random jammed sphere packings with a size distribution has suggested that such quasi-long-range correlations and hyperuniformity are not universal among jammed hard-particle systems. In this paper we show that such systems are indeed hyperuniform with signature quasi-long-range correlations by characterizing the more general local-volume-fraction fluctuations. We argue that the regularity of the void space induced by the constraints of saturation and strict jamming overcomes the local inhomogeneity of the disk centers to induce hyperuniformity in the medium with a linear small-wavenumber nonanalytic behavior in the spectral density, resulting in quasi-long-range spatial correlations. A numerical and analytical analysis of the pore-size distribution for a binary MRJ system in addition to a local characterization of the n-particle loops governing the void space surrounding the inclusions is presented in support of our argument. This paper is the first part of a series of two papers considering the relationships among hyperuniformity, jamming, and regularity of the void space in hard-particle packings.

preprint2011arXiv

Hyperuniformity, quasi-long-range correlations, and void-space constraints in maximally random jammed particle packings. II. Anisotropy in particle shape

We extend the results from the first part of this series of two papers by examining hyperuniformity in heterogeneous media composed of impenetrable anisotropic inclusions. Specifically, we consider maximally random jammed packings of hard ellipses and superdisks and show that these systems both possess vanishing infinite-wavelength local-volume-fraction fluctuations and quasi-long-range pair correlations. Our results suggest a strong generalization of a conjecture by Torquato and Stillinger [Phys. Rev. E. 68, 041113 (2003)], namely that all strictly jammed saturated packings of hard particles, including those with size- and shape-distributions, are hyperuniform with signature quasi-long-range correlations. We show that our arguments concerning the constrained distribution of the void space in MRJ packings directly extend to hard ellipse and superdisk packings, thereby providing a direct structural explanation for the appearance of hyperuniformity and quasi-long-range correlations in these systems. Additionally, we examine general heterogeneous media with anisotropic inclusions and show for the first time that one can decorate a periodic point pattern to obtain a hard-particle system that is not hyperuniform with respect to local-volume-fraction fluctuations. This apparent discrepancy can also be rationalized by appealing to the irregular distribution of the void space arising from the anisotropic shapes of the particles. Our work suggests the intriguing possibility that the MRJ states of hard particles share certain universal features independent of the local properties of the packings, including the packing fraction and average contact number per particle.

preprint2011arXiv

Non-Universality of Density and Disorder in Jammed Sphere Packings

We show for the first time that collectively jammed disordered packings of three-dimensional monodisperse frictionless hard spheres can be produced and tuned using a novel numerical protocol with packing density $ϕ$ as low as 0.6. This is well below the value of 0.64 associated with the maximally random jammed state and entirely unrelated to the ill-defined ``random loose packing'' state density. Specifically, collectively jammed packings are generated with a very narrow distribution centered at any density $ϕ$ over a wide density range $ϕ\in [0.6,~0.74048\ldots]$ with variable disorder. Our results support the view that there is no universal jamming point that is distinguishable based on the packing density and frequency of occurence. Our jammed packings are mapped onto a density-order-metric plane, which provides a broader characterization of packings than density alone. Other packing characteristics, such as the pair correlation function, average contact number and fraction of rattlers are quantified and discussed.

preprint2011arXiv

Phase diagram and structural diversity of the densest binary sphere packings

The densest binary sphere packings have historically been very difficult to determine. The only rigorously known packings in the alpha-x plane of sphere radius ratio alpha and relative concentration x are at the Kepler limit alpha = 1, where packings are monodisperse. Utilizing an implementation of the Torquato-Jiao sphere-packing algorithm [S. Torquato and Y. Jiao, Phys. Rev. E 82, 061302 (2010)], we present the most comprehensive determination to date of the phase diagram in (alpha,x) for the densest binary sphere packings. Unexpectedly, we find many distinct new densest packings.

preprint2010arXiv

Geometrical Ambiguity of Pair Statistics. II. Heterogeneous Media

In a previous paper [Jiao, Stillinger and Torquato, PRE 81, 011105 (2010)], we considered the geometrical ambiguity of pair statistics associated with point configurations. Here we focus on the analogous problem for heterogeneous media (materials). The complex structures of heterogeneous media are usually characterized via statistical descriptors, such as the $n$-point correlation function $S_n$. An intricate inverse problem of practical importance is to what extent a medium can be reconstructed from the two-point correlation function $S_2$ of a target medium. Recently, general claims of the uniqueness of reconstructions using $S_2$ have been made based on numerical studies. Here, we provide a systematic approach to characterize the geometrical ambiguity of $S_2$ for both continuous two-phase heterogeneous media and their digitized representations in a mathematically precise way. In particular, we derive the exact conditions for the case where two distinct media possess identical $S_2$. The degeneracy conditions are given in terms of integral and algebraic equations for continuous media and their digitized representations, respectively. By examining these equations and constructing their rigorous solutions for specific examples, we conclusively show that in general $S_2$ is indeed not sufficient information to uniquely determine the structure of the medium, which is consistent with the results of our recent study on heterogeneous media reconstruction [Jiao, Stillinger and Torquato, PNAS 106, 17634 (2009)]. The uniqueness issue of multiphase media reconstructions and additional structural information required to characterize heterogeneous media are discussed.

preprint2010arXiv

Robust Algorithm to Generate a Diverse Class of Dense Disordered and Ordered Sphere Packings via Linear Programming

We have formulated the problem of generating periodic dense paritcle packings as an optimization problem called the Adaptive Shrinking Cell (ASC) formulation [S. Torquato and Y. Jiao, Phys. Rev. E {\bf 80}, 041104 (2009)]. Because the objective function and impenetrability constraints can be exactly linearized for sphere packings with a size distribution in $d$-dimensional Euclidean space $\mathbb{R}^d$, it is most suitable and natural to solve the corresponding ASC optimization problem using sequential linear programming (SLP) techniques. We implement an SLP solution to produce robustly a wide spectrum of jammed sphere packings in $\mathbb{R}^d$ for $d=2,3,4,5$ and $6$ with a diversity of disorder and densities up to the maximally densities. This deterministic algorithm can produce a broad range of inherent structures besides the usual disordered ones with very small computational cost by tuning the radius of the {\it influence sphere}. In three dimensions, we show that it can produce with high probability a variety of strictly jammed packings with a packing density anywhere in the wide range $[0.6, 0.7408...]$. We also apply the algorithm to generate various disordered packings as well as the maximally dense packings for $d=2,3, 4,5$ and 6. Compared to the LS procedure, our SLP protocol is able to ensure that the final packings are truly jammed, produces disordered jammed packings with anomalously low densities, and is appreciably more robust and computationally faster at generating maximally dense packings, especially as the space dimension increases.

preprint2009arXiv

Optimal Packings of Superballs

Dense hard-particle packings are intimately related to the structure of low-temperature phases of matter and are useful models of heterogeneous materials and granular media. Most studies of the densest packings in three dimensions have considered spherical shapes, and it is only more recently that nonspherical shapes (e.g., ellipsoids) have been investigated. Superballs (whose shapes are defined by |x1|^2p + |x2|^2p + |x3|^2p <= 1) provide a versatile family of convex particles (p >= 0.5) with both cubic- and octahedral-like shapes as well as concave particles (0 < p < 0.5) with octahedral-like shapes. In this paper, we provide analytical constructions for the densest known superball packings for all convex and concave cases. The candidate maximally dense packings are certain families of Bravais lattice packings. The maximal packing density as a function of p is nonanalytic at the sphere-point (p = 1) and increases dramatically as p moves away from unity. The packing characteristics determined by the broken rotational symmetry of superballs are similar to but richer than their two-dimensional "superdisk" counterparts, and are distinctly different from that of ellipsoid packings. Our candidate optimal superball packings provide a starting point to quantify the equilibrium phase behavior of superball systems, which should deepen our understanding of the statistical thermodynamics of nonspherical-particle systems.