Search

Directed motion up a concentration gradient is crucial for the survival and maintenance of numerous biological systems, such as sperms moving towards an egg during fertilization or ciliates moving towards a food source. In these systems, chirality - manifested as a rotational torque - plays a vital role in facilitating directed motion. While systematic studies of active molecules in activity gradients exist, the effect of chirality remains little studied. In this study, we examine the simplest case of a chiral active particle connected to a passive particle in a spatially varying activity field. We demonstrate that this minimal setup can exhibit rich emergent tactic behaviors, with the chiral torque serving as the tuning parameter. Notably, when the chiral torque is sufficiently large, even a small passive particle enables the system to display the desired accumulation behavior. This study provides valuable insights into the design principles of hybrid bio-molecular devices of the future., Comment: 14 pages (6 main + 8 supplementary), 4 figures

Polymer-like structures are ubiquitous in nature and synthetic materials. Their configurational and migration properties are often affected by crowded environments leading to non-thermal fluctuations. Here, we study an ideal Rouse chain in contact with a non-homogeneous active bath, characterized by the presence of active self-propelled agents which exert time-correlated forces on the chain. By means of a coarse-graining procedure, we derive an effective evolution for the center of mass of the chain and show its tendency to migrate towards and preferentially localize in regions of high/low bath activity depending on the model parameters. In particular, we demonstrate that an active bath with non-uniform activity can be used to separate efficiently polymeric species with different lengths and/or connectivity., Comment: 10 pages, 4 figures

We use machine learning algorithms to detect the crystalline phase in undercooled melts in molecular dynamics simulations. Our classification method is based on local conformation and environmental fingerprints of individual monomers. In particular, we employ self-supervised auto-encoders to compress the fingerprint information and a Gaussian mixture model to distinguish ordered states from disordered ones. The resulting identification of crystalline monomers agrees to a large extent with human-defined classifiers such as the stem-length-based classification scheme as developed in our previous work [C. Luo and J.-U. Sommer, Macromolecules 44 (2011), 1523], but does not require any foreknowledge about the structure of semi-crystalline polymers. Because of its local sensitivity, the method allows the resolution of detailed time patterns of crystalline order before an apparent signature of the transition is visible in thermodynamic properties such as for the specific volume. At a pre-transition point, we observe the highest crystallization efficiency using the fraction of monomers being conserved in the crystalline phase as compared to the number of monomers joining that phase., Comment: 41 pages, 19 figures, 5 tables

The force autocorrelation function (FACF), a concept of fundamental interest in statistical mechanics, encodes the effect of interactions on the dynamics of a tagged particle. In equilibrium, the FACF is believed to decay monotonically in time which is a signature of slowing down of the dynamics of the tagged particle due to interactions. Here we analytically show that in odd-diffusive systems, which are characterized by a diffusion tensor with antisymmetric elements, the FACF can become negative and even exhibit temporal oscillations. We also demonstrate that, despite the isotropy, the knowledge of FACF alone is not sufficient to describe the dynamics: the full autocorrelation tensor is required and contains an antisymmetric part. These unusual properties translate into enhanced dynamics of the tagged particle quantified via the self-diffusion coefficient that, remarkably, increases due to particle interactions.

The probability per unit time for a thermally activated Brownian particle to escape over a potential well is in general well-described by Kramers theory. Kramers showed that the escape time decreases exponentially with increasing barrier height. The dynamics slow down when the particle is charged and subjected to a Lorentz force due to an external magnetic field. This is evident via a rescaling of the diffusion coefficient entering as a prefactor in the Kramers escape rate without any impact on the barrier-height-dependent exponent. Here we show that the barrier height can be effectively changed when the charged particle is subjected to an external vortex flow. While the external vortex alone does not affect the mean escape time of the particle, when combined with a magnetic field it effectively pushes the fluctuating particle either radially outside or inside depending on its sign relative to that of the magnetic field. In particular, the effective potential over which the particle escapes can be changed to a flat, a stable, and an unstable potential by tuning the signs and magnitudes of the external vortex and the applied magnetic field. Notably, the last case corresponds to enhanced escape dynamics.

While the behavior of active colloidal molecules is well studied by now for a constant activity, the effect of activity gradients is much less understood. Here we explore one of the simplest molecules in activity gradients, namely active chiral dimers composed of two particles with opposite active torques of the same magnitude. We show analytically that with increasing torque, the dimer switches its behavior from antichemotactic to chemotactic. The origin of the emergent chemotaxis is the cooperative exploration of activity gradient by the two particles. While one of the particles moves into higher activity regions, the other moves towards lower activity region resulting in a net bias in the direction of higher activity. We do a comparative study of chiral active particles with charged Brownian particles under magnetic field and show that despite the fundamental similarity in terms of their odd-diffusive behavior, their dynamics and chemotactic behavior are generally not equivalent. We demonstrate this explicitly in a dimer composed of oppositely charged active particles, which remains antichemotactic for any magnetic field.

While the dynamics of dimers and polymer chains in a viscous solvent is well understood within the celebrated Rouse model, the effect of an external magnetic field on the dynamics of a charged chain is much less understood. Here we generalize the Rouse model for a charged dimer to include the effect of an external magnetic field. Our analytically solvable model allows a fundamental insight into the magneto-generated dynamics of the dimer in the overdamped limit as induced by the Lorentz-force. Surprisingly, for a dimer of oppositely charged particles, we find an enormous enhancement of the dynamics of the dimer center which exhibits even a transient superballistic behavior. This is highly unusual in an overdamped system for there is neither inertia nor any internal or external driving. We attribute this to a significant translation and rotation coupling due to the Lorentz force. We also find that magnetic field reduces the mobility of a dimer along its orientation and its effective rotational diffusion coefficient. In principle, our predictions can be tested by experiments with colloidal particles and complex plasmas.

It is generally believed that collisions of particles reduce the self-diffusion coefficient. Here we show that in odd-diffusive systems, which are characterized by diffusion tensors with antisymmetric elements, collisions surprisingly can enhance the self-diffusion. In these systems, due to an inherent curving effect, the motion of particles is facilitated, instead of hindered by collisions leading to a mutual rolling effect. Using a geometric model, we analytically predict the enhancement of the self-diffusion coefficient with increasing density. This counterintuitive behaviour is demonstrated in the archetypal odd-diffusive system of Brownian particles under Lorentz force. We validate our findings by many body Brownian dynamics simulations in dilute systems.

We consider a rigid assembly of two active Brownian particles, forming an active colloidal dimer, in a gradient of activity. We show analytically that depending on the relative orientation of the two particles the active dimer accumulates in regions of either high or low activity, corresponding to, respectively, chemotaxis and antichemotaxis. Certain active dimers show both chemotactic and antichemotactic behavior, depending on the strength of the activity. Our coarse-grained Fokker-Planck approach yields an effective potential, which we use to construct a nonequilibrium phase diagram that classifies the dimers according to their tactic behavior. Moreover, we show that for certain dimers a higher persistence of the motion is achieved similar to the effect of a steering wheel in macroscopic devices. This work could be useful for designing autonomous active colloidal structures which adjust their motion depending on the local activity gradients.

Thermally activated escape of a Brownian particle over a potential barrier is well understood within Kramers theory. When subjected to an external magnetic field, the Lorentz force slows down the escape dynamics via a rescaling of the diffusion coefficient without affecting the exponential dependence on the barrier height. Here, we study the escape dynamics of a charged Brownian particle from a two-dimensional truncated harmonic potential under the influence of Lorentz force due to an external magnetic field. The particle is driven anisotropically by subjecting it to noises with different strengths along different spatial directions. We show that the escape time can largely be tuned by the anisotropic driving. While the escape process becomes anisotropic due to the two different noises, the spatial symmetry is restored in the limit of large magnetic fields. This is attributed to the Lorentz force induced coupling between the spatial degrees of freedom which makes the difference between two noises irrelevant at high magnetic fields. The theoretical predictions are verified by Brownian dynamics simulations. In principle, our predictions can be tested by experiments with a Brownian gyrator in the presence of a magnetic field.

A Brownian particle performs gyrating motion around a potential energy minimum when subjected to thermal noises from two different heat baths. Here, we propose a magneto-gyrator made of a single charged Brownian particle that is steered by an external magnetic field. Key properties, such as the direction of gyration, the torque exerted by the engine on the confining potential and the maximum power delivered by the microengine can be tuned by varying the strength and direction of the applied magnetic field. Further tunability is obtained by rotating the potential in the plane perpendicular to the direction of the magnetic filed. We show that in this generic scenario, the microengine can be stalled and even reversed by the magnetic field. Finally, we highlight a property of the magneto-gyrator that has no counterpart in the overdamped approximation--the heat loss from the hot to cold bath requires explicit knowledge of the mass of the particle. Consequently, the efficiency of the microengine is mass-dependent even in the overdamped limit., Comment: 14 pages, 4 figures

Different methods for creating Olympic gels are analyzed using computer simulations. First ideal reference samples are obtained from freely interpenetrating semi-dilute solutions and melts of cyclic polymers. The distribution of pairwise concatenations per cyclic molecule is given by a Poisson-distribution and can be used to describe the elastic structure of the gels. Several batches of linear chains decorated with different selectively binding groups at their ends are mixed in the "DNA Origami" technique and network formation is realized. While the formation of cyclic molecules follows mean field predictions below overlap of the precursor molecules, an enhanced ring formation above overlap is found that is not explained by mean field arguments. The "progressive construction" method allows to create Olympic gels with a single reaction step from a concentrated mixture of large compressed rings with a low weight fraction short chains that are below overlap concentration. This method, however, is limited by the difficulty to obtain a sufficiently high degree of polymerization of the large rings.

The bond fluctuation method is used to simulate both non-concatenated entangled and interpenetrating melts of ring polymers. We find that the swelling of interpenetrating rings upon dilution follows the same laws as for linear chains. Knotting and linking probabilities of ring polymers in semi-dilute solution are analyzed using the HOMFLY polynomial. We find an exponential decay of the knotting probability of rings. The correlation length of the semi-dilute solution can be used to superimpose knotting data at different concentrations. A power law dependence $f_{n}\sim\phi R^{2}\sim\phi^{0.77}N$ for the average number $f_{n}$ of linked rings per ring at concentrations larger than the overlap volume fraction of rings $\phi^{*}$ is determined from the simulation data. The fraction of non-concatenated rings displays an exponential decay $P_{OO}\sim\exp(-f_{n})$, which indicates $f_{n}$ to provide the entropic effort for not forming concatenated conformations. Based upon this results we find four different regimes for the conformations of rings in melts that are separated by a critical lengths $N_{OO}$, $N_{C}$ and $N^{*}$. $N_{OO}$ describes the onset of the effect of non-concatenation below which topological effects are not important, $N_{C}$ is the cross-over between weak and strong compression of rings, and $N^{*}$ is defined by the cross-over from a non-concatenation contribution $f_{n}\sim\phi R^{2}$ to an overlap dominated concatenation contribution $f_{n}\sim\phi N^{1/2}$ at $N>N^{*}$. For $N_{OO}

The swelling equilibrium of Olympic gels, which are composed of entangled cyclic polymers, is studied by Monte Carlo Simulations. In contrast to chemically crosslinked polymer networks, we observe that Olympic gels made of chains with a \emph{larger} degree of polymerization, $N$, exhibit a \emph{smaller} equilibrium swelling degree, $Q\propto N^{-0.28}\phi_{0}^{-0.72}$, at the same polymer volume fraction $\phi_{0}$ at network preparation. This observation is explained by a desinterspersion process of overlapping non-concatenated rings upon swelling., Comment: arXiv admin note: substantial text overlap with arXiv:2103.15480

A mean field rate theory description of the homo- and co-polymerization of $f$-functional molecules is developed, which contains the formation of short cyclic structures inside the network. The predictions of this model are compared with Monte-Carlo simulations of cross-linking of star polymers in solution. We find that homo-polymerizations are well captured by mean-field models at concentrations larger than one quarter of the geometrical overlap concentration. All simulation data can be fit using a single geometric parameter for cyclization. The simulation data reveal that within the range of parameters of the present study correlations among multiply connected molecules can be neglected. Thus, mean-field treatments of homopolymerizations are reasonable approximations, if short cycles are properly addressed. Co-polymerization is considered in the case of strict A-B reactions, where all reactive groups of individual molecules are either of type A or B. For these systems we find a clear influence of the local intermixing of A and B groups for all concentrations investigated. In consequence, mean-field models are less appropriate to describe the simulation data. The lack of ring structures containing an odd number of molecules as compared to homopolymerizations at the same extent of reaction allows for the formation of stable AB networks at concentrations one order of magnitude below the geometrical overlap concentration.

Using the bond fluctuation model we study polymer networks obtained by endlinking of symmetric 4-arm star polymers. We consider two types of systems. Solutions of one type (A) of star polymers and solution of two types (A,B) of star polymer where A-type polymers can only crosslink with B-type polymers. We find that network defects in $A$ networks are dominated by short dangling loops close to overlap concentration $c^{*}$. $AB$ networks develop a more perfect network structure, since loop sizes involving an odd number of stars are impossible, and thus, the most frequent dangling loop with largest impact on the phantom modulus is absent. The analysis of the pair-correlation and scattering function reveals that there is an amorphous packing of $A$ and $B$ type stars with a homogenization of $A$ and $B$ concentrations upon cross-linking at intermediate length scales in contrast to the previously suggested crystalline like order of $A$ and $B$ components at $c^{*}$. This result is corroborated by the coincidence of the probabilities of the shortest loop structures (which is impossible upon the previously suggested packing of stars) in both types of networks. Furthermore, we derive the vector order parameters associated with the most frequent network structures based on the phantom model. In particular, for $AB$ networks we can show that there is a dominating cyclic defect with a clearly separated order parameter that could be used to analyze cyclic network defects.

Concatenation and equilibrium swelling of Olympic gels, which are composed of entangled cyclic polymers, is studied by Monte Carlo Simulations. The average number of concatenated molecules per cyclic polymers, $f_n$, is found to depend on the degree of polymerization, $N$, and polymer volume fraction at network preparation, ${\phi}_0$, as $f_n ~ {\phi}_0^{{\nu}/(3{\nu}-1)}N$ with scaling exponent ${\nu} = 0.588$. In contrast to chemically cross-linked polymer networks, we observe that Olympic gels made of longer cyclic chains exhibit a smaller equilibrium swelling degree, $Q ~ N^{-0.28} {\phi}_0^{-0.72}$, at the same polymer volume fraction ${\phi}_0$. This observation is explained by a desinterspersion process of overlapping non-concatenated rings upon swelling, which is tested directly by analyzing the change in overlap of the molecules upon swelling.

Active particles with their characteristic feature of self-propulsion are regarded as the simplest models for motility in living systems. The accumulation of active particles in low activity regions has led to the general belief that chemotaxis requires additional features and at least a minimal ability to process information and to control motion. We show that self-propelled particles display chemotaxis and move into regions of higher activity, if the particles perform work on passive objects, or cargo, to which they are bound. The origin of this cooperative chemotaxis is the exploration of the activity gradient by the active particle when bound to a load, resulting in an average excess force on the load in the direction of higher activity. Using a minimalistic theoretical model, we capture the most relevant features of these active-passive dimers and in particular we predict the crossover between anti-chemotactic and chemotactic behaviour. Moreover we show that merely connecting active particles to chains is sufficient to obtain the crossover from anti-chemotaxis to chemotaxis with increasing chain length. Such an active complex is capable of moving up a gradient of activity such as provided by a gradient of fuel and to accumulate where the fuel concentration is at its maximum. The observed transition is of significance to proto-forms of life enabling them to locate a source of nutrients even in the absence of any supporting sensomotoric apparatus.

We study the motion of a Brownian particle subjected to Lorentz force due to an external magnetic field. Each spatial degree of freedom of the particle is coupled to a different thermostat. We show that the magnetic field results in correlation between different velocity components in the stationary state. Integrating the velocity autocorrelation matrix, we obtain the diffusion matrix that enters the Fokker-Planck equation for the probability density. The eigenvectors of the diffusion matrix do not align with the temperature axes. As a consequence the Brownian particle performs spatially correlated diffusion. We further show that in the presence of an isotropic confining potential, an unusual, flux-free steady state emerges which is characterized by a non-Boltzmann density distribution, which can be rotated by reversing the magnetic field. The nontrivial steady state properties of our system result from the Lorentz force induced coupling of the spatial degrees of freedom which cease to exist in equilibrium corresponding to a single-temperature system.

In systems with overdamped dynamics, the Lorentz force reduces the diffusivity of a Brownian particle in the plane perpendicular to the magnetic field. The anisotropy in diffusion implies that the Fokker-Planck equation for the probabiliy distribution of the particle acquires a tensorial coefficient. The tensor, however, is not a typical diffusion tensor due to the antisymmetric elements which account for the fact that Lorentz force curves the trajectory of a moving charged particle. This gives rise to unusual dynamics with features such as additional Lorentz fluxes and a nontrivial density distribution, unlike a diffusive system. The equilibrium properties are, however, unaffected by the Lorentz force. Here we show that by stochastically resetting the Brownian particle, a nonequilibrium steady state can be created which preserves the hallmark features of dynamics under Lorentz force. We then consider a minimalistic example of spatially inhomogeneous magnetic field, which shows how Lorentz fluxes fundamentally alter the boundary conditions giving rise to an unusual stationary state.

Using Brownian dynamics simulations, the motion of active Brownian particles (ABPs) in the presence of fuel (or 'food') sources is studied. It is an established fact that within confined stationary systems, the activity of ABPs generates density profiles that are enhanced in regions of low activity, which is generally referred to as 'anti-chemotaxis'. We demonstrate that -- contrary to common believes -- in non-stationary setups, emerging here as a result of short fuel bursts, our model ABPs do instead exhibit signatures of chemotactic behavior. In direct competition with inactive, but otherwise identical Brownian particles (BPs), the ABPs are shown to fetch a larger amount of food. From a biological perspective, the ability to turn active would, despite of the absence of sensoric devices, encompass an evolutionary advantage., Comment: 6 pages, 7 figures

The Fokker-Planck equation provides complete statistical description of a particle undergoing random motion in a solvent. In the presence of Lorentz force due to an external magnetic field, the Fokker-Planck equation picks up a tensorial coefficient, which reflects the anisotropy of the particle's motion. This tensor, however, can not be interpreted as a diffusion tensor; there are antisymmetric terms which give rise to fluxes perpendicular to the density gradients. Here, we show that for an inhomogeneous magnetic field these nondiffusive fluxes have finite divergence and therefore affect the density evolution of the system. Only in the special cases of a uniform magnetic field or carefully chosen initial condition with the same symmetry as the magnetic field can these fluxes be ignored in the density evolution., Comment: 7 Pages, 4 Figures

We consider the nonequilibrium dynamics of a charged active Brownian particle in the presence of a space dependent magnetic field. It has recently been shown that the Lorentz force induces a particle flux perpendicular to density gradients, thus preventing a diffusive description of the dynamics. Whereas a passive system will eventually relax to an equilibrium state, unaffected by the magnetic field, an active system subject to a spatially varying Lorentz force settles into a nonequilibrium steady state characterized by an inhomogeneous density and divergence-free bulk fluxes. A macroscopic flux of charged active particles is induced by the gradient of the magnetic field only and does not require additional symmetric breaking such as density or potential gradients. This stands in marked contrast to similar phenomena in condensed matter such as the classical Hall effect. In a confined geometry we observe circulating fluxes, which can be reversed by inverting the direction of the magnetic field. Our theoretical approach, based on coarse-graining of the Fokker-Planck equation, yields analytical results for the density, fluxes, and polarization in the steady state, all of which are validated by direct computer simulation. We demonstrate that passive tracer particles can be used to measure the essential effects of the Lorentz force on the active particle bath, and we discuss under which conditions the effects of the flux could be observed experimentally., Comment: 18 Pages, 6 Figures

In this paper we study a system of entangled chains that bear reversible cross-links in a melt state. The cross-links are tethered uniformly on the backbone of each chain. A slip-link type model for the system is presented and solved for the relaxation modulus. The effects of entanglements and reversible cross-linkers are modelled as discrete form of constraints that influence the motion of the primitive path. In contrast to a non-associating entangled system the model calculations demonstrate that the elastic modulus has a much higher first plateau and a delayed terminal relaxation. These effects are attributed to the evolution of the entangled chains as influenced by tethered reversible linkers. The model is solved for the case when linker survival time $\tau_s$ is greater than the entanglement time $\tau_e$ but less than the Rouse time $\tau_R$.

We study dynamical properties of confined, self-propelled Brownian particles in an inhomogeneous activity profile. Using Brownian dynamics simulations, we calculate the probability to reach a fixed target and the mean first passage time to the target of an active particle. We show that both these quantities are strongly influenced by the inhomogeneous activity. When the activity is distributed such that high-activity zone is located between the target and the starting location, the target finding probability is increased and the passage time is decreased in comparison to a uniformly active system. Moreover, for a continuously distributed profile, the activity gradient results in a drift of active particle up the gradient bearing resemblance to chemotaxis. Integrating out the orientational degrees of freedom, we derive an approximate Fokker-Planck equation and show that the theoretical predictions are in very good agreement with the Brownian dynamics simulations., Comment: 7 pages, 5 figures

Monte Carlo simulations using an explicit solvent model indicate a new pathway for translocation of a polymer chain through a lipid bilayer. We consider a polymer chain composed of repeat units with a given hydrophobicity and a coarse-grained model of a lipid bilayer in the self-organized liquid state. By varying the degree of hydrophobicity the chain undergoes an adsorption transition with respect to the lipid bilayer. Close to the transition point, at a properly balanced hydrophobicity of the chain, the membrane becomes transparent with respect to the chain. At the same time the solvent permeability of the bilayer is strongly increased in the region close to adsorbed chain. Our results indicate that the critical point of adsorption of the polymer chain interacting with the fluctuating lipid bilayer could play a key role for the translocation of molecules though biological membranes., Comment: published in EPL

Via computer simulations, we demonstrate how a densely grafted layer of polymers, a {\it brush}, could be turned into an efficient switch through chemical modification of some of its end-monomers. In this way, a surface coating with reversibly switchable properties can be constructed. We analyze the fundamental physical principle behind its function, a recently discovered surface instability, and demonstrate that the combination of a high grafting density, an inflated end-group size and a high degree of monodispersity are conditions for an optimal functionality of the switch.

We consider adsorption of random copolymer chains onto an interface within the model of Garel et al. Europhysics Letters 8, 9 (1989). By using the replica method the adsorption of the copolymer at the interface is mapped onto the problem of finding the ground state of a quantum mechanical Hamiltonian. To study this ground state we introduce a novel variational principle for the Green's function, which generalizes the well-known Rayleigh-Ritz method of Quantum Mechanics to nonstationary states. Minimization with an appropriate trial Green's function enables us to find the phase diagram for the localization-delocalization transition for an ideal random copolymer at the interface., Comment: 5 pages

We investigate numerically, using the bond-fluctuation model, the adsorption of many random AB--copolymers with excluded volume interactions at the interface between two solvents. We find two regimes, controlled by the total number of polymers. In the first (dilute) regime, the copolymers near the interface extend parallel to it, while in the second regime they extend perpendicular to it. The density at the interface and the density in the bulk depend differently on the total number of copolymers: In the first regime the density at the interface increases more rapidly then in the bulk, whereas the opposite is true in the second regime., Comment: Latex, 14 pages, postscript figures

We investigate numerically using the bond--fluctuation model the adsorption of a random AB--copolymer at the interface between two solvents. From our results we infer several scaling relations: the radius of gyration of the copolymer in the direction perpendicular to the interface ($R_{gz}$) scales with $\chi$, the interfacial selectivity strength, as $R_{gz}=N^{\nu}f(\sqrt{N}\chi)$ where $\nu$ is the usual Flory exponent and $N$ is the copolymer's length; furthermore the monomer density at the interface scales as $\chi^{2\nu}$ for small $\chi$. We also determine numerically the monomer densities in the two solvents and discuss their dependence on the distance from the interface., Comment: Latex text file appended with figures.tar.gz