In [1]:
#FeII
FeII_3d6_4s_a6D = {1:('9/2', '0'), 2:('7/2', '384.79'), 3:('5/2', '667.68'),
4:('3/2', '862.61'), 5:('1/2', '977.05')}
#FeII UV1
FeII_3d6_4p_z6Do = {1:('9/2', '38458.99'), 2:('7/2', '38660.05'), 3:('5/2', '38858.97'),
4:('3/2', '39013.22'), 5:('1/2', '39109.32')}
#FeII UV2
FeII_3d6_4p_z6Fo = {1:('11/2', '41968.07'), 2:('9/2', '42114.82'), 3:('7/2', '42237.06'),
4:('5/2', '42334.84'), 5:('3/2', '42401.32'), 6:('1/2', '42439.85')}
#FeII UV3
FeII_3d6_4p_z6Po = {1:('7/2', '42658.24'), 2:('5/2', '43238.61'), 3:('3/2', '43620.98')}
#FeII UV4
FeII_3d6_4p_z4Fo = {1:('9/2', '44232.54'), 2:('7/2', '44753.81'), 3:('5/2', '45079.90'),
4:('3/2', '45289.82')}
#FeII UV5
FeII_3d6_4p_z4Do = {1:('7/2', '44446.91'), 2:('5/2', '44784.79'), 3:('3/2', '45044.19'),
4:('1/2', '45206.47')}
#MgII
MgII_2p6_3s_2S = {1:('1/2', '0')}
#MgII UV1
MgII_2p6_3p_2Po = {1:('1/2', '35669.31'), 2:('3/2', '35760.88')}
#MgI
MgI_2p6_3s2_1S = {1:('0', '0')}
#MgI UV1
MgI_2p6_3s3p_1Po = {1:('1', '35051.26')}
#MnII 3d5(6S)4s a7S
MnII_3d5_4s_a7S = {1:('0', '0')}
#MnII UV1 3d5(6S)4p z7Po
MnII_3d5_4s_z7Po = {1:('2', '38366.23'), 2:('3', '38543.12'), 3:('4', '38806.69')}
#C II 2s22p
CII_2s2_2p_2Po = {1:('1/2', '0'), 2:('3/2', '63.42')}
#C II UV1 2s2p2
CII_2s_2p2_4P = {1:('1/2', '43003.3'), 2:('3/2', '43025.3'), 3:('5/2', '43053.6')}
#O II 2s22p3
OII_2s2_2p3_4So = {1:('3/2', '0')}
#
OII_2s2_2p3_2Do = {1:('5/2', '26810.55'), 2:('3/2', '26830.57')}
OII_2s2_2p3_2Po = {1:('1/2', '40468.01'), 2:('3/2', '40470.00')}
#Ne IV 2322p3
NeIV_2s2_2p3_4So = {1:('3/2', '0')}
#
NeIV_2s2_2p3_2Do = {1:('5/2', '41234.43'), 2:('3/2', '41278.89')}
NeIV_2s2_2p3_2Po = {1:('1/2', '62434.45'), 2:('3/2', '62440.78')}
#Si I 3s23p2
SiI_3s2_3p2_3P = {1:('0', '0'), 2:('1', '77.115'), 3:('2', '223.157')}
# UV1
SiI_3s2_3p4s_3Po = {1:('0', '39683.163'), 2:('1', '39760.285'), 3:('2', '39955.053')}
# UV2
SiI_3s2_3p4s_1Po = {1:('1', '40991.884')}
# UV3
SiI_3s_3p3_3Do = {1:('1', '45276.188'), 2:('1', '45293.629'), 3:('2', '45321.848')}
# Ti II 3d2(3F)4s
TiII_3d2_4s_a4F = {1:('3/2', '0'), 2:('5/2', '94.1142'), 3:('7/2', '225.7039'),
4:('9/2', '393.4459')}
# 3d2(3F)4p z 4Do
TiII_3d2_4p_z4Do = {1:('1/2', '32532.35'), 2:('3/2', '32602.63'), 3:('5/2', '32698.10'),
4:('7/2', '32767.20')}
# 3d2(3F)4p z 2Fo
TiII_3d2_4p_z2Fo = {1:('5/2', '31207.51'), 2:('7/2', '31490.92')}
# 3d2(3F)4p z 4Fo
TiII_3d2_4p_z4Fo = {1:('3/2', '30836.42'), 2:('5/2', '30958.58'), 3:('7/2', '31113.68'),
4:('9/2', '31301.065')}
# 3d2(3F)4p z 4Go
TiII_3d2_4p_z4Go = {1:('5/2', '29544.45'), 2:('7/2', '29734.62'), 3:('9/2', '29968.33'),
4:('11/2', '30240.94')}
# Ne III 2s22p4
NeIII_2s2_2p4_3P = {1:('2', '0'), 2:('1', '642.876'), 3:('0', '920.550')}
#
NeIII_2s2_2p4_1D = {1:('2', '25840.72')}
#CaII
CaII_3p6_4s_2S = {1:('1/2', '0')}
#CaII UV1
CaII_3p6_4p_2Po = {1:('1/2', '25191.51'), 2:('3/2', '25414.40')}
#Transitions
x_energy = 0.18
y_energy_title = 0.05
x_jvalue = 0.91
xline = [0.2, 0.8]
dyline = 0.05
yline0 = 0.07
yline1 = 0.73
x_term = 0.835
x_term1 = x_term+0.045
x_config = 0.9
x_title = 0.5
y_title=1.12
In [4]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 5
e_config = r'3d$^6$4s'
y_config = 0.01
e_term = r'a $^6$'
e_term1 = r'D'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4s_a6D[i+1][0]
energy = FeII_3d6_4s_a6D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Fe II UV1'
n_J_1 = 5
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^6$4p'
y_config = 0.69
e_term = r'z $^6$'
e_term1 = r'D$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4p_z6Do[i+1][0]
energy = FeII_3d6_4p_z6Do[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
fig.savefig('/Users/Benjamin/Dropbox/Zhu_Projects/Fine Structure Emission/Version 1/Movies/FeII_UV1_0.jpg')
xa = 0.09
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2600.17 2.35E8 2.39E-1',
fontsize=14, rotation='vertical', va='bottom')
fig.savefig('/Users/Benjamin/Dropbox/Zhu_Projects/Fine Structure Emission/Version 1/Movies/FeII_UV1_1.jpg')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
fig.savefig('/Users/Benjamin/Dropbox/Zhu_Projects/Fine Structure Emission/Version 1/Movies/FeII_UV1_2.jpg')
ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2626.45 3.52E7 4.55E-2',
fontsize=14, rotation='vertical', va='bottom')
fig.savefig('/Users/Benjamin/Dropbox/Zhu_Projects/Fine Structure Emission/Version 1/Movies/FeII_UV1_3.jpg')
xa = xa+0.23
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02+dyline,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2586.65 8.94E7 7.17E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2612.65 1.20E8 1.22E-1',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2632.11 6.29E7 8.87E-2',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.035
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
head_width=0.01, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
fig.savefig('/Users/Benjamin/Dropbox/Zhu_Projects/Fine Structure Emission/Version 1/Movies/FeII_UV1_4.jpg')
In [278]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 5
e_config = r'3d$^6$4s'
y_config = 0.01
e_term = r'a $^6$'
e_term1 = r'D'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4s_a6D[i+1][0]
energy = FeII_3d6_4s_a6D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Fe II UV2'
n_J_1 = 6
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^6$4p'
y_config = 0.69
e_term = r'z $^6$'
e_term1 = r'F$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4p_z6Fo[i+1][0]
energy = FeII_3d6_4p_z6Fo[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.09
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2382.76 3.13E8 3.20E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2626.45 3.52E7 4.55E-2',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.23
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02+dyline,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2374.46 4.25E7 3.59E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2396.36 2.59E8 2.79E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2632.11 6.29E7 8.87E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.19
#ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.008, yline1+dyline*3, 0., yline0-yline1-dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*3, 0., yline0-yline1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*3, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2374.46 4.25E7 3.59E-2',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2396.36 2.59E8 2.79E-1',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2632.11 6.29E7 8.87E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa-0.045
xa = xa+0.035
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*4, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*3+0.015, r'51.30$\mu{\rm m}$', fontsize=14, ha='left')
# Note: Can't find the data for 42237.06 -> 0. Why???
Out[278]:
In [383]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 5
e_config = r'3d$^6$4s'
y_config = 0.01
e_term = r'a $^6$'
e_term1 = r'D'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4s_a6D[i+1][0]
energy = FeII_3d6_4s_a6D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Fe II UV3'
n_J_1 = 3
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^6$4p'
y_config = 0.69
e_term = r'z $^6$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4p_z6Po[i+1][0]
energy = FeII_3d6_4p_z6Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.09
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline, 0., yline0-yline1+dyline*2+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2344.21 1.73E8 1.14E-1',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2365.55 5.90E7 4.95E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2381.49 3.10E7 3.51E-2',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.23
ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.02, '*2333.51 1.31E8 8.00E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2349.02 1.15E8 9.50E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.035
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*4, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*3+0.015, r'51.30$\mu{\rm m}$', fontsize=14, ha='left')
Out[383]:
In [382]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 5
e_config = r'3d$^6$4s'
y_config = 0.01
e_term = r'a $^6$'
e_term1 = r'D'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4s_a6D[i+1][0]
energy = FeII_3d6_4s_a6D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Fe II UV4'
n_J_1 = 4
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^6$4p'
y_config = 0.69
e_term = r'z $^4$'
e_term1 = r'D$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4p_z4Do[i+1][0]
energy = FeII_3d6_4p_z4Do[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.09
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2249.87 3.00E6 1.82E-3',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2269.52 4.00E5 3.10E-4',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.23
ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.02, '*2252.52 1.00E6 6.00E-4',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2266.70 1.00E6 8.00E-4',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.035
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[382]:
In [381]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 5
e_config = r'3d$^6$4s'
y_config = 0.01
e_term = r'a $^6$'
e_term1 = r'D'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4s_a6D[i+1][0]
energy = FeII_3d6_4s_a6D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Fe II UV5'
n_J_1 = 4
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^6$4p'
y_config = 0.69
e_term = r'z $^4$'
e_term1 = r'F$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = FeII_3d6_4p_z4Fo[i+1][0]
energy = FeII_3d6_4p_z4Fo[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.09
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2260.78 3.18E6 2.44E-3',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2280.62 4.49E6 4.38E-3',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.23
#ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02+dyline,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.01, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2344.21 1.73E8 1.14E-1',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2365.55 5.90E7 4.95E-2',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2632.11 6.29E7 8.87E-2',
# fontsize=14, rotation='vertical', va='bottom')
ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.02, '*2253.82 4.40E6 3.40E-3',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.80E-3',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.035
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
head_width=0.008, head_length=0.01, color='brown')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[381]:
In [380]:
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'2p$^6$3s'
y_config = 0.01
e_term = r' $^2$'
e_term1 = r'S'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = MgII_2p6_3s_2S[i+1][0]
energy = MgII_2p6_3s_2S[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Mg II UV1'
n_J_1 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2p$^6$3p'
y_config = 0.69
e_term = r' $^2$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = MgII_2p6_3p_2Po[i+1][0]
energy = MgII_2p6_3p_2Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.13
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2796.35 2.60E8 6.08E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2280.62 4.49E6 4.38E-3',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.2
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2803.53 2.57E8 3.03E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.03E-3',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[380]:
In [379]:
#MgI
#MgI_2p6_3s2_1S = {1:('0', '0')}
#MgI UV1
#MgI_2p6_3s3p_1P0 = {1:('1', '35051.26')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'2p$^6$3s$^2$'
y_config = 0.01
e_term = r' $^1$'
e_term1 = r'S'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = MgI_2p6_3s2_1S[i+1][0]
energy = MgI_2p6_3s2_1S[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Mg I UV1'
n_J_1 = 1
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2p$^6$3s3p'
y_config = 0.69
e_term = r' $^1$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = MgI_2p6_3s3p_1Po[i+1][0]
energy = MgI_2p6_3s3p_1Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.25
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2796.35 2.60E8 6.08E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2280.62 4.49E6 4.38E-3',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.2
#ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.01, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2803.53 2.57E8 3.03E-1',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.03E-3',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[379]:
In [378]:
#MnII 3d5(6S)4s a7S
#MnII_3d5_4s_a7S = {1:('0', '0')}
#MnII UV1 3d5(6S)4p z7Po
#MnII_3d5_4s_z7Po = {1:('2', '38366.23'), 2:('3', '38543.12'), 3:('4', '38806.69')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'3d$^5$4s'
y_config = 0.01
e_term = r'a $^7$'
e_term1 = r'S'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = MnII_3d5_4s_a7S[i+1][0]
energy = MnII_3d5_4s_a7S[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.005, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Mn II UV1'
n_J_1 = 3
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^5$4p'
y_config = 0.69
e_term = r'z $^7$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = MnII_3d5_4s_z7Po[i+1][0]
energy = MnII_3d5_4s_z7Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.005, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.13
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2576.88 2.80E8 3.58E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2280.62 4.49E6 4.38E-3',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.15
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2594.50 2.76E8 2.79E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.03E-3',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.15
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*2,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*3, 0., yline0-yline1-dyline*2+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2606.46 2.69E8 1.96E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.03E-3',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.015
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[378]:
In [453]:
#C II 2s22p
#CII_2s2_2p_2Po = {1:('1/2', '0'), 2:('3/2', '63.42')}
#C II UV1 2s2p2
#CII_2s_2p2_4P = {1:('1/2', '43003.3'), 2:('3/2', '43025.3'), 3:('5/2', '43053.6')}
#MnII 3d5(6S)4s a7S
#MnII_3d5_4s_a7S = {1:('0', '0')}
#MnII UV1 3d5(6S)4p z7Po
#MnII_3d5_4s_z7Po = {1:('2', '38366.23'), 2:('3', '38543.12'), 3:('4', '38806.69')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 2
e_config = r'2s$^2$2p'
y_config = 0.01
e_term = r' $^2$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = CII_2s2_2p_2Po[i+1][0]
energy = CII_2s2_2p_2Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'C II UVx'
n_J_1 = 3
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2s2p$^2$'
y_config = 0.69
e_term = r' $^4$'
e_term1 = r'P'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = CII_2s_2p2_4P[i+1][0]
energy = CII_2s_2p2_4P[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.09
#ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2325.40 5.99E1 4.86E-8',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
#ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2324.21 1.40E0 2.27E-9',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.17
#ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02-dyline*1,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline*1, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075-0.01, '*2328.83 6.78E1 2.76E-8',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
#ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*0,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline*2, 0., yline0-yline1-dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075-0.01, '*2327.64 8.49E0 6.90E-9',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
#ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*1,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline*3, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075-0.01, '*2326.11 4.43E1 5.40E-8',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.015
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[453]:
In [ ]:
In [461]:
#O II 2s22p3
#OII_2s2_2p3_4So = {1:('3/2', '0')}
#
#OII_2s2_2p3_2Do = {1:('5/2', '26810.55'), 2:('3/2', '26830.57')}
#OII_2s2_2p3_2Po = {1:('3/2', '40468.01'), 2:('1/2', '40470.00')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'2s$^2$2p$^3$'
y_config = 0.01
e_term = r' $^4$'
e_term1 = r'S$^{\rm o}$'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = OII_2s2_2p3_4So[i+1][0]
energy = OII_2s2_2p3_4So[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'O II'
n_J_1 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2s$^2$2p$^3$'
y_config = 0.69
e_term = r' $^2$'
e_term1 = r'D$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = OII_2s2_2p3_2Do[i+1][0]
energy = OII_2s2_2p3_2Do[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
n_J_2 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
#e_config = r'2s$^2$2p$^3$'
#y_config = 0.69
e_term = r' $^2$'
e_term1 = r'P$^{\rm o}$'
y_term = (n_J_1+n_J_2)*dyline+yline1+0.01+0.1
for i in np.arange(n_J_2):
yline = dyline*(i+n_J_1+1)+yline1+0.1
ax.plot(xline, [yline, yline], 'k')
jvalue = OII_2s2_2p3_2Po[i+1][0]
energy = OII_2s2_2p3_2Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.13
ax.arrow(xline[0]+xa, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '3727.10 1.59E-4',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
ax.arrow(xline[0]+xa, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '3729.86 1.98E-6',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.18
ax.arrow(xline[0]+xa, yline1+dyline*(n_J_1+1)+0.1, 0., yline0-yline1-dyline*(2+n_J_1)+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2470.97 2.12E-2',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
ax.arrow(xline[0]+xa, yline1+dyline*(n_J_1+2)+0.1, 0., yline0-yline1-dyline*(n_J_1+1)-0.1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2471.09 5.22E-2',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.015
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[461]:
In [425]:
#O II 2s22p3
#OII_2s2_2p3_4So = {1:('3/2', '0')}
#
#OII_2s2_2p3_2Do = {1:('5/2', '26810.55'), 2:('3/2', '26830.57')}
#OII_2s2_2p3_2Po = {1:('3/2', '40468.01'), 2:('1/2', '40470.00')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'2s$^2$2p$^3$'
y_config = 0.01
e_term = r' $^4$'
e_term1 = r'S$^{\rm o}$'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = NeIV_2s2_2p3_4So[i+1][0]
energy = NeIV_2s2_2p3_4So[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Ne IV'
n_J_1 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2s$^2$2p$^3$'
y_config = 0.69
e_term = r' $^2$'
e_term1 = r'D$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = NeIV_2s2_2p3_2Do[i+1][0]
energy = NeIV_2s2_2p3_2Do[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
n_J_2 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
#e_config = r'2s$^2$2p$^3$'
#y_config = 0.69
e_term = r' $^2$'
e_term1 = r'P$^{\rm o}$'
y_term = (n_J_1+n_J_2)*dyline+yline1+0.01+0.1
for i in np.arange(n_J_2):
yline = dyline*(i+n_J_1+1)+yline1+0.1
ax.plot(xline, [yline, yline], 'k')
jvalue = NeIV_2s2_2p3_2Po[i+1][0]
energy = NeIV_2s2_2p3_2Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.13
ax.arrow(xline[0]+xa, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2425.139 4.0E-4',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
ax.arrow(xline[0]+xa, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '2422.56 2.7E-4',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.18
ax.arrow(xline[0]+xa, yline1+dyline*(n_J_1+1)+0.1, 0., yline0-yline1-dyline*(2+n_J_1)+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '1601.70 5.3E-1',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.07
ax.arrow(xline[0]+xa, yline1+dyline*(n_J_1+2)+0.1, 0., yline0-yline1-dyline*(n_J_1+1)-0.1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '1601.50 1.33E0',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.015
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[425]:
In [460]:
#Si I 3s23p2
#SiI_3s2_3p2_3P = {1:('0', '0'), 2:('1', '77.115'), 3:('2', '223.157')}
# UV1
#SiI_3s2_3p4s_3Po = {1:('0', '39683.163'), 2:('1', '39760.285'), 3:('2', '39955.053')}
# UV2
#SiI_3s2_3p4s_1Po = {1:('1', '40991.884')}
# UV3
#SiI_3s_3p3_3Do = {1:('1', '45276.188'), 2:('1', '45293.629'), 3:('2', '45321.848')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 3
e_config = r'3s$^2$3p$^2$'
y_config = 0.01
e_term = r' $^3$'
e_term1 = r'P'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = SiI_3s2_3p2_3P[i+1][0]
energy = SiI_3s2_3p2_3P[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Si I UV1'
n_J_1 = 3
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3s$^2$3p4s'
y_config = 0.69
e_term = r' $^3$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = SiI_3s2_3p4s_3Po[i+1][0]
energy = SiI_3s2_3p4s_3Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.11
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0+dyline-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '2515.08 7.39E7 2.10E-1',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.23
ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02-dyline,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*1, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline0+dyline*2, 0., yline1-yline0-0.02-dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.07+0.008, yline1+dyline*2, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07+0.07, yline0+dyline*2, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.07+0.07+0.008, yline1+dyline*3, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.02, '*2507.66 5.47E7 8.59E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2519.96 5.49E7 5.22E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2524.87 2.22E8 7.08E-2 ',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.26
#ax.arrow(xline[0]+xa, yline0+dyline*3, 0., yline1-yline0-0.02-dyline*0,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.008, yline1+dyline*3, 0., yline0-yline1-dyline*0+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline0+dyline*3, 0., yline1-yline0-0.02-dyline*1,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.07+0.008, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.02, '*2516.87 1.68E8 1.59E-1 ',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2529.26 9.04E7 5.20E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.01, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*1, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[460]:
In [478]:
# Ti II 3d2(3F)4s
#TiII_3d2_4s_a4F = {1:('3/2', '0'), 2:('5/2', '94.1142'), 3:('7/2', '225.7039'),
# 4:('9/2', '393.4459')}
# 3d2(3F)4p z 4Do
#TiII_3d2_4p_z4Do = {1:('1/2', '32532.35'), 2:('3/2', '32602.63'), 3:('5/2', '32698.10'),
# 4:('7/2', '32767.20')}
# 3d2(3F)4p z 2Fo
#TiII_3d2_4p_z2Fo = {1:('5/2', '31207.51'), 2:('7/2', '31490.92')}
# 3d2(3F)4p z 4Go
#TiII_3d2_4p_z4Go = {1:('5/2', '29544.45'), 2:('7/2', '29734.62'), 3:('9/2', '29968.33'),
# 4:('11/2', '30240.94')}
# Ne III 2s22p4
#NeIII_2s2_2p4_3P = {1:('2', '0'), 2:('1', '642.876'), 3:('0', '920.550')}
#
#NeIII_2s2_2p4_1D = {1:('2', '25840.72')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 4
e_config = r'3d$^2$4s'
y_config = 0.01
e_term = r'a $^4$'
e_term1 = r'F'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4s_a4F[i+1][0]
energy = TiII_3d2_4s_a4F[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Ti II No5'
n_J_1 = 4
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^2$4p'
y_config = 0.69
e_term = r'z $^4$'
e_term1 = r'D$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4p_z4Do[i+1][0]
energy = TiII_3d2_4p_z4Do[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.20
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02-dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*1, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.07+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*2,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.07+0.07+0.008, yline1+dyline*3, 0., yline0-yline1-dyline*2+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '3073.86 1.71E8 1.21E-1',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.035, '3067.23 3.47E7 4.89E-2',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07+0.07-0.03, yline0+dyline*n_J_0+0.035, '3058.28 1.98E6 4.16E-3 ',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*4, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*3+0.015, r'51.30$\mu{\rm m}$', fontsize=14, ha='left')
Out[478]:
In [477]:
# Ti II 3d2(3F)4s
#TiII_3d2_4s_a4F = {1:('3/2', '0'), 2:('5/2', '94.1142'), 3:('7/2', '225.7039'),
# 4:('9/2', '393.4459')}
# 3d2(3F)4p z 4Do
#TiII_3d2_4p_z4Do = {1:('1/2', '32532.35'), 2:('3/2', '32602.63'), 3:('5/2', '32698.10'),
# 4:('7/2', '32767.20')}
# 3d2(3F)4p z 2Fo
#TiII_3d2_4p_z2Fo = {1:('5/2', '31207.51'), 2:('7/2', '31490.92')}
# 3d2(3F)4p z 4Fo
#TiII_3d2_4p_z4Fo = {1:('3/2', '30836.42'), 2:('5/2', '30958.58'), 3:('7/2', '31113.68'),
# 4:('9/2', '31301.065')}
# 3d2(3F)4p z 4Go
#TiII_3d2_4p_z4Go = {1:('5/2', '29544.45'), 2:('7/2', '29734.62'), 3:('9/2', '29968.33'),
# 4:('11/2', '30240.94')}
# Ne III 2s22p4
#NeIII_2s2_2p4_3P = {1:('2', '0'), 2:('1', '642.876'), 3:('0', '920.550')}
#
#NeIII_2s2_2p4_1D = {1:('2', '25840.72')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 4
e_config = r'3d$^2$4s'
y_config = 0.01
e_term = r'a $^4$'
e_term1 = r'F'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4s_a4F[i+1][0]
energy = TiII_3d2_4s_a4F[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Ti II No2'
n_J_1 = 4
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^2$4p'
y_config = 0.69
e_term = r'z $^4$'
e_term1 = r'F$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4p_z4Fo[i+1][0]
energy = TiII_3d2_4p_z4Fo[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.20
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02-dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*1, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.arrow(xline[0]+xa+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.07+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*2,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.07+0.07+0.008, yline1+dyline*3, 0., yline0-yline1-dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '3242.92 1.47E8 2.32E-1',
fontsize=14, rotation='vertical', va='bottom')
ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.035, '3230.12 2.93E7 6.87E-2',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.07-0.03, yline0+dyline*n_J_0+0.035, '3204.36 1.61E6 3.72E-3 ',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*4, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*3+0.015, r'51.30$\mu{\rm m}$', fontsize=14, ha='left')
Out[477]:
In [476]:
# Ti II 3d2(3F)4s
#TiII_3d2_4s_a4F = {1:('3/2', '0'), 2:('5/2', '94.1142'), 3:('7/2', '225.7039'),
# 4:('9/2', '393.4459')}
# 3d2(3F)4p z 4Do
#TiII_3d2_4p_z4Do = {1:('1/2', '32532.35'), 2:('3/2', '32602.63'), 3:('5/2', '32698.10'),
# 4:('7/2', '32767.20')}
# 3d2(3F)4p z 2Fo
#TiII_3d2_4p_z2Fo = {1:('5/2', '31207.51'), 2:('7/2', '31490.92')}
# 3d2(3F)4p z 4Fo
#TiII_3d2_4p_z4Fo = {1:('3/2', '30836.42'), 2:('5/2', '30958.58'), 3:('7/2', '31113.68'),
# 4:('9/2', '31301.065')}
# 3d2(3F)4p z 4Go
#TiII_3d2_4p_z4Go = {1:('5/2', '29544.45'), 2:('7/2', '29734.62'), 3:('9/2', '29968.33'),
# 4:('11/2', '30240.94')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 4
e_config = r'3d$^2$4s'
y_config = 0.01
e_term = r'a $^4$'
e_term1 = r'F'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4s_a4F[i+1][0]
energy = TiII_3d2_4s_a4F[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
title = 'Ti II No1'
n_J_1 = 4
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3d$^2$4p'
y_config = 0.69
e_term = r'z $^4$'
e_term1 = r'G$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = TiII_3d2_4p_z4Go[i+1][0]
energy = TiII_3d2_4p_z4Go[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.20
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02-dyline*0,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*1, 0., yline0-yline1+dyline*0+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.07+0.008, yline1+dyline*2, 0., yline0-yline1-dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.07, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*2,
# head_width=0.01, head_length=0.02, color='black')
#ax.arrow(xline[0]+xa+0.07+0.07+0.008, yline1+dyline*3, 0., yline0-yline1-dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.035, '3384.73 1.39E8 3.58E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.035, '3230.12 2.93E7 6.87E-2',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.07-0.03, yline0+dyline*n_J_0+0.035, '3204.36 1.61E6 3.72E-3 ',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*4, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*3+0.015, r'51.30$\mu{\rm m}$', fontsize=14, ha='left')
Out[476]:
In [481]:
#CaII
#CaII_3p6_4s_2S = {1:('1/2', '0')}
#CaII UV1
#CaII_3p6_4p_2Po = {1:('1/2', '25191.51'), 2:('3/2', '25414.40')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 1
e_config = r'3p$^6$4s'
y_config = 0.01
e_term = r' $^2$'
e_term1 = r'S'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = CaII_3p6_4s_2S[i+1][0]
energy = CaII_3p6_4s_2S[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Ca II No1'
n_J_1 = 2
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'3p$^6$4p'
y_config = 0.69
e_term = r' $^2$'
e_term1 = r'P$^{\rm o}$'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = CaII_3p6_4p_2Po[i+1][0]
energy = CaII_3p6_4p_2Po[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.13
ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline, 0., yline0-yline1+dyline+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '3934.77 1.47E8 6.82E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2280.62 4.49E6 4.38E-3',
# fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.2
ax.arrow(xline[0]+xa, yline0+dyline*1, 0., yline1-yline0-0.02+dyline*1,
head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.008, yline1+dyline*2, 0., yline0-yline1-dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07, yline1+dyline*2, 0., yline0-yline1+dyline*1+0.02,
# head_width=0.01, head_length=0.02, color='green')
#ax.arrow(xline[0]+xa+0.07+0.06, yline1+dyline*2, 0., yline0-yline1+dyline*2+0.02,
# head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '3969.59 1.40E8 3.03E-1',
fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07-0.03, yline0+dyline*n_J_0+0.02, '*2268.29 3.69E6 3.03E-3',
# fontsize=14, rotation='vertical', va='bottom')
#ax.text(xline[0]+xa+0.07+0.06-0.03, yline0+dyline*n_J_0+0.02, '*2359.83 5.00E7 6.30E-2',
# fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.035
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[481]:
In [486]:
# Ne III 2s22p4
#NeIII_2s2_2p4_3P = {1:('2', '0'), 2:('1', '642.876'), 3:('0', '920.550')}
#
#NeIII_2s2_2p4_1D = {1:('2', '25840.72')}
fig = figure(figsize=(10,8))
ax = fig.add_subplot(111)
ax.set_axis_off()
ax.set_xlim(0,1)
ax.set_ylim(0,1.2)
ax.text(x_energy-0.03, y_energy_title, r'E/hc', fontsize=17, ha='center')
ax.text(x_energy-0.03, y_energy_title-0.05, r'(cm$^{-1}$)', fontsize=17, ha='center')
n_J_0 = 3
e_config = r'2s$^2$2p$^4$'
y_config = 0.01
e_term = r' $^3$'
e_term1 = r'P'
y_term = n_J_0*dyline+yline0+0.01
for i in np.arange(n_J_0):
yline = dyline*(i+1)+yline0
ax.plot(xline, [yline, yline], 'k')
jvalue = NeIII_2s2_2p4_3P[i+1][0]
energy = NeIII_2s2_2p4_3P[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
title = 'Ne III'
n_J_1 = 1
#e_config = r'($1{\rm s}^22{\rm s}^22{\rm p}^63{\rm s}^23{\rm p}^6)$3d$^6$4p'
e_config = r'2s$^2$2p$^4$'
y_config = 0.69
e_term = r' $^1$'
e_term1 = r'D'
y_term = n_J_1*dyline+yline1+0.01
for i in np.arange(n_J_1):
yline = dyline*(i+1)+yline1
ax.plot(xline, [yline, yline], 'k')
jvalue = NeIII_2s2_2p4_1D[i+1][0]
energy = NeIII_2s2_2p4_1D[i+1][1]
ax.text(x_energy, yline-dyline/5, energy, fontsize=16, ha='right')
ax.text(x_jvalue, yline-dyline/5, jvalue, fontsize=16, ha='left')
ax.text(x_config, y_config, e_config, fontsize=20, ha='center')
ax.text(x_term+0.025, y_term, e_term, fontsize=22, ha='left')
ax.text(x_term1+0.005, y_term, e_term1, fontsize=22, ha='left')
ax.text(x_title, y_title, title, fontsize=22, weight='bold', ha='center')
#ax.plot([xline[0]+0.1, xline[0]+0.1], [yline0+dyline, yline1+dyline], 'k')
xa = 0.20
#ax.arrow(xline[0]+xa, yline0+dyline, 0., yline1-yline0-0.02,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline, 0., yline0-yline1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075, '3869.86 1.74E1',
fontsize=14, rotation='vertical', va='bottom')
xa = xa+0.17
#ax.arrow(xline[0]+xa, yline0+dyline*2, 0., yline1-yline0-0.02-dyline*1,
# head_width=0.01, head_length=0.02, color='black')
ax.arrow(xline[0]+xa+0.000, yline1+dyline*1, 0., yline0-yline1+dyline*1+0.02,
head_width=0.01, head_length=0.02, color='green')
ax.text(xline[0]+xa-0.03, yline0+dyline*n_J_0+0.075-0.01, '3968.59 5.40E-2',
fontsize=14, rotation='vertical', va='bottom')
#xa = xa+0.015
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*2, 0., dyline-0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.arrow(xline[0]+xa+0.07+0.06, yline0+dyline*3, 0., -dyline+0.01, lw=1,
# head_width=0.008, head_length=0.01, color='brown')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline+0.015, r'25.99$\mu{\rm m}$', fontsize=14, ha='left')
#ax.text(xline[0]+xa+0.07+0.06+0.015, yline0+dyline*2+0.015, r'35.35$\mu{\rm m}$', fontsize=14, ha='left')
Out[486]:
In [ ]: