Magnetic Stimulation Method with Controllable Induced Field Direction

Information

  • Patent Application
  • 20230128566
  • Publication Number
    20230128566
  • Date Filed
    March 15, 2022
    2 years ago
  • Date Published
    April 27, 2023
    a year ago
Abstract
The disclosure discloses a magnetic stimulation method with a controllable induced field direction, and belongs to the technical field of noninvasive neural regulation. The method includes the following steps: S100, calculating currents i1j, i2j, i3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point Pt, the same below; S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E1, E2, E3; S300, calculating currents that may generate the electric field components E1, E2, E3 at the target point, I1j=E1i1j, I2j=E2i2j, I3j=E3i3j; S400, superimposing the currents of the three energization modes to obtain a resultant current Ij=I1j+I2j+I3j=E1i1j+E2i2j+E3i3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.
Description
TECHNICAL FIELD

The disclosure relates to a magnetic stimulation method with a controllable induced field direction, and belongs to the technical field of noninvasive neural regulation.


BACKGROUND

Magnetic stimulation is a neuromodulation technology. Its basic principle is to make a pulse current into a coil to generate a magnetic field that can penetrate through the scalp, the skull and other obstacle tissues, induce an electric field in the brain, and then induce an induced current to change membrane potentials of neurons, activate the neurons, and produce a series of physiological effects. Due to the non-invasive and safe characteristics of this technique, it is widely used in routine electrophysiological examinations, treatment of neurological diseases and scientific brain researches.


In the coil topology design of the existing magnetic stimulation method, most of them pay attention to the spatial distribution of the amplitude of the induced electric field, and its focusing degree and penetration depth are studied. However, the effect of magnetic stimulation is not only related to the strength of an induced electric field, but also closely related to a direction in which the induced electric field acts on the neurons.


The mechanism of the magnetic stimulation technology is not clear, and effective stimulation in the deep brain is still difficult. The direction and types of neurons in the deep nuclei of the brain are very complex. Different types of neurons have different sensitivity and activation thresholds to the direction of electric fields. At the same time, controlling the intensity and direction of an induced electric field can improve the validness and selectivity of deep neuron stimulation.


For existing coils, such as 8-shaped coils with the consistent electric field directions at a current convergence part, the direction of the induced electric field is changed by means of in-situ rotation, but this method is not flexible for the control of the direction of the electric field, and the direction can only be changed on a plane. When the position of the 8-shaped coil is changed to achieve a change in a direction of an electric field in a three-dimensional space, the 8-shaped coil will inevitably be moved away from an original stimulation target point, and the intensity of the electric field will decrease rapidly, failing to achieve a stimulation effect of neuron activation. Therefore, new stimulation methods and coil structures are needed to achieve direction-controllable magnetic stimulation.


SUMMARY

The disclosure provides a magnetic stimulation method with a controllable induced field direction, so as to solve the problem that if the position of a coil is changed, the intensity of an electric field will decrease rapidly, failing to achieve a stimulation effect of neuron activation in the prior art.


A magnetic stimulation method with a controllable induced field direction includes the following steps:


S100, calculating currents i1j, i2j, i3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point Pt, the same below;


S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E1, E2, E3;


S300, calculating currents that may generate the electric field components E1, E2, E3 at the target point, I1j=E1i1j, I2j=E2i2j, I3j=E3i3j;


S400, superimposing the currents of the three energization modes to obtain a resultant current Ij=I1j+I2j+I3j=E1i1j+E2i2j+E3i3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.


Further, before S100, the method further includes:


S000, constructing a magnetic stimulation coil group.


Further, in S000 specifically includes the following steps:


S001, constructing a three-dimensional space;


S002, placing the target point Pt in the three-dimensional space; setting q lead wires, the end points of ends of which are intersected, above the target point Pt, where all the lead wires form an included angle with an xoy plane and form included angles θ1, θ2, . . . θq with the forward direction of the z axis;


S003, closing the other ends of the q lead wires with the intersected end points by line segments or polygonal connecting lead wires to form a magnetic stimulation coil group which has n coils.


Further, in S100, a current is made into the magnetic stimulation coil group by the following energization modes:


in a first energization mode, currents in all the coils are I11, I12, I13, . . . , I1j, . . . , I1n, thus generating an induced electric field E1;


in a second energization mode, currents in all the coils are I21, I22, I23, . . . , I2j, . . . , I2n, thus generating an induced electric field E2;


in a third energization mode, currents in all the coils are I31, I32, I33, . . . , I3j, . . . , I3n, thus generating an induced electric field E3.


Further, in S200, specifically:


the coil group in each energization mode respectively generates combined induced electric fields E1(a1,b1,c1), E2(a2,b2,c2), E3(a3,b3,c3) at the target point Pt below coil intersections; the three induced electric fields are not coplanar, that is, E1, E2, E3 are linearly independent, where









"\[LeftBracketingBar]"





a
1




b
1




c
1






a
2




b
2




c
2






a
3




b
3




c
3






"\[RightBracketingBar]"



0.




Further, in S300, specifically:


obtaining corresponding unit-direction vector electric fields, e1=E1/|E1|, e2=E2/|E2|, e3=E3/|E3|, and generating currents needing to be made into the unit-direction vector electric fields at the target point,






i
11
=I
11
/|E
1
|,i
12
=I
12
/|E
1
|,i
13
=I
13
/|E
1
|, . . . ,i
1j
=I
1j
/|E
1
|, . . . ,i
1n
=I
1n
/|E
1|;






i
21
=I
21
/|E
2
|,i
22
=I
22
/|E
2
|,i
23
=I
23
/|E
2
|, . . . ,i
2j
=I
2j
/|E
2
|, . . . ,i
2n
=I
2n
/|E
2|;






i
31
=I
31
/|E
3
|,i
32
=I
32
/|E
3
|,i
33
=I
33
/|E
3
|, . . . ,i
3j
=I
3j
/|E
3
|, . . . ,i
3n
=I
3n
/|E
3|;


Further, in S400, specifically: according to the superposition principle of fields, three directional fundamental vector electric fields at the target point are used to combine a unit electric field e in arbitrary direction in a space,






e=λ
1
e
12e23e3,





0≤λi≤1


correspondingly, the currents in all the coils are also superposed,








i
1

=



λ
1



i
11


+


λ
2



i
21


+


λ
3



i
31








i
2

=



λ
1



i
12


+


λ
2



i
22


+


λ
3



i
32








i
3

=



λ
1



i
13


+


λ
2



i
23


+


λ
3



i
33











i
j

=



λ
1



i

1

j



+


λ
2



i

2

j



+


λ
3



i

3

j












i
n

=



λ
1



i
n


+


λ
2



i

2

n



+


λ
3



i

3

n









a ratio i1:i2:i3: . . . :ij: . . . :in of the currents made into the coils in the coil group is adjusted, that is, the direction of the electric field e at the target point is controlled; the sizes ki1, ki2, ki3, . . . :kij, . . . :kin of single currents are overall adjusted, that is, the size ke of the electric field at the target point is adjusted.


The disclosure has the following beneficial effects: The magnetic stimulation method with a controllable induced field direction of the disclosure breaks through the design concept of the existing coil, and can control the direction of the induced electric field at the target point in the three-dimensional space by only changing the current without moving the coil, and then activate, according to an electromagnetic field, neuron mechanism to provide flexible and accurate stimulation with directional characteristics for neuron, so as to provide a means for non-invasive selective neuron stimulation. This helps to improve the accuracy and validness of magnetic stimulation.





BRIEF DESCRIPTION OF FIGURES


FIG. 1A-FIG. 1C are projections of a coil center structure on an xoy plane;



FIG. 2 is a diagram of a coil center structure, which is composed of q wires, q=1, 2, 3, . . . n;



FIG. 3 is an example diagram of making a current;



FIG. 4 illustrates magnetic stimulation coil groups in four rectangular forms;



FIG. 5A-FIG. 5D illustrate an implementation method of controlling a direction of an induced electric field, where FIG. 5A is an energization mode 1; FIG. 5B is an energization mode 2; FIG. 5C is an energization mode 3; FIG. 5D is an energization mode 4;



FIG. 6 is a diagram of electric field distribution on a section under various energization modes;



FIG. 7A-FIG. 7D illustrate electric field distributions on a section when a resultant current is made and a size of a field value of an electric field of a target point, where FIG. 7A is E; FIG. 7B is an x-directional component of E; FIG. 7C is a y-directional component of E; FIG. 7D is a z-directional component of E;



FIG. 8A-FIG. 8E are example diagrams of five magnetic stimulation coil groups, where FIG. 8A illustrates four quadrilateral forms; FIG. 8B illustrates four triangular forms; FIG. 8C illustrates three triangular forms; FIG. 8D illustrates five triangular forms; FIG. 8E illustrates three hexagonal forms.





DETAILED DESCRIPTION

The technical solutions in the examples of the present disclosure will be clearly and completely described below in conjunction with the accompanying drawings in the examples of the present disclosure. Apparently, the described examples are only a part of the examples of the present disclosure, rather than all the examples. Based on the examples in the present disclosure, all other examples obtained by those of ordinary skill in the art without creative work shall fall within the protection scope of the present disclosure.


Referring to FIG. 1A-FIG. 1C, the disclosure provides a magnetic stimulation method with a controllable induced field direction, the magnetic stimulation method with the controllable induced field direction including the following steps:


S100, calculating currents i1j, i2j, i3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point Pt, the same below;


S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E1, E2, E3;


S300, calculating currents that may generate the electric field components E1, E2, E3 at the target point, I1j=E1i1j, I2j=E2i2j, I3j=E3i3j;


S400, superimposing the currents of the three energization modes to obtain a resultant current Ij=I1j+I2j+I3j=E1i1j+E2i2j+E3i3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.


Further, before S100, the method further includes:


S000, constructing a magnetic stimulation coil group.


Further, in S000 specifically includes the following steps:


S001, constructing a three-dimensional space;


S002, placing the target point Pt in the three-dimensional space; setting q lead wires, the end points of ends of which are intersected, above the target point Pt, where all the lead wires form an included angle with an xoy plane and form included angles θ1, θ2, . . . θq with the forward direction of the z axis;


S003, closing the other ends of the q lead wires with the intersected end points by line segments or polygonal connecting lead wires to form a magnetic stimulation coil group which has n coils.


Specifically, the physical principle of a magnetic intervention technology is Maxwell's equations. A time-varying electromagnetic field generated by a coil pulse current can be regarded as a magnetic quasi-static field. An induced electric field in an air domain has the following relationship with the current:







E
=

-



A



t







A
=



μ
0


4

π






Idl
r







E
=


-


μ
0


4

π







I



t






dl
r








where A is a vector magnetic potential; E is the induced electric field; IdI is a current element; r is a distance from the current element to a point to be determined; and μ0 is a vacuum permeability.


It can be known from the above formula that the induced electric field is parallel to the current element IdI, and has an opposite direction to that of the current element IdI. Based on the relationship between the current and the direction of an electric field and the superposition principle of fields, it is easy to construct a coil group structure that is capable of controlling arbitrary direction of an electric field in a space above a target point Pt through intersecting currents. In order to obtain an electric field in arbitrary direction in a three-dimensional space, three non-coplanar fundamental vectors are required, which provide electric field components in the x, y, and z directions of the Cartesian coordinate system.


Further, in S100, a current is made into the magnetic stimulation coil group by the following energization modes:


in a first energization mode, currents in all the coils are I11, I12, I13, . . . , I11, . . . , I1n, thus generating an induced electric field E1;


in a second energization mode, currents in all the coils are I21, I22, I23, . . . , I2j, . . . , I2n, thus generating an induced electric field E2;


in a third energization mode, currents in all the coils are I31, I32, I33, . . . , I3j, . . . , I3n, thus generating an induced electric field E3.


Specifically, an electric field in arbitrary direction is constructed on a two-dimensional plane. FIG. 1 shows a projection of a center structure of a coil group on an xoy plane. An intersection of lead wires is located directly above the target point. The lead wire structures alone or that overlap into a straight line cannot bring electric field components in the x and y directions, and it is impossible to construct, on the plane, a coil structure that controls the direction of the electric field by a current. The lead wire structures that radiate 2, 3, . . . , q lead wires from the center of the target point can bring components in the x and y directions on the plane, thus combining an electric field in arbitrary direction in the xoy plane.


In order to enable the coil to be able to generate an electric field component in the z direction for the electric field in arbitrary direction in the three-dimensional space, the lead wires need to form a certain included angle with the xoy plane. Therefore, the coil structure is a three-dimensional structure rather than a planar structure. As shown in FIG. 2, included angles between the lead wires and the forward direction of the z axis are θ1, θ2, . . . θq, respectively.


The lead wires are closed with polygonal connecting line segments, so that a coil group structure with a controllable induced field direction is obtained. The number q of the intersecting lead wires above the target point corresponds to the number n of coil groups.


The implementation principle of the magnetic stimulation method with the controllable induced field direction is as follows: three energization modes are used to make a current into the coil groups. In a first energization mode, the currents in all coils are I11, I12, I13, . . . , I1j, . . . , I1n; in a second energization mode, the currents in all the coils are I21, I22, I23, . . . , I2j, . . . , I2n; in a third energization mode, the currents in all the coils is I31, I32, I33, . . . , I3j, . . . , I3n (FIG. 3 is taken as an example, m=1, 2, 3 represent the energization modes).


Further, in S200, specifically:


the coil group in each energization mode respectively generates combined induced electric fields E1(a1,b1,c1), E2(a2,b2,c2), E3(a3,b3,c3) at the target point Pt below coil intersections; the three induced electric fields are not coplanar, that is, E1, E2, E3 are linearly independent, where









"\[LeftBracketingBar]"





a
1




b
1




c
1






a
2




b
2




c
2






a
3




b
3




c
3






"\[RightBracketingBar]"



0.




Further, in S300, specifically:


obtaining corresponding unit-direction vector electric fields, e1=E1/|E1|, e2=E2/|E2|, e3=E3/|E3|, and generating currents needing to be made into the unit-direction vector electric fields at the target point,






i
11
=I
11
/|E
1
|,i
12
=I
12
/|E
1
|,i
13
=i
13
/|E
1
|, . . . ,i
1j
=I
1j
/|E
1
|, . . . ,i
1n
=I
1n
/|E
1|;






i
21
=I
21
/|E
2
|,i
22
=I
22
/|E
2
|,i
23
=i
23
/|E
2
|, . . . ,i
2j
=I
2j
/|E
2
|, . . . ,i
2n
=I
2n
/|E
2|;






i
31
=I
31
/|E
3
|,i
32
=I
32
/|E
3
|,i
33
=i
33
/|E
3
|, . . . ,i
3j
=I
3j
/|E
3
|, . . . ,i
3n
=I
3n
/|E
3|;


Further, in S400, specifically: according to the superposition principle of fields, three directional fundamental vector electric fields at the target point are used to combine a unit electric field e in arbitrary direction in a space,






e=λ
1
e
12e23e3,





0≤λi≤1


correspondingly, the currents in all the coils are also superposed,








i
1

=



λ
1



i
11


+


λ
2



i
21


+


λ
3



i
31








i
2

=



λ
1



i
12


+


λ
2



i
22


+


λ
3



i
32








i
3

=



λ
1



i
13


+


λ
2



i
23


+


λ
3



i
33











i
j

=



λ
1



i

1

j



+


λ
2



i

2

j



+


λ
3



i

3

j












i
n

=



λ
1



i
n


+


λ
2



i

2

n



+


λ
3



i

3

n









a ratio i1:i2:i3: . . . :ij: . . . :in of the currents made into the coils in the coil group is adjusted, that is, the direction of the electric field e at the target point is controlled; the sizes ki1, ki2, ki3, . . . :kij, . . . :kin of single currents are overall adjusted, that is, the size ke of the electric field at the target point is adjusted.


As shown in FIG. 8, the coil group can have various forms. A quadrilateral form including four coils in an air domain is taken as an example to describe the implementation of the magnetic stimulation method for controlling the direction of an induced electric field. The schematic diagram of a coil group structure is shown in FIG. 4. The coil group is composed of 4 identical rectangles of a “book-shaped” structure, and edge lengths of the rectangles are a and b. Coils 1 and 2 are coplanar; coils 3 and 4 are coplanar and two coplanar coilsform a wing. An included angle between the plane of each wing coil and a negative direction of the z axis is α. The coil group is of a symmetrical structure, that is, an included angle between two wings is 2α.


In energization mode 1, the currents made into the coil group are I11, I12, I13, I14, and the magnitudes of the currents are equal: I11=I12=I13=I14; the current directions are clockwise, clockwise, counterclockwise, and counterclockwise when seen from the forward direction of the z axis. In energization modes 2 and 3, the magnitudes of the currents made into the coil group are also the same: I21=I22=I23=I24, I31=I32=I33=I34, and the current directions are as shown in FIG. 5.


In the air domain, the direction of the induced electric field is opposite to that of the current. FIG. 6 shows an electric field distribution on a certain section below the coils under the three energization modes. The magnitude of the current is 1 A; the frequency is 2500 Hz; the side lengths a and b are both 70 mm; the included angle is 60°. Table 1 shows electric field values at the target point. The three energization modes only bring electric fields in the x, y, and z directions at the target point.
















Energization
Im1, Im2, Im3, Im4
Ex
Ey
Ez


mode
(A)
(V/m)
(V/m)
(V/m)







1
1, 1, −1, −1
−1.50 × 10−3
0
0


2
1, −1, 1, −1
0
2.44 × 10−4
0


3
1, −1, −1, 1
0
0
8.61 × 10−4











    • Table 1 The induced electric field values of the target point (0, 0, −0.08) under all the energization modes





Note: m=1, 2, 3 are the energization modes; the current is clockwise positive and counterclockwise negative when seen from the forward direction of the z axis, and the same below.


Therefore, the current required to generate a unit fundamental vector electric field at the target can be calculated, as shown in Table 2.
















Energization
ex
ey
ez
im1, im2, im3, im4


mode
(V/m)
(V/m)
(V/m)
(A)







1
1
0
0
−0.67 × 103, −0.67 × 103, 0.67 × 103, 0.67 × 103


2
0
1
0
4.10 × 103, −4.10 × 103, 4.10 × 103, −4.10 × 103


3
0
0
1
1.16 × 103, −1.16 × 103, −1.16 × 103, 1.16 × 103











    • Table 2 Current required to generate a unit fundamental vector electric field at the target point (0, 0, −0.08)





A vector electric field E(2, 1, −1) V/m generated at the target point is taken as an example. Its components in a fundamental vector direction are E1=2, E2=1, E3=−1, Ij=2i1j+i2j−i3j (j=1, 2, 3, 4); a resultant current needing to be made to the coil group is (I1, I2, I3, I4)=(1.61×103, −4.27×103, 6.59×103, −3.93×103) A.



FIG. 7A˜D shows an electric field distribution of a section including a target point when a resultant current is made. An electric field value generated by the target point (0, 0, −0.08) is consistent with an expected value.


In a computational domain with a complex medium and a complex boundary, a finite element calculation method can be used to calculate an induced electric field, or actual measurement can be carried out to calibrate a fundamental vector of an electric field and a current made into it.


The above implementation examples are only used to help understand the method of the disclosure and its core idea. For those skilled in the art, according to the idea of the disclosure, several improvements and modifications can be made in the specific implementation and application scope. These improvements and modifications shall also fall within the protection scope of the disclosure.

Claims
  • 1. A magnetic stimulation method with a controllable induced field direction, the magnetic stimulation method with the controllable induced field direction comprising the following steps: S100, calculating currents i1j, i2j, i3j, j=1, 2, . . . n required to be made to generate a unit-direction vector electric field at a target point Pt, the same below;S200, decomposing a vector electric field E required at the target point to three fundamental vector directions to obtain electric field components E1, E2, E3;S300, calculating currents that may generate the electric field components E1, E2, E3 at the target point, I1j—E1i1j, I2j=E2i2j, I3j=E3i3j;S400, superimposing the currents of three energization modes to obtain a resultant current Ij=I1j+I2j+I3j=E1i1j+E2i2j+E3i3j to be made of each coil of a coil group, that is, generating a required electric field E at the target point for specific directional simulation.
  • 2. The magnetic stimulation method with the controllable induced field direction according to claim 1, wherein before S100, the method further comprises: S000, constructing a magnetic stimulation coil group.
  • 3. The magnetic stimulation method with the controllable induced field direction according to claim 2, wherein S000 specifically comprises the following steps: S001, constructing a three-dimensional space;S002, placing the target point Pt in the three-dimensional space; setting q lead wires, end points of ends of which are intersected, above the target point Pt, wherein all the lead wires form an included angle with an xoy plane and form included angles θ1, θ2, . . . θq with the forward direction of the z axis;S003, closing the other ends of the q lead wires with the intersected end points by line segments or polygonal connecting lead wires to form a magnetic stimulation coil group which has n coils.
  • 4. The magnetic stimulation method with the controllable induced field direction according to claim 1, wherein in S100, a current is made into the magnetic stimulation coil group by the following energization modes: in a first energization mode, currents in all the coils are I11, I12, I13, . . . , I1j, . . . , I1n, thus generating an induced electric field E1;in a second energization mode, currents in all the coils are I21, I22, I23, . . . , I2j, . . . , I2n, thus generating an induced electric field E2;in a third energization mode, currents in all the coils are I31, I32, I33, . . . , I3j, . . . , I3n, thus generating an induced electric field E3.
  • 5. The magnetic stimulation method with the controllable induced field direction according to claim 4, wherein in S200, specifically: the coil group in each energization mode respectively generates combined induced electric fields E1(a1,b1,c1), E2(a2,b2,c2), E3(a3,b3,c3) at the target point Pt below coil intersections; three induced electric fields are not coplanar, that is, E1, E2, E3 are linearly independent, wherein
  • 6. The magnetic stimulation method with the controllable induced field direction according to claim 5, wherein in S300, specifically: obtaining corresponding unit-direction vector electric fields, e1=E1/|E1|, e2=E2/|E2|, e3=E3/|E3|, and generating currents needing to be made into the unit-direction vector electric fields at the target point, i11=I11/|E1|,i12=I12/|E1|,i13=I13/|E1|, . . . ,i1j=I1j/|E1|, . . . ,i1n=I1n/|E1|;i21=I21/|E2|,i22=I22/|E2|,i23=I23/|E2|, . . . ,i2j=I2j/|E2|, . . . ,i2n=I2n/|E2|;i31=I31/|E3|,i32=I32/|E3|,i33=I33/|E3|, . . . ,i3j=I3j/|E3|, . . . ,i3n=I3n/|E3|;
  • 7. The magnetic stimulation method with the controllable induced field direction according to claim 6, wherein in S400, specifically: according to the superposition principle of fields, three directional fundamental vector electric fields at the target point are used to combine a unit electric field e in arbitrary direction in a space, e=λ1e1+λ2e2+λ3e3,0≤λi≤1correspondingly, the currents in all the coils are also superposed,
Priority Claims (1)
Number Date Country Kind
2021112494483 Oct 2021 CN national