人生的意义就是思考其意义分享 http://blog.sciencenet.cn/u/qianlivan 天体物理工作者,家乡云南昆明。

博文

学习日志

已有 3747 次阅读 2013-8-29 13:57 |个人分类:知识|系统分类:教学心得| 天文公式

20140919

重置插图计数器的时候

setcounter{figure}{1}

要放在chapter{}

后面,要不然会被重置为0。喵了个咪的!


20140809

ALFALFA的等效积分时间是40秒。


20140806

脉冲星星表

http://www.atnf.csiro.au/research/pulsar/psrcat/

 


20140725

cadcUrlList.txt

下载CGPS数据

wget --content-disposition -i cadcUrlList.txt


20140712

splatalogue的不同分子(能带)是用sid[]=xx来区分的,不同的xx代表不同的分子(能带)。

CI是125

CO(v=0)是204

C$^{17}$O是226


20140628
GBT谱线搜寻
http://www.cv.nrao.edu/~aremijan/PRIMOS/abstract.html

 



20140213

写公式的时候,表达式的一次方和函数的综量容易混淆。
比如
Z(kT/h\nu)
到底是表示kT/h\nu是函数Z(x)的综量x,还是表示Z\times (kT/h\nu)。
对于后一种情况,可以要求加上指数,
Z(kT/h\nu)^1


20130926

中科院邮箱可以定时发邮件,这个可以作为提醒系统。

上边带(USB)和下边带(LSB)的例子:对于98 GHz~102 GHz,如果本振是100 GHz,那么98 GHz到100 GHz是下边带,100 GHz到102 GHz是上边带。



超导隧道结的光子辅助隧穿




20130925
CARMA观测脚本
http://cedarflat.mmarray.org/observing/scripts/


20130922

IDL求条件最小值

原子参数计算




20130919

Kr\v co et al. 2008

Previous studies have shown that the correlation between HI

and $^{13}$CO center velocities is strong (Li & Goldsmith 2003). Hence, the vH parameter is held fixed for each velocity component. The TF estimates obtained from $^{12}$CO do not carry as much confidence (as discussed in the next section); however, the temperature is still kept as a fixed parameter during our fit. For optically thin spectral lines such as HINSA the temperature cannot sufficiently constrain the shape of the absorption feature to be allowed as a free parameter since a wide range of combinations of temperatures and optical depths would be able to fit the data. Thus the gas temperature must be specified in any fits.


Hence themolecular linewidths for each velocity component
are used as guess values for the HI fit. The HI line widths are
allowed to vary in value during the fit near the initial guess values.

Kr\v co & Goldsmith 2010

This technique directly compares all $^{13}$CO emission components to associated HINSA features by determining the $^{13}$CO center velocities, linewidths, and where available, gas temperatures, and then utilizing those parameters as a template with with to fit the HINSA features.


$^{12}$CO emission at the center velocity for each component gives an estimate of the gas temperature for $^{12}$CO under the assumption that the $^{12}$CO emission is optically thick.





20130913

《射电天文工具》P356
大的光学厚度使临界密度减小了一个因子$\tau$(光学厚度),因此,对于$J=1\to 0$的谱线,临界密度的值为$n^*\approx 50 {\rm cm^{-3}}$。于是,能级布居由碰撞决定,所以对于$J=1\to 0$跃迁的激发温度就是$T_k$。


20130910

《射电天文工具》P342
      如果一个能级的布居由与主要碰撞体H$_2$的碰撞决定,H$_2$密度的一级近似可以通过设置碰撞速率$n({\rm H}_2)<\sigma v>$等于$A$系数而得到,碰撞速率中的括号表示假定为玻尔兹曼分布的速度的平均。定义临界密度为使激发温度$T_{\rm ex}$处在辐射温度和热力学温度$T_k$之间的H$_2$密度,记为$n^*$,即$n^*<\sigma v>=A$。

      在有些情况下,碰撞截面$\sigma$通过计算或者实验已知了。作为一级近似,可以假设$\sigma=10^{-16}{\rm cm^2}$。如果速度v取为$1 {\ \rm km/s}$,临界密度为$n_{\rm H_2}=A\times 10^{10} {\rm cm^{-3}}$。

HD 1-0 跃迁的频率为2.7 THz ,偶极矩是(5.85±0.17) × 10-4 Debye,根据爱因斯坦$A$系数公式
$A_{ul}=\frac{64\pi^4}{3hc^3}\nu^3\vert\mu_{ul}\vert^2$
可以得到$A\sim 7.85\times 10^{-8}$,所以临界密度大约是800 cm$^{-3}$。



20130906

学到一个好用的小工具kpvslice,可以画任意方向的pv图。



20130902

用频率2 MHz的雷达可以探测火星表面下5千米的地层。

行星数据库
JMars
http://jmars.asu.edu/
http://anserver1.eprsl.wustl.edu/
http://www.planetserver.eu/


20130829
《射电天文工具》P206
     电离气体云的辐射,吸收系数

$\kappa_{\nu}=\frac{4}{3}\frac{Z^2 e^6}{c}\frac{N_i N_e}{\nu^2}\sqrt{\frac{1}{2\pi (mkT)^3}}\ln\frac{p_2}{p_1}$

连续谱光深

$\tau_{\nu}=8.235\times 10^{-2}\left(\frac{T_e}{\rm K}\right)^{-1.35}\left(\frac{\nu}{\rm GHz}\right)^{-2.1}\left(\frac{\rm EM}{\rm pc\ cm^{-6}}\right)$

其中辐射量
$\frac{\rm EM}{\rm pc\ cm^{-6}}=\int^{s/\rm pc}_{0}\left(\frac{N_e}{\rm cm^{-3}}\right)^2 {\rm d} \left(\frac{s}{\rm pc}\right)$

《射电天文工具》P303
      复合线光深
$\tau_{\nu}=1.92\times 10^{-3}\left(\frac{T_e}{\rm K}\right)^{-5/2}\left(\frac{\Delta\nu}{\rm GHz}\right)^{-1}\left(\frac{\rm EM}{\rm pc\ cm^{-6}}\right)$

      谱线强度和连续谱强度比

$\frac{T_{\rm L}}{T_c}\left(\frac{\Delta v}{\rm km\cdot s^{-1}}\right)=\frac{6.985\times 10^3}{a(\nu,T_e)}\left[\frac{\nu}{\rm GHz}\right]^{1.1}\left(\frac{T_e}{\rm K}\right)^{-1.15}\frac{1}{1+N({\rm He^+})/N({\rm H^+})}$



$1.92\times 10^{-3}/8.235\times 10^{-2}\times 3\times 10^5=6.99\times 10^3$





https://blog.sciencenet.cn/blog-117333-720645.html

上一篇:旅行杂记(九)从昆明北到王家营
下一篇:观测工作日志
收藏 IP: 159.226.171.*| 热度|

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...
扫一扫,分享此博文

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-5-18 01:19

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部