电子掺杂，空穴掺杂与化学势的关系
在研究超导问题的时候，不可避免的总会遇到体系时电子掺杂还是空穴掺杂，因为这两种不同的掺杂所对应的性质是不同的，最近重新温习一下超导的相关知识，正好也利用一个具体的紧束缚近似模型来说明一下电子和空穴掺杂到底是怎么回事，以及如何的通过自洽方法，通过调整化学势来决定体系到底是电子掺杂还是空穴掺杂{:.info} 半满实空间分析在这里先采用李正中<<固体理论>>第十一章第九节的Hubbard模型来说明一下填充到底如何确定 H=-t\sum_{,\sigma}(C_{i\sigma}^\dagger...
年轻人的旅行,坐标珠海
学习时间久了，就用一场说走就走的旅行来调节一下，所以在几个小时的时间从安排到订票，完成了所有的事项，然后在第二天一大早开启了珠海之旅。{:.success} {:.circle.border.shadow,:width=”600px”,:height=”400px”} 长隆 第一天的行程安排自然就是珠海长隆，朋友找的住宿的地方还是非常棒的，而且离长隆也不是很远，几分钟就可以到达。我因为感觉自己的身体好像并不能玩过山车那么刺激的东西，所以白天的娱乐项目，我尽量就少玩了几个，不过还是去先去升降机做了一下实验，结果就是在失重的一瞬间，整个人不能呼吸，身体还会发抖，跟坐飞机的时候爬升那一段时间完全相同，看来以后和这类项目没有缘分了。在这里留一张鱼的照片吧，娱乐项目我可还没大神到能自拍 白鲸剧场 除了这些娱乐项目之外，有趣的就是白鲸剧场和海豚剧场了。之前去西安的时候，就因为时间原因错过了这些，终于在珠海长隆看到了这些。不得不惊叹于这些动物的智商，在饲养员的训练之下，居然可以完成很多有趣的项目。这里附上照片一张作为留念 {:width=”600px”,:height=”400px”} 烟火 接近300的门票，总算是在最后结束的时候，让我有了一种物超所值的感觉。去年9月份去长沙的时候，总以为橘子洲会有烟花表演，原来因为环境问题早就不是每次去都看得到了，所以也算是留下了点遗憾，虽然之后因为国际交流的原因橘子洲上演了盛大的烟花表演，可自己也没有过去。长隆最后的保留节目就是烟花表演。好多年前，过年的时候在自己所在的城市，在正月十五的时候也是可以看到烟花的，自从上了大学之后，每次这个时间也都在学校了，一直都错过了这些，好在这一次我不会再错过了。烟花表演加上酷炫的无人机表演，配上黑人抬棺的Dj，给长隆一天的游玩画上了句号，也让我觉得不虚此行。完全被烟花表演圈粉了，计划下一次去上海迪士尼看烟花表演，先预定一下计划。 {:.circle.border.shadow,:width=”600px”,:height=”400px”} 海 ...
构建三角形或者平行四边形点阵
平时在做紧束缚模型的时候，都是在n*n的点阵上进行的，但是有时候可能也需要在三角形或者平行四边形样式的点阵上去计算一些性质，正好趁手头空闲就把这个做了一下，还是非常的简单。{:.info} 正方+三角12345678910111213141516171819202122232425262728program e1implicit noneinteger knreal pi,dkcomplex imparameter(kn = 10,pi = 3.1415926535,im = (0.0,1.0))integer m1,m2,m3real t1,t2complex re1,re2call cpu_time(t1) ! 获取当前系统时间open(12,file="squ.dat")open(13,file="tri.dat")do m1 = 1,kn do m2 = m1,m1+10 write(12,*)m1,m2,sqrt(1.0*m1+m2) end doend do!----------------------------do m1 = 1,kn do m2 = 1,m1 write(13,*)m1,m2,sqrt(1.0*m1+m2) end doend doclose(12)close(13)call cpu_time(t2)stopend program e1 {:width=”330px”,:height=”495px”}{:width=”330px”,:height=”495px”} 六边形区域123456789101112131415161718192021222324252627282930313233343536373839! Author:YuXuanLi! E-Mail:yxli406@gmail.com! Graphene BZ plot module pub implicit none complex,parameter::im = (0.0,1.0) real,parameter::pi = 3.14159265359 integer n0 ! n0 边形 real r0 ! 外接圆半径 end module...
超大型系数矩阵对角化
最近在进行一些三维实空间中系统性质的研究，利用紧束缚束缚模型，当三个方向都取开边界时，矩阵会变的非常大，组里的服务器完全不能计算，所以只好寻找一些替代的方法。首先哈密顿量是个大型系数矩阵，这个时候可以通过lanczos 算法将对称的厄密矩阵变成三对角矩阵，这个时候矩阵的本征值求解就变的比较容易了。{:.info} 算法介绍123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127#!/usr/bin/env python# -*- coding: utf-8 -*- # import useful modulesimport matplotlib from math import factorialfrom numpy import *from pylab import *from numpy.polynomial.hermite import * # use LaTeX, choose nice some looking fonts and tweak some settingsmatplotlib.rc('font', family='serif')matplotlib.rc('font', size=16)matplotlib.rc('legend', fontsize=16)matplotlib.rc('legend', numpoints=1)matplotlib.rc('legend', handlelength=1.5)matplotlib.rc('legend',...
Linux中批量执行编译并运行Fortran
最近因为要大量重复的跑一些程序，而且只是参数的小修，所以干脆花点时间整理一个界几个shell脚本，来自动的完成程序的编译及执行。{:.info} 批量编译Fortran并运行12345678910#!/bin/sh #============ get the file name =========== Folder="/home/yxli/te" #要批量编译哪个文件夹下面的Fortranfor file_name in ${Folder}/*.f90do temp_file=`basename $file_name .f90` ifort -mkl $file_name -o $temp_file.out ./$temp_file.out & # 编译成功之后自动运行donerm *out # 删除编译后文件 递归的读取指定文件夹下面的所有Fortran文件并编译执行1234567891011121314151617181920212223242526#!/bin/bash function getdir(){ for element in `ls $1` do dir_or_file=$1"/"$element if [ -d $dir_or_file ] then getdir $dir_or_file else # 下面的全是文件 if [ "${dir_or_file##*.}"x = "f90"x ]||[ "${dir_or_file##*.}"x = "f"x ];then # 筛选处特定后缀的文件 dir_name=`dirname $dir_or_file` # 读取目录 file_name=`basename $dir_or_file .f90` # 读取以f90结尾的文件名 out_file_name="$dir_name/$file_name" # 定义编号成功的文件名称 ifort -mkl...
