|
|
||||||||
1Departments of Mechanical Engineering and 2Biomedical Engineering, Northwestern University, Evanston, Illinois; 3Program in Neuroscience and Cognitive Science, University of Maryland, College Park, Maryland; and 4Institute of Neuroinformatics, University and ETH Zurich, Zurich, Switzerland
Submitted 7 July 2006; accepted in final form 11 May 2007
| ABSTRACT |
|---|
|
|
|---|
| INTRODUCTION |
|---|
|
|
|---|
It is well known that the length of the rat's whiskers varies from long to short along the caudal-rostral dimension (Brecht et al. 1997
; Hartmann et al. 2003
; Neimark et al. 2003
). Thus one plausible mechanism for radial distance encoding is for the rat to compare the identity of whiskers that contacted an object with those that did not. If a whisker of length L touched an object, but a whisker of length L –
L did not, the rat could infer that the object was located at a distance between those two values (after accounting for different whisker base locations). Behavioral studies have shown, however, that rats can determine aperture width with only one whisker remaining on each side of the face (Krupa et al. 2001
). This suggests that cross-whisker comparisons cannot fully explain the rat's distance discrimination capabilities.
A preliminary analysis of the whisker as a cantilever beam suggested that the stiffness properties of the whisker might provide a mechanical explanation for the rat's ability to perform accurate radial distance discriminations. We specifically hypothesized that information about moment at the whisker base is critical for determining radial object distance. To test this hypothesis, we developed two closely related biomechanical models of the whisker. Both models were deliberately developed in analytic form, so that researchers could easily calculate moment at the whisker base during experiments. The analytical models were tested against numerical simulations to quantify limits on their application and together confirmed our hypothesis: by correlating movement to changes in moment at the whisker base the rat could determine the radial distance of an object.
This work continues our characterization of vibrissa dynamics (Hartmann et al. 2003
) and suggests some useful ways to represent the mechanical information encoded in the primary sensory neurons of the trigeminal ganglion (Vg). We interpret these results in the context of our evolving hypothesis that neural responses in Vg can be comprehensively represented using a state-encoding scheme that includes combinations of four and only four mechanical variables at the whisker base: angular position, angular velocity, moment, and rate of change of moment.
| METHODS |
|---|
|
|
|---|
This analysis is based on a total of seven vibrissae obtained from three female Sprague-Dawley rats that had been euthanized in unrelated experiments. All procedures were approved in advance by Northwestern University's Animal Care and Use Committee. Each whisker was grasped firmly at the base and plucked out of the follicle for testing. Visual examination of the whisker revealed that there was a qualitative difference in appearance between approximately the first millimeter of the whisker and the remainder of the whisker. Closer examination under the microscope additionally suggested that this first millimeter is approximately the portion of the whisker that would reside in the follicle, and we therefore used this portion to rigidly attach the whisker to the test stand or load cell during experiments.
Static experiments to determine whisker flexural characteristics
A micromechanical force tester (Mach-1, BioSyntech, Montreal, Canada), was used to impose small vertical displacements on the whisker at known horizontal distances from the base and to measure the associated force. The Mach 1 has a positional accuracy of 1.5 µm, and we used a 50-g load cell to achieve a load resolution of 0.0025 g. This allowed us to characterize force-bending relationships for all but the smallest whiskers. Images of the whiskers as they were deflected during the experiment were acquired with a high-resolution (3,088 x 2,056 pixels) digital camera (Digital EOS Rebel, Canon).
Figure 1A shows the experimental set up used to perform the static force measurements. Whiskers were rigidly fixed at their base to a cylindrical metal test stand using cyanoacrylate (superglue). All whiskers were mounted concave down. A shallow groove 1 mm in length was etched in the top face of the stand. The whisker was placed directly in the groove to ensure that exactly the first millimeter of the whisker was rigidly attached to the stand. As described above, this first millimeter is likely to correspond to the portion of the whisker that would normally reside inside the follicle. Miniature scales (Minitool, Los Gatos, CA) with 100-µm tick-marks were attached both vertically and horizontally to the side of the cylindrical stand. These scales provided an independent measure of displacement that could be compared with the positions given by the Mach-1 micromechanical tester. The inset of Fig. 1A shows a close-up view of the stimulator used to deflect the whisker. The stimulator was custom-machined to a fine taper so that the width that ultimately contacted the whisker was
500 µm.
|
After the whisker had been deflected through a full 1,500 µm, the stimulator was moved back to y = 0 so that it no longer contacted the whisker. The stimulator was moved in the positive x-direction to a different horizontal distance from the base of the whisker. We typically moved the stimulator in 2,000-µm increments in the x-direction, but for some whiskers we moved in 1,000-µm intervals. The stimulator was again positioned carefully just barely above the surface of the whisker, this position was defined as a new y = 0, and the stimulator was lowered to displace the whisker at this new x-location. We continued moving the stimulator further out horizontally from the base of the whisker until we reached the resolution of the force measurement capabilities of the Mach 1 tester.
Dynamic experiments to determine whisker flexural characteristics
In dynamic experiments, the whisker base was mounted directly to the load cell and moved using the Mach-1 tester to hit the tapered stimulator, which was held fixed in position. This experimental setup is shown in Fig. 1B and allowed us to continuously monitor the force at the base of the whisker as it deflected into the stimulator. We lowered each whisker into the stimulator at two different velocities (50 and 500 µm/s) and at five different horizontal locations away from the whisker base (3, 5, 7, 9, and 11 mm). Note that the y-direction is opposite that in Fig. 1A for directional consistency with respect to the whisker.
Analysis of experimental data
Force and displacement data (from the Mach-1 tester), along with the digital images of the whiskers, were imported into MATLAB (v 7.0 2004, The Mathworks, Nattick, MA). As is the convention for load cells, the load measurements from the Mach-1 were provided in grams. These measurements were multiplied by a factor of 9.8 m/s2 to obtain the force in millinewtons (mN). Whisker-stimulator contact forces were always assumed to be normal to the whisker because the contribution of force from friction was assumed to be negligible. The load cell in the Mach-1 tester measured only vertical force, and we therefore divided the measured force by the cosine of the whisker angle at the contact point to obtain the actual force applied.
To extract the geometrical shapes of the whiskers from the high-resolution photographs, the upper and lower outlines of the whisker were located using semiautomated image processing techniques in MATLAB. The shape of the whisker was defined as the average of the upper and lower outlines. For each extraction the averaged points were overlaid on top of the photographed whisker to visually confirm that the averaging technique yielded data points that fell within the upper and lower outlines of the whisker, thus giving an excellent match to the overall shape.
Fundamentals of elasticity: cantilever beam theory
Our goal in this research was to develop an accurate but simple biomechanical model of the rat whisker as a cantilevered beam. Cantilever beam models are derived from elasticity theory (Euler 1744
; Love 1944
; Timoshenko 1970
; Young and Budynas 2001
), which can be used to relate the curvature,
, of a straight cantilever beam to the moment, M, at each point along its length, x
![]() | (1) |
14°), the term
in the denominator is negligible and Eq. 1 can be linearized as
![]() | (2) |
![]() |
In Eq. 2, F is the force exerted normal to the beam at a distance along the whisker, a, from the base of the beam. The linearization assumes that the beam is initially straight and that it deflects only through small angles. This means that the arc length distance a, is essentially the same as a horizontal distance.
If we now assume that the beam is cylindrical with a radius of r then the area moment of inertia I =
r4/4 and Eq. 2 can be solved analytically for y(x)
![]() | (3) |
It is important to note that elasticity theory itself is very general, simply relating curvature to moment. However, for biological materials, Young's modulus (E) is an approximation at best, because these materials are typically anisotropic, heterogeneous, and nonuniform. For a material whose value of E is roughly 5 GPa, the best one might expect is to obtain a value correct to within a few gigapascals.
Comparing results of model 1 with experimental results
Our analysis required a comparison of the shape of the whisker as predicted by our tapered beam model (RESULTS, model 1) with the shape of the real whisker obtained experimentally with high resolution photography. However, real rat whiskers have an inherent curvature. We therefore made the approximation that the deflection of the real whisker (under a force F at a particular arc length location, a) could be expressed as the deflection of a straight tapered cantilever beam (under that same force F, imposed at the same location a), summed with the inherent curved shape of the undeflected whisker (under conditions of zero force). This approximation is schematized in Fig. 2. This analysis is valid as long as the assumptions of the linearized beam model are not violated, namely that the deflections and inherent whisker curvature are sufficiently small. The first part of RESULTS identifies the conditions in which these assumptions are valid.
|
Note that for this model, nodes beyond the point of whisker-stimulator contact deflect linearly, as can be seen mathematically in Eq. 3. In the model, we could therefore assume that the portion of the whisker past the stimulator contact point was translated in the same direction as the last node before the stimulator contact. This portion of the whisker was thus aligned to match the tangent of the deflected whisker at the point of contact.
Numerical simulations
Numerical simulations of whisker bending were performed to identify the limitations on and validate the results of the two analytical models. These numerical simulations also accounted for large angle deflections and inherent whisker curvature. All simulations were performed in MATLAB, and were based on the following principle: if a force
acts at an arc-length a from the base of a beam, the resulting beam shape can be found through repeated application of d
i = (
i x
)/EIi, where d
i is the change in curvature,
i is the vector connecting node i to a, and Ii is the area moment of inertia at node i.
always acts normal to the whisker as long as there is no friction.
| RESULTS |
|---|
|
|
|---|
14°), but the real whisker can bend through very large angles during object contact. The results below account for each of these three complexities and are divided into three parts. In part 1, we develop an analytical model (model 1) to describe the bending of a rat whisker. The model uses the magnitude and location of the imposed force to determine the resultant shape of the whisker after deflection. The model accounts for both whisker taper and inherent whisker curvature and limits on its applicability are tested using numerical simulations.
In part 2, we validate model 1 against experimental data obtained from real rat whiskers, showing an excellent match between theory and experiment. Finally, in part 3, we develop a second analytical model (model 2) that describes the relationship between the rate of change of moment at the whisker base and radial object distance. Numerical simulations are used to show that the inherent curvature of the whisker has a negligible effect on this relationship. We show that measuring changes in moment at the whisker base would permit the rat to extract radial object distance and analyze the consequences of this result for coding in the trigeminal ganglion.
Part 1: Developing an analytical model of a tapered rat whisker with inherent curvature
AN ANALYTICAL EXPRESSION FOR THE DEFORMATION OF A TAPERED WHISKER WITH NO INHERENT CURVATURE.
Expressions for the deflection of a straight cylindrical cantilever beam under a load are readily available in the literature (Young and Budynas 2001
). However, the diameter of a rat whisker decreases approximately linearly with length (Hartmann et al. 2003
; Neimark et al. 2003
). We therefore extended the cylindrical model to account for the taper of the whiskers. The basic derivation for tapered deflections is the same as for the cylindrical case and can be found in APPENDIX A. The analytic solution for the small-angle deflections of a tapered beam was found to be
![]() | (4) |
EFFECTS OF TAPER ON THE SMALL ANGLE APPROXIMATION. The small angle assumption implicit in the linearization of Eq. 1 means that the deformations expressed in Eq. 4 will become inaccurate after a certain bending angle. We used numerical simulations (see METHODS) to explore how linearization impacts the accuracy of the analytical model under large angle deflections and inherent whisker curvature.
Figure 3A shows the difference between the small-angle approximation and the large-angle numerical result for a 200-µN force applied at distances of 10, 20, and 30 mm out along a 60-mm tapered whisker. For the first two locations of applied force, the difference between the small angle approximation (dotted line) and the numerical result (solid line) is negligible. When the force is applied at 30 mm, the small angle approximation clearly diverges from the numerical result as the deflection angle becomes sufficiently large. The effect of taper increases for deflections applied further from the base.
|
To validate the assumptions implicit in the summation, we used numerical simulations to calculate the deflections of an inherently curved whisker through large angles (see METHODS). Before we could compare the analytical results of model 1 and the results of the numerical simulations, however, we noted one additional complexity, as follows: if the whisker is initially straight and deforms only through small angles, the arc length a (the distance as measured along the length of the whisker) differs negligibly from the straight distance from whisker base to point a out along the whisker. This was discussed previously in METHODS. If the whisker is not straight, but instead has an inherent curvature, these two values are different. Thus for the remainder of this paper, it is important to remember that a is always defined as the arc length distance, not the straight distance from base to point of contact distance.
Figure 3B shows the error between deflection profiles found using model 1 and using numerical simulations. The thin solid lines represent straight (black) and inherently curved (gray) whiskers. The inherently curved whisker was chosen to have a constant normalized curvature (ratio between the total arc length and the radius of curvature) of 1. This normalized curvature value is similar to the values found for the whiskers used in this study (data not shown). An increasingly large force was applied at a = 30 mm for both whiskers until the magnitude of deflection at a, y(a), differed by 10% between the two models. The thick solid and dashed lines give the deflected shape of the initially straight (black) and curved (gray) whiskers, as found by using model 1 (dashed) and numerical simulation (solid). It is clear that model 1 yields an accurate description of the deflected whisker up to the force location a, but is less accurate further out.
Figure 3C quantifies the amount of angular deflection that results in 10% error between the two models for a force imposed at any point along the whisker. The inherently curved whisker again had normalized initial curvature of 1, as described for Fig. 3B. The same procedure described for Fig. 3B was repeated for several a values and the resultant deflection angle,
, at which 10% error was reached for each a value was recorded. Figure 3C shows the amount of angular deflection plotted against normalized location of the imposed force, a/L, for an inherently straight (black) and curved (gray) whisker. It is apparent from this figure that imposed force location affects the amount of deflection possible before 10% error results between model 1 and the numerical simulations. As the location of imposed force increases, the amount of deflection before the 10% threshold is reached decreases for both the straight and precurved whiskers. This relationship is steeper for the inherently curved whisker, but both cases show that model 1 is most accurate when forces are applied close to the whisker base.
This analysis has shown that by summing the inherently curved shape of the real whisker with the deformations calculated from Eq. 4, experimenters can obtain an approximation of the deflected whisker shape up to point a with
10% error, provided the force is imposed at a/L < 70%.
Part 2: Validating model 1 against experimental data obtained from real rat vibrissae
Model 1 incorporates the effects of taper and inherent curvature, and we have shown it to be particularly accurate for forces applied close to the base. We used two different methods to determine how well model 1 captured the bending characteristics of a real rat whisker. First, we compared force-displacement curves between model and experiment. Second, we used the model to predict the entire shape of deflected whiskers, and compared this prediction with experimentally-obtained shapes of deflected whiskers.
FORCE DISPLACEMENT CURVES: ANALYTICAL EQUATIONS AND EXPERIMENT. We experimentally quantified bending for a real rat whisker in response to a force imposed at different distances from the base. Figure 4A shows three overlaid images of the E2 whisker bending under the same force (121.8 ± 24.5 µN) imposed at a distance of 7, 8, and 11 mm horizontally from the base. Note that although during experiments the stimulator was positioned at horizontal distances, during all analysis the horizontal distance was converted to arc-length.
|
Figure 4C shows superimposed images of the C3 and
whiskers as they were deflected by approximately the same force (840.3 ± 94 µN) imposed at a horizontal distance of 8 mm. It is clear that the force has a larger effect on the C3 whisker (Dbase = 119 µm, L = 21.50 mm) than on the
whisker (Dbase = 225 µm, L = 66.20 mm). Figure 4D quantifies this effect for seven different whiskers of varying size. In this experiment, whiskers
and
are the longest whiskers, with lengths of 66.2 and 60.3 mm, respectively, whereas E2 and
are the thickest at the base, with base diameters of 232 and 225 µm, respectively. The last four whiskers are all shorter and thinner at the base than
,
, or E2, and therefore require less force to deflect the same amount.
Notably, Fig. 4D shows that at a given horizontal distance away from the base (6 mm in this case) the force-displacement curve follows a linear relationship for each whisker. This relationship can be seen explicitly in Eq. 4. Importantly, this does not mean that for a given force F the whisker will bend linearly along its length, because the proportionality constant between F and y(x) is different at each point x.
CAPTURING THE COMPLETE SHAPE OF A WHISKER: MODEL 1 COMPARED WITH EXPERIMENT. The force-displacement curves in Fig. 4 showed a good match between Eq. 4 and experiment for discrete values of force and displacement. They also serve to quantify the effects of whisker size (base diameter and length) and force location on whisker deflection. However, the curves of Fig. 4 only quantify the relation between force and displacement at point a, where the force is applied. How well can model 1 as described in part 1 characterize the entire shape of the whisker when it contacts an object, purely as a function of whisker length, diameter, and object distance a? To answer this question, model 1 was used to predict the deflection of the whisker everywhere along its length (i.e., at all values of x). These modeling results were then compared with the photographed shape of the whisker (Fig. 1A). Because Fig. 3C shows that the model should remain accurate for relatively large deflections close to the base, it would be surprising if model and experiment were not in good agreement.
We used model 1 to calculate the full shape of the whisker as a function of x analytically while leaving Young's modulus (E) as a free parameter. Experimentally, we took digital photographs to obtain the entire shape of each whisker as it was increasingly deflected by the stimulator. We imported the photographed shape into MATLAB and superimposed the modeling result. The value of E was varied in the model until the best match was found between model and experiment. If E was too large, the model did not deflect enough compared with the experimentally deflected whisker, and if E was too small, the model whisker deflected too much.
The inset of Fig. 5 shows the quantities used to find the best match between model and experiment. The error between the model and the experimental data were found by taking the ratio of the areas between the model and the deflected whisker (area 2) and the area between the deflected whisker and the undeflected whisker (area 1 + area 2). All areas were calculated from the whisker base to the point of contact, a. Normalization to the area between the undeflected whisker and the deflected whisker accounted for any error induced by apparent changes in length due to the small angle approximation, and permitted comparisons of error estimates across whiskers of different lengths. This ratio is referred to as the percent area error, plotted on the y-axis of Fig. 5.
|
, C3, E2, E3,
) were used in the analysis of the complete whisker shape. For each whisker, we averaged over all vertical deflections at each horizontal distance from the base. This amounted to
110 comparisons between experiment and model, >20 values of Young's modulus, for a total of
2,200 comparisons per whisker.
Figure 5 shows the results for the A1 whisker. Plotting error as a function of Young's modulus (E) shows that 1) the smallest error (2.72% for the A1 whisker) is found when the object is closest to the base of the whisker, 2) estimated E for the A1 whisker has a range of
1.5–4.3 GPa, with an average of 2.75 GPa, consistent with the value found for Fig. 4B and with previous estimates (Hartmann et al. 2003
; Neimark et al. 2003
), and 3) the value of the "best" E decreases as the object moves further from the base. These results were representative of all whiskers. A1 does not represent a "best case."
The ranges for E of the other whiskers were mostly similar to that of A1. Table 1 shows geometrical dimensions and average values of E for all seven whiskers. Results for the C3 whisker lay outside the range of results for the other whiskers. For C3, Young's modulus ranged from 4 to 9.5 GPa and had an average value of 6.25 GPa. The C3 whisker was by far the shortest and thinnest of the whiskers and deflections were imposed up to
50% along the whisker length. Most other whiskers only had deflections imposed up to
35% along the length of the whisker. This could help explain the large value of Young's modulus found for the C3 whisker.
|
MODEL 2: AN ANALYTICAL EXPRESSION FOR RADIAL OBJECT DISTANCE AS A FUNCTION OF MOMENT AT THE WHISKER BASE.
Equation 4 describes a relationship between the deflection, y(x), at each point, x, along a tapered whisker and the arc length, a. The value of y(x) is related to a through the force F imposed at point a, the bending stiffness represented by the product EIbase, and the total arc length, L, of the whisker. We asked whether the rat could use the relationship expressed in Eq. 4 to infer information about the radial object distance d, from the whisker base to the contact point. Note that in general the distance d is shorter than arc length distance a, however, assuming a straight whisker and evaluating Eq. 4 at y(d) yields
![]() | (5) |
= y(d)/d (assuming small angle deflections) yields
![]() | (6) |
and LBT is the linear base to tip length of the whisker. Note that L was replaced with LBT to enforce the boundary condition that M = 0 when d = LBT. Solving for the variable d, and taking time derivatives yields
![]() | (7) |
Note that Eq. 7 is expressed in terms of time derivatives. These time derivatives are included primarily for biological plausibility. Recall that
represents the angle that the whisker has rotated since the time of initial contact with the object. M is the moment experienced at the base, which increases as the whisker rotates against the object. In an engineered system, it is easy to set
= 0 at the angle of initial contact and to keep track of its increasing value. In principle, just like the engineered system, the rat could use the absolute position of the whisker (
) combined with an absolute measurement of moment to determine object distance. However, given the well-known difficulty for the nervous system to accurately measure absolute quantities, but its exquisite sensitivity to rates of change, we think it most probable that the rat would use the time derivatives as represented in Eq. 7.
If the whisker moves at constant velocity, then derivatives of moment with respect to
and with respect to time are proportional. If, in contrast, the whisker moves at nonconstant velocity, the rat could keep track of how moment is changing relative to
. Thus most generally, radial distance can be computed as
![]() | (8) |
2° of rotation in the 5 ms after object contact, well within the linear range. It should be noted that the rat will have much less time to compute object distance if contact occurs close to the tip, as the whisker will quickly fold in on itself and/or flick past the object for small
, and
will change accordingly.
Equations 7 and 8 show that, if the rat can keep track of the rate of change of moment and the velocity with which it is "pushing" its whisker against the object, enough information will be present to infer object distance. Taken with the results of previous studies that have described mechanisms for encoding horizontal and vertical position (Ahissar and Arieli 2001
; Ahissar et al. 2000
; Szwed et al. 2006
), these equations effectively show that only three mechanical variables are required to extract 3D spatial information about objects. Those variables are angular position, angular velocity, and rate of change of moment (or curvature). In addition to these three variables, we posit that the rat is sensitive to a fourth variable—moment—so as to remain sensitive to static deflections of its whiskers.
Predicted changes in moment at the whisker base as the whisker rotates against an object
We now use model 2 (Eq. 7) to compute the predicted changes in moment at the whisker base as the whisker is rotated against an object. Figure 6A plots the rate of change of moment at the base of the whisker as a function of contact distance for two different whisking velocities. To highlight the effects of taper (see APPENDIX B), results for the cylindrical whisker are also shown (gray traces). For both tapered and cylindrical whiskers, the steepest change of
is for
, when the imposed force is closer to the vibrissal base. It is clear that
goes to infinity for positions very close to the base. In addition,
goes to zero at the tip of the tapered whisker, meaning that almost no moment is transmitted back to the base when contact is made very near the tip. Instead, the whisker tip might locally deflect and subsequently drag along the object. This suggests that the more distal regions of the whisker may be more sensitive to low-amplitude, high-frequency signals, because these small signals can be amplified by resonance (Andermann et al. 2004
).
|
is larger when the angular velocity
is larger. This is an intuitive result, but the figure makes clear that the rat can obtain the same value of
at the whisker base either by increasing whisking speed or by moving its snout closer to the object. This may suggest the existence of a "sweet spot" or "sweet combination" of object distance and whisking velocity. This location on the whisker would be constrained by the following criteria.
Figure 6B plots the rate of moment change at the base of the whisker as a function of whisker angle,
. As mentioned earlier,
is the angle subtended since initial contact with the object and is interchangeable with time on the x-axis as long as
is constant. Each curve in Fig. 6B represents a different object distance (0.3LBT, 0.6LBT, and 0.9LBT). It is critical to understand that the linear relationship between
and
does not mean that the whisker will bend linearly along its length. The proportionality constant between
(x) and y(x) is different at each point x along the whisker. In addition, imposing a force at position 2x does not make the whisker bend twice as much as if the force were imposed at position x. This can be seen in the uneven spacing of the lines for 0.3LBT, 0.6LBT, and 0.9LBT in Fig. 6B.
Predicted changes in curvature at the whisker base as the whisker rotates into an object
It is clear from Eq. 1 that curvature and moment are directly proportional. The change in curvature at any point along a beam is equal to moment divided by the whisker bending stiffness, EI. Figure 6C plots the angular position of the whisker,
, as a function of normalized contact distance
, for an imposed 0.1-µN-m moment. Simply put, this plot predicts how much the whisker will bend if the whisker is being actuated by a maximum moment of 0.1 µN-m. For a cylindrical beam, the relationship is purely linear. For a tapered beam, much less moment is required for the whisker to deflect past a distal object compared with a more proximal one.
Effects of inherent whisker curvature on moment sensed at the base
Finally, we now show that model 2 holds for all realistic values of whisker curvature (and even much larger curvatures). A mechanical rule of thumb states that if the radius of curvature of a beam is
10 times its maximum cross-sectional depth, many fundamental principles of deformation analysis remain valid (Young and Budynas 2001
). Geometrical analysis showed that the real rat whiskers used in this study exhibited a maximum curvature of 1.7 (units normalized to whisker arc length) along their length. A typical ratio of the radius of curvature to depth was
250. The minimum ratio found along any whisker was
100. Because the minimum value is much >10, fundamental elasticity equations apply.
Numerical simulations were used to compare changes in
profiles for a straight whisker and for a whisker with a large inherent curvature. To be conservative, we modeled the deformation of a whisker bent into the extreme shape of a semi-circle, which has a constant normalized curvature of
3.14. This is roughly twice the maximal curvature found for any of the real whiskers. Base-to-tip length for both models was 60 mm. Figure 6D shows the results of the simulation: the
profiles for the straight and inherently curved whiskers overlap almost exactly. The inherent curvature has negligible effect on the moment that will be sensed at the whisker base.
Effects of whisker or head translations compared with whisker rotations
As established by earlier studies (Szwed et al. 2003
, 2006
), and as schematized in Fig. 7A, cylindrical coordinates are the most natural system to describe whisking movements of the rat. Theta describes the rostral-caudal angle, z is the height of the whisker row, and r is the radial distance out along the whisker. This coordinate system is particularly suited to describe the rotational movements that most typically characterize whisking behavior.
|
Figure 7B shows that the models presented in this paper hold equally well for translation and rotation and can explain the results of the earlier study by Krupa et al. (2001)
. The models also apply to earlier studies that involve small-angle passive displacements of the whiskers in anesthetized rats (Leiser and Moxon 2006
; Lichtenstein et al. 1990
; Shoykhet et al. 2000
; Webber and Stanley 2006
). With knowledge of whisker length and base diameter, approximate Young's modulus (3–4 GPa), the location of the imposed stimulus and its magnitude (which could take the form of a force, rotation or linear deflection), experimenters can now calculate approximately how much moment is experienced at the base of the whisker during passive displacement experiments. As will be shown in DISCUSSION, however, this may not be a very useful calculation to perform for passive displacement experiments.
The variables that this study has found to be important for shape extraction are angular position, angular velocity, moment (or equivalently, curvature), and rate of change of moment. This mechanical analysis suggests that a state-encoding scheme (Paulin 2004
; Paulin and Hoffman 2001
; Paulin et al. 2004
) is a parsimonious and quantitatively rigorous way to represent the responses of Vg neurons. Figure 7C shows an example of a state-encoding scheme using three of the four mechanical variables. Neurons have a certain probability of firing a spike when the whisker is in a particular "state." A state is uniquely defined by whisker position, velocity, moment, and moment-dot; in the example of Fig. 7C only three of the four variables are included. If necessary, velocity could be defined to have two dimensions (rostral-caudal and dorsal-ventral) to account for the directional sensitivity of the cells of velocity information (Jones et al. 2004
). This would result in a higher dimensional space but would otherwise leave the state-encoding representation unchanged.
| DISCUSSION |
|---|
|
|
|---|
WHY DEVELOP AN ANALYTICAL MODEL? This study has developed simple, analytical models for the deformation of rat vibrissae that account for vibrissal curvature as well as taper. The models are well matched by experimental results (Figs. 4 and 5). The advantage of an analytical model over the numerical method also presented in this paper is that it can be solved quickly and exactly, without use of a computer, to obtain a very close approximation to how a real whisker will bend. This is potentially useful to all investigators performing experiments in which the whiskers are deflected by an amount within the confines defined by Fig. 3C. Analytic models also make explicit the dependence of whisker bending properties on mechanical variables. Numerical simulations are required to precisely quantify bending of the whiskers in other cases. It is important to note that the change in curvature and the change in moment at every point along the whisker length are directly proportional, related through the whisker bending stiffness EI.
YOUNG'S MODULUS AND WHISKER STIFFNESS.
Young's modulus (E) for biological materials is an approximation at best. This study found that E approximately equals 3–6 GPa, in line with previous estimates (3–4 GPa, Hartmann et al. 2003
; 9 GPa, Neimark et al. 2003
; 3.5 GPa, Solomon and Hartmann 2006
). A puzzling result of these experiments is that the value for E seemed to decrease as forces were imposed further from the whisker base (Fig. 5). There are at least four possible explanations for this result. First, the result could be taken at face value. The whisker material may vary with length in such a way as to result in lower E values further from the whisker base. Second, it is possible that the equivalent stiffness of the whisker decreases with whisker length. For example, if the whisker tapered parabolically instead of linearly, then the smaller cross sectional area as a function of length would result in an apparently lower E value. Third, Fig. 3C shows that the accuracy of the tapered-beam model decreases as the force is imposed further from the base. It is therefore possible that the decreased accuracy of model 1 is directly responsible for the apparent change in Young's modulus. This is consistent with the increase in error associated with the "best fit" Young's modulus as deflections were imposed at increasing radial distances. Fourth, friction would have the largest effect on the most curved (shortest) whisker, thereby increasing the apparent value for E.
IMPORTANCE OF MOMENT. Almost all previous studies of the vibrissae have focused exclusively on kinematic variables, that is, angular position and its time derivatives. These variables are termed kinematic because they describe the motion of a body without consideration to the forces or moments that affect the motion. During active whisking, however, kinematic variables alone cannot provide a complete representation of all the information transmitted to the rat through its whiskers. Under active whisking conditions, the whisker could well be at the same angular position and yet experience very different moments at the base. Figure 8 shows some of the differences between active and passive whisker displacements. The models presented in this study begin to consider the potentially important role that moment may play in conveying meaningful information to the rat.
|
RELATIVE IMPORTANCE OF WHISKER DIAMETER, LENGTH, CURVATURE, AND TAPER. All equations in this study indicate that moment at the whisker base will depend on the base diameter of the whisker raised to the fourth power. Thus whisker diameter will have the largest influence of any single variable on the moment experienced at the whisker base. There is more tolerance for small deviations in whisker length. The inherent curvature of the whisker plays a relatively small role in determining how the deflected whisker will change shape, whereas in contrast, the taper of the whisker greatly affects how the whisker will bend and the moment transmitted to the base.
Models are highly applicable to natural whisking behaviors
Throughout the METHODS and RESULTS, we have been careful to emphasize the assumptions embedded in the models and the limitations that these assumptions impose. This careful exposition of modeling constraints may leave the impression that the models apply only under very limited conditions. It is therefore important to emphasize that our analysis is in fact very general and that versions of the models will hold even for very complex behaviors.
MODELS CAN APPLY TO A WIDE RANGE OF BOUNDARY CONDITIONS: THE IMPORTANCE OF INSTANTANEOUS MEASUREMENT.
Moment at the whisker base will vary depending on how stiffly or loosely the whisker is held in the follicle, that is, on the boundary conditions in and near the follicle. The rat could presumably change follicular boundary conditions through muscular activation as well as by modulating blood flow to the follicular sinus (Scott 1955
). The models in this study are based on clamped boundary conditions at the whisker base, but more realistic, tissue-like conditions might be modeled with a spring-mounted or a torsional-spring-mounted whisker. It is critical to note, however, that the fundamental results of this study will not change, even if boundary conditions are very different from the clamped condition modeled here. This is because the relationship between moment at the base and radial object distance will remain monotonic regardless of boundary conditions. As long as the rat can learn the monotonic function that relates these variables (M and d), the method proposed here will work for radial distance extraction.
What happens if the rat changes the boundary conditions at the whisker base during the course of a whisk? Equations 7 and 8 show that the rat can determine radial object distance based on the instantaneous rate of change of moment. This means that the rat need only sense distance at a single instant during the whisk, and it does not matter if boundary conditions change before or after that instant. Recent behavioral data (Mitchinson et al. 2007
) have shown that rats often use an exploratory strategy of "minimum impingement," in which they tap, rather than sweep, their whiskers over objects. This suggests that the rat gains a sense of radial object distance in the first few milliseconds immediately after object contact. This strategy is consistent with the one determined to be most effective for radial distance extraction in a hardware model of the whiskers (Solomon and Hartmann 2006
) and also helps avoid measurement complications caused by whisker slip along the object. Finally, we note that regardless of boundary conditions, the amount that the moment will change in a given time interval is directly related to the whisking velocity. We therefore suggest that variations in velocity over the trajectory of the whisk may be of particular behavioral importance to the rat during tasks that require estimates of object distance.
MODELS CAN APPLY TO A WIDE RANGE OF ANGULAR DISPLACEMENTS, VELOCITIES, AND DISTANCES TO OBJECT CONTACT.
Numerous papers have shown that naturalistic rat behaviors use a large range of angular positions, velocities, and distances to object contact (Brecht et al. 1997
; Carvell and Simons 1990
, 1995
; Guic-Robles et al. 1989
; Polley et al. 2005
; Vincent 1912
). It might therefore be asked how the values for these variables presented here fit into these ranges. For example, over what range of angles, whisking amplitudes, and velocities, do the proposed models apply? The short answer is that the fundamental results of the models hold over virtually all distances to contact except very near the tip, all angular velocities, and all angular displacements. Figure 8 shows the broad applicability of the models and the differences between passive displacements and active whisking.
The rat-centered coordinate system for Fig. 8 is defined by
in the top left corner. A value of
= 0 means that the whisker is completely retracted, pointed directly backwards toward the tail of the rat. A value of
= 180° means that the whisker is completely protracted, pointed directly forward toward the snout of the rat. The first row of Fig. 8 shows passive deflection assuming that the whisker behaves as a flexible beam. In this case, pushing a point on the whisker backward or forward causes the whisker to bend and generates a moment at the whisker base. Consistent with the models presented in RESULTS, this figure assumes that the whisker is held rigidly at the base. Assume that the point on the whisker in contact with the stimulator is pushed to some value of
, different from the whisker's rest position. The amount of whisker bending, and hence the moment generated at the base, depends directly on
–
rest, that is, on the position to which the whisker is pushed. This means that there is no way to "decouple" the absolute angular position of the whisker (as measured at the point of stimulator contact) from the moment generated at the base.
The second and third rows of Fig. 8 show that active whisking permits decoupling of the values of absolute whisker position
and the moment generated at the base. In the second row, the whisker is actively protracted forward and behaves as a rigid body until it encounters the object at
90°.