Home | About | Journals | Submit | Contact Us | Français |

**|**Sensors (Basel)**|**v.17(2); 2017 February**|**PMC5335975

Formats

Article sections

- Abstract
- 1. Introduction
- 2. Nano-Positioning Stage
- 3. Capacitive Comb Sensor
- 4. Experiments
- 5. Conclusions
- References

Authors

Related links

Sensors (Basel). 2017 February; 17(2): 257.

Published online 2017 January 28. doi: 10.3390/s17020257

PMCID: PMC5335975

Stefano Mariani, Academic Editor

Jiangsu Provincial Key Laboratory of Advanced Robotics & Collaborative Innovation Center of Suzhou Nano Science and Technology, Soochow University, Suzhou 215123, China; Email: nc.ude.adus@tnehc (T.C.); Email: nc.ude.adus.uts@30092245102 (Y.W.); Email: nc.ude.adus@gnawbp (P.W.); Email: nc.ude.tih@nusnl (L.S.)

Received 2016 November 25; Accepted 2017 January 24.

Copyright © 2017 by the authors.

Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).

This article has been cited by other articles in PMC.

Nano-positioning technology has been widely used in many fields, such as microelectronics, optical engineering, and micro manufacturing. This paper presents a one-dimensional (1D) nano-positioning system, adopting a piezoelectric ceramic (PZT) actuator and a multi-objective topological optimal structure. The combination of a nano-positioning stage and a feedback capacitive comb sensor has been achieved. In order to obtain better performance, a wedge-shaped structure is used to apply the precise pre-tension for the piezoelectric ceramics. Through finite element analysis and experimental verification, better static performance and smaller kinetic coupling are achieved. The output displacement of the system achieves a long-stroke of up to 14.7 μm and high-resolution of less than 3 nm. It provides a flexible and efficient way in the design and optimization of the nano-positioning system.

With continuous developments of the integrated circuit (IC) manufacturing process, ultra-precision machining, and bio-medical engineering, long-stroke and precise nano-positioning technology have been one of the hottest research areas in precision engineering [1,2]. Recent research has focused on the design and manufacturing of nano-positioning stages and high-precision sensors, as well as precise control methods [3,4,5,6]. Compared with traditional bulk mechanical stages, the nano-positioning stages have shorter product lifecycles and higher performance indicators to satisfy the development requirements of microelectronics [7]. The flexible and efficient design methodologies are urgent for meeting the high demands of nano-positioning stages.

Topology optimization, as a kind of conceptual design methodology, could be a choice of priority [8]. Different from the traditional optimizations of sizes and shapes [9], the topology optimization can achieve an optimal structural layout for different conditions without prior information. Meanwhile, the objective function of topology optimization, as design variables, can be designed to achieve the diversification of product design according to different requirements. Related to the initial design and post processing of the nano-positioning stage, the topology optimization can enhance and improve the stage performance significantly [10,11]. Ansola et al. [10] applied an evolutionary structural optimization (ESO) procedure for compliant mechanism design. They introduced the elastic strain energy of the mechanism within the mechanical advantage expressions used as objective functions, in order to improve convergence in the optimization. Chen et al. [11] applied the topology optimization in the size optimization of the stage and enhanced the static and dynamic characteristics significantly. In order to overcome the defects of traditional testing sensors, such as being oversized and low accuracy, high-precision nano-displacement sensors have emerged. Among these, the capacitive sensors have become the first choice of the testing devices, due to their high precision, good sensitivity, and high stability [12,13]. To overcome the drawbacks of limited range of capacitive testing, the capacitive sensors are generally made into offset constant-tooth comb-finger structures to satisfy the range and accuracy demands of nano-positioning [14,15].

The proposed nanopositioners above all have a good performance providing tens of micrometers range of displacement with nanometer resolution and a high bandwidth. However, the structures of the nanopositioners are more complicated, which may cause the fabrication difficulties and high costs. In addition, regarding the comb drive capacitive sensors, researchers are tending to design nanopositioners with a monolithic structure to form a closed-loop control. Here in this paper, we have adopted the topology optimization approach to achieve a wide displacement range and a fine resolution of the nano-positioning stage. In addition, the comb drive capacitive sensors and the nanopositiong are separately designed and connected with each other using a simple connection. The overall structure is relatively simple, and manufacturing cost is very low. Such an idea has proved to be a bold attempt and has also achieved a good result. The goal of this project is to develop an integrated model of micro displacement precision positioning system. The requirements of this system can achieve 10 μm travel, repeat positioning accuracy with 10 nm, and displacement resolution with 5 nm. This paper presents the displacement positioning system of a 1D stage at the nanometre scale by using a microelectromechanical system (MEMS)-based displacement sensor. The 1D nano-positioning stage is optimized by multi-objective topology optimization. A wedge pretension has been applied to verify the accuracy of the stage. An on-chip sensor of comb-finger capacitor is micro-fabricated with a high displacement sensing resolution and is constructed with the stage as a testing device to measure the displacement change of the stage. Experimental results of the travel distance and motion precision of the displacement-controlled manipulation are presented. The performance of the topology stage is superior to the traditional one. Through the combination with the comb sensor, the nano-positioning system can realize the installation of miniaturization and integration, competent for the requirement of precision engineering.

The schematic drawing of the proposed nano-positioning system is shown in Figure 1. It is comprised of a piezoelectric ceramic (PZT) actuator, a lever amplification mechanism, a load stage, and a capacitive comb sensor. The fabrication technique of the stage is wire-electrode cutting. Four multi-leaf springs are designed and machined according to the stiffness requirement. The lever connects the load stage and the PZT actuator by hinges. The topology structure on the lever includes three notches. A preloading screw provides a preload compressive force on the wedge structure for maintaining the feeding accuracy of the PZT actuator. The working principle of the nano-positioning stage is to drive flexible hinges to achieve nanometer motion through the inverse piezoelectric effect of PZT. The whole stage is made of duralumin (Al-Mg alloy, LY11). The size of the comb sensor is 10 mm × 10 mm × 1 mm, including 6 × 70 pairs of differential combs, and its maximum travel distance can achieve 15 μm. The size of the stage, decided by the PZT ceramics, is currently 60 mm × 80 mm × 20 mm, and it can be further reduced to 30 mm × 30 mm × 10 mm. The design objective is to obtain a maximum output displacement of up to 10 μm as well as positioning accuracy of smaller than 5 nm, with high stiffness, good dynamic property and system stability.

In traditional design, both the kinematics and structure must be considered, while the structure design is normally ignored due to the optimal rule [16]. We use mutual energy concept [17] to describe the requirements as follows:

$$\begin{array}{cc}{L}^{k}({u}^{l})=L({t}^{k},{u}^{l})={\displaystyle \underset{{\mathsf{\Gamma}}_{tk}}{\int}t(k){u}^{l}}d\mathsf{\Gamma}& {u}^{l}\end{array}\in V,$$

(1)

$$a({u}^{k},{v}^{l})={\displaystyle \underset{D}{\int}\epsilon {({v}^{l})}^{T}E\epsilon ({u}^{k})}dD,$$

(2)

$$\begin{array}{cc}a({u}^{k},{v}^{l})={L}^{k}({v}^{l})& {u}^{k}\in V,\forall {v}^{l}\end{array}\in V,$$

(3)

where *t*(*k*) is the force at the point *k*, and *u ^{l}* is the displacement at the point

The requirement of kinematics, in order to maximize mutual mean compliance *L ^{q}*(

On the basis of the virtual work principle, it is set as *t*^{3} = *t*^{1} and *t*^{4} = −*t*^{2}. Thus, the design can be described as a multi-objective problem, consisting of four objective functions. In the case, the mutual mean compliance *L*(*t*^{2},*u*^{1}) must be maximized; the mean compliances *L*(*t*^{3},*u*^{1}) and *L*(*t*^{4},*u*^{2}) must be minimized; and the mutual mean compliance *L*(*t*^{5},*u*^{1}) must be minimized. To find the incorporation optimal solution, the multi-objective function is formulated as:

$$\mathrm{max}{f}_{1}=\frac{L({t}^{2},{u}^{1})}{{({w}_{s}L{({t}^{3},{u}^{1})}^{2}+(1-{w}_{s})L{({t}^{4},{u}^{2})}^{2}+{w}_{M}L{({t}^{5},{u}^{1})}^{2})}^{\frac{1}{2}}}.$$

(4)

To be more general:

$$\begin{array}{l}\mathrm{max}{f}_{2}=W\mathrm{log}({L}^{2}({u}^{1}))\\ -\frac{1}{2}(1-W)\mathrm{log}({w}_{s}{L}^{3}{({u}^{3})}^{2}+(1-{w}_{s}){L}^{4}{({u}^{4})}^{2}+{w}_{M}{L}^{5}{({u}^{1})}^{2})\end{array},$$

(5)

where *W*, *w _{s}* and

To get a specified optimization structure, the effect of weighting coefficients *W*, *w _{s}*, and

As the parameters of the multi-objective function are determined, the convergence history of the topology optimization process can be drawn through the finite element method. Figure 3a illustrates the convergence of the objective function (in blue) and the relationship between the objective function and total volume (in green). It is noted from Figure 3b that the mutual compliance *L*^{2}(*u*^{1}) increases as the iteration continues, while the compliances *L*^{3}(*u*^{3}) and *L*^{4}(*u*^{4}) decrease.

(**a**) History of the objective function and total volume; and (**b**) History of mutual compliance and self-compliances.

Finally, the preliminary configuration of the nano-positioning stage is obtained. However, due to the chessboard problem and partial polarization [21], the preliminary configuration produces some “transitional element”. Meanwhile, due to the complex and ambiguous distribution, it is difficult to pick up the element from the topology configuration. Thus, the image-based design method [22] is employed to draw the topology structure. At first, the characteristic function, i.e., pixels, is defined instead of the element density. Then, the distribution of topology optimal configuration for each element is represented by the integer 0 to 27 scale value. By setting an appropriate threshold, the integer scale values are converged, through the method of density penalty [23], to “0” and “1”. The scale values that are larger than the threshold form the real structure, and the ones that are smaller than the threshold form the hole. Sometimes, some fitting curves are necessary in order to make the topology configuration behave well in terms of manufacturability.

In order to verify the superiority of the multi-objective topology optimization, a traditional MPT-1JR micro-stage (designed by Jiangsu Huibo Robotics Technology Co., Ltd., Suzhou, China) with the same size optimization has been applied as well. Both the static and dynamic characteristics are compared.

- (1)Static analysis: On the basis of finite element modelling, static analysis is conducted and the input and output displacements of the hinges at points B and A are studied as shown in Figure 4. The input and output displacements of the points in the core of a circular hinges stand for the mean displacements of the nano-positioning stage. Through the multi-objective topology optimization, the maximum output displacement, amplification ratio, as well as the displacement damage in the PZT are calculated and summarized in Table 2. Compared with the traditional size optimization, it is found that the maximum output displacement and the amplify ratio of the topology design increase, while the input remains approximately the same level. As the goal displacement feed is obtained, the topology design can reduce the damage in PZT because of the reduction of reaction force.Static analysis of the output displacements at points (
**A**) and (**B**) of the hinge with topology design.Besides the kinematics requirement, the structure requirement is considered in the multi-objective topology. The two mean compliances*L*(*t*^{3},*u*^{1}) and*L*(*t*^{4},*u*^{2}) at the input and output points remain sufficiently stiff. Based on the ANSYS simulation results as shown in Table 3, the reduction of the stress values can be clearly observed, that is, the stiffness of the input and output gain the increase. Meanwhile, the stiffness changes of the output and pivot are similar due to the mechanism principles. The vertical displacement of the goal direction is tested. It is found that a small coupling ratio is suitable for the goal requirement and propitious to expand the multi-dimensional stages with reduced cross coupling. - (2)Dynamic characteristics: The stability is dominated by the first resonance mode of the body structure, which is normally a bending mode. In order to obtain high stability, the nano-positioning stage must be designed to raise the frequency of the first resonance mode as high as possible. The mode analysis is performed as shown in Figure 5. The resonance frequency is found to increase from 391.7 to 470.3 Hz through the multi-objective topology optimization. The resonance response of the nano-positioning stage is studied. When the PZT starts or stops operating, because of the characteristics of the PZT ceramic itself, such as creep or hysteresis, the system will result in some frequency disturbance. For a traditional MPT-1JR micro-stage, the frequency response in the PZT ceramics is higher than the one in micro-stage, and oscillation will come up to make the entire system unstable. Through topology optimization, the frequency response can be effectively improved to avoid the emergence of such disturbance.

To overcome the hysteresis and nonlinearity defects of the PZT ceramic and improve the positioning accuracy of the overall system, it is necessary to integrate a displacement sensor with high precision and miniaturization to form a feedback system. As shown in Figure 6a, the capacitive comb sensor consists of fixed comb electrodes, movable comb electrodes, S-type beams and flexible cantilevers. The offset constant-tooth structure with the actuation displacement in the *y-*direction in Figure 6b has been applied in the comb sensor to compensate the nonlinearity and hysteresis effect of piezoelectric ceramic and improve the positioning accuracy. In the figure, *F* is the fixed beam, *M* is the movable beam, *t* is the thickness of comb electrodes, m is the width of electrodes, d is the gap between *F* and *M* electrodes, *L*_{0} is the overlap distance of the electrodes, *L _{m}* and

Schematic drawing of (**a**) the comb sensor and (**b**) the offset constant-tooth structure with readout circuitry

The capacitive comb sensor is microfabricated by a series of MEMS processes as shown in Figure 7.

Process sequence of the micro comb sensor. (**a**) Oxidation and lithography; (**b**) Dry etching; (**c**) Bonding of silicon and glass; (**d**) Sputtering of aluminum; (**e**) Lithography of structure; (**f**) Etched using DRIE.

(a) It is started from a 4-inch, N-type, (100) orientation, double polished silicon (Si) wafer. The back release windows are patterned on the backside of the wafer; (b) a backside deep reactive ion etching (DRIE) process is performed to form a depth of 250 μm; (c) then, a glass wafer is bonded onto the Si wafer to form the support base; (d) the aluminum layer of 1-μm-thick is formed by vacuum evaporation; (e) the aluminum wires and bonding pads are formed by etching processes and then the structure of sensor is patterned on the wafer; and, (f) finally, the device layer including the comb fingers and the flexure beams are etched using DRIE.

With the combination of surface and bulk micromachining processes, the comb sensor demonstrates good detective performance of precision scale and is easy to integrate with the nano-positioning stage. The photo and SEM images of the sensor are shown in Figure 8.

Experiments have been conducted to verify the simulation results. The testing system as shown in Figure 9 includes a Keyence LK-G CCD Laser Displacement Sensor (Keyence, Osaka, Japan), a PTBS200 piezoelectric PZT ceramic (Boshi, Harbin, China), and the fabricated nano-positioning stage. A dual-frequency laser interferometer ZLM700 (JENAer, Jena, Germany) is used to test the resolution of the stage. The stage driven by the piezoelectric ceramic PTBS200/5 × 5/5 is manufactured with a duralumin alloy (LY11), which has a material density of 2.7 g/mm^{3}, Young’s modulus 89 GPa, Poisson ratio of 0.33 and yielding strength of 505 MPa. The size of the stage is 80 mm × 60 mm × 20 mm.

The calibration curve of the comb sensor is provided in Figure 10 according to the relationship between physical change and electrode spacing change. It can be seen that the measurement range of the comb sensor should between −5 μm and 5 μm for a better linearity and the linearity error of the sensor is 3%. The hysteresis curves of the input and output displacements of the stage at points B and A are shown in Figure 11a. As can be seen, the output displacement is up to 14.7 μm at 150 V. On the basis of the output displacement curve, the resolution of the nano-positioning stage is results in being less than 3 nm for an input voltage increment of 1 mV by the Laser Displacement Sensor. To compare the change of the output characteristics and amplifier ratio of the topology optimized stage, a high-order curve fitting has been done as shown in Figure 11b. Referring to Labels (7) and (8), the amplification ratio is about 3.99, which is essentially consistent with the finite element simulation result. It shows that topology optimization can improve the static characteristics of the nano-positioning stage.

(**a**) The input and output displacements of the stage; (**b**) high-order displacement fitting curves of the input and output displacements.

For an applied voltage of 150 V, the displacements along the *x*-direction is 14.7 μm, while along the *y*-direction, it is 0.1 μm. Thus, the cross coupling ratio is 0.6%, which is mainly due to the manufacturing issues, and the actual parts of the topology in the corner and the inflection point have maching errors. In order to control the precision of the nano-positioning stage, the output displacement of the stage is calibrated by the integrated comb sensor before the next closed loop. Figure 12 shows that the comb sensor is suitable for integration and miniaturization in the design of nano-positioning systems.

To test the dynamic characteristics of the system, the PZT power device provides the full value of the step voltage. Then, a laser interferometer is used to measure the corresponding time, as shown in Figure 13. We can see that the response time of rising edge is 2 ms. Therefore, the natural frequency of the system is about 500 Hz.

In order to verify the repetitive positioning accuracy and the stability time of system designed, the displacement response test and the repeat positioning tests of fixed points are carried out, respectively. With the control system commands, the microstage is located at the 10 μm position. An oscilloscope is used to test the response curve of control. The results show that the whole control process of rising is stable. There is no disturbance and overshoot.

With the control system, the microstage is positioned at the centre position (about 7.4 μm), and 10 sets of data are tested. We can see that the repeat positioning accuracy is 10 nm as shown in Table 5. It fully meets the design requirements.

This paper presents a novel nano-positioning system with high-precision, large-stroke and high-integration. The multi-object topology optimization is applied to increase the amplification ratio, reduce the stiffness distribution and improve the frequency property of the nano-positioning stage. Based on the static test, the stiffness at the input and output points is increased and the coupling ratio is reduced. Through the reduction of the whole stiffness, resonant frequency of the nano-positioning system is increased, avoiding the resonance induced oscillation and improving the stability of the system. Combined with the offset constant-tooth comb sensor, the nano-positioning system realizes the installation of integration and miniaturization.

This work is supported by the National Natural Science Foundation of China (No. 61673287, No. 61433010) and the National High-Tech Research and Development Program of China (No. 2015AA042601).

Author Contributions

All of the authors have equally contributed to this work, either in the implementations or in the writing of this article. T.C., L.S. and H.L. managed this work; T.C., Z.Y. and L.S. designed the nano-positioning stage and capacitive comb sensors; Z.Y. and P.W. analyzed and optimized the nano-positioning system; H.L., P.W. and Y.W. conducted fabrication of the capacitive comb sensor; Y.W., Z.Y. and P.W. implemented the experiment; and T.C., H.L. and Y.W. wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

1. Sun L., Wang J., Rong W., Li X., Bao H. A silicon integrated micro nano-positioning XY-stage for nano-manipulation. J. Micromech. Microeng. 2008;18 doi: 10.1088/0960-1317/18/12/125004. [Cross Ref]

2. Hyde H.C., Sandtner W., Vargas E., Dagcan A.T., Robertson J.L., Roux B., Correa A.M., Bezanilla F. Nano-positioning system for structural analysis of functional homomeric proteins in multiple conformations. Structure. 2012;20:1629–1640. doi: 10.1016/j.str.2012.08.022. [PMC free article] [PubMed] [Cross Ref]

3. Zhang D.Y., Ono T. Piezoactuator-integrated monolithic micro-stage with six degrees of freedom. Sens. Actuators A Phys. 2005;122:301–306. doi: 10.1016/j.sna.2005.03.076. [Cross Ref]

4. Zhang Y., Keum H., Park K., Bashir R., Kim S. Micro-Masonry of MEMS Sensors and Actuators. J. Microelectromech. Syst. 2014;23:308–314. doi: 10.1109/JMEMS.2013.2273439. [Cross Ref]

5. Xu Q.S. Design and Development of a Flexure-Based Dual-Stage Nanopositioning System with Minimum Interference Behavior. IEEE Trans. Autom. Sci. Eng. 2012;9:554–563. doi: 10.1109/TASE.2012.2198918. [Cross Ref]

6. Wadikhaye S.P., Yong Y.K., Bhikkaji B., Moheimani S.O.R. Control of a piezoelectrically actuated high-speed serial-kinematic AFM nanopositioner. Smart Mater. Struct. 2014;23 doi: 10.1088/0964-1726/23/2/025030. [Cross Ref]

7. Yong Y.K., Moheimani S.O.R., Kenton B.J., Leang K.K. Invited Review Article: High-speed flexure-guided nanopositioning: Mechanical design and control issues. Rev. Sci. Instrum. 2012;83 doi: 10.1063/1.4765048. [PubMed] [Cross Ref]

8. Jain C., Saxena A. An Improved Material-Mask Overlay Strategy for Topology Optimization of Structures and Compliant Mechanisms. J. Mech. Des. 2010;132 doi: 10.1115/1.4001530. [Cross Ref]

9. Perezzan J.C., Hernandez S. Size and shape optimization of long portal bridges. Trans. Wessex Inst. 1999;5:287–300.

10. Ansola R., Veguería E. A simple evolutionary topology optimization procedure for compliant mechanism design. Finiteelem. Anal. Des. 2007;44:53–62. doi: 10.1016/j.finel.2007.09.002. [Cross Ref]

11. Jong C., Chen K.Y. An Integrated Design Optimization for Monolithic Mechanical Amplifier in PZT Nano-Positioning Stage. Tamkang J. Sci. Eng. 2010;13:165–172.

12. Xu Q.S. New Flexure Parallel-Kinematic Micropositioning System with Large Workspace. IEEE Trans. Robot. 2012;28:478–491. doi: 10.1109/TRO.2011.2173853. [Cross Ref]

13. Oda K., Terao K., Suzuki T., Takao H., Ishimaru I., Oohira F. Vertical comb-drive MEMS mirror with sensing function for phase-shift device. Sens. Actuators A Phys. 2015;181:61–69. doi: 10.1016/j.sna.2012.04.007. [Cross Ref]

14. Maroufi M., Bazaei A., Moheimani S.O.R. A High-Bandwidth MEMS Nanopositioner for On-Chip AFM: Design, Characterization, and Control. IEEE Trans. Control Syst. Technol. 2015;23:504–512. doi: 10.1109/TCST.2014.2345098. [Cross Ref]

15. Xu Q. Design, Fabrication, and Testing of an MEMS Microgripper with Dual-Axis Force Sensor. IEEE Sens. J. 2015;15:6017–6026. doi: 10.1109/JSEN.2015.2453013. [Cross Ref]

16. Jang I.G., Kim I.Y. Computational study of Wolff’s law with trabecular architecture in the human proximal femur using topology optimization. J. Biomech. 2008;41:2353–2361. doi: 10.1016/j.jbiomech.2008.05.037. [PubMed] [Cross Ref]

17. David G.P. Sensitivity analysis of linear elastic systems using domain parameterization and a mixed mutual energy principle. Comput. Methods Appl. Mech. Eng. 1989;77:31–59.

18. Teper R. On the continuity of the concave integral. Fuzzy Sets Syst. 2009;160:1318–1326. doi: 10.1016/j.fss.2008.07.001. [Cross Ref]

19. Kim Y.S., Yoo J.M., Yang S.H., Choi Y.M., Dagalakis N.G., Gupta S.K. Design, fabrication and testing of a serial kinematic MEMS XY stage for multifinger manipulation. J. Micromech. Microeng. 2012;22:85029–85038. doi: 10.1088/0960-1317/22/8/085029. [Cross Ref]

20. Yong Y.K., Liu K., Moheimani S.O.R. Reducing Cross-Coupling in a Compliant XY Nanopositioner for Fast and Accurate Raster Scanning. IEEE Trans. Control Syst. Technol. 2010;18:1172–1179. doi: 10.1109/TCST.2009.2033201. [Cross Ref]

21. Rueda A., Acosta O., Couprie M., Bourgeat P., Fripp J., Dowson N., Romero E., Salvado O. Topology-corrected segmentation and local intensity estimates for improved partial volume classification of brain cortex in MRI. J. Neurosci. Methods. 2010;188:305–315. doi: 10.1016/j.jneumeth.2010.02.020. [PubMed] [Cross Ref]

22. Eisert P., Rurainsky J. Geometry-assisted image-based rendering for facial analysis and synthesis. Signal Process. Image Commun. 2006;21:493–505. doi: 10.1016/j.image.2006.03.003. [Cross Ref]

23. Choi J.S., Yamada T., Izui K., Nishiwaki S., Yoo J. Topology optimization using a reaction–diffusion equation. Comput. Methods Appl. Mech. Eng. 2011;200:2407–2420. doi: 10.1016/j.cma.2011.04.013. [Cross Ref]

Articles from Sensors (Basel, Switzerland) are provided here courtesy of **Multidisciplinary Digital Publishing Institute (MDPI)**

PubMed Central Canada is a service of the Canadian Institutes of Health Research (CIHR) working in partnership with the National Research Council's national science library in cooperation with the National Center for Biotechnology Information at the U.S. National Library of Medicine(NCBI/NLM). It includes content provided to the PubMed Central International archive by participating publishers. |