Programmable self-assembly enables the construction of complex molecular, supramolecular, and crystalline architectures from well-designed building blocks. We introduce a hypergraph-based formalism, Blocks & Bonds (B&B), which generalizes classical chemical graph theory by incorporating directed and multicolored interactions, internal symmetries, and hierarchical organization. Within this framework, we develop the Structure Code (SC), a compact and versatile language for describing self-assembled architectures. We define a Kolmogorov-style structural complexity as the total information content of SC, obtained through its tokenization and Shannon information assignment. Complementing this encoding-based measure, we introduce a much simpler quantity, the compositional complexity, which depends only on the number and cumulative usage of block and bond types in the construction set. A central result of this work is a strong empirical correlation between the token-based structural complexity and the compositional complexity across all examined systems. Owing to this agreement, the compositional complexity emerges as the most practical and broadly applicable measure: it is easy to compute, requires no explicit encoding, and yet closely tracks the actual information content of structurally diverse architectures. Applications to molecular systems (ethylene glycol and glucose), DNA-origami lattices, and crystalline assemblies show that B&B hypergraphs provide a unified, scalable, and information-efficient representation of structural organization, naturally capturing symmetry, modularity, and stereochemistry. This framework establishes a quantitative foundation for complexity-aware classification and inverse design of programmable matter.
Linker-mediated aggregation plays an important role in modern nanoscience. We demonstrate that it departs sharply from classical Smoluchowski kinetics because cluster reactivity evolves during growth. Combining theory with DNA-linked gold-nanoparticle experiments, we establish kinetic critical point controlled by linker abundance. Below threshold, active linkers are depleted and growth arrests; above threshold, clusters accumulate reactive sites, self-accelerate, and cross over to diffusion-limited coarsening. Experiments verify the predicted arrest, accelerated growth, and scaling collapse.
We use in situ transmission electron microscopy with automated tracking to study supported gold nanoparticles (NPs) during high-temperature vacuum annealing. The average mass loss per NP is governed by a flat, nearly size-independent substrate-mediated evaporation profile. On top of this mean shrinkage, individual NPs show significant fluctuations in apparent growth or shrinkage, and NP volume follows a random-walk-like trajectory. To rationalize both the ensemble-mean behavior and the particle-resolved variability, we develop a self-consistent theory that couples substrate-mediated evaporation to collective 2D Ostwald-type mass exchange through a shared adatom field, described in terms of a renormalized screening length and background concentration. In the experimentally relevant regime, the theory predicts an approximately size-independent mean shrinkage rate and clarifies how net mass loss suppresses classical coarsening. Superimposed on this deterministic drift, we quantify stochastic volume trajectories and capture their fluctuation spectrum with a minimal Langevin description consistent with intermittent adatom attachment and detachment events. In addition, we characterize the lateral diffusive motion of NPs, which is responsible for their coalescence. Altogether, our results highlight that stochasticity is intrinsic at the nanoscale and that predicting the evolution of supported NPs at early and intermediate times requires a unified framework combining substrate-mediated evaporation, collective mass exchange, and stochastic fluctuations.
We present a framework for chemical learning based on Competitive Dimerization Networks (CDNs) - systems in which multiple molecular species, e.g., proteins, DNA oligomers, or RNA oligomers, reversibly bind to form dimers. We show numerically that these networks can, in principle, be trained in vitro through directed evolution, enabling the implementation of complex learning tasks such as multiclass classification without digital hardware or prior knowledge of all microscopic association constants. Each molecular species functions analogously to a neuron, with binding affinities acting as tunable synaptic weights. A training protocol involving mutation, selection, and amplification of DNA-based components allows CDNs to robustly discriminate among noisy input patterns. The resulting classifiers exhibit strong output contrast and high mutual information between input and output, especially when guided by a contrast-enhancing loss function. Comparative analysis with in silico gradient descent training reveals closely correlated performance. These results establish CDNs as a promising platform for analog physical computation, bridging synthetic biology and machine learning, and advancing the development of adaptive, energy-efficient molecular computing systems.
Nucleation and growth are studied in a system that undergoes diffusion-controlled condensation under gradual changes in parameters, such as cooling. It is demonstrated that when the Gibbs-Thompson effect becomes negligible, the system falls into a universal regime. i.e., the final droplet size distribution remains invariant under certain rescaling of system parameters. An approximate yet very accurate analytic form is obtained for the final droplet size distribution in this regime.
Recent advances in DNA nanotechnology allow for the assembly of nanocomponents with nanoscale precision, leading to the emergence of DNA-based material fabrication approaches. Yet, transferring these nano- and micron-scale structural arrangements to the macroscale morphologies remains a challenge, which limits the development of materials and devices based on DNA nanotechnology. Here, we demonstrate a materials fabrication approach that combines DNA-programmable assembly with actively driven processes controlled by acoustic fields. This combination provides a prescribed nanoscale order, as dictated by equilibrium assembly through DNA-encoded interactions, and field-shaped macroscale morphology, as regulated by out-of-equilibrium materials formation through specific acoustic stimulation. Using optical and electron microscopy imaging and x-ray scattering, we further revealed the nucleation processes, domain fusion, and crystal growth under different acoustically stimulated conditions. The developed approach provides a pathway for the fabrication of complexly shaped macroscale morphologies for DNA-programmable nanomaterials by controlling spatiotemporal characteristics of the acoustic fields. DNA nanotechnology is useful in preparing nano- and meso-components, but transfer to macroscale arrangements is challenging. Here, the authors report an assembly approach combining DNA programmable assembly with process controlled by acoustic fields to prepare macroscale morphologies.
We explore the potential of nanocrystals (a term used equivalently to nanoparticles) as building blocks for nanomaterials, and the current advances and open challenges for fundamental science developments and applications. Nanocrystal assemblies are inherently multiscale, and the generation of revolutionary material properties requires a precise understanding of the relationship between structure and function, the former being determined by classical effects and the latter often by quantum effects. With an emphasis on theory and computation, we discuss challenges that hamper current assembly strategies and to what extent nanocrystal assemblies represent thermodynamic equilibrium or kinetically trapped metastable states. We also examine dynamic effects and optimization of assembly protocols. Finally, we discuss promising material functions and examples of their realization with nanocrystal assemblies.
Life as we know it relies on the interplay between catalytic activity and information processing carried out by biological polymers. Here we present a plausible pathway by which a pool of prebiotic information-coding oligomers could acquire an early catalytic function, namely sequence- specific cleavage activity. Starting with a system capable of non-enzymatic templated replication, we demonstrate that even non-catalyzed spontaneous cleavage would promote proliferation by generating short fragments that act as primers. Furthermore, we show that catalytic cleavage function can naturally emerge and proliferate in this system. Specifically, a cooperative catalytic network with four subpopulations of oligomers is selected by the evolution in competition with chains lacking catalytic activity. The cooperative system emerges through the functional differentiation of oligomers into catalysts and their substrates. The model is inspired by the structure of the hammerhead RNA enzyme as well as other DNA- and RNA-based enzymes with cleavage activity that readily emerge through natural or artificial selection. We identify the conditions necessary for the emergence of the cooperative catalytic network. In particular, we show that it requires the catalytic rate enhancement over the spontaneous cleavage rate to be at least 10 2 − 10 3 , a factor consistent with the existing experiments. The evolutionary pressure leads to a further increase in catalytic efficiency. The presented mechanism provides an escape route from a relatively simple pairwise replication of oligomers towards a more complex behavior involving catalytic function. This provides a bridge between the information-first origin of life scenarios and the paradigm of autocatalytic sets and hypercycles, albeit based on cleavage rather than synthesis of reactants.
The study of the physico-chemical properties and fatty acid composition of lard lipids of pigs of different weight conditions is of important practical importance in pig breeding for establishing the most optimal terms for slaughtering animals in order to obtain high-quality meat and lard pork. A high concentration of saturated and monounsaturated fatty acids in triglycerides is closely related to their active synthesis and accumulation in the body, even when pigs are kept on low-fat diets The aim of the research was to study the chemical, fatty acid composition and physicochemical properties of pork lard of different breeds. Physico-chemical properties of lard in pigs of the mentioned genotypes depended on breed and sex. There is an inverse correlation between the melting temperature and the iodine number of lard. With an increase in the iodine number, the unsaturation of fats increases, the melting point of lard decreases. The presence of a larger amount of saturated fatty acids in the lard of the Poltava meat breed caused a higher melting point – 42.9. The highest protein was observed in pigs of the meat production direction – 2.68. The smallest amount of protein was recorded in pigs of the Myrhorod breed – 1.98. An inverse-proportional relationship between the content of protein and actual fat in lard was observed. A greater amount of fat caused a lower protein content.
Bond-orientational order in DNA-assembled nanoparticles lattices is explored with the help of recently introduced Symmetry-specific Bond Order Parameters (SymBOPs). This approach provides a more sensitive analysis of local order than traditional scalar BOPs, facilitating the identification of coherent domains at the single bond level. The present study expands the method initially developed for assemblies of anisotropic particles to the isotropic ones or cases where particle orientation information is unavailable. The SymBOP analysis was applied to experiments on DNA-frame-based assembly of nanoparticle lattices. It proved highly sensitive in identifying coherent crystalline domains with different orientations, as well as detecting topological defects, such as dislocations. Furthermore, the analysis distinguishes individual sublattices within a single crystalline domain, such as pair of interpenetrating FCC lattices within a cubic diamond. The results underscore the versatility and robustness of SymBOPs in characterizing ordering phenomena, making them valuable tools for investigating structural properties in various systems.
Even a century after the formulation of Quantum Mechanics (QM), the wave function collapse (WFC) remains a contentious aspect of the theory. Environment-induced decoherence has offered a partial resolution by illustrating how unitary evolution in an open quantum system can lead to effective WFC within its components. However, this approach by itself does not lead to a fully self-consistent reformulation of QM. We introduce a modified set of QM postulates, which exclude both WFC and Born's probability rule. They are replaced with the Maximum Entropy Principle, a weaker postulate that specifies conditional probabilities for mutually compatible observations. Within this formulation, both WFC and Born's rule become emerging properties.
A comprehensive framework of characterizing complex self-assembled structures with a set of orientational order parameters is presented. It is especially relevant in the context of using anisotropic building blocks with various symmetries. Two classes of tensor order parameters are associated with polyhedral nematic and bond orientational order. For the latter, a variation of classical bond order parameters (BOPs) is introduced, which takes advantage of the symmetry of constituent particles and/or expected crystalline phases. These symmetrized BOPs can be averaged over an entire system or assigned locally to an individual bond. By combining that with the bond percolation procedure, one is able to identify coherent domains within a self-assembled structure. As a demonstration of the proposed framework, we apply it to a simulated hybrid system that combines isotropic and patchy particles with octahedral symmetry. Not only does the methodology allow one to identify individual crystalline domains but also it detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.
ABSTRACT Celestial bodies approximated with rigid triaxial ellipsoids in a two-body system can rotate chaotically due to the time-varying gravitational torque from the central mass. At small orbital eccentricity values, rotation is short-term orderly and predictable within the commensurate spin–orbit resonances, while at eccentricity approaching unity, chaos completely takes over. Here, we present the full three-dimensional rotational equations of motion around all three principal axes for triaxial minor planets and two independent methods of numerical solution based on Euler rotations and quaternion algebra. The domains of chaotic rotation are numerically investigated over the entire range of eccentricity with a combination of trial integrations of Euler’s equations of motion and the GALI(k) (Generalized Alignment Index) method. We quantify the dependence of the order–chaos boundaries on shape by changing a prolateness parameter, and find that the main 1:1 spin–orbit resonance disappears for specific moderately prolate shapes already at eccentricities as low as 0.3. The island of short-term stability around the main 1:1 resonance shrinks with increasing eccentricity at a fixed low degree of prolateness and completely vanishes at approximately 0.8. This island is also encroached by chaos on longer time-scales, indicating longer Lyapunov exponents. Trajectories in the close vicinity of the 3:2 spin–orbit resonance become chaotic at smaller eccentricities, but separated enclaves of orderly rotation emerge at eccentricities as high as 0.8. Initial perturbations of rotational velocity in latitude away from the exact equilibrium result in a spectrum of free libration, nutation, and polar wander, which is not well matched by the linearized analysis omitting the inertial terms.
Brownian dynamics is used to study self-assembly in a hybrid system of isotropic particles (IPs), combined with anisotropic building blocks that represent special “designer particles.” Those are modeled as spherical patchy particles (PPs) with binding only allowed between their patches and IPs. In this study, two types of PPs are considered: Octahedral PPs (Oh-PPs) and Square PPs (Sq-PPs), with octahedral and square arrangements of patches, respectively. The self-assembly is additionally facilitated by the simulated annealing procedure. The resultant structures are characterized by a combination of local correlations in cubatic ordering and a symmetry-specific variation of bond orientation order parameters (SymBOPs). By varying the PP/IP size ratio, we detected a sharp crossover between two distinct morphologies in both types of systems. High symmetry phases, NaCl crystal for Oh-PP and square lattice for Sq-PP, are observed for larger size ratios. For the smaller ones, the dominant morphologies are significantly different, e.g., Oh-PPs form a compact amorphous structure with predominantly face-to-face orientation of neighboring PPs. Unusually, for a morphology without a long-range order, it is still possible to identify well organized coherent clusters of this structure, thanks to the adoption of our SymBOP-based characterization.
We present a statistical mechanical description of randomly packed spherical particles, where the average coordination number is treated as a macroscopic thermodynamic variable. The overall packing entropy is shown to have two contributions: geometric, reflecting statistical weights of individual configurations, and topological, which corresponds to the number of topologically distinct states. Both of them are computed in the thermodynamic limit for isostatic and weakly underconstrained packings in 2D and 3D. The theory generalizes concepts of granular and glassy configurational entropies for the case of nonjammed systems. It is directly applicable to sticky colloids and predicts an asymptotic phase behavior of sticky spheres in the limit of strong binding.
In Fall 2020, universities saw extensive transmission of SARS-CoV-2 among their populations, threatening health of the university and surrounding communities, and viability of in-person instruction. Here we report a case study at the University of Illinois at Urbana-Champaign, where a multimodal "SHIELD: Target, Test, and Tell" program, with other non-pharmaceutical interventions, was employed to keep classrooms and laboratories open. The program included epidemiological modeling and surveillance, fast/frequent testing using a novel low-cost and scalable saliva-based RT-qPCR assay for SARS-CoV-2 that bypasses RNA extraction, called covidSHIELD, and digital tools for communication and compliance. In Fall 2020, we performed >1,000,000 covidSHIELD tests, positivity rates remained low, we had zero COVID-19-related hospitalizations or deaths amongst our university community, and mortality in the surrounding Champaign County was reduced more than 4-fold relative to expected. This case study shows that fast/frequent testing and other interventions mitigated transmission of SARS-CoV-2 at a large public university.
Epidemics generally spread through a succession of waves that reflect factors on multiple timescales. On short timescales, super-spreading events lead to burstiness and overdispersion, while long-term persistent heterogeneity in susceptibility is expected to lead to a reduction in the infection peak and the herd immunity threshold (HIT). Here, we develop a general approach to encompass both timescales, including time variations in individual social activity, and demonstrate how to incorporate them phenomenologically into a wide class of epidemiological models through parameterization. We derive a non-linear dependence of the effective reproduction number Re on the susceptible population fraction S. We show that a state of transient collective immunity (TCI) emerges well below the HIT during early, high-paced stages of the epidemic. However, this is a fragile state that wanes over time due to changing levels of social activity, and so the infection peak is not an indication of herd immunity: subsequent waves can and will emerge due to behavioral changes in the population, driven (e.g.) by seasonal factors. Transient and long-term levels of heterogeneity are estimated by using empirical data from the COVID-19 epidemic as well as from real-life face-to-face contact networks. These results suggest that the hardest-hit areas, such as NYC, have achieved TCI following the first wave of the epidemic, but likely remain below the long-term HIT. Thus, in contrast to some previous claims, these regions can still experience subsequent waves.
It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, i.e. constant in time, here we discuss the consequences of dynamic heterogeneity. We integrate the stochastic dynamics of social activity into traditional epidemiological models. The overall epidemic dynamics is condensed to three differential equations and is characterized by an emergent long timescale. Our model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves. Individual waves are suppressed at the state of Transient Collective Immunity (TCI), which subsequently degrades due to the dynamic nature of social activity. Our results also provide a new mechanism for emerging pathogens to transition from a fast-paced epidemic to the endemic state.
It is well recognized that population heterogeneity plays an important role in the spread of epidemics. While individual variations in social activity are often assumed to be persistent, that is, constant in time, here we discuss the consequences of dynamic heterogeneity. By integrating the stochastic dynamics of social activity into traditional epidemiological models, we demonstrate the emergence of a new long timescale governing the epidemic, in broad agreement with empirical data. Our stochastic social activity model captures multiple features of real-life epidemics such as COVID-19, including prolonged plateaus and multiple waves, which are transiently suppressed due to the dynamic nature of social activity. The existence of a long timescale due to the interplay between epidemic and social dynamics provides a unifying picture of how a fast-paced epidemic typically will transition to an endemic state.