1
|
Cherayil BJ. Survival probabilities and first-passage distributions of self-propelled particles in spherical cavities. Phys Rev E 2023; 108:054607. [PMID: 38115486 DOI: 10.1103/physreve.108.054607] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2023] [Accepted: 11/05/2023] [Indexed: 12/21/2023]
Abstract
A model of self-propelled motion in a closed compartment containing simple or complex fluids is formulated in this paper in terms of the dynamics of a point particle moving in a spherical cavity under the action of random thermal forces and exponentially correlated noise. The particle's time evolution is governed by a generalized Langevin equation (GLE) in which the memory function, connected to the thermal forces by a fluctuation-dissipation relation, is described by Jeffrey's model of viscoelasticity (which reduces to a model of ordinary viscous dynamics in a suitable limit). The GLE is transformed exactly to a Fokker-Planck equation that in spherical polar coordinates is in turn found to admit of an exact solution for the particle's probability density function under absorbing boundary conditions at the surface of the sphere. The solution is used to derive an expression (that is also exact) for the survival probability of the particle in the sphere, starting from its center, which is then used to calculate the distribution of the particle's first-passage times to the boundary. The behavior of these quantities is investigated as a function of the Péclet number and the persistence time of the athermal forces, providing insight into the effects of nonequilibrium fluctuations on confined particle motion in three dimensions.
Collapse
Affiliation(s)
- Binny J Cherayil
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, Karnataka, India
| |
Collapse
|
2
|
Abelenda-Núñez I, Ortega F, Rubio RG, Guzmán E. Anomalous Colloidal Motion under Strong Confinement. SMALL (WEINHEIM AN DER BERGSTRASSE, GERMANY) 2023; 19:e2302115. [PMID: 37116105 DOI: 10.1002/smll.202302115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2023] [Revised: 04/13/2023] [Indexed: 06/19/2023]
Abstract
Diffusion of biological macromolecules in the cytoplasm is a paradigm of colloidal diffusion in an environment characterized by a strong restriction of the accessible volume. This makes of the understanding of the physical rules governing colloidal diffusion under conditions mimicking the reduction in accessible volume occurring in the cell cytoplasm, a problem of a paramount importance. This work aims to study how the thermal motion of spherical colloidal beads in the inner cavity of giant unilamellar vesicles (GUVs) is modified by strong confinement conditions, and the viscoelastic character of the medium. Using single particle tracking, it is found that both the confinement and the environmental viscoelasticity lead to the emergence of anomalous motion pathways for colloidal microbeads encapsulated in the aqueous inner cavity of GUVs. This anomalous diffusion is strongly dependent on the ratio between the volume of the colloidal particle and that of the GUV under consideration as well as on the viscosity of the particle's liquid environment. Therefore, the results evidence that the reduction of the free volume accessible to colloidal motion pushes the diffusion far from a standard Brownian pathway as a result of the change in the hydrodynamic boundary conditions driving the particle motion.
Collapse
Affiliation(s)
- Irene Abelenda-Núñez
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Ciudad Universitaria s/n., Madrid, 28040, Spain
| | - Francisco Ortega
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Ciudad Universitaria s/n., Madrid, 28040, Spain
- Unidad de Materia Condensada, Instituto Pluridisciplinar, Universidad Complutense de Madrid, Paseo Juan XXIII 1., Madrid, 28040, Spain
| | - Ramón G Rubio
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Ciudad Universitaria s/n., Madrid, 28040, Spain
| | - Eduardo Guzmán
- Departamento de Química Física, Facultad de Ciencias Químicas, Universidad Complutense de Madrid, Ciudad Universitaria s/n., Madrid, 28040, Spain
- Unidad de Materia Condensada, Instituto Pluridisciplinar, Universidad Complutense de Madrid, Paseo Juan XXIII 1., Madrid, 28040, Spain
| |
Collapse
|
3
|
Sahoo R, Chakrabarti R. Structure and dynamics of an active polymer chain inside a nanochannel grafted with polymers. SOFT MATTER 2023; 19:5978-5988. [PMID: 37497754 DOI: 10.1039/d3sm00618b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
We use computer simulations to investigate the complex dynamics of a polymer, made of active Brownian particles, inside a channel grafted internally with passive polymer chains. Our simulations reveal that this probe-polymer, if passive, exhibits a compact structure when its interaction is repulsive with the grafted chains as it tends to stay within the hollow space created along the axis of the channel. On increasing the attractive interaction, the passive probe-polymer is pulled towards the grafted polymeric region and adopts an extended structure. By contrast, switching on the activity helps the probe-polymer to escape from the local traps caused by the sticky grafted chains. The interplay between the activity of the probe-polymer and its sticky interaction with the grafted chains results in shrinking, followed by swelling as the activity is increased. To elucidate the dynamics we compute the mean square displacement (MSD) of the center of mass of the probe-polymer, which increases monotonically with activity and displays superdiffusive behavior at an intermediate time and enhanced diffusion at a long time period. In addition, compared with the attractive interaction, the active probe-polymer shows faster dynamics when the interaction is repulsive to the grafted polymers. We believe that our current study will provide insights into the structural changes and dynamics of active polymers in heterogeneous media and will be useful in designing polymer-based drug delivery vehicles.
Collapse
Affiliation(s)
- Rajiblochan Sahoo
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| | - Rajarshi Chakrabarti
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| |
Collapse
|
4
|
Yan R, Tan F, Wang J, Zhao N. Conformation and dynamics of an active filament in crowded media. J Chem Phys 2023; 158:114905. [PMID: 36948796 DOI: 10.1063/5.0142559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/04/2023] Open
Abstract
The structural and dynamical properties of active filamentous objects under macromolecular crowding have a great relevance in biology. By means of Brownian dynamics simulations, we perform a comparative study for the conformational change and diffusion dynamics of an active chain in pure solvents and in crowded media. Our result shows a robust compaction-to-swelling conformational change with the augment of the Péclet number. The presence of crowding facilitates self-trapping of monomers and, thus, reinforces the activity mediated compaction. In addition, the efficient collisions between the self-propelled monomers and crowders induce a coil-to-globulelike transition, indicated by a marked change of the Flory scaling exponent of the gyration radius. Moreover, the diffusion dynamics of the active chain in crowded solutions demonstrates activity-enhanced subdiffusion. The center of mass diffusion manifests rather new scaling relations with respect to both the chain length and Péclet number. The interplay of chain activity and medium crowding provides a new mechanism to understand the non-trivial properties of active filaments in complex environments.
Collapse
Affiliation(s)
- Ran Yan
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Fei Tan
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Jingli Wang
- College of Chemistry, Sichuan University, Chengdu 610064, China
| | - Nanrong Zhao
- College of Chemistry, Sichuan University, Chengdu 610064, China
| |
Collapse
|
5
|
Yadav RS, Das C, Chakrabarti R. Dynamics of a spherical self-propelled tracer in a polymeric medium: interplay of self-propulsion, stickiness, and crowding. SOFT MATTER 2023; 19:689-700. [PMID: 36598025 DOI: 10.1039/d2sm01626e] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/17/2023]
Abstract
We employ computer simulations to study the dynamics of a self-propelled spherical tracer particle in a viscoelastic medium, made of a long polymer chain. Here, the interplay between viscoelasticity, stickiness, and activity (self-propulsion) brings additional complexity to the tracer dynamics. Our simulations show that on increasing the stickiness of the tracer particle to the polymer beads, the dynamics of the tracer particle slows down as it gets stuck to the polymer chain and moves along with it. But with increasing self-propulsion velocity, the dynamics gets enhanced. In the case of increasing stickiness as well as activity, the non-Gaussian parameter (NGP) exhibits non-monotonic behavior, which also shows up in the re-scaled self part of the van-Hove function. Non-Gaussianity results owing to the enhanced binding events and the sticky motion of the tracer along with the chain with increasing stickiness. On the other hand, with increasing activity, initially non-Gaussianity increases as the tracer moves through the heterogeneous polymeric environment but for higher activity, the tracer escapes resulting in a negative NGP. For higher values of stickiness, the trapping time distributions of the passive tracer particle broaden and have long tails. On the other hand, for a given stickiness with increasing self-propulsion force, the trapping time distributions become narrower and have short tails. We believe that our current simulation study will be helpful in elucidating the complex motion of activity-driven probes in viscoelastic media.
Collapse
Affiliation(s)
- Ramanand Singh Yadav
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| | - Chintu Das
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| | - Rajarshi Chakrabarti
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| |
Collapse
|
6
|
Kumar P, Chakrabarti R. Dynamics of self-propelled tracer particles inside a polymer network. Phys Chem Chem Phys 2023; 25:1937-1946. [PMID: 36541408 DOI: 10.1039/d2cp04253c] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
The transport of tracer particles through mesh-like environments such as biological hydrogels and polymer matrices is ubiquitous in nature. These tracers can be passive, such as colloids, or active (self-propelled), for example, synthetic nanomotors or bacteria. Computer simulations in principle could be extremely useful in exploring the mechanism of the active transport of tracer particles through mesh-like environments. Therefore, we construct a polymer network on a diamond lattice and use computer simulations to investigate the dynamics of spherical self-propelled particles inside the network. Our main objective is to elucidate the effect of the self-propulsion on the tracer particle dynamics as a function of the tracer size and the stiffness of the polymer network. We compute the time-averaged mean-squared displacement (MSD) and the van-Hove correlations of the tracer. On the one hand, in the case of a bigger sticky particle, the caging caused by the network particles wins over the escape assisted by the self-propulsion. This results an intermediate-time subdiffusion. On the other hand, smaller tracers or tracers with high self-propulsion velocities can easily escape from the cages and show intermediate-time superdiffusion. The stiffer the network, the slower the dynamics of the tracer, and bigger tracers exhibit longer lived intermediate time superdiffusion, since the persistence time scales as ∼σ3, where σ is the diameter of the tracer. At the intermediate time, non-Gaussianity is more pronounced for active tracers. At the long time, the dynamics of the tracer, if passive or weakly active, becomes Gaussian and diffusive, but remains flat for tracers with high self-propulsion, accounting for their seemingly unrestricted motion inside the network.
Collapse
Affiliation(s)
- Praveen Kumar
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| | - Rajarshi Chakrabarti
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai 400076, India.
| |
Collapse
|
7
|
Pial TH, Das S. Specific Ion and Electric Field Controlled Diverse Ion Distribution and Electroosmotic Transport in a Polyelectrolyte Brush Grafted Nanochannel. J Phys Chem B 2022; 126:10543-10553. [PMID: 36454705 DOI: 10.1021/acs.jpcb.2c05524] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022]
Abstract
Controlling ion distribution inside a charged nanochannel is central to using such channels in diverse applications. Here, we show the possibility of using a charged polyelectrolyte (PE) brush-grafted nanochannel for triggering diverse nanoscopic ion distribution and nanofluidic electroosmotic transport by controlling the valence and size of the counterions (that screen the charges of the PE brushes) and the strength of an externally applied axial electric field. We atomistically simulate separate cases of fully charged polyacrylic acid (PAA) brush functionalized nanochannels with Na+, Cs+, Ca2+, Ba2+, and Y3+ counterions screening the PE charges. Four key findings emerge from our simulations. First, we find that the counterions with a greater valence and a smaller size prefer to remain localized inside the brush layer. Second, for the case where there is an added chloride salt with the same cation (as the screening counterions), there are more coions (Cl- ions) in the brush-free bulk than counterions (for counterions Na+, Ca2+, Ba2+, Y3+): this is a manifestation of the overscreening (OS) of the PE brush layer. Contrastingly, the number of Cs+ ions remain higher than the Cl- ions inside the brush-free bulk, ensuring that there is no OS effect for this case. Third, large applied electric field enables a few Na+, Cs+, and Ba2+ counterions to leave the brush layer and to go to the bulk: this makes the OS of the PE brush layer disappear for the cases of PE brushes being screened by the Na+ and Ba2+ ions. On the other hand, no such electric-field-mediated disappearance of OS is observed for the cases of Ca2+ and Y3+ screening counterions; we attribute this to the firm attachment of these counterions to the negatively charged monomers. Free energy associated with a counterion binding to a PE chain corroborates this diversity in the counterion-specific response to the applied electric field. Finally, we demonstrate that such diverse ion distributions, along with specific electric-field-strength-dependent ion properties, lead to (1) electroosmotic (EOS) transport in nanochannels grafted with PAA brushes screened with Cs+ ions to be always counterion dominated, (2) EOS transport in nanochannels grafted with PAA brushes screened with Ca2+ and Y3+ ions to be always coion-dominated, and (3) EOS transport in nanochannels grafted with PAA brushes screened with Na+ and Ba2+ ions to be coion dominated for smaller electric fields and counterion dominated for larger electric fields.
Collapse
Affiliation(s)
- Turash Haque Pial
- Department of Mechanical Engineering, University of Maryland, College Park, Maryland20742, United States
| | - Siddhartha Das
- Department of Mechanical Engineering, University of Maryland, College Park, Maryland20742, United States
| |
Collapse
|
8
|
Theeyancheri L, Sahoo R, Kumar P, Chakrabarti R. In Silico Studies of Active Probe Dynamics in Crowded Media. ACS OMEGA 2022; 7:33637-33650. [PMID: 36188301 PMCID: PMC9520552 DOI: 10.1021/acsomega.2c04709] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/26/2022] [Accepted: 09/07/2022] [Indexed: 06/16/2023]
Abstract
Active systems are made of agents, each of which takes energy from the environment and converts it to directed motion. Therefore, by construction, these systems function out of equilibrium and cannot be described using equilibrium statistical mechanics. Though the most studied aspect has been the collective motion of active particles, the motion at the individual particle level in crowded media is also of prime importance. Examples include the motion of bacteria in hydrogels, single cell migration as a way to search for food or escape from toxic agents, and synthetic active agents transporting through soft crowded media. This review presents an overview of our understanding of single active probe dynamics in crowded media from computer simulations. The active probe is a Janus or a dumbbell-shaped particle, and the medium is made of crowders that are either sticky or repulsive to the probe and could be frozen or mobile. The density and the topology of the crowders also play an important role. We hope our in silico studies will help to elucidate the mechanism of activity-driven transport in crowded media in general and design nanomachines for targeted delivery.
Collapse
|
9
|
Hunter MA, Demir B, Petersen CF, Searles DJ. New Framework for Computing a General Local Self-Diffusion Coefficient Using Statistical Mechanics. J Chem Theory Comput 2022; 18:3357-3363. [PMID: 35657378 DOI: 10.1021/acs.jctc.2c00207] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Widely applicable, modified Green-Kubo expressions for the local diffusion coefficient (Dl) are obtained using linear response theory. In contrast to past definitions in use, these expressions are statistical mechanical results. Molecular simulations of systems with anisotropic diffusion and an inhomogeneous density profile confirm the validity of the results. Diffusion coefficients determined from different expressions in terms of currents and velocity correlations agree in the limit of large systems. Furthermore, they apply to arbitrarily small local regions, making them readily applicable to nanoscale and inhomogeneous systems where knowledge of Dl is important.
Collapse
Affiliation(s)
- Michelle A Hunter
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Baris Demir
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Charlotte F Petersen
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Queensland 4072, Australia
| | - Debra J Searles
- Australian Institute for Bioengineering and Nanotechnology, The University of Queensland, Brisbane, Queensland 4072, Australia.,School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland 4072, Australia
| |
Collapse
|
10
|
Goswami K, Chakrabarti R. Motion of an active particle with dynamical disorder. SOFT MATTER 2022; 18:2332-2345. [PMID: 35244134 DOI: 10.1039/d1sm01816g] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 06/14/2023]
Abstract
We propose a model for investigating the motion of a single active particle in a heterogeneous environment where the heterogeneity may arise due to crowding, conformational fluctuations and/or slow rearrangement of the surroundings. Describing the active particle in terms of the Ornstein-Uhlenbeck process (OUP) and incorporating heterogeneity in a thermal bath using two separate models, namely "diffusing diffusivity" and "switching diffusion", we explore the essential dynamical properties of the particle for its one-dimensional motion. In addition, we show how the dynamical behavior is controlled by dynamical variables associated with active noise such as strength and persistence time. Our model is relevant in the context of single particle dynamics in a crowded environment, driven by activity.
Collapse
Affiliation(s)
- Koushik Goswami
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, Powai 400076, India.
| | - Rajarshi Chakrabarti
- Department of Chemistry, Indian Institute of Technology Bombay, Mumbai, Powai 400076, India.
| |
Collapse
|
11
|
Chen Y, Yan R, Zhao N. Passive and active tracer dynamics in polymer solutions with isotropic-to-nematic phase transition. Phys Chem Chem Phys 2022; 24:7415-7429. [PMID: 35266498 DOI: 10.1039/d2cp00323f] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022]
Abstract
Macromolecular crowding plays a crucial role in determining the dynamics in a living cell. We adopt Langevin dynamics simulations to investigate the anomalous diffusion dynamics of passive and active particles in a solution of polymer chains with tunable stiffness. The solution's anisotropic feature is modulated by changing both the polymer stiffness and volume fraction, where isotropic-to-nematic phase transition is involved. Our results demonstrate the significant impact of polymer flexibility on the dynamics of both passive and active probes. The distinct diffusion mechanism for an active particle is clarified by the interplay between polymer stiffness, crowdedness and activity. Polymer stiffness leads to a global inhibition effect on passive particle diffusion. The diffusion coefficient exhibits an intriguing non-monotonic variation at increasing polymer stiffness, which is due to the fact that the alignment of polymer chains is beneficial for diffusion along the nematic direction but unfavorable for that in the direction perpendicular to it. In sharp contrast, polymer stiffness plays a dominant role in facilitating active particle diffusion. Self-propulsion of the particle can utilize stiffness-induced elastic interactions more efficiently, which promotes its mobility in both directions. Meanwhile, an active particle might have a stronger ability to take advantage of the polymer alignment, contributing substantially enhanced diffusivity. In addition, the diffusion coefficient of an active particle is subject to a tendency of degeneration against varying volume fraction. This counter-intuitive behavior is due to the contrasting factors that increasing crowdedness induces a lower particle speed but a longer persistent motion time.
Collapse
Affiliation(s)
- Ying Chen
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Ran Yan
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| | - Nanrong Zhao
- College of Chemistry, Sichuan University, Chengdu 610064, China.
| |
Collapse
|