Hidden Hyperuniformity in Soft Polymeric Materials

Hyperuniform disordered solids with crystal-like stability

Hyperuniform disordered solids, characterised by unusually suppressed density fluctuations at low wavenumbers (q), are of great interest due to their potentially distinct properties as a unique glass state. From the jamming perspective, there is ongoing debate about the relationship between hyperuniformity and the jamming transition, as well as whether hyperuniformity persists above the jamming point. Here, we successfully generate over-jammed disordered solids exhibiting the strongest class of hyperuniformity, characterised by a power-law density spectrum (qα with α = 4). By decompressing both hyperuniform and conventional over-jammed packings to their respective marginally jammed states, we identify protocol-independent exponents: α ≈ 0.25 for density hyperuniformity and α ≈ 2 for contact-number hyperuniformity, both associated with the jamming transition. Although both marginally jammed and conventional over-jammed packings exhibit marginal stability, we demonstrate that hyperuniform over-jammed packings possess exceptional stability across vibrational, kinetic, thermodynamic, and mechanical properties-similar to crystals. These findings suggest that hyperuniform over-jammed packings offer crucial insights into the ideal disordered solid state and stand out as promising candidates for disordered metamaterials, uniquely combining hyperuniformity with ultrastability.

Colloidal Dispersions of Sterically and Electrostatically Stabilized PbS Quantum Dots: Structure Factors, Second Virial Coefficients, and Film-Forming Properties

Electrostatically stabilized nanocrystals (NCs) and, in particular, quantum dots (QDs) hold promise for forming strongly coupled superlattices due to their compact and electronically conductive surface ligands. However, studies of the colloidal dispersion and interparticle interactions of electrostatically stabilized sub-10 nm NCs have been limited, hindering the optimization of their colloidal stability and self-assembly. In this study, we employed small-angle X-ray scattering (SAXS) experiments to investigate the interparticle interactions and arrangement of PbS QDs with thiostannate ligands (PbS-Sn2S64-) in polar solvents. The study reveals significant deviations from the ideal solution behavior in electrostatically stabilized QD dispersions. Our results demonstrate that PbS-Sn2S64- QDs exhibit long-range interactions within the solvent, in contrast to the short-range steric repulsion characteristic of PbS QDs with oleate ligands (PbS-OA). Introducing highly charged multivalent electrolytes screens electrostatic interactions between charged QDs, reducing the length scale of the repulsive interactions. Furthermore, we calculated the second virial (B2) coefficients from SAXS data, providing insights into how surface chemistry, solvent, and size influence pair potentials. Finally, we explore the influence of long-range interparticle interactions of PbS-Sn2S64- QDs on the morphology of films produced by drying or spin-coating colloidal solutions. The long-range repulsive term of PbS-Sn2S64- QDs promotes the formation of amorphous films, and screening the electrostatic repulsion by the addition of an electrolyte enables the formation of crystalline domains. These findings highlight the critical role of NC-NC interactions in tailoring the properties of functional materials made of colloidal NCs.

Universal Hyperuniform Organization in Looped Leaf Vein Networks

The leaf vein network is a hierarchical vascular system that transports water and nutrients to the leaf cells. The thick primary veins form a branched network, while the secondary veins can develop closed loops forming a well-defined cellular structure. Through extensive analysis of a variety of distinct leaf species, we discover that the apparently disordered cellular structures of the secondary vein networks exhibit a universal hyperuniform organization and possess a hidden order on large scales. Disorder hyperuniform systems lack conventional long-range order, yet they completely suppress normalized infinite-wavelength density fluctuations like crystals. Specifically, we find that the distributions of the geometric centers associated with the vein network loops possess a vanishing static structure factor in the limit that the wave number k goes to 0, i.e., S(k)∼k^{α}, where α≈0.64±0.021, providing an example of class III hyperuniformity in biology. This hyperuniform organization leads to superior efficiency of diffusive transport, as evidenced by the much faster convergence of the time-dependent spreadability S(t) to its longtime asymptotic limit, compared to that of other uncorrelated or correlated disordered but nonhyperuniform organizations. Our results also have implications for the discovery and design of novel disordered network materials with optimal transport properties.

Vibrational properties of disordered stealthy hyperuniform 1D atomic chains

Disorder hyperuniformity is a recently discovered exotic state of many-body systems that possess a hidden order in between that of a perfect crystal and a completely disordered system. Recently, this novel disordered state has been observed in a number of quantum materials including amorphous 2D graphene and silica, which are endowed with unexpected electronic transport properties. Here, we numerically investigate 1D atomic chain models, including perfect crystalline, disordered stealthy hyperuniform (SHU) as well as randomly perturbed atom packing configurations to obtain a quantitative understanding of how the unique SHU disorder affects the vibrational properties of these low-dimensional materials. We find that the disordered SHU chains possess lower cohesive energies compared to the randomly perturbed chains, implying their potential reliability in experiments. Our inverse partition ratio (IPR) calculations indicate that the SHU chains can support fully delocalized states just like perfect crystalline chains over a wide range of frequencies, i.e.ω∈(0,100)cm-1, suggesting superior phonon transport behaviors within these frequencies, which was traditionally considered impossible in disordered systems. Interestingly, we observe the emergence of a group of highly localized states associated withω∼200cm-1, which is characterized by a significant peak in the IPR and a peak in phonon density of states at the corresponding frequency, and is potentially useful for decoupling electron and phonon degrees of freedom. These unique properties of disordered SHU chains have implications in the design and engineering of novel quantum materials for thermal and phononic applications.

Confinement effect of inter-arm interactions on glass formation in star polymer melts

We utilized molecular dynamic simulation to investigate the glass formation of star polymer melts in which the topological complexity is varied by altering the number of star arms (f). Emphasis was placed on how the "confinement effect" of repulsive inter-arm interactions within star polymers influences the thermodynamics and dynamics of star polymer melts. All the characteristic temperatures of glass formation were found to progressively increase with increasing f, but unexpectedly the fragility parameter KVFT was found to decrease with increasing f. As previously observed, stars having more than 5 or 6 arms adopt an average particle-like structure that is more contracted relative to the linear polymer size having the same mass and exhibit a strong tendency for intermolecular and intramolecular segregation. We systematically analyzed how varying f alters collective particle motion, dynamic heterogeneity, the decoupling exponent ζ phenomenologically linking the slow β- and α-relaxation times, and the thermodynamic scaling index γt. Consistent with our hypothesis that the segmental dynamics of many-arm star melts and thin supported polymer films should exhibit similar trends arising from the common feature of high local segmental confinement, we found that ζ increases considerably with increasing f, as found in supported polymer films with decreasing thickness. Furthermore, increasing f led to greatly enhanced elastic heterogeneity, and this phenomenon correlates strongly with changes in ζ and γt. Our observations should be helpful in building a more rational theoretical framework for understanding how molecular topology and geometrical confinement influence the dynamics of glass-forming materials more broadly.

Designer pair statistics of disordered many-particle systems with novel properties

The knowledge of exact analytical functional forms for the pair correlation function g2(r) and its corresponding structure factor S(k) of disordered many-particle systems is limited. For fundamental and practical reasons, it is highly desirable to add to the existing database of analytical functional forms for such pair statistics. Here, we design a plethora of such pair functions in direct and Fourier spaces across the first three Euclidean space dimensions that are realizable by diverse many-particle systems with varying degrees of correlated disorder across length scales, spanning a wide spectrum of hyperuniform, typical nonhyperuniform, and antihyperuniform ones. This is accomplished by utilizing an efficient inverse algorithm that determines equilibrium states with up to pair interactions at positive temperatures that precisely match targeted forms for both g2(r) and S(k). Among other results, we realize an example with the strongest hyperuniform property among known positive-temperature equilibrium states, critical-point systems (implying unusual 1D systems with phase transitions) that are not in the Ising universality class, systems that attain self-similar pair statistics under Fourier transformation, and an experimentally feasible polymer model. We show that our pair functions enable one to achieve many-particle systems with a wide range of translational order and self-diffusion coefficients D, which are inversely related to one another. One can design other realizable pair statistics via linear combinations of our functions or by applying our inverse procedure to other desirable functional forms. Our approach facilitates the inverse design of materials with desirable physical and chemical properties by tuning their pair statistics.

Finite-size effects in the static structure factor S(k) and S(0) for a two-dimensional Yukawa liquid

Finite-size effects in the static structure factor S(k) are analyzed for an amorphous substance. As the number of particles is reduced, S(0) increases greatly, up to an order of magnitude. Meanwhile, there is a decrease in the height of the first peak S_{peak}. These finite-size effects are modeled accurately by the Binder formula for S(0) and our empirical formula for S_{peak}. Procedures are suggested to correct for finite-size effects in S(k) data and in the hyperuniformity index H≡S(0)/S_{peak}. These principles generally apply to S(k) obtained from particle positions in noncrystalline substances. The amorphous substance we simulate is a two-dimensional liquid, with a soft Yukawa interaction modeling a dusty plasma experiment.

Approach to hyperuniformity in a metallic glass-forming material exhibiting a fragile to strong glass transition

We investigate a metallic glass-forming (GF) material (Al₉₀Sm₁₀) exhibiting a fragile-strong (FS) glass-formation by molecular dynamics simulation to better understand this highly distinctive pattern of glass-formation in which many of the usual phenomenological relations describing relaxation times and diffusion of ordinary GF liquids no longer apply, and where instead genuine thermodynamic features are observed in response functions and little thermodynamic signature is exhibited at the glass transition temperature, T_g. Given the many unexpected similarities between the thermodynamics and dynamics of this metallic GF material with water, we first focus on the anomalous static scattering in this liquid, following recent studies on water, silicon and other FS GF liquids. We quantify the "hyperuniformity index" H of our liquid, which provides a quantitative measure of molecular "jamming". To gain insight into the T-dependence and magnitude of H, we also estimate another more familiar measure of particle localization, the Debye-Waller parameter 〈u²〉 describing the mean-square particle displacement on a timescale on the order of the fast relaxation time, and we also calculate H and 〈u²〉 for heated crystalline Cu. This comparative analysis between H and 〈u²〉 for crystalline and metallic glass materials allows us to understand the critical value of H on the order of 10^-3 as being analogous to the Lindemann criterion for both the melting of crystals and the "softening" of glasses. We further interpret the emergence of FS GF and liquid-liquid phase separation in this class of liquids to arise from a cooperative self-assembly process in the GF liquid.

Polymer-Grafted Nanoparticles with Variable Grafting Densities for High Energy Density Polymeric Nanocomposite Dielectric Capacitors

Designing high energy density dielectric capacitors for advanced energy storage systems needs nanocomposite-based dielectric materials, which can utilize the properties of both inorganic and polymeric materials. Polymer-grafted nanoparticle (PGNP)-based nanocomposites alleviate the problems of poor nanocomposite properties by providing synergistic control over nanoparticle and polymer properties. Here, we synthesize "core-shell" barium titanate-poly(methyl methacrylate) (BaTiO₃-PMMA) grafted PGNPs using surface-initiated atom transfer polymerization (SI-ATRP) with variable grafting densities of (0.303 to 0.929) chains/nm² and high molecular masses (97700 g/mL to 130000 g/mol) and observe that low grafted density and high molecular mass based PGNP show high permittivity, high dielectric strength, and hence higher energy densities (≈ 5.2 J/cm³) as compared to the higher grafted density PGNPs, presumably due to their "star-polymer"-like conformations with higher chain-end densities that are known to enhance breakdown. Nonetheless, these energy densities are an order of magnitude higher than their nanocomposite blend counterparts. We expect that these PGNPs can be readily used as commercial dielectric capacitors, and these findings can serve as guiding principles for developing tunable high energy density energy storage devices using PGNP systems.

Comparison of the static structure factor at long wavelengths for a dusty plasma liquid and other liquids

Especially small values of the static structure factor S(k) at long wavelengths, i.e., small k, were obtained in an analysis of experimental data, for a two-dimensional dusty plasma in its liquid state. For comparison, an analysis of S(k) data was carried out for many previously published experiments with other liquids. The latter analysis indicates that the magnitude of S(k) at small k is typically in a range 0.02-0.13. In contrast, the corresponding value for a dusty plasma liquid was found to be as small as 0.0139. Another basic finding for the dusty plasma liquid is that S(k) at small k generally increases with temperature, with its lowest value, noted above, occurring near the melting point. Simulations were carried out for the dusty plasma liquid, and their results are generally consistent with the experiment. Since a dusty plasma has a soft interparticle interaction, our findings support earlier theoretical suggestions that a useful design strategy for creating materials having exceptionally low values of S(0), so-called hyperuniform materials, is the use of a condensed material composed of particles that interact softly at their periphery.

Facile Entropy-Driven Segregation of Imprinted Polymer-Grafted Nanoparticle Brush Blends by Solvent Vapor Annealing Soft Lithography

Polymer-grafted nanoparticles (PGNPs) have attracted extensive research interest due to their potential for enhancing mechanical and electrical properties of both bulk polymer composite materials, as well as thin polymer films incorporating these nanoparticles (NPs). In previous studies, we have shown that an entropic driving force serves to organize low-molecular-mass PGNPs in imprinted blend films of PGNPs with low-molecular-mass homopolymers. In this work, we developed a novel solvent vapor annealing soft lithography (SVA-SL) method to overcome the technical difficulties in processing the high-molecular-mass PGNP blends due to the intrinsically sluggish melt annealing kinetics found in the phase separation of these blend PGNP materials. In particular, we utilized SVA-SL to create nanopatterns in blends of PGNPs having relatively high-molecular-mass-grafted layers but with cores of NPs having greatly different sizes. The minimization of the entropic free energy in the present system corresponded to larger PGNPs partitioning almost exclusively into the "mesa" regions of the imprinted PGNP blend films, as quantified by the estimation of the partition coefficient, K_p. The use of the SVA-SL processing method is important because it allows facile imprint patterning of PGNP materials and large-scale organization of the PGNPs even when the grafted chain lengths are long enough for the chains to be highly entangled, allowing enhanced thermo-mechanical property enhancements of the resulting films and a corresponding extended range of potential nanotech applications.

Hyperuniform jammed sphere packings have anomalous material properties

A spatial distribution is hyperuniform if it has local density fluctuations that vanish in the limit of long length scales. Hyperuniformity is a well known property of both crystals and quasicrystals. Of recent interest, however, is disordered hyperuniformity: the presence of hyperuniform scaling without long-range configurational order. Jammed granular packings have been proposed as an example of disordered hyperuniformity, but recent numerical investigation has revealed that many jammed systems instead exhibit a complex set of distinct behaviors at long, emergent length scales. We use the Voronoi tessellation as a tool to define a set of rescaling transformations that can impose hyperuniformity on an arbitrary weighted point process, and show that these transformations can be used in simulations to iteratively generate hyperuniform, mechanically stable packings of athermal soft spheres. These hyperuniform jammed packings display atypical mechanical properties, particularly in the low-frequency phononic excitations, which exhibit an isolated band of highly collective modes and a band gap around zero frequency.

Hyperuniform Active Chiral Fluids with Tunable Internal Structure

Large density fluctuations observed in active systems and hyperuniformity are two seemingly incompatible phenomena. However, the formation of hyperuniform states has been recently predicted in nonequilibrium fluids formed by chiral particles performing circular motion with the same handedness. Here we report evidence of hyperuniformity realized in a chiral active fluid comprised of pear-shaped Quincke rollers of arbitrary handedness. We show that hyperuniformity and large density fluctuations, triggered by dynamic clustering, coexist in this system at different length scales. The system loses its hyperuniformity as the curvature of particles' motion increases, transforming them into localized spinners. Our results experimentally demonstrate a novel hyperuniform active fluid and provide new insights into an interplay between chirality, activity, and hyperuniformity.

Ultraefficient reconstruction of effectively hyperuniform disordered biphase materials via non-Gaussian random fields

Disordered hyperuniform systems are statistically isotropic and possess no Bragg peaks like liquids and glasses, yet they suppress large-scale density fluctuations in a similar manner as in perfect crystals. The unique hyperuniform long-range order in these systems endow them with nearly optimal transport, electronic, and mechanical properties. The concept of hyperuniformity was originally introduced for many-particle systems and has subsequently been generalized to biphase heterogeneous materials such as porous media, composites, polymers, and biological tissues for unconventional property discovery. Existing methods for rendering realizations of disordered hyperuniform biphase materials reconstruction typically employ stochastic optimization such as the simulated annealing approach, which requires many iterations. Here, we propose an explicit ultraefficient method for reconstructing effectively hyperuniform biphase materials, based on the second-order non-Gaussian random fields where no additional tuning step or iteration is needed. Both the effectively hyperuniform microstructure and the latent material property field can be simultaneously generated in a single reconstruction. Moreover, our method can also incorporate hierarchical uncertainties in the heterogeneous materials, including both uncertainties in the disordered material microstructure and material property variation within each phase. The efficiency and feasibility of the proposed reconstruction method are demonstrated via a wide spectrum of examples spanning from isotropic to anisotropic, effectively hyperuniform to nonhyperuniform, and antihyperuniform systems. Our ultraefficient reconstruction method can be readily incorporated into material design, probabilistic analysis, optimization, and discovery of novel disordered hyperuniform heterogeneous materials.

Long Wavelength Thermal Density Fluctuations in Molecular and Polymer Glass-Forming Liquids: Experimental and Theoretical Analysis under Isobaric Conditions

We establish via an in-depth analysis of experimental data that the dimensionless compressibility (proportional to the dimensionless amplitude of long wavelength thermal density fluctuations) of one-component normal and supercooled liquids of chemically complex nonpolar and weakly polar molecules and polymers follows extremely well a surprisingly simple and general temperature dependence over an exceptionally wide range of pressures and temperatures. A theoretical basis for this behavior is shown to exist in the venerable van der Waals model and its more modern interpretations. Although associated hydrogen-bonding (and to a lesser degree strongly polar) liquids display modestly more complex behavior, rather simple temperature and pressure dependences are also discovered. A new approach to collapse the temperature- and pressure-dependent dimensionless compressibility data onto a master curve is formulated that differs from the empirical thermodynamic scaling approach. As a practical matter, we also find that the dimensionless compressibility scales well as an inverse power law with temperature with an exponent that is system dependent and decreases with pressure. At very high pressures and low temperatures, the thermal liquid behavior appears to approach (but not reach) a repulsion-dominated random close packing limit. All these findings are relevant to our recent theoretical work on the problem of activated relaxation and vitrification of supercooled molecular and polymeric liquids.

Control of Phase Morphology of Binary Polymer Grafted Nanoparticle Blend Films via Direct Immersion Annealing

While the phase separation of binary mixtures of chemically different polymer-grafted nanoparticles (PGNPs) is observed to superficially resemble conventional polymer blends, the presence of a "soft" polymer-grafted layer on the inorganic core of these nanoparticles qualitatively alters the phase separation kinetics of these "nanoblends" from the typical pattern of behavior seen in polymer blends and other simple fluids. We investigate this system using a direct immersion annealing method (DIA) that allows for a facile tuning of the PGNPs phase boundary, phase separation kinetics, and the ultimate scale of phase separation after a sufficient "aging" time. In particular, by switching the DIA solvent composition from a selective one (which increases the interaction parameter according to Timmerman's rule) to an overall good solvent for both PGNP components, we can achieve rapid switchability between phase-separated and homogeneous states. Despite a relatively low and non-classical power-law coarsening exponent, the overall phase separation process is completed on a time scale on the order of a few minutes. Moreover, the roughness of the PGNP blend film saturates at a scale that is proportional to the in-plane phase separation pattern scale, as observed in previous blend and block copolymer film studies. The relatively low magnitude of the coarsening exponent n is attributed to a suppression of hydrodynamic interactions between the PGNPs. The DIA method provides a significant opportunity to control the phase separation morphology of PGNP blends by solution processing, and this method is expected to be quite useful in creating advanced materials.

Structural characterization of many-particle systems on approach to hyperuniform states

The study of hyperuniform states of matter is an emerging multidisciplinary field, impinging on topics in the physical sciences, mathematics, and biology. The focus of this work is the exploration of quantitative descriptors that herald when a many-particle system in d-dimensional Euclidean space R^{d} approaches a hyperuniform state as a function of the relevant control parameter. We establish quantitative criteria to ascertain the extent of hyperuniform and nonhyperuniform distance-scaling regimes as well as the crossover point between them in terms of the "volume" coefficient A and "surface-area" coefficient B associated with the local number variance σ^{2}(R) for a spherical window of radius R. The larger the ratio B/A, the larger the hyperuniform scaling regime, which becomes of infinite extent in the limit B/A→∞. To complement the known direct-space representation of the coefficient B in terms of the total correlation function h(r), we derive its corresponding Fourier representation in terms of the structure factor S(k), which is especially useful when scattering information is available experimentally or theoretically. We also demonstrate that the free-volume theory of the pressure of equilibrium packings of identical hard spheres that approach a strictly jammed state either along the stable crystal or metastable disordered branch dictates that such end states be exactly hyperuniform. Using the ratio B/A, as well as other diagnostic measures of hyperuniformity, including the hyperuniformity index H and the direct-correlation function length scale ξ_{c}, we study three different exactly solvable models as a function of the relevant control parameter, either density or temperature, with end states that are perfectly hyperuniform. Specifically, we analyze equilibrium systems of hard rods and "sticky" hard-sphere systems in arbitrary space dimension d as a function of density. We also examine low-temperature excited states of many-particle systems interacting with "stealthy" long-ranged pair interactions as the temperature tends to zero, where the ground states are disordered, hyperuniform, and infinitely degenerate. We demonstrate that our various diagnostic hyperuniformity measures are positively correlated with one another. The same diagnostic measures can be used to detect the degree to which imperfections in nearly hyperuniform systems cause deviations from perfect hyperuniformity. Moreover, the capacity to identify hyperuniform scaling regimes should be particularly useful in analyzing experimentally or computationally generated samples that are necessarily of finite size.

Circular swimming motility and disordered hyperuniform state in an algae system

Active matter comprises individually driven units that convert locally stored energy into mechanical motion. Interactions between driven units lead to a variety of nonequilibrium collective phenomena in active matter. One of such phenomena is anomalously large density fluctuations, which have been observed in both experiments and theories. Here we show that, on the contrary, density fluctuations in active matter can also be greatly suppressed. Our experiments are carried out with marine algae ([Formula: see text]), which swim in circles at the air-liquid interfaces with two different eukaryotic flagella. Cell swimming generates fluid flow that leads to effective repulsions between cells in the far field. The long-range nature of such repulsive interactions suppresses density fluctuations and generates disordered hyperuniform states under a wide range of density conditions. Emergence of hyperuniformity and associated scaling exponent are quantitatively reproduced in a numerical model whose main ingredients are effective hydrodynamic interactions and uncorrelated random cell motion. Our results demonstrate the existence of disordered hyperuniform states in active matter and suggest the possibility of using hydrodynamic flow for self-assembly in active matter.

Hidden Order Beyond Hyperuniformity in Critical Absorbing States

Disordered hyperuniformity is a description of hidden correlations in point distributions revealed by an anomalous suppression in fluctuations of local density at various coarse-graining length scales. In the absorbing phase of models exhibiting an active-absorbing state transition, this suppression extends up to a hyperuniform length scale that diverges at the critical point. Here, we demonstrate the existence of additional many-body correlations beyond hyperuniformity. These correlations are hidden in the higher moments of the probability distribution of the local density and extend up to a longer length scale with a faster divergence than the hyperuniform length on approaching the critical point. Our results suggest that a hidden order beyond hyperuniformity may generically be present in complex disordered systems.

Stone-Wales defects preserve hyperuniformity in amorphous two-dimensional networks

Disordered hyperuniformity (DHU) is a recently discovered novel state of many-body systems that possesses vanishing normalized infinite-wavelength density fluctuations similar to a perfect crystal and an amorphous structure like a liquid or glass. Here, we discover a hyperuniformity-preserving topological transformation in two-dimensional (2D) network structures that involves continuous introduction of Stone-Wales (SW) defects. Specifically, the static structure factor [Formula: see text] of the resulting defected networks possesses the scaling [Formula: see text] for small wave number k, where [Formula: see text] monotonically decreases as the SW defect concentration p increases, reaches [Formula: see text] at [Formula: see text], and remains almost flat beyond this p. Our findings have important implications for amorphous 2D materials since the SW defects are well known to capture the salient feature of disorder in these materials. Verified by recently synthesized single-layer amorphous graphene, our network models reveal unique electronic transport mechanisms and mechanical behaviors associated with distinct classes of disorder in 2D materials.

Design of nearly perfect hyperuniform polymeric materials

Disordered hyperuniform materials are exotic amorphous systems that simultaneously exhibit anomalous suppression of long-range density fluctuations, comparable in amplitude to that of crystals and quasi-crystalline materials, while lacking the translational order characteristic of simple liquids. We establish a framework to quantitatively predict the emergence of hyperuniformity in polymeric materials by considering the distribution of localized polymer subregions, instead of considering the whole material. We demonstrate that this highly tunable approach results in arbitrarily small long-range density fluctuations in the liquid state. Our simulations also indicate that long-ranged density fluctuation of the whole polymeric material is remarkably insensitive to molecular topology (linear chain, unknotted ring, star, and bottlebrush) and depends on temperature in an apparently near universal fashion. Our findings open the way for the creation of nearly perfect hyperuniform polymeric materials.

Optimized Large Hyperuniform Binary Colloidal Suspensions in Two Dimensions

The creation of disordered hyperuniform materials with extraordinary optical properties (e.g., large complete photonic band gaps) requires a capacity to synthesize large samples that are effectively hyperuniform down to the nanoscale. Motivated by this challenge, we propose a feasible equilibrium fabrication protocol using binary paramagnetic colloidal particles confined in a 2D plane. The strong and long-ranged dipolar interaction induced by a tunable magnetic field is free from screening effects that attenuate long-ranged electrostatic interactions in charged colloidal systems. Specifically, we numerically find a family of optimal size ratios that makes the two-phase system effectively hyperuniform. We show that hyperuniformity is a general consequence of low isothermal compressibilities, which makes our protocol suitable to treat more general systems with other long-ranged interactions, dimensionalities, and/or polydispersity. Our methodology paves the way to synthesize large photonic hyperuniform materials that function in the visible to infrared range and hence may accelerate the discovery of novel photonic materials.

Influence of Branching on the Configurational and Dynamical Properties of Entangled Polymer Melts

We probe the influence of branching on the configurational, packing, and density correlation function properties of polymer melts of linear and star polymers, with emphasis on molecular masses larger than the entanglement molecular mass of linear chains. In particular, we calculate the conformational properties of these polymers, such as the hydrodynamic radius R h , packing length p, pair correlation function g ( r ) , and polymer center of mass self-diffusion coefficient, D, with the use of coarse-grained molecular dynamics simulations. Our simulation results reproduce the phenomenology of simulated linear and branched polymers, and we attempt to understand our observations based on a combination of hydrodynamic and thermodynamic modeling. We introduce a model of "entanglement" phenomenon in high molecular mass polymers that assumes polymers can viewed in a coarse-grained sense as "soft" particles and, correspondingly, we model the emergence of heterogeneous dynamics in polymeric glass-forming liquids to occur in a fashion similar to glass-forming liquids in which the molecules have soft repulsive interactions. Based on this novel perspective of polymer melt dynamics, we propose a functional form for D that can describe our simulation results for both star and linear polymers, covering both the unentangled to entangled polymer melt regimes.

Methodology to construct large realizations of perfectly hyperuniform disordered packings

Disordered hyperuniform packings (or dispersions) are unusual amorphous two-phase materials that are endowed with exotic physical properties. Such hyperuniform systems are characterized by an anomalous suppression of volume-fraction fluctuations at infinitely long-wavelengths, compared to ordinary disordered materials. While there has been growing interest in such singular states of amorphous matter, a major obstacle has been an inability to produce large samples that are perfectly hyperuniform due to practical limitations of conventional numerical and experimental methods. To overcome these limitations, we introduce a general theoretical methodology to construct perfectly hyperuniform packings in d-dimensional Euclidean space R^{d}. Specifically, beginning with an initial general tessellation of space by disjoint cells that meets a "bounded-cell" condition, hard particles of general shape are placed inside each cell such that the local-cell particle packing fractions are identical to the global packing fraction. We prove that the constructed packings with a polydispersity in size are perfectly hyperuniform in the infinite-sample-size limit, regardless of particle shapes, positions, and numbers per cell. We use this theoretical formulation to devise an efficient and tunable algorithm to generate extremely large realizations of such packings. We employ two distinct initial tessellations: Voronoi as well as sphere tessellations. Beginning with Voronoi tessellations, we show that our algorithm can remarkably convert extremely large nonhyperuniform packings into hyperuniform ones in R^{2} and R^{3}. Implementing our theoretical methodology on sphere tessellations, we establish the hyperuniformity of the classical Hashin-Shtrikman multiscale coated-spheres structures, which are known to be two-phase media microstructures that possess optimal effective transport and elastic properties. A consequence of our work is a rigorous demonstration that packings that have identical tessellations can either be nonhyperuniform or hyperuniform by simply tuning local characteristics. It is noteworthy that our computationally designed hyperuniform two-phase systems can easily be fabricated via state-of-the-art methods, such as 2D photolithographic and 3D printing technologies. In addition, the tunability of our methodology offers a route for the discovery of novel disordered hyperuniform two-phase materials.