CN110309524A - A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle - Google Patents
A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle Download PDFInfo
- Publication number
- CN110309524A CN110309524A CN201910214969.1A CN201910214969A CN110309524A CN 110309524 A CN110309524 A CN 110309524A CN 201910214969 A CN201910214969 A CN 201910214969A CN 110309524 A CN110309524 A CN 110309524A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- particle
- data
- oommf
- field
- Prior art date
- Legal status (The legal status 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 status listed.)
- Granted
Links
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of methods of magnetostatic coupling field intensity between calculating magnetic-particle to carry out modeling and relevant calculation to magnetic-particle using the micromagnetics simulation software OOMMF of open source as simulation calculation software.It is modeled respectively first with array structure of the OOMMF to isolated single magnetic-particle structure and more magnetic-particles and obtains non-vector OOMMF tables of data, then the dynamic magnetic spectrum of more magnetic-particles and isolated individual particle is calculated separately according to the OOMMF tables of data result of two kinds of structure simulations, then the frequency values for obtaining corresponding formant in two kinds of structure magnetic spectrums finally by the resonance frequency difference for obtaining two kinds of structures and combine corresponding formula to calculate the interaction field intensity in magnetic moment direction.
Description
Technical field
The invention belongs to micromagnetics to emulate field, be related to a kind of magnetostatic mutually based on the calculating of open source micromagnetics simulation software
The method of effect, and in particular to the calculation method for calculating the magnetostatic coupling between magnetic array unit.
Background technique
In recent years, the periodic array structure of magnet unit due to its unique magnetic property and ultra high density magnetic recording,
The application of sensor etc. and by extensive research and concern.With being constantly progressive for photoetching technique, people already can be with
The patterned magnetic array structure of Nano grade is prepared using the technology, while can also be prepared greatly using templating self-assembly technology
The magnetic array structure material of scale.In the array structure of preparation, the magnetic property of magnetic-particle is not only by the unit grain
The influence of size and pattern, and influenced by intergranular magnetostatic coupling intensity.Crystal grain in magnetic material
It interacts closely related with the macroscopic magnetization of material.Grains heredity in research material is for improving material property and grinding
System exploitation new material is most important.The difference of the grain microstructure due to caused by different preparation methods and preparation condition, so that
It is extremely complex and difficult for carrying out stringent mathematical computations to these interactions.In addition, for special grain structure, at all
It can not go to accurately calculate by formula.
Therefore in view of the defects existing in the prior art, it is really necessary to propose a kind of technical solution to overcome the prior art to exist
Defect.
Summary of the invention
In order to overcome the drawbacks of the prior art, the object of the present invention is to provide a kind of based on the simulation of micromagnetics simulation software
The method of magnetostatic coupling between calculation particle, this method is by the dynamic analog to individual particle system and more particle systems to obtain
Magnetostatic coupling without carrying out cumbersome mathematical computations, while using this method not between by the shape of particle and particle
The limitation of spacing can be good at the magnetostatic coupling power for calculating special construction.
In order to solve technical problem of the existing technology, technical scheme is as follows:
A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle includes the following steps: (1) selected magnetic material
Parameter is expected, by simulation software to single magnetic-particle Modelon Modeling；
(2) simulation calculation is carried out under no externally-applied magnetic field environment, and obtains the vector magnetic moment distribution file of stable state；
(3) the vector magnetic moment distribution file obtained in load step (2) in simulation model is adding weak pulse to rush field H (t) outside
=100exp (- 109T) it is emulated under the magnetic field environment of mT, the time interval that emulation data save is small as far as possible, to improve
Computational accuracy obtains non-vector tables of data after emulation；
(4) the non-vector tables of data for obtaining step (3) is taken out, and includes at least: external magnetic field H (t) data and a certain and H
(t) average magnetization M (t) data in perpendicular direction；(5) fast fourier transform method is utilized, by H (t), M (t) in frequency
H (ω) and M (ω) are transformed in domain respectively；(6) pass through complex permeability formula: χ (ω)=M (ω)/H (ω)=χ ' (ω)-i
χ " (ω) obtains magnetic spectrum χ " (ω)；
(7) the resonant frequency f reflected in magnetic spectrum χ " (ω) is extractedsingle(frequency corresponding to peak value in spectral line)；
(8) more magnetic-particle cell arrays are modeled again by simulation software, and uses two-dimentional boundary condition
(2DPBC) is extended；
(9) magnetic spectrum in array structure is obtained using the data processing method in (2-5) step, and extracts resonant frequency
farray；
(10) then pass through formulaIt is calculated along magnetic moment direction
On interaction field intensity Hu (as unit of mT).
Scheme as a further improvement in step (1), uses OOMMF simulation software modeling and simulating to calculate phase between particle
Interaction intensity.
Scheme as a further improvement, in step (3), using rapid decay weak pulse field H (t)=100exp (-
109T) mT is as excitation external magnetic field.
Scheme as a further improvement, in step (10), the variation of resonant frequency by comparing two kinds of structures goes to calculate
Particles' interaction intensity.
Scheme as a further improvement, using magnetostatic phase between open source micromagnetics simulation software OOMMF simulation calculating particle
When interaction, include the following steps:
(1) magnetic material parameter is selected, by writing building to single magnetic-particle unit for MIF script file in OOMMF
Mould.
(2) OOMMF software loads MIF file, is carried out in a manner of minimizing energy evolution under no externally-applied magnetic field environment
Simulation calculation, and obtain vector magnetic moment distribution file (* .omf file).
(3) in simulation model initialization module load previous step in * .omf file, outside plus weak pulse rush field H (t)=
100exp(-109T) mode that temporal evolution is carried out under the magnetic field environment of mT emulates, and the time interval of simulation stage is 1ps,
Non-vector OOMMF tables of data (* .odt file) is obtained after emulation.
(4) by external magnetic field H (t) data (Oxs_ScriptUZeeman::B* in file in the * .odt file of previous step
Column) and a certain mean magnetizing M (t) data (the Oxs_TimeDriver::m* column in file) with the perpendicular direction H (t) take
Out.
(5) fast fourier transform method is utilized, is transformed to H (ω) and M (ω) respectively in a frequency domain.
(6) pass through formula: χ (ω)=M (ω)/H (ω)=χ ' (ω)-i χ " (ω) obtains magnetic spectrum χ " (ω).
(7) the resonant frequency f reflected in magnetic spectrum χ " (ω) is extractedsingle(frequency corresponding to peak value in spectral line).
(8) more magnetic-particle cell arrays are modeled by OOMMF, and using two-dimentional boundary condition (2DPBC) into
Row extension.
(9) magnetic spectrum in array structure is obtained using the data processing method in (2-5) step, and extracts resonant frequency
farray。
(10) then pass through formulaIt calculates in magnetic moment direction
Interaction field intensity Hu (as unit of mT).
In above-mentioned method, the weak pulse field H (t) of the additional rapid decay in step (3) is applied to and material totality
The perpendicular direction of magnetic moment to make magnetic moment deviate settling position with low-angle under faint pulsed field, and then reduces magnetic moment
Dynamic and magnetostatic coupling error when static state.
In above-mentioned method, external magnetic field H (t) data are located at the Oxs_ScriptUZeeman: in file in step (4):
Bx, By or Bz column；It is determined by the application direction of pulsed field.Mean magnetizing M (t) data are located at the Oxs_ in file
TimeDriver::mx, my, mz column, by actual emulation calculate in select.
In above-mentioned method, formula is existing plural magnetic susceptibility calculation formula in step (6), which can be with table
Levy out the resonant frequency of magnetic moment inside magnet unit.
In above-mentioned method, it is extended in step (8) using two-dimentional boundary condition (2DPBC), to accurately calculate reality
Interaction strength in extensive cell array.
In above-mentioned method, the formula in step (10) is by magnetic moments precession intrinsic frequencySimplify and derives
?.The resonant frequency of magnetic moment has its effective field HeffDirectly determine, and the interacting field between magnetic-particle unit is to change to have
Imitate field HeffThe main reason for, so can reflect the size of interacting field by the change of magnetic moment resonant frequency.
Compared with the calculation method of existing magnetostatic coupling, calculation method of the invention, by strong using computer
Big computing capability, in addition the modeling and simulation computing capability of OOMMF software, can efficiently calculate intergranular interaction field strength
Degree.And it can solve the high complexity and limitation calculated using mathematical modeling mode.This method is not by magnetic material and particle
The freedom degree of the limitation of shape, simulation modeling is larger, and scalability is stronger.
Detailed description of the invention
Fig. 1 is the control interface of OOMMF software.
Fig. 2 is the simulation model for the taper nano dot array of particles established in example, and cone height is all provided with basal diameter
It is set to 130nm, anisotropy axis is along the z-axis direction.
Fig. 3 is to obtain individually isolated taper nano dot particle magnetic spectrum in example after data processing.
Fig. 4 is to show that grain spacing is the taper nano dot array of particles magnetic of 10-190nm in example after data processing
Compose group picture.
Fig. 5 is array junctions under the working frequency (dotted line is shown in figure) and 10 groups of spacing of single isolated unit in example
Resonant frequency (depicted with dots in figure) in structure.
Fig. 6 be in example particle interior zone corresponding to calculated main peak and secondary peak along the z-axis direction on phase interaction
With field intensity Hu.
Specific embodiment
Method as used in the following examples is conventional method unless otherwise specified.
Software as used in the following examples, calculation method etc. can obtain unless otherwise specified from business and internet channel
It arrives.This method is unrelated with magnet unit shape to embody, and selects more special conical nano-magnetic unit as specific real
Apply example.
Currently, establishing simulation model using the theory of micromagnetics is an effective means for studying magnetic systems.Its phase
Research cycle is short, scalability is strong is for the advantages of experiment research, in addition, since the method for experiment has many limitations
Property, so may be implemented more accurate to magnetic material system using the method that micromagnetics is simulated and efficiently study.People are wide
It is general that static and dynamic simulation study is carried out to material system using micromagnetics simulation software OOMMF.However, current OOMMF is soft
Without reference to the calculating for magnetostatic coupling intensity in the Simulation Application of part, thus can not be direct using existing simulation software
Calculate magnetostatic coupling.It is efficient in view of magnetostatic coupling great influence in systems and simulation calculation
Property.The calculation method of magnetostatic coupling intensity of the exploitation based on OOMMF has important value.
Based on current micromagnetics simulation software, the side of magnetostatic coupling field intensity between a kind of calculating magnetic-particle of the present invention
Method includes the following steps:
(1) magnetic material parameter is selected, by simulation software to single magnetic-particle Modelon Modeling；
(2) simulation calculation is carried out under no externally-applied magnetic field environment, and obtains the vector magnetic moment distribution file of stable state；
(3) the vector magnetic moment distribution file obtained in load step (2) in simulation model is adding weak pulse to rush field H (t) outside
=100exp (- 109T) it is emulated under the magnetic field environment of mT, the time interval that emulation data save is small as far as possible, to improve
Computational accuracy obtains non-vector tables of data after emulation；
(4) the non-vector tables of data for obtaining step (3) is taken out, and includes at least: external magnetic field H (t) data and a certain and H
(t) average magnetization M (t) data in perpendicular direction；(5) fast fourier transform method is utilized, by H (t), M (t) in frequency
H (ω) and M (ω) are transformed in domain respectively；(6) pass through complex permeability formula: χ (ω)=M (ω)/H (ω)=χ ' (ω)-i
χ " (ω) obtains magnetic spectrum χ " (ω)；
(7) the resonant frequency f reflected in magnetic spectrum χ " (ω) is extractedsingle(frequency corresponding to peak value in spectral line)；
(8) more magnetic-particle cell arrays are modeled by suture needle software again, and uses two-dimentional boundary condition
(2DPBC) is extended；
(9) magnetic spectrum in array structure is obtained using the data processing method in (2-5) step, and extracts resonant frequency
farray；
(10) then pass through formulaIt is calculated along magnetic moment direction
On interaction field intensity Hu (as unit of mT).
Below by taking OOMMF software as an example, detailed description of the present invention realizes process:
Example:
(1) the hexagonal barium ferrite (BaFe with uniaxial anisotropy is selected12O19) it is artificial material, by writing
MIF file in OOMMF goes to establish individually isolated conical nano-dot cells, is named as Cone_one.mif in this example,
Cone height and basal diameter are disposed as 130nm, and along the z-axis direction, OOMMF software is shown in attached drawing 1, circular cone nanometer to anisotropy axis
Point is referring to a certain unit in attached drawing 2.
(2) it is set as no externally-applied magnetic field environment in Cone_one.mif respective modules and minimizes the evolution side of energy
Formula, carries out simulation calculation using oxsii module loading Cone_one.mif, and the stable state of conical nano-dot cells is obtained after calculating
Magnetic moment distribution file (filename can customize), is saved as Cone_one_130nm.omf file in this example.
(3) modify the Cone_one.mif in (2) step, x-axis direction apply weak pulsed field H (t)=100exp (-
109T) magnetic field of mT is arranged for temporal evolution mode, the Cone_one_130nm.omf in previous step is loaded in initialization
File.The time interval of simulation stage is 1ps, simulation stage 8000.Non-vector OOMMF tables of data is obtained after emulation,
Tables of data, which is customized, saves as ConeDate_one_130nm.odt file.
(4) by external magnetic field H (t) data (Oxs_ in file in ConeDate_one_130nm.odt file
ScriptUZeeman::Bx column) and mean magnetizing M (t) data (Oxs_ in file with the perpendicular direction H (t)
TimeDriver::mx column) it takes out,
(5) Fast Fourier Transform is carried out using matlab software, is transformed to H (ω) and M (ω) respectively in a frequency domain.
(6) pass through formula: χ (ω)=M (ω)/H (ω)=χ ' (ω)-i χ " (ω) obtains individually isolated cone and receives
Rice dot element magnetic spectrum χ " (ω).As shown in figure 3, there are two formants in isolated circular cone nano-dot cells magnetic spectrum, they are located at
Different zones in unit.
(7) resonant frequency f corresponding to each peak of spectral line in magnetic spectrum χ " (ω) is extractedsingle.Calculate each peak value in Fig. 3
Abscissa be required resonant frequency.
(8) conical magnetic Nano dot element array is modeled by OOMMF, model is shown in attached drawing 2, will be between unit
Distance is to be set to 10-190nm (totally 10 groups of spacing) and be extended using two-dimentional boundary condition (2DPBC) with 20nm.
(9) conical nano-dot matrix array structure under 10 groups of spacing is obtained using the same data processing method in (2-5) step
In magnetic spectrum, and extract resonant frequency farray.Magnetic spectrum under 10 groups of spacing in array structure is shown in that Fig. 4, each structures exhibit go out
Two formants, this is consistent with isolated unit magnetic spectrum.The resonant frequency f of extractionarraySee attached drawing 5.
(10) pass through the data combination formula in Fig. 5:It is counted
It calculates.It obtains 10 groups of interacting field Hu (mT is unit) intensity datas along the z-axis direction and sees attached drawing 6, it is seen that corresponding to two peaks
The interaction field intensity in region is different, and difference is larger.This shows that the method in the present invention can be with Thought computation
The interaction of different zones in grain unit.
The above description of the embodiment is only used to help understand the method for the present invention and its core ideas.It should be pointed out that pair
For those skilled in the art, without departing from the principle of the present invention, the present invention can also be carried out
Some improvements and modifications, these improvements and modifications also fall within the scope of protection of the claims of the present invention.
The foregoing description of the disclosed embodiments enables those skilled in the art to implement or use the present invention.
Various modifications to these embodiments will be readily apparent to those skilled in the art, as defined herein
General Principle can be realized in other embodiments without departing from the spirit or scope of the present invention.Therefore, of the invention
It is not intended to be limited to the embodiments shown herein, and is to fit to and the principles and novel features disclosed herein phase one
The widest scope of cause.
Claims (4)
1. a kind of method of magnetostatic coupling field intensity between calculating magnetic-particle, which comprises the steps of:
(1) magnetic material parameter is selected, by simulation software to single magnetic-particle Modelon Modeling；
(2) simulation calculation is carried out under no externally-applied magnetic field environment, and obtains the vector magnetic moment distribution file of stable state；
(3) the vector magnetic moment distribution file obtained in load step (2) in simulation model, outside plus weak pulse rush field H (t)=
100exp(-109T) it is emulated under the magnetic field environment of mT, the time interval that emulation data save is small as far as possible, in terms of improving
Precision is calculated, non-vector tables of data is obtained after emulation；
(4) the non-vector tables of data for obtaining step (3) is taken out, and includes at least: external magnetic field H (t) data and a certain and H (t)
Average magnetization M (t) data in perpendicular direction；
(5) fast fourier transform method is utilized, H (t), H (t) are transformed to H (ω) and M (ω) respectively in a frequency domain；
(6) pass through complex permeability formula: χ (ω)=M (ω)/H (ω)=χ ' (ω)-i χ " (ω) obtains magnetic spectrum χ " (ω)；
(7) the resonant frequency f reflected in magnetic spectrum χ " (ω) is extractedsingle(frequency corresponding to peak value in spectral line)；
(8) more magnetic-particle cell arrays are modeled again by simulation software, and uses two-dimentional boundary condition (2DPBC)
It is extended；
(9) magnetic spectrum in array structure is obtained using the data processing method in (2-5) step, and extracts resonant frequency farray；
(10) then pass through formulaIt is calculated in magnetic moment direction
It interacts field intensity Hu (as unit of mT).
2. method described in claim 1, it is characterised in that: in step (1), use OOMMF simulation software modeling and simulating meter
Calculate particles' interaction intensity.
3. method described in claim 1, it is characterised in that: in step (3), use the weak pulse field H (t) of rapid decay
=100exp (- 109T) mT is as excitation external magnetic field.
4. method described in claim 1, it is characterised in that: in step (10), by the resonant frequency for comparing two kinds of structures
Variation goes to calculate particles' interaction intensity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910214969.1A CN110309524B (en) | 2019-03-21 | 2019-03-21 | Method for calculating strength of magnetostatic interaction field between magnetic particles |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910214969.1A CN110309524B (en) | 2019-03-21 | 2019-03-21 | Method for calculating strength of magnetostatic interaction field between magnetic particles |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110309524A true CN110309524A (en) | 2019-10-08 |
| CN110309524B CN110309524B (en) | 2022-11-11 |
Family
ID=68075682
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910214969.1A Active CN110309524B (en) | 2019-03-21 | 2019-03-21 | Method for calculating strength of magnetostatic interaction field between magnetic particles |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110309524B (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117524386A (en) * | 2024-01-04 | 2024-02-06 | 之江实验室 | Method and device for calculating magnetic alloy permeability based on micromagnetism and machine learning |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120043964A1 (en) * | 2010-05-31 | 2012-02-23 | Oliver Schreck | Method and system for simulating magnetic resonance imaging signals |
| CN103870654A (en) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing |
| CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
-
2019
- 2019-03-21 CN CN201910214969.1A patent/CN110309524B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120043964A1 (en) * | 2010-05-31 | 2012-02-23 | Oliver Schreck | Method and system for simulating magnetic resonance imaging signals |
| CN103870654A (en) * | 2014-03-26 | 2014-06-18 | 西安电子科技大学 | Electromagnetic scattering simulation method based on parallel moment method and physical optics mixing |
| CN107368642A (en) * | 2017-07-13 | 2017-11-21 | 武汉大学 | The multiple dimensioned multiple physical field coupling simulation method of metal increasing material manufacturing |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117524386A (en) * | 2024-01-04 | 2024-02-06 | 之江实验室 | Method and device for calculating magnetic alloy permeability based on micromagnetism and machine learning |
| CN117524386B (en) * | 2024-01-04 | 2024-06-04 | 之江实验室 | Method and device for calculating magnetic alloy permeability based on micromagnetism and machine learning |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110309524B (en) | 2022-11-11 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Crooker et al. | MHD model merging with IMF By: Lobe cells, sunward polar cap convection, and overdraped lobes | |
| Rowland et al. | Implementing quantum logic gates with gradient ascent pulse engineering: principles and practicalities | |
| CN108027409A (en) | Time domain mri | |
| Spalart et al. | Direct simulation and RANS modelling of a vortex generator flow | |
| CN105975713B (en) | A method of electronic circuit electromagnetic radiation is predicted using source method for reconstructing | |
| CN102213658A (en) | Geomechanical magnetic field testing device and method | |
| Wang et al. | New methodologies for computer simulation of paramagnetic resonance spectra | |
| CN102283649A (en) | Selective excitation sequential design method for magnetic resonance imaging in three-dimensional space | |
| Zhou et al. | Differentially rotating strange star in general relativity | |
| CN109856575B (en) | A kind of method, apparatus, equipment and the storage medium of the passive shimming of determining magnetic resonance | |
| CN102288929B (en) | Excitation pulse sequence generator for magnetic resonance system | |
| CN110119586B (en) | Axial conductivity anisotropy transient electromagnetic three-component three-dimensional FDTD forward modeling method | |
| CN110309524A (en) | A kind of method of magnetostatic coupling field intensity between calculating magnetic-particle | |
| Manzin et al. | A micromagnetic solver for large-scale patterned media based on non-structured meshing | |
| CN105241911B (en) | The method and device that low-field nuclear magnetic resonance analyzes fluid is simulated based on LBM | |
| Yang et al. | Demagnetization parameters evaluation of magnetic shields based on anhysteretic magnetization curve | |
| Tejima et al. | Magnetization analysis of stepped giant magneto impedance sensor using assembled domain structure model | |
| CN107219483A (en) | A kind of radial direction kurtosis anisotropic quantitative approach based on diffusion kurtosis imaging | |
| Li et al. | A blended method incorporating a multi-model for pulse wind tunnel aerodynamic identification considering large-scale aircraft | |
| Feix et al. | Is gravitational lensing by intercluster filaments always negligible? | |
| CN1580802A (en) | Space three-shaft motion magnetic field simulating device | |
| Ge et al. | Intelligent suppression of non-maneuvering magnetic interference of aeromagnetic UAV | |
| Washimi et al. | MHD structure of the heliosphere and its response to the 11-year solar cycle variations | |
| CN112444887B (en) | Electromagnetic detection information inversion device and method | |
| Wiesinger et al. | Evaluation of parallel transmit RF-shimming performance for 3 tesla whole-body imaging |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |