Partner: Prof. Bogdan Cichocki, Ph.D. 

Recent publications
1.  Żuk P.J., Cichocki B.^{♦}, Szymczak P.^{♦}, GRPY: An Accurate Bead Method for Calculation of Hydrodynamic Properties of Rigid Biomacromolecules, BIOPHYSICAL JOURNAL, ISSN: 00063495, DOI: 10.1016/j.bpj.2018.07.015, Vol.115, No.5, pp.782800, 2018 Abstract: Two main problems that arise in the context of hydrodynamic bead modeling are an inaccurate treatment of bead overlaps and the necessity of using volume corrections when calculating intrinsic viscosity. We present a formalism based on the generalized RotnePragerYamakawa approximation that successfully addresses both of these issues. The generalized RotnePragerYamakawa method is shown to be highly effective for the calculation of transport properties of rigid biomolecules represented as assemblies of spherical beads of different sizes, both overlapping and nonoverlapping. We test the method on simple molecular shapes as well as real protein structures and compare its performance with other computational approaches. Affiliations:
 
2.  Żuk P.J., Cichocki B.^{♦}, Szymczak P.^{♦}, Intrinsic viscosity of macromolecules within the generalized Rotne–Prager–Yamakawa approximation, JOURNAL OF FLUID MECHANICS, ISSN: 00221120, DOI: 10.1017/jfm.2017.264, Vol.822, pp.R2111, 2017 Abstract: We develop a generalized Rotne–Prager–Yamakawa approximation for the dipolar components of the inverse friction matrix and use it for calculating the intrinsic viscosity of rigidly connected bead conglomerates. Such bead models are commonly used in the calculation of hydrodynamic properties of macromolecules. We consider both the case of nonoverlapping constituent beads as well as overlapping beads of different sizes. We demonstrate the accuracy of the approximation in two test cases and show that it performs well even if the distances between the beads are small or if the beads overlap. Robust performance of this approximation in the case of overlapping beads stems from its correct limiting behaviour at a complete overlap, with one sphere fully immersed in the other. The generalized Rotne–Prager–Yamakawa approximation is thus well suited for evaluation of intrinsic viscosity, which is a key quantity in characterizing molecular conformations of biological macromolecules. Keywords:complex ﬂuids, lowReynoldsnumber ﬂows, mathematical foundations Affiliations:
 
3.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Translational and rotational Brownian displacements of colloidal particles of complex shapes, ARCHIVES OF MECHANICS, ISSN: 03732029, Vol.69, No.3, pp.257267, 2017 Abstract: The exact analytical expressions for the timedependent crosscorrelations of the translational and rotational Brownian displacements of a particle with arbitrary shape were derived by us in [3, 4]. They are in this work applied to construct a method to analyze the Brownian motion of a particle of an arbitrary shape, and to extract accurately the selfdiffusion matrix from the measurements of the crosscorrelations, which in turn allows to gain some information on the particle structure. As an example, we apply our new method to analyze the experimental results of D. J. Kraft et al. for the micrometersized aggregates of the beads [8]. We explicitly demonstrate that our procedure, based on the measurements of the timedependent crosscorrelations in the whole range of times, allows to determine the selfdiffusion (or alternatively the friction matrix) with a much higher precision than the method based only on their initial slopes. Therefore, the analytical timedependence of the crosscorrelations serves as a useful tool to extract information about particle structure from trajectory measurements. Keywords:Brownian motion, Smoluchowski equation, hydrodynamic interactions, selfdiﬀusion matrix, friction coeﬃcients, crosscorrelations of translational and rotational Brownian displacements Affiliations:
 
4.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Note: Brownian motion of colloidal particles of arbitrary shape, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.4942119, Vol.144, pp.07610112, 2016 Abstract: The analytical expressions for the timedependent cross correlations of the translational and rotational Brownian displacements of a particle with arbitrary shape are derived. The reference center is arbitrary, and the reference frame is such that the rotationalrotational diffusion tensor is diagonal. Keywords:Rotational correlation time, Tensor methods, Brownian motion, Colloidal systems, Matrix equations Affiliations:
 
5.  Lisicki M.^{♦}, Cichocki B.^{♦}, Wajnryb E., Nearwall diﬀusion tensor of an axisymmetric colloidal particle, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.4958727, Vol.145, pp.03490418, 2016 Abstract: Hydrodynamic interactions with confining boundaries often lead to drastic changes in the diffusive behaviour of microparticles in suspensions. For axially symmetric particles, earlier numerical studies have suggested a simple form of the nearwall diffusion matrix which depends on the distance and orientation of the particle with respect to the wall, which is usually calculated numerically. In this work, we derive explicit analytical formulae for the dominant correction to the bulk diffusion tensor of an axially symmetric colloidal particle due to the presence of a nearby noslip wall. The relative correction scales as powers of inverse wallparticle distance and its angular structure is represented by simple functions in sines and cosines of the particle’s inclination angle to the wall. We analyse the correction for translational and rotational motion, as well as the translationrotation coupling. Our findings provide a simple approximation to the anisotropic diffusion tensor near a wall, which completes and corrects relations known from earlier numerical and theoretical findings. Affiliations:
 
6.  Liu Y.^{♦}, Bławzdziewicz J.^{♦}, Cichocki B.^{♦}, Dhont J.K.G.^{♦}, Lisicki M.^{♦}, Wajnryb E., Youngf Y.N.^{♦}, Lang P.R.^{♦}, Nearwall dynamics of concentrated hardsphere suspensions: comparison of evanescent wave DLS experiments, virial approximation and simulations, SOFT MATTER, ISSN: 1744683X, DOI: 10.1039/c5sm01624j, Vol.11, pp.73167327, 2015 Abstract: In this article we report on a study of the nearwall dynamics of suspended colloidal hard spheres over a broad range of volume fractions. We present a thorough comparison of experimental data with predictions based on a virial approximation and simulation results. We find that the virial approach describes the experimental data reasonably well up to a volume fraction of ϕ ≈ 0.25 which provides us with a fast and noncostly tool for the analysis and prediction of evanescent wave DLS data. Based on this we propose a new method to assess the nearwall selfdiffusion at elevated density. Here, we qualitatively confirm earlier results [Michailidou, et al., Phys. Rev. Lett., 2009, 102, 068302], which indicate that manyparticle hydrodynamic interactions are diminished by the presence of the wall at increasing volume fractions as compared to bulk dynamics. Beyond this finding we show that this diminishment is different for the particle motion normal and parallel to the wall. Affiliations:
 
7.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Brownian motion of a particle with arbitrary shape, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.4921729, Vol.142, pp.21490219, 2015 Abstract: rownian motion of a particle with an arbitrary shape is investigated theoretically. Analytical expressions for the timedependent crosscorrelations of the Brownian translational and rotational displacements are derived from the Smoluchowski equation. The role of the particle mobility center is determined and discussed. Keywords:Tensor methods, Brownian motion, Matrix theory, Rotational correlation time, Eigenvalues Affiliations:
 
8.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Hydrodynamic radius approximation for spherical particles suspended in a viscous ﬂuid: Inﬂuence of particle internal structure and boundary, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.4871498, Vol.140, pp.16490217, 2014 Abstract: Systems of spherical particles moving in Stokes flow are studied for different particle internal structures and boundaries, including the Navierslip model. It is shown that their hydrodynamic interactions are well described by treating them as solid spheres of smaller hydrodynamic radii, which can be determined from measured singleparticle diffusion or intrinsic viscosity coefficients. Effective dynamics of suspensions made of such particles is quite accurately described by mobility coefficients of the solid particles with the hydrodynamic radii, averaged with the unchanged direct interactions between the particles. Keywords:hydrodynamic radius, Stokes equations Affiliations:
 
9.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Shorttime dynamics and highfrequency rheology of suspensions of spherical coreshell particles with thinshells, COLLOIDS AND SURFACES APHYSICOCHEMICAL AND ENGINEERING ASPECTS, ISSN: 09277757, DOI: 10.1016/j.colsurfa.2012.10.053, Vol.418, pp.2228, 2013 Abstract: Shorttime dynamics and highfrequency rheology for suspensions of nonoverlapping core–shell particles with thin shells were analyzed. In the thinshell limit, the singleparticle scattering coefficients were derived and shown to define a unique effective radius. This result was used to justify theoretically (in the thinshell limit) the accuracy of the annulus approximation with the inner radius equal to the effective hydrodynamic radius of the core–shell particle. The twoparticle virial expansion of the translational and rotational selfdiffusion, sedimentation and viscosity was performed. The virial coefficients were evaluated and shown to be accurately approximated by the effective annulus model, in contrast to the imprecise effective hard sphere model. Keywords:Stokes equations, Core–shell particles, Permeable medium Affiliations:
 
10.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Naegele G.^{♦}, Wajnryb E., Diffusion, sedimentation, and rheology of concentrated suspensions of coreshell particles, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3689322, Vol.136, pp.104902116, 2012 Abstract: Shorttime dynamic properties of concentrated suspensions of colloidal coreshell particles are studied using a precise force multipole method which accounts for manyparticle hydrodynamic interactions. A coreshell particle is composed of a rigid, spherical dry core of radius a surrounded by a uniformly permeable shell of outer radius b and hydrodynamic penetration depth κ−1. The solvent flow inside the permeable shell is described by the BrinkmanDebyeBueche equation, and outside the particles by the Stokes equation. The particles are assumed to interact nonhydrodynamically by a hardsphere nooverlap potential of radius b. Numerical results are presented for the highfrequency shear viscosity, η∞, sedimentation coefficient, K, and the shorttime translational and rotational selfdiffusion coefficients, D t and D r . The simulation results cover the full threeparametric fluidphase space of the composite particle model, with the volume fraction extending up to 0.45, and the whole range of values for κb, and a/b. Manyparticle hydrodynamic interaction effects on the transport properties are explored, and the hydrodynamic influence of the core in concentrated systems is discussed. Our simulation results show that for thin or hardly permeable shells, the coreshell systems can be approximated neither by noshell nor by nocore models. However, one of our findings is that for κ(b − a) ≳ 5, the core is practically not sensed any more by the weakly penetrating fluid. This result is explained using an asymptotic analysis of the scattering coefficients entering into the multipole method of solving the Stokes equations. We show that in most cases, the influence of the core grows only weakly with increasing concentration. Keywords:coreshell particles, suspension, diffusion, sedimentation, effective viscosity Affiliations:
 
11.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Communication: Translational Brownian motion for particles of arbitrary shape, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3689842, Vol.136, pp.07110214, 2012 Abstract: A single Brownian particle of arbitrary shape is considered. The timedependent translational mean square displacement W(t) of a reference point at this particle is evaluated from the Smoluchowski equation. It is shown that at times larger than the characteristic time scale of the rotational Brownian relaxation, the slope of W(t) becomes independent of the choice of a reference point. Moreover, it is proved that in the longtime limit, the slope of W(t) is determined uniquely by the trace of the translationaltranslational mobility matrix μtt evaluated with respect to the hydrodynamic center of mobility. The result is applicable to dynamic light scattering measurements, which indeed are performed in the longtime limit. Keywords:translational and rotational Brownian motion, mean square displacement, particle of arbitrary shape, mobility center Affiliations:
 
12.  Adamczyk Z.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Słowicka A., Wajnryb E., Wasilewska M.^{♦}, Fibrinogen conformations and charge in electrolyte solutions derived from DLS and dynamic viscosity measurements, JOURNAL OF COLLOID AND INTERFACE SCIENCE, ISSN: 00219797, DOI: 10.1016/j.jcis.2012.07.010, Vol.385, pp.244257, 2012 Abstract: Hydrodynamic properties of fibrinogen molecules were theoretically calculated. Their shape was approximated by the bead model, considering the presence of flexible side chains of various length and orientation relative to the main body of the molecule. Using the bead model, and the precise manymultipole method of solving the Stokes equations, the mobility coefficients for the fibrinogen molecule were calculated for arbitrary orientations of the arms whose length was varied between 12 and 18 nm. Orientation averaged hydrodynamic radii and intrinsic viscosities were also calculated by considering interactions between the side arms and the core of the fibrinogen molecule. Whereas the hydrodynamic radii changed little with the interaction magnitude, the intrinsic viscosity exhibited considerable variation from 30 to 60 for attractive and repulsive interactions, respectively. These theoretical results were used for the interpretation of experimental data derived from sedimentation and diffusion coefficient measurements as well as dynamic viscosity measurements. Optimum dimensions of the fibrinogen molecule derived in this way were the following: the contour length 84.7 nm, the side arm length 18 nm, and the total volume 470 nm3, which gives 16% hydration (by volume). Our calculations enabled one to distinguish various conformational states of the fibrinogen molecule, especially the expanded conformation, prevailing for pH < 4 and lower ionic strength, characterized by high intrinsic viscosity of 50 and the hydrodynamic radius of 10.6 nm. On the other hand, for the physiological condition, that is, pH = 7.4 and the ionic strength of 0.15 M NaCl, the semicollapsed conformation dominates. It is characterized by the average angle equal to = 55, intrinsic viscosity of 35, and the hydrodynamic radius of 10 nm. Additionally, the interaction energy between the arms and the body of the molecule was predicted to be 4 kT units, confirming that they are oppositely charged than the central nodule. Results obtained in our work confirm an essential role of the side chains responsible for a highly anisotropic charge distribution in the fibrinogen molecule. These finding can be exploited to explain anomalous adsorption of fibrinogen on various surfaces. Keywords:Bead model of fibrinogen, Charge distribution over fibrinogen, Conformations of fibrinogen molecule, Fibrinogen molecule conformations, Hydrodynamic radius of fibrinogen, Viscosity of fibrinogen solutions Affiliations:
 
13.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Rotational and translational selfdiffusion in concentrated suspensions of permeable particles, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3604813, Vol.134, pp.24490317, 2011 Abstract: In our recent work on concentrated suspensions of uniformly porous colloidal spheres with excluded volume interactions, a variety of shorttime dynamic properties were calculated, except for the rotational selfdiffusion coefficient. This missing quantity is included in the present paper. Using a precise hydrodynamic force multipole simulation method, the rotational selfdiffusion coefficient is evaluated for concentrated suspensions of permeable particles. Results are presented for particle volume fractions up to 45% and for a wide range of permeability values. From the simulation results and earlier results for the firstorder virial coefficient, we find that the rotational selfdiffusion coefficient of permeable spheres can be scaled to the corresponding coefficient of impermeable particles of the same size. We also show that a similar scaling applies to the translational selfdiffusion coefficient considered earlier. From the scaling relations, accurate analytic approximations for the rotational and translational selfdiffusion coefficients in concentrated systems are obtained, useful to the experimental analysis of permeableparticle diffusion. The simulation results for rotational diffusion of permeable particles are used to show that a generalized StokesEinsteinDebye relation between rotational selfdiffusion coefficient and highfrequency viscosity is not satisfied. Keywords:selfdiffusion, permeable particles, concentrated suspensions Affiliations:
 
14.  Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Firstorder virial expansion of shorttime diffusion and sedimentation coefficients of permeable particles suspensions, PHYSICS OF FLUIDS, ISSN: 10706631, DOI: 10.1063/1.3626196, Vol.23, pp.08330318, 2011 Abstract: For suspensions of permeable particles, the shorttime translational and rotational selfdiffusion coefficients, and collective diffusion and sedimentation coefficients are evaluated theoretically. An individual particle is modeled as a uniformly permeable sphere of a given permeability, with the internal solvent flow described by the DebyeBuecheBrinkman equation. The particles are assumed to interact nonhydrodynamically by their excluded volumes. The virial expansion of the transport properties in powers of the volume fraction is performed up to the twoparticle level. The firstorder virial coefficients corresponding to twobody hydrodynamic interactions are evaluated with very high accuracy by the series expansion in inverse powers of the interparticle distance. Results are obtained and discussed for a wide range of the ratio, x, of the particle radius to the hydrodynamic screening length inside a permeable sphere. It is shown that for x≥10, the virial coefficients of the transport properties are wellapproximated by the hydrodynamic radius (annulus) model developed by us earlier for the effective viscosity of porousparticle suspensions. Keywords:Stokes equations, hydrodynamic interactions, diffusion, sedimentation, permeable particles, suspesnion, virial expansion Affiliations:
 
15.  Mongruela A.^{♦}, Lecoq N.^{♦}, Wajnryb E., Cichocki B.^{♦}, Feuillebois F.^{♦}, Motion of a spherocylindrical particle in a viscous fluid in confined geometry, EUROPEAN JOURNAL OF MECHANICS BFLUIDS, ISSN: 09977546, DOI: 10.1016/j.euromechflu.2011.04.005, Vol.30, pp.405408, 2011 Abstract: The motion of a millimeter size spherocylinder particle settling in a very viscous oil in a closed container is measured by laser interferometry, with the goal to model the motion of a particle of this shape in a fluid at microscales. The container is a cylinder with vertical axis and closed at both ends by horizontal plates. The displacement of the particle along the container axis is recorded with a resolution of the order of 100 nm, that is much smaller than the particle–wall separation when in the lubrication regime. The particle friction coefficient, measured as a function of the particle–wall distance, is then used to test the theoretical predictions of an accurate hydrodynamic analysis. The Stokes flow problem is solved by using the hydromultipole method, that is in general appropriate for spheres but is extended here to a nonspherical particle by using a compound of overlapping spheres. The lateral wall effect is negligible but the two parallel horizontal end plane walls are accurately taken into account. The result of the theoretical model is in good quantitative agreement with experiment for the whole settling motion of the spherocylinder, that is for any position between the walls. Keywords:Stokes flows, Suspensions, Sedimentation Affiliations:
 
16.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Shorttime dynamics of permeable particles in concentrated suspensions, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3274663, Vol.132, pp.014503117, 2010 Abstract: We study shorttime diffusion properties of colloidal suspensions of neutral permeable particles. An individual particle is modeled as a solventpermeable sphere of interaction radius a and uniform permeability k, with the fluid flow inside the particle described by the Debye–Bueche–Brinkman equation, and outside by the Stokes equation. Using a precise multipole method and the corresponding numerical code HYDROMULTIPOLE that account for higherorder hydrodynamic multipole moments, numerical results are presented for the hydrodynamic function, H(q), the shorttime selfdiffusion coefficient, Ds, the sedimentation coefficient K, the collective diffusion coefficient, Dc, and the principal peak value H(qm), associated with the shorttime cage diffusion coefficient, as functions of porosity and volume fraction. Our results cover the full fluid phase regime. Generic features of the permeable sphere model are discussed. An approximate method by Pusey to determine Ds is shown to agree well with our accurate results. It is found that for a given volume fraction, the wavenumber dependence of a reduced hydrodynamic function can be estimated by a single master curve, independent of the particle permeability, given by the hardsphere model. The reduced form is obtained by an appropriate shift and rescaling of H(q), parametrized by the selfdiffusion and sedimentation coefficients. To improve precision, another reduced hydrodynamic function, hm(q), is also constructed, now with the selfdiffusion coefficient and the peak value, H(qm), of the hydrodynamic function as the parameters. For wavenumbers qa > 2, this function is permeability independent to an excellent accuracy. The hydrodynamic function of permeable particles is thus well represented in its qdependence by a permeabilityindependent master curve, and three coefficients, Ds, K, and H(qm), that do depend on the permeability. The master curve and its coefficients are evaluated as functions of concentration and permeability. Keywords:Stokes equations, hydrodynamic interactions, selfdiffusion, sedimentation, permeable particles, suspension Affiliations:
 
17.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Highfrequency viscosity of concentrated porous particles suspensions, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3474804, Vol.133, pp.08490619, 2010 Abstract: We determine the highfrequency limiting shear viscosity in colloidal suspensions of rigid, uniformly porous spheres of radius a as a function of volume fraction and inverse porosity parameter x. Our study covers the complete fluidstate regime. The flow inside the spheres is modeled by the Debye–Bueche–Brinkman equation using the boundary condition that fluid velocity and stress change continuously across the sphere surfaces. The manysphere hydrodynamic interactions in concentrated systems are fully accounted for by a precise hydrodynamic multipole method encoded in our HYDROMULTIPOLE program extended to porous particles. A truncated virial expansion is used to derive an accurate and easytouse generalized Saitô formula for. The simulation data are used to test the performance of two simplifying effective particle models. The first model describes the effective particle as a nonporous sphere characterized by a single effective radius dependent on x. In the more refined second model, the porous spheres are modeled as spherical annulus particles with an inner hydrodynamic radius as a function of x, defining the nonporous dry core and characterizing hydrodynamic interactions, and an outer excluded volume radius a characterizing the unchanged direct interactions. Only the second model is in a satisfactory agreement with the simulation data. Keywords:Stokes flow, permeable particles, effective viscosity, lubrication, concentrated suspensions Affiliations:
 
18.  Cichocki B.^{♦}, Wajnryb E., Bławzdziewicz J.^{♦}, Dhont J.K.G.^{♦}, Lang P.R.^{♦}, The intensity correlation function in evanescent wave scattering, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, Vol.132, pp.074704112, 2010  
19.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Highfrequency viscosity and generalized Stokes–Einstein relations in dense suspensions of porous particles, JOURNAL OF PHYSICSCONDENSED MATTER, ISSN: 09538984, DOI: 10.1088/09538984/22/32/322101, Vol.22, pp.32210116, 2010 Abstract: We study the highfrequency limiting shear viscosity, η∞, of colloidal suspensions of uncharged porous particles. An individual particle is modeled as a uniformly porous sphere with the internal solvent flow described by the Debye–Bueche–Brinkman equation. A precise hydrodynamic multipole method with a full account of manyparticle hydrodynamic interactions encoded in the HYDROMULTIPOLE program extended to porous particles, is used to calculate η∞ as a function of porosity and concentration. The secondorder virial expansion for η∞ is derived, and its range of applicability assessed. The simulation results are used to test the validity of generalized Stokes–Einstein relations between η∞ and various shorttime diffusion coefficients, and to quantify the accuracy of a simplifying cell model calculation of η∞. An easytouse generalized Saitˆo formula for η∞ is presented which provides a good description of its porosity and concentration dependence. Keywords:Stokes flow, hydrodynamic interactions, permeable particles, dense suspensions, effective viscosity Affiliations:
 
20.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nägele G.^{♦}, Wajnryb E., Dynamics of permeable particles in concentrated suspensions, PHYSICAL REVIEW E, ISSN: 15393755, DOI: 10.1103/PhysRevE.81.020404, Vol.81, pp.02040414, 2010 Abstract: We calculate shorttime diffusion properties of suspensions of porous colloidal particles as a function of their permeability, for the full fluidphase concentration range. The particles are modeled as spheres of uniform permeability with excluded volume interactions. Using a precise multipole method encoded in the HYDROMULTIPOLE program, results are presented for the hydrodynamic function, H(q), sedimentation coefficient, and selfdiffusion coefficients with a full account of manybody hydrodynamic interactions. While selfdiffusion and sedimentation are strongly permeability dependent, the wavenumber dependence of the hydrodynamic function can be reduced by appropriate shifting and scaling, to a single master curve, independent of permeability. Generic features of the permeable sphere model are discussed. Keywords:Stokes equations, hydrodynamic interactions, permeable particles, concentrated suspensions, selfdiffusion, hydrodynamic function, collective diffusion Affiliations:
 
21.  Cichocki B.^{♦}, EkielJeżewska M.L., Green tensors for Debye–Büche–Brinkman equations generalized for axisymmetric medium, JOURNAL OF MATHEMATICAL PHYSICS, ISSN: 00222488, DOI: 10.1063/1.3498686, Vol.51, pp.103101112, 2010 Abstract: A fluid flow through a nonisotropic porous medium with an axial symmetry is considered. The Green tensors for the corresponding nonisotropic Debye–Büche–Brinkman equations are calculated in terms of single integrals. Shortdistance and farfield limiting behavior is discussed. The exact solution for the Green tensors is found explicitly in the limiting case of an infinite shielding length along the symmetry axis of the system. Keywords:Tensor methods, Porous media, Fluid equations, Integral equations, Exact solutions Affiliations:
 
22.  Cichocki B.^{♦}, EkielJeżewska M.L., Selfdiffusion of a sphere in an effective medium of rods, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.3146786, Vol.130, pp.21490216, 2009 Abstract: Selfdiffusion of a sphere in a network of rods is analyzed theoretically.Hydrodynamic interactions are taken into account according to the model of Dhont et al. [J. Chem. Phys.122, 044905 (2005); Dhont et al., J. Chem. Phys.124, 044907 (2006); Dhont et al., J. Chem. Phys.126, 214501 (2007)] based on the Debye–Bueche–Brinkman equation. The hydrodynamic screening length of the effective medium is assumed to be much larger than the sphere radius and the rod thickness. The selfdiffusion coefficient, given by Dhont et al. in terms of fourdimensional integrals, is in this work expressed in terms of a single integral only and therefore evaluated numerically with a high precision. Moreover, simple expressions for the selfdiffusion coefficient are derived and shown to be independent of the rod length. They can be useful for experimental verification of the model. Keywords:Hydrodynamics, Friction, Hydrological modeling, Self diffusion, Tensor methods Affiliations:
 
23.  Guzowski J.^{♦}, Cichocki B.^{♦}, Wajnryb E., Abade G.C.^{♦}, The shorttime selfdiffusion coefficient of a sphere in a suspension of rigid rods, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.2837296, Vol.128, pp.94502111, 2008  
24.  Cichocki B.^{♦}, Sadlej K., Stokesian dynamics  the BBGKY hierarchy for correlation functions, JOURNAL OF STATISTICAL PHYSICS, ISSN: 00224715, Vol.132, pp.129151, 2008  
25.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Hydrodynamic interactions between spheres in a viscous fluid with a flat free surface or hard wall, JOURNAL OF CHEMICAL PHYSICS, ISSN: 00219606, DOI: 10.1063/1.2724815, Vol.126, pp.184704184707, 2007 Abstract: Hydrodynamic interactions between spheres immersed in a lowReynoldsnumber fluid flow close to a flat free surface or hard wall are investigated. The spheres may have different or equal radii, and may be separated from the boundary or at contact with the free surface. A simple and useful expression is derived for the propagator (Green operator) connecting centers of two spheres. In the derivation, the method of images and the displacement theorems are used. Symmetry of the displacement operators is explicitly shown. The significance of these results in efficient Stokesian and Brownian dynamics simulations is outlined. An example of an application is shown. Keywords:Free surface, Hydrodynamics, Friction, Fluid flows, Mirrors Affiliations:
 
26.  May S.^{♦}, Tokarzewski S., Zachara A., Cichocki B.^{♦}, Efektywna przewodność dwuskładnikowego kompozytu o regularnej dwuwymiarowej strukturze, Prace IPPT  IFTR Reports, ISSN: 22993657, No.24, pp.178, 1992 
List of chapters in recent monographs
1. 305  Cichocki B.^{♦}, EkielJeżewska M.L., Naegele G.^{♦}, Wajnryb E., Drops and Bubbles in Contact with Solid Surface, rozdział: Hydrodynamic Interactions Between Solid Particles at a FluidGas Interface, Ferrari M., Liggieri L., Miller R., CRC Press, Leiden, pp.93104, 2012 
Conference papers
1.  Cichocki B.^{♦}, EkielJeżewska M.L., Wajnryb E., Intrinsic viscosity for Brownian particles of arbitrary shape, Symposium on Microparticles in Stokes Flows in Honor of Francois Feuillebois 65th Birthday, 20110821/0824, Warszawa (PL), DOI: 10.1088/17426596/392/1/012004, Vol.392, pp.01200417, 2012 Abstract: General expressions for the frequencydependent Brownian contribution to the intrinsic viscosity of arbitraryshaped particles have been derived from the Smoluchowski equation. Keywords:intrinsic viscosity, Brownian motion, particle of arbitrary shape Affiliations:

Conference abstracts
1.  EkielJeżewska M.L., Cichocki B.^{♦}, Wajnryb E., Translational and rotational Brownian motion of particles of complex shapes, 30th Marian Smoluchowski Symposium on Statistical Physics, 20170903/0908, Kraków (PL), pp.12, 2017  
2.  EkielJeżewska M.L., Cichocki B.^{♦}, Wajnryb E., Brownian motion of colloidal particles with arbitrary shapes, Dynamics of Interfaces in Complex Fluids and Complex Flows, 20170228/0303, Erlangen (DE), pp.2021, 2017  
3.  Abade G.C.^{♦}, Cichocki B.^{♦}, EkielJeżewska M.L., Nagele G.^{♦}, Wajnryb E., Diffusion, sedimentation, and rheology of concentrated suspensions of coreshell particles, III National Conference of Nano and Micromechanics, 20120704/0706, Warszawa (PL), pp.7980, 2012 Abstract: Shorttime dynamic properties of concentrated suspensions of colloidal coreshell particles have been recently studied [1] using a precise force multipole method which accounts for manyparticle hydrodynamic interactions (HIs). A coreshell particle is composed of a rigid, spherical dry core of radius a surrounded by an uniformly permeable shell of outer radius b and hydrodynamic penetration depth κ1. The solvent flow inside the permeable shell is described by the BrinkmanDebyeBueche equation, and outside the particles by the Stokes equation. The particles are assumed to interact nonhydrodynamically by a hardsphere nooverlap potential of radius b. Numerical results are presented for the highfrequency shear viscosity, sedimentation coefficient and the shorttime translational and rotational selfdiffusion coefficients. The simulation results cover the full threeparametric fluidphase space of the composite particle model, with the volume fraction extending up to 0.45, and the whole range of values for κb, and a/b. Manyparticle hydrodynamic interaction effects on the transport properties are explored, and the hydrodynamic influence of the core in concentrated systems is discussed. Keywords:Stokes equations, BrinkmanDebyeBueche equations, permeable particles, translational and rotational selfdiffusion, sedimentation, effective viscosity Affiliations:
