The present invention relates to a technology for measuring a size distribution of a particle contained in a liquid sample.
In recent years, a development target of a drug moves from a low molecule drug to a bio drug. Since the bio drug is a polymer, the bio drugs are easy to be aggregated, and when the bio drugs are aggregated, toxicity may be generated. For example, the U.S. Food and Drug Administration and the like try to strengthen concentration control regulations of an aggregate. Therefore, there is a need for a technology for quantitatively measuring a size distribution of a desired density with respect to an aggregate in a submicron region of 0.1 to 1 um.
JP-A-2015-083922 (PTL 1) describes a technology for detecting a particle by using an optical measurement. An object of JP-A-2015-083922 (PTL 1) is that “in a single particle detection technology by a scanning molecule counting method which individually detects a single particle by using optical measurement with a confocal microscope or a multiphoton microscope, the single particle can be detected for each type in a sample solution containing a plurality of types of single particles that do not emit light in a specific wavelength band”, and JP-A-2015-083922 (PTL 1) describes a technology that “a technology for detecting a single particle in a sample solution according to the present invention detects light in a specific wavelength band from a light detection region to generate time-series light intensity data while moving a position of a photodetection region of a microscope in a sample solution containing a plurality of types of single particles that do not emit light in different wavelength bands, and individually detects a decrease in light intensity generated when the single particle that does not emit light in the specific wavelength band enters the light detection region in the time-series light intensity data as a signal indicating the existence of each of the single particles” (refer to abstract).
An object of JP-A-2017-102032 (PTL 2) is to “provide a technology capable of reducing a measurement error when measuring a specimen by using light.”, and JP-A-2017-102032 (PTL 2) describes a technology that “an optical measurement method according to the present invention acquires correspondence relationship data that describes a correspondence relationship between an intensity of reflected light when the specimen is irradiated with light and a size of the specimen, and acquires the size of the specimen by using the correspondence relationship data and the intensity of the reflected light. The optical measurement method according to the present invention corrects an inclination of a container by subtracting a component caused by the inclination of the container of the specimen from a detection signal indicating the intensity of the reflected light when the specimen is irradiated with the light.” (refer to abstract).
PTL 1: JP-A-2015-083922
PTL 2: JP-A-2017-102032
In order to obtain a size distribution of an aggregate, it is required to measure a size of a particle. In order to measure the size of the particle by optical measurement, it is considered that it is required not only to specify a plane position of the particle, but also to obtain an optical signal along a depth direction of the particle. The reason is that the particles exist at various depth positions.
In JP-A-2015-083922 (PTL 1), a position of a photodetection region is caused to move faster than Brownian motion of the particle, thereby following Brownian motion of the particle (paragraph 0016 of the same patent literature). However, in the same patent literature, since measurement light is not scanned in the depth direction of the particle (a direction along an optical axis), it is considered that the measurement of the particle size is difficult to be performed even though presence of the particle is detected.
In JP-A-2017-102032 (PTL 2), the particle size is calculated by referring to correspondence relationship data using reflected light intensity. However, in the same patent literature, it is assumed that reflected light from a stationary particle is used, such that it is considered that calculating the size of the particle performing Brownian motion in the depth direction is not assumed.
The present invention has been made in consideration of the above-described circumstances, and an object thereof is to provide an optical measurement technology capable of quantitatively measuring a size distribution of a particle that performs Brownian motion in a sample.
A size distribution measurement device according to the present invention measures a reflected light intensity while scanning a focal point position along an optical axis direction of measurement light, and calculates a size distribution of the particle according to the highest reflected light intensity of the measured reflected light intensities.
According to a size distribution measurement device according to the present invention, it is possible to quantitatively measure a size distribution of a particle that performs Brownian motion in a three-dimensional manner.
In the following, in order to facilitate the understanding of the present invention, a basic principle of measuring a size of a particle will be described first. Next, a specific configuration example of a size distribution measurement device according to embodiments of the present invention will be described.
However, depending on the correspondence relationship between the particle size and the focal point position, the reflected light intensities from particles having different sizes may be the same (refer to a dotted line portion in
A design parameter of an optical system of the size distribution measurement device will be described. In order to measure a particle diameter based upon the reflected light intensity, it is desirable that a spot size of the measurement light is slightly larger than an upper limit size of the particle to be measured. Therefore, on the assumption that a wavelength of the measurement light is defined as X and a numerical aperture of an objective lens is defined as NA, an optical spot diameter Rparticle can be represented by the following formula 1.
In order to detect the reflected light in a state where the focal point is almost positioned on the particle, a z pitch (Δd) is required to be smaller than a resolution Δz in a z direction. That is, the following formula 2 is required to be satisfied. Δz can be defined as a distance at which a reflected signal intensity decreases from a maximum value to a certain threshold value (for example, half of the maximum value) when the focal point position is scanned in the z direction.
[Formula 2]
Δd<Δz (2)
An image acquisition speed can be defined as the number of image acquisitions per second. The number thereof is referred to as a frame rate below. Since the particle in the liquid to be measured performs Brownian motion, in order to obtain a correct particle diameter distribution, it is required to acquire an image at a frame rate that can follow Brownian motion. In the following, the required frame rate will be described by using a formula. Here, for simplifying the description, only Brownian motion in the z direction is considered. An average movement amount L of the particle performing Brownian motion during a time t is represented by the following formula 3.
[Formula 3]
L=√{square root over (2 Dt)} (3)
D in formula 3 is a diffusion constant of the particle, and is represented by the following formula 4 using a Boltzmann constant k, an absolute temperature T, a viscosity μ of a fluid, and a particle diameter a.
In consideration of the required frame rate, it is desirable to interpret that a measurement z position does not move by Δd at regular time intervals, a z position of the particle moves by Δd independently of a component based upon Brownian motion, and the measurement z position is fixed at an origin point. Hereinafter, the embodiment will be interpreted and described as such. Under this interpretation, an average movement amount Lm of the particle disposed at the origin point at time 0 up to m-th measurement (time tm) is represented by the following formula 5. Δt is a reciprocal number of the frame rate (Δt=1/FR), and tm=mΔt.
On the assumption that z resolution is defined as Δz, it is considered that average time Tave at which the particle stays within a z resolution range approximately satisfies the following formula 6.
The following formula 7 is obtained by solving formula 6 with respect to Tave. For simplification of the formula, α=Δz/Δd is set.
In order to detect the reflected light from the particle in a state where the focal point is almost positioned on the particle, an image is required to be acquired a plurality of times while the particle stays within the z resolution range. Such a condition can be represented by the following formula 8. n is a natural number of 1 or more, and is the average number of times the particle stays within the z resolution range.
[Formula 8]
Tave>n Δt (8)
It can be said that formula 8 has a condition that the reflected light from the particle can be continuously detected n times or more on average. Formula 9 can be obtained by substituting formula 7 into formula 8 and rearranging it.
The following formula 10 is obtained by rewriting formula 9 using a relationship of FR=1/Δt.
Since the average number of stays n is considered to be approximately proportional to Δz/Δd, n is represented by the following formula 11 using a proportionality coefficient β of 1 or less.
The following formula 12 is obtained by further rearranging formula 10 using the formula 11. Here, γ=2β/(1−β)2.
Formula 12 is a condition that the frame rate is required to be satisfied in order to correctly measure the particle diameter distribution which is predicted based upon theoretical consideration. A fact that a right side is proportional to a diffusion constant D indicates that as the movement of the particle becomes faster, a high frame rate is required to follow the movement thereof. A fact that the right side is inversely proportional to the z resolution Δz indicates that as the z resolution becomes higher, an amount of deviation from the focal point position of the particle that can be tolerated in detecting the reflected light becomes smaller, such that a higher frame rate is required. A fact that the right side is inversely proportional to the z pitch Δd indicates that as the z pitch becomes finer, a measurement speed in the z direction decreases, such that a higher frame is required. Since it is difficult to predict a value of γ only by theoretical consideration, it is desirable to determine the value thereof based upon simulation or experiment.
According to formula 12, as Δz and Δd become larger, the required frame rate becomes lower. On the other hand, when Δz becomes too large, a probability that a plurality of particles exist within the z resolution range becomes high, and as a result, it becomes difficult to measure a sample having high particle concentration. The z resolution should be determined in consideration of a particle concentration range of a target particle. As described above, there is a restriction that Δd should be smaller than Δz. Even in a range where Δd is smaller than Δz, when Δd becomes too large, a probability that the focal point is almost positioned on the particle decreases, and as a result, the measurement accuracy of the particle diameter deteriorates. Δd should be determined in consideration of the particle diameter measurement accuracy.
In order to verify the validity of the condition of formula 12, the inventors further conduct a simulation study. In the simulation, it is assumed that there is one particle performing Brownian motion in a measurement region, a reflected signal waveform from a particle obtained when the measurement region is measured at the frame rate FR, z resolution Δz, and z pitch Δd is repeatedly calculated 10,000 times, and the obtained waveform is processed with a predetermined algorithm, thereby evaluating a particle detection success rate (number of times one particle is detected÷10,000). The simulation is performed under conditions that simulation is T=300 [K], μ=0.001 [Pa*s], and a size in the z direction of the measurement region is 1 mm.
As a method for avoiding counting the same particle a plurality of times, for example, proposed is an effective analysis method in which the particle is determined only when a reflected signal intensity equal to or greater than a threshold value is obtained at the z position that is continuous for a certain number of times (a first time) or more, and after that, when a reflected signal intensity equal to or lower than the threshold value is obtained at the z position that is continuous for a certain number of times (a second time) or more, it is determined that the measurement of the particle is completed. Here, the first time and the second time may be the same as each other or different from each other. For example, when the threshold value is defined as 0.1 and the first time and the second time are defined as 5, the number of detected particles is 1 at 0.1 fps, 0 at 1 fps, 1 at 10 fps, and 1 at 100 fps. From this example, it can be seen that when the frame rate is insufficient, a correct measurement result cannot be obtained.
A fitting parameter γ is about 4.46. The simulation result is accurately fitted by using formula 13, and it can be seen that the required frame rate is almost proportional to D (inversely proportional to a).
As described above, a large frame rate is desirable in order to cause the detection rate to be close to 1. On the other hand, when the frame rate becomes larger than necessary, noise increases such that the measurement accuracy of the particle diameter deteriorates. The reason is that a frequency band of a signal increases as the frame rate increases. In consideration of performance of a typical detector, it is desirable to set the frame rate to 100 fps or less. When the frame rate is converted to the value of γ, an upper limit of γ is about 10,000.
After a polarization state is adjusted by a λ/2 plate 105, the reference light is reflected by a reference light mirror 106 and incident on the polarization beam splitter 104 again. The measurement light is condensed by an objective lens 108 so that a focal point position of the measurement light matches a measurement position of a sample 110. An XY axis drive mechanism 107 scans the focal point position of the measurement light on an XY plane (a plane perpendicular to a depth direction of the sample 110). A Z axis drive mechanism 109 scans the focal point position of the measurement light along a Z axis direction (an optical axis direction of the measurement light). The reflected light reflected from the sample 110 is incident on the polarization beam splitter 104 again.
The reflected light and the reference light are multiplexed by the polarization beam splitter 104 to form synthesized light. The synthesized light is guided to an interference optical system 112 via a pinhole 111. The synthesized light is split into transmitted light and reflected light by a polarization beam splitter 113.
The reflected light passes through a λ/4 plate 114 of which optical axis is set to about 45 degrees with respect to the horizontal direction, and then is condensed by a condensing lens 115 and bifurcated by a Wollaston prism 116, thereby generating first interference light and second interference light having phase relationships different from each other by 180 degrees. A current differential type photodetector 117 detects the first interference light and the second interference light, and outputs a signal 122 proportional to a difference in intensities therebetween.
The transmitted light is transmitted through a λ/2 plate 118 of which optical axis is set to about 22.5 degrees with respect to the horizontal direction, and then is condensed by a condensing lens 119 and bifurcated by a Wollaston prism 120, thereby generating third interference light and fourth interference light having phase relationships different from each other by 180 degrees. A current differential type photodetector 121 detects the third interference light and the fourth interference light, and outputs a signal 123 proportional to a difference in intensities therebetween.
A signal processing unit 124 calculates a size distribution of a particle contained in the sample 110 based upon the signals 122 and 123. A principle of calculating the size distribution thereof is as described above. Details of a calculation procedure will be described later. A display unit 125 displays a calculation result by the signal processing unit 124.
The signal processing unit 124 acquires a plane image (an observation image on the XY plane) of the sample 110 while scanning the focal point position of the measurement light along the Z axis direction (step 1 in
The signal processing unit 124 specifies individual particles contained in the XY images by comparing the XY images adjacent to each other (step 2 in
The number of the XY image corresponds to the focal point position in the z direction. As described in
The signal processing unit 124 specifies the particle size by referring to correspondence relationship data in which a correspondence relationship between the reflected light intensity and the particle size is defined in advance. Specifically, the size of each particle is acquired by referring to the correspondence relationship data using the maximum reflected light intensity acquired in step 2. The signal processing unit 124 calculates the particle size distribution, and the display unit 125 displays a calculation result thereof. For example, as illustrated in the lower part of
A user inputs a measurement condition to the size distribution measurement device 100. The signal processing unit 124 receives an input of the measurement condition. As the measurement condition, for example, a size range of a particle to be measured, a range of the number of particles contained in the sample 110, a maximum measurement time, or the like can be considered.
The signal processing unit 124 acquires an XY image at each focal point position while scanning the focal point position in the Z axis direction (S1002 to S1003). An initial value of the focal point position is defined as z=0. When tracking individual particles in the XY image, the signal processing unit 124 may appropriately perform image sharpening processing or the like such as gain correction, noise processing, or the like (S1004).
The signal processing unit 124 specifies each particle in the XY image. The signal processing unit 124 calculates a plot (corresponding to the plot in the middle of
When measurement end conditions are reached, the processing proceeds to step S1009, and when the measurement end conditions are not reached, the processing returns to step S1001 and the same processing is repeated (S1008). The measurement end conditions referred to herein are, for example, the range of the number of particles and the maximum measurement time inputted in step S1001. The signal processing unit 124 calculates the size distribution of the particle and displays the calculated size distribution thereof on a screen via the display unit 125 (S1009).
The size distribution measurement device 100 according to the first embodiment acquires the XY image of the sample 110 at each focal point position while scanning the focal point position of the measurement light along the optical axis direction. The size distribution measurement device 100 further specifies each particle on the XY image individually, and obtains the particle size by using the reflected light intensity obtained from the position having the highest reflected light intensity of the respective focal point positions in the Z axis direction. Accordingly, it is possible to accurately obtain the size distribution of the particle that performs Brownian motion along the Z axis direction.
The size distribution measurement device 100 according to the first embodiment scans the measurement light faster than the Brownian motion speed of the particle in the Z axis direction. As a result, the number of particles can be accurately counted by following Brownian motion of the particle in the Z axis direction.
The size distribution measurement device 100 according to the third embodiment also changes the focal point position of the measurement light in each XYZ axis direction, such that the same operations as those of the first and second embodiments can be performed. A slit may be used instead of the pinholes 1702 and 1709. An LED (Light Emitting Diode) can be used instead of the laser light source 1701.
As a reference for switching the measurement method, for example, a ratio of the particle size to the spot size (particle size/spot diameter) can be used. As a result of measuring the particle size by using the reflected light intensity, when the particle size/the spot diameter is equal to or less than a threshold value (for example, 2 or less), the measurement result is adopted. When the particle size/the spot diameter is equal to or greater than the threshold value, the particle size is measured again by using the observation image. In the same manner, as a result of measuring the particle size by using the observation image, when the particle size/the spot diameter is equal to or greater than the threshold value, the measurement result is adopted, and when the particle size/the spot diameter is equal to or less than the threshold value, the particle size is measured again by using the reflected light intensity. Any one of the measurement methods may be used at a boundary portion between the two measurement methods.
For example, the signal processing unit 124 can acquire the observation image of the sample 110 by the same principle as that of a scanning type optical microscope. The reason is that the size distribution measurement device 100 uses an optical system that scans the measurement light.
In
The above embodiments describe the method of calculating the particle size distribution by following Brownian motion of the particle. On the other hand, from a viewpoint of measurement accuracy, it is desirable to prevent Brownian motion of the particle as much as possible. Brownian motion of the particle can be prevented to some extent, depending on a structure of a sample container that stores the sample 110. From the viewpoint of measurement accuracy, it is desirable to prevent an air bubble from being mixed in the sample 110. The air bubble can also be prevented from being mixed therein to some extent, depending on the structure of the sample container. Therefore, in the fifth embodiment and sixth and seventh embodiments of the present invention, a structure example of the sample container suitable for use in the first to fourth embodiments will be described.
The storage hole 1904 is a hole for storing the sample 110, and is formed so as to penetrate the base material 1901. A depth of the storage hole 1904 is, for example, about 300 μm to 1.5 mm. The air bubble storage portion 1902 is a gap portion protruding from a side surface of the storage hole 1904, and communicates with the storage hole 1904. Even though an air bubble is generated when the sample 110 is introduced into the storage hole 1904, the air bubble can be released to the air bubble storage portion 1902 (a gas discharge port). Accordingly, when measurement of the sample 110 in the storage hole 1904 is performed, the air bubble can be prevented from interfering with the measurement thereof.
In order to allow the air bubble to move smoothly from the storage hole 1904 to the air bubble storage portion 1902, an inner wall of the storage hole 1904 is desirably rounded. In
When a hole diameter of the storage hole 1904 is 2.5 to 4 mm and a shape of the air bubble storage portion 1902 in the top view (
As illustrated in
The horizontal reference marker 1905 can be used as a reference when the sample container 1900 is installed horizontally. For example, water containing the air bubble may be contained in the horizontal reference marker 1905, and the sample container 1900 may be installed so that the center of the air bubble is positioned at the center of the horizontal reference marker 1905.
The reference particle storage hole 1906 stores a reference particle at an almost central portion in the depth direction. A position of the reference particle can be used as a reference for the focal point position of the measurement light in the Z axis direction. A procedure of sealing the reference particle will be described later.
An upper surface of the base material 1901 is sealed by a sealing substrate 1908, and a lower surface thereof is sealed by a transmissive substrate 1909. That is, an upper surface and a bottom surface of the storage hole 1904 are also sealed by these substrates. The measurement light is emitted to the sample 110 in the storage hole 1904 via the transmissive substrate 1909. A portion 1910 of a space in the storage hole 1904 that is close to the surface of each substrate is excluded from a measurement target by an influence of interface reflection from the substrate. Therefore, actually measured portion is a region 1911.
The sample container 1900 according to the fifth embodiment includes the air bubble storage portion 1902 connected to the storage hole 1904. As a result, even though the air bubble is generated when the sample 110 is introduced into the storage hole 1904, the air bubble can be released to the air bubble storage portion 1902. As a result, when the measurement of the storage hole 1904 is performed, the air bubble can be prevented from interfering with the measurement thereof, such that the measurement accuracy is improved.
An outflow port 2103 is formed near the well 2102. A bottom portion of the well 2102 and a bottom portion of the outflow port 2103 are connected by a flow path 2104. An upper surface of the well 2102 and an upper surface of the outflow port 2103 are formed on the same plane as an upper surface of the plate 2101. After the sample 110 is introduced into the well 2102, an air bubble can be released through the flow path 2104 and the outflow port 2103.
A bottom surface of the plate 2101 is sealed by a transmissive substrate 2106, which also seals the bottom portion of the well 2102. The measurement light is emitted to the sample 110 in the well 2102 through the transmissive substrate 2106. After the well 2102 stores the sample 110, the upper surface of the well 2102 is covered with a seal 2105 (for example, a sealing tape). Accordingly, the sample 110 is sealed in the well 2102.
In the sample container 2100 according to the sixth embodiment, since an introduction port of the well 2102 is formed in the tapered shape, the sample 110 can be smoothly introduced thereinto. As a result, it is possible to prevent the air bubble from being generated when the sample 110 is introduced thereinto. Even though the air bubble is generated in the well 2102, the air bubble can be released through the flow path 2104 and the outflow port 2103. By the above-described structure, the particle size distribution can be calculated accurately even when the well 2102 of which pitch is specified is used.
The frame member 2301 is a member to which the sample container 2200 is attached. After introducing the sample 110 into the sample container 2200, the sample container 2200 is attached to the frame member 2301. The frame member 2301 includes a wall portion 2302 that covers one end of the flow path 2204. The wall portion 2302 may cover either one end of the flow path 2204 or both ends thereof depending on how the sample container 2200 is handled.
The present invention is not limited to the above-described embodiments, and includes various modifications. For example, the embodiments are described in detail in order to describe the present invention in an easy-to-understand manner, and are not necessarily limited to the one including all the described configurations. It is possible to replace a part of the configuration of one embodiment with a configuration of another embodiment, and it is also possible to add the configuration of another embodiment to the configuration of one embodiment. It is possible to add, delete, and replace another configuration with respect to a part of the configuration of each embodiment.
As the XY axis drive mechanism 107 in the first embodiment, for example, an acousto-optic deflector (AOD), a polygon mirror, a co-vibration type galvanometer mirror (a MEMS mirror), or the like can be used. The same configuration can be used for the XY axis drive in the second and third embodiments. As the Z axis drive mechanism 109 in the first embodiment, for example, a mechanism for driving a stage on which the sample 110 is placed can be used.
As the variable focus lens 126 in the second embodiment, for example, an ultrasonic varifocal lens, a liquid crystal variable focus mirror, a deformable mirror, or the like can be used. The same configuration can be used even when the XYZ axis drive mechanism 127 in the third embodiment scans the focal point in the Z axis direction.
The signal processing unit 124 can be configured by using hardware such as a circuit device or the like on which the function is loaded, and can also be configured by allowing a calculation device (Central Processing Unit: CPU or the like) to execute software on which the function is loaded.
100: size distribution measurement device
110: sample
112: interference optical system
124: signal processing unit
1900: sample container
2100: sample container
2200: sample container
2300: sample container
2301: frame member
2302: wall portion
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/JP2019/000281 | 1/9/2019 | WO | 00 |