The Connection between Biaxial Orientation and Shear Thinning for Quasi-Ideal Rods

The complete orientational ordering tensor of quasi-ideal colloidal rods is obtained as a function of shear rate by performing rheo-SANS (rheology with small angle neutron scattering) measurements on isotropic fd-virus suspensions in the two relevant scattering planes, the flow-gradient (1-2) and the flow-vorticity (1-3) plane. Microscopic ordering can be identified as the origin of the observed shear thinning. A qualitative description of the rheological response by Smoluchowski, as well as Doi–Edwards–Kuzuu theory is possible, as we obtain a master curve for different concentrations, scaling the shear rate with the apparent collective rotational diffusion coefficient. However, the observation suggests that the interdependence of ordering and shear thinning at small shear rates is stronger than predicted. The extracted zero-shear viscosity matches the concentration dependence of the self-diffusion of rods in semi-dilute solutions, while the director tilts close towards the flow direction already at very low shear rates. In contrast, we observe a smaller dependence on the shear rate in the overall ordering at high shear rates, as well as an ever-increasing biaxiality.


Introduction
The connection between the macroscopic rheological properties of semi-flexible or stiff polymer dispersions and the underlying microstructural changes is of high fundamental importance and also plays a huge role in many industrial processes. As a starting point for understanding semi-flexible systems, one can investigate a quasi-ideal stiff polymer suspension and define its rheological behavior, as well as its structure in flow. This benchmark problem could help to disentangle the manifold contributions to the complex flow-response of semi-flexible systems and make further progress in their theoretical description. Shear thinning is one of the main rheological signatures of a complex structural response of a material to shear flow. Generally, this type of non-linear response is probed by steady shear experiments where an increasing shear rate is applied and a stress is measured or vice versa. For a fundamental understanding of shear thinning in terms of the microstructural changes of sheared material, a technique is needed that probes the sample on the relevant length-scales. Moreover, model systems are needed that display very reproducible flow behavior and which can ideally also be modeled with the appropriate microscopic theory. The goal of this paper is exactly this, namely to study the connection between the microstructural and mechanical response of the model methodology, including the explanation of how we obtain all relevant ordering parameters from the experiments in both scattering geometries. We then present the results and discuss the agreements and discrepancies that we find with theory.

Smoluchowski Theory
The Smoluchowski theory is very suitable for the description of nonequilibrium phase behavior as it incorporates in an elegant way the external field, as well as interaction potentials, calculating the relevant equation of motion for a probability density function. In the case of colloidal rods, we are interested in their ODF Ψ(u, t), which is a function of the orientation of a rod u and time t. The N-particle Smoluchowski equation leads to [15]: where D 0 r = 3k B T ln(L/d)/πη s L 3 is the rotational diffusion coefficient in an infinitely-dilute suspension depending on the solvent shear viscosity η s , as well as the length L and thickness d of the rod, andγ is the shear rate. The rotational operator is defined as R = u × ∂/∂u, β = 1/k B T; the shear-field induced torque on the rods is: T(u, t) = −ρ dR du Ψ(u, t)g(R, u, u , t)RU(R, u, u , t) (2) where g(R, u, u , t) is the pair correlation function of the rod under inspection and an additional rod with orientation u with the center-to-center distance R = r − r , U is the pair interaction potential of the two rods and ρ = N/V is the number density. With the tensor: and the shear rate, we can write down the equation of motion for the orientational ordering tensor S from the Smoluchowski equation: The orientational ordering tensor depends on the ODF in the following manner: In Equation (4), the volume fraction of rods in the suspension is given as ϕ = (π/4)d 2 Lρ, and we define the following tensor in addition to the unit tensor I: To solve Equation (4), we use the following closure relation [15] for both of the double contracted forms uuuu : A of the fourth order orientational tensor: where A stands for either S or E and A = (A + A T )/2. It is important to notice that in deriving the equation of motion for S (Equation (4)) from the N-particle Smoluchowski equation for the ODF (Equation (1)), the shear rate-independent expression for the pair correlation function: is used, following Onsager [22] and Doi and Edwards [19]. Here, r and r denote the positions of two neighboring rods, while u and u denote the orientations of their long axis. Assuming a hard-core interaction, this gives the excluded volume: The fact that there is no shear rate dependence in g(R, u, u , t) is one of the main assumptions of the theory, also directly leading to the given form of Equation (9).
The connection of this microscopic description, using the orientational ordering tensor, to the mechanical macroscopic measurables can be made [15] via the deviatoric stress tensor Σ, given as: The shear viscosity can be expressed in terms of the relevant component of the stress tensor using the following constitutive relation, which can be derived by plugging an expansion of the ordering tensor to the third power in shear rate into Equation (10): where C = 8(L/d) 2 /45 ln(L/d) and the concentration-dependent collective rotational diffusion coefficient D coll r is used because the ODF also depends on the distance to the isotropic-nematic transition. This links the shear viscosity with the order parameter tensor and its corresponding eigenvalues.
Finally, also the zero-shear viscosity η 0 as a function of volume fraction ϕ of rods can be calculated from this expression, leading to what is called a mixing rule for this quantity: Equation (12) can be compared, e.g., to the theory by Berry and Russel [31], where the polynomial for the zero-shear viscosity is of second order in ϕ: This mixing rule clearly shows a much stronger concentration dependence of the zero-shear viscosity, as was already pointed out by Dhont and Briels [15].

Doi-Edwards-Kuzuu (DEK) Theory
The DEK theory constitutes a basically similar, but nonetheless very differently interpreted alternative to the theory outlined above [15].
On the one hand, the rotational diffusion coefficient of the dilute system D 0 r used in Smoluchowski theory is replaced by a purely phenomenological function D r . What is taken into account in this function can be made quantitative by the tube model [19]. Thereby, time-dependent tube disengagement of the rods leads to changed rotational relaxation times if the system is in the semi-dilute state. A comparison of the mean square displacement of rods in the dilute state to that in the semi-dilute state directly gives a scaling law for the effective rotational diffusion coefficient. This can be written as: where C is a coefficient with magnitude 10 3 to 10 4 [19]. According to Doi and Edwards [19], one has to take into account this apparent concentration dependence of the rotational diffusion coefficient in the Smoluchowski Equation (1). However, this description is only valid if the surrounding rods under consideration constitute a tube without possessing a preferred direction themselves. In the case of an anisotropic distribution of the surrounding, an effective mean change of the tube has to be considered if the picture of the tube as used by Doi and Edwards is believed to be a valid means of description for the given problem. Then, the function for the rotational diffusion coefficient becomes orientation dependent [19]: Furthermore, if the external perturbation is large, one would have to take tube dilation into account. A description of such a situation is only feasible up to a certain degree of accuracy [19]. Although the DEK approach provides one of the possible non-analytic expressions for such a case, we do not include the corresponding calculation in our treatment, since we want to follow the procedure given by Doi and Edwards [19].
On the other hand, this description is accompanied by neglecting the pair-interaction potential U in Equation (1) entirely. In doing so, the Smoluchowski equation with the alternative description of rotational diffusion, but in absence of interaction, is solvable analytically. Thus, one can give a direct equation for the orientational order-dependent shear viscosity: where the rotational diffusion coefficient is given by Equation (14), the orientational order is represented by the order parameter P 2 , which depends on the highest eigenvalue of S, called λ 1 , as P 2 = 1 2 (3λ 1 − 1), and η 0 is the zero-shear viscosity. It needs to be mentioned here that Equation (16) has been derived by Doi and Edwards [19] for the case of concentrated solutions. It can, however, be used for our purpose, as will be shown in Section 4.
Using the measurements at hand, we are able to review both theoretical approaches from a new perspective.

Scaling
In comparing experiments with theory, one needs to choose the right scaling for shear rate, as well as viscosity. As we indicated in the previous section, the viscosity is enslaved by the ordering in the system. In earlier experiments on worm-like systems [32], indeed, a scaling law for the viscosity was obtained when plotting the viscosity, scaled by the zero-shear viscosity, against the orientational order parameter P 2 , the highest eigenvalue of a projection of S into the measurement plane.
The appropriate scaling of the shear rate can be done by relating it to the rotational diffusion of the rods, D r , such that: where Pe is the so-called Peclet number. The question is if one should consider the collective rotational diffusion D coll r or the self-diffusion D self r , where we have to distinguish between the self-diffusion at infinite dilution D 0 r and the self-diffusion at the volume fraction of interest, as given in Equation (14). The right choice for this scaling will enable us to identify the interdependence of macroscopic flow behavior and microscopic orientational order, which is the focus of our experimental investigation. It should be mentioned that since the relaxation time of the rods can be given as τ r = 1/2D r [33], we will expect a value of τ rγ = 1 at Pe = 2.
The collective rotational diffusion coefficient D coll r goes to zero when approaching the isotropic-nematic spinodal from a dilute, isotropic dispersion by increasing the scaled volume fraction L d ϕ [15]. This number should be compared to its theoretical value at the spinodal point, which can be derived by studying how a perturbation from an isotropic state δS evolves over time, using Equation (4) with the indicated closure relation and setting the shear rate to zero. This procedure leads to: where D coll ). Therefore, the isotropic-nematic spinodal point is located at a scaled volume fraction of L d ϕ I N = 5. It should be mentioned that, depending on which assumptions are made on the interaction potential U between the rods, the values for this number range between L d ϕ I N = 4 and L d ϕ I N = 5. Furthermore, if the effective rod thickness d eff is greater than the actual thickness d, due to the surface charge, the value has to be adapted accordingly. We will, therefore, write down the relevant equations using the effective thickness.
It is to be anticipated that the scaling with D coll r will not suffice, as the self-diffusion is strongly reduced at elevated concentrations. This has been nicely shown in event-driven Brownian dynamics simulations by Tao et al. [13]. From these simulations an alternative expression for D coll r was found for rods of sufficient length: D coll where ν = (1/ √ 2)(L/d eff ) 1/4 and A is a constant, which normalizes the function at ϕ = 0 [13]. This scaling takes the decrease in the rotational self-diffusion into account to some extent. We will use it for the definition of an effective Peclet number Pe eff =γ/D coll r by fitting the correct value of ϕ I N for our system. This fit will be done such that a master curve of the order parameter, as well as scaled viscosity versus Pe eff is obtained, assuming that the viscosity is fully enslaved by the ordering and the ordering is fully enslaved by the collective rotational diffusion. In our theoretical predictions, we will use a corrected Peclet number Pe corr =γ/D coll r , where the theoretical value for ϕ I N will be used for defining a suitable D coll r .

Measurements and Materials
The 1-3 plane rheo-SANS and stress data were acquired simultaneously by mounting an Anton Paar MCR 501 rheometer (Anton Paar GmbH, Graz, Austria), used in the strain-controlled mode with a quartz Couette geometry, in the SANS-1 neutron beam at the SINQ spallation source at the Paul Scherrer Institute (PSI) in Villigen, Switzerland. 1-2 plane scattering experiments were carried out by using a shear cell mounted to the d22 large dynamic range small-angle diffractometer at the Institut Laue-Langevin (ILL) in Grenoble, France [29]. The scattering geometry is shown in Figure 1b. All measurements were performed at 28.5 • C.
In order to find a compromise between intensity and the probed length scale of the rod, we scanned the scalar order parameter over a range of q-values under stationary flow conditions using a detector distance of 6 m and wavelength of the thermal neutrons of 1.3 ± 0.1 nm. The Porod plot for a rod concentration of 11 mg/mL in absence of shear is shown in Figure 2, where the rigid-rod slope I(q) ∼ q −1 is seen at quite low q-values, and it crosses over to a Gaussian chain slope I(q) ∼ q −2 at higher q. We choose the second range for our analysis and evaluate a q-range between 0.032 and 0.046 Å as the intensity is higher in this range. However, an evaluation at smaller scattering vectors does not influence the outcome in the azimuthal averaging profile necessary for the computation of the projected order parameter P 2 , as can be seen from the overlapping of curves in the inset of Figure 2. The measurement data were reduced using the GRASP freeware provided by the ILL (Grenoble, France) and subsequent analysis in Wolfram Mathematica (Wolfram, Oxford, UK) .  While the Anton Paar MCR 501 is inherently a stress-controlled rheometer operating in strain-controlled mode, the rheological measurements presented here were replicated within experimental uncertainties using a strain-controlled ARES G2 rheometer (TA Instruments, New Castle, PA, USA).
The samples for the two different measurements were prepared (see [34]) by using a slightly different Tris buffer. The corresponding concentrations are 90 and 110 mM Tris with suspended filamentous bacteriophage fd in the following concentrations: c ∈ {4.1, 11, 12.43, 15, 20} mg/mL. Respective effective diameters of fd-virus can be considered as d eff (90 mM) ≈ 8.64 nm and d eff (110 mM) ≈ 12 nm [3].

Obtaining the Full Orientation Tensor and the Biaxiality
Intensity profiles, as shown in the inset of Figure 2, are directly proportional to the normalized ODFs f (θ) and f (ψ) [35], which are projections of the overall ODF f (α) into the 1-2 or 1-3 plane, respectively. Thus, the director, which describes the averaged orientation α of the rods with the flow, makes an angle ψ with the flow direction in the 1-3 plane and an angle θ with the flow direction in the 1-2 plane, as can be seen in Figure 1a.
In order to parametrize the orientational ordering by single scalars [36], we fit the intensity profiles with a Maier-Saupe type of ODF: Here, I 0 is the amplitude,C describes the width of the intensity profile and is the second order Legendre polynomial, allowing for the director to point along α max . It has been shown earlier that this function fits the orientational distribution well, from low to high ordering [37][38][39]. This function can be further reduced to a single scalar that quantifies ordering, given by the average orientational order parameter P 2 (α) , which can be calculated using the fitted parameters: (21) In this way, the orientational ordering in this plane is parametrized by P 2 (α) and α max . The procedure is identical for both parametrizations θ and ψ.
To relate the order parameters to the orientational ordering tensor S, we first note that the apparent largest eigenvalues λ 1 (α) of the two plain tensor projections are related to P 2 (α) as P 2 (α) = 1 2 (3λ 1 (α) − 1). In addition to the aforementioned orientation tensor S, we need the traceless orientation tensor Q = (3S − I)/2 in order to describe the relative rod positions in the reference frame of the measurement. The two tensorial quantities are related as follows: where: is the rotation tensor of the laboratory reference frame andQ = 1 3 (2Q + I) is the full orientational ordering tensor in this reference frame. We know from the measurement in the 1-3 direction that the maximum angle ψ max ≡ 0, irrespective of the applied shear rate. Therefore, the rotational tensor only depends on the angle θ. From this observation and the definition of the orientation tensor with the expected symmetry conditions, one finds the following equations for the components of Q: where λ i (α) is defined as the i-th eigenvalue of S measured in the plane parametrized by the angle α, and the biaxiality parameter T can be found as: Since all of the eigenvalues λ i (α) are measured, the traceless tensor Q can be calculated from knowledge of the angle θ, which defines the rotation tensor. This, of course, also gives the full orientational ordering tensorsQ and S. The measured order parameter tensor is evaluated at every shear rate, which corresponds to the effective Peclet number of the system in that instant.

Results
We first study the interdependence between the shear viscosity and the projected order parameter in the 1-3 plane of flow for three different rod concentrations in the dilute to semi-dilute regime. Since the single-particle contribution to the viscosity is negligibly small, we assume that there exists a master curve for the overall comparison between shear thinning and orientational ordering, independent of the rod concentration. For obtaining the master plot, one would ideally scale the viscosity by the zero-shear viscosity η 0 , but for this system of long rods and for this rheometer, it is not feasible to obtain the zero-shear viscosity, as shear thinning sets in immediately. Therefore, we scaled the viscosity with an adjustable η 0 , such that the data collapse on a single curve in a log-lin plot; see Figure 3. For each concentration, we thus find a zero-shear viscosity as the correction value of the shear viscosity in Figure 3. In Figure 4 we compare the thus obtained η 0 , scaled by the solvent shear viscosity, with Smoluchowski, as well as DEK theory and a comparable theory by Berry and Russel [31]. It is seen that neither the Smoluchowski theory nor the theory by Berry and Russel fit the data. However, DEK theory approximates the data reasonably well.
In Figure 3, we also plot the theoretical predictions as full and dashed lines. The discrepancy between theory and experiment in this representation is evident, and neither DEK nor Smoluchowski theory values predict the actual experimental outcome. Fitting the order parameter-dependent corrected viscosity curves with a purely empirical exponential form, suggested by Förster et al. [32] (see the dotted thick line in Figure 3), we find that: describes the data well when a = 5.0. Next, we focus on the relation between shear viscosity and shear rate. A master plot can be obtained for η corr = η/η 0 versus the shear rate, when defining an effective Peclet number Pe eff . The scaling of the shear rate works well when using Pe eff =γ/D coll r , where D coll r is given by Equation (19) and choosing L d ϕ I N = 4.2; see Figure 5. The focus of this scaling lies on the higher shear rate regime, due to the low torques at low shear rates. The same scaling of the shear rate is used when plotting the projected order parameter P 2 (ψ) (in the 1-3 plane) as a function of the imposed shear rate. Again, a collapse of data is achieved by the same Peclet numbers as above; see Figure 5. Note that for the theoretical data to fall on top of each other, we need to use the theoretical L d ϕ I N = 5, giving a corrected Peclet number Pe corr (see Section 2.3), while the scaling law by Tao [13], Equation (19), remains valid in its given form. As the master curve with the right scaling of the shear rate works very well for the orientational ordering, as well as for the viscosity, using L d ϕ I N = 4.2, we now proceed to compare the behavior of rods in the two measurement directions. We apply the same Peclet number scaling to the projected order parameter in the 1-2 plane, as well as in the 1-3 plane; see Figure 6. Before doing so, we need to correct for the disparate effective thicknesses in the two measurements, see Sections 2.3 and 3.1, as well as Equation (19). Again, a data collapse for the curves of different rod concentrations in both measurement directions can be achieved. It is seen that the observed orientational order parameter in the 1-2 plane is higher than the one in the 1-3 plane from the very beginning of shear-induced order, corresponding to the low Peclet number regime. The discrepancy grows non-linearly with Pe eff and shows the largest differences in the high Peclet number regime.
The 1-2 plane measurements provide additional information on the achievable order in the system, as they also give a value for the angle θ the director makes with the flow direction, as plotted in Figure 7b. With the appropriate identified scaling of the shear rate, the orientational order parameters P 2 (θ) and P 2 (ψ) and the orientation of the director θ (given that ψ max ≡ 0), we have the necessary information at hand to identify all of the parameters that describe the orientational ordering, namely θ, the highest eigenvalue λ 1 and the biaxiality T, using Equations (22) and (25) in Section 3.2. These results can be compared to the theoretical prediction from Smoluchowski theory; see Figure 7a,b. It is important to notice that the Peclet numbers of two measurements for different concentrations are not equal, due to the used scaling. Therefore, the full tensorial order can only be found by fitting the relevant data |θ|(Pe eff ), shown in Figure 7b, with an appropriate function and, thus, smoothing out the apparent discrepancies without having the need for averaging. The used fitting function is a Cole-Cole type function [40], providing a sigmoid shape in a log-log representation.
As expected, the angle θ is indirectly proportional to the order parameter and decreases with increasing shear rate. The theoretical prediction starts from a higher mean angle of approximately |θ| = 45 • and has a much steeper decay than the experimental finding. Theory also ultimately predicts a total decay of the angle at very high shear rates. However, measurements cannot be conducted up to values where the mean angle of the rods with the flow direction is effectively zero. Similarly, λ 1 does not reach the high ordering as predicted in theory for high Peclet numbers, while at low Peclet numbers, λ 1 starts to increase at much lower Peclet numbers than predicted by theory; see Figure 7a. The biaxiality parameter T of the system as a function of the corrected Peclet number is plotted in Figure 7c. In experiment, biaxiality is observed as soon as ordering occurs, and it grows together with λ 1 . The theoretical prediction, on the other hand, displays an increasing biaxiality only up to a certain point, corresponding roughly to d 2 λ 1 /dPe 2 corr = 0 and a decay in biaxiality after this point until a certain constant value is reached, which remains up to very high imposed shear rates.

Shear Thinning
As can be seen from Figure 5, there exists a universal scaling for the Peclet number-dependent order parameter if the rotational diffusion that is used to scale the shear rate has the proper concentration dependence. For the shear viscosity versus shear rate data, the same scaling can be applied. In order to have a data collapse here, the shear viscosity needs to be shifted by using the zero-shear viscosity. This procedure leads to a master curve for shear thinning, indicating two important aspects of the rheological behavior of dilute and semidilute suspensions of long, thin, mono-disperse rods: First, we found in Figure 4 that the dependence of the zero-shear viscosity on the relative volume fraction of rods, however few points we have, can be best described by a non-linear polynomial mixing rule of the form: where C = 5800 and r = L/d is the aspect-ratio of the viruses. This mixing rule has been derived by Doi and Edwards [19] using Equation (1), setting U = 0 and inserting an empirical rotational diffusion coefficient, which depends itself on the relative volume fraction of rods: D r = CD 0 r (ρL 3 ) −2 , as briefly described in the theoretical Section 2.2. The given relation stands in marked difference to mixing rules derived via Smoluchowski theory, laid out above, which yields a first order polynomial for η 0 (ϕ)/η s , or the theory by Berry and Russel [31], which yields a second order polynomial. In our experimental investigation, the concentration dependence of the apparent rotational diffusion coefficient is slightly larger than the predicted power of −2 dependence by Doi and Edwards [19]. This result could be compared to numerical simulations [41], where the exponent was found to be slightly smaller than −2, as well as experiments [5], where the dependency was seen to be even less than in the mentioned simulations. One can, therefore, proceed with the criticism of the existing theoretical approaches given here in two different directions, always assuming that the given closure relation in Equation (7) assures a good accuracy, as has been shown by Dhont and Briels [15]. It should, nonetheless, be mentioned here that there remains a possibility that by coming close to the concentration of the isotropic-nematic transition, the closure relation Equation (7) becomes less and less accurate, as it has been shown by Kröger et al. [42] that the used closure does not give as good results as other closures in the nematic phase. Since the discrepancies between theory and experiment are much bigger than the differences in closure accuracy, we will not discuss the associated issues, but instead focus on the fundamental differences between the shown theoretical approaches. Either the problem is seen in the definition of a correct rotational diffusion coefficient, which can be thought of as depending on both, the shear rate and the concentration, or we attribute the huge discrepancies to the definition of a correct pair correlation function, which should be a functional of the imposed field in non-equilibrium. In order to give this reasoning a sound grounding, experiments would be needed that probe the dynamics of rods under shear in order to test the variability of the rotational diffusion coefficient under the important influence parameters. This, however, is not a trivial task.
Since the given difference between Smoluchowski theory as we use it here and the approach by Doi and Edwards is solely attributable to the energy, which is dissipated by ordering of the rods, we single out the used pair correlation function as the weakness of the given theoretical approach. Our conclusion is based on the premise that the rotational diffusion coefficient in the small shear rate regime should not depend on the shear rate, as is assumed in the theory of Doi and Edwards by purely phenomenological reasoning. Instead, given the rotational diffusion coefficient of the infinite dilute regime D 0 r as a state-independent coefficient in Equation (1), a dynamic relation for the pair correlation function is needed, which leads to a higher polynomial order in the mixing rule.
Second, the corrected shear viscosity η corr strongly decreases while the orientational ordering strongly increases with increasing shear rate; see Figure 5. Given the fact that the zero-shear viscosity is much higher than predicted from Smoluchowski theory, this is not surprising. The effect is especially pronounced at very low shear rates, where ordering, as well as the orientation of director θ immediately jumps from 1/3 to λ 1 ≈ 0.4, while θ jumps from 45 • to θ ≈ 25 • ; see Figure 7. This means that shear thinning, as a rheological signature of our system of rigid mono-disperse rods, is most significantly caused by immediate strong ordering of the rods towards the flow direction.
In addition, we observe that the angle θ, which the rod director confined to the 1-2 plane makes with the flow direction, displays a power-law dependence on the shear rate, as can be seen from Figure 7b, which also has been observed in multi-particle collision dynamics simulations [43]. We find that the power-law is independent of the concentration of our suspension if we scale the Peclet number in accordance with Equation (19). The theoretical prediction of a fully-decayed angle at high shear differs from the measurement, which shows that the angle does not reach this plateau also at the highest measurable shear rates, as can be seen from Figure 7b.
Again, the origin of the discrepancy between theory and experiment in Figure 3 should be looked for in the shear-dependent dynamic pair correlation function. Clearly, the non-equilibrium behavior we observe here cannot be accounted for by Equation (8) as has been already pointed out in an earlier paper [15]. Since the expression that Doi and Edwards derive for their empirical rotational diffusion coefficient in this non-equilibrium case has not been used in their outlined analytic approach [19] (see Section 2.2), their prediction fails along the same lines as the Smoluchowski theory.
For high shear rates, so at the moment the tube is fully dilated and there are very little interactions between the rods, the correct scaling of the shear rate should in principle not be done by using D coll r , but rather by using D 0 r . This is seen from evaluating the largest eigenvalue of the orientational ordering tensor as a function of the bare Peclet number, using the rotational diffusion coefficient of the infinite dilute case, shown in the inset of Figure 7a. When doing so, the theoretical values in the large shear rate regime (II in the inset in Figure 7a) are closer to the experiments than what is seen in Figure 7a. This might also explain why in the experiments, λ 1 and θ never reach a plateau and why a quantitative match of theory with experiments of nematic rod-like colloids of the same kind in shear has been found [14]. The reason is that dilation does not play an important role there. Furthermore, this tells us that there exist two distinct regimes in the viscosity versus shear rate curve: a low shear rate regime where the rotational diffusion coefficient strongly depends on the volume fraction, as discussed before, and a shear rate regime that does not show this dependence.

Biaxiality
By using the rheo-SANS measurements in two scattering geometries, we are able to obtain the full ordering tensor Q from which, after diagonalization and rotation, the biaxiality parameter T can be calculated. In Figure 7c, we see a growing biaxiality with growing order, which goes to a plateau value at a Peclet number of approximately Pe e f f ≈ 600, always following the behavior of the largest eigenvalue of S. It comes as no surprise that the dispersions of rods show biaxiality, since order in the gradient direction is forced by the velocity gradient in shear flow, while the exchange of momentum in the vorticity direction is small. The same observation has already been made, however, with rods of a smaller aspect-ratio, in multi-particle collision dynamics simulations [43]. Also for living polymers, which posses characteristics of rods, as well as soft polymers, biaxiality could be seen to occur, although the orientation of the ordering plane there depends on the concentration and is due to hairpin formation, which does not play a role for the very stiff fd-virus [25,27]. That there is order in the vorticity direction at all is, thus, already equivalent to calling the system biaxial. The projection of an order parameter in the 1 − k plane, where k is either two or three, into the perpendicular plane always gives a different result than measuring the order in the plane itself, which can be seen from Figure 6.
The Smoluchowski theory predicts a shear dependence of the biaxiality parameter T displaying a maximum that is roughly located at the Peclet number where the rotational relaxation time of the rods is equal to the reciprocal shear rate, which is also the inflection point of the curve for λ 1 in Figure 7a. The peak and location of the maximum of T slightly depends on the thickness of the rods. This is not the case for the largest eigenvalue λ 1 , as is seen in Figure7c. In contrast to the experimentally-determined shear rate dependence of λ 1 and θ, the experimentally-determined biaxiality is much less sensitive to shear flow at low shear rates than predicted by theory. This suggests that due to the high degree of entanglement, torque is also transmitted in the vorticity direction, leading to an overall flow alignment. The fact that we merely see a monotonic increase of T and no return to a monoaxial distribution might also be connected to the observation that both λ 1 and θ do not reach their plateau at high shear rates and as such to the dilation of the cage. Again, at this high degree of ordering, the shear rate should rather be scaled by the infinite dilute rotational diffusion coefficient.

Conclusions
The complete shear rate-dependent ordering tensor of sheared isotropic dispersions of quasi-ideal colloidal rods was obtained by combining two rheo-SANS geometries. Thus, we obtained the orientational order parameter λ 1 , the angle of the corresponding eigenvector θ and the biaxiality T, which could be related to the shear thinning viscosity. In order to compare results to theory, we scaled the shear rate with an effective rotational diffusion coefficient that approaches zero at the isotropic-nematic transition. While the process of ordering is qualitatively predicted by DEK, as well as Smoluchowski theory, there are some important differences that hint to the need for theoretical improvements.
At very low shear rates, all theoretical predictions underestimate the zero-shear viscosity and corresponding strong shear thinning, as well as the strong shear rate dependence of the orientational order parameter λ 1 and average angle θ, which exemplifies also the interdependence between ordering and shear thinning. We attributed this to the negligence of concentration dependence in the rotational diffusion coefficient for the low shear rate regime, as well as the use of an equilibrium pair-correlation function. While the first of these issues is apparently well-covered by DEK theory, the second has not received enough attention given that the main difficulties in the high shear rate region could be shown to lie in this function.
At high shear rates, we observe that the applied shear rates were not high enough for λ 1 and θ to reach their plateau, while T did not turn back to a monoaxial value of zero. This suggests that for a high degree of ordering a different scaling of the shear rate is needed, probably with the diffusion coefficient at infinite dilution, as the local surrounding of a rod has completely changed.