The disclosure relates to a magnetic stimulation method with a controllable induced field direction, and belongs to the technical field of noninvasive neural regulation.
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.
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
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
1+λ2e2+λ3e3,
0≤λi≤1
correspondingly, the currents in all the coils are also superposed,
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.
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
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:
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.
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
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 (
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
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
1+λ2e2+λ3e3,
0≤λi≤1
correspondingly, the currents in all the coils are also superposed,
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
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
In the air domain, the direction of the induced electric field is opposite to that of the current.
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.
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.
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.
Number | Date | Country | Kind |
---|---|---|---|
2021112494483 | Oct 2021 | CN | national |