-
Effective medium properties of stealthy hyperuniform photonic structures using multiscale physics-informed neural networks
Authors:
Roberto Riganti,
Yilin Zhu,
Wei Cai,
Salvatore Torquato,
Luca Dal Negro
Abstract:
In this article, we employ multiscale physics-informed neural networks (MscalePINNs) for the inverse retrieval of the effective permittivity and homogenization of finite-size photonic media with stealthy hyperuniform (SHU) disordered geometries. Specifically, we show that MscalePINNs are capable of capturing the fast spatial variations of complex fields scattered by arrays of dielectric nanocylind…
▽ More
In this article, we employ multiscale physics-informed neural networks (MscalePINNs) for the inverse retrieval of the effective permittivity and homogenization of finite-size photonic media with stealthy hyperuniform (SHU) disordered geometries. Specifically, we show that MscalePINNs are capable of capturing the fast spatial variations of complex fields scattered by arrays of dielectric nanocylinders arranged according to isotropic SHU point patterns, thus enabling a systematic methodology to inverse retrieve their effective dielectric profiles. Our approach extends the recently developed high-frequency homogenization theory of hyperuniform media and retrieves more general permittivity profiles for applications-relevant finite-size SHU systems, unveiling unique features related to their isotropic nature. In particular, we demonstrate the existence of a transparency region beyond the long-wavelength approximation, enabling effective and isotropic homogenization even without disorder-averaging, in contrast to the case of uncorrelated Poisson random patterns. We believe that the multiscale network approach introduced here enables the efficient inverse design of general effective media and finite-size metamaterials with isotropic electromagnetic responses beyond the limitations of traditional homogenization theories.
△ Less
Submitted 14 May, 2024; v1 submitted 13 May, 2024;
originally announced May 2024.
-
Microstructural and Transport Characteristics of Triply Periodic Bicontinuous Materials
Authors:
Salvatore Torquato,
Jaeuk Kim
Abstract:
3D bicontinuous two-phase materials are increasingly gaining interest because of their unique multifunctional characteristics and advancements in techniques to fabricate them. Due to their complex topological and structural properties, it still has been nontrivial to develop explicit microstructure-dependent formulas to predict accurately their physical properties. A primary goal of the present pa…
▽ More
3D bicontinuous two-phase materials are increasingly gaining interest because of their unique multifunctional characteristics and advancements in techniques to fabricate them. Due to their complex topological and structural properties, it still has been nontrivial to develop explicit microstructure-dependent formulas to predict accurately their physical properties. A primary goal of the present paper is to ascertain various microstructural and transport characteristics of five different models of triply periodic bicontinuous porous materials at a porosity $φ_1=1/2$: those in which the two-phase interfaces are the Schwarz P, Schwarz D and Schoen G minimal surfaces as well as two different pore-channel structures. We ascertain their spectral densities, pore-size distribution functions, local volume-fraction variances, and hyperuniformity order metrics and then use this information to estimate certain effective transport properties via closed-form microstructure-property formulas. Specifically, we estimate the recently introduced time-dependent diffusion spreadability exactly from the spectral density. Moreover, we accurately estimate the fluid permeability of such porous materials from the second moment of the pore-size function and the formation factor, a measure of the tortuosity of the pore space. We also rigorously bound the permeability from above using the spectral density. For the five models with identical cubic unit cells, we find that the permeability, inverse of the specific surface, hyperuniformity order metric, pore-size second moment and long-time spreadability behavior are all positively correlated and rank order the structures in exactly the same way. We also conjecture what structures maximize the fluid permeability for arbitrary porosities and show that this conjecture must be true in the extreme porosity limits by identifying the corresponding optimal structures.
△ Less
Submitted 8 July, 2024; v1 submitted 25 April, 2024;
originally announced May 2024.
-
Stealthy and hyperuniform isotropic photonic bandgap structure in 3D
Authors:
Lukas Siedentop,
Gianluc Lui,
Georg Maret,
Paul M. Chaikin,
Paul J. Steinhardt,
Salvatore Torquato,
Peter Keim,
Marian Florescu
Abstract:
In photonic crystals the propagation of light is governed by their photonic band structure, an ensemble of propagating states grouped into bands, separated by photonic band gaps. Due to discrete symmetries in spatially strictly periodic dielectric structures their photonic band structure is intrinsically anisotropic. However, for many applications, such as manufacturing artificial structural color…
▽ More
In photonic crystals the propagation of light is governed by their photonic band structure, an ensemble of propagating states grouped into bands, separated by photonic band gaps. Due to discrete symmetries in spatially strictly periodic dielectric structures their photonic band structure is intrinsically anisotropic. However, for many applications, such as manufacturing artificial structural color materials or developing photonic computing devices, but also for the fundamental understanding of light-matter interactions, it is of major interest to seek materials with long range non-periodic dielectric structures which allow the formation of {\it isotropic} photonic band gaps. Here, we report the first ever 3D isotropic photonic band gap for an optimized disordered stealthy hyperuniform structure for microwaves. The transmission spectra are directly compared to a diamond pattern and an amorphous structure with similar node density. The band structure is measured experimentally for all three microwave structures, manufactured by 3D-Laser-printing for meta-materials with refractive index up to $n=2.1$. Results agree well with finite-difference-time-domain numerical investigations and a priori calculations of the band-gap for the hyperuniform structure: the diamond structure shows gaps but being anisotropic as expected, the stealthy hyperuniform pattern shows an isotropic gap of very similar magnitude, while the amorphous structure does not show a gap at all. The centimeter scaled microwave structures may serve as prototypes for micrometer scaled structures with bandgaps in the technologically very interesting region of infrared (IR).
△ Less
Submitted 13 March, 2024;
originally announced March 2024.
-
Extraordinary Optical and Transport Properties of Disordered Stealthy Hyperuniform Two-Phase Media
Authors:
Jaeuk Kim,
Salvatore Torquato
Abstract:
Disordered stealthy hyperuniform (SHU) two-phase media are a special subset of hyperuniform structures with novel physical properties due to their hybrid crystal-liquid nature. We have previously shown that the strong-contrast expansion of a linear fractional form of the effective dynamic dielectric constant leads to accurate approximations for disordered two-phase composites when truncated at the…
▽ More
Disordered stealthy hyperuniform (SHU) two-phase media are a special subset of hyperuniform structures with novel physical properties due to their hybrid crystal-liquid nature. We have previously shown that the strong-contrast expansion of a linear fractional form of the effective dynamic dielectric constant leads to accurate approximations for disordered two-phase composites when truncated at the two-point level for distinctly different microstructural symmetries in three dimensions. Here, we further elucidate the extraordinary optical and transport properties of disordered SHU media. Among other results, we prove in detail that SHU layered and transversely isotropic media are perfectly transparent (i.e., no Anderson localization, in principle) within finite wavenumber intervals through the third-order terms. Remarkably, the results for these SHU media imply that there can be no Anderson localization within the predicted perfect transparency interval in practice because the localization length is much larger than any practically large sample size. We further contrast and compare the extraordinary physical properties of SHU layered, transversely isotropic, and fully isotropic media to other model nonstealthy microstructures, including their attenuation characteristics, as measured by the imaginary part of effective dielectric constant, and transport properties, as measured by the time-dependent diffusion spreadability. We demonstrate cross-property relations between them: they are positively correlated as the structures span from nonhyperuniform, nonstealthy hyperuniform, and SHU media. Establishing cross-property relations for SHU media for other wave phenomena (e.g., elastodynamics) and transport properties will also be useful. Cross-property relations are generally useful because they enable one to estimate one property, given a measurement of another.
△ Less
Submitted 12 March, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
Theoretical Prediction of the Effective Dynamic Dielectric Constant of Disordered Hyperuniform Anisotropic Composites Beyond the Long-Wavelength Regime
Authors:
Jaeuk Kim,
Salvatore Torquato
Abstract:
Torquato and Kim [Phys. Rev. X 11, 296 021002 (2021)] derived exact nonlocal strong-contrast expansions of the effective dynamic dielectric constant tensor that treat general three-dimensional (3D) two-phase composites, which are valid well beyond the long-wavelength regime. Here, we demonstrate that truncating this general rapidly converging series at the two- and three-point levels is a powerful…
▽ More
Torquato and Kim [Phys. Rev. X 11, 296 021002 (2021)] derived exact nonlocal strong-contrast expansions of the effective dynamic dielectric constant tensor that treat general three-dimensional (3D) two-phase composites, which are valid well beyond the long-wavelength regime. Here, we demonstrate that truncating this general rapidly converging series at the two- and three-point levels is a powerful theoretical tool for extracting accurate approximations suited for various microstructural symmetries. We derive such closed-form formulas applicable to transverse polarization in layered media and transverse magnetic polarization in transversely isotropic media, respectively. We use these formulas to estimate effective dielectric constant for models of 3D disordered hyperuniform layered and transversely isotropic media: nonstealthy hyperuniform and stealthy hyperuniform (SHU) media. In particular, we show that SHU media are perfectly transparent (trivially implying no Anderson localization, in principle) within finite wave number intervals through the third-order terms. For these two models, we validate that the second-order formulas, which depend on the spectral density, are already very accurate well beyond the long-wavelength regime by showing very good agreement with the finite-difference time-domain simulations. The high predictive power of the second-order formulas implies that higher-order contributions are negligibly small, and thus, it very accurately approximates multiple scattering effects. Therefore, there can be no Anderson localization in practice within the predicted perfect transparency interval in SHU media because the localization length should be very large compared to any practically large sample size. Our predictive theory provides a foundation for the inverse design of novel effective wave characteristics of disordered and statistically anisotropic structures.
△ Less
Submitted 22 December, 2023; v1 submitted 8 December, 2023;
originally announced December 2023.
-
Effective Electromagnetic Wave Properties of Disordered Stealthy Hyperuniform Layered Media Beyond the Quasistatic Regime
Authors:
Jaeuk Kim,
Salvatore Torquato
Abstract:
Disordered stealthy hyperuniform dielectric composites exhibit novel electromagnetic wave transport properties in two and three dimensions. Here, we carry out the first study of the electromagnetic properties of one-dimensional (1D) disordered stealthy hyperuniform layered media. From an exact nonlocal theory, we derive an approximation formula for the effective dynamic dielectric constant tensor…
▽ More
Disordered stealthy hyperuniform dielectric composites exhibit novel electromagnetic wave transport properties in two and three dimensions. Here, we carry out the first study of the electromagnetic properties of one-dimensional (1D) disordered stealthy hyperuniform layered media. From an exact nonlocal theory, we derive an approximation formula for the effective dynamic dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_q,ω)$ of general 1D media that is valid well beyond the quasistatic regime and apply it to 1D stealthy hyperuniform systems. We consider incident waves of transverse polarization, frequency $ω$, and wavenumber $k_q$. Our formula for ${\boldsymbol \varepsilon}_e({k}_q,ω)$, which is given in terms of the spectral density, leads to a closed-form relation for the transmittance $T$. Our theoretical predictions are in excellent agreement with finite-difference time-domain (FDTD) simulations. Stealthy hyperuniform layered media have perfect transparency intervals up to a finite wavenumber, implying no Anderson localization, but non-stealthy hyperuniform media are not perfectly transparent. Our predictive theory provides a new path for the inverse design of the wave characteristics of disordered layered media, which are readily fabricated, by engineering their spectral densities.
△ Less
Submitted 31 July, 2023; v1 submitted 22 May, 2023;
originally announced May 2023.
-
Sub-diffusive wave transport and weak localization transition in three-dimensional stealthy hyperuniform disordered systems
Authors:
Fabrizio Sgrignuoli,
Salvatore Torquato,
Luca Dal Negro
Abstract:
The purpose of this work is to understand the fundamental connection between structural correlations and light localization in three-dimensional (3D) open scattering systems of finite size. We numerically investigate the transport of vector electromagnetic waves scattered by resonant electric dipoles spatially arranged in 3D space by stealthy hyperuniform disordered point patterns. Three-dimension…
▽ More
The purpose of this work is to understand the fundamental connection between structural correlations and light localization in three-dimensional (3D) open scattering systems of finite size. We numerically investigate the transport of vector electromagnetic waves scattered by resonant electric dipoles spatially arranged in 3D space by stealthy hyperuniform disordered point patterns. Three-dimensional stealthy hyperuniform disordered systems (3D-SHDS) are engineered with different structural correlation properties determined by their degree of stealthiness $χ$. Such fine control of exotic states of amorphous matter enables the systematic design of optical media that interpolate in a tunable fashion between uncorrelated random structures and crystalline materials. By solving the electromagnetic multiple scattering problem using Green's matrix spectral method, we establish a transport phase diagram that demonstrates a distinctive transition from a diffusive to a weak localization regime beyond a critical scattering density that depends on $χ$. The transition is characterized by studying the Thouless number and the spectral statistics of the scattering resonances. In particular, by tuning the $χ$ parameter, we demonstrate large spectral gaps and suppressed sub-radiant proximity resonances, facilitating light localization. Moreover, consistently with previous studies, our results show a region of the transport phase diagram where the investigated scattering systems become transparent. Our work provides a systematic description of the transport and weak localization properties of light in stealthy hyperuniform structures and motivates the engineering of novel photonic systems with enhanced light-matter interactions for applications to both classical and quantum devices.
△ Less
Submitted 21 February, 2022; v1 submitted 8 September, 2021;
originally announced September 2021.
-
Gap Sensitivity Reveals Universal Behaviors in Optimized Photonic Crystal and Disordered Networks
Authors:
Michael A. Klatt,
Paul J. Steinhardt,
Salvatore Torquato
Abstract:
Through an extensive series of high-precision numerical computations of the optimal complete photonic band gap (PBG) as a function of dielectric contrast $α$ for a variety of crystal and disordered heterostructures, we reveal striking universal behaviors of the gap sensitivity $\mathcal{S}(α)\equiv dΔ(α)/dα$, the first derivative of the optimal gap-to-midgap ratio $Δ(α)$. In particular, for all ou…
▽ More
Through an extensive series of high-precision numerical computations of the optimal complete photonic band gap (PBG) as a function of dielectric contrast $α$ for a variety of crystal and disordered heterostructures, we reveal striking universal behaviors of the gap sensitivity $\mathcal{S}(α)\equiv dΔ(α)/dα$, the first derivative of the optimal gap-to-midgap ratio $Δ(α)$. In particular, for all our crystal networks, $\mathcal{S}(α)$ takes a universal form that is well approximated by the analytic formula for a one-dimensional quarter-wave stack, $\mathcal{S}_{\text{QWS}}(α)$. Even more surprisingly, the values of $\mathcal{S}(α)$ for our disordered networks converge to $\mathcal{S}_{\text{QWS}}(α)$ for sufficiently large $α$. A deeper understanding of the simplicity of this universal behavior may provide fundamental insights about PBG formation and guidance in the design of novel photonic heterostructures.
△ Less
Submitted 15 July, 2021;
originally announced July 2021.
-
Generation and Structural Characterization of Debye Random Media
Authors:
Zheng Ma,
Salvatore Torquato
Abstract:
In their seminal paper on scattering by an inhomogeneous solid, Debye and coworkers proposed a simple exponentially decaying function for the two-point correlation function of an idealized class of two-phase random media. Such {\it Debye random media}, which have been shown to be realizable, are singularly distinct from all other models of two-phase media in that they are entirely defined by their…
▽ More
In their seminal paper on scattering by an inhomogeneous solid, Debye and coworkers proposed a simple exponentially decaying function for the two-point correlation function of an idealized class of two-phase random media. Such {\it Debye random media}, which have been shown to be realizable, are singularly distinct from all other models of two-phase media in that they are entirely defined by their one- and two-point correlation functions. To our knowledge, there has been no determination of other microstructural descriptors of Debye random media. In this paper, we generate Debye random media in two dimensions using an accelerated Yeong-Torquato construction algorithm. We then ascertain microstructural descriptors of the constructed media, including their surface correlation functions, pore-size distributions, lineal-path function, and chord-length probability density function. Accurate semi-analytic and empirical formulas for these descriptors are devised. We compare our results for Debye random media to those of other popular models (overlapping disks and equilibrium hard disks), and find that the former model possesses a wider spectrum of hole sizes, including a substantial fraction of large holes. Our algorithm can be applied to generate other models defined by their two-point correlation functions, and their other microstructural descriptors can be determined and analyzed by the procedures laid out here.
△ Less
Submitted 20 October, 2020;
originally announced October 2020.
-
Predicting permeability via statistical learning on higher-order microstructural information
Authors:
Magnus Röding,
Zheng Ma,
Salvatore Torquato
Abstract:
Quantitative structure-property relationships are crucial for the understanding and prediction of the physical properties of complex materials. For fluid flow in porous materials, characterizing the geometry of the pore microstructure facilitates prediction of permeability, a key property that has been extensively studied in material science, geophysics and chemical engineering. In this work, we s…
▽ More
Quantitative structure-property relationships are crucial for the understanding and prediction of the physical properties of complex materials. For fluid flow in porous materials, characterizing the geometry of the pore microstructure facilitates prediction of permeability, a key property that has been extensively studied in material science, geophysics and chemical engineering. In this work, we study the predictability of different structural descriptors via both linear regressions and neural networks. A large data set of 30,000 virtual, porous microstructures of different types is created for this end. We compute permeabilities of these structures using the lattice Boltzmann method, and characterize the pore space geometry using one-point correlation functions (porosity, specific surface), two-point surface-surface, surface-void, and void-void correlation functions, as well as the geodesic tortuosity as an implicit descriptor. Then, we study the prediction of the permeability using different combinations of these descriptors. We obtain significant improvements of performance when compared to a Kozeny-Carman regression with only lowest-order descriptors (porosity and specific surface). We find that combining all three two-point correlation functions and tortuosity provides the best prediction of permeability, with the void-void correlation function being the most informative individual descriptor. Moreover, the combination of porosity, specific surface, and geodesic tortuosity provides very good predictive performance. This shows that higher-order correlation functions are extremely useful for forming a general model for predicting physical properties of complex materials. Additionally, our results suggest that neural networks are superior to the more conventional regression methods for establishing quantitative structure-property relationships.
△ Less
Submitted 7 September, 2020;
originally announced September 2020.
-
Nonlocal Effective Electromagnetic Wave Characteristics of Composite Media: Beyond the Quasistatic Regime
Authors:
Salvatore Torquato,
Jaeuk Kim
Abstract:
We derive exact nonlocal expressions for the effective dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ of disordered two-phase composites and metamaterials from first principles. This formalism extends the long-wavelength limitations of conventional homogenization estimates of ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ for arbitrary microstructures so that it can capture sp…
▽ More
We derive exact nonlocal expressions for the effective dielectric constant tensor ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ of disordered two-phase composites and metamaterials from first principles. This formalism extends the long-wavelength limitations of conventional homogenization estimates of ${\boldsymbol \varepsilon}_e({\bf k}_I, ω)$ for arbitrary microstructures so that it can capture spatial dispersion well beyond the quasistatic regime (where $ω$ and ${\bf k}_I$ are frequency and wavevector of the incident radiation). This is done by deriving nonlocal strong-contrast expansions that exactly account for multiple scattering for the range of wavenumbers for which our extended homogenization theory applies, i.e., $0 \le |{\bf k}_I| \ell \lesssim 1$ (where $\ell$ is a characteristic heterogeneity length scale). Due to the fast-convergence properties of such expansions, their lower-order truncations yield accurate closed-form approximate formulas for ${\varepsilon}_e({\bf k}_I,ω)$ that incorporate microstructural information via the spectral density, which is easy to compute for any composite. The accuracy of these microstructure-dependent approximations is validated by comparison to full-waveform simulation methods for both 2D and 3D ordered and disordered models of composite media. Thus, our closed-form formulas enable one to predict accurately and efficiently the effective wave characteristics well beyond the quasistatic regime without having to perform full-blown simulations. Among other results, we show that certain disordered hyperuniform particulate composites exhibit novel wave characteristics. Our results demonstrate that one can design the effective wave characteristics of a disordered composite by engineering the microstructure to possess tailored spatial correlations at prescribed length scales.
△ Less
Submitted 16 March, 2021; v1 submitted 1 July, 2020;
originally announced July 2020.
-
A minimal statistical-mechanical model for multihyperuniform patterns in avian retina
Authors:
Enrique Lomba,
Jean-Jacques Weis,
Leandro Guisández,
Salvatore Torquato
Abstract:
Birds are known for their extremely acute sense of vision. The very peculiar structural distribution of five different types of cones in the retina underlies this exquisite ability to sample light. It was recently found that each cone population as well as their total population display a disordered pattern in which long wave-length density fluctuations vanish. This property, known as hyperuniform…
▽ More
Birds are known for their extremely acute sense of vision. The very peculiar structural distribution of five different types of cones in the retina underlies this exquisite ability to sample light. It was recently found that each cone population as well as their total population display a disordered pattern in which long wave-length density fluctuations vanish. This property, known as hyperuniformity is also present in perfect crystals. In situations like the avian retina in which both the global structure and that of each component display hyperuniformity, the system is said to be multi-hyperuniform. In this work, we aim at devising a minimal statistical-mechanical model that can reproduce the main features of the spatial distribution of photoreceptors in avian retina, namely the presence of disorder, multi-hyperuniformity and local hetero-coordination. This last feature is key to avoid local clustering of the same type of photoreceptors, an undesirable feature for the efficient sampling of light. For this purpose we formulate a simple model that definitively exhibits the required structural properties, namely an equimolar three-component mixture (one component to sample each primary color, red, green, and blue) of non-additive hard disks to which a long-range logarithmic repulsion is added between like particles. A Voronoi analysis of our idealized system of photoreceptors shows that the space-filling Voronoi polygons interestingly display a rather uniform area distribution, symmetrically centered around that of a regular lattice, a structural property also found in human retina. Disordered multi-hyperuniformity offers an alternative to generate photoreceptor patterns with minimal long-range concentration and density fluctuations. This is the key to overcome the difficulties in devising an efficient visual system in which crystal-like order is absent.
△ Less
Submitted 4 May, 2020;
originally announced May 2020.
-
Realizable hyperuniform and nonhyperuniform particle configurations with targeted spectral functions via effective pair interactions
Authors:
Ge Zhang,
Salvatore Torquato
Abstract:
The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function $g_2(r)$ or its corresponding structure factor $S(k)$ is of great fundamental and practical importance. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classic…
▽ More
The capacity to identify realizable many-body configurations associated with targeted functional forms for the pair correlation function $g_2(r)$ or its corresponding structure factor $S(k)$ is of great fundamental and practical importance. A major aim of this paper is to expand our theoretical knowledge of the class of pair correlation functions or structure factors that are realizable by classical disordered ensembles of particle configurations. We first introduce a theoretical formalism that provides a means to draw classical particle configurations from canonical ensembles with certain pairwise-additive potentials that could correspond to targeted analytical functional forms for the structure factor. This formulation enables us to devise an improved algorithm to construct systematically canonical-ensemble particle configurations with such targeted pair statistics. We test the algorithm by targeting several different structure factors across dimensions that are known to be realizable and one hyperuniform target that is known to be nontrivially unrealizable. Subsequently, we also target several families of structure-factor functions that meet the known necessary realizability conditions but were heretofore not known to be realizable. Moreover, we also explore unusual nonhyperuniform targets, including "hyposurficial" and "anti-hyperuniform" examples. In all of these instances, the targeted structure factors were achieved with high accuracy, suggesting that they are indeed realizable by equilibrium configurations with pairwise interactions at positive temperatures. Our findings lead us to the conjecture that any realizable structure factor corresponding to either a translationally invariant equilibrium or nonequilibrium system can be attained by an equilibrium ensemble involving only effective pair interactions.
△ Less
Submitted 5 April, 2020;
originally announced April 2020.
-
Multifunctional Composites for Elastic and Electromagnetic Wave Propagation
Authors:
Jaeuk Kim,
Salvatore Torquato
Abstract:
Composites are ideally suited to achieve desirable multifunctional effective properties since the best properties of different materials can be judiciously combined with designed microstructures. Here we establish cross-property relations for two-phase composite media that link effective elastic and electromagnetic wave characteristics to one another, including the respective effective wave speeds…
▽ More
Composites are ideally suited to achieve desirable multifunctional effective properties since the best properties of different materials can be judiciously combined with designed microstructures. Here we establish cross-property relations for two-phase composite media that link effective elastic and electromagnetic wave characteristics to one another, including the respective effective wave speeds and attenuation coefficients, which facilitate multifunctional material design. This is achieved by deriving accurate formulas for the effective electromagnetic and elastodynamic properties that depend on the wavelengths of the incident waves and the microstructure via the spectral density. Our formulas enable us to explore the wave characteristics of a broad class of disordered microstructures because they apply, unlike conventional formulas, for a wide range of incident wavelengths, i.e., well beyond the long-wavelength regime. This capability enables us to study the dynamic properties of exotic disordered ``hyperuniform'' composites that can have advantages over crystalline ones, such as nearly optimal, direction-independent properties and robustness against defects. We specifically show that disordered ``stealthy'' hyperuniform microstructures exhibit novel wave characteristics, e.g., low-pass filters that transmit waves ``isotropically'' up to a finite wavenumber. Our cross-property relations for the effective wave characteristics can be applied to design multifunctional composites via inverse techniques. Design examples include structural components that require high stiffness and electromagnetic absorption, heat-sinks for CPUs, and sound-absorbing housings for motors that have to efficiently emit thermal radiation and suppress mechanical vibrations, and nondestructive evaluation of the elastic moduli of materials from the effective dielectric response.
△ Less
Submitted 11 May, 2020; v1 submitted 19 August, 2019;
originally announced August 2019.
-
Precise algorithms to compute surface correlation functions of two-phase heterogeneous media and their applications
Authors:
Zheng Ma,
Salvatore Torquato
Abstract:
The quantitative characterization of the microstructure of random heterogeneous media in $d$-dimensional Euclidean space $\mathbb{R}^d$ via a variety of $n$-point correlation functions is of great importance, since the respective infinite set determines the effective physical properties of the media. In particular, surface-surface $F_{ss}$ and surface-void $F_{sv}$ correlation functions (obtainabl…
▽ More
The quantitative characterization of the microstructure of random heterogeneous media in $d$-dimensional Euclidean space $\mathbb{R}^d$ via a variety of $n$-point correlation functions is of great importance, since the respective infinite set determines the effective physical properties of the media. In particular, surface-surface $F_{ss}$ and surface-void $F_{sv}$ correlation functions (obtainable from radiation scattering experiments) contain crucial interfacial information that enables one to estimate transport properties of the media (e.g., the mean survival time and fluid permeability) and complements the information content of the conventional two-point correlation function. However, the current technical difficulty involved in sampling surface correlation functions has been a stumbling block in their widespread use. We first present a concise derivation of the small-$r$ behaviors of these functions, which are linked to the \textit{mean curvature} of the system. Then we demonstrate that one can reduce the computational complexity of the problem by extracting the necessary interfacial information from a cut of the system with an infinitely long line. Accordingly, we devise algorithms based on this idea and test them for two-phase media in continuous and discrete spaces. Specifically for the exact benchmark model of overlapping spheres, we find excellent agreement between numerical and exact results. We compute surface correlation functions and corresponding local surface-area variances for a variety of other model microstructures, including hard spheres in equilibrium, decorated "stealthy" patterns, as well as snapshots of evolving pattern formation processes (e.g., spinodal decomposition). It is demonstrated that the precise determination of surface correlation functions provides a powerful means to characterize a wide class of complex multiphase microstructures.
△ Less
Submitted 31 July, 2018;
originally announced August 2018.
-
Searching for crystal-ice domains in amorphous ices
Authors:
Fausto Martelli,
Nicolas Giovambattista,
Salvatore Torquato,
Roberto Car
Abstract:
We employ classical molecular dynamics simulations to investigate the molecular-level structure of water during the isothermal compression of hexagonal ice (I$h$) and low-density amorphous (LDA) ice at low temperatures. In both cases, the system transforms to high-density amorphous ice (HDA) via a first-order-like phase transition. We employ a sensitive local order metric (LOM) [Martelli et. al.,…
▽ More
We employ classical molecular dynamics simulations to investigate the molecular-level structure of water during the isothermal compression of hexagonal ice (I$h$) and low-density amorphous (LDA) ice at low temperatures. In both cases, the system transforms to high-density amorphous ice (HDA) via a first-order-like phase transition. We employ a sensitive local order metric (LOM) [Martelli et. al., Phys. Rev. B, 97, 064105 (2018)], that can discriminate among different crystalline and non crystalline ice structures and is based on the positions of the oxygen atoms in the first and/or second hydration shell. Our results confirm that LDA and HDA are indeed amorphous, i.e., they lack of polydispersed ice domains. Interestingly, HDA contains a small number of domains that are reminiscent of the unit cell of ice IV, although the hydrogen-bond network (HBN) of these domains differ from the HBN of ice IV. The presence of ice IV-like domains provides some support to the hypothesis that HDA could be the result of a detour on the HBN rearrangement along the I$h$-to-ice IV pressure induced transformation. Both nonequilibrium LDA-to-HDA and I$h$-to-HDA transformations are two-steps processes where a small distortion of the HBN first occurs at low pressures and then, a sudden, extensive re-arrangement of hydrogen bonds at the corresponding transformation pressure follows. Interestingly, the I$h$-to-HDA and LDA-to-HDA transformations occur when LDA and I$h$ have similar local order, as quantified by the site-averaged LOMs. Since I$h$ has a perfect tetrahedral HBN, while LDA does not, it follows that higher pressures are needed to transform I$h$ into HDA than that for the conversion of LDA to HDA. In correspondence with both first-order-like phase transitions, the samples are composed of a large HDA cluster that percolates within the I$h$/LDA samples.
△ Less
Submitted 7 May, 2018;
originally announced May 2018.
-
Structural characterization and statistical-mechanical model of epidermal patterns
Authors:
Duyu Chen,
Wen Yih Aw,
Danelle Devenport,
Salvatore Torquato
Abstract:
In proliferating epithelia of mammalian skin, cells of irregular polygonal-like shapes pack into complex nearly flat two-dimensional structures that are pliable to deformations. In this work, we employ various sensitive correlation functions to quantitatively characterize structural features of evolving packings of epithelial cells across length scales in mouse skin. We find that the pair statisti…
▽ More
In proliferating epithelia of mammalian skin, cells of irregular polygonal-like shapes pack into complex nearly flat two-dimensional structures that are pliable to deformations. In this work, we employ various sensitive correlation functions to quantitatively characterize structural features of evolving packings of epithelial cells across length scales in mouse skin. We find that the pair statistics in direct and Fourier spaces of the cell centroids in the early stages of embryonic development show structural directional dependence, while in the late stages the patterns tend towards statistically isotropic states. We construct a minimalist four-component statistical-mechanical model involving effective isotropic pair interactions consisting of hard-core repulsion and extra short-ranged soft-core repulsion beyond the hard core, whose length scale is roughly the same as the hard core. The model parameters are optimized to match the sample pair statistics in both direct and Fourier spaces. By doing this, the parameters are biologically constrained. Our model predicts essentially the same polygonal shape distribution and size disparity of cells found in experiments as measured by Voronoi statistics. Moreover, our simulated equilibrium liquid-like configurations are able to match other nontrivial unconstrained statistics, which is a testament to the power and novelty of the model. We discuss ways in which our model might be extended so as to better understand morphogenesis (in particular the emergence of planar cell polarity), wound-healing, and disease progression processes in skin, and how it could be applied to the design of synthetic tissues.
△ Less
Submitted 6 December, 2016;
originally announced December 2016.
-
A Cellular Automaton Model for Tumor Dormancy: Emergence of a Proliferative Switch
Authors:
Duyu Chen,
Yang Jiao,
Salvatore Torquato
Abstract:
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 po…
▽ More
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.
△ Less
Submitted 17 October, 2014;
originally announced October 2014.
-
Avian photoreceptor patterns represent a disordered hyperuniform solution to a multiscale packing problem
Authors:
Yang Jiao,
Timothy Lau,
Haralampos Hatzikirou,
Michael Meyer-Hermann,
Joseph C. Corbo,
Salvatore Torquato
Abstract:
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 photorecep…
▽ More
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.
△ Less
Submitted 25 February, 2014;
originally announced February 2014.
-
Photonic band gap in isotropic hyperuniform disordered solids with low dielectric contrast
Authors:
Weining Man,
Marian Florescu,
Kazue Matsuyama,
Polin Yadak,
Geev Nahal,
Seyed Hashemizad,
Eric Williamson,
Paul Steinhardt,
Salvatore Torquato,
Paul Chaikin
Abstract:
We report the first experimental demonstration of a TE-polarization photonic band gap (PBG) in a 2D isotropic hyperuniform disordered solid (HUDS) made of dielectric media with a index contrast of 1.6:1, very low for PBG formation. The solid is composed of a connected network of dielectric walls enclosing air-filled cells. Direct comparison with photonic crystals and quasicrystals permitted us to…
▽ More
We report the first experimental demonstration of a TE-polarization photonic band gap (PBG) in a 2D isotropic hyperuniform disordered solid (HUDS) made of dielectric media with a index contrast of 1.6:1, very low for PBG formation. The solid is composed of a connected network of dielectric walls enclosing air-filled cells. Direct comparison with photonic crystals and quasicrystals permitted us to investigate band-gap properties as a function of increasing rotational isotropy. We present results from numerical simulations proving that the PBG observed experimentally for HUDS at low index contrast has zero density of states. The PBG is associated with the energy difference between complementary resonant modes above and below the gap, with the field predominantly concentrated in the air or in the dielectric. The intrinsic isotropy of HUDS may offer unprecedented flexibilities and freedom in applications (i. e. defect architecture design) not limited by crystalline symmetries.
△ Less
Submitted 12 November, 2013;
originally announced November 2013.
-
Isotropic Band Gaps and Freeform Waveguides Observed in Hyperuniform Disordered Photonic Solids
Authors:
Weining Man,
Marian Florescu,
Eric Paul Williamson,
Yingquan He,
Seyed Reza Hashemizad,
Brian Y. C. Leung,
Devin Robert Liner,
Salvatore Torquato,
Paul M. Chaikin,
Paul J. Steinhardt
Abstract:
Recently, disordered photonic media and random textured surfaces have attracted increasing attention as strong light diffusers with broadband and wide-angle properties. We report the first experimental realization of an isotropic complete photonic band gap (PBG) in a two-dimensional (2D) disordered dielectric structure. This structure is designed by a constrained-optimization method, which combine…
▽ More
Recently, disordered photonic media and random textured surfaces have attracted increasing attention as strong light diffusers with broadband and wide-angle properties. We report the first experimental realization of an isotropic complete photonic band gap (PBG) in a two-dimensional (2D) disordered dielectric structure. This structure is designed by a constrained-optimization method, which combines advantages of both isotropy due to disorder and controlled scattering properties due to low density fluctuations (hyperuniformity) and uniform local topology. Our experiments use a modular design composed of Al2O3 walls and cylinders arranged in a hyperuniform disordered network. We observe a complete PBG in the microwave region, in good agreement with theoretical simulations, and show that the intrinsic isotropy of this novel class of PBG materials enables remarkable design freedom, including the realization of waveguides with arbitrary bending angles impossible in photonic crystals. This first experimental verification of a complete PBG and realization of functional defects in this new class of materials demonstrates their potential as building blocks for precise manipulation of photons in planar optical micro-circuits and has implications for disordered acoustic and electronic bandgap materials.
△ Less
Submitted 11 November, 2013;
originally announced November 2013.
-
Optical cavities and waveguides in hyperuniform disordered photonic solids
Authors:
Marian Florescu,
Salvatore Torquato,
Paul Steinhardt
Abstract:
Using finite difference time domain and band structure computer simulations, we show that it is possible to construct optical cavities and waveguide architectures in hyperuniform disordered photonic solids that are unattainable in photonic crystals. The cavity modes can be classified according to the symmetry (monopole, dipole, quadrupole,etc.) of the confined electromagnetic wave pattern. Owing t…
▽ More
Using finite difference time domain and band structure computer simulations, we show that it is possible to construct optical cavities and waveguide architectures in hyperuniform disordered photonic solids that are unattainable in photonic crystals. The cavity modes can be classified according to the symmetry (monopole, dipole, quadrupole,etc.) of the confined electromagnetic wave pattern. Owing to the isotropy of the band gaps characteristic of hyperuniform disordered solids, high-quality waveguides with freeform geometries (e.g., arbitrary bending angles) can be constructed that have no analogue in periodic or quasiperiodic solids. These capabilities have implications for many photonic applications.
△ Less
Submitted 11 November, 2013;
originally announced November 2013.
-
Quantitative Characterization of the Microstructure and Transport Properties of Biopolymer Networks
Authors:
Yang Jiao,
Salvatore Torquato
Abstract:
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 networ…
▽ More
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.
△ Less
Submitted 21 March, 2012;
originally announced March 2012.
-
Diversity of Dynamics and Morphologies of Invasive Solid Tumors
Authors:
Yang Jiao,
Sal Torquato
Abstract:
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 micro…
▽ More
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.
△ Less
Submitted 4 January, 2012; v1 submitted 3 January, 2012;
originally announced January 2012.
-
Emergent Behaviors from A Cellular Automaton Model for Invasive Tumor Growth in Heterogeneous Microenvironments
Authors:
Yang Jiao,
Salvatore Torquato
Abstract:
Understanding tumor invasion and metastasis is of crucial importance for both fundamental cancer research and clinical practice. In vitro experiments have established that the invasive growth of malignant tumors is characterized by the dendritic invasive branches composed of chains of tumor cells emanating from the primary tumor mass. The preponderance of previous tumor simulations focused on non-…
▽ More
Understanding tumor invasion and metastasis is of crucial importance for both fundamental cancer research and clinical practice. In vitro experiments have established that the invasive growth of malignant tumors is characterized by the dendritic invasive branches composed of chains of tumor cells emanating from the primary tumor mass. The preponderance of previous tumor simulations focused on non-invasive (or proliferative) growth. The formation of the invasive cell chains and their interactions with the primary tumor mass and host microenvironment are not well understood. Here, we present a novel cellular automaton (CA) model that enables one to efficiently simulate invasive tumor growth in a heterogeneous host microenvironment. By taking into account a variety of microscopic-scale tumor-host interactions, including the short-range mechanical interactions between tumor cells and tumor stroma, degradation of extracellular matrix by the invasive cells and oxygen/nutrient gradient driven cell motions, our CA model predicts a rich spectrum of growth dynamics and emergent behaviors of invasive tumors. Besides robustly reproducing the salient features of dendritic invasive growth, such as least resistance and intrabranch homotype attraction, we also predict nontrivial coupling of the growth dynamics of the primary tumor mass and the invasive cells. In addition, we show that the properties of the host microenvironment can significantly affect tumor morphology and growth dynamics, emphasizing the importance of understanding the tumor-host interaction. The capability of our CA model suggests that well-developed in silico tools could eventually be utilized in clinical situations to predict neoplastic progression and propose individualized optimal treatment strategies.
△ Less
Submitted 3 January, 2012; v1 submitted 2 November, 2011;
originally announced November 2011.
-
Phase diagram and structural diversity of the densest binary sphere packings
Authors:
Adam B. Hopkins,
Yang Jiao,
Frank H. Stillinger,
Salvatore Torquato
Abstract:
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)]…
▽ More
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.
△ Less
Submitted 10 August, 2011;
originally announced August 2011.
-
Effects of Random Link Removal on the Photonic Band Gaps of Honeycomb Networks
Authors:
Marian Florescu,
Salvatore Torquato,
Paul J. Steinhardt
Abstract:
We explore the effects of random link removal on the photonic band gaps of honeycomb networks. Missing or incomplete links are expected to be common in practical realizations of this class of connected network structures due to unavoidable flaws in the fabrication process. We focus on the collapse of the photonic band gap due to the defects induced by the link removal. We show that the photonic ba…
▽ More
We explore the effects of random link removal on the photonic band gaps of honeycomb networks. Missing or incomplete links are expected to be common in practical realizations of this class of connected network structures due to unavoidable flaws in the fabrication process. We focus on the collapse of the photonic band gap due to the defects induced by the link removal. We show that the photonic band gap is quite robust against this type of random decimation and survives even when almost 58% of the network links are removed.
△ Less
Submitted 7 November, 2010;
originally announced November 2010.
-
Complete Band Gaps in 2D Photonic Quasicrystals
Authors:
Marian Florescu,
Salvatore Torquato,
Paul J. Steinhardt
Abstract:
We introduce a novel optimization method to design the first examples of photonic quasicrystals with substantial, complete photonic band gaps (PBGs): that is, a range of frequencies over which electromagnetic wave propagation is forbidden for all directions and polarizations. The method can be applied to photonic quasicrystals with arbitrary rotational symmetry; here, we illustrate the results for…
▽ More
We introduce a novel optimization method to design the first examples of photonic quasicrystals with substantial, complete photonic band gaps (PBGs): that is, a range of frequencies over which electromagnetic wave propagation is forbidden for all directions and polarizations. The method can be applied to photonic quasicrystals with arbitrary rotational symmetry; here, we illustrate the results for 5- and 8-fold symmetric quasicrystals. The optimized band gaps are highly isotropic, which may offer advantages over photonic crystals for certain applications.
△ Less
Submitted 20 July, 2010;
originally announced July 2010.
-
Designer disordered materials with large complete photonic band gaps
Authors:
Marian Florescu,
Salvatore Torquato,
Paul J. Steinhardt
Abstract:
We present designs of 2D isotropic, disordered photonic materials of arbitrary size with complete band gaps blocking all directions and polarizations. The designs with the largest gaps are obtained by a constrained optimization method that starts from a hyperuniform disordered point pattern, an array of points whose number variance within a spherical sampling window grows more slowly than the volu…
▽ More
We present designs of 2D isotropic, disordered photonic materials of arbitrary size with complete band gaps blocking all directions and polarizations. The designs with the largest gaps are obtained by a constrained optimization method that starts from a hyperuniform disordered point pattern, an array of points whose number variance within a spherical sampling window grows more slowly than the volume. We argue that hyperuniformity, combined with uniform local topology and short-range geometric order, can explain how complete photonic band gaps are possible without long-range translational order. We note the ramifications for electronic and phononic band gaps in disordered materials.
△ Less
Submitted 20 July, 2010;
originally announced July 2010.
-
Neighbor List Collision-Driven Molecular Dynamics Simulation for Nonspherical Particles. I. Algorithmic Details II. Applications to Ellipses and Ellipsoids
Authors:
Aleksandar Donev,
Salvatore Torquato,
Frank H. Stillinger
Abstract:
In the first part of a series of two papers, we present in considerable detail a collision-driven molecular dynamics algorithm for a system of nonspherical particles, within a parallelepiped simulation domain, under both periodic or hard-wall boundary conditions. The algorithm extends previous event-driven molecular dynamics algorithms for spheres. We present a novel partial-update near-neighbor…
▽ More
In the first part of a series of two papers, we present in considerable detail a collision-driven molecular dynamics algorithm for a system of nonspherical particles, within a parallelepiped simulation domain, under both periodic or hard-wall boundary conditions. The algorithm extends previous event-driven molecular dynamics algorithms for spheres. We present a novel partial-update near-neighbor list (NNL) algorithm that is superior to previous algorithms at high densities, without compromising the correctness of the algorithm. This efficiency of the algorithm is further increased for systems of very aspherical particles by using bounding sphere complexes (BSC). In the second part of this series of papers we apply the algorithm presented in the first part of this series of papers to systems of hard ellipses and ellipsoids. The theoretical machinery needed to treat such particles, including the overlap potentials, is developed in full detail. We describe an algorithm for predicting the time of collision for two moving ellipses or ellipsoids. We present performance results for our implementation of the algorithm. The practical utility of the algorithm is demonstrated by presenting several interesting physical applications, including the generation of jammed packings inside spherical containers, the study of contact force chains in jammed packings, and melting the densest-known equilibrium crystals of prolate spheroids.
△ Less
Submitted 17 May, 2004;
originally announced May 2004.