RU2517762C2 - Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance - Google Patents

Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance Download PDF

Info

Publication number
RU2517762C2
RU2517762C2 RU2012139049/28A RU2012139049A RU2517762C2 RU 2517762 C2 RU2517762 C2 RU 2517762C2 RU 2012139049/28 A RU2012139049/28 A RU 2012139049/28A RU 2012139049 A RU2012139049 A RU 2012139049A RU 2517762 C2 RU2517762 C2 RU 2517762C2
Authority
RU
Russia
Prior art keywords
gradient
diffusion
pulse
pulses
amplitude
Prior art date
Application number
RU2012139049/28A
Other languages
Russian (ru)
Other versions
RU2012139049A (en
Inventor
Артур Михайлович Асланян
Дмитрий Александрович Давыдов
Владимир Леонидович Одиванов
Original Assignee
Общество с ограниченной ответственностью "МИКС"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Общество с ограниченной ответственностью "МИКС" filed Critical Общество с ограниченной ответственностью "МИКС"
Priority to RU2012139049/28A priority Critical patent/RU2517762C2/en
Publication of RU2012139049A publication Critical patent/RU2012139049A/en
Application granted granted Critical
Publication of RU2517762C2 publication Critical patent/RU2517762C2/en

Links

Images

Abstract

FIELD: measurement equipment.
SUBSTANCE: usage: for measurement of substance characteristics by NMR method. Substance: consists in the fact that to determine parameters of self-diffusion of the investigated sample they use a cycle of pulse sequence, made of the specified quantity of gradient pulses, duration, shape, amplitude and intervals between which are permanent, and two radio frequency pulses - 90-degree and 180-degree with the interval t between them, sent in the intervals between the third from the end and the last but one gradient pulse and between the last but one and the last gradient pulse, accordingly. Amplitude of the echo signal is measured at the moment of its maximum - via time t after 180-degree pulse or is produced by averaging by interval of time around this moment. To produce a diffusion drop, the cycles of measurement are repeated to change one of cycle parameters - amplitude of the gradient, duration of gradient pulses or interval between gradient pulses. The repetition period is determined by time of sample relaxation. Positive effect is achieved by setting a quasi-stationary condition in the series of gradient pulses, as a result of which the last pair of pulses in the measurement cycle of the sequence becomes close to equivalence.
EFFECT: higher accuracy of diffusion drop production and determination of self-diffusion coefficient, expanded range of its measurement.
8 dwg

Description

Изобретение относится к области методик измерения характеристик вещества методом ядерного магнитного резонанса (ЯМР) и может быть использовано для повышения точности и расширения диапазона измерения диффузионного затухания, диффузионного спада и определения коэффициента самодиффузии (КСД) веществ.The invention relates to the field of methods for measuring the characteristics of a substance by the method of nuclear magnetic resonance (NMR) and can be used to increase accuracy and expand the measurement range of diffusion attenuation, diffusion decay and determination of the self-diffusion coefficient (CSD) of substances.

Современные методики измерения КСД с использованием ЯМР основаны на регистрации поведения спин-системы ядер атомов исследуемого вещества в неоднородном магнитном поле. Наилучшие характеристики имеют методики с использованием импульсного градиента магнитного поля (ИГМП) [1], имеющие преимущества в чувствительности и диапазоне измерений перед методиками с использованием постоянного градиента. Рассмотрим наиболее распространенную методику определения КСД с ИГМП, в которой измерение диффузионного затухания проводится с использованием импульсной последовательности [2], схема цикла которой показана на рис.1.Modern methods of measuring CSD using NMR are based on recording the behavior of the spin-system of atomic nuclei of the investigated substance in an inhomogeneous magnetic field. The best characteristics are those using a pulsed magnetic field gradient (IGMP) [1], which have advantages in sensitivity and measurement range over techniques using a constant gradient. Consider the most common method for determining CSD with IGMP, in which the diffusion attenuation is measured using a pulse sequence [2], the cycle diagram of which is shown in Fig. 1.

В данном методе используются два одинаковых по амплитуде и длительности градиентных импульса с инвертирующим радиочастотным импульсом между ними. Первый импульс расфазирует спиновый пакет, а второй возвращает его к исходной фазе. Перемещение спинов в пакете в промежутке между импульсами приводит к тому, что часть спинов не восстанавливает свою фазу, за счет чего амплитуда эха уменьшается. Расчет амплитуды эха производится по формуле:This method uses two gradient pulses of the same amplitude and duration with an inverting RF pulse between them. The first pulse dephases the spin packet, and the second returns it to the initial phase. The movement of spins in the packet in the interval between pulses leads to the fact that some of the spins do not restore their phase, due to which the amplitude of the echo decreases. The calculation of the echo amplitude is performed according to the formula

Figure 00000001
Figure 00000001

где А(g) - амплитуда сигнала эха при амплитуде градиента = g, А(0) - амплитуда эха при нулевой амплитуде градиента, γ - гиромагнитное отношение, tg - длительность градиентных импульсов, td - интервал времени диффузии, D - коэффициент самодиффузии. Относительное уменьшение амплитуды эха в присутствии ИГМП называется диффузионным затуханием (ДЗ). КСД измеряется в системе СИ в единицах м2/с, его значение для воды при комнатной температуре 2.5·10-9 м2/с (2500 мкм2/с). Для определения КСД снимается зависимость амплитуды эха от величины градиента, строится график логарифма амплитуды от квадрата градиента. В цикле последовательности также может изменяться интервал диффузии td и длительность градиентных импульсов fg. Зависимость логарифма амплитуды сигнала эха от эффективного параметра

Figure 00000002
называется диффузионным спадом (ДС). В случае свободной диффузии однокомпонентного вещества график должен быть линейным, а КСД определяется по наклону графика ДС. На рис.2 показан ДС для тестового раствора CuSO4 в воде. Диффузионный спад для образца CuSO4 в воде, параметры последовательности: tg= 0.5мс, td= 5мс. По горизонтали - параметр, пропорциональный квадрату g. Величина КСД, вычисленная по наклону ДС = 2700 мкм2/с. По вертикали - условная амплитуда эха.where A (g) is the amplitude of the echo signal at the gradient amplitude = g, A (0) is the amplitude of the echo at zero gradient amplitude, γ is the gyromagnetic ratio, t g is the duration of the gradient pulses, t d is the diffusion time interval, D is the self-diffusion coefficient . The relative decrease in the echo amplitude in the presence of IGMP is called diffusion attenuation (RS). KSD is measured in the SI system in units of m 2 / s; its value for water at room temperature is 2.5 · 10 -9 m 2 / s (2500 μm 2 / s). To determine the CSD, the dependence of the echo amplitude on the gradient value is removed, and a plot of the logarithm of the amplitude on the square of the gradient is constructed. In the sequence cycle, the diffusion interval t d and the gradient pulse duration f g can also vary. The dependence of the logarithm of the amplitude of the echo signal on the effective parameter
Figure 00000002
called diffusion recession (DS). In the case of free diffusion of a single-component substance, the graph should be linear, and the CSD is determined by the slope of the DS graph. Figure 2 shows the DS for a test solution of CuSO 4 in water. Diffusion drop for a CuSO 4 sample in water, sequence parameters: t g = 0.5ms, t d = 5ms. Horizontal - a parameter proportional to the square g. The KSD value calculated by the slope of the DS = 2700 μm 2 / s. Vertical - the conditional amplitude of the echo.

Недостатком метода является его высокая чувствительность к неидентичности градиентных импульсов, которая приводит к неполному восстановлению фаз групп спинов, находящихся в разных частях образца. Измерение КСД маловязких образцов (значение КСД>1000 мкм2/с) затруднений не вызывает, особенно при не очень коротких временах релаксации Т2 (более 10-20 мс). Графики диффузионного спада получаются линейными в широком диапазоне амплитуд, поэтому измеренные значения КСД достаточно точны. Малые значения КСД (менее нескольких сотен мкм2/с) измерить значительно труднее, особенно при наличии коротковременных компонент релаксационного спада. Это связано с тем, что для получения диффузионного спада в достаточном диапазоне амплитуд требуются большие амплитуды и длительности градиентных импульсов и больший интервал диффузии, выбор которого ограничен временем релаксации образца. При использовании таких параметров последовательности получается диффузионный спад с искажением (см. рис.3, ДС для образца вакуумного масла, параметры: tg=4мс, td=10мс, вычисленный КСД = 36.8 мкм2/с).The disadvantage of this method is its high sensitivity to the non-identity of gradient pulses, which leads to incomplete restoration of the phases of the spin groups located in different parts of the sample. The measurement of the CSD of low-viscosity samples (the value of the CSD> 1000 μm 2 / s) is not difficult, especially with not very short relaxation times T2 (more than 10-20 ms). The diffusion decay plots are linear in a wide range of amplitudes, therefore, the measured values of the CSD are quite accurate. Small KSD values (less than a few hundred μm 2 / s) are much more difficult to measure, especially in the presence of short-term components of the relaxation decline. This is due to the fact that obtaining a diffusion decay in a sufficient range of amplitudes requires large amplitudes and durations of gradient pulses and a longer diffusion interval, the choice of which is limited by the relaxation time of the sample. Using these sequence parameters, a diffusion decline with distortion is obtained (see Fig. 3, DS for a vacuum oil sample, parameters: t g = 4ms, t d = 10ms, calculated KSD = 36.8 μm 2 / s).

Такие искажения могут быть связаны с неидентичностью двух градиентных импульсов, в результате чего не происходит полного восстановления фазы. Причин неидентичности может быть несколько.Such distortions can be associated with the non-identity of two gradient pulses, as a result of which complete phase recovery does not occur. There can be several reasons for non-identity.

1. При увеличении произведения амплитуды градиента на длительность импульса увеличивается нагрев опорного резистора в схеме стабилизации тока градиента, что приводит к изменению его сопротивления, которое проявляется в изменении амплитуды второго импульса.1. With an increase in the product of the gradient amplitude by the pulse duration, the heating of the reference resistor in the gradient current stabilization circuit increases, which leads to a change in its resistance, which manifests itself in a change in the amplitude of the second pulse.

2. Поле градиентных катушек оказывает влияние на ферромагнитные полюсные наконечники магнитной системы установки, магнитные параметры которых имеют собственные релаксационные свойства. Это также может приводить к неидентичности импульсов.2. The field of gradient coils affects the ferromagnetic pole pieces of the installation’s magnetic system, the magnetic parameters of which have their own relaxation properties. It can also lead to non-identical pulses.

Попытки уменьшить влияние этих факторов уже предпринимались [3] и приводили к улучшению результатов измерений, но полной компенсации не давали.Attempts to reduce the influence of these factors have already been undertaken [3] and led to improved measurement results, but did not give full compensation.

Целью изобретения является повышение точности получения диффузионного спада и определения коэффициента самодиффузии, расширение диапазона его измерения.The aim of the invention is to increase the accuracy of obtaining diffusion decline and determine the coefficient of self-diffusion, expanding the range of its measurement.

Технический результат достигается тем, что для обеспечения идентичности эффективного действия градиентов предлагается вместо двух подавать серию градиентных импульсов с одинаковой амплитудой, длительностью и интервалами между ними, из которых два последних будут включены в цикл последовательности Хана. В этом случае факторы, влияющие на эффективный вес импульсов, будут приходить в равновесие, в результате чего пара рабочих импульсов будет близка к идентичности, тогда как предварительные импульсы не будут сказываться на результате, поскольку действуют на невозмущенную спиновую систему.The technical result is achieved by the fact that in order to ensure the identity of the effective action of the gradients, it is proposed instead of two to apply a series of gradient pulses with the same amplitude, duration and intervals between them, of which the last two will be included in the Han sequence cycle. In this case, the factors affecting the effective weight of the pulses will come into equilibrium, as a result of which the pair of working pulses will be close to identity, while the preliminary pulses will not affect the result, since they act on the unperturbed spin system.

Заявляемое техническое решение осуществляют следующим образом.The claimed technical solution is as follows.

Исследуемый образец помещается в постоянное магнитное поле внутрь катушки, обеспечивающей воздействие на него переменного магнитного поля в направлении, перпендикулярном вектору постоянного поля. Эта же катушка или другая катушка, также имеющая ось, перпендикулярную вектору постоянного поля, используется для получения радиочастотного сигнала ЯМР. На образец воздействуют радиочастотными импульсами магнитного поля на частоте, близкой к частоте ларморовской прецессии спинов исследуемых ядер:The test sample is placed in a constant magnetic field inside the coil, which provides exposure to an alternating magnetic field in a direction perpendicular to the constant field vector. The same coil or another coil, also having an axis perpendicular to the constant field vector, is used to obtain the radio frequency NMR signal. The sample is exposed to radio frequency magnetic field pulses at a frequency close to the frequency of the Larmor precession of the spins of the nuclei under study:

Figure 00000003
,
Figure 00000003
,

где В0 - индукция магнитного поля, γ - гиромагнитное отношение для данного ядра. Для ядер атомов водорода (протонов) γ=42.58 МГц/Т.where B 0 is the magnetic field induction, γ is the gyromagnetic ratio for a given nucleus. For the nuclei of hydrogen atoms (protons) γ = 42.58 MHz / T.

В импульсной последовательности используются амплитуды и длительности радиочастотных импульсов, поворачивающие вектор намагниченности спин-системы образца на 90 и 180 градусов от исходного состояния относительно вектора постоянного поля. Для получения сигнала спинового эха используется последовательность, состоящая из 90-градусного и 180-градусного импульса с интервалом времени т между ними, амплитуда сигнала эха регистрируется через время 2 т после 90-градусного импульса. Для измерения ДЗ используются одинаковые импульсы градиента магнитного поля амплитудой g и длительностью tg между радиочастотными импульсами и перед регистрацией сигнала эха с интервалом td между ними. Для обеспечения равенства действия градиентных импульсов перед подачей 90-градусного радиочастотного импульса последовательности применяют серию из заданного количества градиентных импульсов амплитудой g и длительностью tg с интервалами между ними, равными td. Последовательность с двумя предварительными импульсами показана на рис.4.In the pulse sequence, the amplitudes and durations of the radio frequency pulses are used, which rotate the magnetization vector of the spin system of the sample by 90 and 180 degrees from the initial state relative to the constant field vector. To obtain a spin echo signal, a sequence consisting of a 90-degree and 180-degree pulse with a time interval t between them is used, the amplitude of the echo signal is recorded 2 t after a 90-degree pulse. To measure the RS, identical pulses of the magnetic field gradient with amplitude g and duration t g between the radio frequency pulses and before recording the echo signal with an interval t d between them are used. To ensure the equal effect of gradient pulses before applying a 90-degree radio frequency pulse, a series of a given number of gradient pulses of amplitude g and duration t g with intervals between them equal to t d is used . The sequence with two preliminary pulses is shown in Fig. 4.

Для получения диффузионного спада циклы последовательности повторяются с изменением одного из параметров последовательности: g, tg, td. Интервалы между циклами выбираются с учетом времени релаксации образца. В качестве диффузионного спада строится график зависимости амплитуды эха от параметра

Figure 00000002
ось ординат градуируется в логарифмическом масштабе.To obtain a diffusion decline, the sequence cycles are repeated with a change in one of the sequence parameters: g, t g , t d . The intervals between cycles are selected taking into account the relaxation time of the sample. As a diffusion decline, a plot of the echo amplitude versus parameter is plotted
Figure 00000002
the ordinate axis is graduated on a logarithmic scale.

На рис.5-8 показаны диффузионные спады, полученные с использованием последовательностей с различным количеством предварительных импульсов. Образец и временные параметры последовательности аналогичны рис.3. Все измерения производились на ЯМР релаксометре с рабочей частотой 18.5 МГц, длительностью радиочастотного 90-градусного импульса 3.6 мкс и длительностью спада свободной индукции около 500 мкс.Figure 5-8 shows the diffusion decays obtained using sequences with different numbers of preliminary pulses. The sample and time parameters of the sequence are similar to Fig. 3. All measurements were performed on an NMR relaxometer with an operating frequency of 18.5 MHz, a duration of a 90-degree radio frequency pulse of 3.6 μs, and a free induction decay of about 500 μs.

Преимуществами заявляемого способа являются:The advantages of the proposed method are:

• Улучшение формы диффузионного спада при исследовании образцов с низкими значениями КСД• Improving the shape of diffusion decline in the study of samples with low KSD values

• Снижение требований к температурной стабильности градиентных импульсов и экранировке магнитной системы• Reducing the requirements for temperature stability of gradient pulses and shielding of the magnetic system

• Расширение диапазона измерения КСД с заданной точностью• Extension of the range of measurement of KSD with a given accuracy

• Возможность выбора количества предварительных импульсов в зависимости от величины измеряемого КСД• Ability to select the number of preliminary pulses depending on the value of the measured KSD

Использование заявляемого способа в методиках измерений методом ЯМР позволит повысить точность и повторяемость результатов измерения параметров самодиффузии при сохранении аппаратных характеристик установки.Using the proposed method in NMR measurement techniques will improve the accuracy and repeatability of the measurement results of self-diffusion parameters while maintaining the hardware characteristics of the installation.

Использованные источникиUsed sources

1. J. Karger, W. Heing, Z. Exp.Techn. Phys., 1971, 19, 453.1. J. Karger, W. Heing, Z. Exp. Technology. Phys., 1971, 19, 453.

2. Е.О. Stejskal, J.E. Tanner, J. Chem. Phys.1965, 42, 288.2. E.O. Stejskal, J.E. Tanner, J. Chem. Phys. 1965, 42, 288.

3. Способ измерения диффузии адсорбированных молекул жидкости. Авт. свид. СССР №649996, опубл. 28.02.1979, бюл. №8.3. A method for measuring the diffusion of adsorbed liquid molecules. Auth. testimonial. USSR No. 649996, publ. 02/28/1979, bull. No. 8.

Claims (1)

Способ определения диффузионного затухания, диффузионного спада и коэффициента самодиффузии, включающий измерение амплитуды сигнала спинового эха, регистрируемой в момент времени 2 т или получаемой усреднением по интервалу времени вокруг этого момента в цикле импульсной последовательности, состоящем из 90-градусного импульса, 180-градусного импульса, подаваемого через время т, включающем два одинаковых градиентных импульса, подаваемых в промежутке между радиочастотными импульсами и между вторым радиочастотным импульсом и регистрацией амплитуды эха, отличающийся тем, что до подачи 90-градусного импульса формируют заданное количество градиентных импульсов, амплитудой, длительностью и формой совпадающих с градиентными импульсами, входящими в цикл измерительной последовательности, интервалы между которыми и между последним из них и первым градиентным импульсом цикла совпадают с интервалом между градиентными импульсами цикла. A method for determining diffusion attenuation, diffusion decay, and self-diffusion coefficient, including measuring the amplitude of the spin echo signal recorded at a time of 2 t or obtained by averaging over the time interval around this moment in a pulse sequence cycle consisting of a 90-degree pulse, a 180-degree pulse, fed through time m, including two identical gradient pulses applied between the radio frequency pulses and between the second radio frequency pulse and the recording echo amplitudes, characterized in that prior to the supply of a 90-degree pulse, a predetermined number of gradient pulses is formed with the amplitude, duration and shape of the pulse coinciding with the gradient pulses included in the measuring sequence cycle, the intervals between which and between the last of them and the first gradient pulse of the cycle coincide with interval between gradient pulses of the cycle.
RU2012139049/28A 2012-09-12 2012-09-12 Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance RU2517762C2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2012139049/28A RU2517762C2 (en) 2012-09-12 2012-09-12 Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2012139049/28A RU2517762C2 (en) 2012-09-12 2012-09-12 Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance

Publications (2)

Publication Number Publication Date
RU2012139049A RU2012139049A (en) 2014-03-20
RU2517762C2 true RU2517762C2 (en) 2014-05-27

Family

ID=50279980

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2012139049/28A RU2517762C2 (en) 2012-09-12 2012-09-12 Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance

Country Status (1)

Country Link
RU (1) RU2517762C2 (en)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU544901A1 (en) * 1974-12-10 1977-01-30 Казанский Ордена Трудового Красного Знамени Государственный Университет Имени В.И.Ульянова-Ленина Method for determining magnitude of impulse gradient magnetic field
US6462542B1 (en) * 2000-07-21 2002-10-08 Schlumberger Technology Corporation Nuclear magnetic resonance measurements and methods of analyzing nuclear magnetic resonance data
US6570382B1 (en) * 1999-12-10 2003-05-27 Schlumberger Technology Corporation Nuclear magnetic resonance method and logging apparatus
RU2378668C2 (en) * 2003-02-27 2010-01-10 Шлюмбергер Текнолоджи Бв Methods of interpreting diffusion-t2 maps obtained using nmr data

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
SU544901A1 (en) * 1974-12-10 1977-01-30 Казанский Ордена Трудового Красного Знамени Государственный Университет Имени В.И.Ульянова-Ленина Method for determining magnitude of impulse gradient magnetic field
US6570382B1 (en) * 1999-12-10 2003-05-27 Schlumberger Technology Corporation Nuclear magnetic resonance method and logging apparatus
US6462542B1 (en) * 2000-07-21 2002-10-08 Schlumberger Technology Corporation Nuclear magnetic resonance measurements and methods of analyzing nuclear magnetic resonance data
RU2378668C2 (en) * 2003-02-27 2010-01-10 Шлюмбергер Текнолоджи Бв Methods of interpreting diffusion-t2 maps obtained using nmr data

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Е.О. Stejskal, J.E. Tanner, J. Chem. Phys., 1965, 42, 288. *

Also Published As

Publication number Publication date
RU2012139049A (en) 2014-03-20

Similar Documents

Publication Publication Date Title
Petrov et al. T2 distribution mapping profiles with phase-encode MRI
US7378845B2 (en) NMR methods for measuring fluid flow rates
Rahmer et al. Rapid acquisition of the 3D MRI gradient impulse response function using a simple phantom measurement
US10557908B2 (en) Magnetic field monitoring of spiral echo train imaging
CN109254255B (en) Nuclear magnetic resonance transverse relaxation time spectrum measuring method and device and readable storage medium
Bartusek et al. Determination of pre-emphasis constants for eddy current reduction
Bottomley Basics of NMR
TWI540330B (en) Method and apparatus for detecting dynamic magnetic field fluctuation
Vashaee et al. A comparison of magnetic resonance methods for spatially resolved T2 distribution measurements in porous media
JP5926260B2 (en) Method and apparatus for adaptively changing the pulse interval in NMR moisture content measurement
Reynaud et al. Fast low‐specific absorption rate B0‐mapping along projections at high field using two‐dimensional radiofrequency pulses
Jarenwattananon et al. Breakdown of Carr-Purcell Meiboom-Gill spin echoes in inhomogeneous fields
JP2014151207A5 (en)
Casanova et al. NMR in inhomogeneous fields
Kittler et al. Pulsed second order field NMR for real time PGSE and single-shot surface to volume ratio measurements
Gurney et al. A simple method for measuring B0 eddy currents
RU2517762C2 (en) Pulse sequence to measure parameters of self-diffusion by method of nuclear magnetic resonance
JP2009085736A (en) Nmr measuring method
Lee et al. Alternating steady state free precession for estimation of current‐induced magnetic flux density: A feasibility study
US9864031B2 (en) Measurement of NMR characteristics of an object containing fast transversal relaxation components
Tan et al. A frequency measurement method using rising-falling edge of square wave for increasing proton magnetometer precision
Renou et al. Radio-frequency pulse calibration using the MISSTEC sequence
RU2696370C1 (en) Method of measuring longitudinal relaxation time in a current medium
Zhu et al. High spatial resolution NMR imaging of polymer layers on metallic substrates
CN108426909B (en) Method and device for synchronously and rapidly measuring nuclear magnetic resonance relaxation time and flow velocity

Legal Events

Date Code Title Description
PD4A Correction of name of patent owner
PC41 Official registration of the transfer of exclusive right

Effective date: 20210517