Next Article in Journal
Molecular Docking-Based Virtual Screening of FDA-Approved Drugs Using Trypanothione Reductase Identified New Trypanocidal Agents
Next Article in Special Issue
How Doping Regulates As(III) Adsorption at TiO2 Surfaces: A DFT + U Study
Previous Article in Journal
Synthesis of Coumarin-Based Photosensitizers for Enhanced Antibacterial Type I/II Photodynamic Therapy
Previous Article in Special Issue
ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel

1
Institute of Theoretical Physics and Mark Kac Center for Complex Systems Research, Jagiellonian University, ul. St. Łojasiewicza 11, 30-348 Kraków, Poland
2
Faculty of Chemistry, Silesian University of Technology, Strzody 9, 44-100 Gliwice, Poland
3
Department of Physics and Astronomy, University of California, Irvine, CA 92697, USA
*
Author to whom correspondence should be addressed.
Molecules 2024, 29(16), 3795; https://doi.org/10.3390/molecules29163795
Submission received: 30 June 2024 / Revised: 23 July 2024 / Accepted: 26 July 2024 / Published: 10 August 2024
(This article belongs to the Special Issue Feature Papers in Computational and Theoretical Chemistry)

Abstract

:
The transport of molecules and particles through single pores is the basis of biological processes, including DNA and protein sequencing. As individual objects pass through a pore, they cause a transient change in the current that can be correlated with the object size, surface charge, and even chemical properties. The majority of experiments and modeling have been performed with spherical objects, while much less is known about the transport characteristics of aspherical particles, which would act as a model system, for example, for proteins and bacteria. The transport kinetics of aspherical objects is an especially important, yet understudied, problem in nanopore analytics. Here, using the Wiener process, we present a simplified model of the diffusion of rod-shaped particles through a cylindrical pore, and apply it to understand the translation and rotation of the particles as they pass through the pore. Specifically, we analyze the influence of the particles’ geometrical characteristics on the effective diffusion type, the first passage time distribution, and the particles’ orientation in the pore. Our model shows that thicker particles pass through the channel slower than thinner ones, while their lengths do not affect the passage time. We also demonstrate that both spherical and rod-shaped particles undergo normal diffusion, and the first passage time distribution follows an exponential asymptotics. The model provides guidance on how the shape of the particle can be modified to achieve an optimal passage time.

1. Introduction

Advances in the field of ion channels in biological membranes, as well as nanopores in solid-state and polymer films, have led scientists to develop an interest in single-molecule and single-particle transport in nanoconfinement [1,2,3,4,5,6,7,8,9,10]. In the resistive pulse technique, the passage of single molecules or particles through a pore can be readily identified as a transient change in the transmembrane current, called the resistive pulse, whereby the amplitude and duration that can be linked to the object size and surface charge [11,12]. The approach provides a high throughput technique to analyze the suspensions of particles and molecules, including mixtures. The application of thin films and pores with a narrow constriction zone enabled extending the resistive pulse technique to DNA sequencing [13,14,15,16,17,18,19,20] and, more recently, protein sequencing [21,22,23], such that current blockage caused by individual nucleotides and amino acids can be identified and differentiated.
In the context of nanopore-based biosensors, there has been a significant focus on the regulation of the transit time of objects through the pore. On the one hand, the technique needs to offer a high throughput characterization of hundreds or even thousands of molecules or particles per second [24,25,26]. On the other hand, however, each molecule needs to spend a sufficient amount of time within the pore to enable detection and characterization. Yet, the majority of resistive pulse models have focused on spherical particles. In this study, we present a model system for the passage of aspherical and spherocylindrical objects, which can be applied to the analysis of specific types of proteins [27], bacteriophages [28], and bacteria [29,30], as well as synthetic particles used, e.g., in drug delivery [31].
A limited number of reports have described the passage of aspherical, rod-like particles through single pores. An early report by Golibersuch [32] provides analytical formulas that predict the magnitude of the current drop when ellipsoidal objects are in a pore. The passage of aspherical proteins through solid-state nanopores allowed researchers to ascertain both the protein shape and the dipole moment [27,33]. Resistive pulse experiments with polymer nanopores with an undulating opening diameter and silica rod-shaped particles tested the early predictions, and demonstrated significant rotations of the particles within the pore [34]. The passage of bacteria of various shapes through silicon nitride pores demonstrated how each type of bacterium could be distinguished based on its resistive-pulse characteristics [30]. Furthermore, the detection of gold nanorods in polymer and silicon nitride nanopores has been documented [35,36].
In this manuscript, we describe the kinetics of aspherical particles as they pass through a single cylindrical pore. We focus on the properties of diffusion as a function of particle size and shape in relation to the pore dimensions. The model we present describes the process of random walk performed by individual particles in a pore. Section 3 describes in detail the model we developed, together with its assumptions and the simplifications applied. This research is inspired by the experimental studies of diffusion [37], performed with submicron (400 nm in diameter) polystyrene particles passing through a single polymer micropore (11 μ m in length and 0.7 μ m in diameter). By balancing the forces acting on the particles, the authors were able to observe the particles’ random walk and determine their diffusion coefficient and electrokinetic velocity.
To describe the diffusion of particles in a pore, we use the Wiener process [38]. The methodology we present is closely connected to the concept of random walks [39,40], which describe the manner in which a particle forms its path from a multitude of random steps. The random walk approach has been used not only to explore diffusion processes [41,42], but also in the context of finance [43], search strategies [44,45], physiology [46], seismology [47], and numerous other systems. In the case of diffusion, the treatment based on the random walk models, along with appropriate extensions, allows for the distinction between normal, sub-, and superdiffusion when the movement is carried by a specific force [48,49] or through a crowded environment [50,51]. Therefore, although this approach is highly simplified, it offers significant insight into the transport properties of confined systems, which are of interest in this context.
In order to accurately describe the diffusion of anisotropic particles in a viscous medium, it is essential to consider the particle orientation in relation to the direction of motion [52,53]. Moreover, for particles that lack spherical symmetry, one must also consider the rotational degrees of freedom, which will play a significant role in the diffusive transport process. In the context of the aforementioned model, we consider the diffusive transport of single aspherical particles in a cylindrical pore. The nanopore has a diameter of one unit and a length of L. The particle is a spherocylinder of diameter d and length l, and is characterized by its length-to-width ratio of f = ( l + d ) / d = 1 + l / d (see Figure 1).
In order to simplify the model, it is common to numerically solve the Langevin equation of a traced particle using a set of effective, macroscopic parameters. For example, multiple collisions with the medium particles are modeled by a random force of specific variance and autocorrelation, whereas the energy dissipation is described by viscosity. Furthermore, these parameters can be combined into diffusion constants (see, e.g., Refs. [54,55,56,57]). Here, we employ a more simplified model, with the objective of developing a framework for simulating a random process that enables a molecule to traverse a channel. This approach allows for a significant reduction in computational time, provides a description that is independent of the microscopic details of the system, and allows for the identification of the fundamental physics of the process. In particular, our model considers the so-called overdamped limit, in which the subsequent positions and orientations of the particle can be described by the Wiener process [38]. This methodology has recently been justified both mathematically [58] and experimentally [59], and has been applied to other problems partially related to the diffusion of anisotropic particles [60].

2. Results and Discussion

2.1. Diffusion Type

The diffusion process through the channel was analyzed in terms of the mean squared displacement (see Appendix A):
x ( t ) x ( 0 ) 2 = x 2 ( t ) = 2 D × t α ,
where x ( t ) denotes the position of the particle’s center after t jumps [61] ( x ( 0 ) = [ 0 , 0 , 0 ] ), D is the diffusion constant, and the exponent α determines the diffusion type (see Appendix A). Note that, for the normal diffusion, α = 1 , while for sub- or superdiffusion, α is smaller or greater than 1, respectively [41,62].
The dependence of x 2 ( t ) on time was recorded during each simulation, and the parameters D and α were determined. A few illustrative examples of the mean squared displacement (MSD) versus time are presented in Figure 2.
As shown in Figure 2, the mean squared displacement obeys the scaling with time given by Equation (1), especially at intermediate times. Moreover, as the particle diameter r increases, the individual MSD curves become smoother (see Figure 2d) due to hindered rotation of thicker particles. The deviations from the scaling observed at the beginning and the end of the simulation are due to boundary conditions. To avoid the effects of reflecting and absorbing ends of the channel, the fitting of the MSD was performed for the time interval [ 0.01 t m i n , t m i n ] , where t m i n denotes the minimal first passage time observed in the simulation of 5000 independent passing events. The fitted values of α and D are shown in Figure 3.
The diffusion constant D for different particle sizes was found to vary between 0.001 and 0.0024 . In general, the diffusion coefficient is higher for smaller and less-elongated particles, which agrees with the intuition supported by the Stokes–Einstein–Smoluchowski–Sutherland relation [63,64,65], as well as previous observations [33,37,66]. The observed fluctuations and irregularities in Figure 3 are related to the fitting methodology. To confirm this, we conducted an additional experiment in which the number of independent trajectories that are averaged in Equation (1) was increased. This resulted in a notable reduction in the observed fluctuations and irregularities. On the other hand, the exponent α that defines the diffusion type for all particles was found to be very close to 1. This observation suggests that the particles perform normal diffusion, regardless of the particle length or thickness, which is in line with the findings reported in [59]. Our simulations, therefore, revealed that although the channel boundaries restrict particle motion by the introduction of trapping events, they do not slow down the transport sufficiently to change the normal diffusion into subdiffusion.

2.2. First Passage Time

Another important characteristic of the diffusive motion is the first passage time, defined as the shortest interval of time required for a particle to leave the pore for the first time. Its distribution is typically given by the exponential law, as follows:
p ( t ) A exp ( λ t ) ,
The distribution allowed us in turn to determine the survival probability S ( t ) :
S ( t ) = 1 0 t p ( t ) d t exp ( λ t ) .
The distribution S ( t ) quantifies the fraction of molecules that remain within the channel at a given time t. The exponent λ defines the rate of the probability decay, which reflects the diminishing probability of finding a tracer in a channel after a specified time t. We determined the first passage time for each particle, and the obtained distribution of survival probability is shown in Figure 4.
The recorded exponential tails (long-time asymptotics) indicate that this process is Markovian, i.e., the imposed geometrical constraints do not introduce memory to the system. Consequently, the next position and orientation of a particle depend only on its current state, regardless of earlier history, which is the case for the standard Wiener process in an empty space.
Our simulations revealed that the parameter λ does not depend on the particle size (see Figure 5). This result supports the previous observation, based on the diffusion exponent α , that the character of the motion does not noticeably depend either on particle length or its thickness.
Conversely, the median of the first passage time shown in Figure 6 indicates a slowdown for thicker particles, despite the almost identical tails of the survival probability for differently shaped particles. This is evident in Figure 5 and Figure 6. This apparent contradiction originates from the exponential tails present in S ( t ) only at long times. At short times t 0 , the survival probability exhibits different dependence, because a particle cannot leave the channel instantaneously. It is thus possible to determine a threshold time, after which an exponential decay of S ( t ) is visible. This time is directly related to the minimal time needed for a particle to pass through the channel. On the other hand, the equivalent information is carried by the median of passing times M t , which is a physical quantity characterizing the transport process (see Figure 6).
Interestingly, there is no additional dependence of the median time M t on the particle length l. It is important to highlight that this phenomenon is not caused by considering rotational and translational movements separately. We checked the scenario in which both translation and rotation were considered simultaneously at each step and rejected if only one of them or both lead to a collision with the channel walls. The obtained results did not show significant differences from the results presented here.
This section showed that normal Markovian diffusion was observed for all particles we considered. The diffusion coefficient is the largest for small spherical particles. At the same time, the first passage time was found to depend on particle thickness only, such that thinner objects pass through the pore faster regardless of their elongation.

2.3. Angle Distribution

The last parameter we analyzed is the particle orientation at the exit time, θ , determined with respect to the long axis of the channel. Taking into account the specific angle distribution on a sphere (see details in Appendix B) the obtained distribution of θ (see Figure 7) confirm that all orientations of a spherical particle ( f = 1 ) occur at equal probabilities. As the aspect ratio of the particle increases, the distribution p ( θ ) becomes non-uniform. As expected, the longest particles we considered align nearly parallel to the channel axis.
To further characterize the angle distribution, we analyze its second moment
σ = 1 π 0 π θ 2 p ( θ ) d θ 1 2 .
In the case of a uniform distribution, the standard deviation is given by σ = π 2 3 1.814 . The measured values of σ for different particles are presented in Figure 8.
Here, as expected, we see a rapid decrease in σ with the increase in the particle length, reflecting the decaying ability of a particle to rotate in the restricted geometry of the channel.

3. Materials and Methods

We assume that the particle’s translations Δ x and rotations Δ φ A ^ are selected randomly. Here, A ^ = A / | A | denotes the unit vector in the direction of A that represents the axis of rotation. The particle’s position and orientation are updated step-by-step according to (5) if and only if the update does not lead to the intersection of the particle’s surface with the channel’s boundaries
Δ x = σ x ξ , Δ φ A ^ = σ φ ( A ^ ) η .
In Equation (5), σ x and σ φ ( A ^ ) denote a set of positive, scalar parameters characterizing random motion. These parameters are controlled by temperature, which measures the strength of random movements, and viscosity, which dissipates the kinetic energy. The components of the displacement vector Δ x and the scalar η that determines the rotation angle are sampled from the normal distribution with a mean value of zero and unit variance. For the translational motion, we assume the same diffusion coefficient in the direction along the particle axis and in the direction that is perpendicular to the axis. The model can also be generalized to take into account possible differences in D in the two directions [52,53]. It is important to note that in the described model, particle rotation and translation properties are directly connected with the distribution of random variables, ξ and η . Thus, the observed diffusion coefficient reflects the width of these distributions and does not depend on particle size or mass, as long as the used distribution is constant. Note that this assumption would not be valid in the case of diffusion in unbounded physical systems. In the case described here, however, the effective diffusion is affected by the geometric constraints because the particle cannot cross the boundary of the pore walls. Consequently, this approach provides understanding of how diffusion properties depend on the relative sizes between the pore and the particles, which we study here.
Specifically, as the particle length increases, its motion becomes restricted due to the geometric confinement provided by the pore. This confinement effectively limits the particle’s ability to move across the pore unless its axis is oriented almost parallel to the channel axis. The direction of the rotation axis A ^ ( | A ^ | = 1 ) is also selected randomly and uniformly on the two-dimensional sphere [58]. The rotation axis is assumed to pass through the geometric origin of the particle. The value of the parameter σ φ ( A ^ ) was adjusted to satisfy the energy equipartition law, i.e., the energy is uniformly distributed among all degrees of freedom. Thus, we use
σ φ ( A ) = σ x A ^ T M A ^ 1 2 ,
where A ^ T M A ^ is the moment of inertia around axis A ^ (see [67] and Appendix C). In Equation (6), M is the inertia tensor of the spherocylinder ( r = d / 2 ):
M = M x x 0 0 0 M y y 0 0 0 M z z ,
where
M x x = M y y = m 1 1 12 l 2 + 1 4 r 2 + 2 m 2 2 5 r 2 + 1 4 l 2 + 3 8 l r ,
and
M z z = 1 2 m 1 + 4 5 m 2 r 2 .
The z-axis is parallel to the spherocylinder axis, while the x and y axes are perpendicular to it. x, y and z-axes cross at the geometric center of the spherocylinder. m 1 and m 2 are masses of the cylinder and the single hemisphere, respectively. We assume the particles have uniform density, thus masses m 1 and m 2 are related to the characteristic dimensions of the spherocylinder
m 1 m 2 = 3 l 2 r .
The position of the geometric center of the particle is described using Cartesian coordinates, wherein the z-axis is parallel to the pore axis, as previously mentioned. The pore entrance is located at z = 0 , while the exit is situated at z = L . Initially, the particle is aligned along the pore axis with its center at x = [ 0 , 0 , 0 ] . The simulation of each particle ends when the particle exits the channel, i.e., when z L . In the pore, the particle is free to move in any direction, as long as its center remains between z = 0 and z = L . The particle is, however, allowed to leave the channel only through z = L . In order to ensure that a spherocylinder does not leave the channel via the entry at z = 0 , the noise-induced displacement leading to the negative value of the z coordinate is rejected. Therefore, there is a reflecting boundary at z = 0 , and an absorbing boundary at z = L .
The random walk of each particle within the pore is modeled according to the following steps:
  • Draw a random translation vector with coordinates selected from independent N ( 0 , σ x ) distributions;
  • Update the particle position until it does not lead to a collision with channel boundaries, or leaving the particle through the pore entry ( z < 0 );
  • Select a random axis of rotation;
  • Calculate the particle’s moment of inertia with respect to rotations around the selected axis and the parameter σ φ ;
  • Draw random angle from the N ( 0 , σ φ ) distribution;
  • Update the particle orientation if it does not lead to a collision with channel boundaries;
  • Increase time;
  • Verify if the particle is still within the channel;
  • If a particle reaches the exit end of the channel ( z > L ), record the first passage time (stopping time) and finish the simulation. Otherwise, return to the first step.
In our simulation, the time unit is defined as the number of random steps, wherein each step consists of one translation and one rotation. Note that the number of steps also includes the rejected jumps.
Our approach of observing subsequent steps and rotations allows us to determine the type of diffusion a given particle follows. During a single simulation, we record the particle position from which the mean squared displacement (MSD) is calculated as a function of the number of steps. We also measure the number of steps a particle needs to leave the channel, the first passage time, and the histogram of the particle orientations inside the channel.
To obtain satisfactory statistics for each specific particle size and channel size, we performed 5000 independent numerical simulations. The parameter σ x was set to 0.05 , which was selected to optimize the time required for a single simulation while ensuring sufficient sampling of traced parameters. In the simulations, the values of length-to-width ratios are chosen to match the experimental conditions [34]. The channel’s length was set to 14.3 , and the width played the role of the unit length. The particle dimensions vary between 1 f 5.6 and 0.03 d 0.9 . The limiting case of f = 1 ( l = 0 ) corresponds to a spherical particle, because the rod connecting two hemispheres is absent in that situation (c.f. Figure 1).

4. Conclusions

The numerical simulations of particles in a channel permitted a systematic examination of the impact of particle shape on diffusion processes within a restricted geometry. The motion of the particles was modeled by the Wiener process, accounting for both the translational and the rotational degrees of freedom. Our analysis demonstrated that all particles that we considered, irrespective of their shape, followed normal diffusion. This finding indicates that the confinement of the particles and their interactions with the channel walls in our system were not able to transform the normal diffusion into subdiffusion. However, the diffusion constant is sensitive to the particle size and shape such that smaller and less elongated molecules are characterized by larger diffusion coefficients. We also found that thicker molecules are characterized by a larger median of the first passage time compared to thinner particles. Surprisingly, however, other diffusion properties, such as the decay of the survival probability, were found the same for all particles. This result might seems counter-intuitive, but it is consistent with the constraints that we imposed in the model. Namely, the same model of random motion was applied to all particles. Finally, as expected, the rotation of longer particles was found to be limited by the channel boundaries. Consequently, the distribution of the direction of the particles as they pass through the channel depends strongly on the particle size.
The numerical study of the noise-driven motion of anisotropic particles in a channel provides a more comprehensive understanding of the interrelationships between the macroscopic properties of the observed diffusive motion and the microscopic dynamics of individual particles. The ability to determine the influence of particle shape on the diffusion within a pore will be beneficial in the analysis of experimental data from biological and artificial nanopores, as well as in the optimization of particle shape for sensing applications. Further studies on the role of shape in the diffusion of single molecules, e.g., viruses, particles and cells, could boost the development of the resistive pulse technique for the identification of aspherical objects from mixtures. The model we have developed could be further modified to account for anomalous diffusion observed in some systems [68,69,70]. Modifications would include other distributions or random displacements, such as Cauchy or other heavy-tailed distributions, external forces or flows, or the inclusion of obstacles [49].

Author Contributions

Conceptualization, all authors; methodology, all authors; software, M.C.; validation, all authors; formal analysis, M.C., A.S. and M.K.; investigation, all authors; resources, all authors; data curation, M.C. and A.S. and M.K.; writing—original draft preparation, M.C., B.D., A.S. and M.K.; writing—review and editing, all authors; visualization, M.C., A.S. and M.K.; supervision, B.D. and Z.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially funded by the Polish National Agency for Academic Exchange of Poland, grant no. PPI/APM/2018/1/00004/U/001, under the Strategic Partnerships program, and by a grant from the Priority Research Area DigiWorld under the Strategic Programme Excellence Initiative at Jagiellonian University.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The numerical data are available at: https://doi.org/10.57903/UJ/QRHLMI. More information will be provided upon reasonable request.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A. Normal and Anomalous Diffusion

The trajectory x ( t ) of Brownian motion can be described in terms of the following stochastic process:
d x d t = ξ ( t ) ,
where ξ ( t ) is the Gaussian white noise ( ξ ( t ) = 0 and ξ ( t ) ξ ( s ) = 2 D δ ( t s ) ). The formal solution of Equation (A1) is
x ( t ) = 0 t ξ ( s ) d s ,
where, without loss of generality, we assumed x ( 0 ) = 0 . The solution x ( t ) is given by the Wiener process [38]. From Equation (A2), one can calculate
x 2 ( t ) = 0 t 0 t ξ ( u ) ξ ( v ) d u d v = 0 t 0 t ξ ( u ) ξ ( v ) d u d v = 0 t 0 t 2 D δ ( u v ) d u d v = 0 t 2 D d u = 2 D t
from which it can be seen that x 2 ( t ) t 1 . For x ( 0 ) = 0 the x 2 ( t ) is equivalent to the mean squared displacement. Equation (A3) shows that for Brownian motion, the mean squared displacement scales linearly in time. Equation (A3) can be generalized into higher number of dimensions
x ( t ) x ( 0 ) 2 = 2 D × t .
For the sake of simplicity, we have omitted the additional multiplicative constant that would otherwise be required due to the increased number of dimensions.
The Brownian motion is a classic example of the so-called normal diffusion. It is known, however, that the MSD does not always scale linearly with time [71,72,73]. Therefore, we can assume that
x ( t ) x ( 0 ) 2 = 2 D × t α ,
where α is the exponent defining the diffusion type. The reference case of α = 1 corresponds to the normal diffusion, while α < 1 represents subdiffusion. Finally, α > 1 stands for superdiffusion. These names indicate whether the MSD scales with time faster (superdiffusion) or slower (subdiffusion) than in the case of normal diffusion. Subdiffusion is typically produced by trapping events, whereas superdiffusion is due to anomalously long jumps [41,42,73].

Appendix B. Renormalization of the Angle Distribution

Let us assume that we measure the angle θ between a vector pointing from the center of the unit sphere to a random (uniformly distributed) point on the sphere and the channel axis. Points on the sphere are selected from the uniform distribution on the unit sphere
p ˜ ( θ , φ ) = 1 4 π ,
where θ ( 0 θ π ) is the polar angle, i.e., the angle between the particle axis and the channel axis when the molecule leaves the channel. φ ( 0 φ 2 π ) stands for the azimuthal angle, thus the angle between the orthogonal projection of the rotation axis onto the plane perpendicular to the particle axis and a fixed direction on that plane. The p ˜ ( θ ) distribution can be calculated as
p ˜ ( θ ) = 0 2 π p ( θ , φ ) sin φ d φ = 1 2 sin θ .
The density of p ˜ ( θ ) , as seen in Equation (A7), is non-uniform. The histogram H ( θ ) of the angle θ can be re-normalized as
p ( θ ) = H ( θ ) b b e b p ( θ ) d θ = 2 H ( θ ) cos ( b b ) cos ( e b ) ,
where b b and e b indicate the beginning and end of the bin containing the angle θ and H ( θ ) is the probability, calculated from the histogram, of observing a given value of θ . The re-normalized distribution p ( θ ) (for a spherical molecule) is uniform (see Figure 7).

Appendix C. Energy of Rotating (Rigid) Molecule

A moving, rigid particle’s kinetic energy K consists of two components: the translational motion of the center of mass ( K T ) and the rotation along the center of mass ( K R ). We will focus first on the rotational part only.
K R = 1 2 ω T M ω ,
where ω = [ ω x x , ω y y , ω z z ] and M is the tensor of inertia. Within the studied model, instead of generating three angles, we select a random axis of rotation A , and the particle is then rotated about this axis. Consequently, the angular velocity ω = φ ˙ A ^ (where A ^ is the unit vector indicating the rotation axis) must be transformed to the principal axes (see Equations (7)–(9)). In the principal axes basis,
ω T = φ ˙ [ A ^ · x ^ , A ^ · y ^ , A ^ · z ^ ] = φ ˙ [ A ^ x , A ^ y , A ^ z ]
and the rotational energy reads
K R = 1 2 φ ˙ 2 [ A ^ x , A ^ y , A ^ z ] M [ A ^ x , A ^ y , A ^ z ] T = 1 2 φ ˙ 2 A ^ T M A ^ ,
showing that A ^ T M A ^ can be interpreted as the inertia moment along the A ^ axis.
The kinetic energy of the translational motion is given as:
K T = 1 2 m v T v ,
where v = v v ^ is the particle velocity. Here, v = | v | and v ^ = v / v . The energy equipartition K T = K R leads to
φ ˙ 2 A ^ T M A ^ = m v 2 .
Therefore, the energy equipartition establishes a relationship between φ ˙ 2 and v 2 , and the model parameters σ x and σ ϕ ( A ^ ) are linked by Equation (6).

References

  1. Hille, B. Ion Channels of Excitable Membranes; Sinauer Associates: Sunderland, MA, USA, 2001. [Google Scholar] [CrossRef]
  2. Xue, L.; Yamazaki, H.; Ren, R.; Wanunu, M.; Ivanov, A.P.; Edel, J.B. Solid-state nanopore sensors. Nat. Rev. Mater. 2020, 5, 931–951. [Google Scholar] [CrossRef]
  3. Hou, X.; Guo, W.; Jiang, L. Biomimetic smart nanopores and nanochannels. Chem. Soc. Rev. 2011, 40, 2385–2401. [Google Scholar] [CrossRef] [PubMed]
  4. Akhtarian, S.; Miri, S.; Doostmohammadi, A.; Brar, S.K.; Rezai, P. Nanopore sensors for viral particle quantification: Current progress and future prospects. Bioengineered 2021, 12, 9189–9215. [Google Scholar] [CrossRef]
  5. Schmeltzer, A.J.; Peterson, E.M.; Harris, J.M.; Lathrop, D.K.; German, S.R.; White, H.S. Simultaneous Multipass Resistive-Pulse Sensing and Fluorescence Imaging of Liposomes. ACS Nano 2024, 18, 7241–7252. [Google Scholar] [CrossRef] [PubMed]
  6. Young, Y.W.; Kappler, M.P.; Call, E.D.; Brown, Q.J.; Jacobson, S.C. Integrated In-Plane Nanofluidic Devices for Resistive-Pulse Sensing. Annu. Rev. Anal. Chem. 2024, 17, 221–242. [Google Scholar] [CrossRef]
  7. Howorka, S.; Siwy, Z. Nanopore analytics: Sensing of single molecules. Chem. Soc. Rev. 2009, 38, 2360–2384. [Google Scholar] [CrossRef] [PubMed]
  8. Namani, S.; Kavetsky, K.; Lin, C.Y.; Maharjan, S.; Gamper, H.B.; Li, N.S.; Piccirilli, J.A.; Hou, Y.M.; Drndic, M. Unraveling RNA Conformation Dynamics in Mitochondrial Encephalomyopathy, Lactic Acidosis, and Stroke-like Episode Syndrome with Solid-State Nanopores. ACS Nano 2024, 18, 17240–17250. [Google Scholar] [CrossRef] [PubMed]
  9. Liu, K.; Pan, C.; Kuhn, A.; Nievergelt, A.P.; Fantner, G.E.; Milenkovic, O.; Radenovic, A. Detecting topological variations of DNA at single-molecule level. Nat. Commun. 2019, 10, 3. [Google Scholar] [CrossRef]
  10. Zhu, J.; Tivony, R.; Bošković, F.; Pereira-Dias, J.; Sandler, S.E.; Baker, S.; Keyser, U.F. Multiplexed Nanopore-Based Nucleic Acid Sensing and Bacterial Identification Using DNA Dumbbell Nanoswitches. J. Am. Chem. Soc. 2023, 145, 12115–12123. [Google Scholar] [CrossRef]
  11. Henriquez, R.R.; Ito, T.; Sun, L.; Crooks, R.M. The resurgence of Coulter counting for analyzing nanoscale objects. Analyst 2004, 129, 478–482. [Google Scholar] [CrossRef]
  12. Kozak, D.; Anderson, W.; Vogel, R.; Chen, S.; Antaw, F.; Trau, M. Simultaneous Size and ζ-Potential Measurements of Individual Nanoparticles in Dispersion Using Size-Tunable Pore Sensors. ACS Nano 2012, 6, 6990–6997. [Google Scholar] [CrossRef]
  13. Bull, R.A.; Adikari, T.N.; Ferguson, J.M.; Hammond, J.M.; Stevanovski, I.; Beukers, A.G.; Naing, Z.; Yeang, M.; Verich, A.; Gamaarachchi, H.; et al. Analytical validity of nanopore sequencing for rapid SARS-CoV-2 genome analysis. Nat. Commun. 2020, 11, 6272. [Google Scholar] [CrossRef] [PubMed]
  14. Goyal, P.; Goyal, P.; Krasteva, P.V.; Van Gerven, N.; Gubellini, F.; Van den Broeck, I.; Troupiotis-Tsaïlaki, A.; Jonckheere, W.; Péhau-Arnaudet, G.; Pinkner, J.S.; et al. Structural and mechanistic insights into the bacterial amyloid secretion channel CsgG. Nature 2014, 516, 250–253. [Google Scholar] [CrossRef]
  15. Wanunu, M. Nanopores: A journey towards DNA sequencing. Phys. Life Rev. 2012, 9, 125–158. [Google Scholar] [CrossRef]
  16. Bayley, H. Sequencing single molecules of DNA. Curr. Opin. Chem. Biol. 2006, 10, 628–637. [Google Scholar] [CrossRef]
  17. Cherf, G.M.; Lieberman, K.R.; Rashid, H.; Lam, C.E.; Karplus, K.; Akeson, M. Automated forward and reverse ratcheting of DNA in a nanopore at 5-Å precision. Nat. Biotechnol. 2012, 30, 344–348. [Google Scholar] [CrossRef]
  18. Deamer, D.; Akeson, M.; Branton, D. Three decades of nanopore sequencing. Nat. Biotechnol. 2016, 34, 518–524. [Google Scholar] [CrossRef]
  19. Hornblower, B.; Coombs, A.; Whitaker, R.D.; Kolomeisky, A.; Picone, S.J.; Meller, A.; Akeson, M. Single-molecule analysis of DNA-protein complexes using nanopores. Nat. Methods 2007, 4, 315–317. [Google Scholar] [CrossRef]
  20. Manrao, E.A.; Derrington, I.M.; Laszlo, A.H.; Langford, K.W.; Hopper, M.K.; Gillgren, N.; Pavlenok, M.; Niederweis, M.; Gundlach, J.H. Reading DNA at single-nucleotide resolution with a mutant MspA nanopore and phi29 DNA polymerase. Nat. Biotechnol. 2012, 30, 349–353. [Google Scholar] [CrossRef]
  21. Dorey, A.; Howorka, S. Nanopore DNA sequencing technologies and their applications towards single-molecule proteomics. Nat. Chem. 2024, 16, 314–334. [Google Scholar] [CrossRef]
  22. Brinkerhoff, H.; Kang, A.S.; Liu, J.; Aksimentiev, A.; Dekker, C. Multiple rereads of single proteins at single–amino acid resolution using nanopores. Science 2021, 374, 1509–1513. [Google Scholar] [CrossRef]
  23. Ouldali, H.; Sarthak, K.; Ensslen, T.; Piguet, F.; Manivet, P.; Pelta, J.; Behrends, J.C.; Aksimentiev, A.; Oukhaled, A. Electrical recognition of the twenty proteinogenic amino acids using an aerolysin nanopore. Nat. Biotechnol. 2020, 38, 176–181. [Google Scholar] [CrossRef]
  24. Yuan, Z.; Liu, Y.; Dai, M.; Yi, X.; Wang, C. Controlling DNA Translocation Through Solid-state Nanopores. Nanoscale Res. Lett. 2020, 15, 80. [Google Scholar] [CrossRef]
  25. Wang, C.; Sensale, S.; Pan, Z.; Senapati, S.; Chang, H.C. Slowing down DNA translocation through solid-state nanopores by edge-field leakage. Nat. Commun. 2021, 12, 140. [Google Scholar] [CrossRef]
  26. Chien, C.C.; Shekar, S.; Niedzwiecki, D.J.; Shepard, K.L.; Drndić, M. Single-Stranded DNA Translocation Recordings through Solid-State Nanopores on Glass Chips at 10 MHz Measurement Bandwidth. ACS Nano 2019, 13, 10545–10554. [Google Scholar] [CrossRef]
  27. Yusko, E.C.; Bruhn, B.R.; Eggenberger, O.M.; Houghtaling, J.; Rollings, R.C.; Walsh, N.C.; Nandivada, S.; Pindrus, M.; Hall, A.R.; Sept, D.; et al. Real-time shape approximation and fingerprinting of single proteins using a nanopore. Nat. Nanotechnol. 2017, 12, 360–367. [Google Scholar] [CrossRef] [PubMed]
  28. Dion, M.B.; Oechslin, F.; Moineau, S. Phage diversity, genomics and phylogeny. Nat. Rev. Microbiol. 2020, 18, 125–138. [Google Scholar] [CrossRef]
  29. Billaudeau, C.; Chastanet, A.; Yao, Z.; Cornilleau, C.; Mirouze, N.; Fromion, V.; Carballido-López, R. Contrasting mechanisms of growth in two model rod-shaped bacteria. Nat. Commun. 2017, 8, 15370. [Google Scholar] [CrossRef] [PubMed]
  30. Hattori, S.; Sekido, R.; Leong, I.W.; Tsutsui, M.; Arima, A.; Tanaka, M.; Yokota, K.; Washio, T.; Kawai, T.; Okochi, M. Contrasting mechanisms of growth in two model rod-shaped bacteria. Sci. Rep. 2020, 10, 15525. [Google Scholar] [CrossRef]
  31. Mitchell, M.J.; Billingsley, M.M.; Haley, R.M.; Wechsler, M.E.; Peppas, N.A.; Langer, R. Engineering precision nanoparticles for drug delivery. Nat. Rev. Drug Discov. 2021, 20, 101–124. [Google Scholar] [CrossRef]
  32. Golibersuch, D. Observation of Aspherical Particle Rotation in Poiseuille Flow via the Resistance Pulse Technique: I. Application to Human Erythrocytes. Biophys. J. 1973, 13, 265–280. [Google Scholar] [CrossRef] [PubMed]
  33. Houghtaling, J.; Ying, C.; Eggenberger, O.M.; Fennouri, A.; Nandivada, S.; Acharjee, M.; Li, J.; Hall, A.R.; Mayer, M. Estimation of Shape, Volume, and Dipole Moment of Individual Proteins Freely Transiting a Synthetic Nanopore. ACS Nano 2019, 13, 5231–5242. [Google Scholar] [CrossRef] [PubMed]
  34. Qiu, Y.; Hinkle, P.; Yang, C.; Bakker, H.E.; Schiel, M.; Wang, H.; Melnikov, D.; Gracheva, M.; Toimil-Molares, M.E.; Imhof, A.; et al. Pores with longitudinal irregularities distinguish objects by shape. ACS Nano 2015, 9, 4390–4397. [Google Scholar] [CrossRef] [PubMed]
  35. Venta, K.E.; Zanjani, M.B.; Ye, X.; Danda, G.; Murray, C.B.; Lukes, J.R.; Drndić, M. Gold Nanorod Translocations and Charge Measurement through Solid-State Nanopores. Nano Lett. 2014, 14, 5358–5364. [Google Scholar] [CrossRef] [PubMed]
  36. Maugi, R.; Hauer, P.; Bowen, J.; Ashman, E.; Hunsicker, E.; Platt, M. A methodology for characterising nanoparticle size and shape using nanopores. Nanoscale 2020, 12, 262–270. [Google Scholar] [CrossRef] [PubMed]
  37. Schiel, M.; Siwy, Z.S. Diffusion and trapping of single particles in pores with combined pressure and dynamic voltage. J. Phys. Chem. C 2014, 118, 19214–19223. [Google Scholar] [CrossRef]
  38. Gardiner, C.W. Handbook of Stochastic Methods for Physics, Chemistry and Natural Sciences; Springer Verlag: Berlin, Germany, 2009. [Google Scholar]
  39. Montroll, E.W.; Weiss, G.H. Random walks on lattices II. J. Math. Phys. 1965, 6, 167. [Google Scholar] [CrossRef]
  40. Montroll, E.W.; Shlesinger, M.F. On the wonderful world of random walks. In Nonequilibrium phenomena II: From Stochastics to Hydrodynamics; Lebowitz, J.L., Montroll, E.W., Eds.; Scientific Research Publishing: Amsterdam, The Netherlands, 1984; pp. 1–121. [Google Scholar]
  41. Metzler, R.; Klafter, J. The random walk’s guide to anomalous diffusion: A fractional dynamics approach. Phys. Rep. 2000, 339, 1–77. [Google Scholar] [CrossRef]
  42. Metzler, R.; Klafter, J. The restaurant at the end of the random walk: Recent developments in the description of anomalous transport by fractional dynamics. J. Phys. A: Math. Gen. 2004, 37, R161–R208. [Google Scholar] [CrossRef]
  43. Meerschaert, M.M.; Scalas, E. Coupled continuous time random walks in finance. Physica A 2006, 370, 114. [Google Scholar] [CrossRef]
  44. Bartumeus, F.; da Luz, M.G.E.; Viswanathan, G.M.; Catalan, J. Animal search strategies: A quantitative random-walk analysis. Ecology 2005, 86, 3078–3087. [Google Scholar] [CrossRef]
  45. Viswanathan, G.M.; Da Luz, M.G.; Raposo, E.P.; Stanley, H.E. The Physics of Foraging: An Introduction to Random Searches and Biological Encounters; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar] [CrossRef]
  46. Collins, J.J.; De Luca, C.J. Random walking during quiet standing. Phys. Rev. Lett. 1994, 73, 764. [Google Scholar] [CrossRef] [PubMed]
  47. Helmstetter, A.; Sornette, D. Diffusion of epicenters of earthquake aftershocks, Omori’s law, and generalized continuous-time random walk models. Phys. Rev. E 2002, 66, 061104. [Google Scholar] [CrossRef] [PubMed]
  48. Lacasta, A.M.; Sancho, J.M.; Romero, A.H.; Sokolov, I.M.; Lindenberg, K. From subdiffusion to superdiffusion of particles on solid surfaces. Phys. Rev. E 2004, 70, 051104. [Google Scholar] [CrossRef]
  49. Kubala, P.; Cieśla, M.; Dybiec, B. Diffusion in crowded environments: Trapped by the drift. Phys. Rev. E 2021, 104, 044127. [Google Scholar] [CrossRef] [PubMed]
  50. Cieśla, M.; Dybiec, B.; Sokolov, I.; Gudowska-Nowak, E. Taming Lévy flights in confined crowded geometries. J. Chem. Phys. 2015, 142, 164904. [Google Scholar] [CrossRef] [PubMed]
  51. Strzelewicz, A.; Krasowska, M.; Cieśla, M. Lévy flights diffusion with drift in heterogeneous membranes. Membranes 2023, 13, 417. [Google Scholar] [CrossRef] [PubMed]
  52. Edwards, S.F.; Evans, K.E. Dynamics of highly entangled rod-like molecules. J. Chem. Soc. Faraday Trans. 2 1982, 78, 113–121. [Google Scholar] [CrossRef]
  53. Tirado, M.M.; Martínez, C.L.; de la Torre, J.G. Comparison of theories for the translational and rotational diffusion coefficients of rod-like macromolecules. Application to short DNA fragments. J. Chem. Phys. 1984, 81, 2047–2052. [Google Scholar] [CrossRef]
  54. Pournin, L.; Weber, M.; Tsukahara, M.; Ferrez, J.A.; Ramaioli, M.; Liebling, T.M. Three-dimensional distinct element simulation of spherocylinder crystallization. Granular Mat. 2005, 7, 119–126. [Google Scholar] [CrossRef]
  55. Pryamitsyn, V.; Ganesan, V. Screening of hydrodynamic interactions in Brownian rod suspensions. J. Chem. Phys. 2008, 128, 134901. [Google Scholar] [CrossRef] [PubMed]
  56. Zhao, T.; Wang, X. Diffusion of rigid rodlike polymer in isotropic solutions studied by dissipative particle dynamics simulation. Polymer 2013, 54, 5241–5249. [Google Scholar] [CrossRef]
  57. Hamid, A.; Shakaib, M.; Molina, J.J.; Qasim, M.D.; Fareed, B.; Uzair, M.; Yamamoto, R. Direct numerical simulations of suspension of disk-shaped particles. Phys. Fluids 2024, 36, 043339. [Google Scholar] [CrossRef]
  58. Straube, A.V.; Höfling, F. Langevin equations and a geometric integration scheme for the overdamped limit of homogeneous rotational Brownian motion. arXiv 2024, arXiv:2403.04501. [Google Scholar]
  59. Li, D.; Liu, Y.; Luo, H.; Jing, G. Anisotropic Diffusion of Elongated Particles in Active Coherent Flows. Micromachines 2024, 15, 199. [Google Scholar] [CrossRef] [PubMed]
  60. Chaki, S.; Olsen, K.S.; Löwen, H. Anisotropic Brownian particles under resetting. arXiv 2024, arXiv:2406.17550. [Google Scholar]
  61. Sokolov, I.M.; Klafter, J. From diffusion to anomalous diffusion: A century after Einstein’s Brownian motion. Chaos 2005, 15, 026103. [Google Scholar] [CrossRef] [PubMed]
  62. Metzler, R.; Jeon, J.H.; Cherstvy, A.G.; Barkai, E. Anomalous diffusion models and their properties: Non-stationarity, non-ergodicity, and ageing at the centenary of single particle tracking. Phys. Chem. Chem. Phys. 2014, 16, 24128–24164. [Google Scholar] [CrossRef] [PubMed]
  63. Einstein, A. Über die von der Molekularkinetischen Theorie der Warme Gefordete Bewegung von in ruhenden Flussigkeiten Suspendierten Teilchen. Ann. Phys. 1905, 17, 549–560. [Google Scholar] [CrossRef]
  64. Sutherland, W. LXXV. A dynamical theory of diffusion for non-electrolytes and the molecular mass of albumin. Philos. Mag. 1905, 9, 781–785. [Google Scholar] [CrossRef]
  65. Von Smoluchowski, M. Zur kinetischen theorie der brownschen molekularbewegung und der suspensionen. Ann. Phys. 1906, 21, 756–780. [Google Scholar] [CrossRef]
  66. Tirado, M.M.; de la Torre, J.G. Rotational dynamics of rigid, symmetric top macromolecules. Application to circular cylinders. J. Chem. Phys. 1980, 73, 1986–1993. [Google Scholar] [CrossRef]
  67. Idema, T. Rotations About an Arbitrary Axis. Available online: https://phys.libretexts.org/Bookshelves/University_Physics/Mechanics_and_Relativity_(Idema)/07%3A_General_Rotational_Motion/7.03%3A_Rotations_About_an_Arbitrary_Axis (accessed on 29 June 2024).
  68. Caspi, Y.; Zbaida, D.; Cohen, H.; Elbaum, M. Anomalous diffusion of high molecular weight polyisopropylacrylamide in nanopores. Macromolecules 2009, 42, 760–767. [Google Scholar] [CrossRef]
  69. Coppens, B.; Belpaire, T.E.; Pešek, J.; Steenackers, H.P.; Ramon, H.; Smeets, B. Anomalous diffusion of nanoparticles in the spatially heterogeneous biofilm environment. iScience 2023, 26, 106861. [Google Scholar] [CrossRef]
  70. Brandani, S.; Hwang, S.; Kärger, J.; Mangano, E. Diffusion anomaly in nanopores as a rich field for theorists and a challenge for experimentalists. Nat. Commun. 2024, 15, 1–3. [Google Scholar] [CrossRef]
  71. Richardson, L.F. Atmospheric diffusion shown on a distance-neighbour graph. Proc. R. Soc. Lond. Ser. A-Contain. Pap. Math. Phys. Character 1926, 110, 709–737. [Google Scholar] [CrossRef]
  72. Jeon, J.H.; Leijnse, N.; Oddershede, L.B.; Metzler, R. Anomalous diffusion and power-law relaxation of the time averaged mean squared displacement in worm-like micellar solutions. New J. Phys. 2013, 15, 045011. [Google Scholar] [CrossRef]
  73. Weigel, A.V.; Simon, B.; Tamkun, M.M.; Krapf, D. Ergodic and nonergodic processes coexist in the plasma membrane as observed by single-molecule tracking. Proc. Natl. Acad. Sci. USA 2011, 108, 6438–6443. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Scheme of a spherocylindrical particle within a cylindrical channel. The length of the channel is designated as L, and its width is of a unit length. The length of the cylindrical portion of the particle is l, resulting in a total length is l + d . The length-to-width ratio is f = ( d + l ) / d = 1 + l / d . The angle between the channel and particle axes is denoted by the variable θ .
Figure 1. Scheme of a spherocylindrical particle within a cylindrical channel. The length of the channel is designated as L, and its width is of a unit length. The length of the cylindrical portion of the particle is l, resulting in a total length is l + d . The length-to-width ratio is f = ( d + l ) / d = 1 + l / d . The angle between the channel and particle axes is denoted by the variable θ .
Molecules 29 03795 g001
Figure 2. The mean squared displacement (MSD) x 2 ( t ) versus time for particles with different lengths and radii: (a) length l = 0 (varying diameter), (b) length l = 4.6 (varying diameter), (c) radius d = 0.03 (varying length), (d) radius d = 0.91 (varying length).
Figure 2. The mean squared displacement (MSD) x 2 ( t ) versus time for particles with different lengths and radii: (a) length l = 0 (varying diameter), (b) length l = 4.6 (varying diameter), (c) radius d = 0.03 (varying length), (d) radius d = 0.91 (varying length).
Molecules 29 03795 g002
Figure 3. The diffusion constant D and the exponent α , which defines the diffusion type as a functions of the particle’s length-to-width ratio f and diameter d. The parameters D and α were extracted from fits of Equation (1) to the data obtained from numerical simulations.
Figure 3. The diffusion constant D and the exponent α , which defines the diffusion type as a functions of the particle’s length-to-width ratio f and diameter d. The parameters D and α were extracted from fits of Equation (1) to the data obtained from numerical simulations.
Molecules 29 03795 g003
Figure 4. Sample survival probability for the particle of length-to-width ratio f = 5.2 and the diameter d = 0.34 . Blue dots correspond to simulation data, while the red line is an exponential fit. The left panel is plotted in the lin–lin scale, while the log–lin scale is used for the right one.
Figure 4. Sample survival probability for the particle of length-to-width ratio f = 5.2 and the diameter d = 0.34 . Blue dots correspond to simulation data, while the red line is an exponential fit. The left panel is plotted in the lin–lin scale, while the log–lin scale is used for the right one.
Molecules 29 03795 g004
Figure 5. The exponent λ characterizing the decay of the survival probability (see Equation (3)) for particles of length-to-width ratio f and diameter d, based on 5000 independent simulations of a single particle diffusion.
Figure 5. The exponent λ characterizing the decay of the survival probability (see Equation (3)) for particles of length-to-width ratio f and diameter d, based on 5000 independent simulations of a single particle diffusion.
Molecules 29 03795 g005
Figure 6. The Median of the first passage time as a function of the particle’s length and radius. The median was calculated from 5000 independent simulations of a spherocylindrical particle of length-to-width ratio f and diameter d passing through the channel.
Figure 6. The Median of the first passage time as a function of the particle’s length and radius. The median was calculated from 5000 independent simulations of a spherocylindrical particle of length-to-width ratio f and diameter d passing through the channel.
Molecules 29 03795 g006
Figure 7. Histogram showing the distribution of particle orientation, i.e., the angle θ between the particle axis and the channel axis, for different molecule’s lengths ( l + d = 0.45 , 0.74 , 1.31 , which correspond to f = 1 , 1.6 , 2.9 ) with the fixed diameter d = 0.45 .
Figure 7. Histogram showing the distribution of particle orientation, i.e., the angle θ between the particle axis and the channel axis, for different molecule’s lengths ( l + d = 0.45 , 0.74 , 1.31 , which correspond to f = 1 , 1.6 , 2.9 ) with the fixed diameter d = 0.45 .
Molecules 29 03795 g007
Figure 8. Parameter σ characterizing the distribution of particle orientation for different length-to-width ratios f and diameters d.
Figure 8. Parameter σ characterizing the distribution of particle orientation for different length-to-width ratios f and diameters d.
Molecules 29 03795 g008
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Cieśla, M.; Dybiec, B.; Krasowska, M.; Siwy, Z.; Strzelewicz, A. Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel. Molecules 2024, 29, 3795. https://doi.org/10.3390/molecules29163795

AMA Style

Cieśla M, Dybiec B, Krasowska M, Siwy Z, Strzelewicz A. Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel. Molecules. 2024; 29(16):3795. https://doi.org/10.3390/molecules29163795

Chicago/Turabian Style

Cieśla, Michał, Bartłomiej Dybiec, Monika Krasowska, Zuzanna Siwy, and Anna Strzelewicz. 2024. "Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel" Molecules 29, no. 16: 3795. https://doi.org/10.3390/molecules29163795

APA Style

Cieśla, M., Dybiec, B., Krasowska, M., Siwy, Z., & Strzelewicz, A. (2024). Numerical Modeling of Anisotropic Particle Diffusion through a Cylindrical Channel. Molecules, 29(16), 3795. https://doi.org/10.3390/molecules29163795

Article Metrics

Back to TopTop