Unraveling the Influence of Topology and Spatial Confinement on Equilibrium and Relaxation Properties of Interlocked Ring Polymers.

We use Langevin dynamics simulations to study linked ring polymers in channel confinement. We address the in- and out-of-equilibrium behavior of the systems for varying degrees of confinement and increasing topological and geometrical complexity of the interlocking. The main findings are three. First, metric observables of different link topologies collapse onto the same master curve when plotted against the crossing number, revealing a universal response to confinement. Second, the relaxation process from initially stretched states is faster for more complex links. We ascribe these properties to the interplay of several effects, including the dependence of topological friction on the link complexity. Finally, we show that transient forms of geometrical entanglement purposely added to the initial stressed state can leave distinctive signatures in force-spectroscopy curves. The insight provided by the findings could be leveraged in single-molecule nanochannel experiments to identify geometric entanglement within topologically linked rings.

Learning how to find targets in the micro-world: the case of intermittent active Brownian particles.

Finding the best strategy to minimize the time needed to find a given target is a crucial task both in nature and in reaching decisive technological advances. By considering learning agents able to switch their dynamics between standard and active Brownian motion, here we focus on developing effective target-search behavioral policies for microswimmers navigating a homogeneous environment and searching for targets of unknown position. We exploit projective simulation, a reinforcement learning algorithm, to acquire an efficient stochastic policy represented by the probability of switching the phase, i.e. the navigation mode, in response to the type and the duration of the current phase. Our findings reveal that the target-search efficiency increases with the particle's self-propulsion during the active phase and that, while the optimal duration of the passive case decreases monotonically with the activity, the optimal duration of the active phase displays a non-monotonic behavior.

Active Brownian particles in a circular disk with an absorbing boundary.

We solve the time-dependent Fokker-Planck equation for a two-dimensional active Brownian particle exploring a circular region with an absorbing boundary. Using the passive Brownian particle as basis states and dealing with the activity as a perturbation, we provide a matrix representation of the Fokker-Planck operator and we express the propagator in terms of the perturbed eigenvalues and eigenfunctions. Alternatively, we show that the propagator can be expressed as a combination of the equilibrium eigenstates with weights depending only on time and on the initial conditions, and obeying exact iterative relations. Our solution allows also obtaining the survival probability and the first-passage time distribution. These latter quantities exhibit peculiarities induced by the nonequilibrium character of the dynamics; in particular, they display a strong dependence on the activity of the particle and, to a less extent, also on its rotational diffusivity.

Survival strategies of artificial active agents.

Artificial cells can be engineered to display dynamics sharing remarkable features in common with the survival behavior of living organisms. In particular, such active systems can respond to stimuli provided by the environment and undertake specific displacements to remain out of equilibrium, e.g. by moving towards regions with higher fuel concentration. In spite of the intense experimental activity aiming at investigating this fascinating behavior, a rigorous definition and characterization of such "survival strategies" from a statistical physics perspective is still missing. In this work, we take a first step in this direction by adapting and applying to active systems the theoretical framework of Transition Path Theory, which was originally introduced to investigate rare thermally activated transitions in passive systems. We perform experiments on camphor disks navigating Petri dishes and perform simulations in the paradigmatic active Brownian particle model to show how the notions of transition probability density and committor function provide the pivotal concepts to identify survival strategies, improve modeling, and obtain and validate experimentally testable predictions. The definition of survival in these artificial systems paves the way to move beyond simple observation and to formally characterize, design and predict complex life-like behaviors.

Analytic Solution of an Active Brownian Particle in a Harmonic Well.

We provide an analytical solution for the time-dependent Fokker-Planck equation for a two-dimensional active Brownian particle trapped in an isotropic harmonic potential. Using the passive Brownian particle as basis states we show that the Fokker-Planck operator becomes lower diagonal, implying that the eigenvalues are unaffected by the activity. The propagator is then expressed as a combination of the equilibrium eigenstates with weights obeying exact iterative relations. We show that for the low-order correlation functions, such as the positional autocorrelation function, the recursion terminates at finite order in the Péclet number, allowing us to generate exact compact expressions and derive the velocity autocorrelation function and the time-dependent diffusion coefficient. The nonmonotonic behavior of latter quantities serves as a fingerprint of the nonequilibrium dynamics.

Geometric Predictors of Knotted and Linked Arcs.

Inspired by how certain proteins "sense" knots and entanglements in DNA molecules, here, we ask if local geometric features that may be used as a readout of the underlying topology of generic polymers exist. We perform molecular simulations of knotted and linked semiflexible polymers and study four geometric measures to predict topological entanglements: local curvature, local density, local 1D writhe, and nonlocal 3D writhe. We discover that local curvature is a poor predictor of entanglements. In contrast, segments with maximum local density or writhe correlate as much as 90% of the time with the shortest knotted and linked arcs. We find that this accuracy is preserved across different knot types and also under significant spherical confinement, which is known to delocalize essential crossings in knotted polymers. We further discover that nonlocal 3D writhe is the best geometric readout of the knot location. Finally, we discuss how these geometric features may be used to computationally analyze entanglements in generic polymer melts and gels.

Topological Friction and Relaxation Dynamics of Spatially Confined Catenated Polymers.

We study catenated ring polymers confined inside channels and slits with Langevin dynamics simulations and address how the contour position and size of the interlocked or physically linked region evolve with time. We show that the catenation constraints generate a drag, or topological friction, that couples the contour motion of the interlocked regions. Notably, the coupling strength decreases as the interlocking is made tighter, but also shorter, by confinement. Though the coupling strength differs for channel and slit confinement, the data outline a single universal curve when plotted against the size of the linked region. Finally, we study how the relaxation kinetics changes after one of the rings is cut open and conclude that considering interlocked circular polymers is key for isolating the manifestations of topological friction. The results ought to be relevant for linked biomolecules in experimental or biological confining conditions.

Optimal navigation strategy of active Brownian particles in target-search problems.

We investigate exploration patterns of a microswimmer, modeled as an active Brownian particle, searching for a target region located in a well of an energy landscape and separated from the initial position of the particle by high barriers. We find that the microswimmer can enhance its success rate in finding the target by tuning its activity and its persistence in response to features of the environment. The target-search patterns of active Brownian particles are counterintuitive and display characteristics robust to changes in the energy landscape. On the contrary, the transition rates and transition-path times are sensitive to the details of the specific energy landscape. In striking contrast to the passive case, the presence of additional local minima does not significantly slow down the active-target-search dynamics.

Emergence of effective temperatures in an out-of-equilibrium model of biopolymer folding.

We investigate the possibility of extending the notion of temperature in a stochastic model for the RNA or protein folding driven out of equilibrium. We simulate the dynamics of a small RNA hairpin subject to an external pulling force, which is time-dependent. First, we consider a fluctuation-dissipation relation (FDR) whereby we verify that various effective temperatures can be obtained for different observables, only when the slowest intrinsic relaxation timescale of the system regulates the dynamics of the system. Then, we introduce a different nonequilibrium temperature, which is defined from the rate of heat exchanged with a weakly interacting thermal bath. Notably, this "kinetic" temperature can be defined for any frequency of the external switching force. We also discuss and compare the behavior of these two emerging parameters, by discriminating the time-delayed nature of the FDR temperature from the instantaneous character of the kinetic temperature. The validity of our numerics are corroborated by a simple four-state Markov model which describes the long-time behavior of the RNA molecule.

Entropic measure unveils country competitiveness and product specialization in the World trade web.

We show how the Shannon entropy function can be used as a basis to set up complexity measures weighting the economic efficiency of countries and the specialization of products beyond bare diversification. This entropy function guarantees the existence of a fixed point which is rapidly reached by an iterative scheme converging to our self-consistent measures. Our approach naturally allows to decompose into inter-sectorial and intra-sectorial contributions the country competitivity measure if products are partitioned into larger categories. Besides outlining the technical features and advantages of the method, we describe a wide range of results arising from the analysis of the obtained rankings and we benchmark these observations against those established with other economical parameters. These comparisons allow to partition countries and products into various main typologies, with well-revealed characterizing features. Our methods have wide applicability to general problems of ranking in bipartite networks.

Target Search of Active Agents Crossing High Energy Barriers.

Target search by active agents in rugged energy landscapes has remained a challenge because standard enhanced sampling methods do not apply to irreversible dynamics. We overcome this nonequilibrium rare-event problem by developing an algorithm generalizing transition-path sampling to active Brownian dynamics. This method is exemplified and benchmarked for a paradigmatic two-dimensional potential with a high barrier. We find that even in such a simple landscape the structure and kinetics of the ensemble of transition paths changes drastically in the presence of activity. Indeed, active Brownian particles reach the target more frequently than passive Brownian particles, following longer and counterintuitive search patterns.

Topological Disentanglement of Linear Polymers under Tension.

We develop a theoretical description of the topological disentanglement occurring when torus knots reach the ends of a semiflexible polymer under tension. These include decays into simpler knots and total unknotting. The minimal number of crossings and the minimal knot contour length are the topological invariants playing a key role in the model. The crossings behave as particles diffusing along the chain and the application of appropriate boundary conditions at the ends of the chain accounts for the knot disentanglement. Starting from the number of particles and their positions, suitable rules allow reconstructing the type and location of the knot moving on the chain Our theory is extensively benchmarked with corresponding molecular dynamics simulations and the results show a remarkable agreement between the simulations and the theoretical predictions of the model.

Dynamical properties of densely packed confined hard-sphere fluids.

Numerical solutions of the mode-coupling theory (MCT) equations for a hard-sphere fluid confined between two parallel hard walls are elaborated. The governing equations feature multiple parallel relaxation channels which significantly complicate their numerical integration. We investigate the intermediate scattering functions and the susceptibility spectra close to structural arrest and compare to an asymptotic analysis of the MCT equations. We corroborate that the data converge in the ß-scaling regime to two asymptotic power laws, viz. the critical decay and the von Schweidler law. The numerical results reveal a nonmonotonic dependence of the power-law exponents on the slab width and a nontrivial kink in the low-frequency susceptibility spectra. We also find qualitative agreement of these theoretical results to event-driven molecular dynamics simulations of polydisperse hard-sphere systems. In particular, the nontrivial dependence of the dynamical properties on the slab width is well reproduced.

Transition path times in asymmetric barriers.

Biomolecular conformational transitions are usually modeled as barrier crossings in a free energy landscape. The transition paths connect two local free energy minima and transition path times (TPT) are the actual durations of the crossing events. The simplest model employed to analyze TPT and to fit empirical data is that of a stochastic particle crossing a parabolic barrier. Motivated by some disagreement between the value of the barrier height obtained from the TPT distributions as compared to the value obtained from kinetic and thermodynamic analyses, we investigate here TPT for barriers which deviate from the symmetric parabolic shape. We introduce a continuous set of potentials, that starting from a parabolic shape, can be made increasingly asymmetric by tuning a single parameter. The TPT distributions obtained in the asymmetric case are very well-fitted by distributions generated by parabolic barriers. The fits, however, provide values for the barrier heights and diffusion coefficients which deviate from the original input values. We show how these findings can be understood from the analysis of the eigenvalues spectrum of the Fokker-Planck equation and highlight connections with experimental results.

Plectoneme dynamics and statistics in braided polymers.

Braids composed of two interwoven polymer chains exhibit a "buckling" transition whose origin has been explained through the onset of plectonemic structures. Here we study, by a combination of simulation and analytics, the dynamics of plectoneme formation and their statistics in steady state. The introduction of an order parameter-the plectonemic fraction-allows us to map out the phase boundary between the straight-braid phase and the plectonemic one. We then monitor the formation and the growth of plectonemes, observing events typical of phase separation kinetics for liquid-gas systems (fusion, fission, and one-dimensional Ostwald ripening) but also of DNA supercoiling dynamics (plectonemic hopping). Finally, we propose a stochastic field theory for the coupled dynamics of twist and local writhe which explains the phenomenology found with Brownian dynamics simulations as well as the power laws underlying the coarsening of plectonemes.

Overtwisting induces polygonal shapes in bent DNA.

By combining analytical results and simulations of various coarse-grained models, we investigate the minimal energy shape of DNA minicircles which are torsionally constrained by an imposed over or undertwist. We show that twist-bend coupling, a cross interaction term discussed in the recent DNA literature, induces minimal energy shapes with a periodic alternation of parts with high and low curvature resembling rounded polygons. We briefly discuss the possible experimental relevance of these findings. We finally show that the twist and bending energies of minicircles are governed by renormalized stiffness constants, rather than the bare ones. This has important consequences for the analysis of experiments involving circular DNA meant to determine DNA elastic constants.

Topological Disentanglement Dynamics of Torus Knots on Open Linear Polymers.

We simulate and study the topological disentanglement occurring when torus knots reach the ends of a semiflexible open polymer (decay into simpler knots or unknotting). Through a rescaling procedure and the application of appropriate boundary conditions, we show that the full unknotting process can be understood in terms of point-like particles representing essential crossings, diffusing on the support [0, 1]. We address the bending and configurational free energy drives on the diffusion process, together with the scaling properties of the effective diffusion and friction coefficients. Agreement with simulations suggests universal features for these two model parameters.

Topologically Linked Chains in Confinement.

We examine how channel confinement affects the equilibrium properties of topologically linked ring polymers and, by contrast, of equivalent unlinked rings, too. By performing extensive simulations of semiflexible rings of different chain length, N, and channel diameter, D, we discover three notable properties purely due to linking. First, upon entering the weak confinement regime, the length of the physically linked portion, lLKThe, becomes independent of chain length. Next, even when confinement is strong enough to pull apart and segregate unlinked rings, lLK stays much larger than in the highly stretched limit. Finally, at fixed N, lLK varies approximately as D0.5, and we provide a simple scaling argument for this power-law behavior. These properties, which may hold for different link topologies, can be tested by current experimental setups on DNA rings confined in microchannels. Moreover, they could be relevant for the efficient in vivo unlinking of newly replicated bacterial chromosomes.

Growth dynamics and complexity of national economies in the global trade network.

We explore the quantitative nexus among economic growth of a country, diversity and specialization of its productions, and evolution in time of its basket of exports. To this purpose we set up a dynamic model and construct economic complexity measures based on panel data concerning up to 1238 exports of 223 countries for 21 years. Key statistical features pertaining to the distribution of resources in the different exports of each country reveal essential in both cases. The parameters entering the evolution model, combined with counterfactual analyses of synthetic simulations, give novel insight into cooperative effects among different productions and prospects of growth of each economy. The complexity features emerging from the analysis of dynamics are usefully compared with gross domestic product per capita (GDPpc) and with an original measure of the efficiency of the economic systems. This measure, whose construction starts from an estimate of bare diversity in terms of Shannon's entropy function, is made fully consistent with the degree of specialization of the products. Comparisons of this measure with the model parameters allow clear distinctions, from multiple perspectives, among developed, emerging, underdeveloped and risky economies.