CN115479921B - Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment - Google Patents
Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment Download PDFInfo
- Publication number
- CN115479921B CN115479921B CN202211156899.7A CN202211156899A CN115479921B CN 115479921 B CN115479921 B CN 115479921B CN 202211156899 A CN202211156899 A CN 202211156899A CN 115479921 B CN115479921 B CN 115479921B
- Authority
- CN
- China
- Prior art keywords
- wavelength
- oscillation
- light
- detection
- brillouin scattering
- 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.)
- Active
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 98
- 238000002474 experimental method Methods 0.000 title claims abstract description 55
- 238000000034 method Methods 0.000 title claims abstract description 32
- 230000010355 oscillation Effects 0.000 claims abstract description 90
- 230000008859 change Effects 0.000 claims description 24
- 238000004088 simulation Methods 0.000 claims description 24
- 239000000463 material Substances 0.000 claims description 23
- 238000002310 reflectometry Methods 0.000 claims description 20
- 239000000523 sample Substances 0.000 claims description 18
- 238000010521 absorption reaction Methods 0.000 claims description 13
- 238000010586 diagram Methods 0.000 claims description 10
- 230000007423 decrease Effects 0.000 claims description 9
- 230000004907 flux Effects 0.000 claims description 8
- 230000008878 coupling Effects 0.000 claims description 6
- 238000010168 coupling process Methods 0.000 claims description 6
- 238000005859 coupling reaction Methods 0.000 claims description 6
- 238000012790 confirmation Methods 0.000 claims description 3
- 238000005094 computer simulation Methods 0.000 abstract description 5
- 230000000694 effects Effects 0.000 abstract description 3
- 238000005516 engineering process Methods 0.000 description 5
- 230000008569 process Effects 0.000 description 5
- 230000003287 optical effect Effects 0.000 description 3
- 101001121408 Homo sapiens L-amino-acid oxidase Proteins 0.000 description 2
- 102100026388 L-amino-acid oxidase Human genes 0.000 description 2
- 101100233916 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) KAR5 gene Proteins 0.000 description 2
- 230000005535 acoustic phonon Effects 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000005693 optoelectronics Effects 0.000 description 2
- 101100012902 Saccharomyces cerevisiae (strain ATCC 204508 / S288c) FIG2 gene Proteins 0.000 description 1
- 239000012472 biological sample Substances 0.000 description 1
- 230000001427 coherent effect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000002105 nanoparticle Substances 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000005086 pumping Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- 238000011896 sensitive detection Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N21/00—Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
- G01N21/17—Systems in which incident light is modified in accordance with the properties of the material investigated
- G01N21/47—Scattering, i.e. diffuse reflection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/12—Timing analysis or timing optimisation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明公开了一种时域布里渊散射实验中优选探测光波长的方法及装置,利用探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论,以此进行计算机仿真实验来优选出探测光波长。本发明优选出时域布里渊散射实验中的探测光波长,可以获得更优的利用时域布里渊散射实验进行检测分析的效果。
The present invention discloses a method and device for optimizing the wavelength of a detection light in a time-domain Brillouin scattering experiment. The method uses the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the wavelength of the detection light changes, and uses this to perform a computer simulation experiment to optimize the wavelength of the detection light. The present invention optimizes the wavelength of the detection light in the time-domain Brillouin scattering experiment, and can obtain a better effect of detection and analysis using the time-domain Brillouin scattering experiment.
Description
技术领域Technical Field
本申请涉及检测分析技术领域,特别是涉及一种时域布里渊散射实验中优选探测光波长的方法及装置。The present application relates to the field of detection and analysis technology, and in particular to a method and device for optimizing the wavelength of detection light in a time-domain Brillouin scattering experiment.
背景技术Background Art
时域布里渊散射技术(Time-domainBrillouin scattering,TDBS)是一种高灵敏度的检测分析技术,其原理是用一束泵浦脉冲激发样品中原子或分子的运动,再用一束可以调节时间差的探测脉冲对某一时刻样品内原子或分子运动变化进行探测,即对声子的探测,在实验上测得的表现形式为样品表面的反射率变化。根据材料的不同、探测光及泵浦光的波长和功率等变化,所测量得到的反射率变化也会因此而改变,其变化一般由一段由电子涨落引起的曲线和一段阻尼振荡组成,其中的阻尼振荡正是研究的核心,包含三个重要的参数:振荡的振幅、振荡的频率、振荡的衰减系数。Time-domain Brillouin scattering (TDBS) is a highly sensitive detection and analysis technology. Its principle is to use a pump pulse to stimulate the movement of atoms or molecules in the sample, and then use a probe pulse with adjustable time difference to detect the change of atomic or molecular movement in the sample at a certain moment, that is, the detection of phonons. The experimentally measured manifestation is the change of reflectivity on the surface of the sample. Depending on the material, the wavelength and power of the probe light and the pump light, the measured reflectivity change will also change. The change is generally composed of a curve caused by electronic fluctuations and a damped oscillation. The damped oscillation is the core of the research, which contains three important parameters: the amplitude of the oscillation, the frequency of the oscillation, and the attenuation coefficient of the oscillation.
TDBS技术具有无损、无接触、灵敏度高等优点,通常能够探测到10-5量级的反射率变化。采用TDBS技术,能够将时域信号转换为频域信号用于测量样品声速(折射率已知的情况下);通过追踪信号的振幅衰减,可以计算出样品声波能量的耗散速率;对纳米尺度的粒子、生物样品等结构进行无创成像;同时TDBS技术还能在电子带间跃迁的研究中发挥作用。TDBS technology has the advantages of being non-destructive, contactless, and highly sensitive, and can usually detect reflectivity changes of the order of 10 -5 . Using TDBS technology, time domain signals can be converted into frequency domain signals for measuring the sound velocity of the sample (when the refractive index is known); by tracking the amplitude attenuation of the signal, the dissipation rate of the sample sound wave energy can be calculated; nanoscale particles, biological samples and other structures can be non-invasively imaged; at the same time, TDBS technology can also play a role in the study of inter-band transitions of electrons.
声子分为光学支和声学支,声学支具有以下的特点:1、振荡频率会随着探测光波长的不同而变化,即布里渊频率其中λpr是探测光波长,npr是在该探测波长下样品的折射率,v是样品的声速;2、声学声子的振荡频率通常只有亚GHz的量级,远小于光学声子的THz的量级;3、通常能够在100ps的时间范围内观测到相干声学声子振荡,而光学声子的可观测范围则小得多。例如在Ge的TDBS实验中,通常测量到的是声学支(CAP),TDBS实验方法可以探测大部分样品中的CAP,在实验中Ge只是作为一个常见的半导体系统进行展示。由于实验中使用的光波长通常在可见光的范围,因此毫米量级的Ge材料可以视为大块材料。Phonons are divided into optical and acoustic branches. The acoustic branch has the following characteristics: 1. The oscillation frequency varies with the wavelength of the detection light, that is, the Brillouin frequency. Where λ pr is the wavelength of the probe light, n pr is the refractive index of the sample at the probe wavelength, and v is the speed of sound of the sample; 2. The oscillation frequency of acoustic phonons is usually only in the sub-GHz range, which is much smaller than the THz range of optical phonons; 3. Coherent acoustic phonon oscillations can usually be observed within a time range of 100ps, while the observable range of optical phonons is much smaller. For example, in the TDBS experiment of Ge, the acoustic branch (CAP) is usually measured. The TDBS experimental method can detect the CAP in most samples. In the experiment, Ge is only displayed as a common semiconductor system. Since the wavelength of light used in the experiment is usually in the visible light range, millimeter-scale Ge materials can be regarded as bulk materials.
在TDBS实验中,为了更高效、更准确地从实验数据中提取有效信息,通常希望振荡信号的振幅足够大,能观测到的信号足够多,同时由于实验中光电器件的分辨率有限,频率越小的振荡分析起来越困难。因此,如何找到一个方法,能够使这三个因素同时满足要求就显得非常重要。In TDBS experiments, in order to extract effective information from experimental data more efficiently and accurately, it is usually hoped that the amplitude of the oscillation signal is large enough and that enough signals can be observed. At the same time, due to the limited resolution of the optoelectronic devices in the experiment, the smaller the frequency of the oscillation, the more difficult it is to analyze. Therefore, it is very important to find a method that can meet these three requirements at the same time.
发明内容Summary of the invention
本发明针对上述问题,提供了一种时域布里渊散射实验中优选探测光波长的方法及装置,利用探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论,以此进行计算机仿真实验来优选出探测光波长。In view of the above problems, the present invention provides a method and device for optimizing the wavelength of detection light in time-domain Brillouin scattering experiments. The method and device use the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the wavelength of the detection light changes, and thus a computer simulation experiment is performed to optimize the wavelength of the detection light.
本发明的第一方面,一种时域布里渊散射实验中优选探测光波长的方法,包括:In a first aspect of the present invention, a method for optimizing the wavelength of a detection light in a time-domain Brillouin scattering experiment comprises:
基于时域布里渊散射实验中反射率的变化公式,得到当探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论;Based on the formula for the change of reflectivity in the time-domain Brillouin scattering experiment, it is concluded that when the wavelength of the probe light changes, the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated;
通过仿真实验来寻找优选的探测波长,具体包括:The optimal detection wavelength is found through simulation experiments, including:
对于一个材料系统,根据所述材料系统的参数进行仿真:For a material system, simulation is performed according to the parameters of the material system:
针对振荡持续时间的考量,利用一个参数τ=1/(αpu×v)来量化震荡的持续时间,得到探测波长不同时会引起参数τ变化的结论,其中αpu是泵浦光的吸收系数,v是介质的声速;In consideration of the duration of oscillation, a parameter τ = 1/(α pu × v) is used to quantify the duration of oscillation, and it is concluded that different detection wavelengths will cause the parameter τ to change, where α pu is the absorption coefficient of the pump light and v is the sound velocity of the medium;
针对振荡振幅的考量,找到振幅最大值时的探测波长;Considering the oscillation amplitude, find the detection wavelength when the amplitude is maximum;
针对振荡频率的考量,找到振荡频率随探测波长的增加而下降的规律;Considering the oscillation frequency, we found that the oscillation frequency decreases as the detection wavelength increases;
将振荡信号的振幅A0、频率f和振荡持续时间τ三者做乘积并绘图展现出来,通过图示找出优选的探测波长。The amplitude A 0 , frequency f and oscillation duration τ of the oscillation signal are multiplied and plotted, and the optimal detection wavelength is found through the plot.
本发明的进一步技术方案为:所述方法还包括对找出的优选探测波长进行多次仿真实验进行验证。A further technical solution of the present invention is: the method further comprises performing multiple simulation experiments to verify the found optimal detection wavelength.
本发明的进一步技术方案为:所述基于时域布里渊散射实验中反射率的变化公式,反射率的具体变化公式为:A further technical solution of the present invention is: the specific change formula of the reflectivity in the time-domain Brillouin scattering experiment is:
其中,i表示虚数,k=2π/λpr表示探测光波数,λpr表示探测光波长,表示介质介电常数对光子能量的导数,v表示介质的声速;acv表示相对形变势耦合常数,npr表示介质对于探测光的复折射率,αpu表示泵浦光的吸收系数,Rpu表示介质对于泵浦光的反射率,F表示泵浦光通量,ρ表示介质的密度,t表示时间,c是光速,表示约化普朗克常数;ωpu表示pump光的角频率。in, i represents an imaginary number, k = 2π/λ pr represents the detection light wave number, λ pr represents the detection light wavelength, represents the derivative of the dielectric constant of the medium with respect to the photon energy, v represents the sound velocity of the medium; a cv represents the relative deformation potential coupling constant, n pr represents the complex refractive index of the medium for the probe light, α pu represents the absorption coefficient of the pump light, R pu represents the reflectivity of the medium for the pump light, F represents the pump light flux, ρ represents the density of the medium, t represents the time, c is the speed of light, represents the reduced Planck constant; ω pu represents the angular frequency of the pump light.
本发明的进一步技术方案为:所述材料系统的参数包括介质介电常数对光子能量的导数、介质声速、相对形变势耦合常数、泵浦光的吸收系数、介质对于泵浦光的反射率、介质的密度、介质对于探测光的复折射率。A further technical solution of the present invention is: the parameters of the material system include the derivative of the dielectric constant of the medium with respect to photon energy, the sound velocity of the medium, the relative deformation potential coupling constant, the absorption coefficient of the pump light, the reflectivity of the medium to the pump light, the density of the medium, and the complex refractive index of the medium to the detection light.
本发明的进一步技术方案为:当所述材料系统为Ge时,优选的时域布里渊散射实验探测波长范围为600-820nm。A further technical solution of the present invention is: when the material system is Ge, the preferred detection wavelength range of the time-domain Brillouin scattering experiment is 600-820nm.
本发明的第二方面,提供了一种时域布里渊散射实验中优选探测光波长的装置,包括:A second aspect of the present invention provides a device for optimally detecting light wavelength in a time-domain Brillouin scattering experiment, comprising:
探测光波长关系确认模块,用于基于时域布里渊散射实验中反射率的变化公式,得到当探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论;The detection light wavelength relationship confirmation module is used to obtain the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the detection light wavelength changes based on the reflectivity change formula in the time-domain Brillouin scattering experiment;
仿真实验模块,用于通过仿真实验来寻找优选的探测波长,具体包括:The simulation experiment module is used to find the optimal detection wavelength through simulation experiments, including:
对于一个材料系统,根据所述材料系统的参数进行仿真:For a material system, simulation is performed according to the parameters of the material system:
针对振荡持续时间的考量,利用一个参数τ=1/(αpu×v)来量化震荡的持续时间,得到探测波长不同时会引起参数τ变化的结论,其中αpu是泵浦光的吸收系数,v是介质的声速;In consideration of the duration of oscillation, a parameter τ = 1/(α pu × v) is used to quantify the duration of oscillation, and it is concluded that different detection wavelengths will cause the parameter τ to change, where α pu is the absorption coefficient of the pump light and v is the sound velocity of the medium;
针对振荡振幅的考量,找到振幅最大值时的探测波长;Considering the oscillation amplitude, find the detection wavelength when the amplitude is maximum;
针对振荡频率的考量,找到振荡频率随探测波长的增加而下降的规律;Considering the oscillation frequency, we found that the oscillation frequency decreases as the detection wavelength increases;
将振荡信号的振幅A0、频率f和振荡持续时间τ三者做乘积并绘图展现出来,通过图示找出优选的探测波长。The amplitude A 0 , frequency f and oscillation duration τ of the oscillation signal are multiplied and plotted, and the optimal detection wavelength is found through the plot.
本发明提供的一种时域布里渊散射实验中优选探测光波长的方法及装置,利用探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论,以此进行计算机仿真实验来优选出探测光波长。本发明优选出时域布里渊散射实验中的探测光波长,可以获得更优的利用时域布里渊散射实验进行检测分析的效果。The present invention provides a method and device for optimizing the wavelength of a detection light in a time-domain Brillouin scattering experiment. The method and device utilize the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the wavelength of the detection light changes, and thus a computer simulation experiment is performed to optimize the wavelength of the detection light. The present invention optimizes the wavelength of the detection light in the time-domain Brillouin scattering experiment, and can obtain a better effect of detection and analysis using the time-domain Brillouin scattering experiment.
附图说明BRIEF DESCRIPTION OF THE DRAWINGS
图1(a)是本发明实施例中泵浦波长为800nm,泵浦光通量为0.3mJ/cm2,探测光波长分别在539nm,652nm及1127nm下的TDBS仿真结果图;FIG. 1( a ) is a diagram showing the TDBS simulation results in an embodiment of the present invention when the pump wavelength is 800 nm, the pump light flux is 0.3 mJ/cm 2 , and the probe light wavelengths are 539 nm, 652 nm, and 1127 nm respectively;
图1(b)是本发明实施例中提取阻尼振荡振幅,其中探测光子能量在1-4eV范围内得到的绘制图;FIG1( b ) is a plot obtained by extracting the damped oscillation amplitude in an embodiment of the present invention, wherein the detected photon energy is in the range of 1-4 eV;
图1(c)是本发明实施例中可观测到振荡信号的时间范围以及振荡信号频率的仿真结果图;FIG1( c ) is a simulation result diagram of the time range in which an oscillation signal can be observed and the frequency of the oscillation signal in an embodiment of the present invention;
图1(d)是本发明实施例中振荡信号的振幅、频率与观测范围三者的乘积随探测波长的变化效果图;FIG. 1( d ) is a diagram showing the variation of the product of the amplitude, frequency and observation range of the oscillation signal with the detection wavelength in an embodiment of the present invention;
图2(a)是本发明另一实施例中泵浦波长为800nm,泵浦光通量为0.3mJ/cm2,探测光子能量分别在625,730,870,920nm下的TDBS实验结果图;FIG. 2( a ) is a diagram showing the TDBS experimental results in another embodiment of the present invention, with a pump wavelength of 800 nm, a pump light flux of 0.3 mJ/cm 2 , and detection photon energies of 625, 730, 870, and 920 nm;
图2(b)是本发明另一实施例中振荡信号的振幅、频率与观测范围三者的乘积随探测波长的变化效果图;FIG2( b ) is a diagram showing the variation of the product of the amplitude, frequency and observation range of the oscillation signal with the detection wavelength in another embodiment of the present invention;
图3是本发明实施例中时域布里渊散射实验中优选探测光波长的装置结构示意图。FIG. 3 is a schematic diagram of the structure of a device for optimally detecting light wavelength in a time-domain Brillouin scattering experiment according to an embodiment of the present invention.
具体实施方式DETAILED DESCRIPTION
下面结合附图和实施例对本发明作进一步的详细说明。可以理解的是,此处所描述的具体实施例仅仅用于解释本发明,而非对本发明的限定。另外还需要说明的是,为了便于描述,附图中仅出示了与本发明相关的部分而非全部结构。The present invention will be further described in detail below in conjunction with the accompanying drawings and embodiments. It is to be understood that the specific embodiments described herein are only used to explain the present invention, rather than to limit the present invention. It should also be noted that, for ease of description, only parts related to the present invention, rather than all structures, are shown in the accompanying drawings.
在更加详细地讨论示例性实施例之前应当提到的是,一些示例性实施例被描述成作为流程图描绘的处理或方法。虽然流程图将各步骤描述成顺序的处理,但是其中的许多步骤可以被并行地、并发地或者同时实施。此外,各步骤的顺序可以被重新安排。当其操作完成时所述处理可以被终止,但是还可以具有未包括在附图中的附加步骤。所述处理可以对应于方法、函数、规程、子例程、子程序等等。It should be mentioned before discussing the exemplary embodiments in more detail that some exemplary embodiments are described as processes or methods depicted as flow charts. Although the flow charts describe the steps as sequential processes, many of the steps therein can be implemented in parallel, concurrently or simultaneously. In addition, the order of the steps can be rearranged. The process can be terminated when its operation is completed, but can also have additional steps not included in the accompanying drawings. The process can correspond to a method, function, procedure, subroutine, subprogram, etc.
本发明实施例针对一种时域布里渊散射实验中优选探测光波长的方法及装置,提供了如下实施例:The embodiment of the present invention provides the following embodiments for a method and device for optimizing the wavelength of detection light in a time-domain Brillouin scattering experiment:
基于本发明的实施例1Embodiment 1 based on the present invention
本发明实施例1的一种时域布里渊散射实验中优选探测光波长的方法,包括:A method for optimizing the wavelength of a detection light in a time-domain Brillouin scattering experiment according to Embodiment 1 of the present invention comprises:
基于时域布里渊散射实验中反射率的变化公式,得到当探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论;Based on the formula for the change of reflectivity in the time-domain Brillouin scattering experiment, it is concluded that when the wavelength of the probe light changes, the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated;
具体地,TDBS实验中的阻尼振荡信号受到多维度的影响,诸如探测波长、材料复折射率随入射光波长的变化、材料介电常数随入射光子能量的导数等等。反射率的变化可以表达为:Specifically, the damped oscillation signal in the TDBS experiment is affected by multiple dimensions, such as the detection wavelength, the change of the complex refractive index of the material with the wavelength of the incident light, the derivative of the dielectric constant of the material with the energy of the incident photon, etc. The change of reflectivity can be expressed as:
其中,i表示虚数,k=2π/λpr表示探测光波数,λpr表示探测光波长,表示介质介电常数对光子能量的导数,v表示介质的声速;acv表示相对形变势耦合常数,npr表示介质对于探测光的复折射率,αpu表示泵浦光的吸收系数,Rpu表示介质对于泵浦光的反射率,F表示泵浦光通量,ρ表示介质的密度,t表示时间,c是光速,表示约化普朗克常数;ωpu表示pump光的角频率。从公式中可以看出,当探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制,因此可以通过计算机仿真来寻找最佳的探测波长。in, i represents an imaginary number, k = 2π/λ pr represents the detection light wave number, λ pr represents the detection light wavelength, represents the derivative of the dielectric constant of the medium with respect to the photon energy, v represents the sound velocity of the medium; a cv represents the relative deformation potential coupling constant, n pr represents the complex refractive index of the medium for the probe light, α pu represents the absorption coefficient of the pump light, R pu represents the reflectivity of the medium for the pump light, F represents the pump light flux, ρ represents the density of the medium, t represents the time, c is the speed of light, represents the reduced Planck constant; ω pu represents the angular frequency of the pump light. It can be seen from the formula that when the wavelength of the detection light changes, the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated, so the optimal detection wavelength can be found through computer simulation.
通过仿真实验来寻找优选的探测波长,具体包括:The optimal detection wavelength is found through simulation experiments, including:
对于一个材料系统,根据所述材料系统的参数进行仿真:For a material system, simulation is performed according to the parameters of the material system:
针对振荡持续时间的考量,利用一个参数τ=1/(αpu×v)来量化震荡的持续时间,得到探测波长不同时会引起参数τ变化的结论,其中αpu是泵浦光的吸收系数,v是介质的声速;In consideration of the duration of oscillation, a parameter τ = 1/(α pu × v) is used to quantify the duration of oscillation, and it is concluded that different detection wavelengths will cause the parameter τ to change, where α pu is the absorption coefficient of the pump light and v is the sound velocity of the medium;
针对振荡振幅的考量,找到振幅最大值时的探测波长;Considering the oscillation amplitude, find the detection wavelength when the amplitude is maximum;
针对振荡频率的考量,找到振荡频率随探测波长的增加而下降的规律;Considering the oscillation frequency, we found that the oscillation frequency decreases as the detection wavelength increases;
将振荡信号的振幅A0、频率f和振荡持续时间τ三者做乘积并绘图展现出来,通过图示找出优选的探测波长。The amplitude A 0 , frequency f and oscillation duration τ of the oscillation signal are multiplied and plotted, and the optimal detection wavelength is found through the plot.
优选地,对于一个新的材料系统,只需要查到其介质介电常数对光子能量的导数、介质声速、相对形变势耦合常数、泵浦光的吸收系数、介质对于泵浦光的反射率、介质的密度、介质对于探测光的复折射率等参数即可进行仿真。Preferably, for a new material system, it is only necessary to find out the derivative of the dielectric constant of the medium with respect to the photon energy, the sound velocity of the medium, the relative deformation potential coupling constant, the absorption coefficient of the pump light, the reflectivity of the medium to the pump light, the density of the medium, the complex refractive index of the medium to the detection light and other parameters to perform simulation.
在一个具体实施例中,参见图1,是给定介质为Ge大块材料,泵浦光波长和通量分别为800nm和0.3mJ/cm2的仿真结果。在入射光波长大于300nm的范围,Ge对入射光的吸收系数呈下降趋势,因此采用越大的探测光波长,所能探测的距离就越远,在时域上的表现就是振荡的持续时间越长。用一个参数τ=1/(αpu×v)来量化振荡的持续时间,如图1(c)所示,可观测到振荡信号的时间范围τ=1/(αpu×v),以及振荡信号的频率,在探测波长小于550nm时,τ的数值仅在1ps的级别,此时很难观测到完整的振荡,图1(a)也佐证了这点,图1(a)泵浦波长为800nm,泵浦光通量为0.3mJ/cm2,探测光波长分别在539nm,652nm及1127nm下的TDBS仿真结果图,表现形式为反射率变化随时间延迟的阻尼振荡,探测波长为539nm所对应的曲线,很难观测到一个完整的振荡;而随着探测波长增加到1127nm,在200ps的范围内仍能观测到振荡,因此,由于探测波长不同引起的τ的变化,是实验中重点考量的一个因素。In a specific embodiment, see Figure 1, which is a simulation result given that the medium is a bulk Ge material, the pump light wavelength and flux are 800nm and 0.3mJ/cm2 respectively. In the range of incident light wavelength greater than 300nm, the absorption coefficient of Ge to incident light shows a downward trend, so the larger the detection light wavelength, the longer the detection distance, which is manifested in the time domain as a longer duration of oscillation. A parameter τ=1/(α pu ×v) is used to quantify the duration of the oscillation. As shown in Figure 1(c), the time range of the oscillation signal τ=1/(α pu ×v) and the frequency of the oscillation signal can be observed. When the detection wavelength is less than 550nm, the value of τ is only at the level of 1ps. At this time, it is difficult to observe a complete oscillation. Figure 1(a) also confirms this point. Figure 1(a) shows the TDBS simulation results when the pump wavelength is 800nm, the pump light flux is 0.3mJ/ cm2 , and the detection light wavelengths are 539nm, 652nm and 1127nm respectively. The expression is a damped oscillation of the reflectivity change with time delay. It is difficult to observe a complete oscillation in the curve corresponding to the detection wavelength of 539nm; as the detection wavelength increases to 1127nm, oscillation can still be observed within the range of 200ps. Therefore, the change of τ caused by different detection wavelengths is a key factor to be considered in the experiment.
另一个需要重点考量的因素是振荡的振幅,因此对310-1240nm范围内的波长进行了仿真,并提取了对应的振幅,绘制出图1(b)的曲线,图1(b)提取阻尼振荡的振幅A0,其中探测光子能量在1-4eV的范围内,注意到620nm处有振幅的最大值,且随着波长增加或减小,振幅都快速地减小,这是由于620nm处对振幅的贡献最大,而随着波长增加或减少,的贡献急剧下降。Another important factor to consider is the amplitude of the oscillation. Therefore, the wavelengths in the range of 310-1240nm were simulated and the corresponding amplitudes were extracted to plot the curve in Figure 1(b). Figure 1(b) extracts the amplitude A0 of the damped oscillation, where the detection photon energy is in the range of 1-4eV. It is noted that the maximum amplitude is at 620nm, and the amplitude decreases rapidly as the wavelength increases or decreases. This is because the maximum amplitude at 620nm The contribution to the amplitude is the largest, and as the wavelength increases or decreases, The contribution has dropped sharply.
另外,振荡的频率对实验数据的提取也有一定的影响,如图1(c)所示,在310-1240nm的探测波长范围内,振荡的频率总体呈现随波长的增加而下降的趋势。在TDBS实验中,由于噪声及光电器件精度的影响,振荡频率越大意味着在同一个时间尺度下能够观测到的振荡数量越多,这对于分析信号是有利的。因此,将振荡信号的振幅A0、频率f和振荡持续时间τ三者做乘积并以图1(d)的形式展示出来,图1(d)振荡信号的振幅A0、频率f与观测范围τ三者的乘积随探测波长的变化,反映了三个重要指标对于实验结果的贡献,灰度框部分是TDBS实验建议最佳波长选取范围。从图中可以看到,虽然在1100nm后有较高的数值,但考虑振荡的振幅过小,不纳入考虑范围。而灰度区域覆盖的波长范围,无论是振荡的振幅、频率,还是持续时间,都符合实验结果的要求,因此可认为在Ge的TDBS实验中,最佳的探测波长选取范围为600-820nm。In addition, the frequency of oscillation also has a certain influence on the extraction of experimental data. As shown in Figure 1(c), within the detection wavelength range of 310-1240nm, the frequency of oscillation generally shows a downward trend with the increase of wavelength. In the TDBS experiment, due to the influence of noise and the accuracy of optoelectronic devices, the larger the oscillation frequency, the more oscillations can be observed at the same time scale, which is beneficial for analyzing the signal. Therefore, the amplitude A0, frequency f and oscillation duration τ of the oscillation signal are multiplied and displayed in the form of Figure 1(d). Figure 1(d) shows that the product of the amplitude A0, frequency f and observation range τ of the oscillation signal changes with the detection wavelength, reflecting the contribution of the three important indicators to the experimental results. The gray box part is the optimal wavelength selection range recommended for the TDBS experiment. As can be seen from the figure, although there are higher values after 1100nm, the amplitude of the oscillation is too small and is not taken into consideration. The wavelength range covered by the grayscale area, whether it is the oscillation amplitude, frequency, or duration, meets the requirements of the experimental results. Therefore, it can be considered that in the Ge TDBS experiment, the optimal detection wavelength selection range is 600-820nm.
优选地,所述方法还包括对找出的优选探测波长进行多次仿真实验进行验证。Preferably, the method further comprises performing multiple simulation experiments to verify the found preferred detection wavelength.
具体地,参见图2,图2(a)是在相同介质和泵浦条件下的TDBS实验结果图,曲线由背景热信号和阻尼振荡组成。图2(a)泵浦波长为800nm,泵浦光通量为0.3mJ/cm2,探测光子能量分别在625,730,870,920nm下的TDBS实验结果图,通过减去背景信号可以提取出阻尼振荡的形式,对实验数据处理后提取上述提到的三个参数并作乘积绘制出如图2(b)所示的曲线,图2(b)与图1(d)类似,纵轴由实验数据提取出的三个参数的乘积组成,横轴是对应的探测波长。其具有和图1(d)一样的趋势,证明了该仿真对实验有指导意义。Specifically, see Figure 2. Figure 2(a) is a TDBS experimental result diagram under the same medium and pumping conditions. The curve consists of background thermal signals and damped oscillations. Figure 2(a) is a TDBS experimental result diagram with a pump wavelength of 800nm, a pump light flux of 0.3mJ/ cm2 , and detection photon energies of 625, 730, 870, and 920nm. By subtracting the background signal, the damped oscillation form can be extracted. After processing the experimental data, the three parameters mentioned above are extracted and multiplied to draw a curve as shown in Figure 2(b). Figure 2(b) is similar to Figure 1(d). The vertical axis consists of the product of the three parameters extracted from the experimental data, and the horizontal axis is the corresponding detection wavelength. It has the same trend as Figure 1(d), proving that the simulation has guiding significance for the experiment.
基于本发明的实施例2Embodiment 2 based on the present invention
本发明实施例2所提供的一种时域布里渊散射实验中优选探测光波长的装置300可执行本发明任意实施例所提供的时域布里渊散射实验中优选探测光波长的方法,具备执行方法相应的功能模块和有益效果,该系统可以由软件和/或硬件(集成电路)的方式实现,并一般可集成于服务器或终端设备中。The device 300 for preferentially detecting the wavelength of light in a time-domain Brillouin scattering experiment provided in Example 2 of the present invention can execute the method for preferentially detecting the wavelength of light in a time-domain Brillouin scattering experiment provided in any embodiment of the present invention, and has functional modules and beneficial effects corresponding to the execution method. The system can be implemented by software and/or hardware (integrated circuit) and can generally be integrated into a server or terminal device.
图3是本发明实施例2中的一种时域布里渊散射实验中优选探测光波长的装置300的结构示意图。参照图3,本发明实施例的时域布里渊散射实验中优选探测光波长的装置300具体可以包括:FIG3 is a schematic diagram of the structure of a device 300 for preferentially detecting light wavelengths in a time-domain Brillouin scattering experiment in Embodiment 2 of the present invention. Referring to FIG3 , the device 300 for preferentially detecting light wavelengths in a time-domain Brillouin scattering experiment in the embodiment of the present invention may specifically include:
探测光波长关系确认模块310,用于基于时域布里渊散射实验中反射率的变化公式,得到当探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论;The detection light wavelength relationship confirmation module 310 is used to obtain the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the detection light wavelength changes based on the reflectivity change formula in the time-domain Brillouin scattering experiment;
仿真实验模块320,用于通过仿真实验来寻找优选的探测波长,具体包括:The simulation experiment module 320 is used to find the optimal detection wavelength through simulation experiments, specifically including:
对于一个材料系统,根据所述材料系统的参数进行仿真:For a material system, simulation is performed according to the parameters of the material system:
针对振荡持续时间的考量,利用一个参数τ=1/(αpu×v)来量化震荡的持续时间,得到探测波长不同时会引起参数τ变化的结论,其中αpu是泵浦光的吸收系数,v是介质的声速;In consideration of the duration of oscillation, a parameter τ = 1/(α pu × v) is used to quantify the duration of oscillation, and it is concluded that different detection wavelengths will cause the parameter τ to change, where α pu is the absorption coefficient of the pump light and v is the sound velocity of the medium;
针对振荡振幅的考量,找到振幅最大值时的探测波长;Considering the oscillation amplitude, find the detection wavelength when the amplitude is maximum;
针对振荡频率的考量,找到振荡频率随探测波长的增加而下降的规律;Considering the oscillation frequency, we found that the oscillation frequency decreases as the detection wavelength increases;
将振荡信号的振幅A0、频率f和振荡持续时间τ三者做乘积并绘图展现出来,通过图示找出优选的探测波长。The amplitude A 0 , frequency f and oscillation duration τ of the oscillation signal are multiplied and plotted, and the optimal detection wavelength is found through the plot.
除了上述2个模块以外,装置300还可以包括其他部件,然而,由于这些部件与本公开实施例的内容无关,因此在这里省略其图示和描述。In addition to the above two modules, the device 300 may also include other components. However, since these components are irrelevant to the content of the embodiment of the present disclosure, their illustration and description are omitted here.
一种时域布里渊散射实验中优选探测光波长的装置300的具体工作过程参照上述时域布里渊散射实验中优选探测光波长的方法实施例1的描述,不再赘述。The specific working process of the device 300 for optimally detecting the wavelength of light in a time-domain Brillouin scattering experiment is described in the above-mentioned method embodiment 1 for optimally detecting the wavelength of light in a time-domain Brillouin scattering experiment, and will not be repeated here.
综合上述各实施例提供的时域布里渊散射实验中优选探测光波长的方法及装置,利用探测光波长发生变化时,阻尼振荡的振幅、频率以及衰减系数均会被调制的结论,以此进行计算机仿真实验来优选出探测光波长。本发明优选出时域布里渊散射实验中的探测光波长,可以获得更优的利用时域布里渊散射实验进行检测分析的效果。Based on the above-mentioned methods and devices for optimizing the wavelength of the detection light in the time-domain Brillouin scattering experiment provided by the embodiments, the conclusion that the amplitude, frequency and attenuation coefficient of the damped oscillation will be modulated when the wavelength of the detection light changes is used to perform a computer simulation experiment to optimize the wavelength of the detection light. The present invention optimizes the wavelength of the detection light in the time-domain Brillouin scattering experiment, and can obtain a better effect of detection and analysis using the time-domain Brillouin scattering experiment.
注意,上述仅为本发明的较佳实施例及所运用技术原理。本领域技术人员会理解,本发明不限于这里所述的特定实施例,对本领域技术人员来说能够进行各种明显的变化、重新调整和替代而不会脱离本发明的保护范围。因此,虽然通过以上实施例对本发明进行了较为详细的说明,但是本发明不仅仅限于以上实施例,在不脱离本发明构思的情况下,还可以包括更多其他等效实施例,而本发明的范围由所附的权利要求范围决定。Note that the above are only preferred embodiments of the present invention and the technical principles used. Those skilled in the art will understand that the present invention is not limited to the specific embodiments described herein, and that various obvious changes, readjustments and substitutions can be made by those skilled in the art without departing from the scope of protection of the present invention. Therefore, although the present invention has been described in more detail through the above embodiments, the present invention is not limited to the above embodiments, and may include more other equivalent embodiments without departing from the concept of the present invention, and the scope of the present invention is determined by the scope of the appended claims.
Claims (6)
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211156899.7A CN115479921B (en) | 2022-09-22 | 2022-09-22 | Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202211156899.7A CN115479921B (en) | 2022-09-22 | 2022-09-22 | Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115479921A CN115479921A (en) | 2022-12-16 |
CN115479921B true CN115479921B (en) | 2024-09-20 |
Family
ID=84394430
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202211156899.7A Active CN115479921B (en) | 2022-09-22 | 2022-09-22 | Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115479921B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102679907A (en) * | 2012-06-01 | 2012-09-19 | 哈尔滨工业大学深圳研究生院 | High-precision differential interference measuring system and method based on LED light source |
CN110672566A (en) * | 2019-10-18 | 2020-01-10 | 南昌航空大学 | A Method for Detecting the Structure of Acoustic Photonic Crystals in Water by Stimulated Brillouin Scattering |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4100574B2 (en) * | 2005-12-01 | 2008-06-11 | 国立大学法人 東京大学 | Optical fiber characteristic measuring apparatus and optical fiber characteristic measuring method |
US20130330711A1 (en) * | 2012-06-06 | 2013-12-12 | National Taiwan University | Sensor for detection of a target of interest |
EP3140620A1 (en) * | 2014-05-05 | 2017-03-15 | Filippo Bastianini | Apparatus for interrogating distributed optical fibre sensors using a stimulated brillouin scattering optical frequency-domain interferometer |
KR101727091B1 (en) * | 2016-06-30 | 2017-05-04 | 한국광기술원 | Apparatus and method of distributed fiber sensor using optical frequency domain reflectometry based on Brillouin dynamic grating |
KR102429715B1 (en) * | 2020-11-19 | 2022-08-04 | 중앙대학교 산학협력단 | Optical frequency-domain reflectometry device and method of brillouin dynamic grating based on the modulation of a single light source |
-
2022
- 2022-09-22 CN CN202211156899.7A patent/CN115479921B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102679907A (en) * | 2012-06-01 | 2012-09-19 | 哈尔滨工业大学深圳研究生院 | High-precision differential interference measuring system and method based on LED light source |
CN110672566A (en) * | 2019-10-18 | 2020-01-10 | 南昌航空大学 | A Method for Detecting the Structure of Acoustic Photonic Crystals in Water by Stimulated Brillouin Scattering |
Also Published As
Publication number | Publication date |
---|---|
CN115479921A (en) | 2022-12-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Huang et al. | Low reflection effect by 3D printed functionally graded acoustic black holes | |
Zarei et al. | Laser ultrasonics for nondestructive testing of composite materials and structures: a review | |
Rogers et al. | Optical generation and characterization of acoustic waves in thin films: Fundamentals and applications | |
Laureti et al. | Comparative study between linear and non-linear frequency-modulated pulse-compression thermography | |
Fahim et al. | Model-based damage reconstruction in composites from ultrasound transmission | |
Ronchi et al. | Discrimination of nano-objects via cluster analysis techniques applied to time-resolved thermo-acoustic microscopy | |
Xu et al. | Surface/sub-surface crack-scattered nonlinear rayleigh waves: a full analytical solution based on elastodynamic reciprocity theorem | |
Larsen et al. | Finite-difference time-domain analysis of time-resolved terahertz spectroscopy experiments | |
Xu et al. | Research on the detection of the hole in wood based on acoustic emission frequency sweeping | |
Moura et al. | Depth evaluation of surface-breaking cracks using laser-generated transmitted Rayleigh waves | |
Ruiz-Veloz et al. | Thermally corrected solutions of the one-dimensional wave equation for the laser-induced ultrasound | |
CN115479921B (en) | Method and device for optimizing detection light wavelength in time domain Brillouin scattering experiment | |
Sun et al. | Laser ultrasonic spatially resolved acoustic spectroscopy for grain size study based on improved variational mode decomposition (IVMD) | |
Golato et al. | Multimodal sparse reconstruction in guided wave imaging of defects in plates | |
Wei et al. | A method for reconstructing two-dimensional surface and internal temperature distributions in structures by ultrasonic measurements | |
Demirkiran et al. | Analysis of microcantilevers excited by pulsed-laser-induced photoacoustic waves | |
Bouzzit et al. | Study of Rayleigh waves interaction with a spherical ball in contact with a plane surface for the development of new NDT method for ball bearings | |
Jiang et al. | Acoustic emission source localization technique based on least squares support vector machine by using FBG sensors | |
Zhao et al. | Dynamic surface acoustic response to a thermal expansion source on an anisotropic half space | |
Qi et al. | Detection of parallel double crack in pressure vessel based on optical fibre ultrasonic sensing | |
Ong et al. | Non-contact measurement techniques for structural health monitoring | |
Chen et al. | Novel optical fiber dynamic light scattering measurement system for nanometer particle size | |
Bai et al. | A numerical study of shock waves generated through laser ablation of explosives | |
Munian et al. | Ultrasonic guided waves through Coalesence of matrix cracks and delaminations in laminated composites | |
Zhang et al. | Defect modes and localisation of quasi-Lamb waves along a sidewall of corrugated aluminium plates |
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 |