Source author record

Alexei V. Tkachenko

Alexei V. Tkachenko appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
11topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

14 published item(s)

preprint2022arXiv

Chaos over Order: Mapping 3D Rotation of Triaxial Asteroids and Minor Planets

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 3D rotational equations of motion around all three principle 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$) 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.

preprint2021arXiv

Controlling morphology in hybrid isotropic/patchy particle assemblies

Brownian Dynamics is used to study self-assembly in a hybrid system of istotropic 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 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.

preprint2021arXiv

Geometric and Topological Entropies of Sphere Packing

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 packings in 2D and 3D, and the result is further expanded to the case of "floppy" particle clusters. The theory is directly applicable to sticky colloids, and in addition, generalizes concepts of granular and glassy configurational entropies for the case of non-jammed systems.

preprint2021arXiv

Stochastic social behavior coupled to COVID-19 dynamics leads to waves, plateaus and an 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, i.e. 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 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 the 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 the endemic state.

preprint2021arXiv

Symmetry-specific orientational order parameters for complex structures

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, respectively. For the latter, a variation of classical bond order parameters is introduced that takes advantage of the symmetry of constituent particles, and/or expected crystalline phases. These Symmetrized Bond Order Parameters (SymBOPs) can be averaged over an entire system, or assigned locally to an individual bond. By combining that with 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 it also detects coherent clusters of a peculiar compact amorphous structure that is not space-filling and lacks any long-range order.

preprint2021arXiv

Time-dependent heterogeneity leads to transient suppression of the COVID-19 epidemic, not herd immunity

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.

preprint2020arXiv

Ligation of random oligomers leads to emergence of autocatalytic sequence network

The emergence of longer information-carrying and functional nucleotide polymers from random short strands was a major stepping stone at the dawn of life. But the formation of those polymers under temperature oscillation required some form of selection. A plausible mechanism is template-based ligation where theoretical work already suggested a reduction in information entropy. Here, we show how nontrivial sequence patterns emerge in a system of random 12mer DNA sequences subject to enzyme-based templated ligation reaction and temperature cycling. The strands acted both as a template and substrates of the reaction and thereby formed longer oligomers. The selection for templating sequences leads to the development of a multiscale ligation landscape. A position-dependent sequence pattern emerged with a segregation into mutually complementary pools of A-rich and T-rich sequences. Even without selection for function, the base pairing of DNA with ligation showed a dynamics resembling Darwinian evolution.

preprint2015arXiv

Chromatic patchy particles: effects of specific interactions on liquid structure

We study the structural and thermodynamic properties of patchy particle liquids, with a special focus on the role of "color", i.e. specific interactions between individual patches. A possible experimental realization of such "chromatic" interactions is by decorating the particle patches with single-stranded DNA linkers. The complementarity of the linkers can promote selective bond formation between predetermined pairs of patches. By using MD simulations, we compare the local connectivity, the bond orientation order, and other structural properties of the aggregates formed by the "colored" and "colorless" systems. The analysis is done for spherical particles with two different patch arrangements (tetrahedral and cubic). It is found that the aggregated (liquid) phase of the "colorless" patchy particles is better connected, denser and typically has stronger local order than the corresponding "colored" one. This, in turn, makes the colored liquid less stable thermodynamically. Specifically, we predict that in a typical case the chromatic interactions should increase the relative stability of the crystalline phase with respect to the disordered liquid, thus expanding its region in the phase diagram.

preprint2015arXiv

Onset of autocatalysis of information-coding polymers

Self-replicating systems based on information-coding polymers are of crucial importance in biology. They also recently emerged as a paradigm in material design on nano- and micro-scales. We present a general theoretical and numerical analysis of the problem of spontaneous emergence of autocatalysis for heteropolymers capable of template-assisted ligation driven by cyclic changes in the environment. Our central result is the existence of the first order transition between the regime dominated by free monomers and that with a self-sustaining population of sufficiently long chains. We provide a simple, mathematically tractable model supported by numerical simulations, which predicts the distribution of chain lengths and the onset of autocatalysis in terms of the overall monomer concentration and two fundamental rate constants. Another key result of our study is the emergence of the kinetically-limited optimal overlap length between a template and each of its two substrates. The template-assisted ligation allows for heritable transmission of the information encoded in chain sequences thus opening up the possibility of long-term memory and evolvability in such systems.

preprint2013arXiv

A simple analytical formula for the free-energy of ligand-receptor mediated interactions

Recently \1, we presented a general theory for calculat- ing the strength and properties of colloidal interactions mediated by ligand-receptor bonds (such as those that bind DNA-coated colloids). In this communication, we derive a surprisingly simple analytical form for the inter- action free energy, which was previously obtainable only via a costly numerical thermodynamic integration. As a result, the computational effort to obtain potentials of in- teraction is significantly reduced. Moreover, we can gain insight from this analytic expression for the free energy in limiting cases. In particular, the connection of our general theory to other previous specialised approaches is now made transparent. This important simplification will significantly broaden the scope of our theory.

preprint2013arXiv

DNA-Programmed Mesoscopic Architecture

We study the problem of the self-assembly of nanoparticles (NPs) into finite mesoscopic structures with a programmed local morphology and complex overall shape. Our proposed building blocks are NPs directionally-functionalized with DNA. The combination of directionality and selectivity of interactions allows one to avoid unwanted metastable configurations which have been shown to lead to slow self-assembly kinetics even in much simpler systems. With numerical simulations, we show that a variety of target mesoscopic objects can be designed and self-assembled in near perfect yield. They include cubes, pyramids, boxes and even an Empire State Building model. We summarize our findings with a set of design strategies that lead to the successful self-assembly of a wide range of mesostructures.

preprint2013arXiv

Precursors of order in aggregates of patchy particles

We study computationally the local structure of aggregated systems of patchy particles. By calculating the probability distribution functions of various rotational invariants we can identify the precursors of orientation order in amorphous phase. Surprisingly, the strongest signature of local order is observed for 4-patch particles with tetrahedral symmetry, not for 6-patch particles with the cubic one. This trend is exactly opposite to their known ability to crystallize. We relate this anomaly to the observation that a generic aggregate of patchy systems has coordination number close to 4. Our results also suggest a significant correlation between rotational order in the studied liquids with the corresponding crystalline phases, making this approach potentially useful for a broader range of patchy systems.

preprint2009arXiv

A Generalized Theory of DNA Looping and Cyclization

We have developed a generalized semi-analytic approach for efficiently computing cyclization and looping $J$ factors of DNA under arbitrary binding constraints. Many biological systems involving DNA-protein interactions impose precise boundary conditions on DNA, which necessitates a treatment beyond the Shimada-Yamakawa model for ring cyclization. Our model allows for DNA to be treated as a heteropolymer with sequence-dependent intrinsic curvature and stiffness. In this framework, we independently compute enthlapic and entropic contributions to the $J$ factor and show that even at small length scales $(\sim \ell_{p})$ entropic effects are significant. We propose a simple analytic formula to describe our numerical results for a homogenous DNA in planar loops, which can be used to predict experimental cyclization and loop formation rates as a function of loop size and binding geometry. We also introduce an effective torsional persistence length that describes the coupling between twist and bending of DNA when looped.

preprint2007arXiv

Effects of Sequence Disorder on DNA Looping and Cyclization

Effects of sequence disorder on looping and cyclization of the double-stranded DNA are studied theoretically. Both random intrinsic curvature and inhomogeneous bending rigidity are found to result in a remarkably wide distribution of cyclization probabilities. For short DNA segments, the range of the distribution reaches several orders of magnitude for even completely random sequences. The ensemble averaged values of the cyclization probability are also calculated, and the connection to the recent experiments is discussed.