## Abstract

*In-situ* micro/nano characterization is an indispensable methodology for material research. However, the precise *in-situ* SEM twisting of 1D material with large range is still challenge for current techniques, mainly due to the testing device’s large size and the misalignment between specimen and the rotation axis. Herein, we propose an *in-situ* twist test robot (iTRo) to address the above challenges and realize the precise *in-situ* SEM twisting test for the first time. Firstly, we developed the iTRo and designed a series of control strategies, including assembly error initialization, triple-image alignment (TIA) method for rotation axis alignment, deformation-based contact detection (DCD) method for sample assembly, and switch control for robots cooperation. After that, we chose three typical 1D material, i.e., magnetic microwire Fe_{74}B_{13}Si_{11}C_{2}, glass fiber, and human hair, for twisting test and characterized their properties. The results showed that our approach is able to align the sample to the twisting axis accurately, and it can provide large twisting range, heavy load and high controllability. This work fills the blank of current *in-situ* mechanical characterization methodologies, which is expected to give significant impact in the fundamental nanomaterial research and practical micro/nano characterization.

## Introduction

Since the mid-nineteenth century, mechanical test has been an irreplaceable methodology for material study owing to its advantages in offering intuitionistic result. Benefiting from the micro/nano technologies, the *in-situ* SEM mechanical characterization has also become possible in recent decades. Compared with traditional methods, *in-situ* SEM test is able to provide a complete picture of the evolution of micro/nano specimen as it is loaded and fails owing to its real-time imaging ability at high magnification. At current stage, scientists have successfully realized the *in-situ* indention^{1,2,3}, bending^{4} and tension (compression)^{5, 6} test of micro/nano materials, which greatly promotes our fundamental understanding on the micro/nano world.

Twisting test plays an indispensable role among all the three widely accepted mechanical test approaches in material research, i.e., stretching (compressing), bending and twisting. From the view of fundamental study, the grain boundary sliding and crystal twinning phenomenon of material can be observed much easily under twisting^{7, 8}. On the other hand, from the view of application, many products and components are subjected to torsional forces during their operation, making twisting test be a powerful approach for device evaluation. To twist the small sample, a commonly used approach is to upgrade the traditional torsion test instrument in accuracy^{9, 10}. However, limited to the inherent drawbacks of the traditional instruments in structure and working mechanism, it is hardly possible to align the sample perfectly into the rotation axis of the rotation table, so that the effective gauge need to be lengthened to reduce the influence of the vertical misalignment^{11,12,13}. It not only affect the measurement accuracy, but also make it nearly impossible to integrate within SEM for *in-situ* twisting test. Recently, a nanotube-based torsional NEMS device is reported for nano twisting test^{14,15,16}, which overcomes the drawback of the traditional instrument in both size and accuracy. However, this system only works for very few specific samples due to its limit rotation range, torsion load, and inefficient operation mode.

The rapid progress of robotics, especially micro/nano robotics, provides new opportunities for the precise and efficient operation at small scale^{17,18,19,20}. Nowadays, robot-aided *in-situ* SEM characterization has been regarded as one of the most powerful mechanical test methods at small scale owing to its real-time imaging ability at nanometer resolution^{21,22,23}. However, current nanorobotic manipulation systems for SEM are mainly subjected to linear motion. Therefore, although they are capable for the *in-situ* stretching (compressing) and bending test effectively, they cannot meet the requirement of twisting test needing precise rotational operation.

The rotational operation inside SEM raises many serious challenges, mainly including the high manipulation flexibility in narrow space; accurate sample rotational positioning in 3D while only 2D SEM surface image information is available; system cooperation for alignment, assembly and twisting operation at microscopy environment. Herein, in this paper, we develop a nanorobotic manipulation system, named *in-situ* twisting test robot (iTRo), to deal with the above challenges. Firstly, the iTRo with six degree-of-freedoms (DOFs) and two independent manipulation units is developed and integrated within SEM. Then, a series of operation and control strategies are proposed to conduct the assembly error calibration, sample alignment, sample assembly and twisting operation. Lastly, torsion test of magnetic microwire, glass fiber, and human hair is implemented and the experimental results are discussed.

## Results

### Development of iTRo

The iTRo for SEM mainly consists of two independent 3-DOF manipulators, i.e., left manipulator (LM) and right manipulator (RM) on the same basement stage, and a SEM vision feedback system. LM is composed of three independent nanopositioners: two linear positioners and one rotary positioner; and RM is composed by three independent linear positioners. On the end of LM and RM, two sample holders are designed to hold the two ends of sample for twisting test. Additionally, a specific connection port is also designed to control the robot by PC outside SEM’s chamber (Fig. 1a,b and Movie 1).

Here, four steps are designed for a fully *in-situ* twisting test procedure, including robot initialization, LM alignment, RM alignment and twisting operation, which are automated by a switch control strategy (Fig. 1c). In each step, the operation of each nanopositioner is controlled automatically by real-time SEM image processing.

### Robot Initialization

As illustrated in Fig. 2a and Supplementary Fig. 1a, ten coordinates are established for iTRo. The world coordinate {U} is established on the SEM stage and the coordinate {O} is established on the sample rotation stage of SEM. The coordinates {A}, {B}, {C}, {D}, {E} and {F} are established on the six nanopositioners; {P} is established on the endpoint of the sample and {Q} is established on the endpoint of the RM sample holder; {M} is established on the microscopy image system. According to the structure of iTRo, it can be seen that the rotatory assembly errors may come from two joints in coordinates {O} and {A}, which are indicated by *Φ* and *Ψ*, respectively. Because the rotational errors would significantly affect the operation accuracy and even lead to operation failure, they must be calibrated to initialize the robot before twisting test.

In the robot initialization process, the sample’s position in {M} and {P} are obtained from SEM image and robot firstly. After that, the rotation matrix \({R}_{LM\_{\rm{calibration}}}\) between {M} and {P} is calculated based on robot kinematic model (Details are available in methods section). As a result, the initial assembly errors, *Φ* and *Ψ*, are calibrated as *Φ* = 0.0327 rad and *Ψ* = 0.1324 rad finally. After such calibration, we adjust the robot to make *Φ* and *Ψ* infinitely closed to 0 degree. In this case, a standardized rotation matrix \({}_{M}{}^{P}R\) between {P} and {M} corresponding to {U} can be simply represented as a standard matrix:

where *f*
_{
x
} and *f*
_{
y
} are the digital resolutions of the microscope in the *X*
_{
M
}- and *Y*
_{
M
}- axes. Similarly, the rotation matrix between {Q} and {M} corresponding to {U} can also be calibrated and standardized as \({}_{M}{}^{Q}R\).

### Sample alignment on LM

To implement twisting test precisely, the misalignment \({\rm{\Delta }}{P}_{err}={[\begin{array}{ccc}{\rm{\Delta }}{x}_{errP} & {\rm{\Delta }}{y}_{errP} & {\rm{\Delta }}{z}_{errP}\end{array}]}^{{\rm{T}}}\) between sample’s rotation axis and LM’s rotation axis must be small enough (Supplementary Fig. 1). However, in actual practice, it’s nearly impossible to make Δ*P*
_{
err
} be acceptable by manually operation. Herein, we use a triple image alignment (TIA) method to calibrate the offset error and then implement the alignment automatically.

As illustrated in Supplementary Fig. 1b, *P* (*P*
_{
x
}, *P*
_{
y
}, *P*
_{
z
}) is defined as a certain point on the sample’s own rotation axis. Because *Z*
_{
P
}- axis is parallel to the rotation axis, the value of Δ*z*
_{
errP
} would not affect the twisting accuracy. Hence, Δ*z*
_{
errP
} is simplified ignored in sample alignment, i.e., to define Δ*z*
_{
errP
} = 0. At initial condition, the sample’s position *P*
_{
o
} (*P*
_{
x
}, *P*
_{
y
}) in {M} is obtained from SEM image firstly. After that, we rotate the sample clockwise by α degrees, and then anticlockwise by 2α degrees. During this process, the coordinate differences along *Y*
_{
M
}- axis, i.e., Δ*y*
_{
f
} and Δ*y*
_{
b
} (unit: pixel), are calculated by image processing (Supplementary Fig. 1c,d). After mathematical deduction, the offset of sample to LM’s rotation axis can be represented as (details are available in methods section, Supplementary Table 1):

The rule in equation (2) indicates that the sample can be aligned to LM’s rotation axis theoretically regardless of the initial offset.

In this experiment, we set α = 90°, and *f*
_{
x
} = 0.2908, *f*
_{
y
} = 0.2885 corresponding to microscopy magnification 500 × . The values of Δ*y*
_{
f
} and Δ*y*
_{
b
} are directly calculated from the feedback SEM image. By substituting these values into equation (2), we obtain the initial offset \({\rm{\Delta }}{P}_{err}={[\begin{array}{ccc}-69.5 & -257.9 & 0\end{array}]}^{{\rm{T}}}{\rm{\mu }}{\rm{m}}\), which is then used to adjust LM to eliminate the eccentric error.

As the comparison results shown in Fig. 2b, the sample almost maintains on the same straight line while it’s at angle θ + α, θ, and θ − α. To further evaluate the alignment quality, we define two more parameters, i.e., peak value *e*
_{
P
} and standard deviation *e*
_{
SD
} as:

where *i* indicates the *i*th images captured during LM’s rotation in 360°, *j* is the total number of captured image, *y*[*i*] is sample’s fluctuation displacement along *y* axis shown in the *i*th image. *e*
_{
p
} represents the sample’s largest error from the mean value and *e*
_{
SD
} represents the dispersion of sample’s error about *y*- axis while LM rotates in one revolution (360°). During LM alignment experiment, we capture the SEM image every 30°, so that *j* = 12. The results show *e*
_{
p
} and *e*
_{
SD
} is 258 μm and 197 μm respectively at the initial condition. After TIA process in one iteration, *e*
_{
p
} reduces to 22 μm and *e*
_{
SD
} reduces to 11 μm, which improves 91.5% and 94.4% respectively (Fig. 2c). As to human hair, the sample’s fluctuation is already much less than sample’s diameter (~50 μm), which accuracy is acceptable for twisting test. Please be noted that the accuracy can be improved further by repeating TIA process for multiple times.

### Sample assembly on RM

After the sample is aligned onto the rotation axis of LM, the next challenge is to assemble the other end of it onto RM. Assume the sample holder of RM is located under the sample at the initial condition (Supplementary Fig. 1b), and the offset error between them is defined as \({\rm{\Delta }}{Q}_{err}={[\begin{array}{ccc}{\rm{\Delta }}{x}_{errQ} & {\rm{\Delta }}{y}_{errQ} & {\rm{\Delta }}{z}_{errQ}\end{array}]}^{{\rm{T}}}\). SEM could provide clear 2D surface image with nanometer resolution, thereby the relative position of sample and RM stage along *X*
_{
M
}- and *Y*
_{
M
}- axis can be seen directly. In actual practice, the displacement error on *X*
_{
M
}-*Y*
_{
M
} plane would not affect the sample assembly accuracy, thereby we ignore Δ*x*
_{
errQ
} and Δ*y*
_{
errQ
} during RM alignment process, i.e., define \({\rm{\Delta }}{x}_{errQ}={\rm{\Delta }}{y}_{errQ}=0\). Yet, for the distance difference along Z_{M}- axis (depth direction), it’s hard to know from SEM image directly and must be calibrated.

To calibrate Δz_{
errQ
}, we propose a deformation-based contact detection (DCD) approach. According to the imaging principle of SEM, the electrons on the sample surface would change when sample deforms, which would then cause a brightness change on sample surface. Here, DCD approach is designed based on the above principle: we move RM stage until the highlighted area on sample surface changes, which means RM stage contact and deform the sample. The advantage of DCD approach lies in it only requires a single microscope information to determine the contact state in 3D space.

Here, we first define the highlighted area function *S*
_{
i
} of a binaried SEM image,

where *f*(•) is a nonlinear function to exert a binary image processing, *I*
_{
i
}(*x*, *y*) is a M × N gray image captured from SEM vision system in the *i*th step. The DCD rule is defined as:

when

where *l*
_{step} is a constant step size of the *Z*
_{
M
}- axis linear positioner on RM, which is set as 6.67 μm in this paper, ε is a given judgement criterion, which is defined as 4% of *S*
_{0}. As the above rule shows, if the highlighted area change Δ*S* satisfies with equation (7), Δ*z*
_{errQ} is of a very small value, less than 0.5*l*
_{step}, and then can be treated as zero. (More details are available in methods section)

As shown in Fig. 2d, corresponding to the black line in Fig. 2e, the total highlighted area is 20210 at the initial condition (step 1), and changes to 26119 (step 4) as RM moves up (along *Z*
_{
M
}- axis) gradually with *l*
_{step}. In this condition, equation (7) is satisfied (ΔS = 5909 > 20210 × 4%) and then the step size of the positioner is halved (3.33 μm). Next, RM stage starts to move down for one step, so that the error can be controlled within 3.33 μm, which can be ignored compared with the sample’s size. After such alignment, the end of the sample is fixed on the sample holder of RM by glue (Aron alpha, Toagosei Co. Ltd.) finally. (Data are available in Supplementary Table 2).

###
*In-situ* SEM twisting test

Three types of 1D materials, magnetic microwire, glass fiber and human hair are chosen for the *in-situ* torsion test inside SEM, which are on behalf of crystal, amorphous and organic material, respectively.

Firstly, a magnetic microwire (Fe_{74}B_{13}Si_{11}C_{2}) is aligned and assembled on iTRo following the above processes. Then, we twist the sample in one direction continuously with angular velocity 9.5 deg/s until it’s broken at twisting angle θ = 565° (Fig. 3a,b,c,d). As shown in Fig. 3e, the fracture plane is perpendicular to the rotation axis, clearly indicating it’s a typical ductile failure caused by shear stress. In addition, taking into consideration of the gauge length *L* = 1400 μm and sample radius *R* = 8 μm, we can obtain the shear strain is approximate *R* = θ × *R*/*L* = 5.6%, which is a rational value for ferrous alloy.

Secondly, a single-mode optical fiber (G.652.D) with diameter approximate 100 µm is taken as for twisting test. After the fiber is aligned on iTRo with gauge length 1000 μm, the fiber is twisted with a constant angular velocity 9.5 deg/s until is broken at 85°. It can be seen that the fracture plane is perpendicular to the maximum tensile stress (at 45°) approximately (Fig. 4), i.e., helical fracture, which agrees well with the brittle property of glass.

For human hair, we twist it forward-backward for several circles rather than in one individual direction considering its high elasticity. The length between the two ends is approximate 680 μm, as shown in Fig. 5a. The sample is twisted by rotating LM with angular velocity 9.5 deg/s and rotation range ± 180° by 25 circles (Movie 1). Figure 5b,c show images of the sample from 0° to 180° and from 0° to −180° with every 60° in the first twisting circle, clearly demonstrating the torsional deformation on the hair surface (More results are available in Supplementary Fig. 2).

Unlike the traditional materials (crystal and amorphous), as a kind of biomaterial, human air shows more complicate behaviors during twisting. Firstly, we investigate the growth of typical cracks on hair surface. As shown in Fig. 5d, four cracks (I, II, III, IV) on the hair surface are observed by SEM and the images at circle 0th, 5th, 10th, 15th, 20th and 25th are selectively given. Besides, the shape and size change of the cracks with time are also quantitatively given in Fig. 5e,f. The results demonstrate a similar behavior of sample under twisting at macro scale, that the cracks expand gradually with increasing of twisting cycle. For instance, the area of crack IV increased from 10 μm^{2} to 18 μm^{2} after 25 circles, approximate 80% (Data are given in Supplementary Table 3). This phenomenal agrees with the behavior of material under shear stress causing by torsion.

In addition, as shown in Fig. 5g (Data are available in Supplementary Table 4), the slop angel of crack V and VI is 48° and 121° at the initial condition. However, when the hair is twisted from +180° to 0°, the slop angle of crack V and VI increases to 52° and 124° respectively. On contrary, while it is twisted from −180° to 0°, the slop angle of crack V and VI reduces to 39° and 119°, respectively. The above angles are measured in each twisting circle and the results show convinced statistical significance. It indicates the deformation of hair lags the applied extra load, which is an obvious slow rebound effect, which agree with torsional-hysteresis characteristics of human hair in published results^{12}.

We also find that the diameter of hair decreases slightly with increase of twisting circles. As shown in Fig. 5h (Data are given in Supplementary Table 5), taken sample b as an example, its diameter reduces from 47 μm to 45 μm after 25 circles. Human hair is a kind of composite biomaterial with complicated multi-layer and porous structure, and has been reported with distinct viscoelasticity^{24}. During twisting, the soft components in the hair would be compressed, so to make the diameter become a little smaller. Besides, because of the slow rebound effect of hair, the diameter could not recover to its initial condition immediately. Therefore, the diameter decrease slightly with increasing of twisting circle. This explanation is also verified by a further experiment (Supplementary Fig. 3), that after twisting for 25 circles, the diameter of the hair can recovery (increase ~1.5 μm) after 25 mins standing.

## Discussion

Stretching (or compressing), bending and twisting are three widely accepted mechanical characterization methodologies in material research. However, unlike the other two types of test, *in-situ* SEM twisting test is still very challenge at small scale, since the rotational positioning raise higher demands to the flexibility, accuracy of the assembly, alignment and manipulation. In the proposed iTRo (Fig. 1), we design two independent manipulation units with 6-DOFs in total to ensure the flexibility and employ the nanopositioner to ensure the positioning accuracy. As for SEM, it could only provide the 2D grayscale images. Hence, most of the traditional robotic sensing and positioning techniques are not available, especially for rotational positioning requiring 3D manipulation. Therefore, taking into consideration of the specific features of SEM imaging system, we propose a TIA strategy for sample alignment on LM and a DCD approach for sample alignment on RM (Fig. 2), by which the sample can be assembled precisely. In addition, we also develop a switch control strategy to automate the torsion test process and to guarantee the efficiency (Fig. 1c). As a result, we are able to conduct the *in-situ* torsion test inside SEM effectively (Figs 3, 4 and 5, and Movie 1).

To demonstrate the versatility of iTRo, we test three typical materials: magnetic microwire, glass fiber and human hair. As the results show, the metal microwire’s ductile failure, glass fiber’s helical fracture, and human hair’s crack area extension, diameter change, and slow rebound effect are investigated successfully, which verify the effectiveness and efficiency of iTRo, including the loading ability, twisting range, accuracy and controllability. Specifically, we find some unique behaviors of human hair during the twisting operation. For instance, we can conclude that the cross section of hair is not of a circular shape but more like an irregular ellipse based on the multi-view images (Fig. 5b,c and Supplementary Fig. 2), which is also proofed by the cross-section image of hair (Supplementary Fig. 4). In addition, we also find an obvious slow rebound effect in human hair: 1) the twist deformation of hair lags the twisting angle, which is demonstrated by the slop change of the crack in twisting process (Fig. 5b, Supplementary Table 4); 2) the diameter of hair reduce gradually with increasing of twisting circle (Fig. 5c, Supplementary Table 5), but is able to recovery again after standing for a few minutes (Supplementary Fig. 3). Human hair is a kind of composite biomaterial with complicated multi-layer and porous structure, and has been reported with distinct viscoelasticity^{24}. During twisting, the soft components in the hair would be compressed, and the porous structure would also act as damper to slow down the deformation. Thus, the hair shows the slow rebound effect during twisting.

The iTRo provides us the possibility to study the micro/nano material’s behavior under torsion load under SEM *in-situ*, which fills the blanks in the *in-situ* micro/nano material characterization. Except for 1D material, this system could also be used to test 2D material and 3D structures under torsion load. In addition, this system can be easily expended to study the material under composite load, including stretching (compressing), bending and twisting, owing to its high manipulation flexibility and accuracy.

Torque sensing is also a big challenge in micro/nano twisting test^{11, 25,26,27}. However, this paper mainly focuses on solving the operation problems in *in-situ* SEM twisting test, i.e., precise assembly and positioning, especially the misalignment problem in existing twisting test approach^{11, 13}. Up to now, we have been able to investigate the sample’s behavior during twisting and study the fracture plane after twisting. As future work, we will cooperate with other researches in sensing field to integrate iTRo with torsion sensors, and then to investigate the micro/nano material’s property more deeply.

## Methods

### Nanorobotic manipulators

LM and RM are assembled on the same basement stage, and they are independent to each other. LM is mainly composed of three nanopositioners: two linear positioners (ECS3030, Attocube Inc.) and one rotary positioner (ECR3030, Attocube Inc.); and RM is composed by three linear positioners (ECS3030, Attocube Inc.). The travel range, resolution, repeatability, maximum drive velocity and maximum load of the linear positioner are 20 mm, 1 nm, 50 nm, 4.5 mm/s and 90 N respectively. The travel range, resolution, repeatability, maximum drive velocity and maximum load of the rotary positioner are 360° endless, 1 μ°, 5% over full range, 10°/s and 20 N respectively.

### Sample preparation

The magnetic microwire consists of two layer initially, the inner metal core (Fe_{74}B_{13}Si_{11}C_{2}) and the coated shell (glass). In the twisting experiment, we twist the microwire directly without removing the glass shell. Because the glass is more brittle than metal, so it is teared off before metal broken during twisting. Thus, only the inner metal core is left after fracture. The hair is taken from the experimenter, i.e., the author of this paper. The diameter of the hair is approximate 40–60 μm regarding to various sample. In all experiments, the hair is tested without any chemical or physical treatment. In the cross section observation (Supplementary Fig. 4), the hair is cut by a scissor.

### Center line recognition

Considering the nature of twisting test, the center line of the sample should be taken as the reference during the experiment. Taken a human hair (diameter ~50 μm) as example, its image is captured from the microscopy firstly (Supplementary Fig. 5a). Then, Gaussian smoothing filter is adopted to reduce the noise. After that, the edge of sample is recognized by Canny detection, as the red dots shown in Fig. S5B. To find out the center line, the middle point of the sample *y*
_{center} is calculated firstly by:

where *TP* is the height of the full screen, *y*
_{
m
} is the distance between sample’s upper edge and image’s upper edge, *y*
_{
n
} is the distance between sample’s bottom edge and image’s bottom edge, respectively. In addition, *k* is the number of pixels extracted from the top/bottom edge of the sample, which is defined as *k* = 30 in our experiment.

After the middle point is calculated, a straight line is drawn to indicate the center line of the sample, as the yellow line shown in Supplementary Fig. 5b. This approach is able to ensure the recognition accuracy as high as to two pixels, which is adequate for the automatic robot manipulation.

In material study, the sample for twisting test is usually prepared to a symmetrical structure. Hence, the above center line recognition process can work for most cases in actual practice. Moreover, note that the application scope of TIA approach is not limited to symmetrical sample, since it’s based on a to-be-aligned point. For the unsymmetrical sample, the reference point on twisting axis can be pointed out by operator first. After that, the centerline can be set following the same method mentioned above.

### Details for calibrating the initial assembly error

Based on the robot kinematic model, the rotation matrix between {P} and {M} corresponding to {U} can be represented as:

where *f*
_{
x
} and *f*
_{
y
} are the digital resolutions of the microscope in the *X*
_{
M
}- and *Y*
_{
M
}- axes.

Define we have *n* image positions \({S}_{Mi}={[\begin{array}{ccc}{x}_{Mi} & {y}_{Mi} & {z}_{Mi}\end{array}]}^{{\rm{T}}}\) of the sample measured from the microscopy image and the corresponding sample positions \({S}_{Pi}={[\begin{array}{ccc}{x}_{Pi} & {y}_{Pi} & {z}_{Pi}\end{array}]}^{{\rm{T}}}\) obtained from the robot movement, i = 1,2, …*n*. Then, the initial rotatory errors *Φ* and *Ψ* can be calculated by solving the following matrix,

Through calculating the inverse matrix of \([\begin{array}{cccc}{S}_{P1} & {S}_{P2} & \mathrm{...} & {S}_{Pn}\end{array}]\) by the Moore-Penrose matrix inverse method, the rotation matrix \({R}_{LM\_\mathrm{calibration}}\) can be obtained by the following equation:

Note that the value of n depends on the required calibration accuracy of the system, the larger the value, the greater the precision. In this manuscript, we chose *n* = 6 under the magnification of 500×, and then \({R}_{LM\_{\rm{c}}{\rm{a}}{\rm{l}}{\rm{i}}{\rm{b}}{\rm{r}}{\rm{a}}{\rm{t}}{\rm{i}}{\rm{o}}{\rm{n}}}\) is calculated as:

As a result, the initial assembly error *Φ* and *Ψ* are calibrated as *Φ* = 0.0327 rad and *Ψ* = 0.1324 rad, respectively.

### Process of the TIA approach for LM alignment

As illustrated in Supplementary Fig. 1, at the initial condition, define the position of the sample in the rotation axis coordinate {\(\widehat{O}\)} is \({P}_{0}={[\begin{array}{ccc}{x}_{{P}_{0}} & {y}_{{P}_{0}} & {z}_{{P}_{0}}\end{array}]}^{{\rm{T}}}\), and in the microscopy coordinate {M} is \({P}_{0\_M}={[\begin{array}{ccc}{x}_{0} & {y}_{0} & {z}_{0}\end{array}]}^{{\rm{T}}}\). The relationship between *P*
_{0} and *P*
_{0_M
} can be represented as:

After that, if the sample is rotated clockwise by α degrees, we will have two new position of the sample, i.e., \({P}_{1}={[\begin{array}{ccc}{x}_{{P}_{1}} & {y}_{{P}_{1}} & {z}_{{P}_{1}}\end{array}]}^{{\rm{T}}}\) in {\(\widehat{O}\)} and \({P}_{1\_M}={[\begin{array}{ccc}{x}_{1} & {y}_{1} & {z}_{1}\end{array}]}^{{\rm{T}}}\) in {M}. The relationship between *P*
_{1} and *P*
_{1_M
} can be represented as:

Similarly, after the sample is rotated anticlockwise by 2α degrees, the relationship between *P*
_{2} and *P*
_{2_M
} can be represented as:

Based on the robotic kinetic mode, the relationship between *P*
_{0} and *P*
_{1}, *P*
_{0} and *P*
_{2} can be represented as:

where *R*(*α*) is the rotation matrix:

As illustrated in the microscopy image (Supplementary Fig. 1), we get:

and

where the coordinate difference along *Y*
_{
M
}- axis, i.e., Δ*y*
_{
f
} and Δ*y*
_{
b
}, can be calculated based on image processing. During the rotational movement, Δ*z*
_{
f
} and Δ*x*
_{
b
} equal to zero. By solving the above equations (13)–(20), the initial eccentric position *P*
_{
0
}, or say Δ*P*
_{
err
} can be obtained as equation (2):

Based on the above process, the offset Δ*P*
_{
err
} of the sample from the rotation axis of LM can be calculated easily by rotating the robot for α and −α from the initial condition^{28}. Details for the alignment process are illustrated in Supplementary Table 1.

### Process of the DCD approach for RM alignment

The process of DCD approach for RM alignment is illustrated in Supplementary Table 2. During the experiment, the robot moves upward with a constant step size and the highlighted area change Δ*S* is calculated dynamically in each step. Once the value of Δ*S* is larger than the criterion ε, which means the robot touches and deforms the sample, the robot would move backward for half-step and then stop. In this paper, the step size is chosen as 6.67 μm. Note that if a smaller step size is chosen, the sample alignment accuracy of Δ*z*
_{
errQ
} can be improved further.

The calculation of highlighted area *S*
_{
i
} is conducted by image processing. Firstly, the captured grey-scale SEM image (1024 × 800) is converted to binary image. In this paper, we adopt quasi-auto threshold for every sample and then a threshold value is achieved by the image process system. To judge the highlighted area by the same standard, the threshold must be served as a fixed value in the following alignment process. After binaryzation, the highlighted area in SEM image is recognized by contour extraction using Suzuki algorithm^{29}. The algorithm chooses external contour of each highlighted area and store nearly all the points in vector. The distance between two close contour points \(C{P}_{1}\,({{\rm{x}}}_{CP\_1},{{\rm{y}}}_{CP\_1})\) and \(C{P}_{2}\,({{\rm{x}}}_{CP\_2},{{\rm{y}}}_{CP\_2})\) should not exceed 1 pixel, which means:

At last the algorithm will calculate the close highlighted area by Green’s theorem:

where *H* is a positively oriented, piecewise smooth, simple closed curve in a plane, and let *G* be the region bounded by *H*.

Lastly, the highlighted area *S*
_{
i
} is calculated by counting all the circled pixels. For the sample of human hair, the value *ε* is defined as 4% of the initial highlighted area *S*
_{
0
} (unit: pixel). Note that for other types of samples with totally different imaging properties, the value of *ε* may need slight adjustment experimentally to achieve higher accuracy. In this paper, we didn’t chose the amount of the highlighted pixels in a binary image, but added one more step to find the contours in the binary image and took the circled areas as the highlighted area. Compared with accounting highlight pixels in a binary image, the contour border extraction method of highlighted area could eliminate the random noise and detect the deformation more exactly (Supplementary Fig. 6).

### Data and materials availability

All data needed to evaluate the conclusions in the paper are present in the paper and/or the Supplementary Materials. Additional data related to this paper may be requested from Y. Shen (yajishen@cityu.edu.hk).

## References

- 1.
Gouldstone, A., Van Vliet, K. J. & Suresh, S. Nanoindentation: Simulation of defect nucleation in a crystal.

*Nature***411**, 656–656, doi:10.1038/35079687 (2001). - 2.
Minor, A., Morris, J. Jr. & Stach, E. Quantitative

*in situ*nanoindentation in an electron microscope.*Appl. Phys. Lett.***79**, 1625–1627, doi:10.1063/1.1400768 (2001). - 3.
Lee, C., Wei, X., Kysar, J. W. & Hone, J. Measurement of the elastic properties and intrinsic strength of monolayer graphene.

*science***321**, 385–388, doi:10.1126/science.1157996 (2008). - 4.
Chen, C., Pei, Y. & De Hosson, J. T. M. Effects of size on the mechanical response of metallic glasses investigated through

*in situ*TEM bending and compression experiments.*Acta Mater.***58**, 189–200, doi:10.1016/j.actamat.2009.08.070 (2010). - 5.
Yu, M.-F.

*et al*. Strength and breaking mechanism of multiwalled carbon nanotubes under tensile load.*Science***287**, 637–640, doi:10.1126/science.287.5453.637 (2000). - 6.
Kiener, D., Hosemann, P., Maloy, S. & Minor, A.

*In situ*nanocompression testing of irradiated copper.*Nat. Mater.***10**, 608–613, doi:10.1038/nmat3055 (2011). - 7.
Yuan, X.

*et al*. Effect of plastic deformation on the Σ2 grain boundary plane distribution in WC–Co cemented carbides.*Int. J. Refract. Met. Hard Mater.***47**, 38–43, doi:10.1016/j.ijrmhm.2014.05.012 (2014). - 8.
Klassen-Neklyudova, M. V.

*Mechanical twinning of crystals*. (Springer Science & Business Media, 2012). - 9.
Liu, D.

*et al*. Anomalous plasticity in the cyclic torsion of micron scale metallic wires.*Phys. Rev. Lett.***110**, 244301, doi:10.1103/PhysRevLett.110.244301 (2013). - 10.
Fu, S., Wang, L., Chen, G., Yu, D. & Chen, X. A tension-torsional fatigue testing apparatus for micro-scale components.

*Rev. Sci. Instrum.***87**, 015111, doi:10.1063/1.4939856 (2016). - 11.
Walter, M. & Kraft, O. A new method to measure torsion moments on small-scaled specimens.

*Rev. Sci. Instrum.***82**, 035109, doi:10.1063/1.3557824 (2011). - 12.
Liu, D., He, Y., Hu, P. & Ding, H. Characterizing Torsional Properties of Microwires Using an Automated Torsion Balance.

*Experimental Mechanics*, 1–15 (2016). - 13.
Liu, D., Peng, K. & He, Y. Direct measurement of torsional properties of single fibers.

*Meas. Sci. Technol.***27**, 115017, doi:10.1088/0957-0233/27/11/115017 (2016). - 14.
Cohen-Karni, T., Segev, L., Srur-Lavi, O., Cohen, S. R. & Joselevich, E. Torsional electromechanical quantum oscillations in carbon nanotubes.

*Nat Nano***1**, 36–41, doi:10.1038/nnano.2006.57 (2006). - 15.
Nagapriya, K.

*et al*. Torsional stick-slip behavior in WS 2 nanotubes.*Phys. Rev. Lett.***101**, 195501, doi:10.1103/PhysRevLett.101.195501 (2008). - 16.
Garel, J.

*et al*. Ultrahigh torsional stiffness and strength of boron nitride nanotubes.*Nano Lett.***12**, 6347–6352, doi:10.1021/nl303601d (2012). - 17.
Dong, L., Tao, X., Zhang, L., Zhang, X. & Nelson, B. J. Nanorobotic spot welding: controlled metal deposition with attogram precision from copper-filled carbon nanotubes.

*Nano Lett.***7**, 58–63, doi:10.1021/nl061980+ (2007). - 18.
Li, G., Xi, N., Yu, M. & Fung, W.-K. Development of augmented reality system for AFM-based nanomanipulation.

*IEEE/ASME Transactions on mechatronics***9**, 358–365, doi:10.1109/TMECH.2004.828651 (2004). - 19.
Fukuda, T., Arai, F. & Dong, L. Assembly of nanodevices with carbon nanotubes through nanorobotic manipulations.

*Proceedings of the IEEE***91**, 1803–1818, doi:10.1109/JPROC.2003.818334 (2003). - 20.
Yajing Shen, Masahiro Nakajima, Zhenhai Zhang & Toshio Fukuda. Dynamic Force Characterization Microscopy Based on Integrated Nanorobotic AFM and SEM System for Detachment Process Study.

*IEEE-ASME Transactions on Mechatronics***20**(6), 3009–3017, doi:10.1109/TMECH.2015.2413779 (2015). - 21.
Shi, C.

*et al*. Recent advances in nanorobotic manipulation inside scanning electron microscopes.*Microsystems & Nanoengineering***2**(2016). - 22.
Zimmermann, S.

*et al*. Automated mechanical characterization of 2-D materials using SEM based visual servoing.*International Journal of Optomechatronics***7**, 283–295, doi:10.1080/15599612.2013.879501 (2013). - 23.
Shang, W., Lu, H., Wan, W., Fukuda, T. & Shen, Y. Vision-based Nano Robotic System for High-throughput Non-embedded Cell Cutting.

*Scientific reports***6**, doi:10.1038/srep22534 (2016). - 24.
Chen, N. & Bhushan, B. Morphological, nanomechanical and cellular structural characterization of human hair and conditioner distribution using torsional resonance mode with an atomic force microscope.

*Journal of microscopy***220**, 96–112, doi:10.1111/jmi.2005.220.issue-2 (2005). - 25.
Lu, W.-Y. & Song, B. Quasi-static torsion characterization of micro-diameter copper wires.

*Experimental Mechanics***51**, 729–737, doi:10.1007/s11340-010-9377-8 (2011). - 26.
Song, B. & Lu, W.-Y. An improved experimental technique to characterize micro-diameter copper wires in torsion.

*Experimental Mechanics***55**, 999–1004, doi:10.1007/s11340-014-9978-8 (2015). - 27.
Behlow, H.

*et al*. Direct measurement of shear properties of microfibers.*Rev. Sci. Instrum.***85**, 095118, doi:10.1063/1.4895679 (2014). - 28.
Yajing Shen, Wenfeng Wan, Haojian Lu, Toshio Fukuda & Wanfeng Shang. Automatic Sample Alignment Under Microscopy for 360° Imaging Based on the Nanorobotic Manipulation System.

*IEEE Transactions on Robotics***33**(1), 220–226, doi:10.1109/TRO.2016.2623337 (2017). - 29.
Suzuki, S. Topological structural analysis of digitized binary images by border following.

*Computer Vision, Graphics, and Image Processing***30**, 32–46, doi:10.1016/0734-189X(85)90016-7 (1985).

## Acknowledgements

This work was supported in part by the National Natural Science Foundation of China 61403323, ShenZhen (China) Basic Research Project No. JCYJ20160329150236426, Natural Science Foundation of Shaanxi Province of China (2016JQ5062) and the RGC General Research Fund 11278716.

## Author information

### Affiliations

### Contributions

H. Lu performed the twisting test experiment and analyzed the data. W. Shang proposed the alignment strategy and control system. X. Wei, Z. Yang and T. Fukuda provided valuable advices in robot developed. Y. Shen conceived the project, developed the nanorobotic system, designed the experiment and analyzed the data. H. Lu and W. Shang contributed equally to this work. All the authors drafted the manuscript.

### Corresponding author

## Ethics declarations

### Competing Interests

The authors declare that they have no competing interests.

## Additional information

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Electronic supplementary material

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Lu, H., Shang, W., Wei, X. *et al.* Nanorobotic System iTRo for Controllable 1D Micro/nano Material Twisting Test.
*Sci Rep* **7, **3077 (2017). https://doi.org/10.1038/s41598-017-03228-4

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.