WO2021132009A1 - 測定装置および測定方法 - Google Patents

測定装置および測定方法 Download PDF

Info

Publication number
WO2021132009A1
WO2021132009A1 PCT/JP2020/047122 JP2020047122W WO2021132009A1 WO 2021132009 A1 WO2021132009 A1 WO 2021132009A1 JP 2020047122 W JP2020047122 W JP 2020047122W WO 2021132009 A1 WO2021132009 A1 WO 2021132009A1
Authority
WO
WIPO (PCT)
Prior art keywords
physical field
measured
physical
magnetic resonance
measurement
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.)
Ceased
Application number
PCT/JP2020/047122
Other languages
English (en)
French (fr)
Inventor
エンスト デイヴィッド ヘルブスレブ
憲和 水落
義治 芳井
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Sumida Corp
Kyoto University NUC
Original Assignee
Sumida Corp
Kyoto University NUC
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 Sumida Corp, Kyoto University NUC filed Critical Sumida Corp
Priority to CN202080083340.9A priority Critical patent/CN114746764B/zh
Priority to US17/773,530 priority patent/US11808831B2/en
Priority to EP20904834.7A priority patent/EP4083647B1/en
Publication of WO2021132009A1 publication Critical patent/WO2021132009A1/ja
Anticipated expiration legal-status Critical
Ceased legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/02—Measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/032—Measuring direction or magnitude of magnetic fields or magnetic flux using magneto-optic devices, e.g. Faraday or Cotton-Mouton effect
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/543—Control of the operation of the MR system, e.g. setting of acquisition parameters prior to or during MR data acquisition, dynamic shimming, use of one or more scout images for scan plane prescription
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/10—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects by using electron paramagnetic resonance
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/24—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux
    • G01R33/26—Arrangements or instruments for measuring magnetic variables involving magnetic resonance for measuring direction or magnitude of magnetic fields or magnetic flux using optical pumping
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/28—Details of apparatus provided for in groups G01R33/44 - G01R33/64
    • G01R33/32—Excitation or detection systems, e.g. using radio frequency signals
    • G01R33/323—Detection of MR without the use of RF or microwaves, e.g. force-detected MR, thermally detected MR, MR detection via electrical conductivity, optically detected MR
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/443—Assessment of an electric or a magnetic field, e.g. spatial mapping, determination of a B0 drift or dosimetry
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/44—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using nuclear magnetic resonance [NMR]
    • G01R33/48—NMR imaging systems
    • G01R33/54—Signal processing systems, e.g. using pulse sequences ; Generation or control of pulse sequences; Operator console
    • G01R33/56—Image enhancement or correction, e.g. subtraction or averaging techniques, e.g. improvement of signal-to-noise ratio and resolution
    • G01R33/5608—Data processing and visualization specially adapted for MR, e.g. for feature analysis and pattern recognition on the basis of measured MR data, segmentation of measured MR data, edge contour detection on the basis of measured MR data, for enhancing measured MR data in terms of signal-to-noise ratio by means of noise filtering or apodization, for enhancing measured MR data in terms of resolution by means for deblurring, windowing, zero filling, or generation of gray-scaled images, colour-coded images or images displaying vectors instead of pixels
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N24/00—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects
    • G01N24/006—Investigating or analyzing materials by the use of nuclear magnetic resonance, electron paramagnetic resonance or other spin effects using optical pumping
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R33/00—Arrangements or instruments for measuring magnetic variables
    • G01R33/20—Arrangements or instruments for measuring magnetic variables involving magnetic resonance
    • G01R33/60—Arrangements or instruments for measuring magnetic variables involving magnetic resonance using electron paramagnetic resonance

Definitions

  • the present invention relates to a measuring device and a measuring method.
  • a magnetic measuring device performs magnetic measurement by optical detected magnetic resonance (ODMR) using electron spin resonance (see, for example, Patent Document 1).
  • a medium having a spin sublevel level and an optical transition level is irradiated with a high-frequency magnetic field (microwave) and light for excitation between spin sublevels and excitation between optical transitions, respectively.
  • Changes in the number of occupied numbers due to magnetic resonance between sub-levels are detected with high sensitivity by optical signals.
  • nitrogen in a diamond structure and electrons in a lattice defect (NVC: Nitrogen Vacancy Center) are initialized by photoexcitation by irradiation with a high-frequency magnetic field of about 2.87 GHz, and then three spin subs in the ground state.
  • a method for measuring a DC magnetic field using NVC there is a measuring method using a Ramsey Pulse Sequence.
  • the Ramsey pulse sequence (a) excitation light is applied to the NVC, (b) the first ⁇ / 2 pulse of microwaves is applied to the NVC, and (c) a predetermined time interval from the first ⁇ / 2 pulse.
  • a second ⁇ / 2 pulse of microwave is applied to the NVC at tt, (d) the measurement light is applied to the NVC to measure the amount of light emitted from the NVC, and (e) the magnetic flux density is determined based on the measured amount of light emitted. Can be estimated.
  • a measuring method of an alternating magnetic field using NVC there is a measuring method using a spin echo pulse sequence (Spin Echo Pulse Sequence).
  • spin echo pulse sequence (a) excitation light is applied to the NVC, (b) the first ⁇ / 2 pulse of the microwave is applied to the NVC at a phase of 0 degrees of the magnetic field to be measured, and (c) the microwave.
  • a ⁇ pulse is applied to the NVC at a phase of 180 degrees of the magnetic field to be measured, (d) a second ⁇ / 2 pulse of the microwave is applied to the NVC at a phase of 360 degrees of the magnetic field to be measured, and (e) the measurement light is applied to the NVC.
  • Is irradiated to measure the amount of light emitted from NVC and (f) the magnetic flux density can be estimated based on the measured amount of light emitted.
  • the magnetic field to be measured is estimated by using the electron spin quantum operation based on the Rabi vibration in the color center.
  • the time interval between the first ⁇ / 2 pulse and the ⁇ pulse and the time interval between the ⁇ pulse and the second ⁇ / 2 pulse can hold the information of the physical field to be measured. If the coherence time is longer than the coherence time, the information of the physical field to be measured disappears, so that it is difficult to accurately measure a long-period physical field such as an alternating magnetic field.
  • An object of the present invention is to obtain a measuring device and a measuring method for accurately measuring a long-period physical field by utilizing quantum operations in a predetermined quantum system.
  • the measuring device is a magnetic resonance member arranged in an AC physical field to be measured and capable of quantum manipulation in a predetermined quantum system, a coil for applying a microwave magnetic field to the magnetic resonance member, and the coil.
  • a high-frequency power supply that conducts a microwave current
  • an irradiation device that irradiates the magnetic resonance member with light in the DC physical field measurement sequence
  • an AC physical field to be measured from the magnetic resonance member in the DC physical field measurement sequence.
  • a detection device that detects the corresponding physical event and a DC physical field measurement sequence are executed a predetermined number of times, and in each of the DC physical field measurement sequences, the high-frequency power supply and the irradiation device are controlled, and the physical detected by the detection device. It includes a measurement control unit that specifies the detected value of an event, and a calculation unit that calculates the measurement results for a specific period of the AC physical field to be measured based on the detected values corresponding to a plurality of DC physical field measurement sequences.
  • a DC physical field measurement sequence is executed a predetermined number of times, and in each of the DC physical field measurement sequences, a physical event corresponding to the AC physical field to be measured is detected and the physical event thereof is detected.
  • the detected value of the event is specified, and
  • the measurement result for a specific period of the AC physical field to be measured is calculated based on the detected value corresponding to the plurality of DC physical field measurement sequences.
  • the DC physical field measurement sequence consists of a magnetic resonance member arranged in the AC physical field to be measured and capable of quantum manipulation in a predetermined quantum system, a coil that applies a microwave magnetic field to the magnetic resonance member, and the coil.
  • a high-frequency power supply that conducts a microwave current, an irradiation device that irradiates the magnetic resonance member with light, and a detection device that detects a physical event corresponding to the AC physical field to be measured from the magnetic resonance member. Is executed.
  • a measuring device and a measuring method for accurately measuring a long-period physical field can be obtained by utilizing quantum operations in a predetermined quantum system in a magnetic resonance member.
  • FIG. 1 is a diagram showing a configuration of a measuring device according to an embodiment of the present invention.
  • FIG. 2 is a timing chart for explaining the operation of the measuring device shown in FIG.
  • FIG. 3 is a timing chart for explaining the measurement sequence in FIG.
  • FIG. 1 is a diagram showing a configuration of a measuring device according to an embodiment of the present invention.
  • FIG. 2 is a timing chart for explaining the operation of the measuring device shown in FIG.
  • the measuring device shown in FIG. 1 includes a magnetic resonance member 1.
  • the magnetic resonance member 1 is a member arranged in the AC physical field to be measured and capable of quantum operation in a predetermined quantum system, and here includes a color center capable of electron spin quantum operation.
  • the number of color centers included in the magnetic resonance member 1 may be one or a plurality (that is, an ensemble).
  • the quantum system is electron spin and is capable of quantum manipulation based on Rabi oscillation.
  • a nuclear spin, a flux qubit, or the like may be used as this quantum system.
  • the magnetic resonance member 1 is made capable of forming a magnetic flux qubit or a member including a nuclear spin capable of quantum manipulation.
  • the AC physical field to be measured is an AC magnetic field having a specific cycle, but may be another physical field (electric field, temperature field, etc.) that alternates at a specific cycle.
  • the AC physical field to be measured may be an AC physical field having a single frequency or an AC physical field having a plurality of frequency components.
  • an optical detection magnetic resonance member is used as the magnetic resonance member 1, and in the magnetic resonance member 1, a physical event corresponding to the AC physical field to be measured (in the case of the optical detection magnetic resonance member, fluorescence). Light emission) is optically detected.
  • the physical event may be a change in electrical characteristics (such as a change in the resistance value of the magnetic resonance member 1) or may be electrically detected.
  • the optical detection magnetic resonance member used as the magnetic resonance member 1 is a plate material such as diamond having NVC and is fixed to the support plate 1a.
  • the measuring device shown in FIG. 1 includes a coil 2, a high frequency power supply 3, an irradiation device 4, and a detection device 5.
  • the coil 2 applies a microwave magnetic field to the magnetic resonance member 1.
  • the microwave frequency is set according to the type of the magnetic resonance member 1. For example, when the magnetic resonance member 1 is a diamond having NVC, the coil 2 applies a microwave magnetic field of about 2.87 GHz.
  • the high-frequency power source 3 conducts a microwave current (that is, a current for generating the above-mentioned microwave magnetic field) through the coil 2.
  • the irradiation device 4 irradiates the magnetic resonance member 1 with light (in this embodiment, excitation light having a predetermined wavelength and measurement light having a predetermined wavelength) in the DC physical field measurement sequence.
  • the detection device 5 detects a physical event corresponding to the AC physical field to be measured from the magnetic resonance member 1 in the DC physical field measurement sequence.
  • the detection device 5 is a light receiving device that detects fluorescence (the above-mentioned physical event) emitted from the magnetic resonance member 1 when the measurement light is irradiated.
  • the detection device 5 includes an electrode pair provided on the magnetic resonance member 1, measures the current value flowing from the electrode pair to the magnetic resonance member 1, and changes according to the AC physical field to be measured.
  • the resistance value of (the above-mentioned physical event) may be detected.
  • the measuring device shown in FIG. 1 includes an arithmetic processing unit 11.
  • the arithmetic processing unit 11 includes, for example, a computer, executes a program on the computer, and operates as various processing units.
  • the arithmetic processing unit 11 operates as the measurement control unit 21 and the arithmetic unit 22.
  • the measurement control unit 21 sets the DC physical field measurement sequence SQi (in a specific period within one cycle) for each cycle PEj of the AC physical field to be measured (here, the AC magnetic field to be measured). ) Execute a predetermined number of times, control the high frequency power supply 3 and the irradiation device 4 in each of the DC physical field measurement sequences SQ1, ..., SQn, and specify the detected value of the physical event detected by the detection device 5. .. In FIG. 2, a series of DC physical field measurement sequences SQ1, ..., SQn are started from the phase zero of the AC physical field to be measured, but a series of DC physical field measurement sequences SQ1, ..., SQn may be started from any unknown phase.
  • the time length of each DC physical field measurement sequence SQi and the time interval between the DC physical field measurement sequences are constant. Further, the time interval between the DC physical field measurement sequences may be zero or a predetermined time length other than zero.
  • the irradiation device 4 is provided with a laser diode or the like as a light source
  • the detection device 5 is provided with a photodiode or the like as a light receiving element
  • the measurement control unit 21 is obtained by amplifying the output signal of the light receiving element or the like. Based on the output signal of the detection device 5, the above-mentioned fluorescence light amount is specified as the above-mentioned detection value.
  • FIG. 3 is a timing chart for explaining the measurement sequence in FIG.
  • the above-mentioned DC physical field measurement sequence SQi is a Ramsey pulse sequence
  • the above-mentioned microwave includes two ⁇ / 2 pulses
  • the above-mentioned physical event is the two ⁇ /. It corresponds to the phase change of the electron spin of the color center (here, NVC) in the free aging motion at the time interval tt between two pulses.
  • the time interval tt between the two ⁇ / 2 pulses may be set according to (a) the effective lateral relaxation time (free induction decay time) T 2 * of the magnetic resonance member 1 or (b). )
  • the frequency of the AC physical field to be measured may be set to be included in the effective sensitivity frequency range (range of the upper limit frequency f or less corresponding to the time interval tt) obtained by the time interval tt.
  • the sensitivity of the magnetic field to be measured in the measurement sequence SQi becomes preferable.
  • the effective lateral relaxation time (free induction decay time) T 2 * has a value peculiar to the color center used for the measurement. Further, when the frequency of the AC physical field to be measured is lowered (that is, the time of one cycle is lengthened), the number of measurement sequences SQi per cycle is changed, or the interval between the measurement sequence SQi is changed. To.
  • the sensitivity H is expressed by the product of dBmin (the minimum value of the detectable magnetic field strength) and the square root of the measurement time Tm. Then, dBmin depends on T 2 * , the number of NV centers to be observed, and the like.
  • the DC physical field measurement sequence SQi described above is a Ramsey pulse sequence
  • the sensitivity H is substantially constant in the measurable frequency region regardless of the frequency of the alternating magnetic field to be measured. Specifically, the lower the frequency component in the magnetic field to be measured, the longer the measurement time per cycle, and the larger the number of measured values per specific period (here, one cycle), so that the sensitivity is substantial. It becomes constant. Further, as the number of measured values per specific period increases, the number of measured values used for curve fitting described later also increases, and a more accurate fitting curve can be obtained.
  • the calculation unit 22 calculates the measurement result of the AC physical field to be measured based on the detected values corresponding to the above-mentioned multiple DC physical field measurement sequences.
  • the calculation unit 22 has the above-mentioned multiple (n) times of the DC physical field measurement sequences SQ1 to SQn for each of the plurality of cycles of the AC physical field to be measured. Based on the detected value corresponding to, the measured value BMi of the measured AC physical field for one cycle is calculated, and (b) the average of the measured values of the measured AC physical field for one cycle, which is the same number as the multiple cycles, is calculated. The calculation is performed, and the noise in the measurement result for one cycle is attenuated.
  • the components other than the frequency (and the frequency that is an integral multiple of the frequency) that is the reciprocal of the one cycle are attenuated, so that the noise that is randomly generated at the high frequency is attenuated. Even if the frequency is an integral multiple of that, if it occurs in a random phase, that component will also be attenuated.
  • the frequency is an integral multiple of that, if it occurs in a random phase, that component will also be attenuated.
  • shotkey noise such temporal averaging is effective for noise attenuation, and if the magnetic resonance member 1 is an ensemble, a spatially averaged detection value can be obtained. Effective for noise attenuation.
  • the calculation unit 22 calculates the magnetic flux density from the above-mentioned fluorescence detection value according to the following equation.
  • Si is the detection value of fluorescence in the i-th measurement sequence SQi (that is, the corresponding phase i) in one cycle
  • BMi is the i-th measurement sequence SQi in the first cycle corresponding to the detected value Si.
  • a and b are constants
  • tt is the time interval between two ⁇ / 2 pulses (free aging motion time)
  • ⁇ is. , Magnetic rotation ratio (constant).
  • the measurement control unit 21 sets the DC physical field measurement sequence (here, Ramsey pulse sequence) SQ1, ..., SQn into one cycle of the alternating alternating physics field to be measured (here, alternating magnetic field).
  • a predetermined number of times n (n> 1) is executed per (in a specific period within one cycle (here, the entire area of one cycle)), and the calculation unit 22 performs each DC physics based on the detection value of the detection device 5.
  • the DC physical field measurement value BMi corresponding to the field measurement sequence SQi is derived.
  • the length of one cycle of the AC physical field to be measured here, the AC magnetic field
  • the measurement control unit 21 and the calculation unit 22 continuously and repeatedly execute this operation over a plurality of cycles PE1, ..., PEm (m> 1) of the AC physical field to be measured (here, an AC magnetic field), and perform the plurality of cycles.
  • the measured values for PE1, ..., PEm [BM1 (1), ..., BMn (1)], ..., [BM1 (m), ..., BMn (m)] are specified.
  • the calculation unit 22 calculates the average value BMAi of the measured values BMi (1), ..., BMi (m) for a plurality of cycles for each phase i in one cycle according to the following equation, and the calculation result [ BMA1, ..., BMAN] is used as the measurement result.
  • BMAi (BMi (1) + ... + BMi (m)) / m
  • the calculation unit 22 further performs curve fitting on the measurement results [BMA1, ..., BMAN], and derives a fitting curve (that is, a measurement waveform for one cycle).
  • the optimum value of the parameter in the functional form is derived by the least squares method based on the above measurement result.
  • the fitting curve is obtained by performing interpolation (interpolation or extrapolation) according to a predetermined method based on the above measurement result.
  • the calculation unit 22 derives the measurement result of the maximum amplitude of the AC physical field to be measured based on the fitting curve derived as described above. For example, as shown in FIG. 2, a measured value (estimated value) of the maximum amplitude Bac of the alternating magnetic field to be measured of a sine wave can be obtained.
  • the measurement control unit 21 (a) first irradiates the magnetic resonance member 1 with excitation light of a predetermined wavelength by the irradiation device 4, and the magnetic resonance member 1. After aligning the electron spin states of 1, (b) the coil 2 and the high-frequency power supply 3 apply a microwave magnetic field of the first ⁇ / 2 pulse to the magnetic resonance member 1, and (c) then for a predetermined time. After the interval tt elapses, a second ⁇ / 2 pulse microwave magnetic field is applied to the magnetic resonance member 1 by the coil 2 and the high frequency power supply 3, and (d) then, the irradiation device 4 emits the measurement light for projection measurement. At the same time as irradiating, the detection device 5 receives the fluorescence emitted by the magnetic resonance member 1 and detects the received light amount (detected light amount).
  • the direction of the electron spin changes in proportion to the time integration of the magnetic flux density of the external magnetic field (here, the alternating magnetic field to be measured). Therefore, from the detected light amount, the external magnetic field (here, the alternating magnetic field to be measured) is used. ) Magnetic flux density can be derived.
  • the time width tw of the ⁇ / 2 pulse is set to the time (about several tens of nanoseconds) for rotating the electron spin by ⁇ / 2, and is specified in advance from the period of the Rabi vibration of the magnetic resonance member 1.
  • the irradiation time of the excitation light and the irradiation time of the measurement light are about several microseconds to several tens of microseconds, respectively.
  • the above-mentioned time interval tt is set to several hundred microseconds or less. Further, the shorter the time interval between the excitation light irradiation and the first ⁇ / 2 pulse and the time interval between the second ⁇ / 2 pulse and the measurement light irradiation, the better.
  • the Ramsey pulse sequence (DC physical field measurement sequence SQi) is continuously executed in one cycle of the alternating magnetic field to be measured.
  • the irradiation of the measurement light in a certain measurement sequence SQi and the irradiation of the excitation light in the next measurement sequence SQ (i + 1) are performed separately, but both may be performed together.
  • the i-th measurement sequence SQi is executed in a certain phase Pi (may be an unknown phase) in a certain period PEj of the alternating magnetic field to be measured, and the i-th measurement sequence is also executed in the next period PE (j + 1). SQi is executed in the same phase Pi.
  • the measurement result obtained as described above may be stored as data in a storage device (not shown), transmitted to an external device, or displayed on a display device.
  • the length of the above-mentioned one measurement sequence SQi is set to half a cycle or less of the AC physical field to be measured, and for example, a low frequency AC physical field of about 1 kHz or less (particularly about 100 Hz or less) is a measurement target. Will be done.
  • the magnetic field strength is usually measured accurately.
  • the alternating magnetic field to be measured has a frequency component of about 1 kHz or less in such a case, the physical field can be measured satisfactorily.
  • the magnetic resonance member 1 is a member arranged in the AC physical field to be measured and capable of quantum operation in a predetermined quantum system.
  • the coil 2 applies a microwave magnetic field to the magnetic resonance member 1.
  • the high frequency power supply 3 conducts a microwave current through the coil 2.
  • the irradiation device 4 irradiates the magnetic resonance member 1 with light, and the detection device 5 detects a physical event corresponding to the AC physical field to be measured from the magnetic resonance member 1.
  • the measurement control unit 21 executes the DC physical field measurement sequence a predetermined number of times, controls the high frequency power supply 3 and the irradiation device 4 in each of the DC physical field measurement sequences, and the physical event detected by the detection device 5. Identify the detected value of.
  • the calculation unit 22 calculates the measurement results for a specific period (here, one cycle) of the AC physical field to be measured based on the detected values corresponding to the plurality of DC physical field measurement sequences.
  • a long-period physical field such as an alternating magnetic field can be accurately measured by utilizing the quantum operation in a predetermined quantum system in the magnetic resonance member 1.
  • the measured value of the physical field to be measured can be accurately obtained.
  • the above-mentioned multiple times of the DC physical field measurement sequence are uniformly executed (the length of the time interval tt in the DC physical field measurement sequence and the time interval of the DC physical field measurement sequence are measured AC. (When it is set to be constant regardless of the phase of the physical field), the above-mentioned multiple DC physical field measurement sequences are performed without synchronizing the multiple DC physical field measurement sequences with the AC physical field to be measured. May be good. In the case of the spin echo pulse sequence, synchronization with the alternating magnetic field to be measured is required.
  • the waveform of the alternating magnetic field to be measured is assumed to be a sine wave, but instead, the waveform of the alternating magnetic field to be measured is another waveform (triangle wave, saw wave, square wave, etc.). It may be a composite of a plurality of waveforms, etc.), as long as the period of the alternating magnetic field to be measured is known. Further, the alternating magnetic field to be measured may have a direct current component. Further, the same applies when a physical field other than the alternating magnetic field is the measurement target.
  • the time interval tt in the DC physical field measurement sequence may be changed according to the phase within one cycle of the AC physical field to be measured.
  • the phase of the AC physical field to be measured may be specified based on the phase of the external AC signal.
  • the measured AC physical field is generated by a physical phenomenon caused by an external AC signal
  • the measured AC physical field may be synchronized with the external AC signal.
  • the above-mentioned DC physical field measurement sequence is repeatedly executed with the time interval tt constant, and the waveform of the AC physical field to be measured is specified from the measurement results obtained by the multiple DC physical field measurement sequences.
  • the phase may be specified.
  • the measurement result for one cycle is obtained as an example of the specific period, but the measurement result for a cycle longer than one cycle (for example, for two cycles) may be obtained. Further, in the above embodiment, the measurement result for a period shorter than one cycle (that is, a specific phase range) may be obtained. In that case, for example, in each cycle of the above-mentioned plurality of cycles, only the detected value corresponding to the period may be acquired for each cycle and only the measured value corresponding to the period may be derived. However, it is preferable that the specific period is one cycle.
  • the calculation unit 22 estimates the actual rotation angle of the electron spin from the above-mentioned detected value, measured value, or measured result, and uses the detected value, measured value, or measured result as the electron spin. It may be converted to a value corresponding to the actual rotation angle of, and the above-mentioned curve fitting may be performed based on the converted value. At that time, the above-mentioned curve fitting may be performed directly by using a fitting function in consideration of the 2 ⁇ spin phase period without performing the conversion.
  • the average value BMi of the measured values BMi (1), ..., BMi (m) is regarded as the part of the phase i in the measurement result, but for a plurality of (m) cycles for a certain phase i. If there is no part of the measured values BMi (1), ..., BMi (m) and only the remaining part of BMi (1), ..., BMi (m) is derived, the remaining part The average value is calculated as the above-mentioned average value BMi.
  • BMAi (BMi (1) + ... + BMi (L-1) + BMi (L + 1) + ... + BMi (m)) / (m-1)
  • the measured values of a plurality of cycles in a specific period are averaged for each phase to obtain the measurement result, but when the time averaging is not required, the measured values for one cycle are specified.
  • the measured value of the period may be the measurement result (without averaging).
  • the present invention is applicable to, for example, a magnetic measuring device and a magnetic measuring method.

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Condensed Matter Physics & Semiconductors (AREA)
  • Engineering & Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Signal Processing (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Radiology & Medical Imaging (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Artificial Intelligence (AREA)
  • Power Engineering (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Measuring Magnetic Variables (AREA)

Abstract

磁気共鳴部材１は、被測定交流物理場内に配置され所定の量子系での量子操作の可能な部材である。コイル２および高周波電源３は、その磁気共鳴部材１にマイクロ波の磁場を印加する。照射装置４は、その磁気共鳴部材１に光を照射し、検出装置５は、その磁気共鳴部材１から、被測定交流物理場に対応する物理的事象を検出する。そして、測定制御部２１は、直流物理場測定シーケンスを所定複数回実行し、直流物理場測定シーケンスのそれぞれにおいて、検出装置５により検出された物理的事象の検出値を特定する。演算部２２は、その複数回の直流物理場測定シーケンスに対応する検出値に基づいて、被測定交流物理場の特定期間分の測定結果を演算する。

Description

測定装置および測定方法
　本発明は、測定装置および測定方法に関するものである。
　ある磁気計測装置は、電子スピン共鳴を利用した光検出磁気共鳴（ＯＤＭＲ：Optically Detected Magnetic Resonance）で磁気計測を行っている（例えば特許文献１参照）。
　ＯＤＭＲでは、スピンサブレベルの準位と光学遷移準位とをもつ媒質に高周波磁場（マイクロ波）と光をそれぞれスピンサブレベル間の励起と光遷移間の励起のために照射することで、スピンサブレベル間の磁気共鳴による占有数の変化などが光信号により高感度で検出される。通常、基底状態の電子が緑光で励起された後、基底状態に戻る際に赤光を発する。例えば、ダイヤモンド構造中の窒素と格子欠陥（ＮＶＣ：Nitrogen　Vacancy　Center）中の電子は、２．８７ＧＨｚ程度の高周波磁場の照射により、光励起により初期化された後では、基底状態中の３つのスピンサブレベルの中で一番低いレベル（ｍｓ＝０）から、基底状態中のそれより高いエネルギー軌道のレベル（ｍｓ＝±１）に遷移する。その状態の電子が緑光で励起されると、無輻射で基底状態中の３つのサブレベルの中で一番低いレベル（ｍｓ＝０）に戻るため発光量が減少し、この光検出より、高周波磁場により磁気共鳴が起こったかどうかを知ることができる。ＯＤＭＲでは、このような、ＮＶＣといった光検出磁気共鳴材料が使用される。
　そして、ＮＶＣを使用した直流磁場の測定方法として、ラムゼイパルスシーケンス（Ramsey Pulse Sequence）を使用した測定方法がある。ラムゼイパルスシーケンスでは、（ａ）励起光をＮＶＣに照射し、（ｂ）マイクロ波の第１のπ／２パルスをＮＶＣに印加し、（ｃ）第１のπ／２パルスから所定の時間間隔ｔｔでマイクロ波の第２のπ／２パルスをＮＶＣに印加し、（ｄ）測定光をＮＶＣに照射してＮＶＣの発光量を測定し、（ｅ）測定した発光量に基づいて磁束密度を見積もれる。
　また、ＮＶＣを使用した交流磁場の測定方法として、スピンエコーパルスシーケンス（Spin Echo Pulse Sequence）を使用した測定方法がある。スピンエコーパルスシーケンスでは、（ａ）励起光をＮＶＣに照射し、（ｂ）マイクロ波の第１のπ／２パルスを被測定磁場の位相０度でＮＶＣに印加し、（ｃ）マイクロ波のπパルスを被測定磁場の位相１８０度でＮＶＣに印加し、（ｄ）マイクロ波の第２のπ／２パルスを被測定磁場の位相３６０度でＮＶＣに印加し、（ｅ）測定光をＮＶＣに照射してＮＶＣの発光量を測定し、（ｆ）測定した発光量に基づいて磁束密度を見積もれる。
　このように、ラムゼイパルスシーケンスやスピンエコーパルスシーケンスでは、カラーセンタにおけるラビ振動に基づく電子スピン量子操作を利用して、被測定磁場が推定される。
特開２０１２－１１０４８９号公報
　しかしながら、スピンエコーパルスシーケンスでは、第１のπ／２パルスとπパルスとの時間間隔、およびπパルスと第２のπ／２パルスとの時間間隔が、被測定物理場の情報を保持できるスピンコヒーレンス時間よりも長くなると、被測定物理場の情報が消えてしまうため、長周期な、交流磁場などの物理場を正確に測定することは困難である。
　本発明は、所定の量子系での量子操作を利用して、長周期な物理場を正確に測定する測定装置および測定方法を得ることを目的とする。
　本発明に係る測定装置は、被測定交流物理場内に配置され所定の量子系での量子操作の可能な磁気共鳴部材と、その磁気共鳴部材にマイクロ波の磁場を印加するコイルと、そのコイルにマイクロ波の電流を導通させる高周波電源と、直流物理場測定シーケンスにおいて、その磁気共鳴部材に光を照射する照射装置と、直流物理場測定シーケンスにおいて、その磁気共鳴部材から、被測定交流物理場に対応する物理的事象を検出する検出装置と、直流物理場測定シーケンスを所定複数回実行し、直流物理場測定シーケンスのそれぞれにおいて、高周波電源および照射装置を制御し、検出装置により検出された物理的事象の検出値を特定する測定制御部と、複数回の直流物理場測定シーケンスに対応する検出値に基づいて、被測定交流物理場の特定期間分の測定結果を演算する演算部とを備える。
　本発明に係る測定方法は、（ａ）直流物理場測定シーケンスを所定複数回実行し、直流物理場測定シーケンスのそれぞれにおいて、被測定交流物理場に対応する物理的事象を検出してその物理的事象の検出値を特定し、（ｂ）その複数回の直流物理場測定シーケンスに対応するその検出値に基づいて、被測定交流物理場の特定期間分の測定結果を演算する。そして、直流物理場測定シーケンスは、被測定交流物理場内に配置され所定の量子系での量子操作の可能な磁気共鳴部材と、その磁気共鳴部材にマイクロ波の磁場を印加するコイルと、そのコイルにマイクロ波の電流を導通させる高周波電源と、その磁気共鳴部材に光を照射する照射装置と、その磁気共鳴部材から、被測定交流物理場に対応する物理的事象を検出する検出装置とを使用して実行される。
　本発明によれば、磁気共鳴部材における所定の量子系での量子操作を利用して、長周期な物理場を正確に測定する測定装置および測定方法が得られる。
図１は、本発明の実施の形態に係る測定装置の構成を示す図である。 図２は、図１に示す測定装置の動作について説明するタイミングチャートである。 図３は、図２における測定シーケンスについて説明するタイミングチャートである。
　以下、図に基づいて本発明の実施の形態を説明する。
　図１は、本発明の実施の形態に係る測定装置の構成を示す図である。図２は、図１に示す測定装置の動作について説明するタイミングチャートである。図１に示す測定装置は、磁気共鳴部材１を備える。磁気共鳴部材１は、被測定交流物理場内に配置され所定の量子系での量子操作の可能な部材であって、ここでは、電子スピン量子操作の可能なカラーセンタを含む。なお、磁気共鳴部材１に含まれるカラーセンタは１つでも複数（つまり、アンサンブル）でもよい。この実施の形態では、この量子系は、電子スピンであって、ラビ振動に基づく量子操作が可能である。なお、電子スピンの代わりに、この量子系として、原子核スピン、磁束量子ビットなどを使用してもよい。その場合には、磁気共鳴部材１は、量子操作が可能な原子核スピンを含んだものや、磁束量子ビットを形成可能なものとされる。また、この実施の形態では、被測定交流物理場は、特定周期の交流磁場であるが、特定周期で交番する他の物理場（電場、温度場など）であってもよい。なお、被測定交流物理場は、単一周波数の交流物理場でもよいし、複数の周波数成分を有する交流物理場でもよい。
　また、この実施の形態では、磁気共鳴部材１として、光検出磁気共鳴部材が使用され、磁気共鳴部材１において、被測定交流物理場に対応する物理的事象（光検出磁気共鳴部材の場合、蛍光発光）が光学的に検出される。なお、物理的事象は、電気特性の変化（磁気共鳴部材１の抵抗値の変化など）であってもよく、電気的に検出されてもよい。
　ここでは、磁気共鳴部材１として使用される光検出磁気共鳴部材は、ＮＶＣを有するダイヤモンドなどの板材であって、支持板１ａに固定されている。
　また、図１に示す測定装置は、コイル２、高周波電源３、照射装置４、および検出装置５を備える。
　コイル２は、磁気共鳴部材１にマイクロ波の磁場を印加する。マイクロ波の周波数は、磁気共鳴部材１の種類に応じて設定される。例えば、磁気共鳴部材１が、ＮＶＣを有するダイヤモンドである場合、コイル２は、２．８７ＧＨｚ程度のマイクロ波磁場を印加する。高周波電源３は、コイル２にマイクロ波の電流（つまり、上述のマイクロ波の磁場を生成するための電流）を導通させる。
　照射装置４は、直流物理場測定シーケンスにおいて、磁気共鳴部材１に光（この実施の形態では、所定波長の励起光と所定波長の測定光）を照射する。検出装置５は、直流物理場測定シーケンスにおいて、磁気共鳴部材１から、被測定交流物理場に対応する物理的事象を検出する。
　この実施の形態では、検出装置５は、測定光の照射時において磁気共鳴部材１から発せられる蛍光（上述の物理的事象）を検出する受光装置である。なお、検出装置５は、磁気共鳴部材１に設けられた電極対を備え、電極対から磁気共鳴部材１に流れる電流値を測定して、被測定交流物理場に応じて変化する磁気共鳴部材１の抵抗値（上述の物理的事象）を検出するようにしてもよい。
　さらに、図１に示す測定装置は、演算処理装置１１を備える。演算処理装置１１は、例えばコンピュータを備え、プログラムをコンピュータで実行して、各種処理部として動作する。この実施の形態では、演算処理装置１１は、測定制御部２１および演算部２２として動作する。
　測定制御部２１は、例えば図２に示すように、被測定交流物理場（ここでは、被測定交流磁場）の１周期ＰＥｊごとに、直流物理場測定シーケンスＳＱｉを（１周期以内の特定期間において）所定複数回実行し、直流物理場測定シーケンスＳＱ１，・・・，ＳＱｎのそれぞれにおいて、高周波電源３および照射装置４を制御し、検出装置５により検出された物理的事象の検出値を特定する。なお、図２では、被測定交流物理場の位相ゼロから、一連の直流物理場測定シーケンスＳＱ１，・・・，ＳＱｎが開始されているが、一連の直流物理場測定シーケンスＳＱ１，・・・，ＳＱｎは、任意の未知の位相から開始させてもよい。
　なお、この実施の形態では、各直流物理場測定シーケンスＳＱｉの時間長および直流物理場測定シーケンス間の時間間隔は一定である。また、直流物理場測定シーケンス間の時間間隔は、ゼロでも、ゼロ以外の所定時間長でもよい。
　例えば、照射装置４は、レーザーダイオードなどを光源として備え、検出装置５は、フォトダイオードなどを受光素子として備え、測定制御部２１は、受光素子の出力信号に対して増幅などを行って得られる検出装置５の出力信号に基づいて、上述の蛍光の光量を上述の検出値として特定する。
　図３は、図２における測定シーケンスについて説明するタイミングチャートである。
　この実施の形態では、上述の所定の直流物理場測定シーケンスＳＱｉとして、例えば図３に示すようなラムゼイパルスシーケンスが適用される。ただし、これに限定されるものではない。
　この実施の形態では、上述の直流物理場測定シーケンスＳＱｉは、ラムゼイパルスシーケンスとされるため、上述のマイクロ波として２つのπ／２パルスを含み、上述の物理的事象は、その２つのπ／２パルスの間の時間間隔ｔｔにおける自由歳差運動でのカラーセンタ（ここでは、ＮＶＣ）の電子スピンの位相変化に対応する。そして、その２つのπ／２パルスの間の時間間隔ｔｔは、（ａ）磁気共鳴部材１の実効横緩和時間（自由誘導減衰時間）Ｔ２ ＊に応じて設定されてもよいし、（ｂ）被測定交流物理場の周波数が当該時間間隔ｔｔによって得られる有効感度周波数範囲（時間間隔ｔｔに対応する上限周波数ｆ以下の範囲）に含まれるように設定されてもよい。
　上述の時間間隔ｔｔが実効横緩和時間（自由誘導減衰時間）Ｔ２ ＊の約半分に等しいとき、測定シーケンスＳＱｉにおける被測定磁場の感度が好ましくなる。なお、実効横緩和時間（自由誘導減衰時間）Ｔ２ ＊は、当該測定に使用されるカラーセンタに固有な値を有する。また、被測定交流物理場の周波数が低くされる（つまり、１周期の時間が長くされる）場合、１周期あたりの測定シーケンスＳＱｉの数が変更されるか、測定シーケンスＳＱｉ間の間隔が変更される。
　ここで、感度Ｈは、ｄＢｍｉｎ（検出可能な磁場強度の最小値）と測定時間Ｔｍの平方根との積で表現される。そして、ｄＢｍｉｎは、Ｔ２ ＊や観測するＮＶセンタの数などに依存する。この実施の形態では、上述の直流物理場測定シーケンスＳＱｉは、ラムゼイパルスシーケンスとされるため、測定可能な周波数領域において、感度Ｈは、被測定交流磁場の周波数に拘わらず略一定となる。具体的には、被測定磁場における周波数成分が低周波数になるほど、１周期あたりの測定時間が長くなり、特定期間（ここでは、１周期）あたりの測定値の数が多くなるので、感度が実質的に一定となる。また、特定期間あたりの測定値の数が増えると、後述のカーブフィッティングに使用される測定値も多くなり、より正確なフィッティング曲線が得られる。
　演算部２２は、上述の複数回の直流物理場測定シーケンスに対応する検出値に基づいて、被測定交流物理場の測定結果を演算する。
　この実施の形態では、演算部２２は、例えば図２に示すように、（ａ）被測定交流物理場の複数周期のそれぞれについて、上述の複数（ｎ）回の直流物理場測定シーケンスＳＱ１～ＳＱｎに対応する検出値に基づいて１周期分の被測定交流物理場の測定値ＢＭｉを演算し、（ｂ）その複数周期と同数の、１周期分の被測定交流物理場の測定値の平均を演算し、これにより、１周期分の測定結果におけるノイズを減衰させる。具体的には、このように平均化することで、その１周期の逆数である周波数（およびその整数倍の周波数）の成分以外は減衰するため、高周波でランダムに発生するノイズは減衰する。なお、その整数倍の周波数の成分であっても、ランダムな位相で発生する場合には、その成分も減衰する。例えば、ショットキーノイズの場合、このような時間的な平均化がノイズの減衰に有効であり、また、磁気共鳴部材１がアンサンブルであれば空間的な平均化がなされた検出値が得られ、ノイズの減衰に有効である。
　例えば、ｍ周期における第ｊ周期（ｊ＝１，・・・，ｍ）の測定値をＢＭ１（ｊ），・・・，ＢＭｎ（ｊ）とすると、［ＢＭ１（１）～ＢＭ１（ｍ）の平均値］，・・・，［ＢＭｎ（１）～ＢＭｎ（ｍ）の平均値］が平均化後の測定結果とされる。
　具体的には、この実施の形態では、演算部２２は、検出装置５で得られる上述の蛍光の検出値に基づいて、１周期分の被測定交流磁場の強度の測定値ＢＭｉ（ｉ＝１，・・・，ｎ）を演算する。
　例えば、演算部２２は、次式に従って、上述の蛍光の検出値から磁束密度を計算する。
　Ｓｉ＝［（ａ＋ｂ）／２］＋［（ａ－ｂ）／２］×ｃｏｓ（γ・ＢＭｉ・ｔｔ）
　ここで、Ｓｉは、１周期におけるｉ番目の測定シーケンスＳＱｉ（つまり、対応する位相ｉ）での蛍光の検出値であり、ＢＭｉは、検出値Ｓｉに対応する１周期におけるｉ番目の測定シーケンスＳＱｉ（つまり、対応する位相ｉ）での磁束密度であり、ａ，ｂは、定数であり、ｔｔは、２つのπ／２パルスの間の時間間隔（自由歳差運動時間）であり、γは、磁気回転比（定数）である。なお、ａおよびｂは、ＢＭｉまたはｔｔを変化させたときのＳｉの最大値および最小値であり、例えば実験で、既知で一定のＢＭｉに対して、ｔｔを変化させたときのＳｉを測定することで、特定可能である。なお、ａは、Ｂ＝０としたときの検出光量として求めてもよい。
　次に、当該実施の形態に係る磁場測定装置の動作について説明する。
　例えば図２に示すように、測定制御部２１は、直流物理場測定シーケンス（ここでは、ラムゼイパルスシーケンス）ＳＱ１，・・・，ＳＱｎを、被測定交流物理場（ここでは交流磁場）の１周期あたり（１周期内の特定期間（ここでは１周期の全域）に）所定の複数回数ｎ（ｎ＞１）を実行し、演算部２２は、検出装置５の検出値に基づいて、各直流物理場測定シーケンスＳＱｉに対応する直流物理場測定値ＢＭｉを導出する。なお、ここで、被測定交流物理場（ここでは交流磁場）の１周期の長さは既知である。
　測定制御部２１および演算部２２は、この動作を、被測定交流物理場（ここでは交流磁場）の複数周期ＰＥ１，・・・，ＰＥｍ（ｍ＞１）にわたって連続して繰り返し実行し、複数周期ＰＥ１，・・・，ＰＥｍ分の測定値［ＢＭ１（１），・・・，ＢＭｎ（１）］，・・・，［ＢＭ１（ｍ），・・・，ＢＭｎ（ｍ）］を特定する。
　そして、演算部２２は、１周期内の各位相ｉについて、複数周期分の測定値ＢＭｉ（１），・・・，ＢＭｉ（ｍ）の平均値ＢＭＡｉを次式に従って演算し、その演算結果［ＢＭＡ１，・・・，ＢＭＡｎ］を測定結果とする。
　ＢＭＡｉ＝（ＢＭｉ（１）＋・・・＋ＢＭｉ（ｍ））／ｍ
　ここでは、演算部２２は、さらに、その測定結果［ＢＭＡ１，・・・，ＢＭＡｎ］に対してカーブフィッティングを行い、フィッティング曲線（つまり、１周期分の測定波形）を導出する。
　ここで、当該フィッティング曲線の関数形（例えば正弦波など）が既知である場合には、例えば、上述の測定結果に基づいて、その関数形内のパラメータの最適値を最小２乗法などで導出して、フィッティング曲線の関数を得る。また、当該フィッティング曲線の関数形が未知である場合には、例えば、上述の測定結果に基づいて、所定の方法に従って補間（内挿や外挿）を行うことで、フィッティング曲線を得る。
　さらに、演算部２２は、上述のように導出したフィッティング曲線に基づいて、被測定交流物理場の最大振幅の測定結果を導出する。例えば図２に示すように、正弦波の被測定交流磁場の最大振幅Ｂａｃの測定値（推定値）が得られる。
　ここで、各直流物理場測定シーケンスＳＱｉとして実行されるラムゼイパルスシーケンスについて説明する。
　１回のラムゼイパルスシーケンスでは、例えば図３に示すように、測定制御部２１は、（ａ）まず、照射装置４で、所定波長の励起光を磁気共鳴部材１に照射して、磁気共鳴部材１の電子スピンの状態を揃え、（ｂ）その後、コイル２および高周波電源３で、第１のπ／２パルスのマイクロ波磁場を磁気共鳴部材１に印加し、（ｃ）その後、所定の時間間隔ｔｔ経過後、コイル２および高周波電源３で、第２のπ／２パルスのマイクロ波磁場を磁気共鳴部材１に印加し、（ｄ）その後、照射装置４で、射影測定用の測定光を照射するとともに、磁気共鳴部材１の発する蛍光を検出装置５で受光しその受光光量（検出光量）を検出する。
　なお、時間間隔ｔｔにおいては、外部磁場（ここでは被測定交流磁場）の磁束密度の時間積分に比例して電子スピンの向きが変化するため、検出光量から、外部磁場（ここでは被測定交流磁場）の磁束密度が導出可能となっている。
　さらに、π／２パルスの時間幅ｔｗは、電子スピンをπ／２だけ回転させる時間（数十ナノ秒程度）とされ、当該磁気共鳴部材１のラビ振動の周期などから予め特定される。また、例えば、励起光の照射時間および測定光の照射時間は、それぞれ数マイクロ秒から数十マイクロ秒程度である。また、例えば、上述の時間間隔ｔｔは、数百マイクロ秒以下とされる。また、励起光照射と第１のπ／２パルスとの時間間隔、および第２のπ／２パルスと測定光照射との時間間隔は、それぞれ、短いほど良い。
　この実施の形態では、例えば図３に示すように、被測定交流磁場の１周期の中で、ラムゼイパルスシーケンス（直流物理場測定シーケンスＳＱｉ）が連続的に実行される。なお、図３では、ある測定シーケンスＳＱｉにおける測定光の照射と次の測定シーケンスＳＱ（ｉ＋１）における励起光の照射とが別々に行われているが、両者をまとめて行ってもよい。
　ここで、被測定交流磁場のある周期ＰＥｊにおいて、ｉ番目の測定シーケンスＳＱｉは、ある位相Ｐｉ（未知の位相でもよい）で実行され、次の周期ＰＥ（ｊ＋１）においても、ｉ番目の測定シーケンスＳＱｉは、同一の位相Ｐｉで実行される。
　なお、上述のようにして得られた測定結果は、データとして図示せぬ記憶装置に記憶されたり、外部装置に送信されたり、表示装置で表示されたりするようにしてもよい。
　また、例えば、上述の１測定シーケンスＳＱｉの長さは、被測定交流物理場の半周期以下とされ、例えば、約１ｋＨｚ以下（特に、約１００Ｈｚ以下）の低周波数の交流物理場が測定対象とされる。
　なお、上述のスピンエコーパルスシーケンスの場合、電子スピンの横緩和時間Ｔ２が１ミリ秒程度であり、被測定交流磁場が約１ｋＨｚ以下の周波数成分を有するときには、通常、正確な磁場強度の測定が困難であるが、本実施の形態では、そのような場合で被測定交流磁場が約１ｋＨｚ以下の周波数成分を有していても、物理場の測定が良好に行われる。
　以上のように、上記実施の形態によれば、磁気共鳴部材１は、被測定交流物理場内に配置され所定の量子系での量子操作の可能な部材である。コイル２は、その磁気共鳴部材１にマイクロ波の磁場を印加する。高周波電源３は、そのコイル２にマイクロ波の電流を導通させる。直流物理場測定シーケンスにおいて、照射装置４は、その磁気共鳴部材１に光を照射し、検出装置５は、その磁気共鳴部材１から、被測定交流物理場に対応する物理的事象を検出する。そして、測定制御部２１は、直流物理場測定シーケンスを所定複数回実行し、直流物理場測定シーケンスのそれぞれにおいて、高周波電源３および照射装置４を制御し、検出装置５により検出された物理的事象の検出値を特定する。演算部２２は、複数回の直流物理場測定シーケンスに対応する検出値に基づいて、被測定交流物理場の特定期間分（ここでは、１周期分）の測定結果を演算する。
　これにより、磁気共鳴部材１における所定の量子系での量子操作を利用して、交流磁場などの長周期な物理場が正確に測定される。具体的には、上述のように、低周波領域での感度Ｈが良好に得られるため、測定対象の物理場の測定値が正確に得られる。
　また、上述のように、上述の複数回の直流物理場測定シーケンスが均一に実行される場合（直流物理場測定シーケンスにおける時間間隔ｔｔの長さおよび直流物理場測定シーケンスの時間間隔が被測定交流物理場の位相に拘わらず一定と設定される場合）、複数回の直流物理場測定シーケンスを被測定交流物理場に同期させずに、上述の複数回の直流物理場測定シーケンスを行うようにしてもよい。なお、スピンエコーパルスシーケンスの場合、被測定交流磁場との同期が必要となる。
　なお、上述の実施の形態に対する様々な変更および修正については、当業者には明らかである。そのような変更および修正は、その主題の趣旨および範囲から離れることなく、かつ、意図された利点を弱めることなく行われてもよい。つまり、そのような変更および修正が請求の範囲に含まれることを意図している。
　例えば、上記実施の形態では、一例として、被測定交流磁場の波形を正弦波と仮定しているが、その代わりに、被測定交流磁場の波形は他の波形（三角波、のこぎり波、矩形波、複数の波形を合成したものなど）であってもよく、被測定交流磁場の周期が既知であればよい。また、被測定交流磁場は、直流成分を有していてもよい。さらに、交流磁場以外の物理場が測定対象である場合でも同様である。
　また、上記実施の形態において、被測定交流物理場の１周期内の位相に応じて、直流物理場測定シーケンスにおける時間間隔ｔｔを変化させてもよい。その場合、被測定交流物理場と同期していることが既知である外部交流信号があれば、その外部交流信号の位相に基づいて被測定交流物理場の位相を特定してもよい。例えば、外部交流信号に起因する物理現象によって被測定交流物理場が発生している場合、被測定交流物理場は、外部交流信号に同期していることがある。あるいは、例えば、時間間隔ｔｔを一定にして上述の直流物理場測定シーケンスを繰り返し実行し、複数回の直流物理場測定シーケンスで得られた測定結果から、被測定交流物理場の波形などを特定してその位相を特定するようにしてもよい。
　さらに、上記実施の形態では、特定期間の一例として、１周期分の測定結果を得ているが、１周期分より長い周期（例えば２周期分など）の測定結果を得るようにしてもよい。さらに、上記実施の形態において、１周期より短い期間（つまり、特定の位相範囲）の測定結果を得るようにしてもよい。その場合、例えば、上述の複数周期における各周期において、１周期ごとにその期間に対応する検出値だけ取得してその期間に対応する測定値だけを導出するようにしてもよい。ただし、特定期間を１周期とするのが好ましい。
　さらに、上記実施の形態において、被測定物理場の強度（振幅）が大きく、上述の時間間隔ｔｔの期間において磁気共鳴部材１におけるカラーセンタの電子スピンが１周（２π）以上回転してしまう場合、その電子スピンの実際の回転角を２πで除算したときの余りの角度に対応する検出値、測定値および測定結果が得られる。そのため、その場合には、演算部２２は、上述の検出値、測定値、または測定結果から電子スピンの実際の回転角を推定して、その検出値、測定値、または測定結果を、電子スピンの実際の回転角に対応する値へ変換し、変換後の値に基づいて、上述のカーブフィッティングを行うようにしてもよい。なお、その際、当該変換を行わずに、直接的に、２πスピン位相周期を考慮したフィッティング関数を使用して、上述のカーブフィッティングを行うようにしてもよい。
　さらに、上記の実施の形態では、ある位相ｉについての複数（ｍ）周期分の測定値ＢＭｉ（１），・・・，ＢＭｉ（ｍ）のすべてが導出される場合、位相ｉについてのｍ個の測定値ＢＭｉ（１），・・・，ＢＭｉ（ｍ）の平均値ＢＭｉが、測定結果のうちのその位相ｉの部分とされているが、ある位相ｉについての複数（ｍ）周期分の測定値ＢＭｉ（１），・・・，ＢＭｉ（ｍ）の一部がなく、ＢＭｉ（１），・・・，ＢＭｉ（ｍ）の残りの部分のみが導出される場合、その残りの部分の平均値が、上述の平均値ＢＭｉとして計算される。つまり、例えば、ｍ個の測定値ＢＭｉ（１），・・・，ＢＭｉ（Ｌ－１），ＢＭｉ（Ｌ），ＢＭｉ（Ｌ＋１），・・・，ＢＭｉ（ｍ）のうち（１＜Ｌ－１，Ｌ，Ｌ＋１＜ｍ）、Ｌ番目（Ｌ番目の周期）の測定値ＢＭｉが得られない場合（例えば、対応する検出値が適切な値として得られない場合、Ｌ番目の測定が行われていない場合など）には、残りの（ｍ－１）個の測定値の平均値が、次式のようにして、上述の平均値ＢＭｉとして計算される。なお、２個以上の測定値が得られない場合も同様である。
　ＢＭＡｉ＝（ＢＭｉ（１）＋・・・＋ＢＭｉ（Ｌ－１）＋ＢＭｉ（Ｌ＋１）＋・・・＋ＢＭｉ（ｍ））／（ｍ－１）
　さらに、上記実施の形態では、複数周期分の、特定期間の測定値を位相ごとに平均化して測定結果を得ているが、時間的平均化が必要ない場合には、１周期分の、特定期間の測定値を（平均化せずに）測定結果としてもよい。
　本発明は、例えば、磁気測定装置および磁気測定方法に適用可能である。

Claims (4)

  1. 　被測定交流物理場内に配置され所定の量子系での量子操作の可能な磁気共鳴部材と、
    　前記磁気共鳴部材にマイクロ波の磁場を印加するコイルと、
    　前記コイルに前記マイクロ波の電流を導通させる高周波電源と、
    　直流物理場測定シーケンスにおいて、前記磁気共鳴部材に光を照射する照射装置と、
    　前記直流物理場測定シーケンスにおいて、前記磁気共鳴部材から、前記被測定交流物理場に対応する物理的事象を検出する検出装置と、
    　前記直流物理場測定シーケンスを所定複数回実行し、前記直流物理場測定シーケンスのそれぞれにおいて、前記高周波電源および前記照射装置を制御し、前記検出装置により検出された前記物理的事象の検出値を特定する測定制御部と、
    　前記複数回の前記直流物理場測定シーケンスに対応する前記検出値に基づいて、前記被測定交流物理場の特定期間分の測定結果を演算する演算部と、
    　を備えることを特徴とする測定装置。
  2. 　前記演算部は、前記被測定交流物理場の１周期ごとに前記直流物理場測定シーケンスを前記所定複数回実行し、前記被測定交流物理場の複数周期について得られた、前記複数周期と同数の前記特定期間分の測定結果の平均を演算することで、前記特定期間分の測定結果におけるノイズを減衰させることを特徴とする請求項１記載の測定装置。
  3. 　前記磁気共鳴部材は、電子スピン量子操作の可能なカラーセンタを含み、
    　前記直流物理場測定シーケンスは、前記マイクロ波として２つのπ／２パルスを含み、
    　前記物理的事象は、ラビ振動において、前記２つのπ／２パルスの間の時間間隔における自由歳差運動での前記カラーセンタの電子スピンの位相変化に対応し、
    　前記２つのπ／２パルスの間の時間間隔は、（ａ）前記磁気共鳴部材の実効横緩和時間に応じて設定されるとともに、（ｂ）前記被測定交流物理場の周波数が当該時間間隔によって得られる有効感度周波数範囲に含まれるように設定されること、
    　を特徴とする請求項１記載の測定装置。
  4. 　（ａ）直流物理場測定シーケンスを所定複数回実行し、前記直流物理場測定シーケンスのそれぞれにおいて、被測定交流物理場に対応する物理的事象を検出して前記物理的事象の検出値を特定し、
    　（ｂ）前記複数回の前記直流物理場測定シーケンスに対応する前記検出値に基づいて、前記被測定交流物理場の特定期間分の測定結果を演算し、
    　前記直流物理場測定シーケンスは、前記被測定交流物理場内に配置され所定の量子系での量子操作の可能な磁気共鳴部材と、前記磁気共鳴部材にマイクロ波の磁場を印加するコイルと、前記コイルに前記マイクロ波の電流を導通させる高周波電源と、前記磁気共鳴部材に光を照射する照射装置と、前記磁気共鳴部材から、前記被測定交流物理場に対応する前記物理的事象を検出する検出装置とを使用して実行されること、
    　を特徴とする測定方法。
PCT/JP2020/047122 2019-12-24 2020-12-17 測定装置および測定方法 Ceased WO2021132009A1 (ja)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202080083340.9A CN114746764B (zh) 2019-12-24 2020-12-17 测量装置和测量方法
US17/773,530 US11808831B2 (en) 2019-12-24 2020-12-17 Measurement apparatus and measurement method
EP20904834.7A EP4083647B1 (en) 2019-12-24 2020-12-17 Measurement device and measurement method

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
JP2019233583A JP7387107B2 (ja) 2019-12-24 2019-12-24 測定装置および測定方法
JP2019-233583 2019-12-24

Publications (1)

Publication Number Publication Date
WO2021132009A1 true WO2021132009A1 (ja) 2021-07-01

Family

ID=76576084

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/JP2020/047122 Ceased WO2021132009A1 (ja) 2019-12-24 2020-12-17 測定装置および測定方法

Country Status (5)

Country Link
US (1) US11808831B2 (ja)
EP (1) EP4083647B1 (ja)
JP (1) JP7387107B2 (ja)
CN (1) CN114746764B (ja)
WO (1) WO2021132009A1 (ja)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP7194327B2 (ja) * 2018-07-03 2022-12-22 スミダコーポレーション株式会社 磁場測定装置および磁場測定方法
JP7403272B2 (ja) * 2019-10-11 2023-12-22 Ｔｉａｎｍａ Ｊａｐａｎ株式会社 磁気光学式計測装置
DE102020123993A1 (de) * 2020-09-15 2022-03-17 Endress+Hauser SE+Co. KG Quantensensor
JP7616569B2 (ja) * 2021-08-03 2025-01-17 スミダコーポレーション株式会社 測定装置および測定方法
JP7749184B2 (ja) * 2021-11-19 2025-10-06 スミダコーポレーション株式会社 デジタイズ装置およびデジタイズ方法
US20260017556A1 (en) * 2022-08-05 2026-01-15 Okinawa Institute Of Science And Technology School Corporation Quantum error correction system and process for quantum error correction
WO2024203140A1 (ja) * 2023-03-29 2024-10-03 京セラ株式会社 磁気検出器

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001178701A (ja) * 1999-12-24 2001-07-03 Toshiba Medical System Co Ltd Ｍｒｉ装置及びｍｒイメージング方法
JP2011180570A (ja) * 2010-02-02 2011-09-15 Kyoto Univ 蛍光顕微鏡装置
JP2012110489A (ja) 2010-11-24 2012-06-14 Sumitomo Electric Ind Ltd 磁気計測装置、及び、磁気センサ素子の製造方法
JP2016205954A (ja) * 2015-04-21 2016-12-08 ルネサスエレクトロニクス株式会社 磁気計測装置
US20170343695A1 (en) * 2016-05-31 2017-11-30 Lockheed Martin Corporation Magneto-Optical Detecting Apparatus and Methods
US20180136291A1 (en) * 2016-11-08 2018-05-17 Linh M. Pham Methods and apparatus for optically detecting magnetic resonance
JP2020008298A (ja) * 2018-07-03 2020-01-16 スミダコーポレーション株式会社 磁場測定装置および磁場測定方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2016118756A1 (en) * 2015-01-23 2016-07-28 Lockheed Martin Corporation Apparatus and method for high sensitivity magnetometry measurement and signal processing in a magnetic detection system
WO2017127096A1 (en) * 2016-01-21 2017-07-27 Lockheed Martin Corporation Diamond nitrogen vacancy sensor with dual rf sources
US10228429B2 (en) * 2017-03-24 2019-03-12 Lockheed Martin Corporation Apparatus and method for resonance magneto-optical defect center material pulsed mode referencing
US10408890B2 (en) * 2017-03-24 2019-09-10 Lockheed Martin Corporation Pulsed RF methods for optimization of CW measurements
US10459041B2 (en) * 2017-03-24 2019-10-29 Lockheed Martin Corporation Magnetic detection system with highly integrated diamond nitrogen vacancy sensor
US10838021B2 (en) * 2017-07-11 2020-11-17 Lockheed Martin Corporation Apparatus and method for simultaneous ramsey vector magnetometry
CN108646203B (zh) * 2018-05-04 2020-05-15 中国科学技术大学 一种纳米尺度的微波磁场测量方法
JP7209176B2 (ja) * 2018-10-26 2023-01-20 スミダコーポレーション株式会社 磁場発生源検出装置および磁場発生源検出方法
JP7265699B2 (ja) * 2019-02-22 2023-04-27 スミダコーポレーション株式会社 測定装置
CN110095740B (zh) * 2019-05-15 2020-07-24 中国科学院地质与地球物理研究所 电子自旋磁场测量方法及系统

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001178701A (ja) * 1999-12-24 2001-07-03 Toshiba Medical System Co Ltd Ｍｒｉ装置及びｍｒイメージング方法
JP2011180570A (ja) * 2010-02-02 2011-09-15 Kyoto Univ 蛍光顕微鏡装置
JP2012110489A (ja) 2010-11-24 2012-06-14 Sumitomo Electric Ind Ltd 磁気計測装置、及び、磁気センサ素子の製造方法
JP2016205954A (ja) * 2015-04-21 2016-12-08 ルネサスエレクトロニクス株式会社 磁気計測装置
US20170343695A1 (en) * 2016-05-31 2017-11-30 Lockheed Martin Corporation Magneto-Optical Detecting Apparatus and Methods
US20180136291A1 (en) * 2016-11-08 2018-05-17 Linh M. Pham Methods and apparatus for optically detecting magnetic resonance
JP2020008298A (ja) * 2018-07-03 2020-01-16 スミダコーポレーション株式会社 磁場測定装置および磁場測定方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See also references of EP4083647A4

Also Published As

Publication number Publication date
EP4083647B1 (en) 2025-06-18
JP7387107B2 (ja) 2023-11-28
JP2021103093A (ja) 2021-07-15
EP4083647A1 (en) 2022-11-02
US11808831B2 (en) 2023-11-07
CN114746764B (zh) 2025-03-25
EP4083647A4 (en) 2023-12-27
US20220373628A1 (en) 2022-11-24
CN114746764A (zh) 2022-07-12

Similar Documents

Publication Publication Date Title
JP7387107B2 (ja) 測定装置および測定方法
JP7194327B2 (ja) 磁場測定装置および磁場測定方法
Sheng et al. Emission of electromagnetic pulses from laser wakefields through linear mode conversion
US20180275224A1 (en) Generation of magnetic field proxy through rf frequency dithering
WO2021054141A1 (ja) 測定装置および測定方法
JP7465467B2 (ja) 磁場測定装置および磁場測定方法
JP7369344B2 (ja) 磁場測定装置および磁場測定方法
US10036797B2 (en) Electron spin resonance apparatus
JP7114406B2 (ja) 電子スピン共鳴測定装置及び方法
JP7265699B2 (ja) 測定装置
JP6454241B2 (ja) 電子スピン共鳴装置
JP2022161376A (ja) 電子スピン共鳴測定装置
CN117907908A (zh) 一种磁场强度测量方法及装置
Romalis et al. Vapor cell Rydberg atom electrometry with time-separated fields
JP5256136B2 (ja) 電磁波測定装置、及び電磁波測定方法
JP7812074B2 (ja) 測定装置および測定方法
JP7025010B2 (ja) 歪み検出装置、歪み検出方法及び歪み検出プログラム
JP7007648B2 (ja) 自発分極検出装置、自発分極検出方法及び自発分極検出プログラム
JP2025149785A (ja) 測定装置および測定方法
CN116736197A (zh) 一种切片锁频的快速量子调控磁测量方法及系统
Lisin et al. Measurement of g-factors of ground and excited optical states of Er3+ in YLiF4 and LuLiF4 in zero dc magnetic field by photon echo method
Masuyama et al. Highly sensitive macro-scale diamond magnetometer operated by dynamical decoupling sequence with coplanar waveguide resonator
CN117741529A (zh) 一种基于双点测磁共振谱拟合的量子测量方法
WO2018174911A1 (en) Generation of magnetic field proxy through rf frequency dithering
CN117740760A (zh) 一种空气电离产生剩余电流的太赫兹波形探测系统

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 20904834

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

ENP Entry into the national phase

Ref document number: 2020904834

Country of ref document: EP

Effective date: 20220725

WWG Wipo information: grant in national office

Ref document number: 202080083340.9

Country of ref document: CN

WWG Wipo information: grant in national office

Ref document number: 2020904834

Country of ref document: EP