Probing the rotation curve of the outer accretion disk in FU Orionis objects with long-wave

合集下载

托福阅读curve

托福阅读curve

托福阅读curveTOEFL Reading CurveThe TOEFL Reading section is known for its challenging texts and critical thinking questions, making it a crucial component of the overall TOEFL test. One aspect that test-takers often ponder is the "curve" or scaling system used to calculate their final scores. However, it is important to note that the TOEFL Reading section does not have a predetermined curve like some other standardized tests; instead, it follows a predetermined scoring rubric.The TOEFL Reading section is designed to assess a test-taker's ability to comprehend and analyze academic texts commonly encountered in universities and colleges. This section consists of three to four passages, each followed by a set of 10 to 14 questions. The questions can be in various formats, such as multiple-choice, drag-and-drop, or matching.While the TOEFL Reading section does not have a set curve, the scoring rubric takes into account the difficulty level of the passages and the performance of the test-takers. Each question in the section is worth one point, resulting in a total score of 30 or 40 points, depending on the number of passages. The raw scores are then converted to a scaled score ranging from 0 to 30, which is the score ultimately reported to test-takers.The scaled score is based on the statistical performance of all test-takers who have taken the same version of the TOEFL Reading section. The performance of the test-takers helps determine the difficulty level of the questions, and subsequently, the scaled score. Therefore, the curve is not predetermined and can vary from test to test.To perform well on the TOEFL Reading section, it is essential to develop solid reading comprehension skills. Some effective strategies include active reading, identifying the main ideas and supporting details, understanding the organization and structure of the passage, and practicing with sample questions.In conclusion, the TOEFL Reading section does not have a specific curve, but instead relies on a predetermined scoring rubric that considers the difficulty level of the questions and the test-takers' performance. Developing strong reading skills and practicing with sample questions are key to performing well on this section of the TOEFL test.。

推理正确率分布曲线 英文

推理正确率分布曲线 英文

推理正确率分布曲线英文English:The distribution curve of inference accuracy represents the spread of correct inference rates across a population or dataset. This curve typically takes the shape of a bell curve, with the majority of individuals or data points clustering around the mean accuracy rate, and fewer individuals or data points at the extremes of high and low accuracy. The central tendency of this curve is indicative of the average accuracy level, while the spread or standard deviation reflects the variability or consistency of inference accuracy within the population or dataset. Factors such as cognitive ability, prior knowledge, experience, and task complexity can influence the shape and parameters of the accuracy distribution curve. In educational settings, understanding the distribution of inference accuracy can inform instructional strategies, curriculum design, and assessment practices, allowing educators to tailor interventions and support to meet the diverse needs of learners. Additionally, in fields like psychology and neuroscience, analyzing individual differences in inference accuracy distribution can shed light on cognitive processes,learning mechanisms, and neural substrates underlying reasoning and decision-making.Translated content:推理正确率分布曲线代表了在人群或数据集中正确推理率的分布情况。

石油英语词汇(E)

石油英语词汇(E)

石油英语词汇(E1)--------------------------------------------------------------------------------E & ST 适用性试验E by N 东偏北E by S 东偏南E long 东经e 弹性e 当量E 地;接地;地线e 电子E 东E 恩氏度E 工程E 空E 能e 偏心度E 气体膨胀系数e 误差e 效率E 早第三纪;下第三系e 自然对数的底E&P 勘探及生产E-field ratio telluric method 大地电流场比法E-H tee 超高频T形接头E-layer E电离层E-log 电测井E-mail 电子信箱E-phase method 电场相位法E-Z tree 简易树e.a. 弹性轴E.B.B. 最高质量e.c. 漆包线e.c.d.r. 外部临界阻尼电阻e.r 在途中e.r. 蒸发率e.s. 接地开关E.S. 静电的e.s. 静电学e.s.c.g.s 厘米-克-秒静电制单位E.S.D 超硬铝E.S.H.P 总有效马力E.U. 能量单位E.U. 熵单位Ea 地;接地;地线ea 每EAC 接触放电加工EAC 循环进位EAD 放电涂覆处理EAE 扩展的运算单元EAEG 欧洲地球物理勘探工作者协会EAG 电磁波传播衰减增益eagle 鹰;鹰金元eaglestone 鹰石EAM 电子记帐机EAPG 欧洲石油地学家和工程师协会EAPG 欧洲石油地质学家协会ear muffs 护耳套;耳机上防噪橡皮护圈ear plugs 耳塞ear protectors 护耳器ear sensitivity 听觉灵敏度ear 耳;耳柄;耳状物;吊环earing 出耳子;形成花边Earlachitina 厄尔几丁虫属Earlougher type curve 埃洛弗尔解释图版early breakthrough 早期突破early detection 早期检测early development 早期开发early diagenesis 早期成岩作用early event time 事件最早时间early failure 过早损坏;早期失效early field life 油田开发初期early finish date 最早结束时刻early pressure information 测压早期资料early pressure maintenance scheme 早期保持压力方案early production period 开采初期early production system 早期开采系统early sand control 早期防砂early setting cement 快凝水泥early simulation 早期模拟early stage 早期early start 最早开始early strength concrete 早强混凝土early strength 初期强度early time data 早期数据early time treatment 早期处理early transient regime 早期不稳定阶段early warning rada 预警雷达early warning 早期报警early water breakthrough 早期见水early 早先的early-stage waterflooding 早期注水early-time portion 早期段earmark 记号;打上记号earn 挣得earnest money 定金earnest 热心;真实;认真的;真挚的earning block 获利区块earning per share 每股收益earning performance 盈利状况earning power 收益能力;赚钱能力earning rate 收益率earning well 收益井earning 赚;所得earnings statement 收益表EAROM 电改写只读存储器earphone unit 头带受话机earphone 耳机earpiece 耳机earshot 听力范围earth acceleration 重力加速度earth albedo 地球反照率earth alkali metal 碱土金属earth anchor 地锚earth anchorage 地锚固定earth attraction 地球引力earth auger 土螺旋钻earth axis 地轴earth bar 接地棒earth bearing strength 地耐力earth bus 接地母线earth cable 接地电缆earth capacitance 对地电容earth circuit 接地电路earth clamp 接地夹earth coal 土褐煤earth color 矿物颜料earth conductivity 大地电导率earth connection 接地earth contraction 地球收缩earth coordinates 大地坐标系earth core 地核earth coupling 大地耦合earth crust 地壳earth current 大地电流earth curvature 地球曲率earth curve 地壳弯曲earth ditch 土沟earth dynamics 地球动力学earth electric field 地电场earth electricity 地电earth electrode 接地电极earth ellipsoid 地球椭圆体earth embankment lay pipe 土堤敷设管线earth embankment 土堤earth evolution model 地球演化模式earth excavation 挖土方earth exploration satellite 地球勘测卫星earth fall 土塌earth fault 接地故障earth fill 填土方earth filling 填土处理earth flow 泥流earth grabbing bucket 抓斗earth gravitational model 地球重力模型earth gravity field 地球重力场earth history 地史学earth holography 大地全息摄影术earth inductor 地磁感应器earth interior 地球内部earth lead 接地导线earth leakage 接地漏电earth magnetic effect 地磁场影响earth magnetic field 地磁场earth magnetism 地磁earth mantle 地幔earth medium 大地介质earth movement 地壳运动earth nucleus 地核Earth observation satellite 地球观测卫星earth orbital imagery 地球轨道成象earth orbital operation 绕地球轨道运行earth pad 土基earth pillar 土柱earth pitch 软沥青earth plate 接地板earth potential 地电位earth pressure 土压earth pulsation 地壳脉动earth radiation 大地辐射earth rammer 夯土机earth resistance 接地电阻;土抗力earth resources information 地球资源信息earth resources observation system 地球资源观测系统earth resources technology satellite 地球资源技术卫星earth resources terrestrial satellite 地球资源卫星earth resources test satellite 地球资源试验卫星earth resources 地球资源earth respone function 地响应函数earth return circuit 接地回线earth return 接地回线earth revolution 地球公转earth rotation 地球自转earth satellite vehicle 地球卫星运载火箭earth satellite 地球卫星earth science 地学earth screw 取土样的麻花钻earth shaking 极其重大的earth shell 地壳earth shield 接地屏蔽earth shock 地震earth silicon 硅石earth slide 滑坡earth station 地面站earth storage 土油池earth strain 地应变earth stress 地应力earth surface 地面earth switch 接地开关earth temperature 地温earth terminal 接地端子earth terrain camera 地表摄象机earth test 接地试验earth tester 接地电阻测试器earth thermometer 地温计earth tide correction 固体潮校正earth tide effect 固体潮效应earth tide gravimeter 固体潮重力仪earth tide perturbation 固体潮摄动earth tide table 固体潮值表earth tide 固体潮earth transmission characteristics 地层透射特性earth tremor 地颤;地震预兆earth tripolite 硅藻土earth wave 地波earth wax 地蜡earth zone 地球带earth 地earth's compliance factor 地球顺从系数earth's magnetic dip angle 地磁倾角earth's magnetic field 地磁场earth's magnetism 地磁earth's shadow 地球阴影earth's sphere 地球earth's spheroid 地球椭球体earth's spin vector 地球自转矢量earth-air current 大地-大气电流earth-attenuation 地层衰减earth-baseplate resonance 地面-底板共振earth-creep 土滑earth-filled bag 土袋earth-filtering effect 地层滤波效应earth-forming element 构成地球的元素earth-ionosphere waveguide 地-电离层波导earth-leakage protection 接地漏电保护earth-loss 地层损耗earth-magnetic navigation 地磁导航earth-orbital photography 地球轨道摄影术earth-orbiting satellite 地球轨道卫星earth-probing radar 地质探测雷达earth-resistivity 地电阻率earth-scraper 铲土机earth-shift 地层移动earthenware pipe 瓦管earthenware 陶器earthern pit 土质泥浆池earthing cable 接地电缆earthing device 接地装置earthing wire 地线earthing 接地earthly heat 地热earthly 地球的earthmover 土方机械;大型挖土机earthmoving equipment 运土机earthmoving 运土earthometer 兆欧计earthquake acceleration 地震加速度earthquake centre 震中earthquake counter measure 防震措施earthquake country 震区earthquake depth 震源深度earthquake epicenter 震中earthquake focus 震源earthquake intensity 地震烈度earthquake origin 震源earthquake period 地震周期earthquake precursors 地震前兆earthquake prediction 地震预报earthquake proof construction 防震建筑earthquake record 地震记录earthquake region 震区earthquake scale 地震震级earthquake seismology 天然地震学earthquake strength 地震强度earthquake wave 地震波earthquake 天然地震earthquake-resistant 抗震的earthrise 地出earthwork 土方工程;土工;土方earthy bauxite 土状铝土矿earthy element 土族元素earthy material 泥质材料earthy spring 泥泉earthy water 硬水easamatic 简易自动式的EASC 期望年度脱销费用ease of ignition 易点燃性ease of money 放松银根ease off 修正;放松ease os extinguishability 易熄灭性ease 安逸;容易;减轻EASE 电子分析和模拟设备easel 框;绘图桌easement 缓和easer 辅助炮眼easily hydrated clay 易水化粘土easily-worn parts 易损零件easing the bit in 钻头轻压慢转钻进地层East African Graben 东非地堑East African Rift Valley 东非裂谷east by north 东偏北east by south 东偏南east digging 易挖的east longitude 东经East Pacific plate 东太平洋板块east 东;东方的easterly trade wind 东贸易风Eastern hemisphere 东半球Eastern oil 美国东部石油公司easting 东横坐标;东向分量;朝东方向;东航Eastman survey instrument 伊斯门测斜仪Eastman whipstock turbine 伊斯门造斜涡轮Easto circulating sub 伊士托循环接头eastonite 富美黑云母easy bend 慢弯管easy cleavage 明显解理;明显劈理easy curve 平缓曲线easy drilling 在易破碎岩层中钻进;轻快钻进easy grade 平缓坡度easy gradient 平缓梯度easy instruction automatic computer 教学用自动计算机easy push fit 轻推配合easy running fit 轻转配合easy running 平滑运转easy slide fit 轻滑配合easy 容易的;平缓的easy-to-drill formation 易钻地层easy-to-read 易读的easy-to-use 易于使用的eat away 侵蚀eat up 消耗eater 食者;腐蚀物eating 侵蚀作用;食物;食用的EATT 电磁波传播衰减eaves 屋檐EB 电子束ebb 退落ebb-current flow structure 落潮流构造ebb-reflux 退潮ebb-tidal delta 退潮三角洲ebbing 沉陷ebbtide 落潮EBCDIC 扩充二-十进制交换码EBIT 扩展基本信息带EBM 电子束加工ebonite 硬橡皮;胶木ebony 乌木;黑檀的eboulement 崩塌ebp 终沸点EBR 电子束记录ebsemble correlation function 总体相关函数ebullated bed reactor 沸腾床反应器ebullated bed 沸腾床ebullated dryer 沸腾干燥器ebullator 循环泵ebullience 沸腾;起泡ebulliency =ebullienceebulliometer 沸点测定计ebullioscope 沸点计ebullioscopy 沸点测定法ebullition 起泡EBW 电子束焊EC CAP 电解质电容器EC test 电解腐蚀试验ec 大地电流EC 弹性系数EC 地壳EC 电解的;电解质的EC 电子俘获EC 电子计算机EC 电子控制的ec 例如EC 涡流;涡电流;杂散电流EC 误差控制EC 误差校正EC 有效浓度ec- 出;出自EC-GC 电子俘获气相色谱EC-GLC 电子俘获气液色谱ECB 环境协调委员会ecboline 麦角碱ECC 错误检查和校正eccentered gun 偏心射孔器eccentered screw pump 偏心螺杆泵eccentering arm 偏心臂eccentralizer 偏心器eccentric adjuster 偏心调节装置eccentric annulus 偏心环空eccentric anomaly 偏心异常eccentric axis 偏心轴eccentric bit 偏心钻头eccentric cam 偏心凸轮eccentric clip 偏心夹环eccentric disc 偏心圆盘eccentric distance 偏心距离eccentric drive 偏心传动eccentric gas lift valve 偏心气举阀eccentric gear 偏心齿轮eccentric gearing 偏心传动eccentric injection mandrel 偏心配水器eccentric load 偏心负荷eccentric mandrel 偏心杆eccentric orifice with flange taps 法兰取压偏心孔板eccentric orifice 偏心孔板eccentric pattern 螺旋状排列eccentric pin 偏心销eccentric production mandrel 偏心配产器eccentric pump 偏心泵eccentric rebel tool 偏心变向器eccentric reducer 偏心大小头eccentric rod 偏心杆eccentric rotating 偏心旋转eccentric rotor sliding vane compressor 偏心转子滑叶压缩机eccentric shaft 偏心轴eccentric sliding vane pump 偏心滑叶泵eccentric stabilizer 偏心稳定器eccentric tongs 偏心钳eccentric underreaming bit 管下扩眼偏心钻头eccentric water distributor 偏心配水器eccentric wear 偏磨eccentric weight 偏心锤eccentric yokes 偏心接箍eccentric 偏心轮;偏心装置;偏心的;反常的eccentrically loaded 偏心载荷的eccentricity curve 偏心率曲线eccentricity tester 径向跳动检查仪eccentricity 偏心;偏心距ECCM 电子反干扰设备Ecculiomphalus 松旋螺属ECD 当量循环密度ECD 电导检测器ECD 电子俘获探测器ECDM 电解放电加工ecf 高程校正系数ECG 电解磨削echelette 光栅echelon faults 雁行断层echelon folding 雁行褶皱作用echelon folds 雁行褶皱echelon fracture 雁行式裂缝echelon pattern 雁行构造型式echelon structure 雁行构造echelon 梯队;阶梯光栅;排成梯队echelon-sections 剖面阶梯形排列Echinatisporis 棘刺孢属echinenone 海胆酮echini echinus的复数Echinochara 刺轮藻属echinochrome 海胆色素Echinocypris 棘星介属echinoderm 棘皮动物Echinodermata 棘皮动物门echinoid 海胆类Echinoidea 海胆类Echinosphaerites 刺海林檎属Echinosporites 刺纹单缝孢属echinozoan 海胆的echinus 海胆Echiuroidae NFDAB纲echo amplifier 回声放大器echo amplitude 反射波振幅echo arrival 回波初至echo check 回送检验echo deep sounding 回声测深echo depth sounder 回声测深仪echo depth sounding sonar 回声测深仪echo depth sounding 回声测深echo depth-recorder 回声深度记录器echo distortion 回波失真echo effect 回波效应echo impulse 回波脉冲echo liquid level instrument 液面回声探测仪echo locator 回声勘定器echo ranger 回声测距仪echo ranging sonar 回声测距声呐echo ranging 回声测距echo repeater unit 回波中继续置echo sounder 回声测深仪echo sounding machine 回声测深仪echo suppressor 回波抑制器echo time 回声时间echo trouble 反射故障echo wave 回波echo 回声echo-fathom 回声测深echo-image 回波图象echo-pulse 回波脉冲echo-ranging system 回声测距系统echo-resonator 回波谐振器echo-sounding device 回声侧深仪echo-sounding instrument 回声测深仪echo-sounding receiver 回声测深接收机echo-sounding 回声探测echo-strenghth indicator 回波强度指示器echoed signal 回波信号echogram 音响测深图echograph 回声深度记录器echoing characteristics 回波特性echoing 回波现象echolocation 回声勘定echometer 回声测距仪echometry 测回声术ecidioclimate 微生态气候ECL 发射极耦合逻辑eclipse factor 阴影率eclipse 〔日ecliptic coordinate 黄道坐标ecliptic 黄道eclogite 榴辉岩ECM 电解加工ECM 电色谱ECM 电子干扰措施ECM 欧洲共同市场ECMA 欧洲计算机厂家协会ECO 电子耦合振荡器ECO 设计变动命令书eco-activist 生态活动家eco-atmosphere 生态大气eco-catastrophe 生态灾难ecochronology 生态年代学ecocide 生态灭绝ecoclimate 生态气候ecocline 生态差型ecocrisis 生态危机ecocycle rule 生态循环规律ecofactor 生态因素ecogenesis 生态发生;生态种发生ecogeography 生态地理学ecography 描述生态学ecogroup 生态群ecologic adaptation 生态适应ecologic age 生态年龄ecologic amplitude 生态幅度ecologic assemblage 生态组合ecologic balance 生态平衡ecologic botany 生态植物学ecologic change 生态变化ecologic character 生态性状ecologic community 生态群落ecologic complex 生态复合体ecologic distribution 生态分布ecologic disturbance 生态失调ecologic effect 生态效应ecologic environment 生态环境ecologic equilibrium 生态平衡ecologic facies 生态相ecologic food chain 生态食物链ecologic niche 生态境ecologic potential 生态本能ecologic reef 生态礁ecologic regime 生态状况ecologic replacement 生态更替ecologic setting 生态背景ecologic speciation 生态性物种形成ecologic succession 生态序列ecologic threshold 生态临界ecologic tolerance 生态耐性ecologic 生态的ecological charges 生态费用ecological fabric 生态织物ecological =ecologicecology of reservoir 油藏生态学ecology 生态学economatrix 数理经济学econometric analysis model 计量经济分析模型econometric model 经济计量模型econometrics 计量经济学economic algae 经济藻类economic analysis 经济分析economic and socal development strategy 经济与社会发展战略economic and technical development zone 技术经济开发区economic and technological development zone 经济开发区Economic and Trade Arbitration Commission of CCPIT 中国国际贸易促进会经贸仲裁委员会economic assessment 经济评价economic balance 经济平衡economic barometer 经济观测指标economic behaviour 经济行为economic coefficient 经济系数economic conditions 经济情况Economic Contract Law of the People's Republic of China 中华人民共和国经济合同法economic contraction 经济萎缩economic cooperation 经济合作Economic Council of Arab League 阿拉伯联盟经济理事会economic cycle 经济周期economic data 经济资料economic depletion 经济耗竭economic depth 经济深度economic dispatch 经济分配economic entity 经济实体economic evaluation model 经济评价模型economic evaluation 经济评价economic factors 经济因素economic feasibility 经济可行性economic flow rate 有经济价值的产量economic forecast 经济预测economic gain 经济收益economic geography 经济地理economic geology 经济地质学economic growth rate 经济增长率economic haul 经济运距economic index 经济指数economic indicators 经济指标;经济指数economic integration 经济一体化economic interests 经济利益economic internal rate of return 经济内部收益率economic leverage 经济杠杆economic life 经济寿命;经济开采期限economic limit rate 经济极限产量economic limit 经济极限economic loss 经济损失economic mathematical model 经济数学模型economic model 经济模式economic movement mechanism 经济运行机制economic net present value 经济净现值economic order quantity 经济定购量economic order 经济秩序economic phenomenon 经济现象economic pipe size 经济管径economic policy 经济政策economic potential 经济潜力economic production life 经济开采寿命economic projection 经济预测economic prosperity 经济繁荣economic rationality 经济合理性economic recession 经济衰退economic retrechment 经济紧缩economic return 经济收益economic sanction 经济制裁economic scale 经济规模economic selenology 经济月质学economic storage 经济库容economic structure 经济结构economic system reform 经济体制改革economic system 经济体制;经济制度economic tanker size 经济船型economic target system 经济指标体系economic target 经济指标economic thickness 经济厚度economic trend 经济趋势economic value 经济价值;工业价值economic water-oil ratio 经济水-油比economic worth 经济价值economic yardstick 经济标准;经济指标economic yield 经济产量economic zone 经济区economic 经济的economical 经济的economically recoverable oil 经济上有开采价值的石油economics of energy 能源经济学economics of management 管理经济学economics of scale 规模经济学economics 经济学economill 轻便钻头economist 经济学家economizer bank 预热管;节热器排管economizer 节油器economy 经济;经济体系econtone 生态混合群落;生态过渡区;溶泌区ECOSOC 经济及社会理事会ecosphere 生态层;生物域;大气层ecostratigraphic classification 生态地层分类ecostratigraphic unit 生态地层单位ecostratigraphic 生态地层学的Ecostratigraphy 生态地层学ecosystem 生态系统ecotelemetry 生态遥测术ecotonal community 交错区群落ecotope 生态区ecotopic 适应特殊生态的ecotype 生态类型ecoulement 重力滑动ECP 套管外封隔器ECP 误差校正剖析ECP 岩心穿透深度ECP 有效岩心穿透深度ECP 有效岩心孔隙度ECPL 工程零部件清单ECPWS 工程更改建议说明书ECR 工程更改申请ECR 工程控制室ECR 设备更改申请ecronic 河口湾ECS 环境控制系统ECS 扩展磁心存储器ECS 末端电池转换开关ectexis 泌出混合岩化作用ECTL 射极耦合晶体管逻辑ectype 复制品;副本ed. 版ed. 编辑edaphic control 底土控制edaphic 土壤的;土壤圈的EDC 程序超过磁鼓容量EDC 工程设计改变EDC 工程图纸改变EDC 估计完工日期eddy conductivity 涡动传导性eddy current brake 涡流闸eddy current coefficient 涡流系数eddy current constant 涡流常数eddy current coupling 电磁离合器eddy current damping 涡流阻尼eddy current inspection 涡流探伤eddy current seismometer 涡流地震检波器eddy current test 涡流探伤eddy current thickness meter 涡流测厚仪eddy current 涡流;涡电流eddy diffusion 涡流扩散eddy diffusivity 涡流扩散系数eddy effect 涡流效应eddy flux 涡流通量eddy marking 涡流痕迹eddy motion 涡流运动eddy resistance 涡流阻力eddy shedding 涡流分离eddy stress 涡流应力eddy thermal conductivity 涡流热传导eddy viscosity 涡流粘度eddy wind 旋风eddy 涡流eddy-current loss 涡流损耗eddy-current type geophone 涡流式检波器eddying flow 涡流edel metal 贵金属Edenian 艾登阶edenite 浅闪石edetate 乙二胺四乙酸盐edetic acid 乙二胺四乙酸edeyen 沙丘沙漠EDF 经验判别函数EDF 能量密度函数edge angle 棱角edge aquifer 边部水体edge away 楔出;尖灭edge connector 插头座edge contact 边缘接触edge crack 边缘裂隙edge damage 破边;边伤edge effect 边缘效应edge enhancement 边缘增强edge flare 卷边对接edge fog 边缘模糊edge gradient 边缘梯度edge joint 边缘连接;卷边接头;边缘焊接头edge lease 油气田边缘租地edge orientation 晶棱定向edge orifice 锐边孔板edge pressure 边缘压力edge radius 棱角半径edge reflection 边界反射edge resolution 边缘分辨率edge seam 边缘线状裂纹;倾斜地层edge stress 边缘应力edge tool 削边刀edge value 边界值;油水边界异常值edge water drive reservoir 边缘水驱油藏edge water drive 边水驱edge water flood 边缘注水edge water incursion 边水侵入edge water limit 边水界限edge water line 边水线edge water pressure 边水压力edge water 边水edge wave 边缘波edge weld 角接焊edge well 边井edge zone 边缘地带edge 边EDGE 电子数据收集设备edge-frame 内图廓edge-to-edge contact 边-边接触edge-to-surface contact 边-面接触edgefold 折边edger 轧边机;磨边器;弯曲膜瞠edgewise breccia 竹叶角砾岩edgewise conglomerate 竹叶砾岩Edgeworth expansion 埃奇沃斯展式edging 边缘;滚压边缘修饰;镶边;轧边;磨边Ediacara fauna 埃迪卡拉动物群edible clay 可食粘土edible oil 食用油edible 适合食用的edifice 大建筑物edingtonite 钡沸石EDIPS 地球资源观测系统数字图象处理系统EDIS 勘探资料解释系统Edison effect 爱迪生效应edit control character 编辑控制符edit mode 编辑方式edit routine 编辑程序edit 初步整理;编排;校定edit-modify 编辑修改Editia 美脊介属editing operation 编辑操作editing subroutine 编辑子程序editing 编辑edition 版本editor in chief 主编editor 编辑editorial committee 编委会editorial 编辑的editorship 编者职位EDM 电子距离测量EDM 放电加工EDMI 电磁波测距仪Edmond's balance 爱德蒙气体比重天平EDP 电子数据处理EDP 勘探数据处理EDPC 电子数据处理中心EDPE 电子数据处理设备EDPM 电子数据处理机EDPS 电子数据处理系统EDR 估算损害比EDR 实验数据记录EDS 编辑数据选择EDS 能量色散X-射线探测仪EDST 长期钻柱测试EDTA reagent 乙二胺四乙酸试剂EDTA titration 乙二胺四乙酸滴定法EDTA 乙二胺四醋酸EDTCC 电子数据传输通信中心education 教育educt 离析物eduction column 气举管eduction pipe 气举管;排泄管eduction tube 气举管;排泄管eduction valve 排气阀;泄流阀eduction 推断;析出eductor 气举管;喷射器edulcoration 冼净EDVAC 电子离散变量自动计算机Edward balance 爱德华天然气比重天平edwardite 独居石EE 电机工程EE 电气工程师EE 电眼ee 允许误差;错误不在此限EEC 欧洲经济共同体EECL 射极-射极耦合逻辑EELS 电子发射器定位系统EERI 地震工程研究学会EEROM 电可擦只读存储器EEZ 专属经济区EF 经验公式EF 侵蚀因素EF 倾角指示器EF 震源EFD 工艺流程图EFDL 射极跟随器二极管逻辑EFDTL 射极跟随器二极管-晶体管逻辑eff 效率eff 有效的effect of anisotropy 各向异性效应effect of dragging 拖曳作用effect of inertia 惯性作用effect of relative project 相关项目效果effect of wellbore storage 井筒储存效应effect 作用effective absorption coefficient 有效吸收系数effective address 有效地址effective age 有效使用期effective anisotropy 有效各向异性effective antenna length 有效天线长度effective aperture 有效孔径effective area 有效面积effective array length 有效组合长度effective atmosphere 有效大气压effective atomic number 有效原子序数effective attenuation 有效衰减effective bit-weight 有效钻压effective brake area 有效制动面积effective capacitance 有效电容effective capacity 有效容量effective capture cross-section 有效俘获截面effective clearance 有效间隙effective competition 有效竞争effective compressibility 有效压缩系数effective concentration 有效浓度effective constant 有效常数effective constituent 有效成分effective contact radius 有效接触半径effective core penetration 有效岩心穿透深度effective core porosity 有效岩心孔隙度effective coross-section 有效截面effective current 有效电流effective data transfer rate 有效数据传输率effective date 生效日期effective decline rate 有效递减率effective demand 有效需求effective depth 有效深度effective diameter 有效直径effective discharge 有效排量effective dose 有效剂量effective field intensity 有效场强effective film 有效界膜;膜的有效厚度effective grain-size 有效粒径effective half-life 有效半衰期effective head 有效压头effective heat duty 有效热负荷effective heating time 有效加热时间effective height 有效高度effective instruction 有效指令effective length 有效长度effective life 有效寿命;有效使用期effective management 有效管理effective matrix parameter 有效骨架参数effective molecular weight 有效分子量effective multiplication factor 有效放大系数;有效增殖系数effective noise temperature 等效噪声温度effective output 有效产量effective pay factor 生产层有效因素effective pay thickness 产层有效厚度effective pay 有效产层effective percentage modulation 有效调制深度effective perforation 有效炮眼effective permeability 有效渗透率effective placement 有效充填effective porosity 有效孔隙度effective power 有效功率;水力功率effective pressure head 有效压头effective pressure 有效压力effective price 有效价格effective profile 有效剖面effective radiation 有效辐射effective radius 有效半径effective range 有效范围;有效测程effective recourse 有效追索effective resistance 有效电阻effective rock volume 有效岩石体积effective service life 有效使用期effective shot depth 激发有效深度effective size of grain 有效粒径effective sound pressure 有效声压effective source rock 有效生油岩effective stress 有效应力effective surface 有效表面effective susceptibility 有效磁化率effective temperature 有效温度effective terrestrial radiation 有效大地辐射;有效地面辐射effective thickness 有效厚度effective thread 有效螺纹effective time 有效时间effective transmission rate 有效传输率effective up-date rate 有效更新率effective value 有效值effective velocity 有效速度effective vibration length 有效振动距离effective viscosity 有效粘度effective volatility 有效挥发度effective water resistivity 水的有效电阻率effective water saturation 有效含水饱和度effective wavelength 有效波长effective word 有效字effective work 有效功;有效工作effective 有效的effectiveness 效率effector 效应器官;操纵器;格式控制字符effectuation 有效化;完成effervesce 起泡effervescence 起泡effervescent 起泡的effervescing steel 沸腾钢efficacy 效力efficiency coefficient method 功效系数法efficiency contract 效率合约efficiency curve 效率曲线efficiency diagram 效率图efficiency diode 增效二极管;阻尼二极管efficiency distribution coefficient 有效分布系数efficiency estimate 有效估计量efficiency extraction 有效提取efficiency factor 效率因子efficiency of conversion 转换效率efficiency of core penetration 岩心穿透深度efficiency of displacement 驱替效率efficiency of estimator 估计量的有效性efficiency of generator 发电机效率efficiency of imbibition 自吸效率efficiency of labour 劳动生产率efficiency of management 管理效率efficiency of pump 泵效率efficiency of rectification 整流效率efficiency ratio 效率比efficiency statistic 有效统计量efficiency test 有效试验efficiency value 效率值efficiency 效率efficient departmentalization 有效分权efficient 有效的effictiveness 效果efflatum 火山喷出物efflorescence 风化effluence 射出effluent channel 排泄道effluent clarity 流出液清洁度effluent concentration 流出物浓度effluent control 流出物控制effluent face 流出面effluent gases 排放气effluent holding reservoir 污水储池effluent oil recovery 污油回收effluent seepage 渗流effluent specification 排放规范effluent standard 排放标准effluent stream 外排流effluent treatment 流出物处理effluent 流出effluent-disposal standard 污水处理标准effluent-end saturation 流出端饱和度effluve 高压放电effluvia effluvium的复数effluvium 无声放电;磁素;臭气;散出efflux coefficient 排出系数efflux cup method 流杯法efflux time 排出时间efflux velocity 排出速度efflux 流出;射流;泄漏;时间消逝effort 努力;成果;工作effuse 喷出effuser 喷管;扩散器;集气管effusiometer 气体扩散计effusion 喷发effusive breccia 喷出角砾岩effusive eruption 溢流喷发effusive mass 喷发体effusive period 喷发期effusive rock 喷发岩effusive 射流的effy 效率EFL 错误频率极限eg exempli gratia 例如EG 甘醇egg end 半球形的底板egg shaped 蛋形的egg 蛋;卵形物eggbeater PDC bit 打蛋器型PDC钻头EGL 地面海拔高度EGMBE 乙二醇一丁醚EGP 处理事故的砾石充填EGR 废气再循环egress and ingress 出入egress hole 出口孔egress of heat 热传导egress pressure 出口压力egress 出口EGS 欧洲地球物理学会EGT 废气温度Egyptian General Petroleum Corp. 埃及石油总公司EH control 电动-液压控制EHF 极高频EHP horse-power) 电马力EHP 有效马力EHPH 马力-小时ehrwaldite 玻基二辉岩EHT 极高压EHV 极高压EI 电子撞击EI 工程索引EI 能量指数EI 侵蚀指数Ei-function Ei函数EIA 地址错误EIA 电子工业协会EIA 工程工业协会EIA 美国能源情报局EIC 工业情报中心eicosane 二十烷eidograph 缩放仪eidophor 大图象投射器Eifelian 艾斐尔阶eigen space 本征空间eigen 本征的;固有的eigen- 本征eigenfrequency spectrum 固有频谱eigenfrequency 本征频率eigenfunction 本征函数eigenmode 正则型;本征型eigenperiod 本征周期eigenstate 特征状态eigentone 本征音;固有振动频率eigenvalue extraction 特征值析取eigenvalue 本征值eigenvector 特征矢量eigenvibration 本征振动eigenwavefront 特征波前eigenwert 特征值eight digit number 八位数eight-digit binary number 八位二进制数eight-point mooring system 八点系泊系统eighty-board 井架工作平台eighty-column card 八十列卡片eighty-column puncher 八十列穿孔机eigram 双字母组eigth bit analog 八位模拟eikonal equation 程函方程eikonogen 影源eikonometer 光象测定器EIL 符号部分出错einstein 爱因斯坦Einstein's equation 爱因斯坦方程einsteinium 锿Einthoven galvanometer 弦线电流计EIO 操作码出错EIP 电激发极化法EIPS 特优犁钢EIR 环境影响报告EIS 环境影响报告EIU 内外加厚的EIW 电感应焊ejaculation 突发eject 喷出ejecta 喷出物ejected electron 发射的电子ejected photoelectron 发射的光电子ejected rock 喷出岩ejection curve 退出曲线ejection efficency 退出效率ejection nozzle 喷嘴ejection test 喷射试验ejection 喷出ejective fold 隔挡褶皱ejector air pump 喷射空气泵ejector condenser 喷射冷凝器ejector pin 推顶杆ejector plate 推顶杆板ejector priming 喷射泵启动ejector pump 喷射泵ejector return pin 复位杆ejector rod 推顶柱ejector sleeve 推顶套ejector vacuum pump 喷射真空泵ejector 射流泵ejector-type through-tubing tool 喷射式过油管下井仪eka- 准eka-element 准元素eka-silicon 准硅EKB 方钻杆补心海拔高度ekerite 钠闪花岗岩ekistics 城市与区域计划学Ekman dredge 艾克曼采泥器eksedofacies 风化环境相ektexic 泌出变熔作用ekzema 盐穹EL 弹性极限EL 电测井EL 电致发光el 高地;海拔;正视图EL 设备表elaborate 精心制成的elaboration 精心制作elaeometer 验油比重计elapse 过去elapsed time 经过的时间elastance 倒电容elastic absorption 弹性吸收elastic acoustical reactance 弹性声阻抗elastic after effect 弹性后效elastic after-working 弹性后效elastic anisotropy 弹性各向异性elastic beam 弹性梁elastic behavior 弹性特性elastic bitumen 弹性沥青elastic body 弹性体elastic boundary 弹性边界elastic break-down 弹性失效elastic buffer 弹性缓冲器elastic coefficient 弹性系数elastic collision 弹性碰撞elastic compaction 弹性挤压作用elastic connector 弹性连接器elastic constant 弹性常数elastic coupling 弹性联轴节elastic deformation 弹性变形elastic demand 弹性需求elastic discontinuity 弹性不连续性elastic drive 弹性驱动elastic effect 弹性效应elastic elongation 弹性伸长elastic energy 弹性能elastic equilibrium 弹性平衡elastic extension 弹性延伸elastic fatigue 弹性疲劳elastic fibre 弹性纤维elastic flow 弹性流elastic fluid 弹性流体elastic force 弹性力elastic formation 弹性地层elastic gel 弹性凝胶elastic half-space 弹性半度空间elastic heterogeneity 弹性非均匀性elastic hysteresis 弹性滞后elastic impedance 弹性阻抗elastic imperfection 弹性不完整elastic instability 弹性不稳定性elastic isotropy 弹性各向同性elastic joint 弹性连轴节elastic lag 弹性滞后;弹性惯性elastic limit 弹性极限elastic liquid 弹性液体elastic medium 弹性介质elastic migration theory 弹性波偏移理论elastic mineral pitch 弹性沥青elastic module 弹性系数elastic modulus 弹性模量elastic nylon 弹性锦纶elastic peak 弹性峰值elastic plan 弹性计划elastic potential theory 弹性势理论elastic properties 弹性elastic range 弹性范围。

The Cores of Dark Matter Dominated Galaxies theory vs. observations

The Cores of Dark Matter Dominated Galaxies theory vs. observations

a r X i v :a s t r o -p h /9708176v 2 25 M a y 1998to be published in July 20,1998issue of the Astrophysical JournalPreprint typeset using L A T E X style emulateapj v.04/03/99THE CORES OF DARK MATTER-DOMINATED GALAXIES:THEORY VERSUS OBSERVATIONSAndrey V.Kravtsov and Anatoly A.KlypinAstronomy Department,New Mexico State University,Las Cruces,NM 88003-0001,USAandJames S.Bullock and Joel R.PrimackPhysics Department,University of California,Santa Cruz,CA 95064to be published in July 20,1998issue of the Astrophysical JournalABSTRACTWe use the rotation curves of a sample of dark matter dominated dwarf and low-surface brightness (LSB)late-type galaxies to study their radial mass distributions.We find that the shape of the rotation curves is remarkably similar for all (both dwarf and LSB)galaxies in the sample,suggesting a self-similar density distribution of their dark matter (DM)halos.This shape can be reproduced well by a density profile with a shallow central cusp [ρ(r )∝1/r γ,γ≈0.2−0.4]corresponding to a steeply rising velocity curve [v (r )∝r g ,g ≈0.9−0.8].We further show that the observed shapes of the rotation curves are well matched by the average density profiles of dark matter halos formed in very high resolution simulations of the standard cold dark matter model (CDM),the low-density CDM model with cosmological constant (ΛCDM),and the cold+hot dark matter model with two types of neutrino (CHDM).This is surprising in light of several previous studies,which suggested that the structure of simulated dark matter halos is inconsistent with the dynamics of dwarf galaxies.We discuss possible explanations for this discrepancy and show that it is most likely due to the systematic differences at small radii between the analytic model proposed by Navarro,Frenk,&White,with γNFW =1,and the actual central density profiles of the dark matter halos.We also show that the mass distributions in the hierarchically formed halos are on average consistent with the shape of rotation curves of dark matter dominated galaxies.However,the scatter of the individual profiles around the average is substantial.Finally,we show that the dark matter halos in our hierarchical simulations and the real galaxies in our sample exhibit a similar decrease in their characteristic densities with increasing characteristic radial scales and show increase in their maximum rotation velocities with increase in the radii at which their maximum velocities occur.Subject headings:cosmology:theory –dark matter:halos —galaxies:kinematics and dynamics —galaxies:structure1.INTRODUCTIONThe amount of luminous matter (stars and gas)in many spiral and irregular galaxies is not sufficient to explain the amplitude and shape of their rotation curves (RCs).This discrepancy is usually interpreted as evidence for the presence of an extended dark matter (DM)halo surround-ing the visible regions of galaxies (e.g.,Casertano &van Gorkom 1991;Persic,Salucci,&Stel 1996,and references therein).The extent of the dark matter halos,estimated using satellite dynamics,is 1∼0.2−0.5h −1Mpc (Zarit-sky &White 1994;Carignan et al.1997;Zaritsky et al.1997).However,the dynamical contribution of the dark matter can be substantial even in the very inner regions of galaxies:the observed rotation velocities of some dwarf and low-surface brightness (LSB)galaxies imply that DM constitutes a dominant fraction (up to ∼95%)of dynami-cal mass within the last measured point of their RCs (e.g.,Carignan &Freeman 1988;Martimbeau,Carignan,&Roy 1994;de Blok &McGaugh 1997).These dark matter dom-inated galaxies offer a unique opportunity for probing di-rectly the density structure of DM halos which can be then compared with predictions of theoretical models.The detailed structure of DM halos formed via dissipa-tionless hierarchical collapse in CDM-like models was re-cently studied using high-resolution N -body simulations (Dubinski &Carlberg 1991;Navarro,Frenk,&White 1996,1997,hereafter NFW96and NFW97).The halo density profiles were found to be cuspy (coreless)and well fitted by the following two-parameter profile (NFW96):ρNFW (r )=ρs2KRAVTSOV ET AL.B95),who pointed out that shapes of the density profilesof four dwarf galaxies analyzed by Moore(1994)are re-markably similar and are wellfitted by the following phe-nomenological density profile:ρbρB(r)=.(3)(r/r0)γ[1+(r/r0)α](β−γ)/αNote thatρ(r≪r0)∝r−γ,ρ(r≫r0)∝r−β,andαcharacterizes the sharpness of the change in logarithmicslope.This family includes both cuspy profiles of the typeproposed by NFW96(α,β,γ)=(1,3,1)and the so-calledmodified isothermal profile(α,β,γ)=(2,2,0),which isthe most widely used model for the halo density distri-bution in analyses of observed rotation curves.It is alsoconvenient to make directfits with an analytic model sim-ilar to(3)for the velocity profile:(r/r t)gV(r)=V tln r/r for large r,which has an approximate slope of b∼0.34for values of r near a typical virial radius.CORES OF DARK MATTER-DOMINATED GALAXIES3Table1The sample of dwarf and LSB galaxiesr0V0ρ0r t V t r max V max DistanceGalaxy M B h−1kpc km s−1108h3M⊙kpc−3h−1kpc km s−1h−1kpc km s−1h−1Mpc Reference (1)(2)(3)(4)(5)(6)(7)(8)(9)(10)(11)NOTES.–Col.(2)M B,blue absolute magnitude;col.(3)bestfit r0(see eq.[3];thefitting procedure is described in§2.2); col.(4)V0=V(r0);col.(5)bestfitρ0(see eq.[3]);col.(6)bestfit r t(see eq.[4]);col.(7)bestfit V t(see eq.[4]);col.(8)r max; col.(9)V max=V(r max);col.(10)distance to galaxy adopted in this study;REFERENCES.–(1)Carignan&Freeman1988;(2)Carignan&Beaulieu1989;(3)Lake et al.1990;(4)Meurer et al. 1996;(5)Martimbeau et al.1994;(6)Cˆo t´e et al.1991;(7)Jobin M.&Carignan C.1990;(8)van Zee et al.1997;(9)de Blok et al.1996.4KRAVTSOV ET AL.constrained for all of the galaxies.However,a fair number of galaxies in the sample do show the turnover and thus can be used to constrain α.The plausible value of the parameter α=2was determined using rotation curves of these galaxies.We generalize this value to all of the galax-ies (which in no way contradicts the data,but is not,of course,a strict procedure)and thus currently we can only talk about a plausible range of αvalues as a “universal fit”(if any such universal value exist at all).For example,our results will not change drastically if we use α=1.5instead of α=2.However,α=1gives a poorer fit to the data.We fix the parameter γto 0.2:the value which best fits most of the observed rotation curves.The corresponding best-fit slopes of the profile (4)are (a,b,g )=(1.50,0.34,0.9).Note that g =1−γ/2.With parameters α,β,and γ(a,b,g )fixed,we fitted the data for the remaining free parameters of the profile (3):ρ0and r 0(V t and r t in eq.[4]).Our fits thus have the same number of free parameters as do profiles (1)and (2).Note that while the particular set of the parameters (α,β,γ)=(2,3,0.2)used in the paper didn’t result from a strict fitting proce-dure,it was motivated by all possible constraints of avail-able data.The only theoretically suggested value is that of βbut it also seems to be favored by data (Burkert 1995).Hopefully,as new RC observations come along,they can be used to pinpoint the parameters αand βwith a betteraccuracy.Fig.1.—Rotation curves of (a)and (b)LSB galaxies (symbols)normalized to the best fit values of r 0and rota-tional velocities v 0at r 0predicted by density model (eq.[3]).The solid line on both panels shows the analytic rotation curve corresponding to the density profile (eq.[3])with (α,β,γ)=(2,3,0.2).The rotation curves for different dwarf and LSB galaxies have virtually identical shapes,which is very well matched over the entire observed range of scales by the ana-lytic model.Note,that the RC of NGC 2915extends outside the scale of the plot:the outer part of this RC can be seen in Figure 2.Figure 1shows rotation curves of dwarf (a)and LSB (b)galaxies normalized to their best fit values of r 0and to the rotational velocities v 0at r 0,predicted by analytic profile (3).The best fit values of r 0,ρ0,V 0=V (r 0),r t ,and V t are given in Table 1for each galaxy in the sam-ple.The formal errors of each of these values are less than about ∼2−5%.Figure 1a shows that all of the dwarf galaxies have rotation curves of virtually identical shape 3with a remarkably small scatter.The rotation curves of the two dwarf galaxies,DDO154and NGC2915,cannot be described by a smooth density distribution model in their outer parts.The RC of DDO154shows a decrease in ro-tational velocity in the three outermost observed points.Conversely,the RC of NGC2915has a sharp upturn at∼>5h −1kpc (or r/r 0∼>4.5).This upturn can be seen in Figure 2.The explanation of this peculiar behavior is not clear (see,however,Burkert &Silk 1997),but it is obvious that it cannot be explained by any smooth model for the mass distribution.Note,however,that apart from the pe-culiar outer regions,the rotation curves of both DDO154and NGC2915have the same shape as the rest of the galax-ies.The shape of the galaxies’rotation curves is well matched by the rotation curve corresponding to density profile (3)with (α,β,γ)=(2,3,0.2)or correspondingly to RC (4)with (a,b,g )=(1.5,0.34,0.9).This result is in per-fect agreement with Burkert (1995)who showed similarity of rotation curves for four dwarf galaxies (two of which,DDO154and DDO170,were included in our sample).As was mentioned above,the similar fit by ρB (r )(eq.[2])proposed by Burkert (1995)is equally good.Note,how-ever,that our profile does not have any flat core,whereas ρB (r )predicts such a core at r ≪r b .The fact that both profiles fit the data equally well is easy to understand if we notice that ρB (r )predicts a flat density distribution at the scales well below the observational resolution (∼<1kpc).Thus,ρB (r )and profile (3)can be virtually identical in the range of scales resolved in observations and thus provide an equally good fit to the data.Figure 1b shows that the rotation curves of dark mat-ter dominated LSB galaxies are also well described by the same analytic density profile.The larger amplitude of scatter in the case of LSB galaxies can be explained by the larger observational errors associated with a given point of a rotation curve and thus most likely reflects ob-servational uncertainties rather than intrinsic scatter of the halo properties.Most of the LSB galaxies in our sam-ple are located at considerably larger distances than dwarf galaxies.Therefore,the dwarf galaxies have been observed with considerably higher resolution and smaller observa-tional errors than LSB galaxies.The estimated errors are typically 10−20%(de Blok et al.1996),especially in the inner regions of galaxies (∼<10kpc).One important issue is whether subtraction of baryon component (stars and gas)in the galaxies from our sample will affect results of the rotation curve analysis presented above.As we mentioned in §2.1,the combined contribu-tion of stars and gas is ∼<15%for all of our galaxies (∼<10%in most cases).It is clear that ideally one has to subtract3Bythe shape of a rotation curve we mean its particular functional form.For example,the shape of the rotation curve described by equation(4)is x g /(1+x a )(g +b )/a (where x ≡r/r t ).By saying that the RC shape is similar for all our galaxies,we mean that all their rotation curves can be described by this functional form with fixed values of parameters a ,b ,and g .CORES OF DARK MATTER-DOMINATED GALAXIES5contributions of both gas and stars from the observed RC in order to get the mass distribution of the dark matter to an accuracy of better than10%.However,it is well known that this is not an easy thing to do.For stars,the exact mass-to-light ratio is not known and we cannot convert vis-ibleflux into the stellar mass without making additional assumptions(hence,“maximum disk controversy”).In the case of gas,we know exactly how to convert the21-cmflux into the mass of gas(no mass-to-light uncertainty).This conversion,however,relies on other assumptions which can easily lead to uncertainties as high as10%in the contri-bution of gas.For example,we need to know distance to the galaxy in order to make this conversion.We also need to have a reliable way to estimate the profile of molecu-lar gas to recover the dark matter density profile.The distances to the dwarf galaxies in our sample are very un-certain(often by a factor of two or more)and so is the conversion.It is not clear whether it makes any sense to subtract the gas with uncertainties of its contribution this big.After all,the motivation for use of dark matter domi-nated galaxies for this kind of analysis is to avoid dubious or uncertain correction procedures which are very unlikely to result in a better determination of the shape of DM density distribution.The distance uncertainty is not so severe for LSB galaxies and so conversion could,in princi-ple,have been done in this case.We have not done this for one simple reason:subtraction of the gas component could change any given point of rotation curve by at most10% (in practice less than that).However,the errors associated with each point of RC are of the order of10−20%(which combine both observational errors and assymetries in the rotation curves between receding and approaching sides) and it seems unlikely that correction due to gas subtrac-tion would improve or systematically change the answer (unless the observational errors are significantly overesti-mated).We have tested the effect of gas subtraction on the RC shapes by performing RC shape analysis for two galaxies(dwarf NGC5585and LSB F583−1)with and without subtraction of gas.The two galaxies have been se-lected to have a clearly visible turnover of the RC and to have a fairly high fraction of gas inside the last measured point of the rotation curve.This fraction is8%for NGC 5585(Cˆo te et al.1991)and5%for F583−1(de Blok et al. 1996).The results offitting the(α,β,γ)=(2.0,3.0,0.2) model to the RC with and without gas subtraction result in very similar results:the difference in the bestfit param-eters is∼<10%and corrected and uncorrected RCs have virtually identical shape.Note also that the dwarf galaxies that we used have on average a higher(or at least as high) fraction of gas(typically∼6−10%)than LSB galaxies (typically3−7%).Therefore,if gas would introduce sys-tematic differences in the shape of RC,we could expect that the scatter in Figure1to be larger for dwarf galax-ies(some galaxies have much more gas than the others). Yet,the shape of rotation curves for dwarf galaxies is very uniform.We have repeated thefitting procedure described above using the analytic profiles(1)and(2).As was mentioned above,ρB(r)results in afit that is equally good to thefit by profile(3)shown in Figure1.However,the analytic profile proposed by NFW failed to produce a reasonable fit to the data,as was indeed pointed out in NFW96(see theirfig.12).The major difficulty with this profile,as was noted before by Flores&Primack1994and B95,is that the inner slope of the density distribution(γ=1)is con-siderably steeper than implied by the rotation curves.The finite spatial extent of the data and incorrect inner slope of the profile(1)lead to implausible solutions of theχ2-minimization procedure(the values of r s increase without convergence).The observed similarity of the shapes of the rotation curves for seventeen different galaxies,selected solely on the basis of their dark matter content,and the remarkably small amount of scatter,implies that their matter distri-butions are self-similar in terms of the density structure. Of course,this includes both stellar and gaseous matter as well as DM.Both stellar and gaseous masses are uncer-tain because of uncertainties in the distance,mass-to-light ratio,and atomic-to-molecular gas ratio of each galaxy. To the extent that we can neglect the stellar and gaseous components(a subject that we intend to address in a sub-sequent paper),the self-similar rotation curves of these DM-dominated galaxies imply that they all have the same density structure.The question we now ask is whether the disagreement between this density structure andρNFW(r) indicates a failure of CDM-type models?PARISON WITH THEORETICAL MODELS3.1.Numerical simulationsWe have used the new Adaptive Refinement Tree(ART) N-body code(see Kravtsov,Klypin,&Khokhlov1997 for details)to simulate the evolution of collisionless dark matter in the three cosmological structure formation mod-els:(1)standard cold dark matter model(CDM:Ω0=1, h=0.5,σ8=0.7);(2)a low-density CDM model with cos-mological constant(ΛCDM:ΩΛ=0.7,h=0.7,σ8=1.0); and(3)a cold+hot dark matter model with two types of neutrino(CHDM;Ω0=1andΩν=0.2;h=0.5;σ8=0.7;cf.Primack et al.1995).HereΩ0,ΩΛ,andΩνare the present-epoch values of the density of matter,vac-uum energy(as measured by the cosmological constant), and massive neutrinos,respectively.The rmsfluctuation in spheres of radius8h−1Mpc,σ8,was chosen to conform with the local abundance of galaxy clusters,forΛCDM and CHDM models it is also in agreement with measure-ments of the cosmic microwave background anisotropy by the COBE satellite.The simulations followed trajecto-ries of1283cold dark matter particles in a box of size of L box=7.5h−1Mpc.In the CHDM simulation,two addi-tional equal-mass“massive neutrino”species were evolved, which brings the number of particles in the simulation to 3×1283.To test for the possible effects of thefinite box size,we have run an additional simulation of theΛCDM model with the box size twice as large:L box=15h−1 Mpc=21.43Mpc.We will denote the twoΛCDM simu-lations asΛCDM7.5andΛCDM15according to their box sizes.We have used a2563uniform grid covering the entire computational volume andfiner refinement meshes con-structed recursively and adaptively inside the high-density regions.The comoving cell size corresponding to a refine-ment level L is∆x L=∆x0/2L,where∆x0=L box/256 is the size of the uniform grid cell(L=0corresponds to the uniform grid).The increase of spatial resolution corre-sponding to each successive refinement level was accompa-6KRAVTSOV ET AL.nied by the decrease of the integration time step by a factor of2.The simulations were started at redshift z i=40in the CDM andΛCDM7.5simulations and at z i=30in the CHDM andΛCDM15simulations.Particle trajecto-ries were integrated with the time step of∆a0=0.0015 on the zeroth-level uniform grid in the case of the CDM andΛCDM runs and with∆a0=0.006in the CHDM run. The time step on a refinement level L is∆a L=∆a0/2L. The time step for the highest refinement level corresponds to∼>40,000time steps over the Hubble time.Six refine-ment levels were introduced in the highest density regions corresponding to a cell size of∆x6=0.46h−1kpc.The dynamic range of the simulations is thus256×26=16,384. Note,that the resolution is constant in comoving coordi-nates which means that actual physical resolution is higher at earlier epochs(the halos were resolved with six refine-ment levels as early as z≈1).The refinement criterion was based on the local overdensity of dark matter particles. Regions with overdensity higher thanδ=n th(L)23(L+1) were refined to the refinement level L.Here,n th(L)is the threshold number of particles per mesh cell of level L estimated using the cloud-in-cell method(Hockney& Eastwood1981).We have used values n th=5at all levels in the CDM andΛCDM runs;for the CHDM run we have used n th=10at the levels L=0,1and n th=5for all of the higher levels.These values of threshold were sug-gested by results of the tests presented in Kravtsov et al. (1997);they ensure that refinements are introduced only in the regions of high-particle density,where the two-body relaxation effects are not important.For the dark matter halos used in our analysis the spa-tial resolution is equal to≈0.5−2h−1kpc(corresponding to the6th to4th refinement levels).For each of the ana-lyzed halos,we have taken into account only those regions of the density and circular velocity profiles that correspond to scales at least twice as large as the formal resolution. The mass resolutions(particle mass)of our simulations are listed in Table2,and are in the range of∼(1−10)×107h−1 M⊙.Therefore a typical halo of mass∼1011h−1M⊙in our simulations contains several thousands of particles. These simulations are comparable in spatial and mass resolution,as well as in the box size,to those of NFW96,97. There is,however,a significant difference:our simulations are direct simulations of all DM halos in a given compu-tational volume,whereas NFW96,97simulate with high resolution a handful of individual halos.The fact that we analyze a statistically large sample consisting of dozens of galaxy-size halos in each simulation allows us to make con-clusions about average halo properties and estimate the amount of cosmic scatter.A summary of the numerical simulations is given in Table2.The parameters listed in this table are defined in the text above.3.2.Tests of numerical effectsThere are several effects which can affect the halo den-sity profiles at scales larger than some particular scale re-lated to this effect.We have tested the reliability of the simulated density and velocity profiles by comparing re-sults of the simulations with different resolutions and time steps.Specifically,the tests were used to determine the range of numerical parameters for which the convergence of density profiles was found at scales larger than two for-mal resolution elements(formal resolution is equal to the size of the refinement mesh cell).Tests presented in Kravtsov et al.(1997)show that the density profiles are not affected by the force reso-lution down to a scale of about one resolution element (a similar conclusion was reached by NFW96).To test the effects of the time step we have used a set of643-particle simulations of the CDM model with parameters identical to those described in the previous section.These test simulations were started from identical initial con-ditions,but evolved with different time steps:∆a0= 0.006,0.003,0.0015,parison of the density profiles for the same halos in these simulations shows that for halos of all masses,the profiles converge for runs with ∆a0∼<0.0015(the value used in our CDM andΛCDM simulations)at all scales,down to the resolution limit. We further use two1283-particle simulations of theΛCDM model with the box size of15h−1Mpc and with time steps of∆a0=0.006and∆a0=0.0015.The comparison shows that the most massive halos(virial mass M vir>1013h−1 M⊙)have systematically shallower central(r∼<10−20h−1 kpc)density profiles in the∆a0=0.006run as compared to the halos from the∆a0=0.0015run.However,the dif-ference decreases with decreasing halo mass and for M vir ∼<5×1012h−1M⊙the density profiles from the two runs are identical within statistical noise at scales larger than one resolution element.This mass dependence is due to the different accuracy of numerical integration in objects of different masses.The accuracy depends on the average displacement of particles during a single time step:for the integration to be accurate,the displacement should be∼<10−20%of the resolution element.Particles in-side more massive halos have considerably higher veloci-ties(v∼>300−400km/s)and thus average displacements that are larger than the displacements of particles inside small halos(v∼<200km/s).In this study we focus on the mass distribution of the small halos(M∼<1×1012 M⊙),for which the tests indicate convergence of the den-sity profiles for time steps∆a0≤0.006.The time step of all simulations presented in this paper,except for the CHDM simulation,is four times smaller than the above value(see Table2).As an additional test,we have com-pared average RC shapes for CDM halos in the7.5h−1 Mpc box simulation shown in Figure2a and for halos in the same mass range(∼<1×1012M⊙)from an identical simulation(identical initial conditions and simulation pa-rameters)with time step∆a0=0.006.We have found that average RC shapes and the scatter in these two sim-ulations are indistinguishable.The mass resolution in our simulations(particle mass) is∼(0.6−5)×107h−1M⊙for L box=7.5h−1Mpc runs,and1.3×108h−1M⊙for the test L box=15h−1 MpcΛCDM run(see Table2).Therefore,halos of mass M vir=1012h−1M⊙and M vir=1011h−1M⊙(the range of masses used in our comparison with the data)are re-solved with∼100,000and∼10,000particles,respec-tively.For reference,there are∼>100−200particles inside the innermost point(2formal resolutions)of the rotation curve used in thefitting procedure described -parison of the average velocity profiles in theΛCDM7.5 andΛCDM15simulations(the latter has eight times worse mass resolution than the former)shows that there are no systematic differences between profiles in these two runsCORES OF DARK MATTER-DOMINATED GALAXIES 7Table 2Parameters of the Numerical SimulationsModelΩ0ΩΛΩνhσ8z i∆a 0L box Particle mass ×10−3h −1Mpc×107h −1M ⊙acold particles;bhot particles.(see Fig.2).The force resolution can introduce errors in rotational velocities.To estimate this effect,we assume that the fi-nite force resolution results in a flat core (ρ=const)at scales smaller than the resolution element h r in an oth-erwise ideal NFW halo (Eq.[1]).This results in the ve-locity profile v soft (r )/v s =F (x )/xF (1).The error is ∼18%at r ≈h r ,and ∼<5%at r ∼>2h r (see Fig.5).Thus,the ve-locity profiles of simulated halos should not be significantly affected at scales r ∼>2h r ,which is where we perform the fit to analytic models.To test whether the box size of our simulations (7.5h −1Mpc)is large enough not to miss all important tidal effects,we have compared the density and velocity profiles of halos formed in ΛCDM 7.5and ΛCDM 15simulations.We have not found any systematic differences between halo density profiles in these simulations.The average profiles of ha-los are identical within the statistical noise (see Figure 2).We have also used another indirect way of testing for the proper simulation of the tidal fields.Tidal torques from the surrounding large-scale structure presumably play a major role in the acquisition of the angular momentum,J =|J |,by the galaxy-size halos (Peebles 1969;Doroshke-vich 1970).Therefore,we can test if the tidal effects were simulated properly by comparing the so-called spin pa-rameter for the halos in our runs with previous results based on the larger-box simulations.The spin param-eter,λ,of a halo is defined as λ≡J |E |1/2/(GM 5/2),where J is the angular momentum of the halo,E is its total energy,and M is the halo virial mass.We have found that the distributions of λis very nearly log-normal 4,P (λ)=(1/λ√−2φ(r )≈2.15v max8KRAVTSOV ET AL.most dwarf and LSB galaxies are measured only to radii of ∼<10−30h−1kpc and often are still rising at the last mea-sured point.Therefore,the mass distribution in the outerparts of the galactic halos(and often maximum rotation velocity)is poorly constrained.To avoid any bias in the fitting procedure we considered only the inner30h−1kpc of the simulated halos.We then normalized each rotation profile to its bestfit values of r0and rotational velocity v0 at r0and computed the average of these normalized pro-files over all halos considered in each cosmological model (∼50−60).In Figure2we compare the average nor-malized dark matter velocity profiles for halos formed in CDM,ΛCDM,and CHDM models,shown by solid lines, with corresponding profiles of the dwarf galaxies from our sample,shown with different symbols(the symbols are as in Fig.1).The average profile from the larger-boxΛCDM15 simulation is shown with a dashed line.This profile does not extend to values of r/r0which are as low as for the ΛCDM7.5profile(due to worse spatial resolution).How-ever,for values of r/r0,where the two profiles overlap, they are indistinguishable.The dotted lines show the2σenvelope representing the scatter of individual halo profiles around the average.It should be noted that the scatter in the inner regions of the halo velocity profiles is substantial. This scatter possibly reflects physical differences between individual halos:our tests show that it is unlikely that the scatter can be attributed to the statistical noise associated with thefinite mass resolution.The mass resolution of our simulations is very high(see Table2):the number of dark matter particles inside the smallest scale,r min,of rota-tion curve used in thefitting procedure is∼>200for large (∼1012M⊙)halos and∼>100for smaller(∼1011M⊙) halos.Figure2shows that on average the velocity profiles of halos formed in hierarchical structure formation models and observed dark matter halos are in good agreement. It also shows that both cold dark matter halos and halos of dark matter dominated galaxies exhibit a certain self-similarity of the mass distribution in their inner regions.It was noted previously(e.g.,B95,NFW96)that hierarchi-cal formation of the halos should also result in well-defined scaling properties of the mass distribution.It is thus inter-esting to compare the scaling properties of galaxy halos in our sample with those of the DM halos formed in the three hierarchical models studied in this paper.Figure3shows the plot of the best-fit parameters r0andρ0of the model density distribution(3)for the dwarf(solid circles)and LSB(open circles)galaxies together with corresponding parameters of DM halos formed in CDM(a),ΛCDM(b), and CHDM(c)simulations.As before,the values of the remaining parameters of the profile(3)werefixed to(α,β,γ)=(2,3,0.2).For both galaxies and simulated halos, the parameters r0andρ0are clearly correlated:the ha-los that are compact are systematically denser.DM halos in all models are fairly consistent with the observational points,except possibly for the CDM model that appears to form halos somewhat denser than observed.Note that the absence of halos at r0∼<2h−1kpc is due to ourfinite nu-merical resolution rather than the generic failure of these models to produce very compact halos.The characteristic density of the DM halos correlates strongly with halo mass in a way that reflects the mass dependence of the epoch of halo formation(NFW96):low-mass small halos collapse at systematically higher redshifts(when the universe was denser)and are therefore denser than the larger higher-mass halos.Thus,the correlation observed in Figure3is likely to reflect the different formation epochs of individualhalos.Fig.2.—Average normalized dark matter velocity profiles for halos formed in(top panel)CDM,(middle panel)ΛCDM, and(bottom panel)CHDM models with corresponding profiles of the dwarf galaxies from our sample.The dotted lines show the2σenvelope representing scatter of individual halo profiles around the average.It should be noted that although the ve-locity profiles of the hierarchically formed dark matter halos are on average consistent with the shape of observed rotation curves,the scatter in the inner regions of the halo velocity profiles is substantial.This scatter possibly reflects real phys-ical differences between individual halos.The average profile from the larger-boxΛCDM15simulation(with2times worse spatial and8times worse mass resolutions)is shown with a dashed line in the middle panel.This profile does not extend to values of r/r0which are as low as for theΛCDM7.5pro-file(due to worse spatial resolution).However,for values of r/r0,where the two profiles overlap,they are indistinguish-able.This suggests that the shape is not affected by thefinite size of the simulation box and mass resolution.The peculiar upturn in the rotation curve of NGC2915is discussed in§2.2.A similar correlation can be observed in the r max−v max plane,shown in Figure4(values of r max and v max for each galaxy are given in Table1).The maximum point in a galaxy’s DM velocity profile and the corresponding radius is a nice set of physical parameters for comparison with simulations.Ideally,such a comparison would not force any pre-supposedfit to either the data or the the simu-lated profiles.Unfortunately,most of the galaxy rotation。

氮气吸脱附曲线英语

氮气吸脱附曲线英语

氮气吸脱附曲线英语Here's a sample text written in English, following the given requirements:Okay, let's talk about the nitrogen adsorption-desorption curve. It's basically a graph that shows how nitrogen molecules interact with a solid surface. You know, when you expose a material to nitrogen gas, it starts adsorbing the gas molecules on its surface. And when you change the pressure or temperature, those molecules can desorb back into the gas phase.Now, this curve is really useful in understanding the porous structure of materials. Like, if you see a steeprise at low pressures, that usually means there are lots of micropores in the material. And the plateau region at higher pressures tells you about the mesopores and macropores.But here's the cool part: the shape of the curve canalso reveal the type of interaction between the nitrogen and the material. Like, if the desorption curve doesn't follow the same path as the adsorption curve, that suggests some sort of interaction between the gas and the solid.And speaking of interactions, did you know that the nitrogen adsorption-desorption curve can also be used to calculate the surface area of a material? Yeah, it's pretty amazing. By measuring the amount of nitrogen adsorbed at different pressures, you can estimate the total surface area of the material, even down to the nanometer scale.So in a nutshell, this curve is a powerful tool for characterizing porous materials. It gives you insights into their pore structure, surface area, and even the nature of interactions between the solid and the gas. And all this from just a simple graph!。

尼康3D扫描仪ModelMaker H120产品介绍说明书

尼康3D扫描仪ModelMaker H120产品介绍说明书

Premium portable metrologyModelMakerHandheld scannersMCAxArticulated armsNIKON METROLOGY I VISION BEYOND PRECISIONUltra-fast high-definition3D scanning MODELMAKER H120More than two decades since theinception of the ModelMaker productline, the cutting-edge ModelMakerH120 firmly pushes the ever-exactingboundaries of handheld laser scanning.Incorporating blue laser technology,ultra-fast frame rate, speciallydeveloped Nikon optics and the abilityto measure the most challengingmaterials this represents the nextgeneration of portable laser scanning.The H120 makes no compromisesin addressing the market needsby efficiently deliveringthe most detailedand accurate datain a fraction of thetime of competingtechnologies.UNCOMPROMISING PERFORMANCEBy combining a frame rate of 450 Hz, a stripe width of 120 mm and a resolution of 35 μm, users benefit from high productivity and detailed measurements with a single sensor. Without relying on interpolation techniques to artificially boost data density, the ModelMaker H120guarantees fast data collection over a large area without compromising on small details – offering great flexibility in a single solution even when cycle time is critical, no matter the type of parts measured. Furthermore, the superior accuracy of the ModelMaker H120 ensures it stands far apart from similar technology, further pushing the traditionally accepted boundaries of handheld laser scanners.MEASURE THE MOST CHALLENGING MATERIALSThe 4th generation of Nikon’s patented Enhanced Sensor Performance (ESP4) provides faster-then-ever real-time dynamic adjustment of the laser intensity for every point. Users can confidently scan across parts with strong colour transitions and varying reflectivity from any direction with no loss in scanner speed and no need for prior part preparation. ModelMaker scanners also benefit from intelligent reflection control which allows users to measure very shiny or polished materials while unwanted reflections are filtered out.IMMEDIATE PRODUCTIVITYSimple system set-up, immediate boot-up and no need for scanner warm-up combined with the structural rigidity, thermal stability and absolute encoder technology of the MCAx arms allows users to switch on and start confidently collecting accurate data straightaway.EXTREMEL Y LOW NOISE DATABy combining specially-developed Nikon optics and low-speckle blue laser technology, the ModelMaker H120 achieves super low-noisemeasurements and can cleanly resolve details such as sharp edges and even surface scratches and abrasions which other scanners simply cannot.ENHANCED USER EXPERIENCEInnovative features such as thermal compensation, an integrated locking connector, contrasting full field of view projector, excellent touch probe clearance and a compact size give the user all the feedback and assurance he needs to concentrate purely on the measurement task.Integrated lockBest in class accuracyCompact and lightweightLow noise blue laserUp to 450,000 points per secondFull FOV indicatorMODELMAKER MMDxINTUITIVE SCANNING AND INSPECTION SOFTWAREScanning technology optimized for your applicationThe ModelMaker MMDx range of handheld laser scanners is ideally suited for portable 3D inspection and reverse engineering applications. With choices of scanner models for high detail, all-round scanning or high productivity, users can select the best hardware for their needs.MMDx incorporates 3rd generation Enhanced Sensor Performance (ESP3) to scan almost any sample materials and surface finishes without user interaction.The digital camera technology offers a measuring accuracy down to 7 microns and benefits from a true non-interpolated resolution of more than a thousand points per stripe, allowing freeform surfaces and features to be scanned accurately and efficiently.Featuring high frame rates and laser stripes up to 200 mm, the MMDx range provides the ultimate in scanning productivity. The scanner’s digital cameras benefit from a true (non-interpolated) resolution of over 1000 points per stripe, providing optimum resolution for scanning freeform surfaces and features efficiently.Weighing around 400 g and featuring an angled laser plane for comfort while scanning, MMDx scanners are optimized for ergonomic use. Set-up time and portability is optimized through the use of isolated thermal zones, temperature compensation and on-board processing – which means no external controller or extraneous cabling.ModelMaker scanners and MCAx arms seamlesslyinteract with Focus software for scan and tactile probe dataacquisition and inspection processing. It is a total solution that tightlyintegrates hardware and software to guarantee smooth and error-free operation.Focus software is specifically designed to easily control data flows with minimal user interaction.Users can complete handheld data acquisition and inspection jobs in Focus without compromising performance.Alternatively, through the Nikon Metrology API, the ModelMaker scanners and MCAx arms can be used directly in many 3rd party inspection and reverse engineering software applications, including PolyWorks ®, Metrolog ® and Geomagic ®.Scan rateProductivity Resolution Accuracy H120• • • • • •• • • • • • •• • • • •• • • • •MMDx50• •• •• • • •• • • •MMDx100• •• •• •• •MMDx200• •• • • •••APPLICATIONSThe combined solution of ModelMaker scanners and MCAx arms delivers high-productivity and precise non-contact and contact metrology. Used to optimise production workflow through rapid, reliable and accurate analysis of product dimensions – both freeform and geometric – it has proven to be an invaluable tool across many industries and workplaces from the shop floor to the metrology lab.Able to robustly measure almost any material and with the flexibility to inspect parts of sizes ranging from a few millimetres to several metres and more allows the solution to span many industries including Automotive, Aerospace, Power Generation and Consumer products, and well as Universities, Research Institutes and scanning service providers – especially for components such as tools and dies, body-in-white / sheet metal parts, castings, injection moulded, soft or fragile materials and additive manufactured parts.The ModelMaker handheld laser scanners paired with MCAx portable articulated co-ordinate measuring arms and Focus software allow you to reduce measurement times by rapidly diagnosing production issues in all areas of manufacture. This enables delivery of your products faster and with greater confidence by meeting the highest quality standards.Key benefits for your application• High accuracy and fast data throughput saves time and money • Optimized for hard-to-scan surfaces• Designed for use under all shop floor or field conditions• Extreme temperature stability and zero warm-up time• Quick and easy plug-and-play setup• Enhanced ergonomics for stress-free usage• Short learning curve• Seamless transition between scanning and touch-probing • Compatible with all major brands of point cloud software Uses within your process• Fast & accurate multi-sensor 3D inspection• Part-to-CAD inspection: First article inspection against CAD model • Inspection of geometric features• Gap-and-flush inspection• Reverse engineering: from concept studio clay to class A surfaces •Digitizing for additive manufacturingMCAxAccurate and portable multi-sensor measurement• Tactile probing performance from 0.023 mm and scanning system accuracy from 0.028 mm ensures the highest standard of measurement results • Available in six lengths between 2.0 m and 4.5 m to suit small to large measurement tasks• Absolute encoder technology means no referencing or warm-up period is required• Advanced carbon fiber construction for strength and thermally stability in all environments• Automatic probe recognition and repeatable probe and scanner mounting allows immediate switching between measurement tools• The ergonomic wrist features haptic feedback whilst the arm provides audio and visual notifications • Low friction handling positions for reduced user stress and fatigue• Counterbalance for effortless control infinite rotation of all principle axes for unrestricted use • Integrated lock secures the arm easily and safely • Quickly and easily attaches to a variety of stands / tripods or vacuum mount• Supports a wide variety of fixed and touch-trigger probes in many lengths and stylus configurations • MCAx++ and MCAx+ include a certified length standard for performance verification in the field • Certified performance according to ASME B89.4.22. VDI/VDE 2617-9 certification is also availableThe MCAx Manual Coordinate measuring Arm is a precise, reliable and easy-to-use portable 7-axis measuring arm. It is the perfect partner for the ModelMaker H120 and MMDx laser scanners and Focus Handheld scanning and inspection software due to its high precision, repeatability and stability. This total solution’s accuracy, capability and portability make it feel perfectly at home in the metrology lab, on the shop floor and in-the-field.The arm can be equipped with a wide range of probing systems aside from laser scanning, such as a large choice of probes for a variety of tasks including touch-trigger measurements and continuous scanning. Its flexibility makes this measurement arm the perfect solution for the widest range of measurement tasks. The MCAx range of 7-axis articulated arms is available in six different sizes and in three accuracy levels giving users the ability to specify thebest system for their needs.Infinite joint rotationRotating gripsAbsolute encodersCarbon constructionCounterbalanceIntegrated lockMagnetic mountSPECIFICATIONSComplies with 21 CFR 1040.10 and 1040.11, Laser Notice No. 50 dated June 24, 20071 Typical values are 30% better than published values.2 L aser scanner Accuracy is determined by scanning a plane from various directions, each time using the entire scanner field of view. The result is the maximum 1σ deviation of the scan data to fitted plane features.3T he Scanning performance test indicates the performance of the laser scanner combined with a MCAx arm. The test is performed by scanning a highly accurate reference plate in 5 different orientations of the articulated arm and laser scanner. The 5 resulting point clouds are merged together and a best-fit plane is constructed through this combined point cloud. For each of the points, the deviation distance to the best-fit plane is calculated. The result of the test is the 2σ value of all of the deviations.4T he Point repeatability test (or SPAT) is the reference test to determine measurement arm repeatability with a ball probe. The probe is placed in a conical socket and points are measured from multiple approach directions and is tested different zones of the arm measurement volume. The result is the maximum of the X, Y or Z range divided by two.5T he Volumetric accuracy test most accurately represents the reasonable expectations for probing performance in practical measuring applications since it involves measuring a certified length standard many times in several locations and orientations and compares the resulting measurements to the actual length. It is the most appropriate test for determining machine accuracy and repeatability. The result is the maximum deviation of the measuring distance less the theoretical length.Probing and scanning specifications are achieved under stable environmental conditions with the MCAx arm mounted on a base plate or magnetic base. A 15 mm diameter, 50 mm long, steel ball probe connected to both probeports is used for the probing performance tests. Probing specifications are based on a subset of ASME B89.4.22:2004. Probing certification to VDI/VDE 2617-9 is also available.Working temperature Storage temperature 0 – 50˚ C -30 – 70˚ CRelative humidity Operational elevation10 – 90% non-condensing 0 – 2000 m Universal worldwide voltage CE Compliance 110 – 240 V AC (50 – 60 Hz)YesH -120_E N _0118– C o p y r i g h t N i k o n M e t r o l o g y N V 2018. A l l r i g h t s r e s e r v e d . T h e m a t e r i a l s p r e s e n t e d h e r e a r e s u m m a r y i n n a t u r e , s u b j e c t t o c h a n g e a n d i n t e n d e d f o r g e n e r a l i n f o r m a t i o n o n l y.。

FVCOM使用手册(完整版)

FVCOM使用手册(完整版)

An Unstructured Grid, Finite-Volume Coastal Ocean Model(无规则网格的有限体积海岸海洋模型)FVCOM User Manual(FVCOM用户手册)FVCOM软件用户许可协议 (3)第一章序言 (4)第二章:模型公式 (6)2.1 直角坐标系下的原始方程 (7)2.2 -坐标下的控制方程 (12)2.3 二维(垂直积分)方程 (13)2.4 湍流闭合模型 (15)2.4.1 水平扩散系数 (15)2.4.2 垂直旋转粘性和热扩散系数 (16)2.5 球面坐标系下的原始方程 (24)第三章有限体积离散法 (27)3.1 不规则三角网格的设计 (27)3.2 笛卡尔坐标下的离散方法 (29)3.2.1 二维外部模式 (29)3.2.2 三维内模式 (37)3.3 外部与内部模式的输运一致性 (44)3.4 干/湿处理方法 (46)3.4.1 标准 (48)3.4.2 Isplit的上限 (52)3.5 球坐标系下的有限体积离散方法 (57)3.6 岸边界条件的微元处理 (63)第四章:外部强迫 (66)4.1 风应力、热通量和降水/蒸发 (66)4.2 潮汐强迫 (67)4.3 增加海岸或江河流量的方法 (69)4.3.1 TCE方法 (69)4.3.2 MCE方法 (72)4.4 水平分辨率和时间步长的规范 (74)4.5 通过底部输入地下水 (77)4.5.1 简单盐平衡地下水通量形式 (77)4.5.2 地下水输入的完全格式 (78)第五章:开边界处理 (79)5.1 开边界处理的初始设定 (79)5.2 普遍辐射开边界条件 (82)5.3 新的有限体积开边界条件模块 (87)第六章:数据同化方法 (97)6.1 推导方法 (100)6.2 OI方法 (102)6.3 Kalman筛选 (104)6.3.1减小序列Kalman筛选(RRKF) (106)6.3.2 集合Kalman筛选(EnKF) (109)6.3.3 集合平方根Kalman过滤(EnSRF) (111)6.3.4. 集合变换 Kalman筛选 (ETKF) (113)6.3.5确认实验 (114)第七章:FVCOM沉积模块 (120)7.1 控制方程 (121)7.2 简单测试情况 (122)第八章:FVCOM生物模块 (123)8.1灵活生物模块(FBM) (124)8.1.1 FBM流程图 (124)8.1.2 FBM中的方程和函数 (126)8.2 提前选择生物模块 (157)8.2.1 养分-浮游植物-浮游动物(NPZ模型) (158)8.2.2 磷限制低养分层食物网模型 (160)8.2.3. The Multi-Species NPZD Model (168)8.2.3 多物种NPZD模型 (168)8.2.4 水质量模型 (171)第九章:示踪-追踪模型 (174)第十章:三维拉格朗日粒子追踪 (175)第十二章:代码平行 (193)12.1 区域分解 (194)12.2 区域设置 (195)12.3 数据交换 (196)12.4数据收集 (197)12.5 执行 (198)第十三章:模型代码描述和总说明 (199)13.1 在使用FVCOM前的用户应知 (199)13.3 数值稳定的标准 (206)13.4子程序和函数描述 (207)第14章模式安装,编译和运行 (231)14.1 获得FVCOM (232)14.2a 编译METIS库 (233)14.2b 编译FVCOM (233)14.3a 运行FVCOM(连续) (238)14.3b 运行FVCOM(平行) (239)第十五章:模型设置 (240)15.1 FVCOM运行时间控制变量文件casename_run.dat (240)15.2 FVCOM输入文件 (253)15.3特殊设置的必需输入文件 (256)15.4 原始输入文件的输入文件格式 (257)15.5 建立和使用FVCOM模块 (268)第十六章:FVCOM测试例子 (292)第十七章:不规则三角形网格产生 (319)17.1数据准备 (320)17.2 网格产生 (324)感谢 (347)参考文献 (348)FVCOM Software Users’ License AgreementFVCOM软件用户许可协议All users should read this agreement carefully. A user, who receives any version of the source code of FVCOM, must accept all the terms and conditions of this agreement and also agree that this agreement is like any written negotiated agreement signed by you. You may be required to have another written agreement directly with Dr. Changsheng Chen at SMAST/UMASS-D and Dr. Robert C. Beardsley at WHOI 所有用户须仔细阅读此协议。

Autodesk VR AR技术在工具路径和探头路径定义和可视化中的应用说明书

Autodesk VR AR技术在工具路径和探头路径定义和可视化中的应用说明书

MFG124360Exploring Toolpath and Probe Path Definition and Visualization in VR/ARZhihao CuiAutodeskDescriptionVisualizing and defining 3D models on a 2D screen has always been a challenge for CAD and CAM users. Tool paths and probe paths add other levels of complexity to take into consideration, as the user cannot fully appreciate the problem on a 2D viewer. Imagine yourself trying to define a tool axis on a complex shape—it’s very hard to take every single aspect of the shape into account, except by guessing, calculating, and retrying repeatedly. With augmented reality (AR) and virtual reality (VR) technologies, the user gains the ability to inspect and define accurate 3D transformations (position and rotation) for machine tools in a much more natural way. We will demonstrate one potential workflow to address this during the class, which includes how to export relevant models from PowerMill software or PowerInspect projects; how to reconstruct, edit, and optimize models in PowerShape software and 3ds Max software; and eventually how to add simple model interactions and deploy them in AR/VR environments with game engines like Stingray or Unity.SpeakerZhihao is a Software Engineer in Advanced Consulting team within Autodesk. His focus for AR and VR technologies is in manufacturing industry and he wishes to continuously learn and contribute to it.Data PreparationThe first step of the journey to AR or VR is generating the content to be visualized. Toolpath and probe path need to be put into certain context to be meaningful, which could be models for the parts, tools, machines or even the entire factory.PowerMill ExportFigure 1 Typical PowerMill ProjectPartsExporting models of the part is relatively simple.1. Choose the part from the Explorer -> Models -> Right click on part name -> ExportModel…2. Follow the Export Model dialog to choose the name with DMT format1.1DGK file is also supported if additional CAD modification is needed later. See Convert PowerMill part mesh on page 7Figure 2 PowerMill – Export ModelToolTool in PowerMill consists three parts – Tip, Shank and Holder.To export the geometry of the tool, type in the macro commands shown in Figure 3, which would generate STL files2 contains the corresponding parts. Three lines of commands3 are used instead of exporting three in one file (See Figure 11), or one single mesh would be created instead of three which will make the coloring of the tool difficult.EDIT TOOL ; EXPORT_STL TIP "powermill_tool_tip.stl"EDIT TOOL ; EXPORT_STL SHANK "powermill_tool_shank.stl"EDIT TOOL ; EXPORT_STL HOLDER "powermill_tool_holder.stl"Figure 3 PowerMill Macro - Export ToolToolpathToolpath is the key part of the information generated by a CAM software. They are created based on the model of the part and various shapes of the tool for different stages (e.g. roughing, polishing, etc.). Toolpaths are assumed to be fully defined for visualization purposes in this class, and other classes might be useful around toolpath programming, listed on page 15. Since there doesn’t exist a workflow to directly stream data into AR/VR environment, a custom post-processor4 is used to extract minimal information needed to describe a toolpath, i.e. tool tip position, normal direction and feed rate (its format is described in Figure 17).The process is the same way as an NC program being generated for a real machine to operate. Firstly, create an NC program with the given post-processor shown in Figure 4. Then grab and drop the toolpath onto the NC program and write it out to a text file shown in Figure 5.2DDX file format can also be exported if geometry editing is needed later in PowerShape3 The macro is also available in addition al material PowerMill\ExportToolMesh.mac4 The file is in additional material PowerMill\simplepost_free.pmoptzFigure 4 PowerMill Create NC ProgramFigure 5 PowerMill Insert NC ProgramPowerInspect ExportFigure 6 Typical PowerInspect OMV ProjectCADCAD files can be found in the CAD tab of the left navigation panel. The model can be re-processed into a generic mesh format for visualization using PowerShape, which is discussed in Section Convert PowerMill part mesh on page 7.Figure 7 Find CAD file path in PowerInspectProbeDefault probe heads are installed at the following location:C:\Program Files\Autodesk\PowerInspect 2018\file\ProbeDatabaseCatalogueProbes shown in PowerInspect are defined in Catalogue.xml file and their corresponding mesh files are in probeheads folder. These files will be used to assemble the probe mentioned in Section Model PowerInspect probe on page 9.Probe toolAlthough probe tool is defined in PowerInspect, they cannot be exported as CAD geometries to be reused later. In Model PowerInspect probe section on page 9, steps to re-create the probe tool will be introduced in detail based on the stylus definition.Probe pathLike toolpath in PowerMill, probe path can be exported using post processor5 to a generic MSR file format, which contains information of nominal and actual probing points, measuring tolerance, etc.This can be achieved from Run tab -> NC Program, which is shown in Figure 8.Figure 8 Export Probe Path from PowerInspect5 The file is in additional materialPowerInspect\MSR_ResultsGenerator_1.022.pmoptzModelling using PowerShapeConvert PowerMill part meshDMT or DGK files can be converted to mesh in PowerShape to FBX format, which is a more widely adopted format.DMT file contains mesh definition, which can be exported again from PowerShape after color change and mesh decimation if needed (discussed in Section Exporting Mesh in PowerShape on page 10).Figure 9 PowerShape reduce meshDGK file exported from PowerMill / PowerInspect is still parametric CAD model not mesh, which means further editing on the shape is made possible. Theoretically, the shape of the model won’t be changed since the toolpath is calculated based on the original version, but further trimming operations could be carried here to keep minimal model to be rendered on the final device. For example, not all twelve blades of the impeller may be needed to visualize the toolpath defined on one single surface. It’s feasible to remove ten out of the twelve blades and still can verify what’s going on with the toolpath defined. After editing the model, PowerShape can convert the remaining to mesh and export to FBX format as shown below.Figure 10 Export FBX from PowerShapeModel PowerMill toolImport three parts of the tool’s STL files into PowerShape, and change the color of individual meshes to match PowerMill’s color scheme for easier recognition.Figure 11 PowerShape Model vs PowerMill assembly viewBefore exporting, move the assembled tool such that the origin is at the tool tip and oriented z-axis upwards, which saves unnecessary positional changes during AR/VR setup. Then follow Figure 10 to export FBX file from PowerShape to be used in later stages.Model PowerInspect probeTake the example Probe OMP400. OMP400.mtd file6 contains where the mesh of individual components of the probe head are located and their RGB color. For most of the probe heads, DMT mesh files will be located in its subfolder. They can be dragged and dropped into PowerShape in one go to form the correct shape, but all in the same color (left in Figure 14). To achieve similar looking in PowerInspect, it’s better to follow the definition fi le, and import each individual model and color it according to the rgb value one by one (right in Figure 14).<!-- Head --><machine_part NAME="head"><model_list><dmt_file><!-- Comment !--><path FILE="probeheads/OMP400/body.dmt"/><rgb R="192"G="192"B="192"/></dmt_file>Figure 12 Example probe definition MTD fileFigure 13 PowerShape apply custom colorFigure 14 Before and after coloring probe headFor the actual probe stylus, it’s been defined in ProbePartCatalogue.xml file. For theTP20x20x2 probe used in the example, TP20 probe body, TP20_STD module and M2_20x2_SS stylus are used. Construct them one by one in the order of probe body, module and stylus, and each of them contains the definition like the below, which is the TP20 probe body.6C:\Program Files\Autodesk\PowerInspect 2018\file\ProbeDatabaseCatalogue<ProbeBody name="TP20"from_mounting="m8"price="15.25"docking_height="0"to_mounting="AutoMagnetic"length="17.5"><Manufacturer>Renishaw</Manufacturer><Geometry><Cylinder height="14.5"diameter="13.2"offset="0"reference_length="14.5" material="Aluminium"color="#C8C8C8"/><Cylinder height="3.0"diameter="13.2"offset="0"reference_length="3.0" material="Stainless"color="#FAFAFA"/></Geometry></ProbeBody>Figure 15 Example TP20 probe body definitionAlmost all geometries needed are cylinder, cone and sphere to model a probing tool. Start with the first item in the Geometry section, and use the parameters shown in the definition to model each of the geometries with solid in PowerShape and then convert to mesh. To make the result look as close as it shows in PowerInspect, color parameter can also be utilized (Google “color #xxx” to convert the hex color).Figure 16 Model PowerInspect ToolU nlike PowerMill tool, PowerInspect probe’s model origin should be set to the probe center instead of tip, which is defined in the MSR file. But the orientation should still be tuned to be z-axis facing upwards.DiscussionsExporting Mesh in PowerShapeIn PowerShape, there are different ways that a mesh can be generated and exported. Take the impellor used in PowerMill project as an example, the end mesh polycount is 786,528 if it’s been converted from surfaces to solid and then mesh with a tolerance set to 0.01. However, if the model was converted straight from surface to mesh, the polycount is 554,630, where the 30% reduce makes a big impact on the performance of the final AR/VR visualization.Modifying the tolerance could be another choice. For visualization purposes, the visual quality will be the most impactable factor of choosing the tolerance value. If choosing the value is set too high, it may introduce undesired effect that the simulated tool is clipped into the model in certain position. However, setting the tolerance too small will quickly result in a ridiculous big mesh, which will dramatically slow down the end visualization.Choosing the balance of the tolerance here mainly depends on what kind of end devices will the visualization be running on. If it will be a well-equipped desktop PC running VR, going towards a large mesh won’t ne cessarily be a problem. On the other hand, if a mobile phone is chosen for AR, a low polycount mesh will be a better solution, or it can be completely ignored as a placeholder, which is discussed in Section On-machine simulation on page 12.Reading dataSame set of model and paths data can be used in multiple ways on different devices. The easiest way to achieve this is through game engines like Stingray or Unity 3D, which has built-in support for rendering in VR environment like HTC Vive and mobile VR, and AR environment like HoloLens and mobile AR.Most of the setup in the game engine will be the same for varies platform, like models and paths to be displayed. Small proportion will need to be implemented differently for each platform due to different user interaction availability. For example, for AR when using HoloLens, the user will mainly control the application with voice and gesture commands, while on the mobile phones, it will make more sense to offer on-screen controls.For part and tool models, FBX files can be directly imported into the game engines without problem. Unit of the model could be a problem here, where export from PowerShape is usually in millimeter but units in game engines are normally in meters. Unit change in this case could result in a thousand times bigger, which may cause the user seeing nothing when running the application.For toolpath data, three sets of toolpath information are exported from PowerMill with the given post-processor, i.e. tool tip position, tool normal vector and its feed rate. They can be read line by line, and its positions can be used to create toolpath lines. And together with the normal vector and feed rates, an animation of the tool head can be created.Position(x,y,z) Normal(i,j,k) Feed rate33.152,177.726,52.0,0.713,-0.208,0.67,3000.0Figure 17 Example toolpath output from PowerMillFor probe path data, similar concept could be applied with an additional piece of information7–actual measured point, which means not only the nominal probe path can be simulated ahead of time, but also the actual measured result could be visualized with the same setup.7 See page 14 for MSR file specification.STARTG330 N0 A0.0 B0.0 C0.0 X0.0 Y0.0 Z0.0 I0 R0G800 N1 X0 Y0 Z25.0I0 J0 K1 O0 U0.1 L-0.1G801 N1 X0.727 Y0.209 Z27.489 R2.5ENDFigure 18 Example probe path output from PowerInspectUse casesOn-machine simulationWhen running a new NC program with a machine tool, it’s common to see the machine operator tuning down the feed rate and carefully looking through the glass to see what is happening inside the box. After several levels of collision checking in CAM software and machine code simulator, why would they still not have enough confidence to run the program?Figure 19 Toolpath simulation with AR by Hans Kellner @ AutodeskOne potential solution to this problem is using AR on the machine. Since how the fixture is used nowadays is still fairly a manual job constrained by operator’s experience, variations of fixtures make it a very hard process to verify ahead of machining process. Before hitting the start button for the NC program, the operator could start the AR simulation on the machine bed, with fixtures and part held in place. It will become an intuitive task for the operator to check for collisions between part of the virtual tool and the real part and fixtures. Furthermore, a three-second in advance virtual simulation of the tool head can be shown during machining process to significantly increase the confidence and therefore leave the machine always running at full speed, which ultimately increases the process efficiency.Toolpath programming assistanceProgramming a toolpath within a CAM software can sometimes be a long iterative try and error process since the user always imagines how the tool will move with the input parameters. Especially with multi-axis ability, the user will often be asked to provide not only the basic parameters like step over values but also coordinate or direction in 3D for the calculation to start. Determining these 3D values on a screen becomes increasingly difficult when othersurfaces surround the places needed to be machined. Although there are various ways to let the user to navigate to those positions through hiding and sectioning, workarounds are always not ideal and time-consuming. As shown in Figure 20, there’s no easy and intuitive way to analyze the clearance around the tool within a tight space, which is one of the several places to be considering.Figure 20 Different angles of PowerMill simulation for a 5-axis toolpath in a tight spaceTaking the example toolpath in PowerMill, a user will need to recalculate the toolpath after each modification of the tool axis point, to balance between getting enough clearance8and achieving better machining result makes the user and verify the result is getting better or worth. However, this workflow can be changed entirely if the user can intuitively determine the position in VR. The tool can be attached to the surface and freely moved by hand in 3D, which would help to determine the position in one go.Post probing verificationProbing is a common process to follow a milling process on a machine tool, making sure the result of the manufacturing is within desired tolerance. Generating an examination report in PowerInspect is one of the various ways to control the quality. However, what often happens is that if an out of tolerance position is detected, the quality engineer will go between the PC screen and the actual part to determine what is the best treatment process depending on different kind of physical appearance.8 Distance between the tool and the partFigure 21 Overlay probing result on to a physical partOverlaying probing result with AR could dramatically increase the efficiency by avoiding this coming back and forth. Same color coded probed point can be positioned exactly at the place of occurrence, so that the surrounding area can be considered separately. The same technique could also be applied to scanning result, as shown in Figure 22.Figure 22 Overlaying scanning result on HoloLens by Thomas Gale @ AutodeskAppendixReference Autodesk University classesPowerMillMFG12196-L: PowerMILL Hands on - Multi Axis Machining by GORDON MAXWELL MP21049: How to Achieve Brilliant Surface Finishes for CNC Machining by JEFF JAJE MSR File format9G330 Orientation of the probeG800 Nominal valuesG801 Measured valuesN Item numberA Rotation about the X axisB Rotation about the Y axisC Rotation about the Z axisXYZ Translations along the X, Y and Z axes (these are always zero)U Upper toleranceL Lower toleranceO OffsetI and R (in G330) just reader valuesR (in G801) probe radius9 Credit to Stefano Damiano @ Autodesk。

抽运-检测型非线性磁光旋转铷原子磁力仪的研究

抽运-检测型非线性磁光旋转铷原子磁力仪的研究

抽运-检测型非线性磁光旋转铷原子磁力仪的研究缪培贤;杨世宇;王剑祥;廉吉庆;涂建辉;杨炜;崔敬忠【摘要】报道了一种抽运-检测型的非线性磁光旋转铷原子磁力仪.其原理是线偏振光通过处于外磁场环境中被极化的原子介质后,由于原子对线偏振光中左、右圆偏成分不同的吸收和色散,导致光的偏振方向会产生与磁场相关的转动.分析了该磁力仪的工作原理,并测试了它对不同磁场大小的响应.测试结果表明,磁力仪测量范围为100—100000 nT,极限灵敏度为0.2 pT/Hz1/2,磁场分辨率为0.1 pT.进一步研究了不同磁场下原子系综极化态的横向弛豫时间,讨论了原子磁力仪高磁场采样率的获得方法.本文的原子磁力仪在5000—100000 nT的磁场测量范围内磁场采样率可实现1—1000 Hz范围内可调,能够测量低频的微弱交变磁场.本文的研究内容为大磁场测量范围、高灵敏度、高磁场采样率的原子磁力仪研制提供了重要参考.%We report a rubidium atomic magnetometer based on pump-probe nonlinear magneto-optical rotation. The rubid-ium vapor cell is placed in a five-layer magnetic shield with inner coils that can generate uniform magnetic fields along the direction of pump beam, and the cell is also placed in the center of a Helmholtz coil that can generate an oscillating magnetic field perpendicular to the direction of pump beam. The atoms are optically pumped by circularly polarized pump beam along a constant magnetic field in a period of time, then the pump beam is turned off and a π/2 pulse of oscillating magnetic fiel d for 87Rb atoms is applied. After the above process, the individual atomic magnetic moments become phase coherent, resulting in a transverse magnetization vector precessing at the Larmor frequency in the mag-netic field. The linearly polarized probingbeam is perpendicular to the direction of magnetic field, and can be seen as a superposition of the left and right circularly polarized light. Because of the different absorptions and dispersions of the left and right circularly polarized light by rubidium atoms, the polarization direction of probing beam rotates when probing beam passes through rubidium vapor cell. The rotation of the polarization is subsequently converted into an electric signal through a polarizing beam splitter. Finally, the decay signal related to the transverse magnetization vector is measured. The Larmor frequency proportional to magnetic field is obtained by the Fourier transform of the decay signal. The value of magnetic field is calculated from the formula: B =(2π/γ)f , where γ and f are the gyromagnetic ratio and Larmor frequency, respectively. In order to measure the magnetic field in a wide range, the tracking lock mode is proposed and tested. The atomic magnetometer can track the magnetic field jump of 1000 nT or 10000 nT, indicating that the atomic magnetometer has strong locking ability and can be easily locked after start-up. The main performances in different magnetic fields are tested. The results show that the measurement range of the atomic magnetometer is from 100 nT to 100000 nT, the extreme sensitivity is 0.2 pT/Hz1/2, and the magnetic field resolution is 0.1 pT. The transverse relaxation times of the transverse magnetization vector in different magnetic fields are obtained, and the relaxation time decreases with the increase of the magnetic field. When the measurement range is from 5000 nT to 100000 nT, the magnetic field sampling rate of the atomic magnetometer can be adjusted in a range from 1 Hz to 1000 Hz. Theatomic magnetometer in high sampling rate can measure weak alternating magnetic field at low frequency. This paper provides an important reference for developing the atomic magnetometer with large measurement range, high sensitivity and high sampling rate.【期刊名称】《物理学报》【年(卷),期】2017(066)016【总页数】11页(P47-57)【关键词】原子磁力仪;非线性磁光旋转;灵敏度;磁场采样率【作者】缪培贤;杨世宇;王剑祥;廉吉庆;涂建辉;杨炜;崔敬忠【作者单位】兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州 730000;兰州空间技术物理研究所, 真空技术与物理重点实验室, 兰州730000【正文语种】中文报道了一种抽运-检测型的非线性磁光旋转铷原子磁力仪.其原理是线偏振光通过处于外磁场环境中被极化的原子介质后,由于原子对线偏振光中左、右圆偏成分不同的吸收和色散,导致光的偏振方向会产生与磁场相关的转动.分析了该磁力仪的工作原理,并测试了它对不同磁场大小的响应.测试结果表明,磁力仪测量范围为100—100000 nT,极限灵敏度为0.2 pT/Hz1/2,磁场分辨率为0.1 pT.进一步研究了不同磁场下原子系综极化态的横向弛豫时间,讨论了原子磁力仪高磁场采样率的获得方法.本文的原子磁力仪在5000—100000 nT的磁场测量范围内磁场采样率可实现1—1000 Hz范围内可调,能够测量低频的微弱交变磁场.本文的研究内容为大磁场测量范围、高灵敏度、高磁场采样率的原子磁力仪研制提供了重要参考.高灵敏度的原子磁力仪在生物医学[1,2]、惯性导航[3,4]、军事磁异反潜[5]、基础物理研究等[6−9]领域具有重要的应用.目前国际上出现了Mz和Mx模式的光泵磁力仪、相干布居囚禁磁力仪、非线性磁光旋转(nonlinear magneto-optical rotation,NMOR)磁力仪、无自旋交换弛豫(spin-exchange relaxation free,SERF)磁力仪等多种原子磁力仪[10],其中SERF磁力仪灵敏度已达到fT/Hz1/2量级[11−13].近年来,国内有多家单位开展了原子磁力仪的研究.例如浙江大学研制了铷光泵磁力仪,零磁场附近灵敏度达到0.5 pT/Hz1/2[14];北京大学详细讨论了铯光泵磁力仪的参数优化问题,得到最优的灵敏度为2.5 pT/Hz1/2[15];国防科学技术大学研制了NMOR铷原子磁力仪,测量范围为±60 nT,灵敏度达到1 pT/Hz1/2[16],后来经过进一步优化实验条件,灵敏度达到0.2 pT/Hz1/2[17].总体而言,国内原子磁力仪的研制还处于起步阶段,在灵敏度、测量范围、磁场采样率等指标上还有很大的提升空间[18].本文系统地研究了抽运-检测型的NMOR铷原子磁力仪,测试结果表明,磁力仪测量范围为100—100000 nT,极限灵敏度为0.2 pT/Hz1/2,磁场分辨率为0.1 pT,磁场采样率最高可达1000 Hz.研究的NMOR铷原子磁力仪用两束激光完成外磁场中原子系综极化态的制备与探测,圆偏振抽运光与外磁场平行,线偏振探测光与外磁场垂直.铷原子磁力仪采用87Rb原子D1线跃迁制备极化态原子介质,即基态52S1/2到第一激发态的52P1/2的跃迁,对应波长为795 nm.基态52S1/2的两个精细能级分别是52S1/2(Mj=−1/2)和52S1/2(Mj=+1/2),795 nm的左旋圆偏振光(σ+光子)可被处于52S1/2(Mj= −1/2)基态的87Rb原子吸收,使得87Rb原子跃迁到52P1/2(Mj=+1/2)激发态上,激发态87Rb原子通过辐射光子后跃迁到52S1/2(Mj=−1/2)或52S1/2(Mj=+1/2)基态上,左旋圆偏振光持续作用将使铷泡内绝大部分87Rb原子最终处于52S1/2(Mj=+1/2)基态上.同理,右旋圆偏振光(σ−光子)持续作用将使铷泡内绝大部分87Rb原子最终处于52S1/2(Mj=−1/2)态上.这样,圆偏振的抽运光完成了原子系综极化态的制备.这里引入二能级磁共振的经典物理图像来解释NMOR铷原子磁力仪的工作原理[19].经过抽运光作用后,极化态的87Rb原子磁矩与外磁场B近似平行或反平行.在与外磁场垂直的平面内施加角频率ω约等于拉莫尔进动频率ω0的激励磁场B′[19],原子磁矩将在实验室坐标系中做复杂的运动,而在以角频率ω旋转的转动坐标系中,原子磁矩绕B′做进动.由于铷泡内原子间频繁的碰撞,在激励磁场的作用下使大部分铷原子磁矩绕外磁场进动的相位角趋于一致,原子系综呈现出绕外磁场进动的宏观磁化强度[20].原子磁矩在旋转坐标系中进动π角度时,相当于在外磁场B量子化轴方向上原子发生了磁共振跃迁.如果激励磁场持续作用,87Rb原子将在两个基态能级间来回跃迁.本文NMOR铷原子磁力仪要求原子磁矩在旋转坐标系中进动π/2角度,即原子系综宏观磁化强度进动到与外磁场B垂直的平面内,然后关闭激励磁场.线偏振光可以看作是左、右圆偏振光的矢量叠加,当线偏振的探测光穿过铷泡时,由于原子对线偏振光中左、右圆偏成分不同的吸收和色散,导致线偏振光的偏振方向会随着原子磁矩绕外磁场的拉莫尔进动而相对原来偏振方向做摆动,用差分探测方式探测偏振光偏振方向的摆动即可获得原子磁矩拉莫尔进动自由弛豫信号,并由此信号傅里叶变换出拉莫尔进动频率.由外磁场B与拉莫尔进动频率f的依赖关系可获得外磁场大小[18]:其中γ是旋磁比.对于87Rb原子,γ/2π的值为6.99583 Hz/nT[18].NMOR铷原子磁力仪要求探测光不能过于破坏原子系综的极化态,显然探测光的频率不能等于87Rb原子的D1线跃迁频率.我们在实验中设定探测光频率相对于87Rb原子的D1线跃迁频率红失谐4 GHz.研制的NMOR铷原子磁力仪如图1所示.铷泡为Φ25 mm×50 mm的圆柱型气室,气室中充有100 Torr的氮气缓冲气体,采用交流无磁加热使铷泡工作在100◦C.待测外磁场B方向与抽运光方向平行,与探测光方向垂直.实验时抽运激光被扩束为10 mm×30 mm的长方形光斑,光强为20µW/mm2;探测光为直径2mm的圆斑,进入铷泡前光功率为100µW.原子磁力仪具体工作过程是:795 nm抽运激光经过声光调制器AOM和1/4玻片形成圆偏振光,扩束后作用在铷泡上,将87Rb原子磁矩抽运在与外磁场平行的方向上;抽运激光作用一段时间后关闭,用信号源给亥姆霍兹线圈输入特定时长的正弦交变信号以产生原理部分描述的激励磁场,驱动87Rb原子磁矩在与外磁场垂直的平面内绕外磁场B做拉莫尔进动;红失谐的探测激光经过偏振片,成为线偏振光穿过铷泡,用偏振分光棱镜(PBS)、光电探测器、差分放大电路、美国NI公司的PCI-5922数据采集卡和计算机中编写的Labview程序实现铷原子拉莫尔进动信号的提取及处理,得到外磁场大小.计算机可设定数字信号处理(DSP)模块的时序组合,实现磁场采样率的设定.DSP给声光调制器AOM、信号源和PCI-5922数据采集卡输入电平触发信号,分别控制作用于铷泡的抽运激光开或关、正弦交变磁场开或关以及PCI-5922数据采集卡的采集触发.图1中铷泡、铷泡加热模块、亥姆霍兹线圈被置于五层坡莫合金的磁屏蔽筒内,磁屏蔽筒内含有可产生精密待测磁场的线圈.本文系统地研究了NMOR铷原子磁力仪的测量范围、灵敏度、分辨率、磁场采样率这些性能指标.在具体介绍这些内容之前,有必要先描述原子磁力仪的时序控制过程及跟踪式锁频过程.首先介绍原子磁力仪时序控制过程.图2显示了NMOR铷原子磁力仪在关闭抽运光后不同时长激励磁场的作用效果,外磁场环境为10000 nT.在原理部分描述到,如果抽运光作用结束后激励磁场持续作用,87Rb原子将在两个基态能级间来回跃迁.图2(a)激励磁场作用10 ms,反映了该物理过程.图2(a)中插图显示了0.5 Ms时间内的测试结果,一个包络终止代表着87Rb原子在外磁场量子化轴方向上两个基态能级间的一次跃迁.将激励磁场作用时间设定为0.1 Ms,即原子系综的宏观磁化强度进动到与外磁场垂直的平面内,测试结果如图2(b)所示,由自由弛豫过程中的正弦信号可傅里叶变换出拉莫尔进动频率.图3(a)显示了NMOR铷原子磁力仪工作时的时序示意图;图3(b)显示在10000 nT磁场环境下获得的实测数据,原子磁力仪的工作周期T=10ms,抽运激光作用时长t1=3 ms,激励磁场作用时长t2=0.1Ms,该时序磁场采样率为100 Hz;图3(c)是图3(b)中的部分曲线的放大.其次介绍原子磁力仪跟踪式锁频过程,该过程在Labview程序中完成.Labview程序在每一个原子磁力仪工作周期内能够获得拉莫尔进动频率和外磁场数值,将前一个工作周期中获得的拉莫尔进动频率设定为下一个工作周期中信号源的输出频率,即实现了跟踪式锁频.本文描述的原子磁力仪跟踪式锁频方法与Mz光泵磁力仪不同,即使激励磁场振荡频率偏离拉莫尔进动频率很远,只要特定时长激励磁场的作用能够使原子系综横向磁化强度矢量不为零,本文描述的原子磁力仪就能够实现跟踪式锁频.为了验证跟踪式锁频能力,设计这样的实验:设定原子磁力仪工作时序为T=100Ms,t1=30ms,t2=0.1Ms.设定激励磁场振荡频率为70 kHz,对应约10000 nT的测量磁场.保持激励磁场振荡频率不改变,改变线圈电流,使测量磁场从5000nT增加至15000 nT.图4(a)显示激励磁场关闭后磁力仪获得的自由弛豫正弦信号最大振幅随着扫描磁场的变化,可以看出在10000 nT附近自由弛豫正弦信号振幅最大.从原理上讲,只要横向磁化矢量不为零,铷泡中的铷原子就能够对线偏振光中左、右圆偏成分实现吸收和色散,通过差分探测获得与磁场相关的自由弛豫正弦振荡信号.横向磁化矢量越大,会使自由弛豫正弦振荡信号的振幅越大.在工作原理部分我们重点描述了激励磁场振荡角频率ω约等于拉莫尔进动角频率ω0的情况,实际上当ω与ω0相差较大时,在转动坐标系中原子磁矩会感受一有效磁场(有效磁场的描述详见参考文献[19])的作用,且在转动坐标系中磁矩进动角频率ω1为[19]可以分析,设定ω0=ω时特定时长的激励磁场作用满足π/2的脉冲效果,使横向磁化矢量最大;而后因外界磁场改变导致ω0与ω相差较大时,在特定时长内激励磁场的作用效果ω1t2可能会出现3π/2+δ,5π/2+δ′等脉冲效果,其中δ或δ′的绝对值小于等于π/2,在转动坐标系中该脉冲效果使原子磁矩在与外磁场垂直平面内的投影矢量的模达到最大值,即横向磁化矢量达到极大值,因此图4(a)中在10000 nT两侧出现若干峰值也不难理解.图4(b)显示在上述扫描磁场过程中磁力仪输出的磁场值,在自由弛豫正弦信号振幅最小时易出现与外磁场无关的数据,图4(b)中若干跳点输出磁场值用(1)式换算成频率,发现该频率正好等于铷泡交流无磁加热的输出频率.图4(b)的实验结果表明,如果该原子磁力仪在跟踪式锁频模式下工作,在很宽的磁场范围内磁力仪能够实现瞬时锁定.设定磁场线圈电流使磁屏蔽筒内磁场在10000 nT 和9000 nT,或者50000 nT和40000 nT之间来回跃变,采用跟踪式锁频模式,实验结果如图4(c)和图4(d)所示,表明该原子磁力仪对1000 nT或10000 nT的跃变磁场能够实现瞬时锁定,分别对应着7 kHz或70 kHz的频率跃变.上述实验结果表明本文描述的原子磁力仪跟踪式闭环锁定可行,而且具有很强的闭环锁定能力.接下来详细介绍NMOR铷原子磁力仪的各项性能指标.1)磁场测量范围本文的NMOR铷原子磁力仪用精密电流源给磁屏蔽筒中的磁场线圈通入逐渐增加的电流I来检验磁场测量范围,采用跟踪式锁频模式测量外磁场B的大小,测试结果如图5所示.原子磁力仪可响应100—100000 nT范围内的磁场.图5中数据线性拟合结果为从表达式(3)可知,当线圈电流I为零时,磁屏蔽筒内有约27 nT的剩余磁场.2)灵敏度和分辨率本文采用磁场噪声功率谱密度(@1 Hz)来表征原子磁力仪的灵敏度.值得注意的是,目前一些文献采用功率谱或者均方根幅度谱来表征原子磁力仪的灵敏度,从物理意义上来说是不准确的.功率谱密度使测量独立于信号持续时间和采样数量,通过功率谱密度测量可检测信号的本底噪声.若采用功率谱或均方根幅度谱,我们在实验中发现随着采样时间的延长会得到更优的灵敏度指标,显然用于表征原子磁力仪的灵敏度指标不合理.首先分析500 nT外磁场环境下如何获得磁力仪的灵敏度指标.图6(a)显示了截取的自由弛豫正弦信号,代表经过铷泡的线偏振探测光偏振方向的摆动.图6(b)是图6(a)中数据的快速傅里叶变换(FFT),分析出的拉莫尔进动频率为3.5 kHz,对应着约500 nT的外磁场.图6(c)表示300 s时间内采集的磁场数据,磁场采样频率为10 Hz,磁场波动小于10 pT.图6(c)中插图部分显示了4 s时间内的磁场数据,原子磁力仪的磁场分辨率为0.1 pT.图6(d)是由图6(c)中磁场数据处理得到的噪声功率谱密度,用1 Hz频点附近11个数据的平均值代表原子磁力仪的灵敏度,得到灵敏度指标为0.2 pT/Hz1/2.本研究采用美国安捷伦科技公司的B2912 A型精密电流源产生待测磁场,电流源精度为10−6,当电流源输出的量程值分别为1MA,10MA,100MA,1 A时,分别对应着1 nA,10 nA,100 nA,1µA的电流分辨率.原子磁力仪测量的磁场由电流源产生,因此电流源的噪声将反映在磁力仪灵敏度指标测试中.图7显示了磁力仪灵敏度指标和线圈电流与外磁场大小的依赖关系.当I>100 MA时,磁力仪灵敏度约为12pT/Hz1/2,对应电流分辨率为1µA;当10 MA<I<100 MA时(图中阴影部分),磁力仪灵敏度约为1 pT/Hz1/2,对应电流分辨率为100 nA;当1MA<I<10 MA时,磁力仪灵敏度约为0.2 pT/Hz1/2,对应电流分辨率为10 nA;特殊地,当I<1 MA时,在50 nT磁场环境中磁力仪的灵敏度依旧为0.2 pT/Hz1/2,此时对应电流分辨率为1 nA.综上所述,本文的NMOR铷原子磁力仪的极限灵敏度为0.2 pT/Hz1/2.图7中线圈电流I与外磁场B在1 MA附近呈现非严格的线性关系,这是由磁屏蔽筒内的剩余磁场导致的,可参考表达式(3).3)横向弛豫时间对磁场大小的依赖关系原子系综宏观磁化强度被激励磁场作用至与外磁场垂直的平面内,该横向磁化强度将呈指数形式衰减,衰减函数的时间常数为横向弛豫时间T2,即信号幅度衰减至e−1倍所需的时间[20].本文中用y=A exp(−t/T2)函数来拟合出T2.图8(a)显示了500 nT磁场下的弛豫信号,此时原子磁力仪的工作周期T=100ms,抽运激光作用时长t1=30 Ms,激励磁场作用时长t2=5 Ms.以激励磁场关闭时为时间零点,将弛豫信号中的波峰随时间的变化曲线绘制在图8(b)中,通过指数拟合得到横向弛豫时间T2为5.946 Ms.图8(c)显示了横向弛豫时间随磁场的变化,可以看出随着磁场的增加,横向弛豫时间逐渐减小,这是由于铷泡所在区域磁场梯度的增加导致了原子系综宏观磁化强度的弛豫加快.图8(c)的实验结果对Labview程序编写时自由弛豫信号截取时长的设定具有重要参考意义.4)磁场采样率磁场采样率S是原子磁力仪的一项重要指标.目前国内光泵磁力仪磁场采样率大都小于20 Hz,而国外已出现磁场采样率为100 Hz、甚至1000 Hz的原子磁力仪[18].例如美国Geometrics公司推出的G-824 A型航空铯磁力仪的采样率达到了1000 Hz,而美国限制出口该磁力仪[18].本文的NMOR铷原子磁力仪通过设定工作周期T、抽运激光作用时长t1、激励磁场作用时长t2,可实现磁场采样率S在1—1000 Hz范围内可调.实验中当以1000 Hz磁场采样率测量10000 nT附近的恒磁场时,90%的数据落在(10000±0.1)nT以内.高磁场采样率的磁力仪可用于测量环境中低频的交变磁场,图9显示了原子磁力仪测量(10000±100)nT范围内频率为100 Hz交变磁场的实验结果,测量时激励磁场振荡频率固定为70 kHz.图9(a)是原子磁力仪采集的原始数据,随着磁场的波动原始信号的最大振幅也跟着波动;图9(b)是原子磁力仪时序示意图,设定工作周期T=1 Ms,抽运激光作用时长t1=0.3 ms,激励磁场作用时长t2=0.1 ms;图9(c)显示了测量的磁场数据.NMOR铷原子磁力仪的拉莫尔进动频率是由自由弛豫正弦信号的快速傅里叶变换曲线拟合得到,因此磁场采样率S的设定需要考虑与拉莫尔进动频率相适应,必须保证有足够多的数据能够精确拟合出拉莫尔进动频率.本文原子磁力仪在5000—100000 nT待测磁场范围内实现磁场采样率S在1—1000 Hz范围内可调,在100—5000 nT待测磁场范围内可设定S≤20 Hz.另外,本文描述的原子磁力仪在高磁场采样率条件下无法使用跟踪式锁频,这是因为跟踪式锁频步骤是在Labview程序中实现,而在程序流程中计算机与信号源通讯需要时间,采用跟踪式锁频测量时S≤20 Hz.信号源输出频率为定值时磁场采样率S可在1—1000 Hz范围内可调,参考图4(a)的实验结果,适用于测量稳定磁场附近小于1000 nT的磁场波动.本文详细地描述了NMOR铷原子磁力仪的工作原理和测量方法,系统地研究了测量范围、灵敏度、分辨率、横向弛豫时间、磁场采样率等性能指标.实验结果表明原子磁力仪测量范围为100—100000 nT,极限灵敏度为0.2 pT/Hz1/2,磁场分辨率为0.1 pT,制备的铷原子极化态横向弛豫时间在毫秒量级,磁场采样率最高可达1000 Hz.本文用噪声功率谱密度讨论原子磁力仪的灵敏度指标时考虑了精密电流源的电流噪声,该做法对磁力仪的灵敏度指标标定具有借鉴意义.本文原子磁力仪的若干性能指标在国内以及国际上都具有先进性.除了上述列出的性能指标外,磁力仪的空间分辨率也是磁力仪的一项重要指标,而本研究采用Φ25 mm×50 mm的圆柱型气室,体积较大,下一步可研究微型原子气室的原子磁力仪.本研究的原子磁力仪在生物医学、基础物理研究方面具有潜在的应用前景.本文所描述的原子磁力仪实验装置是在浙江工业大学林强教授及其团队老师吴彬、郑文强、程冰,以及浙江科技学院李曙光副教授的帮助下搭建完成的,上述研究人员在作者搭建原子磁力仪过程中给予了诸多技术资料、技术协助和有益讨论.作者本人现场参观了浙江工业大学的原子磁力仪装置,从中获得启发,完成了本文的研究内容.作者对林强教授团队表示由衷的感谢.We report a rubidiuMatoMicMagnetoMeter based on puMp-probe nonlinearMagneto-op tical rotation.The rubidiuMvapor cell is p laced in a five-layer Magnetic shield With inner coils that can generate uniforMMagnetic fields along the direction of puMp beam,and the cell is also p laced in the center of a Helmholtz coil that can generate an oscillating Magnetic field perpendicular to the direction of puMp beam.The atoMs are op tically puMped by circularly polarized puMp beaMalong a constant magnetic field in a period of time,then the puMp beaMis turned off and aπ/2 pulse of oscillating magnetic field for87Rb atoMs is app lied.A fter the above p rocess,the individual atoMic magnetic moments becoMe phase coherent,resu lting in AtransverseMagnetization vector precessing at the LarMor frequency in theMagnetic field.The linearly polarized probing beaMis perpendicular to the direction ofmagnetic field,and can be seen as a superposition of the left and right circularly polarized light.Because of the diff erent absorptions and dispersions of the left and right circularly polarized light by rubidiuMatoMs,the polarization direction of p robing beaMrotateswhen probing beaMpasses through rubidiuMvapor cell.The rotation of the polarization is subsequently converted into an electric signal through a polarizing beaMsplitter.Finally,the decay signal related to the transverseMagnetization vector isMeasured.The LarMor frequency p roportional to Magnetic field isobtained by the Fourier transforMof the decay signal.The value ofmagnetic field is calculated froMthe formula:B=(2π/γ)f,where γ and f are the gyromagnetic ratio and LarMor frequency,respectively.In order toMeasure theMagnetic field in a Wide range,the tracking lock Mode is p roposed and tested.The atoMicMagnetoMeter can track themagnetic field juMp of 1000 nT or 10000 nT,indicating that the atoMicmagnetometer has strong locking ability and can be easily locked after start-up.The Main perforMances in diff erent Magnetic fields are tested.The results shoWthat the MeasureMent range of the atoMic magnetometer isfroM100 nT to 100000 nT,the extreme sensitivity is 0.2 pT/Hz1/2,and the magnetic field resolution is 0.1 pT.The transverse relaxation tiMes of the transverse Magnetization vector in diff erent Magnetic fields are obtained,and the relaxation tiMe decreases With the increase of the Magnetic field.When the MeasureMent range is froM5000 nT to 100000 nT,themagnetic field saMp ling rate of the atoMicmagnetometer can be ad justed in a range froM1 Hz to 1000 Hz.The atoMic MagnetoMeter in high saMp ling rate can Measure weak alternating Magnetic field at loWfrequency.This paper provides an iMportant reference for developing the atoMic MagnetoMeter With large measurement range,high sensitivity and high saMp ling rate.【相关文献】[1]Xu S,C raWford C W,Rochester S,Yashchuk V,Budker D,Pines A 2008 Phys.Rev.A 78 013404[2]Maser D,Pandey S,Ring H,Ledbetter MP,Knappe S,K itching J,Budker D 2011Rev.Sci.Instrum.82 086112[3]Kornack T W,Ghosh R K,RoMalis MV 2005 Phys.Rev.Lett.95 230801[4]Meyer D,Larsen M2014 Gyroscopy and Navigation 5 75[5]C leMT R 1998 Nav.Eng.J.110 139[6]Savukov IM,Seltzer S J,RoMalis MV 2005 Phys.Rev.Lett.95 063004[7]Budker D,RoMalis MV 2007 Nat.Phys.3 227[8]Savukov I M,RoMalis MV 2005 Phys.Rev.Lett.94 123001[9]Yashchuk V V,G ranwehr J,K iMball D F,Rochester S M,Trabesinger A H,U rban JT,Budker D,Pines A 2004 Phys.Rev.Lett.93 160801[10]Liu G B,Sun X P,Gu S H,Feng JW,Zhou X 2012 Physics 41 803(in Chinese)[刘国宾,孙献平,顾思洪,冯继文,周欣2012物理41 803][11]A llred J C,LyMan R N,Kornack T W,RoMalis MV 2002 Phys.Rev.Lett.89 130801[12]KoMinis I K,Kornack T W,A llred J C,RoMalis MV 2003 Nature 422 596[13]Dang H B,Maloof A C,RoMalis MV 2010 Appl.Phys.Lett.97 151110[14]Li S G,Zhou X,Cao X C,Sheng J T,Xu Y F,Wang Z Y,Lin Q 2010 Acta Phys.Sin.59 877(in Chinese)[李曙光,周翔,曹晓超,盛继腾,徐云飞,王兆英,林强2010物理学报59 877][15]Gu Y,Shi R Y,Wang Y H 2014 Acta Phys.Sin.63 110701(in Chinese)[顾源,石荣晔,王延辉2014物理学报63 110701][16]D ing Z C,Li Y Y,Wang Z G,Yang K Y,Yuan J 2015 sers 42 0408003(in Chinese)[丁志超,李莹颖,汪之国,杨开勇,袁杰2015中国激光42 0408003][17]Wang Z G,Luo H,Fan Z F,Xie Y P 2016 Acta Phys.Sin.65 210702(in Chinese)[汪之国,罗晖,樊振方,谢元平2016物理学报65 210702][18]Dong H B,Zhang C D 2010 Chin.J.Eng.Geophys.7 460(in Chinese)[董浩斌,张昌达2010工程地球物理学报7 460][19]Wang Y Q,Wang Q J,Fu J S,Dong T Q 1986 The Theory of FrequencyStandards(Beijing:Science Press)pp168–173(in Chinese)[王义遒,王庆吉,傅济时,董太乾1986量子频标原理 (北京:科学出版社)第168—173页][20]Ek lund E J 2008 Ph.D.D issertation(USA:University of California Irvine)PACS:07.55.Ge,32.60.+i,32.80.Xx,42.50.Gy DOI:10.7498/aps.66.160701†Corresponding author.E-Mail:*******************。

圆锥曲线论 英文

圆锥曲线论 英文

圆锥曲线论英文Conic sections are a fundamental topic in mathematics that deals with the properties and equations of curves formed by the intersection of a plane with a cone. These curves include the circle, ellipse, parabola, and hyperbola. In this document, we will explore the characteristics and equations of these conic sections.1. Circle:A circle is a conic section formed when a plane intersects a cone at a right angle to its axis. It is defined as the set of all points in a plane that are equidistant from a fixed center point. The equation of a circle with center (h, k) and radius r is given by (x h)^2 + (y k)^2 = r^2.2. Ellipse:An ellipse is formed when a plane intersects a cone at an angle that is less than a right angle. It is defined as the set of all points in a plane for which the sum of the distances from two fixed points (called foci) is constant. The equation of an ellipse with center (h, k), major axis length 2a, and minor axis length 2b is given by ((x h)^2 / a^2) + ((y k)^2 / b^2) = 1.3. Parabola:A parabola is formed when a plane intersects a cone parallel to one of its generating lines. It is defined as the set of all points in a plane that are equidistant from a fixed point (called the focus) and a fixed line (called the directrix). The equation of a parabola with vertex (h, k) and focal length p is given by (x h)^2 = 4p(y k).4. Hyperbola:A hyperbola is formed when a plane intersects a cone at an angle greater than a right angle. It is defined as the set of all points in a plane for which the absolute value of the difference of the distances from two fixed points (called foci) is constant. The equation ofa hyperbola with center (h, k), transverse axis length 2a, and conjugate axis length 2b is given by ((x h)^2 / a^2) ((y k)^2 / b^2) = 1.These equations provide a mathematical representation of the conic sections and allow us to analyze their properties. By manipulating these equations, we can determine important characteristics such as the shape, size, orientation, and position of the conic sections.In addition to their geometric properties, conic sections have various applications in different fields. For example, circles are commonly used in geometry, physics, and engineering to represent objects with rotational symmetry. Ellipses are used in astronomy to describe the orbits of planets and satellites. Parabolas are used in physics to model the trajectory of projectiles and in engineering to design reflectors and antennas. Hyperbolas are used in physics and engineering to describe the behavior of waves and particles.In conclusion, conic sections are a fascinating topic in mathematics with diverse applications in various fields. Understanding the properties and equations of circles, ellipses, parabolas, and hyperbolas allows us to analyze and solve problems involving these curves. By studying conic sections, we gain valuable insights into the fundamental principles of geometry and their practical applications.。

A new paradigm for the universe (preliminary version)

A new paradigm for the universe (preliminary version)

a rXiv:as tr o-ph/31133v13Nov23A new paradigm for the universe (preliminary version)Colin Rourke Email:cpr@ URL:/~cpr Abstract A new paradigm for the structure of galaxies is proposed.The main hypothesis is that a normal galaxy contains a hypermassive black hole at its centre which generates the spiral arms.The paradigm gives satisfactory explanations for:•The rotation curve of a galaxy.•The spherical bulge at the centre of a normal galaxy.•The spiral structure and long-term stability of a normal galaxy.•The age and orbits of globular clusters.•The origin and prevalence of solar systems.•(Highly speculative)the origin of life.The paradigm is compatible with direct observations but not with many of the current interpretations of these observations.It is also incompatible with large swathes of current cosmological theory and in particular with the expanding universe and hot big bang theories.A tentative new explanation for the observed redshift of distant objects is given which is compatible with a static model for the universe.This paper is being circulated in a very preliminary form in the hope that others will work both on interpreting observational data in the light ofthe new paradigm and on the obvious gaps in underlying theory.AMS Classification 85A15;85A40,83F05,83C40,83C57Keywords Galaxy,new paradigm,spiral structure,black hole,rotation curve,redshift1IntroductionThis paper is concerned with the nature of normal galaxies such as the Milky Way or the Andromeda nebula.It is an attempt to formulate a theory which explains satisfatorily both the spiral structure and the rotation curve.Current explanations for these have an unnatural ad hoc nature—in particular the hy-pothesis of exactly the correct amount of unobserved dark matter to fit the rotation curve.The main hypothesis is that the centre of a normal galaxy contains a hyper-massive black hole1of mass around1015solar masses.The centre generates the spiral arms by a process,which will be outlined,whereby matter is ejected from the centre and condenses into solar systems.This implies that young stars in a galaxy are moving outwards as well as around the centre.This gen-eral outward movement has not been observed and the reason for this is that the frequency shift due to the outward motion is cancelled by the gravitational frequency shift from the gravitationalfield of the centre.As will be seen later, in a normal galaxy stars move outwards at near escape velocity,so the two opposing frequency shifts are almost the same.Also the motion of stars is far from Keplerian,being strongly controlled by inertial drag effects from the (rotating)centre.The result is that the outward progress takes a very long time—commensurate with the lifetime of a star and hence the outward velocity is rather smaller than(about one tenth of)the observed rotational velocity. The general picture which emerges is of a structure stable over an extremely long timescale(at least1012years)with stars born and aging on their outward journey from the centre and returning to the centre to be recycled with new matter to form new solar systems.This timescale is incompatible with current estimates for the age of the universe and entails the abandonment of the big bang theory in its current form.The tentative suggestion is that galaxies have a natural lifetime of perhaps1016years with the universe considerably older than this.The paper is organised as follows.Section2contains a discussion of galactic rotation curves.This is not the most logical place to start,but it provides by far the clearest evidence for the main hypothesis.The new paradigm gives a natural explanation of the observed rotation curves for galaxies.Section3 discusses another strong piece of evidence—the spherical bulge at the centre of normal galaxies.Section4describes the main proposal,namely the generator for spiral arms and the way in which energy feeds from the central black hole. In section5the full dynamic of a spiral galaxy is discussed and the way in which the arms are formed in differing types of spiral galaxies.Section6is concerned with a further strong piece of evidence,namely the age and orbits of globular clusters.Section7discusses the formation of solar systems,which appear in a natural way in the new paradigm and tentatively suggests how life starts on suitable planets.Section8concerns direct evidence including evidence from our own galaxy near Sagittarius A∗and section9contains the suggested new explanation for redshift.Finally section10contains speculative material on thenature and long-term evolution of galaxies.Three estimates for the mass of a normal galaxy are given(in sections2,3and 9)and several observational predictions are made.It is worth remarking that,in contrast to other suggested alternatives to current mainstream cosmology,this paper does not propose any new physics apart from one mild hypothesis(with supporting plausibility arguments)about the behaviour of light in negatively curved space-time;note that this hypothesis is not needed for the main content of the paper—the new model for galaxies. Indeed the whole paper can be seen as a strong supporting evidence for the correctness of standard Einsteinian relativity.The paper is a very preliminary study of the subject and has obvious gaps, notably in the interpretation of observations and in the theory of the spiral arm generator.It is being circulated in the hope that others will help complete the work.2The rotation curveThe observed rotation curves for galaxies are quite striking.Essentially the curve(of tangential velocity against distance from the centre)comprises two approximately straight lines with a short transition region.Thefirst line passes through the origin and the second is horizontal.For a typical example,enter “galactic rotation curve”into Google,hit“I’m feeling lucky”and look at the figure on the left near the bottom.2Galactic rotation curves are so characteristic(and simple to describe)that there must be some strong structural reason for them.They are very far indeed from the curve you would get with a standard Keplerian model of rotation with any reasonable mass distribution.The current explanation involves a fortuitous arrangement of“dark matter”(ie matter for whose existence there is no other evidence)and begs several questions not least of which are the stability and prevalence of this arrangement.This explanation is strained to the limits by several observations which show that the horizontal straight line section of the rotation curve extends far outside the limits of the main visible part of galaxies.The explanation given here involves no dark matter.Essentially the rotation curve is a consequence of inertial drag due to rotation of the hypermassive cen-tral black hole.Inertial drag is one of the stranger consequences of general relativity.It is a true embodiment of Mach’s principle understood in the sense that the matter in the universe determines the concept of inertia.For a discus-sion of this effect see Misner,Thorne and Wheeler[8;section21.12].Only the most basic properties of inertial drag will be needed:(1)A rotating body causes the local concept of“inertial frame”to rotate inthe same sense as its own rotation.(2)This effect drops offproportionately with1rω+1×0r +1=ωK3The effect described here does not depend on the mass,so“heavy”is strictly un-necessary.However if the mass is not huge,the effect is negligible in practical terms. Zero size is for simplicity and is also immaterial.4A calculation of K can be deduced from the discussion in Misner,Thorne and Wheeler(op cit).See section21.12and in particular equations21.155-6:the rotation effect for a body at distance r is43.(1)The rotation of the inertial frames tends to increase v.The accelerationisωfor r small and isωKK+r so the decelleration due to this effect is:1K+r )=vr+KAdding the two effects:dvK+r −vK+r=2KωrThus:ddr +v=2rωKrThere are two asymptotes.For r small v≈rωand the curve is roughly a straight line through the origin.And for r large the curve approaches the horizontal line v=2ωK.A rough graph is given infigure1where K=ω=1.1210203040Figure1This curve is a good,but not perfect,fit to observed galactic rotation curves. However our analysis so far has been very simple minded and has ignored all effects except inertial drag from the centre.In this sense the result is quite remarkable.The salient features of the rotation curve are entirely explained as the effect of inertial drag.The most important effect that has been ignored,for a standard spiral galaxy with bilateral symmetry,is the gravitational attraction of the arms.This causes a“flywheel effect”:stars in the arms will tend to rotate with the local rotation of the arms.This extends the region where rotation is roughly plate-like and causes the rotation curve to be nearer to a straight-line through the origin of somewhat smaller slope than the asymptote.Theflywheel effect breaks down as r increases to the point where the forces required to maintain it are too great and the rotation curve turns fairly sharply towards the horizontal asymptote. This modifies the curve offigure1to something likefigure2.2110203040Figure2Finally note that perturbations to tangential velocity die out like1/r and the limiting horizontal asymptote is highly stable.Perturbation in theflywheel effect(due to non-uniform mass in the arms)will result influctuations in the horizontal portion of the rotation curve as illustrated infigure3.Figure3is a perfectfit for observed rotation curves.Several comments need to be made.(1)The analysis in this section has totally ignored radial velocity.Indeed it is a quite remarkable effect of inertial drag that the tangential velocity is controlled asymptotically independently of radial velocity or acceleration.In section5the effects on radial acceleration caused by the rotation are discussed and the full dynamic of a spiral galaxy is derived in outline.This will explain the familiar spiral structure.21203040Figure3(2)An unholy mixture of relativistic and Newtonian dynamics has been used. For largish r this is justified since inertial drag is then the dominant relativistic effect.However for small r the approximations used are probably very coarse and the analysis needs to be done properly in a fully relativistic setting which may well significantly alter the theoretical rotation curve.(3)Nevertheless this crude discussion does allow afirst estimate of the central mass of a normal galaxy.The observed rotation curves show roughly plate-like rotation extending about to the edge of the central spherical paring withfigure3,this would give the radius of the bulge as about10times the Schwarzschild radius.From looking at a gallery of galaxy photographs,the radius of the bulge is commonly about10−1times the galactic radius so the Schwarzschild radius is about10−2times the galactic radius.Thus for the Milky Way with a diameter of1018km the Schwarzschild radius is5×1015km giving a mass of roughly1.7×1015solar masses.(4)The control of tangential velocity that has been examined in this section applies most strongly to motion in the galactic plane(and indeed to stars mov-ing outward at the start of their lives).The control gets progessively weaker the further out one goes and random motions and local gravitational effects change the motion.Thus one would expect that the rotation curves obtained infigures 1to3apply best in the galactic plane,with significant variations outside it. This is indeed what is observed.As remarked earlier,the effect described in this section is independent of mass. However for rotating bodies of small mass the effect is unobservably small.For example the sun has K≈2km andω=2π/25days.Thus2Kωis4km per4 days or.04km per hour.Observational prediction Galaxies with poor bilateral symmetry will have smallerflywheel effect and the curve will be nearer tofig1thanfigs2or3.3The spherical bulgeThe next piece of evidence for the existence of a hypermassive black hole at the centre of galaxies is so obvious and commonplace that,like many common-place observations,it is easily overlooked.Normal galaxies have a pronounced spherical bulge at their centres.No satisfactory explanation for this has been proposed.If a galaxy is a rotating disc composed of stars and gas with,per-haps,a massive but not hypermassive black hole at the centre(say107)solar masses,then there is no reason to expect the formation of a spherical bulge. One might see a pronounced cluster at the centre,but why should this extend to great distances on either side of the plane of the galaxy?However if there is a hypermassive black hole at the centre then a spherical bulge is exactly what would be seen,because of gravitational lensing effects. The bulge is not real,but an artifact of the distortion of light caused by the black hole.A graphic demonstration of this effect can be found on the web at: http://www.photon.at/~werner/bh/The visible size of the bulge can be used to give a second rough estimate of the mass of the central black hole similar to that given in section2.The lensing effects extend to roughly20times the Schwarzschild radius.However the expected size of the bulge depends heavily on the actual shape of the galaxy, thus for a given mass of centre,a thin disc will not give such a large bulge as a thicker one.A guess from looking at pictures is that the Schwarzschild radius of the black hole at the centre is commonly about3×10−3times the galactic diameter.Thus for the Milky Way with a diameter of1018km the Schwarzschild radius is3×1015km giving a mass of roughly1015solar masses. Observational prediction Repeated(and distorted)observations of distant ob-jects due to the same lensing effect that causes the apparent bulge.4The generator for spiral armsHere is a sketch of the proposed nature of a normal galaxy.The centre contains all but a small proportion(less than.1%)of the mass.The remainder(1011–1012solar masses)is the visible part of the galaxy.The centre comprises two parts.A central black hole and a surrounding rapidly rotating sphere of matter some of which is in plasma form,which I shall call the generator.The generator is mostly concentrated in the equatorial band forming a rotating toroidal belt which I shall call simply the belt.The belt has a complex electro-magnetic structure similar to that modelled(see eg Williams[12])for(lighter)blackholes and used to explain Active Galactic Nuclei.5The generator is fed from two sources.Dying stars and debris fall into it and get torn apart by the huge tidal forces and broken down into atoms or smaller particles and energy feeds directly into it from the central black hole both in the form of Penrose-process energy and directly from the gravitationalfield by tidal effects.The result becomes highly unstable as it builds up and it forms sharp bulges which explode outwardsflinging elementary particles,energy and heavier particles out into space to condense into solar systems and form the familiar arms.The generator is highly massive and stratified,with a plasma of small particles at the inside,where the input of energy from the black hole is greatest,and with layers of heavier particles and dust as one moves outwards.The thickness of the generator implies that the polar radiation observed for so-called Active Galactic Nuclei does not escape and explains why this is not observed for normal galaxies.The explosions do not occur in random places:most normal galaxies have a pronounced bilateral symmetry with two main opposing arms(eg M101,M83 etc).Why does this happen?The suggestion is that this is simply a stable solution.Once two arms have formed,then the gravitational pull of these arms will form bulges at the roots of the arms and encourage explosions there which feed the arms.The bilateral symmetry arises because the bulges are tidal bulges which always have bilateral symmetry.This tendency to bilateral structure is weak and looking at a gallery of galaxies you canfind many examples where it fails to form or where other weak arms have formed as well as the two main arms.One general observed property of spiral arms is worth commenting on.The roots of the arms are offset.A very clear example is M83(Southern Pinwheel) but any gallery of galaxies shows the same phenomenon repeated.This property is related to the apparent size of the central black hole and the belt rotation.A full explanation will have to wait for a good mathematical model for the generator.But this offset,which is clearly a real phenomenon,can be used to explain the general rotation of galaxies:The offset jets contain a deal of lighter particles,eg photons,which are radiated away from the galaxy and this “wind”radiates angular momentum away from the system and causes rotation in exactly the same way that a pinwheel ter(in section10)it will be seen how the rotation stabilises and this will explain the near uniformity of rotation across different galaxies as observed in rotation curves.5The full dynamicThis section is in very preliminary form and contains a similar unholy mixture of relativistic and Keplerian dynamics to that used in section2.Using the new paradigm,the full dynamics of a galaxy are capable of being modelled accurately with a fully relativistic treatment(for example by using the models for orbits in Kerr spacetime given in[4])and this will be done in the next version.What follow are heuristic arguments which establish the scene very roughly.A fully relativistic model may well give a significantly different picture.Most of the visible part of a galaxy lies within100Schwarzschild radii of the central black hole and the part of pre-stars’(and hence stars’)orbits which determine the overall picture is the actual ejection from the generator,which is very close indeed to the central black hole.Thus the dynamics are highly relativised and under strong central control including in particular the inertial dragfields used in section2to explain the rotation curve.The orbit of a star takes it from the centre to the outside and back in roughly its lifetime.In other words the orbit diameter is roughly the radius of the galaxy(say1018km)and its period is roughly1010years.A purely Newtonian estimate of the mass of the centre which supports such an orbit is1010solar masses.6The discrepancy from the estimates in sections2and3of roughly 1015solar masses is entirely due to relativistic effects.The strongest effect is the swan-neck effect caused by the extreme curvature near the centre.A good illustration is to be found in Misner,Thorne and Wheeler[8;figure21.12 and the text in section23.8].This effect causes radial distance to be grossly distorted near the black hole so that the real distance a star travels is far larger than the non-relativised distance.Since distance appears as d3in the calculation of M,this is highly effective in correcting the discrepancy.Moreover the effect is arbitrarily large(depending where the motion starts)and so is capable of correcting all the discrepancy.There are two other effects which are worth mentioning briefly:(1)Theflywheel effect revistedThis is not a relativistic effect.As was seen in section2,the tangential motion of stars in the galactic plane is stabilised by gravitational attraction within thet2”where M is the ratio of centralmasses,d the ratio of orbital diameters and t the ratio of orbital periods.Both d and t are1010approx.arms.This also has an effect on radial motion.The outward gravitational pull of the arms has a tendency to stabilise the outward velocity,in other words to make it more uniform than it would otherwise be.(2)The slingshot effectThis is another effect of inertial drag and is caused by the decrease of inertial drag as r increases.It is only significant for large r.The best way to think of this is to image that inertial drag causes plate-like rotation of a particle. If the effect were to stop suddenly,the tangential velocity would throw the particle outwards.As the effect decreases there is a corresponding outwards acceleration.To quantify this let v=v rot+v inert where v rot is the tangential velocity due to rotation of the local inertial frame and v inert is the tangential velocity measured in the local inertial frame.The effective outward acceleration is v2inert/r(the familiar“centrifugal acceleration”).For large r,v rot≈v inert≈Kω(see section2)and the corresponding outward acceleration is K2ω25ωand this is the apparentfixed frame as far asthe shape of the galaxy is concerned.Working in this frame,˙r=q say(the apparent outward velocity)and˙θ=0for r≤10K and˙θ=−ωr for r≥10K.Henceθ=constant(zero say)for r≤10K andθ=15+2Kωlog(r)+C)for r≥10K where C is a constant such thatθ=0for r=10K ie C= 2ωK(1−log(10K)).Without the log term this is a standard“logarithmic spiral”.The effect of the log term is to decrease the pitch of the spiral as r increases and the outwarddecelleration has the same effect.Rather than sketching this curve,two typical examples of galaxies are reproduced infigure4,both of whichfit the curve extremely well.Thefirst NGC1365is a typical barred galaxy.The transition between the two intervals for r is very clear.The pitch change in the arms is obscurred by the angle of view.The second M51is a typical complete open spiral with the plate region coinciding with the central bulge and very clear decrease in spiral pitch as in the theoretical curve.Figure4:NGC1365and M51The various observed shapes of spiral galaxies can be explained by adjusting the basic picture.The key variables are the ratio of the radius of the plate region r0=10K to the overall visible disc radius,R say,which determines the overall extent of the spiral and the size of q compared to the other constants, which determines the pitch of the spiral.For a standard spiral galaxy such as M51,M83or M101,r0is one or two orders of magnitude smaller than R and the arms wrap right round the centre.For a bar galaxy,R is not much bigger than r0and the central rotating plate forms most of the galaxy with the arms moving out only a little further.Theflywheel effect dominates for this type of galaxy.Indeed for some bar spiral galaxies,theflywheel effect has captured stellar material in advance of the arms to form apparently forward pointing arms as well as the usual trailing arms.A good example of a barred galaxy with very full data can be found in[7](NGC5383)showing constant rotation along the whole extent of the bar.The data there also show outward motion along the arms and hence directly support the main hypothesis of this paper. As in section2,the analysis of this section applies only to stars moving outward in the galactic plane.The control gets progessively weaker the further out one goes and random motions and local gravitational effects change the motion. Furthermore,the orbits being examined here are in Kerr spacetime and theseare known to be chaotic Hartl[4].7As stars age their orbits move out of the plane andfill the whole of the galactic halo which leads naturally to the subject of globular clusters.6Globular clustersAnother strong piece of evidence for the main hypothesis is provided by globular clusters which are known to be very old objects and have caused some heartache about the age of the universe according to the big-bang theory.In the new paradigm there is a very natural explanation for their age.The inner arms of a galaxy comprise young stars and star formation regions. As the stars move out along the arms they age and mature.Thus for example our sun was formed about2×109years ago and has moved about37Hartl proves this for spinning particles.The spin of stars may not significant,but local random effects are likely to cause similar chaotic effects.stripped down into small elementary particles.However some of the heavier particles will survive this destruction.Here is a tentative suggestion for how this happens.As supposed in section4,the belt has a complex stratified structure with a plasma of small particles at the inside,where the input of tidal energy from the black hole is greatest,and with layers of heavier particles as one moves outwards.The belt may even be surrounded by an orbiting region of cold dust. There is observational evidence for significant amounts of dust near the centre. As stars feed into the belt,part of their mass will mix with the outer layers and thus some of the heavier atoms resulting from fusion will survive.Thus, when the explosions which create the spiral arms occur,the outgoing matter will comprise a good deal of small particles from the inner layers(which condense to form the observed mix of light elements)8together with a sprinkling of heavier matter.When the resulting clouds condense into stars,the heavier matter will condense around them to form solar systems.Some tentative remarks on life.The timescale for the galaxy is huge.Probably several orders of magnitude greater than current estimates of the age of the universe.This is plenty of time for life to have arisen many times over on suitable planets.When these planets are destroyed(as their sun falls into the centre)some of the molecules may survive and and be recycled into new planets. Thus in a steady state planets will start out seeded with molecules which will help to start life over again.Indeed standard selection processes over a galactic timescale will favour lifeforms which can arise easily from the debris left over from their ultimate demise in the galactic centre.This might explain how life arose on earth rather more quickly than totally random processes can explain. Observational predictions Solar systems around most if not all stars. (Tentative)pre-life molecules on all planets.8Direct evidenceMany observations have been referred to in previous sections,all of which give indirect evidence for the main hypothesis.This section is concerned with direct evidence.There is direct evidence for a galactic centre of mass greater than1010solar masses[6].This is considerably smaller than the hypothesised mass for a normal galaxy,but it does show that hypermassive black holes are likely to exist.A recent paper[11]gave an upper limit for the mass of the centre of the Milky Way of around4×106solar masses.Since this is direct contradiction with the main hypothesis of this paper,it is necessary to examine this evidence in detail. The team observed a star“S2”near Sgr A∗over ten years.The motionfits a roughly Keplerian orbit of period15years or so around Sgr A∗.Their data is reproduced infigure5,but with another suggested motion superimposed,which will be explained.Figure5The suggested Keplerian orbit is not wellfitted by the last few observations. S2is moving very fast(as would be expected)as it passes near to Sgr A∗, but continues to move very fast as it moves away.Read carefully the dates on the observations.In a Keplerian orbit it would slow down symmmetrically to the acceleration on approach.The last four observations arefitted perfectly by a straight line exactly parallel to the galactic equator(exactly along the belt with the correct sense of rotation).The suggestion is that these observations do not show a closed orbit but rather they show S2falling into the centre and being captured by the belt.The web site referred to in reference[11],gives some other evidence.Click on the link“Proper motions”on the left and look at the image halfway down the page.There is another star“S1”which also shows clear evidence of acceleration towards the galactic equator(the straight line in figure5)which wouldfit with imminent capture by the belt.Observational predictions S2is either continuing to move fast along the belt line or has disappeared from view.9Stars near the galactic centre are likely。

Exact conserved quantities on the cylinder I conformal case

Exact conserved quantities on the cylinder I conformal case
DCPT-02/69 EMPG-02-20 May 2003
arXiv:hep-th/0211094v3 23 Jul 2003
Exact conserved quantities on the cylinder I: conformal case.
ቤተ መጻሕፍቲ ባይዱ
D. Fioravanti
a b
a
and M. Rossi
it became well known that these local integrals of motion represent one series of local commuting charges of conformal field theories without extended symmetries [18, 19] (cfr. [8, 1, 20] for details concerning the matrix Lax formulation). In addition, this series is exactly that which is still conserved, after suitable modification, if the conformal field theory is perturbed by its Φ(1,3) operator. Actually, applying the investigations of (2) the present paper and of [8] to the setup elaborated in [19] about the A2 KdV theory, it is only a matter of calculations to extend the subsequent results to the only other series [18], characterising the only other integrable perturbations. The very intriguing feature of the conformal formulation in terms of integrable hierarchies is that the conformal monodromy matrix already contains the perturbing field and the screening operator [19]. A similar discretisation should be also allowed in conformal field theories with extended symmetry algebra (e.g. for W algebras it is of basic importance the setting of [21]; cfr. [22] for development in a peculiar case). We remark that similar results have been obtained in pioneering works [1, 2] by Bazhanov, Lukyanov and Zamolodchikov from a different starting point and via a different approach. They define directly in the continuous field theory the transfer matrix as an operator series acting on Virasoro modules (for 0 < β 2 < 1/2). On the other hand, we intend to create a bridge between conformal field theories and the powerful algebraic Bethe Ansatz formulation of lattice KdV [8, 9] – based on this interesting generalisation of the Yang-Baxter algebra –, showing that the continuum field theory limit can be recovered in the braided case as well. As particular consequence, in the second paper we will be able to formulate the Φ(1,3) perturbation of conformal field theory within the same framework of [8] by merging two KdV theories. In any case, the monodromy matrix of [1, 2] is different from the na¨ ıve continuous limit [8] of our monodromy matrices, which, unlike in [1, 2], are not solutions of the usual Yang-Baxter equation. In section 2 we summarise the main results on (m)KdV theories obtained in [8]. In section 3, starting from the Bethe equations of [8], we write the nonlinear integral equation for each vacuum of the left and right quantum KdV equations. In section 4 we use this equation to calculate up to quadratures the energy of the vacua. In section 5 we apply the nonlinear integral equation to find exact expressions for the continuous limit of the vacuum eigenvalues of the transfer matrix. From the asymptotic expansion of these eigenvalues we obtain the local abelian charges, mentioned before as characterising the Φ(1,3) perturbation. In addition, we have solved exactly the theory in two cases: at 1 the free fermion point β 2 = 2 and in the limit of infinite twist (˜ ω → +∞). Some results are also compared with results and conjectures in [1, 2]. In Section 6 we summarise our work describing also possible further applications.

芯片测试仪设备操作指南说明书

芯片测试仪设备操作指南说明书

MechanicsOscillationsElliptical Oscillation of a String PendulumDESCRIPTION OF ELLIPTICAL OSCILLATIONS OF A STRING PENDULUM AS THE SU-PERIMPOSITION OF TWO COMPONENTS PERPENDICULAR TO ONE ANOTHER.UE1050121 06/15 MEC/UDFig. 1: Experiment set-upGENERAL PRINCIPLESDepending on the initial conditions, a suitable suspended string pendulum will oscillate in such a way that the bob’s motion describes an ellipse for small pendulum deflections. If the motion is resolved into two perpendicu-lar components, there will be a phase difference between those components.This experiment will investigate the relationship by measuring the oscillations with the help of two perpendicularly mounted dynamic force sensors. The amplitude of the components and their phase difference will then be evaluated. The phase shift between the oscillations will be shown directly by displaying the oscillations on a dual-channel oscilloscope.Three special cases shed light on the situation:a) If the pendulum swings along the line bisecting the two force sensors, the phase shift φ = 0°.b) If the pendulum swings along a line perpendicular to that bisecting the two force sensors, the phase shift φ = 180°.c) If the pendulum bob moves in a circle, the phase shift φ =oscillation directions of the string pendulum under in-vestigationLIST OF EQUIPMENT1 SW String Pendulum Set 1012854 (U61025)1 SW Stand Equipment Set 1012849 (U61022)1 SW Sensors Set @230 V 1012850 (U61023-230) or1 SW Sensors Set (@115 V 1012851 (U61023-115) 1 USB Oscilloscope 2x50 MHz 1017264 (U112491)SET-UP∙Screw the stand rods with both external and internal threads into the outer threaded sockets of the base plate. ∙Extend both rods by screwing rods with external thread only onto the ends of them.∙Attach double clamps near the top of both stand rods and turn them to point inwards so that the slots are vertical and facing one another.∙Attach both springs from the spring module to the lugs on the cross bar (angled side).∙Hang the large loop of string from the lug on the flat side.Fig. 3 Assembly of spring module∙Connect the springs and vector plate to the hook of a dynamic force sensor with a small loop of string and care-fully pull everything taut.∙Attach the force sensor with the screw tightened by hand. ∙Attach the second force sensor in the same way.Fig. 4 Attachment of dynamic force sensors to spring module∙Pull the string through the eyelet of the spring module (in the middle of the metal disc).∙Thread the end of the string through the two holes of the length adjustment slider.Fig. 5 Set up of string3B Scientific GmbH, Rudorffweg 8, 21031 Hamburg, Germany, ∙Clamp the cross bar into the slots of the two double clamps, suspend a weight from the end of the string and set up the height of the pendulum using the length ad-justment slider.Fig. 6 Attachment of cross bar in double clamp ∙ Connect the force sensors to the inputs for channels A and B of the MEC amplifier board.∙ Connect outputs A and B of the MEC control unit to channels CH1 and CH2 of the oscilloscope.EXPERIMENT PROCEDURE∙Set the oscilloscope time base time/div to 1 s, select a vertical deflection for channels CH1 and CH2 of 50 mV DC and set the trigger to “Edge” mode, “Normal” sweep, “Source CH1” and “Slope +”.∙Slightly deflect the string pendulum and allow it to oscil-late in a plane which bisects the alignment of the two force sensors (oscillation path a in Fig. 2). Observe the oscilloscope trace and save it.∙Slightly deflect the string pendulum and allow it to oscil-late in a plane which is perpendicular to the one which bi-sects the two force sensors (oscillation path b in Fig. 2). Observe the oscilloscope trace and save it.∙Slightly deflect the string pendulum and allow it to oscil-late in a circle (oscillation path c in Fig. 2). Observe the oscilloscope trace and save it.SAMPLE MEASUREMENT AND EVALUA-TIONWhen the pendulum is oscillating in the plane of the bisecting angle between the sensors, the two sensors will experience symmetric loading (oscillation path a in Fig. 2). The signals from the two force sensors will be in phase, i.e. the phase shift between them will be φ= 0° (Fig. 7).Fig. 7: Oscillation components for a string pendulum swingingalong the line bisecting the two force sensorsWhen the pendulum is oscillating in the plane perpendicular to the bisecting angle between the sensors, the two sensors will experience asymmetric loading (oscillation path b in Fig. 2). The signals from the two force sensors will be wholly out of phase, i.e. the phase shift between them will be φ= 180° (Fig. 8).Fig. 8: Oscillation components for a string pendulum swingingalong the line perpendicular to that bisecting the two force sensorsThe circular oscillation is a superimposition of the oscillations along the plane of the bisecting angle between the sensors and the angle perpendicular to it with a phase shift of φ = 90°(Fig. 9).Fig. 9: Oscillation components for a string pendulum describ-ing a circle。

PROBING THE PARETO FRONTIER FOR

PROBING THE PARETO FRONTIER FOR

Department of Computer Science,University of British ColumbiaTechnical Report TR-2008-01,January2008(revised May2008)PROBING THE PARETO FRONTIER FORBASIS PURSUIT SOLUTIONSEWOUT V AN DEN BERG AND MICHAEL P.FRIEDLANDER∗Abstract.The basis pursuit problem seeks a minimum one-norm solution of an underdetermined least-squares problem.Basis pursuit denoise(BPDN)fits the least-squares problem only approximately, and a single parameter determines a curve that traces the optimal trade-offbetween the least-squares fit and the one-norm of the solution.We prove that this curve is convex and continuously differentiable over all points of interest,and show that it gives an explicit relationship to two other optimization problems closely related to BPDN.We describe a root-finding algorithm forfinding arbitrary points on this curve;the algorithm is suitable for problems that are large scale and for those that are in the complex domain.At each iteration,a spectral gradient-projection method approximately minimizes a least-squares problem with an explicit one-norm constraint.Only matrix-vector operations are required.The primal-dual solution of this problem gives function and derivative information needed for the root-finding method. Numerical experiments on a comprehensive set of test problems demonstrate that the method scales well to large problems.Key words.basis pursuit,convex program,duality,root-finding,Newton’s method,projected gradient,one-norm regularization,sparse solutionsAMS subject classifications.49M29,65K05,90C25,90C061.Basis pursuit denoise.The basis pursuit problem aims tofind a sparse solution of the underdetermined system of equations Ax=b,where A is an m-by-n matrix and b is an m-vector.Typically,m n,and the problem is ill-posed.The approach advocated by Chen et al.[15]is to solve the convex optimization problemx 1subject to Ax=b.(BP)minimizexIn the presence of noisy or imperfect data,however,it is undesirable to exactlyfit the linear system.Instead,the constraint in(BP)is relaxed to obtain the basis pursuit denoise(BPDN)problemx 1subject to Ax−b 2≤σ,(BPσ)minimizexwhere the positive parameterσis an estimate of the noise level in the data.The case σ=0corresponds to a solution of(BP)—i.e.,a basis pursuit solution.There is now a significant body of work that addresses the conditions under which a solution of this problem yields a sparse approximation to a solution of the underdetermined system;see Cand`e s,Romberg,and Tao[11],Donoho[24],and Tropp[48],and references therein.The sparse approximation problem is of vital importance to many applications in signal processing and statistics.Some important applications include image reconstruction,such as MRI[36,37]and seismic[31,32] images,and model selection in regression[26].In many of these applications,the data sets are large,and the matrix A is available only as an operator.In compressed sensing[10–12,23],for example,the matrices are often fast operators such as Fourier or wavelet transforms.It is therefore crucial to develop algorithms for the sparse reconstruction problem that scale well and work effectively in a matrix-free context.∗Department of Computer Science,University of British Columbia,Vancouver V6T1Z4,B.C., Canada({ewout78,mpf}@cs.ubc.ca).Friedlander is corresponding author.This work was supported by the NSERC Collaborative Research and Development Grant334810-05.May30,200812 E.van den BERG and M.P.FRIEDLANDERWe present an algorithm,suitable for large-scale applications,that is capable of finding solutions of (BP σ)for any value of σ≥0.Our approach is based on recasting (BP σ)as a problem of finding the root of a single-variable nonlinear equation.At each iteration of our algorithm,an estimate of that variable is used to define a convex optimization problem whose solution yields derivative information that can be used by a Newton-based root finding algorithm.1.1.One-norm regularization.The convex optimization problem (BP σ)is only one possible statement of the one-norm regularized least-squares problem.In fact,the BPDN label is typically applied to the penalized least-squares problem (QP λ)minimize xAx −b 22+λ x 1,which is the problem statement proposed by Chen et al.[14,15].A third formulation,(LS τ)minimize x Ax −b 2subject to x 1≤τ,has an explicit one-norm constraint and is often called the Lasso problem [46].The parameter λis related to the Lagrange multiplier of the constraint in (LS τ)and to the reciprocal of the multiplier of the constraint in (BP σ).Thus,for appropriate parameter choices of σ,λ,and τ,the solutions of (BP σ),(QP λ),and (LS τ)coincide,and these problems are in some sense equivalent.However,except for special cases—such as A orthogonal—the parameters that make these problems equivalent cannot be known a priori.The formulation (QP λ)is often preferred because of its close connection to convex quadratic programming,for which many algorithms and software are available;some examples include iteratively reweighted least squares [7,section 4.5]and gradient projection [27].For the case where an estimate of the noise-level σis known,Chen et al.[15,section 5.2]argue that the choice λ=σ√2log n has important optimality properties.However,this argument hinges on the orthogonality of A .We focus on the situation where σis approximately known—such as when we can estimate the noise levels inherent in an underlying system or in the measurements taken.In this case it is preferable to solve (BP σ),and here this is our primary goal.An important consequence of our approach is that it can also be used to efficiently solve the closely related problems (BP )and (LS τ).Our algorithm also applies to all three problems in the complex domain,which can arise in signal processing applications.1.2.Approach.At the heart of our approach is the ability to efficiently solve a sequence of (LS τ)problems using a spectral projected-gradient (SPG)algorithm [5,6,18].As with (QP λ),this problem is parameterized by a scalar;the crucial difference,however,is that the dual solution of (LS τ)yields vital information on how to update τso that the next solution of (LS τ)is much closer to the solution of (BP σ).Let x τdenote the optimal solution of (LS τ).The single-parameter functionφ(τ)= r τ 2with r τ:=b −Ax τ(1.1)gives the optimal value of (LS τ)for each τ≥0.As we describe in section 2,its derivative is given by −λτ,where λτ≥0is the unique dual solution of (LS τ).Importantly,this dual solution can easily be obtained as a by-product of the minimization of (LS τ);thisPROBING THE PARETO FRONTIER FOR BASIS PURSUIT SOLUTIONS3 is discussed in section2.1.Our approach is then based on applying Newton’s method tofind a root of the nonlinear equationφ(τ)=σ,(1.2)which defines a sequence of regularization parametersτk→τσ,where xτσis a solutionof(BPσ).In other words,τσis the parameter that causes(LSτ)and(BPσ)to share the same solution.There are four distinct components to this paper.Thefirst two are related to the root-finding algorithm for(BPσ).The third is an efficient algorithm for solving (LSτ)—and hence for evaluating the functionφand its derivativeφ .The fourth gives the results of a series of numerical experiments.Pareto curve(section2).The Pareto curve defines the optimal trade-offbetween the two-norm of the residual r and the one-norm of the solution x.The problems (BPσ)and(LSτ)are two distinct characterizations of the same curve.Our approach uses the functionφto parameterize the Pareto curve byτ.We show that for all points of interest,φ—and hence the Pareto curve—is continuously differentiable.We are also able to give an explicit expression for its derivative.This surprising result permits us to use a Newton-based root-finding algorithm tofind roots of the nonlinear equation (1.2),which correspond to points on the Pareto curve.Thus we canfind solutions of (BPσ)for anyσ≥0.Rootfinding(section3).Each iteration of the root-finding algorithm for(1.2) requires the evaluation ofφandφ at someτ,and hence the minimization of(LSτ). This is a potentially expensive subproblem,and the effectiveness of our approach hinges on the ability to solve this subproblem only approximately.We present rate-of-convergence results for the case whereφandφ are known only approximately.This is in contrast to the usual inexact-Newton analysis[22],which assumes thatφis known exactly.We also give an effective stopping rule for determining the required accuracy of each minimization of(LSτ).Projected gradient for Lasso(section4).We describe an SPG algorithm for solving (LSτ).Each iteration of this method requires an orthogonal projection of an n-vector onto the convex set x 1≤τ.In section4.2we give an algorithm for this projection with a worst-case complexity of O(n log n).In many important applications,A is a Fourier-type operator,and matrix-vector products with A and A T can be obtained with O(n log n)cost.The projection cost is typically much smaller than the worst-case, and the dominant cost in our algorithm consists of the matrix-vector products,as it does in other algorithms for basis pursuit denoise.We also show how the projection algorithm can easily be extended to project complex vectors,which allows us to extend the SPG algorithm to problems in the complex domain.Implementation and numerical experiments(sections5and6).To demonstrate the effectiveness of our approach,we run our algorithm on a set of benchmark problems and compare it to other state-of-the-art solvers.In sections6.1and6.2we report numerical results on a series of(BPσ)and(BP)problems,which are normally considered as distinct problems.In section6.3we report numerical results on a series of(LSτ)problems for various values ofτ,and compare against the equivalent(QPλ)formulations.In section6.4we show how to capitalize on the smoothness of the Pareto curve to obtain quick and approximate solutions to(BPσ).1.3.Assumption.We make the following blanket assumption throughout:Assumption1.1.The vector b∈range(A),and b=0.This assumption is only needed in order to simplify the discussion,and it implies that(BPσ)is feasible for all4 E.van den BERG and M.P.FRIEDLANDERσ≥0.In many applications,such as compressed sensing[10–12,23],A has full row rank,and therefore this assumption is satisfied automatically.1.4.Related work.Homotopy approaches.A number of approaches have been suggested for solving (BPσ),many of which are based on repeatedly solving(QPλ)for various values ofλ. Notable examples of this approach are Homotopy[41,42]and Lars[26],which solve (QPλ)for essentially all values ofλ.In this way,they eventuallyfind the value ofλthat recovers a solution of(BPσ).These active-set continuation approaches begin withλ= A T b ∞(for which the corresponding solution is xλ=0),and gradually reduceλin stages that predictably change the sparsity pattern in xλ.The remarkable efficiencyof these continuation methods follows from their ability to systematically update the resulting sequence of solutions.(See Donoho and Tsaig[25]and Friedlander and Saunders[28]for discussions of the performance of these methods.)The computational bottleneck for these methods is the accurate solution at each iteration of a least-squares subproblem that involves a subset of the columns of A.In some applications(such as the seismic image reconstruction problem[32])the size of this subset can become large, and thus the least-squares subproblems can become prohibitively expensive.Moreover, even if the correct valueλσis known a priori,the method must necessarily begin with λ= A T b ∞and traverse all critical values ofλdown toλσ.Basis pursuit denoise as a cone program.The problem(BPσ)withσ>0can be considered as a special case of a generic second-order cone program[8,Ch.5]. Interior-point(IP)algorithms for general conic programs can be very effective if the matrices are available explicitly.Examples of general-purpose software for cone programs include SeDuMi[45]and MOSEK[39].The software package 1-magic[9] contains an IP implementation specially adapted to(BPσ).In general,the efficiency of IP implementations relies ultimately on their ability to efficiently solve certain linear systems that involve highly ill-conditioned matrices.Basis pursuit as a linear program.The special caseσ=0corresponding to(BP) can be reformulated and solved as a linear program.Again,IP methods are known to be effective for general linear programs,but many IP implementations for general linear programming,such as CPLEX[16]and MOSEK,require explicit matrices.The solver PDCO[44],available within the SparseLab package,is capable of using A as an operator only,but it often requires many matrix-vector multiplications to converge, and as we report in section6.2,it is not generally competitive with other approaches. We are not aware of simplex-type implementations that do not require A explicitly.Sampling the Pareto curve.A common approach for obtaining approximate solu-tions to(BPσ)is to sample various points on the Pareto curve;this is often accomplished by solving(QPλ)for a decreasing sequence of values ofλ.As noted by Das and Den-nis[19],and more recently by Leyffer[35],a uniform distribution of weightsλcan lead to an uneven sampling of the Pareto curve.In contrast,by instead parameterizing the Pareto curve byσorτ(via the problem(BPσ)or(LSτ)),it is possible to obtain a more uniform sample of the Pareto curve;see section6.4.Projected gradient.Our application of the SPG algorithm to solve(LSτ)follows Birgin et al.[5]closely for minimization of general nonlinear functions over arbitrary convex sets.The method they propose combines projected-gradient search directions with the spectral step length that was introduced by Barzilai and Borwein[1].A nonmonotone linesearch is used to accept or reject steps.The key ingredient of Birgin et al.’s algorithm is the projection of the gradient direction onto a convex set,which in our case is defined by the constraint in(LSτ).In their recent report,FigueiredoPROBING THE PARETO FRONTIER FOR BASIS PURSUIT SOLUTIONS5Fig.2.1:A typical Pareto curve(solid line)showing two iterations of Newton’s method.The first iteration is available at no cost.et al.[27]describe the remarkable efficiency of an SPG method specialized to(QPλ). Their approach builds on the earlier report by Dai and Fletcher[18]on the efficiency of a specialized SPG method for general bound-constrained quadratic programs(QP s).2.The Pareto curve.The functionφdefined by(1.1)yields the optimal value of the constrained problem(LSτ)for each value of the regularization parameterτ.Its graph traces the optimal trade-offbetween the one-norm of the solution x and the two-norm of the residual r,which defines the Pareto curve.Figure2.1shows the graph ofφfor a typical problem.The Newton-based root-finding procedure that we propose for locating specific points on the Pareto curve—e.g.,finding roots of(1.2)—relies on several important properties of the functionφ.As we show in this section,φis a convex and differentiable function ofτ.The differentiability ofφis perhaps unintuitive,given that the one-norm constraint in(LSτ)is not differentiable.To deal with the nonsmoothness of the one-norm constraint,we appeal to Lagrange duality theory.This approach yields significant insight into the properties of the trade-offcurve.We discuss the most important properties below.2.1.The dual subproblem.The dual of the Lasso problem(LSτ)plays a prominent role in understanding the Pareto curve.In order to derive the dual of(LSτ), wefirst recast(LSτ)as the equivalent problemr 2subject to Ax+r=b, x 1≤τ.(2.1) minimizer,xThe dual of this convex problem is given bymaximizeL(y,λ)subject toλ≥0,(2.2)y,λwhere{ r 2−y T(Ax+r−b)+λ( x 1−τ)}L(y,λ)=infx,ris the Lagrange dual function,and the m-vector y and scalarλare the Lagrange multipliers(e.g.,dual variables)corresponding to the constraints in(2.1).We use the6 E.van den BERG and M.P.FRIEDLANDERseparability of the infimum in r and x to rearrange terms and arrive at the equivalent statementL(y,λ)=b T y−τλ−supr {y T r− r 2}−supx{y T Ax−λ x 1}.We recognize the suprema above as the conjugate functions of r 2and ofλ x 1, respectively.For an arbitrary norm · with dual norm · ∗,the conjugate function of f(x)=α x for anyα≥0is given byf∗(y):=supx {y T x−α x }=0if y ∗≤α,∞otherwise;(2.3)see Boyd and Vandenberghe[8,section3.3.1].With this expression of the conjugate function,it follows that(2.2)remains bounded if and only if the dual variables y and λsatisfy the constraints y 2≤1and A T y ∞≤λ.The dual of(2.1),and hence of (LSτ),is then given bymaximizey,λb T y−τλsubject to y 2≤1, A T y ∞≤λ;(2.4)the nonnegativity constraint onλis implicitly enforced by the second constraint.Importantly,the dual variables y andλcan easily be computed from the optimal primal solutions.To derive y,first note that from(2.3),supr{y T r− r 2}=0if y 2≤1.(2.5)Therefore,y=r/ r 2,and we can without loss of generality take y 2=1in(2.4). To deriveλ,note that as long asτ>0,λmust be at its lower bound,as implied by the constraint A T y ∞≤λ.Hence,we takeλ= A T y ∞.(If r=0orτ=0,the choice of y orλ,respectively,is arbitrary.)The dual variable y can then be eliminated, and we arrive at the following necessary and sufficient optimality conditions for the primal-dual solution(rτ,xτ,λτ)of(2.1):Axτ+rτ=b, xτ 1≤τ(primal feasibility);(2.6a)A T rτ ∞≤λτ rτ 2(dual feasibility);(2.6b)λτ( xτ 1−τ)=0(complementarity).(2.6c)2.2.Convexity and differentiability of the Pareto curve.LetτBPbe the optimal objective value of the problem(BP).This corresponds to the smallest value ofτsuch that(LSτ)has a zero objective value.As we show below,φis nonincreasing,and thereforeτBPis thefirst point at which the graph ofφtouches the horizontal axis.Our assumption that0=b∈range(A)implies that(BP)is feasible,and thatτBP>0. Therefore,at the endpoints of the interval of interest,φ(0)= b 2>0andφ(τBP)=0.(2.7) As the following result confirms,the function is convex and strictly decreasing overthe intervalτ∈[0,τBP].It is also continuously differentiable on the interior of this interval—this is a crucial property.PROBING THE PARETO FRONTIER FOR BASIS PURSUIT SOLUTIONS 7Theorem 2.1.(a)The function φis convex and nonincreasing.(b)For all τ∈(0,τBP ),φis continuously differentiable,φ (τ)=−λτ,and theoptimal dual variable λτ= A T y τ ∞,where y τ=r τ/ r τ 2.(c)For τ∈[0,τBP ], x τ 1=τ,and φis strictly decreasing.Proof .(a)The function φcan be restated asφ(τ)=inf x f (x,τ),(2.8)wheref (x,τ):= Ax −b 2+ψτ(x )and ψτ(x ):= 0if x 1≤τ,∞otherwise .Note that by (2.3),ψτ(x )=sup z {x T z −τ z ∞},which is the pointwise supremum of an affine function in (x,τ).Therefore it is convex in (x,τ).Together with the convexity of Ax −b 2,this implies that f is convex in (x,τ).Consider any nonnegative scalars τ1and τ2,and let x 1and x 2be the corresponding minimizers of (2.8).For any β∈[0,1],φ(βτ1+(1−β)τ2)=inf x f (x,βτ1+(1−β)τ2)≤f βx 1+(1−β)x 2,βτ1+(1−β)τ2≤βf (x 1,τ1)+(1−β)f (x 2,τ2)=βφ(τ1)+(1−β)φ(τ2).Hence,φis convex in τ.Moreover,φis nonincreasing because the feasible set enlarges as τincreases.(b)The function φis differentiable at τif and only if its subgradient at τis unique[43,Theorem 25.1].By [4,Proposition 6.5.8(a)],−λτ∈∂φ(τ).Therefore,to prove differentiability of φit is enough show that λτis unique.Note that λappears linearly in (2.4)with coefficient −τ<0,and thus λτis not optimal unless it is at its lower bound,as implied by the constraint A T y ∞≤λ.Hence,λτ= A T y τ ∞.Moreover,convexity of (LS τ)implies that its optimal value is unique,and so r τ≡b −Ax τis unique.Also, r τ >0because τ<τBP (cf.(2.7)).As discussed in connection with(2.5),we can then take y τ=r τ/ r τ 2,and so uniqueness of r τimplies uniqueness of y τ,and hence uniqueness of λτ,as required.The continuity of the gradient follows from the convexity of φ.(c)The assertion holds trivially for τ=0.For τ=τBP , x τBP 1=τBP by definition.It only remains to prove part (c)on the interior of the interval.Note that φ(τ)≡ r τ >0for all τ∈[0,τBP ).Then by part (b),λτ>0,and hence φis strictly decreasing for τ<τBP .But because x τand λτboth satisfy the complementarity condition in (2.6),it must hold that x τ 1=τ.2.3.Generic regularization.The technique used to prove Theorem 2.1does not in any way rely on the specific norms used in the objective and regularization functions,and it can be used to prove similar properties for the generic regularized fitting problemminimize x Ax −b s subject to x p ≤τ,(2.9)8 E.van den BERG and M.P.FRIEDLANDERwhere1≤(p,s)≤∞define the norms of interest,i.e., x p=(i|x i|p)1/p.Moregenerally,the constraint here may appear as Lx p,where L may be rectangular. Such a constraint defines a seminorm,and it often arises in discrete approximations of derivative operators.In particular,least-squares with Tikhonov regularization[47], which corresponds to p=s=2,is used extensively for the regularization of ill-posed problems;see Hansen[29]for a comprehensive study.In this case,the Pareto curve defined by the optimal trade-offbetween x 2and Ax−b 2is often called the L-curve because of its shape when plotted on a log-log scale[30].If we define¯p and¯s such that1/p+1/¯p=1and1/s+1/¯s=1,then the dual of the generic regularization problem is given bymaximizey,λb T y−τλsubject to y ¯s≤1, A T y ¯p≤λ.As with(2.4),the optimal dual variables are given by y=r/ r ¯p andλ= A T y ¯s. This is a generalization of the results obtained by Dax[21],who derives the dual for p and s strictly between1and∞.The corollary below,which follows from a straightfoward modification of Theorem2.1,asserts that the Pareto curve defined for any1≤(p,s)≤∞in(2.7)has the properties of convexity and differentiability.Corollary2.2.Letθ(τ):= rτ s,where rτ:=b−Axτ,and xτis the optimal solution of(2.9).(a)The functionθis convex and nonincreasing.(b)For allτ∈(0,τBP),θis continuously differentiable,θ (τ)=−λτ,and the optimal dual variableλτ= A T yτ ¯p,where yτ=rτ/ rτ ¯s.(c)Forτ∈[0,τBP], xτ p=τ,andθis strictly decreasing.3.Rootfinding.As we briefly outlined in section1.2,our algorithm generatesa sequence of regularization parametersτk→τσbased on the Newton iterationτk+1=τk+∆τk with∆τk:=σ−φ(τk)/φ (τk),(3.1)such that the corresponding solutions xτk of(LSτk)converge to xσ.For valuesofσ∈(0, b 2),Theorem2.1implies thatφis convex,strictly decreasing,and continuously differentiable.In that case it is clear thatτk→τσsuperlinearly for all initial valuesτ0∈(0,τBP)(see,e.g.,Bertsekas[3,proposition1.4.1]).The efficiency of our method,as with many Newton-type methods for large problems,ultimately relies on the ability to carry out the iteration described by(3.1) with only an approximation ofφ(τk)andφ (τk).Although the nonlinear equation(1.2) that we wish to solve involves only a single variableτ,the evaluation ofφ(τ)involves the solution of(LSτ),which can be a large optimization problem that is expensive to solve to full accuracy.For systems of nonlinear equations in general,inexact Newton methods assume that the Newton system analogous to the equationφ (τk)∆τk=σ−φ(τk)is solved only approximately,with a residual that is a fraction of the right-hand side.A constant fraction yields a linear convergence rate,and a fraction tending to zero yields aPROBING THE PARETO FRONTIER FOR BASIS PURSUIT SOLUTIONS9 superlinear convergence rate(see,e.g.,Nocedal and Wright[40,theorem7.2]).However, the inexact-Newton analysis does not apply to the case where the right-hand side(i.e., the function itself)is known only approximately,and it is therefore not possible to know a priori the accuracy required to achieve an inexact-Newton-type convergence rate.This is the situation that we are faced with if(LSτ)is solved approximately.As we show below,with only approximate knowledge of the function valueφthis inexact version of Newton’s method still converges,although the convergence rate is sublinear. The rate can be made arbitrarily close to superlinear by increasing the accuracy with which we computeφ.3.1.Approximate primal-dual solutions.In this section we use the duality gap to derive an easily computable expression that bounds the accuracy of the computed function value ofφ.The algorithm for solving(LSτ)that we outline in section4maintains feasibility of the iterates at all iterations.Thus,an approximate solution¯xτand its corresponding residual¯rτ:=b−A¯xτsatisfy¯xτ 1≤τ,and ¯rτ 2≥ rτ 2>0,(3.2).We where the second set of inequalities holds because¯xτis suboptimal andτ<τBPcan thus construct the approximations¯yτ:=¯rτ/ ¯rτ 2and¯λτ:= A T¯yτ ∞to the dual variables that are dual feasible,i.e.,they satisfy(2.6b).The value of the dual problem(2.2)at any feasible point gives a lower bound on the optimal value rτ 2,and the value of the primal problem(2.1)at any feasible point gives an upper bound on the optimal value.Therefore,b T¯yτ−τ¯λτ≤ rτ 2≤ ¯rτ 2.(3.3) We use the duality gapδτ:= ¯rτ 2−(b T¯yτ−τ¯λτ)(3.4) to measure the quality of an approximate solution¯xτ.By(3.3),δτis necessarily nonnegative.Let¯φ(τ):= ¯rτ 2be the objective value of(LSτ)at the approximate solution¯xτ. The duality gap at¯xτprovides a bound on the difference betweenφ(τ)and¯φ(τ).If we additionally assume that A is full rank(so that its condition number is bounded), we can also useδτto provide a bound on the difference between the derivativesφ (τ)),and¯φ (τ).From(3.3)–(3.4)and from Theorem2.1(b),for allτ∈(0,τBP¯φ(τ)−φ(τ)<δand|¯φ (τ)−φ (τ)|<γδτ(3.5)τfor some positive constantγthat is independent ofτ.It follows from the definition ofφ and from standard properties of matrix norms thatγis proportional to the condition number of A.3.2.Local convergence rate.The following theorem establishes the local con-vergence rate of an inexact Newton method for(1.2)whereφandφ are known only approximately.10 E.van den BERG and M.P.FRIEDLANDER Theorem 3.1.Suppose that A has full rank,σ∈(0, b 2),and δk :=δτk →0.Then if τ0is close enough to τσ,the iteration (3.1)—with φand φ replaced by ¯φand ¯φ—generates a sequence τk →τσthat satisfies |τk +1−τσ|=γδk +ηk |τk −τσ|,(3.6)where ηk →0and γis a positive constant.Proof .Because φ(τσ)=σ∈(0, b 2),equation (2.7)implies that τσ∈(0,τBP ).By Theorem 2.1we have that φ(τ)is continuously differentiable for all τclose enough to τσ,and so by Taylor’s theorem,φ(τk )−σ= 10φ (τσ+α[τk −τσ])dα·(τk −τσ)=φ (τk )(τk −τσ)+10 φ (τσ+α[τk −τσ])−φ (τk ) ·dα(τk −τσ)=φ (τk )(τk −τσ)+ω(τk ,τσ),where the remainder ωsatisfiesω(τk ,τσ)/|τk −τσ|→0as |τk −τσ|→0.(3.7)By (3.5)and because (3.2)holds for τ=τk ,there exist positive constants γ1and γ2,independent of τk ,such that φ(τk )−σφ (τk )−¯φ(τk )−σ¯φ (τk ) ≤γ1δk and |φ (τk )|−1<γ2.Then,because ∆τk = σ−¯φ(τk ) /¯φ (τk ),|τk +1−τσ|=|τk −τσ+∆τk |= −¯φ(τk )−σ¯φ (τk )+1φ (τk ) φ(τk )−σ−ω(τk ,τσ) ≤ φ(τk )−σφ (τk )−¯φ(τk )−σ¯φ (τk ) + ω(τk ,τσ)φ (τk )=γ1δk +γ2|ω(τk ,τσ)|=γ1δk +ηk |τk −τσ|,where ηk :=γ2|ω(τk ,τσ)|/|τk −τσ|.With τk sufficiently close to τσ,(3.7)implies that ηk <1.Applythe above inequality recursively ≥1times to obtain|τk + −τσ|≤γ1 i =1(γ1) −i δk +i −1+(ηk ) |τk −τσ|,and because δk →0and ηk <1,it follows that τk + →τσas →∞.Thus τk →τσ,as required.By again applying (3.7),we have that ηk →0.Note that if (LS τ)is solved exactly at each iteration,such that δk =0,then Theorem 3.1shows that the convergence rate is superlinear,as we expect of a standard Newton iteration.In effect,the convergence rate of the algorithm depends on the rate at which δk →0.If A is rank deficient,then the constant γin (3.6)is infinite;we thus expect that ill-conditioning in A leads to slow convergence unless δk =0,i.e.,φis evaluated accurately at every iteration.。

Pocket Surf III I ABSoLute MoBILItY 二维表面粗糙度测量仪说明书

Pocket Surf III I ABSoLute MoBILItY 二维表面粗糙度测量仪说明书
1 Rz = 5 (Rz1 + Rz2 + ... + Rz5) The Rz definition is identical to the definition in DIN4768: 1990. The ten point height Rz as well as the parameter symbol Ry of ISO 4287:1984 has been canceled. Maximum roughness depth Rmax is the the largest single roughness depth with the evalution length. (DIN EN ISO 4288; Rmax is also called Rz1max ).
• Selectable traverse length 1, 3 or 5 cut-offs of 0.8 mm/.030“
• Operates in any position – horizontal, vertical, and upside down
• Four switchable probe positions – axial (folded) or at 90°, 180° or 270°
Ra Mean roughness DIN EN ISO 4287, ASME B46.1
Roughness average Ra is the arithmetic average of the absolute values of the roughness profile ordinates.
Meets ASME-B46.1, ISO, DIN standards and MIL specifications

圆极化与线计划的设置

圆极化与线计划的设置

圆极化与线计划的设置英文回答:Circular polarization and linear polarization are two different ways of describing the orientation of electromagnetic waves. In circular polarization, the electric field vector of the wave rotates in a circular pattern as the wave propagates. This rotation can be clockwise or counterclockwise. In linear polarization, the electric field vector of the wave oscillates in a straight line.Circular polarization can be achieved by combining two perpendicular linearly polarized waves with a phase difference of 90 degrees. This can be done using a device called a quarter-wave plate or a combination of a half-wave plate and a linear polarizer. The resulting wave will have a rotating electric field vector.Linear polarization can be achieved by using a linearpolarizer, which is a device that only allows waves with a specific orientation of the electric field vector to pass through. The polarizer absorbs or blocks waves with orientations perpendicular to the desired polarization.The choice between circular polarization and linear polarization depends on the specific application. Circular polarization is often used in satellite communication to minimize signal loss due to the rotation of the satellite and to improve signal reception in areas with high levelsof interference. It is also used in some optical systems to eliminate the effects of birefringence. Linear polarization, on the other hand, is commonly used in antennas, optical filters, and polarization-sensitive detectors.中文回答:圆极化和线极化是描述电磁波方向的两种不同方式。

花样滑冰专业术语

花样滑冰专业术语

花样滑冰专业术语花样滑冰 figure skating 花样滑冰鞋 figure skate 冰场 ice arena; rink 起滑脚 starting foot 冰刀套 skate guard 花样 pattern起跳 take off奥运会花样滑冰四个项目男子单人滑 men女子单人滑 ladies双人滑 pairs冰舞 ice dancing(包括:规定舞 compulsory dance创编舞 original dance 和自由舞 freedance)男子单人滑和女子单人滑两个项目自由滑 free skating 短节目 short program花样滑冰技术动作特种圆形 advanced figure 半周、半圆 half-circle 旋转 spin燕式旋转 arabesque spin 倒滑;退滑 back skating 直立滑行 ride冰上表演 ice show艺术印象 artistic impression附加动作 additional move 前一周半跳 axel-paulsen(源自阿克塞尔?保尔森,19世纪挪威花样溜冰选手)旋转轴 axle of revolution 倒滑压步 back cross over 后内刃 back in 后内括弧形 back in bracket 后内变刃形 back in change 后内圆形 back in circle 后内外勾形 back in counter 后内结环形 back in loop后内勾形 back in rocker 后内环绕,后内螺旋形 back in spiral 后辅刃back out后内括弧形 back out bracket 后外变刃形 back out change 后外圆形 back out circle 后外勾形 back out counterback out loop 后外结环形后外环绕,后外螺旋形 back out spiral 环绕图形 spiral figure单脚蹲转 Haines spin单脚直立旋转 one-foot upright spin单脚括弧形 paragraph bracket 单脚结环形 paragraph loop 单个动作individual part 单人旋转 solo spin抱身 hand-to-body grip 连接动作 connecting move 连接步 connecting step 步法 footwork技术水平 technical merit 侧翻举 cartwheel lift变刃 change of edge借跳 partner-assisted jump 旋转动作 spinning movement 滑脚 employed foot;gliding foot;tracting foot滑步 glide滑弧线 curved stroke滑区 skating area编排 combination硬转身 forced turn集体动作 group move外刃 outside;outer edge内刃 inner edge;inside右脚外刃 right outside右脚内刃 right inside用刀刃蹬冰 stroke from the edge 用刀齿蹬冰 stroke from the point of the skate握臂 hand-to-arm grip 携手 hand-to-hand grip 稳健 sureness撑竿式跳 pole vaule jump 横一字 spread eagle镜子滑 mirror skating刀齿旋转 pirouette;toe-scratch 刀齿向下 toe of the skate pointing downward刀齿蹬冰 toe push双人旋转 pair spin分腿 split分腿举 split lift平刃旋转 flat-foot spin 规定图形滑 skating of prescribed movements 同脚变刃的前进跳跃 forward jump 前刃变后刃半周跳 Mohawk评分 assign marks人工冷冻冰场artificially frozen rink “6”字形 figure six“8”字形 curve eight浮脚 free foot双人动作 pair move转体 turn转体半周 half turn转体两周 double turn圆形滑 concentric stroking 同足后内结环一周跳 one-foot Salchow jump 异足后外结环一周跳 half loop jump 后内点冰一周跳 flip jump;toe Salchow 后内点冰“3”字跳 half toe Salchow。

胸部评估.

胸部评估.
Kusmols breath
Sighing respiration
Tidal breath Ataxic breath Inhibitory breath
Palpation

Thoracic expansion
Vocal
fremitus(触觉语颤) Pleural friction fremitus(胸膜摩擦感) Confirm the inspection

Vertical lines
Anterior middle line(前正中线) Mid-clavicular lines(锁骨中线) Spinal line(后正中线) axillary lines (anterior,middle, posterior) (腋前、中、后线) Scapular lines(肩胛线)
叩诊的位置
叩诊肺前界
叩诊肺下界(肩胛下线);
肩胛下线叩诊肺下界移动度:
(三)percussion
1.叩诊方法 direct percussion indirect percussion:the palmar surface of the left distal phalanx of the middle finger serves as the pleximeter and is firmly placed on the chest wall in an interspace;parallel to the ribs.
3.胸部叩诊音的分类
The normal percussion note varies with the thick-ness of the chest wall and the force applied by the examiner. 1 The clear, long, low-pitched sound elicited over the normal lung is termed resonance.

微积分术语中英文对照

微积分术语中英文对照

微积分术语中英文对照A、B:Absolute convergence :绝对收敛Absolute extreme values :绝对极值Absolutemaximum and minimum :绝对极大与极小Absolute value:绝对值Absolute value function :绝对值函数Acceleration:加速度Antiderivative :原函数,反导数Approximate integration:近似积分(法)Approximation :逼近法bydifferentials:用微分逼近linear :线性逼近法by Simpson’s Rule :Simpson法则逼近法by the Trapezoidal Rule :梯形法则逼近法Arbitrary constant :任意常数Arc length :弧长Area :面积under a curve :曲线下方之面积between curves :曲线间之面积in polar coordinates:极坐标表示之面积ofasector of a circle:扇形之面积of a surface of a revolution :旋转曲面之面积Asymptote :渐近线horizontal :水平渐近线slant:斜渐近线vertical :垂直渐近线Averagespeed :平均速率Average velocity :平均速度Axes,coordinate:坐标轴Axesofellipse :椭圆之对称轴Binomial series:二项式级数Binomial theorem:二项式定理C:Calculus :微积分differential:微分学integral:积分学Cartesian coordinates:笛卡儿坐标一般指直角坐标Cartesian coordinatessystem:笛卡儿坐标系Cauch’s Mean ValueTheorem:柯西中值定理Chain Rule:链式法则Circle :圆Circular cylinder :圆柱体,圆筒Closedinterval:闭区间Coefficient :系数Composition of function :复合函数Compound interest :复利Concavity:凹性Conchoid:蚌线Conditionally convergent:条件收敛Cone:圆锥Constant function :常数函数Constant of integration:积分常数Continuity :连续性atapoint :在一点处之连续性of afunction:函数之连续性on an interval:在区间之连续性from the left:左连续from the right:右连续Continuous function :连续函数Convergence:收敛interval of:收敛区间radius of :收敛半径Convergent sequence :收敛数列series :收敛级数Coordinates:坐标Cartesian :笛卡儿坐标cylindrical:柱面坐标polar :极坐标rectangular :直角坐标spherical:球面坐标Coordinateaxes:坐标轴Coordinate planes :坐标平面Cosine function:余弦函数Criticalpoint:临界点Cubicfunction:三次函数Curve :曲线Cylinder:圆筒,圆柱体,柱面Cylindrical Coordinates :圆柱坐标D:Decreasing function :递减函数Decreasing sequence :递减数列Definiteintegral:定积分Degree of a polynomial:多项式之次数Density :密度Derivative :导数of a composite function:复合函数之导数of a constant function :常数函数之导数directional :方向导数domainof:导数之定义域of exponential function :指数函数之导数higher :高阶导数partial:偏导数of a powerfunction :幂函数之导数of apowerseries:羃级数之导数of a product:积之导数of a quotient :商之导数asa rate of change :导数当作变化率right—hand :右导数second:二阶导数as theslope of a tangent:导数看成切线之斜率Determinant :行列式Differentiablefunction :可导函数Differential :微分Differential equation :微分方程partial:偏微分方程Differentiation:求导法implicit :隐求导法partial :偏微分法term by term :逐项求导法Directional derivatives :方向导数Discontinuity :不连续性Disk method:圆盘法Distance:距离Divergence:发散Domain:定义域Dotproduct :点积Double integral:二重积分change of variable in:二重积分之变数变换in polar coordinates :极坐标二重积分E、F、G:Ellipse :椭圆Ellipsoid :椭圆体Epicycloid :外摆线Equation:方程式Even function:偶函数Expected Valued :期望值Exponential Function:指数函数Exponents , laws of:指数率Extremevalue:极值ExtremeValueTheorem :极值定理Factorial :阶乘First DerivativeTest :一阶导数试验法First octant :第一卦限Focus :焦点Fractions:分式Function :函数FundamentalTheorem ofCalculus:微积分基本定理Geometricseries:几何级数Gradient:梯度Graph :图形Green Formula :格林公式H:Half-angle formulas :半角公式Harmonic series:调和级数Helix:螺旋线Higher Derivative:高阶导数Highermathematics高等数学Horizontal asymptote:水平渐近线Horizontalline :水平线Hyperbola:双曲线Hyperboloid :双曲面I:Implicit differentiation:隐求导法Implicit function :隐函数Improper integral :反常积分, 广义积分Increasing,Decreasing Test :递增或递减试验法Increment :增量Increasing Function :增函数Indefinite integral :不定积分Independentvariable:自变量Indeterminatefrom :不定型Inequality :不等式Infinite point:无穷极限点Infinite series :无穷级数Inflectionpoint :反曲点Instantaneous velocity:瞬时速度Integer:整数Integral :积分Integrand :被积函数Integration:积分Integration by part:分部积分法Intercepts :截距Intermediatevalue of Theorem:中值定理Interval:区间Inversefunction :反函数Inversetrigonometricfunction :反三角函数Iterated integral :逐次积分L:Laplace transform :Laplace 变换Law ofsines:正弦定理Law ofCosines :余弦定理Leastupper bound :最小上界Left-handderivative:左导数Left—handlimit:左极限Lemniscate :双钮线Length :长度Levelcurve :等高线L'Hospital's rule: 洛必达法则Limacon :蚶线Limit :极限Linear approximation:线性近似Linear equation :线性方程式Linear function:线性函数Linearity :线性Linearization:线性化Line in theplane :平面上之直线Line in space:空间之直线Local extreme :局部极值Local maximumand minimum :局部极大值与极小值Logarithm:对数Logarithmic function:对数函数M、N、O:Maximumandminimum values :极大与极小值Mean Value Theorem:均值定理Multipleintegrals :重积分Multiplier :乘子Natural exponential function:自然指数函数Natural logarithmfunction:自然对数函数Natural number:自然数Normal line:法线Normal vector:法向量Number:数Octant :卦限Oddfunction:奇函数One-sided limit:单边极限Open interval :开区间Optimization problems :最佳化问题Order :阶Ordinary differential equation:常微分方程Origin :原点Orthogonal :正交的P、Q:Parabola:拋物线Parabolic cylinder:抛物柱面Paraboloid:抛物面Parallelepiped:平行六面体Parallel lines :平行线Parameter :参数Partial derivative:偏导数Partial differential equation:偏微分方程Partialfractions:部分分式Partial integration:部分积分Partition :分割Period:周期Periodic function:周期函数Perpendicular lines :垂直线Piecewise defined function:分段定义函数Plane:平面Point of inflection:反曲点Polar axis :极轴Polar coordinate:极坐标Polarequation:极方程式Pole :极点Polynomial:多项式Positive angle :正角Point—slopeform:点斜式Power function:幂函数Product:积Quadrant :象限QuotientLaw of limit :极限的商定律Quotient Rule:商定律R:Radius of convergence :收敛半径Rangeof a function :函数的值域Rate of change :变化率Rational function:有理函数Rationalizing substitution :有理代换法Rational number:有理数Real number:实数Rectangular coordinates:直角坐标Rectangular coordinate system :直角坐标系Relativemaximum and minimum :相对极大值与极小值Revenue function :收入函数Revolution ,solid of :旋转体Revolution,surface of:旋转曲面Riemann Sum:黎曼和Right-handderivative:右导数Right—hand limit:右极限Root :根S:Saddlepoint :鞍点Scalar:纯量Secant line:割线Secondderivative :二阶导数SecondDerivative Test :二阶导数试验法Second partialderivative:二阶偏导数Sector :扇形Sequence:数列Series:级数Set :集合Shell method:剥壳法Sinefunction:正弦函数Singularity :奇点Slant,Obliqueasymptote :斜渐近线Slope:斜率Slope-intercept equationof a line :直线的斜截式Smooth curve:平滑曲线Smoothsurface :平滑曲面Solid of revolution :旋转体Space :空间Speed :速率Spherical coordinates :球面坐标Squeeze Theorem:夹挤定理Stepfunction :阶梯函数Strictlydecreasing:严格递减Strictly increasing :严格递增Substitutionrule: 替代法则Sum:和Surface :曲面Surface integral :面积分Surface of revolution :旋转曲面Symmetry :对称T:Tangent function:正切函数Tangentline:切线Tangentplane :切平面Tangentvector :切向量Taylor'sformula :泰勒公式Totaldifferential :全微分Trigonometric function:三角函数Trigonometric integrals :三角积分Trigonometric substitutions:三角代换法Tripe integrals :三重积分V、X、Z:Valueof function:函数值Variable :变量Vector :向量Velocity:速度Verticalasymptote :垂直渐近线Volume :体积X—axis:x轴X —coordinate :x坐标X -intercept :x截距Zerovector:函数的零点Zeros of a polynomial :多项式的零点。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

Received date/Accepted date Abstract Studies of the Spectral Energy Distribution of Young Stellar Objects suggest that the outer disk of FU Orionis objects might be self-gravitating. In this paper we propose a method to test directly whether, in these objects, significant deviations from Keplerian rotation occur. In a first approach, we have used a simplified model of the disk vertical structure that allows us to quickly bring out effects related to the disk self-gravity. We find that the often studied optical and near-infrared line profiles are produced too close to the central object to provide significant evidence for non-Keplerian rotation. Based on parameters relevant for the case of FU Ori, we show that high-resolution long-wavelength spectroscopy, of the far-infrared H2 pure rotational lines (sometimes observed in “passive” protostellar disks) and sub-mm CO lines, should be well suited to probe the rotation curve in the outer disk, thus measuring to what extent it is affected by the disk self-gravity. The results of the present exploratory paper should be extended soon to a more realistic treatment of the disk vertical structure. Key words. accretion, accretion disks – gravitation – stars: pre-main sequence
1. Introduction
Very young (Class 0 or Class I) protostellar sources are thought to be characterized by a fairly high mass accretion rate (Calvet et al., 2000), but they appear to be still deeply embedded in their protostellar envelopes, so that their disk remains hidden to direct observations. On the other hand, in the more evolved T Tauri stars the mass accretion rate is generally considered to be modest, so that their disks are heated by the combined contribution of viscous dissipation and irradiation from the central star, making it more difficult to extract detailed information about the accretion process. In this context, FU Orionis objects (Hartmann & Kenyon, 1996) are a rather small but remarkable class of pre-main sequence stars, because they are ideal “laboratories” to test the process of disk accretion during the early stages of star formation. In fact, the disk of FU Orionis objects, differently from that of Class I objects, can be studied directly from its optical emission, and, differently from that of T Tauri stars, is likely to be the site of “active” accretion. The distinctive feature of FU Orionis objects is that they undergo a violent outburst phase. During the outburst, they can increase their luminosity by more than 4 magnitudes in a period of a few years, and then slowly return to a quiescent state on a much longer timescale. It is commonly believed that such outbursts are the result of a significant increase of the mass accretion rate in the disk
which usually surrounds Young Stellar Objects. Many possible mechanisms to trigger the outburst phase have been discussed in the literature, among which are a tidal interaction with a companion star (Bonnell & Bastien, 1992), the onset of a thermal instability in a partially ionized disk, when the outer disk is already accreting at a sufficiently high rate (Bell et al., 1995), or the onset of thermal instability induced by the presence of a satellite within the disk (Clarke & Syer, 1996). The high accretion rate makes the emission of these systems dominated by the accretion luminosity, with only a minor contribution from the central star. The modeling of the Spectral Energy Distribution (SED) of FU Orionis objects provides an estimate of the ˙ (as discussed below in Sect. 2.1) which product M⋆ M 2 turns out to be of the order of 10−4 M⊙ /yr. For a stellar mass M⋆ ≈ M⊙ , this would correspond to an accretion rate as high as ≈ 10−4 M⊙ /yr. “Standard” models based on the presence of a Keplerian accretion disk lead to a satisfactory fit of the available photometric data only for wavelengths shorter than 10µm; the luminosity at longer wavelengths is usually much larger than expected. A flaring disk model, in which the outer disk is illuminated by the inner disk, could in principle explain the far-infrared excess, but the required amount of flaring turns out to be often too large (Kenyon & Hartmann, 1991); therefore, the long-wavelength part of the SED is generally at-
2
G. Lodato and G. Bertin: Probing the rotation curve of FU Orionis obnfalling envelope, which is heated by the accretion disk luminosity. Recently we have shown the v
相关文档
最新文档