CN107831094B - 基于碱金属原子弛豫率变化测量气体扩散常数的方法 - Google Patents
基于碱金属原子弛豫率变化测量气体扩散常数的方法 Download PDFInfo
- Publication number
- CN107831094B CN107831094B CN201711033721.2A CN201711033721A CN107831094B CN 107831094 B CN107831094 B CN 107831094B CN 201711033721 A CN201711033721 A CN 201711033721A CN 107831094 B CN107831094 B CN 107831094B
- Authority
- CN
- China
- Prior art keywords
- magnetic field
- alkali metal
- relaxation rate
- gradient
- gas diffusion
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
Landscapes
- Physics & 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)
- General Physics & Mathematics (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analysing Materials By Optical Means (AREA)
Abstract
本发明提供了一种基于碱金属原子弛豫率变化测量气体扩散常数的方法,包括步骤:用二次函数y=a*(x+b)2+c对测量得到的横向弛豫率随纵向磁场梯度变化规律进行拟合,得到拟合常数a的数值;通过公式计算得到气体扩散常数D。将磁场梯度引起的碱金属原子弛豫这一弊端转化为可利用的资源,充分利用了磁力仪中碱金属原子横向弛豫率受磁场梯度影响这一特点,提出了一种利用磁力仪中碱金属原子横向弛豫率随纵向磁场梯度变化情况实现气体扩散常数测量的方案,达到了“变废为宝”的效果。
Description
技术领域
本发明涉及测量技术领域,特别地,涉及一种利用磁力仪中碱金属原子横向弛豫率随磁场梯度变化关系实现气体扩散常数测量的方法。
背景技术
扩散作为一种普遍现象,遍布于化学工业、食品安全以及建筑材料等各个领域,扩散过程较为复杂。扩散系数作为扩散的重要物性参数,对于气体扩散常数的测量,有利于加深对扩散过程的了解。然而,目前还没有一种统一的方法用于气体扩散常数的测量。
对于气体扩散系数的实验测定,国内外学者给出了不同的测量方法,主要包括:激光全息干涉法([1]何茂刚,郭盈,钟秋,等.激光全息干涉法测量二元气体扩散系数[J].工程热物理学报,2010,V31(3):369-372.[2]王操.高温高压下二元气体扩散系数的激光测试系统与仿真[D].华中科技大学,2012.)、气相色谱分析法([3]Karaiskakis G,GavrilD.Determination of diffusion coefficients by gas chromatography.[J].Cheminform,2004,35(35):147.)、膜池法([4]Villet M C,Gavalas G R.Measurement ofconcentration-dependent gas diffusion coefficients in membranes from apsuedo-steady state permeation run[J].Journal of Membrane Science,2007,297(1–2):199-205.)和Stefan扩散管法([5]徐洪峰,王国香,史继诚,等.二元气体扩散系数的一种测定方法[J].大连交通大学学报,2012,33(5):90-92.)等。其中,Stefan扩散管法因其实验装置简单、操作方便、实验数据精确度高等优点,是迄今为止测量气体扩散系数最常用的方法,但也存在技术上的不足。
发明内容
本发明目的在于提供一种利用磁力仪中碱金属原子横向弛豫率随磁场梯度变化关系实现气体扩散常数测量的全新的方法,以解决背景技术中的技术问题。
为实现上述目的,本发明提供了基于碱金属原子弛豫率变化测量气体扩散常数的方法,包括步骤:
A、泵浦光沿z轴方向传输,探测光沿x轴方向传输,所述泵浦光和探测光均经过加热至至工作温度的铯原子气室;
B、在z轴方向上对铯原子气室施加静磁场,在x轴方向上对铯原子气室施加交变磁场;
C、改变梯度线圈中电流大小和方向,使得纵向磁场梯度在一定范围内变化,利用自由感应衰减法测量不同磁场梯度下铯原子的横向弛豫率;
D、用二次函数y=a*(x+b)2+c对测量得到的横向弛豫率随纵向磁场梯度变化规律进行拟合,得到拟合常数a的数值;
E、通过公式计算得到气体扩散常数D,
其中,R为原子气室半径,γ为碱金属原子的旋磁比。
优选的,步骤A中,所述铯原子气室的温度为60℃。
优选的,步骤A中,所述铯原子气室内充有He作为缓冲气体,N2作为淬火气体。
优选的,所述静磁场的强度为10μT。
优选的,所述交变磁场的强度为1μT。
优选的,所述纵向磁场梯度在-20nT/mm~20nT/mm范围内变化。
本发明具有以下有益效果:
本发明将磁场梯度引起的碱金属原子弛豫这一弊端转化为可利用的资源,提出了一种利用磁力仪中碱金属原子横向弛豫率随纵向磁场梯度变化情况实现气体扩散常数测量的方案。该方案充分利用了磁力仪中碱金属原子横向弛豫率受磁场梯度影响这一特点,达到了“变废为宝”的效果。
为清楚阐述本发明所述方法的作用原理,下面对碱金属原子横向弛豫率(1/T2)随磁场梯度变化作简要介绍:
对于充有缓冲气体的原子气室,碱金属原子密度矩阵演化方程为:
(1)式中,第一项表示自由原子哈密顿量的演化,还包括碱金属原子与外磁场的相互作用;最后一项表示碱金属原子的空间扩散作用;其余项表示碱金属原子间自旋交换碰撞、碱金属原子与缓冲气体间自旋破坏碰撞,以及光泵浦作用等。为重点研究磁场梯度对碱金属原子横向弛豫率的影响,可将密度矩阵演化方程简化为:
利用微扰法求解得到碱金属原子横向弛豫率:
当原子气室内缓冲气体压强较大时,碱金属原子横向弛豫率可近似为:
当气室内缓冲气体压强较低时,碱金属原子横向弛豫率为:
对于实际原子气室,碱金属原子的弛豫来源还包括自旋交换弛豫、自旋破坏弛豫,以及光泵浦作用引起的弛豫等。由于上述弛豫作用与磁场梯度无关,在特定的实验环境下,对于特定的原子气室,可将上述几项弛豫作用视为常数。因此,碱金属原子横向弛豫率随气室内磁场梯度变化关系满足公式(6):
根据碱金属原子横向弛豫率表达式可知,碱金属原子横向弛豫率随纵向磁场梯度的变化满足二次关系式。
本申请实验中,沿z轴方向给碱金属原子气室主动施加一个磁场梯度,通过改变梯度线圈中电流值和通电方向实现磁场梯度大小和方向的调节。利用自由感应衰减法测量得到不同纵向磁场梯度下铯原子的系列横向弛豫率。用二次函数y=a*(x+b)2+c对测量得到的横向弛豫率随纵向磁场梯度变化规律进行拟合,得到拟合常数a、b、c。由于常数a依赖于原子气室半径R、气体扩散常数D以及碱金属原子的旋磁比γ,根据下列公式(7)进而得到气体扩散常数D,
利用本方案实现气体扩散常数的测量,对于加深对传质过程的了解具有重要的理论和工程实践意义。
除了上面所描述的目的、特征和优点之外,本发明还有其它的目的、特征和优点。下面将参照图,对本发明作进一步详细的说明。
附图说明
构成本申请的一部分的附图用来提供对本发明的进一步理解,本发明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中:
图1是本发明优选实施例的磁场梯度线圈结构示意图;
图2是本发明优选实施例的气体扩散常数测量实验装置图;
其中,1.泵浦光激光器,2.探测光激光器,3.偏振片,4.衰减片,5.扩束准直系统,6.1/4波片,7.五层磁屏蔽系统,8.加热装置,9.铯原子气室,10.主线圈,11.副线圈,12.磁场梯度线圈,13.偏振分光片PBS,14.平衡探测器,15.信号处理系统;
图3是本发明优选实施例的铯原子横向弛豫率随纵向磁场梯度变化曲线及二次拟合结果(300Torr He,50Torr N2);
图4是本发明优选实施例的铯原子横向弛豫率随纵向磁场梯度变化曲线及二次拟合结果(50Torr N2)。
具体实施方式
以下结合附图对本发明的实施例进行详细说明,但是本发明可以根据权利要求限定和覆盖的多种不同方式实施。
参见图1、图2,本发明规定沿z轴方向的磁场为纵向磁场,发明中利用一对绕向相反的同轴线圈(后文中称之为磁场梯度线圈)产生纵向磁场梯度,通过改变梯度线圈中电流大小和方向实现磁场梯度大小和方向的调节。磁场梯度线圈如图1中所示,线圈匝数为n,电流值为I,线圈半径为R,间距为2d,产生纵向磁场梯度为其中,κ为线圈结构常数,满足:
将铯原子气室9置于五层磁屏蔽系统7中,加热装置8对铯原子气室9进行加热。主线圈10产生纵向稳态磁场B0,副线圈11产生横向交变磁场B1cos(ωt),梯度线圈12产生纵向磁场梯度。泵浦光激光器1产生D1线线偏振光,经偏振片3和1/4波片6后转化为圆偏振光沿z轴传输进入到气室9对铯原子进行泵浦,衰减片4实现泵浦光光强的调节,扩束准直系统5改善光束质量。探测光激光器2产生D1线光经偏振片3后变为线偏光,经准直系统5和衰减片4后沿x方向传输,通过原子气室9后探测光偏振面发生一定角度的偏转,透射光经由偏振分光片13分解为两束正交的线偏振光,平衡探测器14探测两束光光强,最后将信号送入信号处理系统15进行处理。
本发明一种基于磁力仪中碱金属原子横向弛豫率随纵向磁场梯度变化关系实现气体扩散常数测量的方法,主要包括以下步骤:
1、打开泵浦光和探测光激光器,使其分别稳定在相应的频率上
打开泵浦光激光器1使其波长稳定在894nm(铯原子D1线对应波长),泵浦光经偏振片3和1/4波片6后变为圆偏振光,泵浦光经利用衰减片4将泵浦光功率调节至1mW,经准直扩束5后沿z轴方向传输对铯原子进行泵浦。
打开探测光激光器2使其波长稳定在894nm,利用衰减片4将探测光功率调节至30μW,探测光经偏振片3经准直扩束5后沿x轴方向传输经铯原子气室后透射。
2、利用加热装置将铯原子气室加热至其工作温度。实验中所加梯度场的调节范围是根据具体的实验环境确定的,例如磁屏蔽系统内剩余梯度场的大小、梯度线圈匝数、大小以及通电电流值等;实验中只要保证外加梯度场在一定范围内变化即可观测到相应的碱金属原子横向弛豫率随磁场梯度的二次变化关系。
球形铯原子气室9采用Pyrex玻璃吹制而成,半径R为10mm,气室内充有He作为缓冲气体,N2作为淬火气体。利用加热装置8将铯原子气室9加热至60℃。
3、对主线圈通稳恒电流,使其产生实验所需静磁场B0
对主线圈10通稳恒电流使其产生沿z轴方向10μT的静磁场。
4、对副线圈施加交变电流,使其产生横向交变磁场B1cos(ωt)
对副线圈11通交变电流使其产生大小约为1μT的横向交变磁场。
5、改变梯度线圈中通电电流值和方向实现纵向磁场梯度大小和方向的调节;
改变梯度线圈12中电流大小和方向使得纵向磁场梯度在-20nT/mm~20nT/mm范围内变化,利用自由感应衰减法测量不同磁场梯度下铯原子的横向弛豫率。磁场变化范围是根据实验中具体的梯度线圈尺寸及通电电流大小确定的,并非一定要在这个范围内变化。
6、测量不同纵向磁场梯度下铯原子对应的横向弛豫率,得到铯原子横向弛豫率随纵向磁场梯度变化规律。
例如,利用二次函数y=a*(x+b)2+c对6中测量结果进行拟合,得到图3、图4中a值分别为4.5969*105、7.9418*105。根据已知的铯原子半径R=10mm,铯原子旋磁比γ=3.5Hz/nT,根据公式得到图3、4中铯原子气室中气体扩散常数分别为0.1218cm2/s、0.7050cm2/s。
实验证实,本方案可行性较高,且实验结果与理论吻合。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (6)
2.根据权利要求1所述的方法,其特征在于,步骤A中,所述铯原子气室的温度为60℃。
3.根据权利要求1所述的方法,其特征在于,步骤A中,所述铯原子气室内充有He作为缓冲气体,N2作为淬火气体。
4.根据权利要求1所述的方法,其特征在于,所述静磁场的强度为10μT。
5.根据权利要求1所述的方法,其特征在于,所述交变磁场的强度为1μT。
6.根据权利要求1所述的方法,其特征在于,所述纵向磁场梯度在-20nT/mm~20nT/mm范围内变化。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711033721.2A CN107831094B (zh) | 2017-10-30 | 2017-10-30 | 基于碱金属原子弛豫率变化测量气体扩散常数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711033721.2A CN107831094B (zh) | 2017-10-30 | 2017-10-30 | 基于碱金属原子弛豫率变化测量气体扩散常数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107831094A CN107831094A (zh) | 2018-03-23 |
CN107831094B true CN107831094B (zh) | 2020-04-07 |
Family
ID=61650801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711033721.2A Active CN107831094B (zh) | 2017-10-30 | 2017-10-30 | 基于碱金属原子弛豫率变化测量气体扩散常数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107831094B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111220934A (zh) * | 2020-01-14 | 2020-06-02 | 杭州昕磁科技有限公司 | 基于脉冲泵浦磁力仪的梯度检测系统 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4414509A (en) * | 1980-11-26 | 1983-11-08 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Low energy electron magnetometer using a monoenergetic electron beam |
US8427146B2 (en) * | 2009-06-26 | 2013-04-23 | Seiko Epson Corporation | Magnetic sensor for measuring a magnetic field using optical pumping method |
CN105651649A (zh) * | 2016-01-27 | 2016-06-08 | 东南大学 | 一种适用于原子磁强计的原子密度实时在线测量方法 |
Family Cites Families (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
HUP0103605A3 (en) * | 1998-06-17 | 2002-11-28 | Medi Physics Inc Princeton | Resilient containers for hyperpolarized gases |
US6653832B2 (en) * | 2001-03-09 | 2003-11-25 | Battelle Memorial Institute | Method for high resolution magnetic resonance analysis using magic angle technique |
KR20050008788A (ko) * | 2002-06-05 | 2005-01-21 | 파마시아 앤드 업존 캄파니 엘엘씨 | 고처리량 스크리닝을 위한 불소 nmr의 용도 |
CN1279370C (zh) * | 2003-07-03 | 2006-10-11 | 石油大学(北京) | 测量井下套管内流体特性的探测器及探测系统 |
US7521928B2 (en) * | 2006-11-07 | 2009-04-21 | Trustees Of Princeton University | Subfemtotesla radio-frequency atomic magnetometer for nuclear quadrupole resonance detection |
CN201181274Y (zh) * | 2008-03-19 | 2009-01-14 | 苏州特尔纳米技术有限公司 | 电子显微镜下微纳样本的无线操作机构 |
WO2014008226A1 (en) * | 2012-07-02 | 2014-01-09 | Millikelvin Technologies Llc | Techniques, systems and machine readable programs for magnetic resonance |
CN204154603U (zh) * | 2014-10-28 | 2015-02-11 | 延边大学 | 透皮扩散池电磁搅拌器 |
US10139347B2 (en) * | 2015-09-23 | 2018-11-27 | Halliburton Energy Services, Inc. | Measurement of noble gas adsorption via laser-induced breakdown spectroscopy for wettability determination |
CN106597338B (zh) * | 2016-12-28 | 2019-03-29 | 北京航空航天大学 | 一种基于电子共振相频分析测量原子横向弛豫时间的方法 |
CN107192633A (zh) * | 2017-07-10 | 2017-09-22 | 北京航空航天大学 | 一种serf态下在线测量原子磁强计气室内碱金属密度的方法 |
-
2017
- 2017-10-30 CN CN201711033721.2A patent/CN107831094B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4414509A (en) * | 1980-11-26 | 1983-11-08 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Low energy electron magnetometer using a monoenergetic electron beam |
US8427146B2 (en) * | 2009-06-26 | 2013-04-23 | Seiko Epson Corporation | Magnetic sensor for measuring a magnetic field using optical pumping method |
CN105651649A (zh) * | 2016-01-27 | 2016-06-08 | 东南大学 | 一种适用于原子磁强计的原子密度实时在线测量方法 |
Non-Patent Citations (5)
Title |
---|
Determination of Diffusion Coefficients by Gas Chromatography;Karaiskakis G 等;《ChemInform》;20040806;第35卷(第35期);摘要 * |
Measurement of concentration-dependent gas diffusion coefficients in membranes from a psuedo-steady state permeation run;Villet M C 等;《Journal of Membrane Science》;20070705;第297卷(第1-2期);第199-205页 * |
NEW METHOD FOR PREDICTION OF BINARY GAS-PHASE DIFFUSION COEFFICIENTS;Fuller E N 等;《Industrial & Engineering Chemistry》;19660531;第58卷(第5期);第18-27页 * |
Relaxation of excited states of an emitter near a metal nanoparticle: An analysis based on superradiance theory;Protsenko I E 等;《Journal of Experimental and Theoretical Physics》;20140228;第119卷(第2期);第227-241页 * |
Statistical analysis of formation and relaxation of atom clusters based on data of molecular dynamic modeling of gas-phase nucleation of metallic nanoparticles;Korenchenko A E 等;《High Temperature》;20160630;第54卷(第6期);第243–248页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107831094A (zh) | 2018-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Martynov et al. | Exploring the sensitivity of gravitational wave detectors to neutron star physics | |
CN108693488B (zh) | 一种基于双抽运光束的无自旋交换弛豫原子自旋磁场测量装置 | |
CN106886000B (zh) | 一种利用核磁共振实现磁场幅度稳定的装置及方法 | |
CN108717168B (zh) | 一种基于光场幅度调制的标量磁场梯度测量装置及方法 | |
CN105301541B (zh) | 原子磁强计的磁线圈x、y轴非正交角的测量装置与方法 | |
Santos-da-Costa et al. | Two new tests to the distance duality relation with galaxy clusters | |
CN113777106B (zh) | 一种原子磁强计的碱金属原子数密度空间分布均匀性的测试系统及方法 | |
JP2018004462A (ja) | 磁場計測装置、磁場計測装置の調整方法、および磁場計測装置の製造方法 | |
CN110568381B (zh) | 基于双光束三轴矢量原子磁强计光磁非正交角原位测量方法 | |
CN105929458A (zh) | 航空磁场矢量检测装置及监测方法 | |
CN111025201A (zh) | 一种原子磁力仪探头光路结构 | |
Wei et al. | In-situ measurement of the density ratio of K-Rb hybrid vapor cell using spin-exchange collision mixing of the K and Rb light shifts | |
CN107831094B (zh) | 基于碱金属原子弛豫率变化测量气体扩散常数的方法 | |
Audet et al. | Electron injector for compact staged high energy accelerator | |
Vlasova et al. | High-sensitive absorption measurement in ultrapure quartz glasses and crystals using time-resolved photothermal common-pass interferometry and its possible prospects | |
Parks et al. | A simple pendulum laser interferometer for determining the gravitational constant | |
Fu et al. | Analysis on the effect of electron spin polarization on a hybrid optically pumped k-rb-21 ne co-magnetometer | |
Nasyrov et al. | Simple method for characterization of anti-relaxation coating of optical cells | |
Zavattini et al. | Optical polarimetry for fundamental physics | |
Yang et al. | Optimizing noble gas pressure for enhanced self-compensation in spin-exchange relaxation-free comagnetometers | |
Zhu et al. | Magnetic effect in the test of the weak equivalence principle using a rotating torsion pendulum | |
Qiang et al. | Investigation of spatial spin polarization of 87Rb and 129Xe atomic ensemble pumped by Gaussian beam in a miniature nuclear magnetic resonance gyroscope | |
Shi et al. | Optimized design and calibration of the triaxis induction magnetometer with crosstalk and nonorthogonality compensation | |
Ruan et al. | Effects of pump laser intensity on the cell temperature working point in a K-Rb-21 Ne spin-exchange relaxation-free co-magnetometer | |
Motuz et al. | Theoretical assessment of the OTDR detector noise on plasma current measurement in tokamaks |
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 |