Cluster-assisted accretion for massive stars
吸积模型应用到活动星系核所遇到的困难-星际物理与化学
Additional heating to the corona is often assumed. But the mechanism is unclear!
• Temperarure inpedendent on M Ti~1011-1012K，近维里温度 Te~109K • Low L: Eddington ratio≾0.01
Narayan et al. 1998
Comparison of the Spectra
log fν soft
hard log hν
Accretion Energy Radiation Low T, high L, geometrically thin disk ≫ 1: multi-color blackbody Soft state spectra
•
Advection-Dominated Accretion Flow (ADAF)
Spectrum from Standard Disk
Multi-colour Blackbody
spectra (optically thick)
kT ~1 keV ( BH binaries)
Log n f(n)
kT ~ 10 eV (AGN)
Log n
Spectrum from ADAF
活动星系核（AGN)
• 活动星系核：星系的核心有剧烈电磁辐射 • 提供活动星系核巨大能源 的发动机—黑洞吸积过程 星系 • 低吸积率的黑洞：正常星系 • 高吸积率的黑洞：活动星系核
活动星系核 类星体 黑洞
吸积的基本模式
•
The standard accretion disk (Shakura & Sunyaev 1973)
COSMIC掩星L2波段截断高度的分布规律统计
COSMIC掩星L2波段截断高度的分布规律统计文章统计了COSMIC掩星数据L2波段截断高度的分布规律，得到了其截断高度的分布特点，为以后L2波段数据的误差分析提供了参考。
标签：无线电掩星；COSMIC；L2波段；截断高度引言在SAC-C和CHAMP无线电掩星计划中，LEO掩星接收机采用的是闭环（Close-Loop，CL）跟踪技术。
闭环技术容易在10km以下产生追踪误差。
而且在低对流层，闭环技术不能跟踪上升掩星[1]。
因此，在2001年，Sokolovskiy 提出了开环追踪技术（Open-Loop），能够有效地克服追踪误差的问题[2]。
开环技术能够跟踪信号到较低的高度，而且不容易发生失锁现象。
但是开环技术只能够消除L1的追踪误差的问题，在低对流层，大部分掩星L2波段信号质量都不好，不能用于电离层的校准，低于一定高度L2波段的信号不可用[3]，因此我们需要对L2信号进行截断，文章主要分析L2信号截断高度的分布规律，了解不同纬度中L2信号截断高度的分布特点，对以后做误差分析提供一定程度的参考。
1 数据介绍为了分析L2截断高度的分布规律，选取了COSMIC网站中COSMIC掩星计划的atmprf干廓线数据，干廓线数据中包含了“Znid”属性，即为L2的截断高度，时间选用的是2011年1、4、7、10四个月份的前五天，选取了全球作为研究区域，研究L2的截断高度随纬度的变化。
表1是对所选区域掩星数量（atmprf 文件）的统计。
我们选取的廓线总计24108条，可以看出在中纬度地区即-60- -30N，30-60N 区域内掩星数量要多于热带地区和两极地区，两极地区最少。
2 分析与讨论我们将全球区域分为6个部分，0-30S、30-60S、60-90S、0-30N、30-60N、60-90N。
分别统计在不同区域内L2截断高度的分布规律，并计算了每个部分L2截断高度的平均值，如图1所示。
在图1中，上面三幅图是北半球L2截断高度的分布图，可以看出在北半球0-15N范围内L2的截断高度主要分布在18-20km，L2质量不好，主要因为在热带低对流层，信号在穿过大气时会产生超折射效应以及多路径效应[4]，影响了L2信号的传播，所以L2质量不好，截斷高度较高。
天文学名词英汉对照表
天文学名词英汉对照表【推荐][名词委审定]汉英天文学名词(定义版, 1998)Ⅰ1054超新星||supernova of 1054(CM T au); 公元1054年出现在金牛座ζ星附近的超新星。
160分钟振荡||160-minute oscillation; 在太阳光球上观测到的周期约为 160分钟的速度起伏。
5分钟振荡||five-minute oscillation; 在太阳光球上观测到的周期约为 5分钟的速度起伏。
Ap星||Ap star; 光谱型B5—F5的特殊主序星，具有异常强的和可变的锰、硅、铕、铬、锶谱线和较强且变化的磁场。
与金属线星类似，自转速度较小，但多半是单星。
BN天体||BN object; 美国天文学家贝克林（Becklin ）和诺伊格鲍尔（Neugebauer）在猎户星云中发现的一个点状红外源。
被认为是恒星刚形成阶段的候选者。
Be星||Be star; 光谱中出现(或曾出现)氢的巴耳末发射线，光度级为Ⅱ—Ⅴ，主要为Ｂ型的恒星。
B型星||B star; 光谱型为B的恒星。
光谱主要特征为中性氦吸收线和氢吸收线。
CCD 摄谱仪||CCD spectrograph; 用CCD(电荷耦合器件)作为辐射接收器的摄谱仪。
CCD 照相机||CCD camera; 用CCD(电荷耦合器件)作为辐射接收器的照相装置。
CCD测光||CCD photometry; 利用CCD 进行的二维光度测量。
CCD天文学||CCD astronomy; 用CCD (电荷耦合器件)作为辐射接收器和探测器的实测天文学。
CH星||CH star; 光谱中CH分子的G带(430.3nm)特强，422.6nm中性钙线与CN分子带偏弱的G3 —K4型星族Ⅱ巨星。
E冕||E corona; 由日冕自身的高次电离原子辐射形成的日冕成分。
F 冕||F corona; 由行星际尘埃云散射太阳光球辐射形成的日冕成分。
G型星||G star; 光谱型为G的恒星。
AR5060——具有复杂结构的大黑子但没有大耀斑发生的太阳活动区(英文)
AR5060——具有复杂结构的大黑子但没有大耀斑发生的太阳活动区(英文)李治凯;栾蒂;李志【期刊名称】《天文研究与技术》【年(卷),期】1989()S1【摘要】AR5060是No.Ⅳ联测期中的第二个目标活动区。
它从1988年6月25日东边缘初现到7月8日转出西边缘消失的14天中,黑子群一直保持最复杂的FKC、EKC型和最复杂的BGD磁型。
6月29日黑子群面积发展到3000面积单位,是第22周以来第一群最大的黑子(更大的是1989年3月的AR5395,面积达3600单位)。
该活动区的黑子群发生过强烈的运动和磁性重联。
似乎具备发生强烈大耀斑的位形特征和动力学条件,可是在这期间,全球耀斑监测所观测到的120多个耀斑(据SGD)中,亚耀斑占81％,1级耀斑占15％,2级耀斑只有3个占4％,而且这3个2级耀斑的X射线级别只达到M6.5,M9.2,M3.9,没有一个达到X级。
在AR5060活动区耀斑活动高峰期的6月28日,29日,30日和7月1日这四天中,云南天文台26CM太阳望远镜观测到其中一个2B/M6.5耀斑(1988年6月29日0737UT)、几个1级耀斑和其它许多亚耀斑。
从黑子群和色球单色照片上作耀斑发生点同黑子相对位置的比较,结果是出乎意料的,在结构复杂、运动剧烈的黑子群内部发生的都是小耀斑,而3个2B/M级耀斑都发生在黑子群以外只有卫星黑子浮现和消失的时期和地点。
【总页数】7页(P74-80)【关键词】AR5060;太阳活动区;面积单位;重联;云南天文台;动力学条件;太阳望远镜;发生点;位形;相对位置【作者】李治凯;栾蒂;李志【作者单位】中国科学院云南天文台【正文语种】中文【中图分类】P1【相关文献】1.太阳活动区卫星黑子与高能耀斑 [J], 洪琴芳2.2.9AR5060——具有复杂结构的大黑子但没有大耀斑发生的太阳活动区 [J], 李治凯3.太阳活动的节律及大耀斑期的预测(英文) [J], 罗葆荣4.AR5060活动区的黑子群及其耀斑活动(英文) [J], 张桂清;包克寒;薛义生;段书林5.AR5060活动区1988年7月2日的耀斑及其磁场特征(英文) [J], 李京因版权原因，仅展示原文概要，查看原文内容请购买。
被动型星载氢钟H型气体电离模型及参数优化
基 于平 面盘绕馈 能 天线的 Ｈ 型 电 离 系统 ，并建 立该 电 离方 法的数 学模 型 。首 先 ，由 Ｍａ ｘ ｗｅ ｌ ｌ 方程 出发 ，推导 出 了电离泡 中的 电场 和磁 场 的分布模 型 ，得到 气体 电 离机 制 ；再 由电磁 场 的
分布方程推 导 出氢气击 穿判 据表 达 式 ，建 立 了气体 电离条件 的数 学模 型 ；最后 ，利 用 电 离击 穿条件判 据方程优化 了电离装 置的设计参 数 。数 值仿 真 与测 试表 明该模 型 是有 效 的 ，采 用优 化 的设计 参数 可以使 电离功耗 和 电磁 干扰 特性得 到改善 ，且 电离直流 电源 总功 率可 小于 ２ ｗ。 关键 词 气体 电 离模 型 电磁 场分 布 电离参数 优 化 被 动 型 星载 氢钟
ＤＯＩ ：１ ０ ． ３ ７ ８ ０ ／ ｊ ． ｉ ｓ ｓ ｎ ． １ ０ ０ ０ — ７ ５ ８ Ｘ． ２ ０ １ ３ ． ０ ５ ． ０ ０ ９
１ 引 言
被 动 型氢 原子 钟是 实现 氢 原子钟 小 型化 的主要 途径 ，它体积 、质量 相对 较小 ，且 相对 于其 他小 型原 子钟 具有 无可 比拟 的中 、长期稳 定 度 ，适 用 于 卫 星定 位 、导 航 、空 间探 测及 天 文 观 测等 系 统 ， 可 以提供超 稳 定度 的标 准频 率信 号 ｌ 】 ］ 。氢 气 电离装 置是 被 动型 星载 氢原 子钟 的重要 物 理部件 ，它将 经 过 提纯器 过 滤后 的纯 净 的氢气 电离 成 为氢原 子 ，以便 经 过选 态 以后产 生高 能态 氢原 子并 注入 到储 存 泡 中 ，从 而产生 １ ． ４ ２ ０ ＧＨｚ 跃 迁微 波信 号 。 电离 电 源的功 率不 仅影 响 了小 型氢 钟 的星载 质量 和体 积 ，还决 定 了 电离 电磁 能对 整个 电子 系统 的 电磁干 扰特 性 （ 电离 振荡 器功 率越 大 ， 电磁 辐射 干扰 越强 ） ｌ ２ ］ 。因此 ，优化 电离参 数从 而 降低 电 离
黑洞的准正模式(quasinormal modes)
Quasi-Normal Modes of Stars and Black HolesKostas D.KokkotasDepartment of Physics,Aristotle University of Thessaloniki,Thessaloniki54006,Greece.kokkotas@astro.auth.grhttp://www.astro.auth.gr/˜kokkotasandBernd G.SchmidtMax Planck Institute for Gravitational Physics,Albert Einstein Institute,D-14476Golm,Germany.bernd@aei-potsdam.mpg.dePublished16September1999/Articles/Volume2/1999-2kokkotasLiving Reviews in RelativityPublished by the Max Planck Institute for Gravitational PhysicsAlbert Einstein Institute,GermanyAbstractPerturbations of stars and black holes have been one of the main topics of relativistic astrophysics for the last few decades.They are of partic-ular importance today,because of their relevance to gravitational waveastronomy.In this review we present the theory of quasi-normal modes ofcompact objects from both the mathematical and astrophysical points ofview.The discussion includes perturbations of black holes(Schwarzschild,Reissner-Nordstr¨o m,Kerr and Kerr-Newman)and relativistic stars(non-rotating and slowly-rotating).The properties of the various families ofquasi-normal modes are described,and numerical techniques for calculat-ing quasi-normal modes reviewed.The successes,as well as the limits,of perturbation theory are presented,and its role in the emerging era ofnumerical relativity and supercomputers is discussed.c 1999Max-Planck-Gesellschaft and the authors.Further information on copyright is given at /Info/Copyright/.For permission to reproduce the article please contact livrev@aei-potsdam.mpg.de.Article AmendmentsOn author request a Living Reviews article can be amended to include errata and small additions to ensure that the most accurate and up-to-date infor-mation possible is provided.For detailed documentation of amendments, please go to the article’s online version at/Articles/Volume2/1999-2kokkotas/. Owing to the fact that a Living Reviews article can evolve over time,we recommend to cite the article as follows:Kokkotas,K.D.,and Schmidt,B.G.,“Quasi-Normal Modes of Stars and Black Holes”,Living Rev.Relativity,2,(1999),2.[Online Article]:cited on<date>, /Articles/Volume2/1999-2kokkotas/. The date in’cited on<date>’then uniquely identifies the version of the article you are referring to.3Quasi-Normal Modes of Stars and Black HolesContents1Introduction4 2Normal Modes–Quasi-Normal Modes–Resonances7 3Quasi-Normal Modes of Black Holes123.1Schwarzschild Black Holes (12)3.2Kerr Black Holes (17)3.3Stability and Completeness of Quasi-Normal Modes (20)4Quasi-Normal Modes of Relativistic Stars234.1Stellar Pulsations:The Theoretical Minimum (23)4.2Mode Analysis (26)4.2.1Families of Fluid Modes (26)4.2.2Families of Spacetime or w-Modes (30)4.3Stability (31)5Excitation and Detection of QNMs325.1Studies of Black Hole QNM Excitation (33)5.2Studies of Stellar QNM Excitation (34)5.3Detection of the QNM Ringing (37)5.4Parameter Estimation (39)6Numerical Techniques426.1Black Holes (42)6.1.1Evolving the Time Dependent Wave Equation (42)6.1.2Integration of the Time Independent Wave Equation (43)6.1.3WKB Methods (44)6.1.4The Method of Continued Fractions (44)6.2Relativistic Stars (45)7Where Are We Going?487.1Synergism Between Perturbation Theory and Numerical Relativity487.2Second Order Perturbations (48)7.3Mode Calculations (49)7.4The Detectors (49)8Acknowledgments50 9Appendix:Schr¨o dinger Equation Versus Wave Equation51Living Reviews in Relativity(1999-2)K.D.Kokkotas and B.G.Schmidt41IntroductionHelioseismology and asteroseismology are well known terms in classical astro-physics.From the beginning of the century the variability of Cepheids has been used for the accurate measurement of cosmic distances,while the variability of a number of stellar objects(RR Lyrae,Mira)has been associated with stel-lar oscillations.Observations of solar oscillations(with thousands of nonradial modes)have also revealed a wealth of information about the internal structure of the Sun[204].Practically every stellar object oscillates radially or nonradi-ally,and although there is great difficulty in observing such oscillations there are already results for various types of stars(O,B,...).All these types of pulsations of normal main sequence stars can be studied via Newtonian theory and they are of no importance for the forthcoming era of gravitational wave astronomy.The gravitational waves emitted by these stars are extremely weak and have very low frequencies(cf.for a discussion of the sun[70],and an im-portant new measurement of the sun’s quadrupole moment and its application in the measurement of the anomalous precession of Mercury’s perihelion[163]). This is not the case when we consider very compact stellar objects i.e.neutron stars and black holes.Their oscillations,produced mainly during the formation phase,can be strong enough to be detected by the gravitational wave detectors (LIGO,VIRGO,GEO600,SPHERE)which are under construction.In the framework of general relativity(GR)quasi-normal modes(QNM) arise,as perturbations(electromagnetic or gravitational)of stellar or black hole spacetimes.Due to the emission of gravitational waves there are no normal mode oscillations but instead the frequencies become“quasi-normal”(complex), with the real part representing the actual frequency of the oscillation and the imaginary part representing the damping.In this review we shall discuss the oscillations of neutron stars and black holes.The natural way to study these oscillations is by considering the linearized Einstein equations.Nevertheless,there has been recent work on nonlinear black hole perturbations[101,102,103,104,100]while,as yet nothing is known for nonlinear stellar oscillations in general relativity.The study of black hole perturbations was initiated by the pioneering work of Regge and Wheeler[173]in the late50s and was continued by Zerilli[212]. The perturbations of relativistic stars in GR werefirst studied in the late60s by Kip Thorne and his collaborators[202,198,199,200].The initial aim of Regge and Wheeler was to study the stability of a black hole to small perturbations and they did not try to connect these perturbations to astrophysics.In con-trast,for the case of relativistic stars,Thorne’s aim was to extend the known properties of Newtonian oscillation theory to general relativity,and to estimate the frequencies and the energy radiated as gravitational waves.QNMs werefirst pointed out by Vishveshwara[207]in calculations of the scattering of gravitational waves by a Schwarzschild black hole,while Press[164] coined the term quasi-normal frequencies.QNM oscillations have been found in perturbation calculations of particles falling into Schwarzschild[73]and Kerr black holes[76,80]and in the collapse of a star to form a black hole[66,67,68]. Living Reviews in Relativity(1999-2)5Quasi-Normal Modes of Stars and Black Holes Numerical investigations of the fully nonlinear equations of general relativity have provided results which agree with the results of perturbation calculations;in particular numerical studies of the head-on collision of two black holes [30,29](cf.Figure 1)and gravitational collapse to a Kerr hole [191].Recently,Price,Pullin and collaborators [170,31,101,28]have pushed forward the agreement between full nonlinear numerical results and results from perturbation theory for the collision of two black holes.This proves the power of the perturbation approach even in highly nonlinear problems while at the same time indicating its limits.In the concluding remarks of their pioneering paper on nonradial oscillations of neutron stars Thorne and Campollataro [202]described it as “just a modest introduction to a story which promises to be long,complicated and fascinating ”.The story has undoubtedly proved to be intriguing,and many authors have contributed to our present understanding of the pulsations of both black holes and neutron stars.Thirty years after these prophetic words by Thorne and Campollataro hundreds of papers have been written in an attempt to understand the stability,the characteristic frequencies and the mechanisms of excitation of these oscillations.Their relevance to the emission of gravitational waves was always the basic underlying reason of each study.An account of all this work will be attempted in the next sections hoping that the interested reader will find this review useful both as a guide to the literature and as an inspiration for future work on the open problems of the field.020406080100Time (M ADM )-0.3-0.2-0.10.00.10.20.3(l =2) Z e r i l l i F u n c t i o n Numerical solutionQNM fit Figure 1:QNM ringing after the head-on collision of two unequal mass black holes [29].The continuous line corresponds to the full nonlinear numerical calculation while the dotted line is a fit to the fundamental and first overtone QNM.In the next section we attempt to give a mathematical definition of QNMs.Living Reviews in Relativity (1999-2)K.D.Kokkotas and B.G.Schmidt6 The third and fourth section will be devoted to the study of the black hole and stellar QNMs.In thefifth section we discuss the excitation and observation of QNMs andfinally in the sixth section we will mention the more significant numerical techniques used in the study of QNMs.Living Reviews in Relativity(1999-2)7Quasi-Normal Modes of Stars and Black Holes 2Normal Modes–Quasi-Normal Modes–Res-onancesBefore discussing quasi-normal modes it is useful to remember what normal modes are!Compact classical linear oscillating systems such asfinite strings,mem-branes,or cavitiesfilled with electromagnetic radiation have preferred time harmonic states of motion(ωis real):χn(t,x)=e iωn tχn(x),n=1,2,3...,(1) if dissipation is neglected.(We assumeχto be some complex valuedfield.) There is generally an infinite collection of such periodic solutions,and the“gen-eral solution”can be expressed as a superposition,χ(t,x)=∞n=1a n e iωn tχn(x),(2)of such normal modes.The simplest example is a string of length L which isfixed at its ends.All such systems can be described by systems of partial differential equations of the type(χmay be a vector)∂χ∂t=Aχ,(3)where A is a linear operator acting only on the spatial variables.Because of thefiniteness of the system the time evolution is only determined if some boundary conditions are prescribed.The search for solutions periodic in time leads to a boundary value problem in the spatial variables.In simple cases it is of the Sturm-Liouville type.The treatment of such boundary value problems for differential equations played an important role in the development of Hilbert space techniques.A Hilbert space is chosen such that the differential operator becomes sym-metric.Due to the boundary conditions dictated by the physical problem,A becomes a self-adjoint operator on the appropriate Hilbert space and has a pure point spectrum.The eigenfunctions and eigenvalues determine the periodic solutions(1).The definition of self-adjointness is rather subtle from a physicist’s point of view since fairly complicated“domain issues”play an essential role.(See[43] where a mathematical exposition for physicists is given.)The wave equation modeling thefinite string has solutions of various degrees of differentiability. To describe all“realistic situations”,clearly C∞functions should be sufficient. Sometimes it may,however,also be convenient to consider more general solu-tions.From the mathematical point of view the collection of all smooth functions is not a natural setting to study the wave equation because sequences of solutionsLiving Reviews in Relativity(1999-2)K.D.Kokkotas and B.G.Schmidt8 exist which converge to non-smooth solutions.To establish such powerful state-ments like(2)one has to study the equation on certain subsets of the Hilbert space of square integrable functions.For“nice”equations it usually happens that the eigenfunctions are in fact analytic.They can then be used to gen-erate,for example,all smooth solutions by a pointwise converging series(2). The key point is that we need some mathematical sophistication to obtain the “completeness property”of the eigenfunctions.This picture of“normal modes”changes when we consider“open systems”which can lose energy to infinity.The simplest case are waves on an infinite string.The general solution of this problem isχ(t,x)=A(t−x)+B(t+x)(4) with“arbitrary”functions A and B.Which solutions should we study?Since we have all solutions,this is not a serious question.In more general cases, however,in which the general solution is not known,we have to select a certain class of solutions which we consider as relevant for the physical problem.Let us consider for the following discussion,as an example,a wave equation with a potential on the real line,∂2∂t2χ+ −∂2∂x2+V(x)χ=0.(5)Cauchy dataχ(0,x),∂tχ(0,x)which have two derivatives determine a unique twice differentiable solution.No boundary condition is needed at infinity to determine the time evolution of the data!This can be established by fairly simple PDE theory[116].There exist solutions for which the support of thefields are spatially compact, or–the other extreme–solutions with infinite total energy for which thefields grow at spatial infinity in a quite arbitrary way!From the point of view of physics smooth solutions with spatially compact support should be the relevant class–who cares what happens near infinity! Again it turns out that mathematically it is more convenient to study all solu-tions offinite total energy.Then the relevant operator is again self-adjoint,but now its spectrum is purely“continuous”.There are no eigenfunctions which are square integrable.Only“improper eigenfunctions”like plane waves exist.This expresses the fact that wefind a solution of the form(1)for any realωand by forming appropriate superpositions one can construct solutions which are “almost eigenfunctions”.(In the case V(x)≡0these are wave packets formed from plane waves.)These solutions are the analogs of normal modes for infinite systems.Let us now turn to the discussion of“quasi-normal modes”which are concep-tually different to normal modes.To define quasi-normal modes let us consider the wave equation(5)for potentials with V≥0which vanish for|x|>x0.Then in this case all solutions determined by data of compact support are bounded: |χ(t,x)|<C.We can use Laplace transformation techniques to represent such Living Reviews in Relativity(1999-2)9Quasi-Normal Modes of Stars and Black Holes solutions.The Laplace transformˆχ(s,x)(s>0real)of a solutionχ(t,x)isˆχ(s,x)= ∞0e−stχ(t,x)dt,(6) and satisfies the ordinary differential equations2ˆχ−ˆχ +Vˆχ=+sχ(0,x)+∂tχ(0,x),(7) wheres2ˆχ−ˆχ +Vˆχ=0(8) is the homogeneous equation.The boundedness ofχimplies thatˆχis analytic for positive,real s,and has an analytic continuation onto the complex half plane Re(s)>0.Which solutionˆχof this inhomogeneous equation gives the unique solution in spacetime determined by the data?There is no arbitrariness;only one of the Green functions for the inhomogeneous equation is correct!All Green functions can be constructed by the following well known method. Choose any two linearly independent solutions of the homogeneous equation f−(s,x)and f+(s,x),and defineG(s,x,x )=1W(s)f−(s,x )f+(s,x)(x <x),f−(s,x)f+(s,x )(x >x),(9)where W(s)is the Wronskian of f−and f+.If we denote the inhomogeneity of(7)by j,a solution of(7)isˆχ(s,x)= ∞−∞G(s,x,x )j(s,x )dx .(10) We still have to select a unique pair of solutions f−,f+.Here the information that the solution in spacetime is bounded can be used.The definition of the Laplace transform implies thatˆχis bounded as a function of x.Because the potential V vanishes for|x|>x0,the solutions of the homogeneous equation(8) for|x|>x0aref=e±sx.(11) The following pair of solutionsf+=e−sx for x>x0,f−=e+sx for x<−x0,(12) which is linearly independent for Re(s)>0,gives the unique Green function which defines a bounded solution for j of compact support.Note that for Re(s)>0the solution f+is exponentially decaying for large x and f−is expo-nentially decaying for small x.For small x however,f+will be a linear com-bination a(s)e−sx+b(s)e sx which will in general grow exponentially.Similar behavior is found for f−.Living Reviews in Relativity(1999-2)K.D.Kokkotas and B.G.Schmidt 10Quasi-Normal mode frequencies s n can be defined as those complex numbers for whichf +(s n ,x )=c (s n )f −(s n ,x ),(13)that is the two functions become linearly dependent,the Wronskian vanishes and the Green function is singular!The corresponding solutions f +(s n ,x )are called quasi eigenfunctions.Are there such numbers s n ?From the boundedness of the solution in space-time we know that the unique Green function must exist for Re (s )>0.Hence f +,f −are linearly independent for those values of s .However,as solutions f +,f −of the homogeneous equation (8)they have a unique continuation to the complex s plane.In [35]it is shown that for positive potentials with compact support there is always a countable number of zeros of the Wronskian with Re (s )<0.What is the mathematical and physical significance of the quasi-normal fre-quencies s n and the corresponding quasi-normal functions f +?First of all we should note that because of Re (s )<0the function f +grows exponentially for small and large x !The corresponding spacetime solution e s n t f +(s n ,x )is therefore not a physically relevant solution,unlike the normal modes.If one studies the inverse Laplace transformation and expresses χas a com-plex line integral (a >0),χ(t,x )=12πi +∞−∞e (a +is )t ˆχ(a +is,x )ds,(14)one can deform the path of the complex integration and show that the late time behavior of solutions can be approximated in finite parts of the space by a finite sum of the form χ(t,x )∼N n =1a n e (αn +iβn )t f +(s n ,x ).(15)Here we assume that Re (s n +1)<Re (s n )<0,s n =αn +iβn .The approxi-mation ∼means that if we choose x 0,x 1, and t 0then there exists a constant C (t 0,x 0,x 1, )such that χ(t,x )−N n =1a n e (αn +iβn )t f +(s n ,x ) ≤Ce (−|αN +1|+ )t (16)holds for t >t 0,x 0<x <x 1, >0with C (t 0,x 0,x 1, )independent of t .The constants a n depend only on the data [35]!This implies in particular that all solutions defined by data of compact support decay exponentially in time on spatially bounded regions.The generic leading order decay is determined by the quasi-normal mode frequency with the largest real part s 1,i.e.slowest damping.On finite intervals and for late times the solution is approximated by a finite sum of quasi eigenfunctions (15).It is presently unclear whether one can strengthen (16)to a statement like (2),a pointwise expansion of the late time solution in terms of quasi-normal Living Reviews in Relativity (1999-2)11Quasi-Normal Modes of Stars and Black Holes modes.For one particular potential(P¨o schl-Teller)this has been shown by Beyer[42].Let us now consider the case where the potential is positive for all x,but decays near infinity as happens for example for the wave equation on the static Schwarzschild spacetime.Data of compact support determine again solutions which are bounded[117].Hence we can proceed as before.Thefirst new point concerns the definitions of f±.It can be shown that the homogeneous equation(8)has for each real positive s a unique solution f+(s,x)such that lim x→∞(e sx f+(s,x))=1holds and correspondingly for f−.These functions are uniquely determined,define the correct Green function and have analytic continuations onto the complex half plane Re(s)>0.It is however quite complicated to get a good representation of these func-tions.If the point at infinity is not a regular singular point,we do not even get converging series expansions for f±.(This is particularly serious for values of s with negative real part because we expect exponential growth in x).The next new feature is that the analyticity properties of f±in the complex s plane depend on the decay of the potential.To obtain information about analytic continuation,even use of analyticity properties of the potential in x is made!Branch cuts may occur.Nevertheless in a lot of cases an infinite number of quasi-normal mode frequencies exists.The fact that the potential never vanishes may,however,destroy the expo-nential decay in time of the solutions and therefore the essential properties of the quasi-normal modes.This probably happens if the potential decays slower than exponentially.There is,however,the following way out:Suppose you want to study a solution determined by data of compact support from t=0to some largefinite time t=T.Up to this time the solution is–because of domain of dependence properties–completely independent of the potential for sufficiently large x.Hence we may see an exponential decay of the form(15)in a time range t1<t<T.This is the behavior seen in numerical calculations.The situation is similar in the case ofα-decay in quantum mechanics.A comparison of quasi-normal modes of wave equations and resonances in quantum theory can be found in the appendix,see section9.Living Reviews in Relativity(1999-2)K.D.Kokkotas and B.G.Schmidt123Quasi-Normal Modes of Black HolesOne of the most interesting aspects of gravitational wave detection will be the connection with the existence of black holes[201].Although there are presently several indirect ways of identifying black holes in the universe,gravitational waves emitted by an oscillating black hole will carry a uniquefingerprint which would lead to the direct identification of their existence.As we mentioned earlier,gravitational radiation from black hole oscillations exhibits certain characteristic frequencies which are independent of the pro-cesses giving rise to these oscillations.These“quasi-normal”frequencies are directly connected to the parameters of the black hole(mass,charge and angu-lar momentum)and for stellar mass black holes are expected to be inside the bandwidth of the constructed gravitational wave detectors.The perturbations of a Schwarzschild black hole reduce to a simple wave equation which has been studied extensively.The wave equation for the case of a Reissner-Nordstr¨o m black hole is more or less similar to the Schwarzschild case,but for Kerr one has to solve a system of coupled wave equations(one for the radial part and one for the angular part).For this reason the Kerr case has been studied less thoroughly.Finally,in the case of Kerr-Newman black holes we face the problem that the perturbations cannot be separated in their angular and radial parts and thus apart from special cases[124]the problem has not been studied at all.3.1Schwarzschild Black HolesThe study of perturbations of Schwarzschild black holes assumes a small per-turbation hµνon a static spherically symmetric background metricds2=g0µνdxµdxν=−e v(r)dt2+eλ(r)dr2+r2 dθ2+sin2θdφ2 ,(17) with the perturbed metric having the formgµν=g0µν+hµν,(18) which leads to a variation of the Einstein equations i.e.δGµν=4πδTµν.(19) By assuming a decomposition into tensor spherical harmonics for each hµνof the formχ(t,r,θ,φ)= mχ m(r,t)r Y m(θ,φ),(20)the perturbation problem is reduced to a single wave equation,for the func-tionχ m(r,t)(which is a combination of the various components of hµν).It should be pointed out that equation(20)is an expansion for scalar quantities only.From the10independent components of the hµνonly h tt,h tr,and h rr transform as scalars under rotations.The h tθ,h tφ,h rθ,and h rφtransform asLiving Reviews in Relativity(1999-2)13Quasi-Normal Modes of Stars and Black Holes components of two-vectors under rotations and can be expanded in a series of vector spherical harmonics while the components hθθ,hθφ,and hφφtransform as components of a2×2tensor and can be expanded in a series of tensor spher-ical harmonics(see[202,212,152]for details).There are two classes of vector spherical harmonics(polar and axial)which are build out of combinations of the Levi-Civita volume form and the gradient operator acting on the scalar spherical harmonics.The difference between the two families is their parity. Under the parity operatorπa spherical harmonic with index transforms as (−1) ,the polar class of perturbations transform under parity in the same way, as(−1) ,and the axial perturbations as(−1) +11.Finally,since we are dealing with spherically symmetric spacetimes the solution will be independent of m, thus this subscript can be omitted.The radial component of a perturbation outside the event horizon satisfies the following wave equation,∂2∂t χ + −∂2∂r∗+V (r)χ =0,(21)where r∗is the“tortoise”radial coordinate defined byr∗=r+2M log(r/2M−1),(22) and M is the mass of the black hole.For“axial”perturbationsV (r)= 1−2M r ( +1)r+2σMr(23)is the effective potential or(as it is known in the literature)Regge-Wheeler potential[173],which is a single potential barrier with a peak around r=3M, which is the location of the unstable photon orbit.The form(23)is true even if we consider scalar or electromagnetic testfields as perturbations.The parameter σtakes the values1for scalar perturbations,0for electromagnetic perturbations, and−3for gravitational perturbations and can be expressed asσ=1−s2,where s=0,1,2is the spin of the perturbingfield.For“polar”perturbations the effective potential was derived by Zerilli[212]and has the form V (r)= 1−2M r 2n2(n+1)r3+6n2Mr2+18nM2r+18M3r3(nr+3M)2,(24)1In the literature the polar perturbations are also called even-parity because they are characterized by their behavior under parity operations as discussed earlier,and in the same way the axial perturbations are called odd-parity.We will stick to the polar/axial terminology since there is a confusion with the definition of the parity operation,the reason is that to most people,the words“even”and“odd”imply that a mode transforms underπas(−1)2n or(−1)2n+1respectively(for n some integer).However only the polar modes with even have even parity and only axial modes with even have odd parity.If is odd,then polar modes have odd parity and axial modes have even parity.Another terminology is to call the polar perturbations spheroidal and the axial ones toroidal.This definition is coming from the study of stellar pulsations in Newtonian theory and represents the type offluid motions that each type of perturbation induces.Since we are dealing both with stars and black holes we will stick to the polar/axial terminology.Living Reviews in Relativity(1999-2)K.D.Kokkotas and B.G.Schmidt14where2n=( −1)( +2).(25) Chandrasekhar[54]has shown that one can transform the equation(21)for “axial”modes to the corresponding one for“polar”modes via a transforma-tion involving differential operations.It can also be shown that both forms are connected to the Bardeen-Press[38]perturbation equation derived via the Newman-Penrose formalism.The potential V (r∗)decays exponentially near the horizon,r∗→−∞,and as r−2∗for r∗→+∞.From the form of equation(21)it is evident that the study of black hole perturbations will follow the footsteps of the theory outlined in section2.Kay and Wald[117]have shown that solutions with data of compact sup-port are bounded.Hence we know that the time independent Green function G(s,r∗,r ∗)is analytic for Re(s)>0.The essential difficulty is now to obtain the solutions f±(cf.equation(10))of the equations2ˆχ−ˆχ +Vˆχ=0,(26) (prime denotes differentiation with respect to r∗)which satisfy for real,positives:f+∼e−sr∗for r∗→∞,f−∼e+r∗x for r∗→−∞.(27) To determine the quasi-normal modes we need the analytic continuations of these functions.As the horizon(r∗→∞)is a regular singular point of(26),a representation of f−(r∗,s)as a converging series exists.For M=12it reads:f−(r,s)=(r−1)s∞n=0a n(s)(r−1)n.(28)The series converges for all complex s and|r−1|<1[162].(The analytic extension of f−is investigated in[115].)The result is that f−has an extension to the complex s plane with poles only at negative real integers.The representation of f+is more complicated:Because infinity is a singular point no power series expansion like(28)exists.A representation coming from the iteration of the defining integral equation is given by Jensen and Candelas[115],see also[159]. It turns out that the continuation of f+has a branch cut Re(s)≤0due to the decay r−2for large r[115].The most extensive mathematical investigation of quasi-normal modes of the Schwarzschild solution is contained in the paper by Bachelot and Motet-Bachelot[35].Here the existence of an infinite number of quasi-normal modes is demonstrated.Truncating the potential(23)to make it of compact support leads to the estimate(16).The decay of solutions in time is not exponential because of the weak decay of the potential for large r.At late times,the quasi-normal oscillations are swamped by the radiative tail[166,167].This tail radiation is of interest in its Living Reviews in Relativity(1999-2)。
超大质量双黑洞系统中的吸积盘Formationand
的而FRI射电星系中的吸积盘则很可能是非标准的
•观测：大约 7%的射电源属于X-形态源 [Leahy & Parma 92] •理论：喷流等离子体遗迹可探测时间大约一千万年，7%的探测概率意味 着低亮度的FRII射电星系寿命大约一亿年。这和测得的AGN寿命一致
rL
1
1 m
2/3
rs
m=
if 104
The whole inner region is destroyed by waves
if 104 and H/r<0.1
Opens a gap in the accretion disk The outer boundary of the gap:
evolution merging
Galaxy + galaxy
Normal galaxy
星系与活动星系的大统一模型（？）
gas: starbursts galaxy Evolution Short-lived AGNs
Martini 2004 106 yr tQ a few 108 yr
1.3、Gravitational Wave Astrophysics:
•观测：喷流遗迹与寄住星系盘垂直 •理论：在吸积盘与超大质量双黑洞发生相互作用前，吸积盘因角动量守恒 而与星系盘共面，因而产生的喷流与吸积盘和星系盘垂直
•观测：活动的喷流随机取向 •理论：新喷流方向垂直于双黑洞轨道面，而星系并合形成的超大质量双黑 洞相对于星系盘取向随机
Random distribution of jet orientations (Liu, Zhang, Wu, 2005, in preparation)
UBVI CCD photometry of two old open clusters NGC 1798 and NGC 2192
a r X i v :a s t r o -p h /9901140v 1 12 J a n 1999Mon.Not.R.Astron.Soc.000,1–8(1998)Printed 1February 2008(MN L A T E X style file v1.4)UBV I CCD photometry of two old open clusters NGC 1798and NGC 2192Hong Soo Park and Myung Gyoon Lee ⋆Department of Astronomy,Seoul National University,Seoul 151-742,KoreaAccepted 1998??.Received 1998??;in original form 1998July ??ABSTRACT We present UBV I CCD photometry of two open clusters NGC 1798and NGC 2192which were little studied before.Color-magnitude diagrams of these clusters show several features typical for old open clusters:a well-defined main-sequence,a red giant clump,and a small number of red giants.The main sequence of NGC 1798shows a distinct gap at V ≈16.2mag.From the surface number density distribution we have measured the size of the clusters,obtaining 8′.3(=10.2pc)for NGC 1798and 7′.3(=7.5pc)for NGC 2192.Then we have determined the reddening,metallicity,and distance of these clusters using the color-color diagrams and color-magnitude diagrams:E (B −V )=0.51±0.04,[Fe/H]=−0.47±0.15dex and (m −M )0=13.1±0.2(d =4.2±0.3kpc)for NGC 1798,and E (B −V )=0.20±0.03,[Fe/H]=−0.31±0.15dex and (m −M )0=12.7±0.2(d =3.5±0.3kpc)for NGC 2192.The ages of these clusters have been estimated using the morphological age indicators and the isochrone fitting with the Padova isochrones:1.4±0.2Gyrs for NGC 1798and 1.1±0.1Gyrs for NGC 2192.The luminosity functions of the main sequence stars in these clusters are found to be similar to other old open clusters.The metallicity and distance of these clusters are consistent with the relation between the metallicity and galactocentric distance of other old open clusters.Key words:Hertzsprung–Russell (HR)diagram –open clusters and associations:general –open clusters and associations:individual:NGC 1798,NGC 2192–stars:luminosity function.1INTRODUCTIONOld open clusters provide us with an important informationfor understanding the early evolution of the Galactic disk.There are about 70known old open clusters with age >1Gyrs (Friel 1995).These clusters are faint in general so thatthere were few studies about these clusters until recently.With the advent of CCD camera in astronomy,the numberof studies on these clusters has been increasing.However,there are still a significant number of old open clusters forwhich basic parameters are not well known.For example,metallicity is not yet known for about 30clusters amongthem.Recently Phelps et al.(1994)and Janes &Phelps (1994)presented an extensive CCD photometric survey of potentialold open clusters,the results of which were used in the studyon the development of the Galactic disk by Janes &Phelps(1994).In the sample of the clusters studied by Phelps et al.⋆corresponding author,E-mail:mglee@astrog.snu.ac.krthere are several clusters for which only the non-calibrated photometry is available.We have chosen two clusters among them,NGC 1798and NGC 2192,to study the characteristics of these clus-ters using UBV I CCD photometry.These clusters are lo-cated in the direction of anti-galactic centre.To date there is published only one photometric study of these clusters,which was given by Phelps et al.(1994)who presented non-calibrated BV CCD photometry of these clusters.From the instrumental color-magnitude diagrams of these clus-ters Phelps et al.estimated the ages of these clusters using the morphological age indicators,obtaining the values of 1.5Gyrs for NGC 1798and 1.1Gyrs for NGC 2192.However,no other properties of these clusters are yet known.In this paper we present a study of NGC 1798and NGC 2192based on UBV I CCD photometry.We have es-timated the basic parameters of these clusters:size,red-dening,metallicity,distance,and age.Also we have derived the luminosity function of the main sequence stars in these clusters.Section 2describes the observations and data re-duction.Sections 3and 4present the analysis for NGC 1798c 1998RAS2Hong Soo Park and Myung Gyoon Leeand NGC2192,respectively.Section5discusses the results. Finally Section6summarizes the primary results.2OBSER V ATIONS AND DATA REDUCTION 2.1ObservationsUBV I CCD images of NGC1798and NGC2192were obtained using the Photometrics512CCD camera at the Sobaeksan Observatory61cm telescope in Korea for several observing runs between1996November and1997October. We have used also BV CCD images of the central region of NGC1798obtained by Chul Hee Kim using the Tek1024 CCD camera at the Vainu Bappu Observatory2.3m tele-scope in India on March4,1998.The observing log is given in Table1.The original CCD images wereflattened after bias sub-traction and several exposures for eachfilter were combined into a single image for further reduction.The sizes of the field in a CCD image are4′.3×4′.3for the PM512CCD im-age,and10′.6×10′.6for the Tek1024CCD image.The gain and readout noise are,respectively,9electrons/ADU and 10.4electrons for the PM512CCD,and9electrons/ADU and10.4electrons for the Tek1024CCD.Figs.1and2illustrate grey scale maps of the V CCD images of NGC1798and NGC2192made by mosaicing the images of the observed regions.It is seen from thesefigures that NGC1798is a relatively rich open cluster,while NGC 2192is a relatively poor open cluster.2.2Data ReductionInstrumental magnitudes of the stars in the CCD images were obtained using the digital stellar photometry reduction program IRAF⋆/DAOPHOT(Stetson1987,Davis1994). The resulting instrumental magnitudes were transformed onto the standard system using the standard stars from Lan-dolt(1992)and the M67stars in Montgomery et al.(1993) observed on the same photometric nights.The transforma-tion equations areV=v+a V(b−v)+k V X+Z V,(B−V)=a BV(b−v)+k BV X+Z BV,(U−B)=a UB(u−b)+k UB X+Z UB,and(V−I)=a V I(v−i)+k V I X+Z V I,where the lower case symbols represent instrumental mag-nitudes derived from the CCD images and the upper case symbols represent the standard system values.X is the air-mass at the midpoint of the observations.The results of the transformation are summarized in Table2.The data ob-tained on non-photometric nights were calibrated using the photometric data for the overlapped region.The total number of the measured stars is1,416for NGC1798and409for NGC2192.Tables3and4list the⋆IRAF is distributed by the National Optical Astronomy Ob-servatories,which are operated by the Association of Universities for Research in Astronomy,Inc.under contract with the National Science Foundation.photometry of the bright stars in the C-regions of NGC1798 and NGC2192,respectively.The X and Y coordinates listed in Table3and4are given in units of CCD pixel(=0′′.50). The X and Y values are increasing toward north and west, respectively.We have divided the entire region of thefields into sev-eral regions,as shown in Figs.1and2,for the analysis of the data.The C-region represents the central region of the cluster,and the F-regions(F,Fb,Fir,and Fi regions)rep-resent the controlfield regions,and the N-region represents the intermediate region between the central region and the field region.The radius of the C-region is300pixel for NGC 1798and NGC2192.The ratio of the areas of the C-region, N-region,Fb-region,and(Fi+Fir)-regions for NGC1798is 1:1.50:1.00:1.07,and the ratio of the areas of the C-region, N-region,and F-region for NGC2192is1:1.26:0.98.3ANALYSIS FOR NGC17983.1The Size of the ClusterWe have investigated the structure of NGC1798using star-counts.The centre of the cluster is estimated to be at the po-sition of(X=710pixel,Y=1110pixel),using the centroid method.Fig.3illustrates the projected surface number den-sity profile derived from counting stars with V<19.5mag in the entire CCDfield.The magnitude cutofffor starcounts was set so that the counts should be free of any photometric incompleteness problem.Fig.3shows that most of the stars in NGC1798are concentrated within the radius of250pixel (=125′′),and that the outskirts of the cluster extend out to about500pixel(=250′′)from the center.The number density changes little with radius beyond500pixel,show-ing that the outer region of the observedfield can be used as a controlfield.Therefore we have estimated the approxi-mate size of NGC1798for which the cluster blends in with thefield to be about500′′in diameter,which corresponds to a linear size of10.2pc for the distance of NGC1798as determined below.3.2Color-Magnitude DiagramsFigs.4and5show the V−(B−V)and V−(V−I)color-magnitude diagrams(CMDs)of the measured stars in the observed regions in NGC1798.Thesefigures show that the C-region consists mostly of the members of NGC1798with some contamination of thefield stars,while the F-regions consist mostly of thefield stars.The N-region is intermediate between the C-region and the F-region.The distinguishable features seen in the color-magnitude diagrams of the C-region are:(a)There is a well-defined main sequence the top of which is located at V≈16 mag;(b)There is seen a distinct gap at V≈16.2mag in the main sequence,which is often seen in other old open clusters (e.g.M67);(c)There is a poorly defined red giant branch and these is seen some excess of stars around(B−V)=1.3 and V=15.6mag on this giant branch,which is remarked by the small box in thefigures.This may be a random excess of stars.However,the positions of the stars in the CMDs are consistent with the positions of known red giant clump in other old open clusters.Therefore most of these stars arec 1998RAS,MNRAS000,1–8UBV I CCD photometry of two old open clusters NGC1798and NGC21923probably red giant clump stars;and(d)There are a small number of stars along the locus of the red giant branch. 3.3Reddening and MetallicityNGC1798is located close to the galactic plane in the anti-galactic centre direction(b=4◦.85and l=160◦.76)so that it is expected that the reddening toward this cluster is significant.We have estimated the reddening for NGC1798 using two methods as follows.First we have used the mean color of the red giant clump.Janes&Phelps(1994)estimated the mean color and magnitude of the red giant clump in old open clusters to be (B−V)RGC=0.87±0.02and M V,RGC=0.59±0.09,when the difference between the red giant clump and the main sequence turn-offof the clusters,δV,is smaller than one. The mean color of the red giant clump in the C-region is estimated to be(B−V)RGC=1.34±0.01((V−I)RGC= 1.47±0.01,and(U−B)RGC=1.62±0.04),and the cor-responding mean magnitude is V RGC=15.57±0.05.δV is estimated to be0.8±0.2,which is the same value derived by Phelps et al.(1994).From these data we have derived a value of the reddening,E(B−V)=0.47±0.02.Secondly we have used the color-color diagram to es-timate the reddening and the metallicity simultaneously. We havefitted the mean colors of the stars in the C-region with the color-color relation used in the Padova isochrones (Bertelli et al.1994).This process requires iteration,because we need to know the age of the cluster as well as the red-dening and metallicity.We have iterated this process until all three parameters are stabilized.Fig.6illustrates the results offitting in the(U−B)−(B−V)color-color diagram.It is shown in thisfigure that the stars in NGC1798are reasonablyfitted by the color-color relation of the isochrones for[Fe/H]=−0.47±0.15 with a reddening value of E(B−V)=0.55±0.05. The error for the metallicity,0.15,was estimated by com-paring isochrones with different metallicities.As a reference the mean locus of the giants for solar abundance given by Schmidt-Kaler(1982)is also plotted in Fig.6.Finally we derive a mean value of the two estimates for the reddening, E(B−V)=0.51±0.04.3.4DistanceWe have estimated the distance to NGC1798using two methods as follows.First we have used the mean magnitude of the red giant clump.We have derived a value of the ap-parent distance modulus(m−M)V=14.98±0.10from the values for the mean magnitudes of the red giant clump stars described above.Secondly we have used the the zero-age main sequence (ZAMS)fitting,following the method described in Vanden-Berg&Poll(1989).VandenBerg&Poll(1989)presented the semi-empirical ZAMS as a function of the metallicity[Fe/H] and the helium abundance Y:V=M V(B−V)+δM V(Y)+δM V([Fe/H])whereδM V(Y)=2.6(Y−0.27)andδM V([Fe/H])=−[Fe/H](1.444+0.362[Fe/H]).Before the ZAMSfitting,we subtracted statistically the contribution due to thefield stars in the CMDs of the C-region using the CMDs of the Fb-region for BV photometry and the CMDs of the Fi+Fir region for V I photometry. The size of the bin used for the subtraction is∆V=0.25 and∆(B−V)=0.1.The resulting CMDs are displayed in Fig.7.We used the metallicity of[Fe/H]=–0.47as derived above and adopted Y=0.28which is the mean value for old open clusters(Gratton1982).Using this method we have obtained a value of the apparent distance modulus(m−M)V=14.5±0.2.Finally we calculate a mean value of the two estimates,(m−M)V=14.7±0.2.Adopting the extinction law of A V=3.2E(B−V),we derive a value of the intrinsic distance modulus(m−M)0=13.1±0.2.This corresponds to a distance of d=4.2±0.3kpc.3.5AgeWe have estimated the age of NGC1798using two methods as follows.First we have used the morphological age index (MAI)as described in Phelps et al.(1994).Phelps et al. (1994)and Janes&Phelps(1994)presented the MAI–δV relation,MAI[Gyrs]=0.73×10(0.256δV+0.0662δV2).From the value ofδV derived above,0.8±0.2mag,we obtain a value for the age,MAI=1.3±0.2Gyrs.Secondly we have estimated the age of the cluster us-ing the theoretical isochrones given by the Padova group (Bertelli et al.1994).Fitting the isochrones for[Fe/H]=–0.47to the CMDs of NGC1798,as shown in Fig.8,we estimate the age to be1.4±0.2Gyrs.Both results agree very well.3.6Luminosity FunctionWe have derived the V luminosity functions of the main sequence stars in NGC1798,which are displayed in Fig.9.The Fb-region was used for subtraction of thefield star contribution from the C-region and the magnitude bin size used is0.5mag.This controlfield may not be far enough from the cluster to derive thefield star contribution.If so,we might have oversubtracted thefield contribution,obtaining flatter luminosity functions than true luminosity functions. However,the fraction of the cluster members in thisfield must be,if any,very low,because the surface number density of this region is almost constant with the radius as shown in Fig.3.The luminosity function of the C-region in Fig. 9(a)increases rapidly up to V≈16.5mag,and stays almost flat for V>16.5mag.The luminosity functions of the N-region and the(R+Fir)-region are steeper than that of the C-region.A remarkable drop is seen at V=16.2mag(M V= 1.5mag)in the luminosity function of the C-region based on smaller bin size of0.2mag in Fig.9(b).This corresponds to the main sequence gap described above.4ANALYSIS FOR NGC21924.1The Size of the ClusterWe have investigated the structure of NGC2192using star-counts.We could not use the centroid method to estimatec 1998RAS,MNRAS000,1–84Hong Soo Park and Myung Gyoon Leethe centre of this cluster,because this cluster is too sparse. So we have used eye-estimate to determine the centre of the cluster to be at the position of(X=465pixel,Y=930 pixel).Fig.10illustrates the projected surface number den-sity profile derived from counting stars with V<18mag in the entire CCDfield.The magnitude cutofffor starcounts was set so that the counts should be free of any photomet-ric incompleteness problem.Fig.10shows that most of the stars in NGC2192are concentrated within the radius of200 pixel(=100′′),and that the outskirts of the cluster extend out to about440pixel(=220′′)from the centre.Therefore the approximate size of NGC2192is estimated to be about 440′′in diameter,which corresponds to a linear size of7.5pc for the distance of NGC2192as determined below.4.2Color-Magnitude DiagramsFigs.11and12show the V−(B−V)and V−(V−I)color-magnitude diagrams of the measured stars in the observed regions in NGC2192.The distinguishable features seen in the color-magnitude diagrams of the C-region are:(a)There is a well-defined main sequence the top of which is located at V≈14mag;(b)There are a group of red giant clump stars at(B−V)=1.1and V=14.2mag,which are remarked by the small box in thefigures;and(c)There are a small number of stars along the locus of the red giant branch. 4.3Reddening and MetallicityNGC2192is located11degrees above the galactic plane in the anti-galactic centre direction(b=10◦.64and l= 173◦.41)but higher than NGC1798so that it is expected that the reddening toward this cluster is significant but smaller than that of NGC1798.We have estimated the red-dening for NGC2192using two methods as applied for NGC 1798.First we have used the mean color of the red giant clump.The mean color of the red giant clump in the C-region is estimated to be(B−V)RGC=1.08±0.01((V−I)RGC=1.07±0.01,and(U−B)RGC=0.61±0.02),and the corresponding mean magnitude is V RGC=14.20±0.05.δV is estimated to be0.6±0.2,which is similar to the value derived by Phelps et al.(1994).From these data we have derived a value of the reddening,E(B−V)=0.19±0.03.Secondly we have used the color-color diagram to es-timate the reddening and the metallicity simultaneously. We havefitted the mean colors of the stars in the C-region with the color-color relation used in the Padova isochrones (Bertelli et al.1994).Fig.13illustrates the results offit-ting in the(U−B)−(B−V)color-color diagram.It is shown in thisfigure that the stars in NGC2192are rea-sonablyfitted by the color-color relation of the isochrones for[Fe/H]=−0.31±0.15dex with a reddening value of E(B−V)=0.21±0.01.The error for the metallicity,0.15, was estimated by comparing isochrones with different metal-licities.As a reference the mean locus of the giant for solar abundance given by Schmidt-Kaler is also plotted in Fig.13. Finally we derive a mean value of the two estimates for the reddening,E(B−V)=0.20±0.03.4.4DistanceWe have estimated the distance to NGC2192using two methods as for NGC1798.First we have used the mean magnitude of the red giant clump.We have derived a value of the apparent distance modulus(m−M)V=13.61±0.10 from the values for the mean magnitudes of the red giant clump stars described previously.Secondly we have used the ZAMSfitting.Before the ZAMSfitting,we subtracted statistically the contribution due to thefield stars in the CMDs of the C-region using the CMDs of the F-region.The size of the bin used for the subtraction is∆V=0.25and∆(B−V)=0.1.The resulting CMDs are displayed in Fig.14.We used the metallicity of[Fe/H]=–0.31as derived before and adopted Y=ing this method we have obtaineda value of the apparent distance modulus(m−M)V=13.1±0.2.Finally we calculate a mean value of the two estimates,(m−M)V=13.3±0.2.Adopting the extinction law of A V=3.2E(B−V),we derive a value of the intrinsic distance modulus(m−M)0=12.7±0.2.This corresponds to a distance of d=3.5±0.3kpc.4.5AgeWe have estimated the age of NGC2192using two methods as follows.First we have used the morphological age index. From the value ofδV derived above,0.6±0.2mag,we obtain a value for the age,MAI=1.1±0.2Gyrs.Secondly we have estimated the age of the cluster using the theoretical isochrones given by the Padova group(Bertelli et al.1994). Fitting the isochrones for[Fe/H]=–0.31to the CMDs of NGC2192,as shown in Fig.15,we estimate the age to be 1.1±0.1Gyrs.Both results agree very well.4.6Luminosity FunctionWe have derived the V luminosity functions of the main sequence stars in NGC2192,which are displayed in Fig.16.The F-region was used for subtraction of thefield star contribution from the C-region.The luminosity function of the C-region in Fig.16(a)increases rapidly up to V≈14 mag,and stays almostflat for V>15mag.The luminosity function of the N-region is steeper than that of the C-region. Fig.16(b)displays a comparison of the luminosity functions of NGC1798,NGC2192,and NGC7789which is another old open cluster of similar age(Roger et al.1994).Fig.16(b) shows that the luminosity functions of these clusters are similar in that they are almostflat in the faint part.The flattening of the faint part of the luminosity functions of old open clusters has been known since long,and is believed to be due to evaporation of low mass stars(Friel1995).5DISCUSSIONWe have determined the metallicity and distance of NGC 1798and NGC2192in this study.We compare them with those of other old open clusters here.Fig.17illustrates the radial metallicity gradient of the old open clusters com-piled by Friel(1995)and supplemented by the data in Wee&Lee(1996)and Lee(1997).Fig.17shows that thec 1998RAS,MNRAS000,1–8UBV I CCD photometry of two old open clusters NGC1798and NGC21925 mean metallicity decreases as the galactocentric distance in-creases.The positions of NGC1798and NGC2192we haveobtained in this study are consistent with the mean trendof the other old open clusters.The slope we have deter-mined for the entire sample including these two clusters is∆[Fe/H]/R GC=−0.086±0.011dex/kpc,very similar tothat given in Friel(1995),∆[Fe/H]/R GC=−0.091±0.014dex/kpc.There are only four old open clusters located beyondR GC=13kpc in Fig.17.These four clusters follow the meantrend of decreasing outward.However,the number of theclusters is not large enough to decide whether the metallictykeeps decreasing outward or it stops decreasing somewherebeyond R GC=13kpc and stays constant.Further studiesof more old open clusters beyond R GC=13kpc are neededto investigate this point.6SUMMARY AND CONCLUSIONSWe have presented UBV I photometry of old open clustersNGC1798and NGC2192.From the photometric data wehave determined the size,reddening,metallicity,distance,and age of these clusters.The luminosity functions of themain sequence stars in these clusters are similar to those ofthe other old open clusters.The basic properties of theseclusters we have determined in this study are summarizedin Table5.ACKNOWLEDGMENTSProf.Chul Hee Kim is thanked for providing the BV CCDimages of NGC1798.This research is supported in part bythe Korea Science and Engineering Foundation Grant No.95-0702-01-01-3.REFERENCESBertelli,G.,Bressan,A.,Chiosi,C.,Fagotto,F.,&Nasi,E.1994,A&AS,106,275Davis,L.E.,1994,A Reference Guide to the IRAF/DAOPHOTPackageFriel,E.D.1995,ARA&A,33,381Gratton,R.G.,1982,ApJ,257,640Janes,K.,&Phelps,R.L.1994,AJ,108,1773Landolt,A.U.,1992,AJ,104,340Lee,M.G.,1997,113,729Montgomery,K.A.,Marschall,L.A.,&Janes,K.A.1993,AJ,106,181Phelps,R.L.,Janes,K.,&Montgomery,K.A.,1994,AJ,107,1079Roger,C.M.,Paez,E.,Castellani,V.,&Staniero,O.1994,A&A,290,62Schmidt-Kaler,T.1982,in Landolt-Bornstein VI,2b(Berlin:Springer)Stetson,P.B.,1987,PASP,99,191VandenBerg,D.A.,&Poll,H.E.,1989,AJ,98,1451Wee,S.O.,&Lee,M.G.,1996,Jour.Korean Astro.Soc.,29,181This paper has been produced using the Royal AstronomicalSociety/Blackwell Science L A T E X stylefile.c 1998RAS,MNRAS000,1–86Hong Soo Park and Myung Gyoon LeeTable2.Transformation coefficients for the standard stars.Date Color a k Z rms n(stars) 96.11.11V0.003–0.101–6.0400.00919(B−V) 1.090–0.155–0.4670.01317(U−B) 1.008–0.154–1.7110.03416 97.02.11V0.028–0.198–6.0500.00934(B−V) 1.150–0.118–0.5930.01035(U−B) 1.079–0.324–1.5750.03226(V−I)0.983–0.1300.3020.00829 97.02.13V–0.007–0.176–6.0050.01034(B−V) 1.160–0.133–0.5880.01334(U−B) 1.079–0.349–1.5260.02320(V−I)0.986–0.0970.2710.01532 97.03.17V-0.019–0.225–6.0680.01548(B−V) 1.221–0.093–0.7810.01844(U−B) 1.008–0.309–1.4580.02735(V−I)0.956–0.1570.3140.02052Table 3.UBV I photometry of the bright stars in the C-regionof NGC1798.ID X[px]Y[px]V(B−V)(U−B)(V−I) 1723.21110.815.8600.7790.3240.9143680.51092.015.380 1.3250.627 1.5254746.01160.315.9050.7410.2200.9545741.81165.514.380 1.4830.880 1.6336684.61030.115.8140.8640.342 1.0427745.31073.315.3760.9200.237 1.0849707.21038.814.777 1.4860.723 1.653 11549.41153.715.305 1.5290.915 1.666 12547.91162.714.846 1.6270.947 1.764 13602.31068.015.325 1.5240.847 1.660 14598.11080.815.776 1.3340.575 1.490 15525.81116.815.510 1.0020.346 1.169 18833.91144.214.976 1.4940.866 1.578 21981.91125.114.711 1.510 1.031 1.440 26835.71098.015.706 1.3180.468 1.407 27656.51022.115.653 1.4160.180 1.571 28659.8955.915.216 1.2160.384 1.385 29656.4976.813.6670.6610.3100.777 30674.2822.512.0330.4450.0500.532 33668.7966.315.687 1.3130.641 1.446 34764.3879.915.460 1.3370.508 1.408 38749.01008.115.359 1.3220.744 1.426 39774.4965.915.545 1.3040.590 1.402 43571.8930.612.875 1.970 1.264 2.252 45512.8908.715.810 1.3480.474 1.485 46642.1949.515.622 1.3290.590 1.477 48732.61211.315.821 1.3530.643 1.492 50678.81210.515.297 1.3680.770 1.501 51742.11244.414.7990.747-0.046 1.216 52673.11250.215.180 1.3060.780 1.457 56728.21183.415.665 1.4050.803 1.545 57614.91229.413.085 2.1770.635 3.581 62563.01338.415.9090.7060.6380.778 65913.91203.013.686 1.2490.851 1.354 67792.71314.215.6010.9750.470 1.275 68932.61268.414.4660.7170.1970.763 70801.11338.715.561 1.3270.785 1.442 72827.61315.513.5260.9670.521 1.161 78444.01094.814.605 1.5790.825 1.723 82684.3822.615.8280.9510.128 1.189 97598.51330.315.580 1.3360.710 1.483Table4.UBV I photometry of the bright stars in the C-region of NGC2192.ID X[px]Y[px]V(B−V)(U−B)(V−I) 4502.4879.614.9830.4320.1140.424 5537.5891.712.9690.9240.4750.930 7535.1931.315.1700.598-0.2010.483 9434.8947.613.6760.570-0.0560.648 10562.3991.914.037 1.0960.612 1.043 11412.9884.814.165 1.0780.608 1.072 12598.4778.715.1840.464-0.0300.477 13595.9930.014.4800.5020.1720.549 14625.8971.814.7430.7200.0890.741 16614.3684.214.3970.7660.1750.758 17463.3689.514.7660.4660.0120.480 18553.3710.715.4920.447-0.1350.428 19591.7723.814.4660.5450.0930.566 21623.0867.712.634 1.1020.713 1.108 22648.0929.015.3380.3940.1150.399 26461.1711.314.162 1.0860.670 1.091 29583.9802.315.0940.600-0.0160.644 32720.4889.515.3260.3700.0710.410 37395.9778.315.2640.3960.0310.379 38405.5778.514.8900.4810.1010.409 40258.21042.915.4310.4240.1200.484 42376.6883.915.1870.4550.1360.470 44180.1908.714.3200.2870.1120.295 47367.1938.314.8970.4030.2520.459 51183.3998.714.202 1.1140.552 1.108 55475.81052.514.6370.3880.1670.465 56389.11055.113.6060.6520.2560.777 58481.41074.114.3470.4500.1150.520 59514.91085.115.0970.4290.1370.517 60478.51089.614.2360.5280.1590.582 61498.51127.914.6090.5090.1350.578 63278.71154.414.372 1.0700.599 1.053 64527.41158.715.4180.4460.2450.510 65343.41167.814.266 1.0120.618 1.035 66362.21169.314.2660.8790.4690.879 74346.91039.814.5810.4520.1240.510 77688.91124.215.4200.3900.2900.400 83678.91079.715.4950.4640.1070.504 84624.21042.813.5010.3850.2170.389 88579.61028.914.8250.4560.1620.530 Table5.Basic properties of NGC1798and NGC2192.Parameter NGC1798NGC2192RA(2000)5h11m40s6h15m11sDEC(2000)47◦40′37′′39◦51′1′′l160◦.76173◦.41b4◦.8510◦.64Age1.4±0.2Gyrs1.1±0.1GyrsE(B−V)0.51±0.040.20±0.04[Fe/H]−0.47±0.15dex−0.31±0.15dex(m−M)013.1±0.212.7±0.2distance4.2±0.3kpc3.5±0.3kpcR GC12.5kpc11.9kpcz355pc646pcdiameter10.2pc(8′.3)7.5pc(7′.3)c 1998RAS,MNRAS000,1–8UBV I CCD photometry of two old open clusters NGC1798and NGC21927 Table1.Observing log for NGC1798and NGC2192.Date Target Filter Seeing Telescope Condition96.11.11NCC1798UBV2′′.3SAO a-61cm Photometric97.01.12NCC1798UBV I2′′.2SAO-61cm Non-photometric97.02.12NCC1798,NGC2192UBV I3′′.3SAO-61cm Photometric97.02.13NCC1798,NGC2192UBV I2′′.5SAO-61cm Non-photometric97.03.17NCC2192UBV I2′′.2SAO-61cm Photometric97.10.21NCC1798BV I2′′.3SAO-61cm Non-photometric97.10.24NCC2192BV I2′′.7SAO-61cm Non-photometric97.03.04NCC1798BV2′′.7VBO b-2.3m Non-photometrica Sobaeksan Astronomical Observatoryb Vainu Bappu Observatoryc 1998RAS,MNRAS000,1–8。
群星脉冲星代码
群星脉冲星代码脉冲星（pulsar）是一种自然界中的奇特天体，在1983年获得了诺贝尔物理学奖。
它们是由质量较大的恒星在耗尽核心燃料后，通过引力坍缩而形成的。
脉冲星通常具有非常强烈的磁场，并且以极高的自转速度旋转，所以它们被称为“自转的磁星”。
下面是一个简单的群星脉冲星的代码示例：```pythonimport numpy as np#定义脉冲星类class Pulsar:def __init__(self, mass, radius, period):self.mass = mass # 质量（单位：太阳质量）self.radius = radius # 半径（单位：太阳半径）self.period = period # 自转周期（单位：秒）def calculate_pulse_period(self):#计算脉冲周期（单位：秒）G=6.67某10某某(-11)#万有引力常数c=3某10某某8#光速pulse_period = 2 某 np.pi 某 np.sqrt((self.radius某某3) / (G 某 self.mass)) / creturn pulse_perioddef display_info(self):#打印脉冲星的信息print("Mass: {} M☉".format(self.mass))print("Radius: {} R☉".format(self.radius))print("Period: {} seconds".format(self.period))print("Pulse Period: {}seconds".format(self.calculate_pulse_period())#创建一个脉冲星对象并打印信息pulsar = Pulsar(1.4, 10, 0.01)pulsar.display_info。
全氘甲烷大气示踪剂测量技术
第39卷第5期原子能科学技术Vol.39,No.5 2005年9月Atomic Energy Science and Technology Sep.2005全氘甲烷大气示踪剂测量技术陈绍华,邢丕峰,罗顺忠,羊衍秋(中国工程物理研究院核物理与化学研究所,四川绵阳 621900)摘要:文章通过系统研究甲烷中痕量全氘甲烷的分析方法、分析条件和检测能力等,证明了脉冲进样四极质谱(QMS)法用于大气扩散实验中全氘甲烷的浓度测量是可行的。
同时研究建立了痕量全氘甲烷脉冲进样QMS测量技术。
在50mL/min He流洗、分流比1∶100、电子倍增器检测、慢扫描和SIM模式测量条件下,CD4的最小检测量(体积分数)为511×10-9。
结合天然甲烷与全氘甲烷同位素分离技术,对甲烷中CD4的最小检测量可达到5×10-10,对空气中CD4的探测下限约为815×10-15。
研究表明:气相色谱2四极质谱(GC2QMS)法对CD4的探测下限约是脉冲进样QMS法的40倍。
关键词:测量技术;示踪剂;全氘甲烷;气相色谱2四极质谱法中图分类号:O657163 文献标识码:A 文章编号:100026931(2005)0520409206Measurement T echnologyof Atmospheric T racer Deuterated MethanesC H EN Shao2hua,XIN G Pi2feng,L UO Shun2zhong,YAN G Yan2qiu(I nstitute of N uclear Physics and Chemist ry,China A cadem y of Engineering Physics,P.O.B ox9192220,M iany ang621900,China)Abstract:Analyzing met hods,analyzing condition and testing ability of t race deuterated met hanes CD4in met hane were st udied.The result shows t hat t he measuring met hod of CD4concent ration in at mosp here diff usion test s by p ulse injecting quadrupole mass spec2 t rometer(QMS)met hod is feasible.Pulse injection QMS measuring technology of t race CD4was established.On such condition as rinsing2He gases flow50mL/min,splitting ratio1∶100,electron multiplier testing,slow scanning and SIM modeling,t he lowest detectable limit of CD4in met hane is5.1×bined wit h t he isotope separate technology of CD4and met hane,t he lowest detectable limit of CD4in met hane is5.0×10-10,and for CD4in air,which is8.5×10-15.It is shown t hat t he lowest detectable limit of CD4by gas chromatograp hy2quadrupole mass spect rometer met hod is40times of p ulse injection QMS met hod.K ey w ords:measurement technology;tracer;deuterated methanes;gas chromatography2 quadrupole mass spectro meter收稿日期:2003212222;修回日期:2004212209基金项目:中国工程物理研究院核科学研究基金资助项目(2001077)作者简介:陈绍华(1969—),男,四川安县人,副研究员,博士,分离理论与设备专业 全氘甲烷(12CD4和13CD4)是目前已知的最有发展前途的非放射性远程大气扩散示踪剂之一。