red blood cells ( rbcs ) have been a model system and an inspiration for the biophysics of the membrane for decades@xcite . among the numerous works in the literature ,
many studies have been devoted to adhesion of cells or vesicles and capsules on _ flat and solid _
surfaces @xcite , where , in the case of elastic capsules only , a buckling instability may occur due to the interplay between bending and stretching energy@xcite .
buckling has also been observed in freestanding capsules due to osmotic pressure@xcite or shear forces by external flow @xcite , and in vesicles due to asymmetric lipid distributions @xcite or shear forces by external flow @xcite .
the influence of intercellular adhesion @xcite is less known despite its relevance in dense environments such as blood at physiological hematocrit or even tissues .
theoretical studies on aggregation between two rbcs based on minimizing the free energy of the membrane@xcite indicate a significant change of the geometry of contact zones between cells when varying the interaction energy . some of these shapes have been observed experimentally @xcite but not quantified as a function of the dextran or fibrinogen concentrations .
a variety of membrane shapes has been predicted theoretically for rbc doublets@xcite , depending mainly on the nondimensionalized reduced adhesion strength @xmath0 and the reduced volume @xmath1 where @xmath2 and @xmath3 are the enclosed volume and surface area of the membrane , a parameter that is indeed physiologically related to the hydration state and age of rbcs and can vary over a significant range .
examples include male - female ( convex - concave ) shapes and more or less pronounced sigmoid shapes obtained by varying @xmath4 and @xmath0 as depicted in fig .
[ 3dshapes](a ) .
the 3d reconstructions of our experimental confocal images show rbc doublets that look very similar to the theoretical predictions , too ( fig.[3dshapes](b ) ) .
experimental images also show that the amplitude of the deformation of contact zones depends non monotonically on macromolecular ( dextran ) concentration .
there are direct physiological implications of aggregation since rbcs in blood in stasis
are known to form aggregates in the form of so called _
rouleaux_. this clustering process is reversible and typical hydrodynamic shear forces in physiological flow are usually sufficient to break _ rouleaux _ , at least in larger vessels . in the 1960 s , fibrinogen could be identified as one of the main plasma proteins causing rbc aggregation or clustering @xcite .
it is known that the number and size of aggregates increase with fibrinogen concentration in the physiological range ( from approximately @xmath5 for healthy adults up to @xmath6 in acute inflammatory phases ) .
the level of rbc aggregation is used in one of the most fundamental standard hematological blood tests , the erythrocyte sedimentation rate ( esr ) which is carried out worldwide many thousand times each day .
the more and the larger the aggregates , the faster the sedimentation of the rbcs : by simply measuring visually the sedimentation front in a standardized glass capillary one gets a robust and quick , if unspecific , indication on the inflammatory state of the patient .
of course , the tendency of the rbcs to form aggregates may also increase the risk of thrombosis and cardiovascular diseases , especially in combination with stenosis .
the molecular mechanisms of macromolecular induced rbc aggregation have been the subject of controversial studies @xcite .
two models have been proposed : one is based on the physical effect of depletion @xcite and the other on physisorption or bridging @xcite .
in addition to the specific adhesion mechanism , the strength of rbc aggregation also depends on their physical properties , such as deformability @xcite , surface charge @xcite and reduced volume @xcite .
it can also be induced by other macromolecules such as dextran , a widely used molecule in laboratory experiments or as a plasma expander and in veterinary medicine .
afm based single cell force spectroscopy , rheology and sedimentation of rbcs @xcite as well as theoretical models @xcite depict that an increasing dextran concentration leads first to an increase of the interaction strength among rbcs up to a specified macromolecule concentration . beyond that maximum the interaction strength decreases again , and cell  cell aggregation vanishes .
this leads to a characteristic bell  shaped adhesion - energy versus concentration curve . for the case of fibrinogen
, it was found the interaction energy increases with the concentration , but data from the literature are typically limited to lower concentrations than for the case of dextran @xcite .     and adhesion strengths @xmath0 . from left to right : male - female shapes , strong sigmoid shapes and low sigmoid shapes.(b ) 3d - confocal imaging of rbc - doublets .
the experimental results qualitatively reproduce the theoretical predictions.[3dshapes ] ]    we present here a systematic experimental study in order to analyze quantitatively the evolution of the contact zone as a function of aggregating molecule concentration , in a large domain of parameter space by using either dextran or fibrinogen .
we analyze the case of rbc doublets as well as that of larger _
rouleaux_. a systematic numerical investigation supports the experimental results .
the numerical results show that in the case of two adhering cells , the interface remains flat up to a critical value of interaction energy above which it bifurcates into a deformed , buckled state .
our analytical model reveals a forward bifurcation with critical parameters close to those of the full numerical problem .
the analysis of the interfacial shape of two cells _ rouleaux _ , here called doublets , was performed by confocal imaging of doublets during free sedimentation , i.e. while not in contact with surfaces .
typical images are shown in fig .
[ experimentaldoublets ] for doublets in dextran 70kda solutions ( see s. i. for the fibrinogen case ) .
there is a significant evolution of the interfacial shape with polymer concentration : at low and high concentrations the interface is only weakly deformed while in the intermediate range where interaction energy is higher , the interface is strongly bent .
qualitatively , three different kinds of interfacial shapes were identified : parachute , male - female and sigmoid . in the sigmoid ( s - shape interface ) , cells were mostly laterally displaced , with a non - monotonous curvature of the contact zone , while in other shapes the two cells had a common symmetry axis .
the parachute shape was defined by at least one of the cells having a convex and a concave side .
the male - female shape was defined when the interface had a bulge , with rather flat non - interacting large membrane parts .
cells with a flat interface were also classified under this category .
similar shapes were found for fibrinogen solutions . as can be seen in fig .
[ experimentaldoublets](a ) , all shapes can be found for all polymer concentrations , probably due to dispersity of cell properties in the sample , with a probability that depends on concentration .    a quantitative analysis of the distribution of the different interfacial shapes for the different concentrations of dextran
is shown in fig .
[ experimentaldoublets](b ) , revealing an increase of the occurrence of sigmoid shapes in the intermediate range of polymer concentrations , i.e. at larger interaction energy favoring slightly larger contact areas . at low and high concentrations of dextran ,
doublets are mostly in the parachute ( about 60 % ) and male - female shapes . for fibrinogen ( fig .
[ experimentaldoublets](c ) ) , mostly parachute shapes were observed in the physiological range of the protein s concentration .
the case of echinocyte rbc doublets , obtained after sedimentation and contact with glass , was also studied ( see s. i. ) .
quantification of the interfacial shape deformation was done by fitting it with a heuristically chosen sinus function ( see s.i . ) , whose amplitude was taken as a measure of the interfacial curvature ( fig .
s2 ) . a bell shape dependency of the deformation amplitude on polymer concentration
is seen ( fig .
[ amplitudesdoublets ] ) , similar to the dependency found in afm single cell force spectroscopy @xcite or aggregation index ( see s.i . ) . for the case of dextran
, both doublets with discocytic and with echinocytic cells were analysed .
the amplitude for discocytic cells can reach much larger values than echinocytes that have a close - to - spherical shape and probably a higher membrane stiffness @xcite .      to compute the interaction energies corresponding to the experimental doublets and gain further insight on the variability of shapes we performed numerical simulations ( see s.i . for details of the method ) .
based on several successful studies ( especially equilibrium shapes and shape dynamics under external flow ) which have shown that 2d and 3d models capture the same essential features @xcite , we restrict ourselves here to a 2d model . in 2d the notion of in - plane shear elasticity looses its meaning and the vesicles and inextensible capsule model ( often evoked to model rbc ) are equivalent .    for reduced areas @xmath7 ( where @xmath3 is the area of the 2d vesicle and @xmath8 its perimeter ) in the range ( @xmath9 ) with different combinations ( cells with same or different reduced areas ) and a membrane stiffness of @xmath10 j , we varied the interaction energy @xmath11 in the range 0 - 8.23 @xmath12j / m@xmath13 .
the resulting interfacial shapes for different combinations of reduced area @xmath14 , interaction energy @xmath11 and curvature energy @xmath15 ( fig .
[ numericdoublets ] ) were in qualitative agreement with the experimentally obtained shapes , namely male - female , s - shapes and parachute shapes .
figure [ numericdoublets](a ) shows the evolution of the deformation as a function of the interaction energy . as for the analysis of experimental data we fit the interfacial shape with a sine function .
the mean amplitude of the fitted curve was compared to experimental results ( fig .
[ amplitudesdoublets ] ) , by heuristically choosing a correspondence between concentration and interaction energy : we assume a linear increase of the interaction energy with concentration up to 6.9 @xmath12j / m@xmath13 at 60 mg / ml , followed by a symmetric linear decrease , in order to get a good compromise between an agreement on the width and the maximum value of bell shaped curves . this piecewise linear variation of the interaction energy is a simplification which nevertheless yields a qualitative agreement for the overall behavior .
interestingly , the numerical results also show that the transition from a flat to a deformed interface only occurs at some finite interaction energy @xmath11 ( fig .
[ numericdoublets](b ) ) . in principle
, this method could allow the determination of the experimental interaction energy for a given concentration . however , the numerical results depend quite sensitively on the reduced areas @xmath14 of the two cells ( fig . [ numericdoublets](b ) ) , a parameter that varies slightly between cells in experiments , even in the same sample .
this dispersity is probably also the reason why the bell shape in fig .
[ amplitudesdoublets](a ) of the experimental data does not show this sharp transition .     and @xmath16 and variation of the energy [ fig - model ] ]    the evolution of the amplitude of contact zone s deformation in fig .
[ numericdoublets](b ) suggests that a bifurcation akin to buckling occurs at a critical value of the interaction energy , triggering a change of shape from flat to s - shape .
to get some insight into the involved mechanism , this behavior can actually be recovered by an analytical model through a minimization of the energy of the system .
let us suppose that the shape of an rbc in a symmetric doublet ( in 2d ) is approximated as shown in fig .
[ fig - model](a ) , with two semi - circular caps of radius @xmath17 , a roughly straight portion where curvature is negligible ( length @xmath18 ) and a contact zone with a sinusoidal shape of curvilinear length @xmath19 and amplitude @xmath20 described by the equation :    @xmath21    for @xmath22 .
the corresponding interaction energy ( per unit length in the third dimension ) is , at order 4 in @xmath20 :    @xmath23    where @xmath11 is the interaction energy per unit area .
the curvature energy of the system is concentrated in the circular caps @xmath24 and the deformed contact zone @xmath25 given by : @xmath26    the total energy @xmath27 therefore depends on @xmath20 , @xmath17 and @xmath18 which are related through the constant area @xmath3 and constant perimeter @xmath8 constraints .    in the following
we define a dimensionless energy @xmath28 , a rescaled interaction energy @xmath29 , and rescale all lengths by @xmath8 : @xmath30 and so on , and introduce the reduced area @xmath7 .
the total energy @xmath31 can then be derived and expanded at order 4 in @xmath20 ( see supplemental material ) .
figure [ fig - model](b ) shows the energy due to deformation for @xmath32 and different interaction energies @xmath33 . below a critical value @xmath34
, the minimum of energy corresponds to @xmath35 ( flat interface ) , while for @xmath36 , the energy is minimal for a finite value of @xmath37 ( see s. i. ) .
this defines the bifurcation threshold @xmath34 corresponding to @xmath38 , which is an increasing function of @xmath14 ( fig .
[ fig - model](c ) ) : @xmath39    the bifurcation diagram of the deformation amplitude ( fig .
[ fig - model](d ) ) reveals a classical supercritical bifurcation , owing to the competition between interaction energy and curvature energy through the constant perimeter and area constraints .
orders of magnitude of the bifurcation threshold and amplitude of the deformation are in good agreement with experiments and numerical simulations .
we have @xmath40 .
if we take @xmath41 m and @xmath42 j as in simulations , the critical interaction energy is @xmath43j / m@xmath13 for @xmath32 and @xmath44j / m@xmath13 for @xmath45 , which is in reasonable agreement with the simulations of fig . [ numericdoublets](b ) , while the amplitude of the deformation @xmath20 saturates at @xmath46 m in this model , rather close to the value @xmath47 m found in the simulations given the simplicity of the model .
interestingly , the energy of the equilibrium configuration decreases linearly with @xmath33 as shown in fig .
[ fig - model](d ) , with no sign of the shape bifurcation due to the smallness of the non linear terms arising from the deformation in the curvature and interaction energy terms .
this result can be of interest in the modeling of clustering dynamics in flow for instance , where the dissociation rate due to hydrodynamic stresses can be assumed to simply scale like @xmath48 despite the possibly complex shapes of contact zones .
a characterization of larger clusters was performed with _ rouleaux _ of at least 7 cells .
the fast sedimentation and arbitrary orientation of the rouleaux only allowed us to image them when settled on the ( bsa treated ) cover slip . since fibrinogen led to significant inter - connections between _ rouleaux _ which eventually led to a percolated gel like structure at 45% hematocrit , we only focused on rouleaux in linear configuration in dextran solutions ( fig .
[ imagelargecluster]a ) . at 10mg
/ ml the cells in a _
rouleau _ were typically laterally shifted and we did not observe large clusters . at higher concentrations of dextran
rouleaux _ were observed and as shown in fig .
[ imagelargecluster ] , only the first three interfaces exhibited significant deformation ( fig .
[ imagelargecluster]a ) while the interfaces far from the tip were always in a male - female shape with very low deformation , mostly below the experimental resolution limit .
again , a large variety of shapes was observed for these first three interfaces ( fig .
[ imagelargecluster]b ) .    as for doublets
, these interfaces can be sorted into different categories and the distribution of shapes is shown separately for the first , second and third interface starting from the doublet tip in fig .
[ amplitudeslargecluster]a .
unlike the case of doublets , we observed mostly male - female shapes and no parachute shapes at all , likely due to the constraints induced by the next attaching cell . for concentrations in the lower to medium - range ,
s - shape contact zones exist more frequently for the first contact zone than for the next ones , while male - female shapes are indeed prevalent everywhere except at very low dextran concentrations where rbcs tend to stay very close to their free discoid shape due to the low interaction energy , in agreement with the theoretical results on low - aggregation behavior mentioned in @xcite .
the amplitudes of the interfacial deformations ( determined by the same method as for doublets ) are shown in fig .
[ amplitudeslargecluster]a .
again , a bell shape behavior for all three interfaces was observed .
interestingly , the amplitude of the first contact zone was very similar to the values measured for echinocyte doublets ( see s.i . ) , too .
indeed , the aggregation into larger clusters gives less degree of freedom to rbcs due to higher constraints on the membrane , leading to a higher effective stiffness of rbcs .
numerical simulations were performed for 7 cells and different interaction energies @xmath11 , a constant reduced area @xmath49 and membrane stiffness @xmath10 j ( fig . [ simulationslargecluster ] ) .
the results agree qualitatively with the experimental results and the amplitudes of the interfacial deformations of the first and second interface were in good agreement with the experimental data using the same energy scale as for doublets ( see fig .
[ amplitudeslargecluster ] ) . however , the amplitudes between experiment and theory differed significantly for the third shape . here
the variability of the experimentally observed shapes plays a role , as well as the blurring of the contact zone due to resolution in experiments , which tends to increase the amplitude of the fitted sine function .
red blood cells aggregate due to the presence of plasma proteins or biomimicking polymers .
this is an example of adhesion between two very soft interfaces , different from the case of adhesion on a solid substrate , where we observe a buckling of the cell  cell interface .
the nature of this instability is revealed by our numerical simulations and explained by an analytical model .
the critical threshold of buckling and the final geometry of the aggregate depends strongly on the physical parameters , i.e. on the reduced volume of the cells . for large reduced volumes that correspond to more spherical cells , the critical interaction energy for buckling increases . in the case of doublets of two cells having identical reduced volumes close to that of a typical rbc , a sigmoid deformation of the interface takes place , while for two different reduced volumes male - female shapes can be obtained due to the asymmetry of properties that forces a different mode .
we also characterized long linear clusters experimentally and numerically and have found that only the three interfaces at the two ends are considerably deformed .
it would undoubtedly be interesting to check how these contact shapes obtained in quasi - static conditions evolve in pathological conditions , or under flow with the additional constraint of hydrodynamic stresses .
human blood withdrawal from healthy donors as well as blood preparation and manipulation were performed according to regulations and protocols that were approved by the ethic commission of the rztekammer des saarlandes ( reference no 24/12 ) .
we obtained informed consent from the donors after the nature and possible consequences of the studies were explained .
venous blood of three healthy donors ( one female , 26 years , blood group 0 negative ; one male , 43 years , blood group 0 negative ; one male , 27 years , blood group b positive ) was drawn into conventional edta tubes ( s - monovette ; sarstedt , nmbrecht , germany ) and all measurements were performed within 4 hours after the blood had been drawn . after washing the rbcs three times ( 704 g , 3 min ) with phosphate - buffered solution ( pbs , life technologies , waltham , ma , usa ) 0.1 @xmath12l cellmask ( life technologies , waltham , ma , usa ) was added for fluorescence labeling to 200 @xmath12l rbcs in 1ml pbs for 10min . after washing two times ( 704 g , 3min )
the supernatant was removed and the rbcs were resuspended in solutions of dextran ( dextran 70 from leuconostoc mesenteroides , sigma - aldrich , st louis , usa ) or ( washed ) fibrinogen ( fibrinogen from human plasma , sigma - aldrich , st louis , usa ) in pbs .
the fibrinogen as delivered contains approximately 60% fibrinogen , 25% sodium chloride and 15% of sodium citrate .
in order to remove the salt the solutions were ultrafiltrated with a 10kda filter ( vivaspin turbo 4 , satorius ag , goettingen , germany ) at 3500 g in a centrifuge ( hermle z326k , hermle labortechnik gmbh , wehingen , germany ) for 10 hours by adding 2ml of pbs every hour .
the remaining protein solution was freeze - dried and the proteins were resuspended in pbs .
fibrinogen was used in the measurements within 24 hours after freeze - drying and the osmolality of this solutions was verified to be in the order of pbs ( + /-
5% ) .
cylindrical wells of diameter 5 mm , height 1.6 mm ( @xmath12-slide 18 well - flat from ibidi , munich , germany ) were used as sample chambers .
rbcs sitting on the bottom slide of the well changed to echinocytes due to the glass effect . for the measurements in which we wanted to preserve the discoyte shape of sedimented cells , the wells were covered with 1 mg / ml bovine serum albumine ( bsa , sigma - aldrich , st louis , usa ) in pbs for 30 min at 37@xmath50c .
the bsa - pbs solution was removed and after a drying time of 10 min at 37@xmath50c the rbc suspensions were added .
all observations were performed with the confocal microscope described in the previous section at a temperature of 23 @xmath50c with an hematocrit of 0.2 % for the imaging of rbc doublets and 0.4 % for clusters larger than 6 rbcs .
images of doublets of rbc in their normal discoid shape could be recorded in free volume during the slow sedimentation process which prevented any contact with the glass surface of the cover slip .
this procedure did allow to study the aggregates without any other influence but we should also note that it was also not possible to investigate the interfacial shape of sedimented cells in discocytic state . in this configuration , they were horizontally oriented with their interface and the rbc membrane was very thin .
the optical 2-d sections showed a good signal only when the membrane was oriented along the axis of the point spread function of the laser focus .
however , images of sedimented rbcs were also taken on the non - treated cover slips , i.e. for the case of echinocytes which were often oriented with their interface perpendicular to the cover slip , most likely due to their rounder shape . in total ,
30 cells for discocytes and 15 for echinocytes were evaluated . for clusters that were larger than 6 cells
, measurements could be only be performed on sedimented cells . in this case
the cover slip was always treated with bsa to preserve the discocytic shape , and contact zones were always oriented perpendicular to the cover slip , thereby allowing a good characterization .
the raw images were contrast enhanced and the position of the interfacial membrane contours was manually determined by visual inspection . the point cloud that indicated the interface
was extracted with matlab ( the mathworks gmbh , ismaning , germany ) and stored into a one dimensional array to allow for a fitting procedure with origin ( originlab corporation , northampton , ma , usa ) .
the interface was fitted with a heuristically chosen sinus function which yielded reasonable agreement . for sigmoid
shapes the wave length of the sinus was left as a free parameter and for the male - female and parachute shape the wave length was limited to a range of @xmath51 around the double of the cell diameter .
this accounts only for the images shown in fig.[3dshapes ] . in the preparation step for the wells
celltak ( celltak , bc bioscience , san jose , usa ) was added as specified by the manufacturer in order to immobilize the cells on the bsa treated cover slip .
a stack of 50 images was taken at a distance of 0.1 @xmath12 m in z - direction , while the x - y area was 12 x 12 @xmath12 m . a surface rendering with imaris ( bitplane ag , zuerich , switzerland ) allowed to look through the cell membranes attached to the bottom onto the contact zones .
10    n.  mohandas and p.  g. gallagher .
red cell membrane : past , present , and future .
, 112(10):39393948 , 2008 .
u.  seifert and r.  lipowsky .
adhesion of vesicles .
, 42:47684771 , oct 1990 .
u.  s. schwarz and s.  a. safran .
physics of adherent cells .
, 85:13271381 , aug 2013 .
s.  komura , k.  tamura , and t.  kato .
buckling of spherical shells adhering onto a rigid substrate .
, 18:343 , 2005 .
s.  knoche and j.  kierfeld .
osmotic buckling of spherical capsules .
, 10:83588369 , 2014 .
s.  s. datta , s .- h .
kim , j.  paulose , a.  abbaspourrad , d.  r. nelson , and d.  a. weitz . delayed buckling and guided folding of inhomogeneous capsules .
, 109:134302 , 2012 .
d.  barths - biesel . motion and deformation of elastic capsules and vesicles in flow .
, 48:2552 , 2016 .
woo , , and a.  wallqvist .
spontaneous buckling of lipid bilayer and vesicle budding induced by antimicrobial peptide magainin 2 : a coarse - grained simulation study .
, 115:81228129 , 2011 .
v.  narsimhan , a.  p. spann , and e.  s.  g. shaqfeh .
pearling , wrinkling , and buckling of vesicles in elongational flows . , 777:126 , 2015 .
ionic strength dependence of localized contact formation between membranes : nonlinear theory and experiment .
, 77(2):817  828 , 1999 .
p.  ziherl and s.  svetina .
flat and sigmoidally curved contact zones in vesicle  vesicle adhesion . , 104(3):761765 , 2007 .
s.  svetina and p.  ziherl .
morphology of small aggregates of red blood cells . , 73:8491 , 2008 .
d.  tilley , w.  t. coakley , r.  k. gould , s.  e. payne , and l.  a. hewison .
real time observations of polylysine , dextran and polyethylene glycol induced mutual adhesion of erythrocytes held in suspension in an ultrasonic standing wave field . , 14(8):499507 , 1987 .
e.  w. merrill , e.  r. gilliland , t.  s. lee , and e.  w. salzman .
blood rheology : effect of fibrinogen deduced by addition . , 18(4):437446 , 1966 .
r.  e. wells , t.  h. gawronski , p.  j. cox , and r.  d. perera .
influence of fibrinogen on flow properties of erythrocyte suspensions .
, 207(5):10351040 , 1964 .
s.  chien , s.  usami , h.  m. taylor , j.  l. lundberg , and m.  i. gregersen .
effects of hematocrit and plasma proteins on human blood rheology at low shear rates .
, 21(1):8187 , 1966 .
mechanical deformability of biological membranes and the sphering of the erythrocyte . , 13(3):209  217 , 1973 .
r.  fhraeus .
the influence of the rouleau formation of the erythrocytes on the rheology of the blood .
, 161(2):151165 , 1958 .
e.  w. merrill , g.  c. cokelet , a.  britten , and r.  e. wells .
non - newtonian rheology of human blood - effect of fibrinogen deduced by `` subtraction '' .
, 13(1):4855 , 1963 .
a.  s. popel and p.  c. johnson .
microcirculation and hemorheology .
, 37:43 , 2005 .    j.  j. bishop , a.  s. popel , m.  intaglietta , and p.  c. johnson .
rheological effects of red blood cell aggregation in the venous network : a review of recent studies .
, 38(2 - 3):26374 , 2001 .
z.  marton , g.  kesmarky , j.  vekasi , a.  cser , r.  russai , b.  horvath , and k.  toth .
red blood cell aggregation measurements in whole blood and in fibrinogen solutions by different methods . , 24(2):7583 , 2001 .
b.  neu and h.  j. meiselman .
depletion - mediated red blood cell aggregation in polymer solutions .
, 83(5):24822490 , 2002 .
s.  asakura and f.  oosawa . on interaction between two bodies immersed in a solution of macromolecules . ,
( 22):12551256 , 1954 .
h.  baumler and e.  donath .
does dextran indeed significantly increase the surface - potential of human red - blood - cells .
, 120(2):113122 , 1987 .
e.  evans and d.  needham .
attraction between lipid bilayer membranes in concentrated solutions of nonadsorbing polymers : comparison of mean - field theory with measurements of adhesion energy .
, 21(6):18221831 , 1988 .
p.  steffen , c.  verdier , and c.  wagner .
quantification of depletion - induced adhesion of red blood cells .
, 110(1):018102 , 2013 .
h.  bumler , b.  neu , e.  donath , and h.  kiesewetter .
basic phenomena of red blood cell rouleaux formation . , 36(5 - 6):439442 , 1998 .
n.  maeda , m.  seike , s.  kume , t.  takaku , and t.  shiga .
fibrinogen - induced erythrocyte aggregation : erythrocyte - binding site in the fibrinogen molecule .
, 904(1):8191 , 1987 .
d.  e. brooks .
the effect of neutral polymers on the electrokinetic potential of cells and other charged particles : ii . a model for the effect of adsorbed polymer on the diffuse double layer . , 43(3):687699 , 1973 .
s.  chien and k.  m. jan .
red cell aggregation by macromolecules : roles of surface adsorption and electrostatic repulsion .
, 1(4 - 5):385409 , 1973 .
s.  chien .
shear dependence of effective cell volume as a determinant of blood viscosity .
, 168(3934):977979 , 1970 .
jan and s.  chien .
role of surface electric charge in red blood cell interactions .
, 61(5):638654 , 1973 .
s.  svetina , d.  kuzman , r.  e. waugh , p.  ziherl , and b.  zeks .
the cooperative role of membrane skeleton and bilayer in the mechanical behaviour of red blood cells .
, 62:107113 , 2004 .
m.  brust , o.  aouane , m.  thibaud , d.  flormann , c.  verdier , l.  kaestner , m.  w. laschke , h.  selmi , a.  benyoussef , t.  podgorski , g.  coupier , c.  misbah , and c.  wagner .
the plasma protein fibrinogen stabilizes clusters of red blood cells in microcapillary flows .
, 4:4348 , 2014 .
d.  flormann , e.  kuder , p.  lipp , c.  wagner , and l.  kaestner .
is there a role of c - reactive protein in red blood cell aggregation ?
, 37(4):474482 , 2015 .
o.  baskurt , m.  hardeman , m.  rampling , and h.  meiselman . .
ios press , 2007 .
badr kaoui , george biros , and chaouqi misbah
. why do red blood cells have asymmetric shapes even in a symmetric flow ?
, 103(18):188101 , 2009 .
dmitry  a. fedosov , matti peltomaki , and gerhard gompper .
deformation and dynamics of red blood cells in flow through cylindrical microchannels .
, 10:42584267 , 2014 .
o.  aouane , m.  thibaud , a.  benyoussef , c.  wagner , and c.  misbah .
vesicle dynamics in a confined poiseuille flow : from steady state to chaos . , 90:033011 , sep 2014 .
canham . the minimum energy of bending as a possible explanation of the biconcave shape of the human red blood cell .
, 26(1):6181 , 1970 .
w.  helfrich .
elastic properties of lipid bilayers : theory and possible experiments . , 28(11 - 12):693703 , 1973 .
g.  ghigliotti , t.  biben , and c.  misbah .
rheology of a dilute two - dimensional suspension of vesicles . ,
653:489518 , jun 2009 .
b.  neu and h.j .
depletion - mediated red blood cell aggregation in polymer solutions .
, 83(5):24822490 , 2002 .
yaling liu and wing  kam liu .
rheology of red blood cell aggregation by computer simulation . , 220(1):139154 , 2006 .
m.  thibaud and c.  misbah .
rheology of a vesicle suspension with finite concentration : a numerical study . , 88:062707 , dec 2013 .
the research leading to this results has received fundings from the german research foundation ( dfg , sfb 1027 ) , the centre national de la recherche scientifique ( cnrs ) and the german french university ( dfh / ufa ) .
t. p. and c. m. acknowledge support from cnes and labex tec21 .
d. f. performed experiments and data analysis , o. a. ran numerical simulations , l. k. supervised the experiments , c.r . programmed the analysis software , s.s . compared 3d analysis , c. m. supervised numerical simulations , t. p. designed the research , made the analytical model , supervised data analysis , c. w. designed the research , supervised experiments and numerical simulations .
a simple way to quantitatively characterize the clustering of the rbcs is based on the quantification of sedimented aggregates within a certain area on a microscope slide .
chien and jan @xcite defined a microscopic aggregation index ( mai ) as the ratio of the total number of cells to the total number of `` particles '' ( individual cells and aggregates ) , representing the average number of cells per particle . here , we investigate rbc aggregation by confocal microscopy and define the mai as the the ratio of the number of aggregated cells at the bottom of the well after sedimentation to the total number of cells , i.e. the fraction of aggregated cells .
all observations were performed at a temperature of 23 @xmath50c with a hematocrit of 10 % .
this concentration ensured an almost fully covered bottom surface after 5 min of sedimentation for most sample solutions , with an average coverage of @xmath52 cells / mm@xmath13 .
imaging was performed using a confocal microscope ( tsc sp5 ii ; leica microsystems , mannheim , germany ) with a laser of 633 nm as the excitation wavelength .
rbcs were classified as part of a cluster if the membranes of two cells were optically indistinguishable over a distance of at least 4 pixels , which corresponds to 1 @xmath12 m in our image with a size of 144 x 144 @xmath12m@xmath13 .
white - light microscopy was used as a control to exclude any possible influence of cellmask on aggregation ( data not shown ) .
fibrinogen and dextran solutions were compared to pbs as negative control . due to the increase of viscosity with increase in dextran or fibrinogen concentration , the sedimentation slows down , which leads to a lower total number of cells per image at higher viscosities in comparison to rbcs in pbs .      in dextran solutions , the number of _ rouleaux _ and the mai exhibit a bell shaped curve with a maximum around a concentration of 50 mg / ml and almost no _ rouleaux _ at 120 mg / ml .
[ microindex](b ) ) . for fibrinogen , a monotonous increase of the mai is seen up to a concentration of 12mg / ml , which corresponds to the physiological range .
both results are reminiscent of the dependencies of interaction energy , in - flow clustering and rheology on dextran concentration obtained in previous studies @xcite .
another qualitative observation is that the dextran and fibrinogen induced _ rouleaux _ look slightly different with more interconnections in the fibrinogen case .
all measurements of contact zone shapes of rbc clusters were performed at 0.2% v / v hematrocrit .
the procedure for extracting interfacial shapes for further analysis is exemplified in fig .
[ figcontour ] .
the contact zones of rbc clusters are first isolated .
the images are then thresholded and skeletonized to extract the interface shape on which a sine function can be fitted to measure the amplitude of deformation .
when the doublets were allowed to sediment on the bsa treated cover slip , the cells remained in discocyte state for several minutes but they were mostly oriented with their interface parallel to the coverslip which did prevent a characterization of the thin interfacial membrane that was then oriented perpendicular to the z - axis of the point spread function .
this was different for the non - treated cover slips : due to the so called `` glass effect '' , the rbcs transformed into the echinocyte state and we observed only male - female and sigmoid shapes ( fig . [ figechinocytes](a ) ) .
a quantitative analysis of the distribution of the different interfacial shapes for the different concentrations of dextran is shown in fig . [ figechinocytes](b ) .
for all different cases we observed an increase in the number of sigmoid interfaces in the intermediate range of polymer concentrations , i.e. at larger interaction strengths .
the sigmoid shape allows for the largest contact zone , favorized by large interaction strength . at low and high concentrations of dextran we observe only few sigmoid shapes and more parachute than male - female shapes for the cells in discocytic state . with fibrinogen , mostly parachute shapes were observed .
this is most likely related to the low ( but physiological ) concentrations of the protein . for cells in echinocyte state
no parachute shape was observes which is to be expected in view of the spherical shape of a single echinocyte .      based on several successful studies ( especially equilibrium shapes and shape dynamics under external flow ) that have shown that 2d and 3d models capture the same essential features @xcite , have restricted ourselves here to a 2d model . in 2d
the notion of in - plane shear elasticity looses its meaning and the vesicles and inextensible capsule model ( often evoked to model rbc ) are equivalent .
vesicles consist of closed bilayer membranes enclosing an inner fluid and suspended in an outer fluid . using the typical reduced area ( equivalent of the surface - to - volume ratio in three dimensions ) of the rbc
, the minimum energy shape of an unstressed vesicle is the same as for an rbc despite the lack of a cytoskeleton .
indeed , the stress free shape is dictated by the curvature energy and the rigidity of the cell .
to describe the mechanical behavior of the vesicle , we employ ( for each cell ) the two dimensional helfrich free energy@xcite ( per unit length in the 3rd dimension ) : @xmath53 where @xmath54 is the mean curvature , @xmath55 is the arclength element along the membrane contour , @xmath56 is the membrane bending rigidity modulus , and @xmath57 is the spontaneous curvature that can be neglected in case of symmetrical membranes ; actually in 2d the spontaneous curvature plays no role . in order to fulfill the condition of membrane inextensibility we supplemented the curvature energy with a term @xmath58 , where @xmath59 is a lagrange multiplier that depends on the position along the membrane .
the fluid inside and outside the membrane is described by the stokes equations ( the inertial effect is small and is neglected ) with boundary conditions at the membrane and far away from the membrane .
far away from the membrane the flow vanishes ( in the absence of any imposed external flow ) . at the membrane
we require velocity continuity , and force balance .
in addition , the velocity field along the membrane has to be divergence - free in order to enforce membrane incompressibility .
this condition represents an implicit equation that allows one to determine the still unknown lagrange multiplier @xmath59 ; the details of how this is achieved in practice can be found in @xcite .
the hydrodynamics stress jump across is counterbalanced by the membrane force .
this force is composed of the bending force obtained as a functional derivative of the helfrich free energy supplied with a tension - like force to fulfill the condition of inextensibility of the membrane .
@xmath60\mathbf{n }   - c \zeta \mathbf{n } + \frac{\partial{\zeta}}{\partial{l}}\mathbf{t}\]]where @xmath61 and @xmath62 are the normal and tangent unit vectors . @xmath63
designates the vector position of a given point on a given cell ; the curvature as well as the lagrange multiplier are defined for each cell .
in addition to this force we had to consider the force due to cell - cell interaction ( due to the presence of proteins or dextran ) . for this purpose
we refer to the work of neu and meiselman @xcite .
neu and meiselman have presented a model that takes into account the strong electrostatic repulsion due to the negative charge of rbc membranes , and the weak osmotic forces of attraction due to the depletion effect induced by macromolecules surrounding rbcs @xcite .
this model contains several parameters depending on the physiochemical properties and the concentration of the polymers .
often , the interaction energy between cells is described by a morse potential ( or alternatively a mie type potential ) , having an exponential form , which has proven to adequately describe experimental data @xcite . as is the case in many physical problems , a morse potential can be well approximated by a lennard - jones potential , and we adopted this form in the present study for practical reasons .
the interaction potential per unit ( length of the 1st membrane @xmath64 length of the 2nd membrane @xmath64 length in the 3rd ( invariant ) dimension ) is @xmath65\ ] ] where @xmath66 is a microscopic parameter with dimension [ energy / length@xmath67 .            in this formulation
, @xmath66 is a microscopic parameter of interaction between surface elements .
given the range of the interaction potential ( typically a few times @xmath78 ) , an integration of the potential over planar parallel surfaces of extension much larger than this range must be performed in order to derive the macroscopic interaction energy per unit area @xmath11 .
the zero force position is therefore @xmath81 , which is smaller than @xmath82 .
the interaction energy per unit surface between infinite plates ( _ i.e. _ whose extension is much larger than @xmath83 ) is then :        recently the rbc - rbc adhesion energy has been quantified using the atomic force microscopy - based single cell force spectroscopy  @xcite . from these data , in principle , we can determine the adhesion energy @xmath11 and equilibrium distance @xmath83 from which the two parameters entering the lennard - jones potential can be determined , namely the minimal energy @xmath66 and the range of the potential @xmath78 .
however , here we chose to adapt the results from the numerical simulations , namely the amplitude of the interfacial deformation , directly to the experimental data and use this to deduce the interaction energy . in this procedure , due to numerical constraints ,
the zero force length ( i.e. the length corresponding to minimal potential ) was fixed in our simulations to 489 nm .
this was probably significantly larger than the experimental value but still in the order of the optical resolution limit .
in addition , note that as long as @xmath78 is small compared to the lateral extension of the membrane , the product @xmath85 is the most quantity here .
the stokes equations were solved using the boundary integral technique as presented in  @xcite    @xmath86 dl(\mathbf{x } ) } \\+ \frac{(1-\lambda)}{2 \pi ( 1 + \lambda)}\int\limits_{\sum_{i}{\partial \omega_{i}}}\mathbf{u}(\mathbf{x } ) \cdot \mathbf{t}(\mathbf{x},\mathbf{x_0 } ) \cdot \mathbf{n}(\mathbf{x})dl(\mathbf{x } ) \ ] ]    where @xmath87 is the velocity of a vector position @xmath88 belonging to the membrane @xmath89 , @xmath90 is the viscosity contrast between the inner and outer fluids with the respective viscosities @xmath91 and @xmath92 , and @xmath93 and @xmath94 are the green s functions associated with the velocity and the stress fields respectively .    to match the average values of a human rbc , the reduced area ( @xmath14 ) of the vesicle typically lies in the range ( @xmath95 ) , the typical radius is @xmath96 , the membrane stiffness is @xmath97 .
here we did not investigate any dynamic effects and we therefore fixed the viscosity ratio simply to unity , even if typically a factor of 5 was assumed .
we considered doublets and sets of @xmath98 cells , separated by an initial surface - to - surface interdistance of 600 nm .
we provide here the details of the analytical model showing buckling of contact zones in symmetric rbc doublets .
we suppose that the shape of an rbc in a symmetric doublet ( in 2d ) is approximated as shown in fig .
[ fig - model ] , with two semi - circular caps of radius @xmath17 , a roughly straight portion where curvature is negligible ( length @xmath18 ) and a contact zone with a sinusoidal shape of curvilinear length @xmath19 and amplitude @xmath20 described by the equation :                      in the following we define a dimensionless energy @xmath28 , a rescaled interaction energy @xmath106 , and rescale all lengths by @xmath8 : @xmath30 and so on , and use the definition of the reduced area @xmath7 .    the constant perimeter and area constraints can then be written as the following ( dropping the * for simplicity ) : @xmath107 and can be used to eliminate @xmath17 and @xmath18 as in equations [ rexpr ] and [ l1expr ] .
if we solve @xmath111 to find the equilibrium configurations we get the two solutions whose expressions are given by [ aeq ] , which defines the bifurcation threshold @xmath34 corresponding to @xmath38 ( eq . [ betac ] )