Zhaohua Tian's Web

之前写文章一直不注重逻辑，以为做科研就是会推导，会计算。当自己开始写文章之后，才发现，做科研绝不仅仅只是推导和计算，还要有物理思想，也要在算之前学会分析某个问题是否有意义。写文章的过程，是完整的向读者展现自己的科研成果，而不是为了炫技。无论是在题目、摘要，还是正文，目的都是要让别人以最小的代价，知道你做了什么，你想强调什么，你为什么做的牛逼。最近一边写文章，一边被老板吐槽，幡然醒悟自己以前的一些科研习惯是多么的不好。现在对一些需要注意的点进行总结。 目前对写文章的一些总结： 摘要部分 简洁、通俗易懂，少用一些不常见的名词来描述你的工作，这样别人不知道你做的是什么。 第一句一般是背景，也可以是你做的东西。 要有条理，有层次，每一句话都是有意义的。 最后一句要说意义。 介绍部分 基本框架可以写成： 大背景是啥 问题是啥 别人怎么解决 解决的不好在哪？ 然后我们xx 不要天花乱坠的说别人的工作，一定要注意和自己的工作进行联系和对话，不要跑偏 介绍我们的工作的时候，可以和摘要类似，说法上稍微有点变化就行。正文部分 即使是理论文章，公式也不能太多，将最关键的公式写成行间公式，普通公式就写成 ...

教科书结果重复继续使用MEEP来作计算，这次将先会重复一下教科书的结果“Photonic Crystals: Molding the Flow of Light “：P128 Figure5 上图是要计算的结构，是一个纳米梁波导，其基本参数为: 周期设置为a=1 圆孔半径设置为r=0.25 宽度设置为w=1 高度设置为h=0.4 所用材料折射率为:3.4641 下面是MEEP的具体实现，首先是导入相应的包和设置基本参数 12345678910111213141516171819#### import meep as mpfrom meep import mpbresolution = 25 # pixels/a 分辨率设置# 三维结构设置# z周期a = 1 # units of um# 圆孔半径r = 0.25 # units of um# 高度h = 0.4 # units of um# 纳米梁宽度w = 1 # units of um# 归一化处理r = r/a # units of "a"h ...

Quantum #Optics学习知识应该细嚼慢咽，进行深度思考与优化。 真正的大师永远怀着一颗学徒的心。 每次涉及到新的领域的时候，就会发现，自己的一些基础概念不过关，这些都需要自己去仔细的理解、消化、思考。 以这个文件为目录，我将会对一些经典的重要的知识点进行一个非常详尽的说明与概括，并写好相应的参考文献。 量子光学 与量子电动力学的联系与区别 真空自发辐射模型 Lamb位移 无限深势的本征态 简谐振子的本征态 压缩态、相干态、Fock态、热态的定义、光子数分布 量子光场的时空间分布的物理意义的理解，以及与经典光场分布的联系 量子拍频效应 非线性光学产生压缩态的具体方程 多光子、多端口分束器原理 Bell不等式 薛定谔方程方法 主方程方法 输入、输出方法 量子信息 发展历史 量子计算的模型 几种量子算法 量子模拟 量子随机行走 量子隐形传态 拓扑与拓扑光子学 发展历史 量子霍尔效应、分数量子霍尔效应 获得诺贝尔奖的理论？ 一维SSH模型 二维SSH模型 光学拓扑的理论分类 PT-Symmetry系统 发展历史 经典文献重复 光学PT示例 量子PT文献 Waveguide-Q ...

写作一直是自己的弱项，现在想要进行一个非常系统化的训练，不然文章完全写不出来。根据我自己的感受，写好文章需要具备的素质是： 能够自如的介绍清楚技术方面的内容，比如能够介绍清楚方程的定义、推导以及求解等，能够说清楚我们解决问题的思路以及结果。 能够系统的总结和概括自己文章的主要内容、创新点 能够对研究的动机、背景做出非常系统、简洁的描述 相关语句比较流畅，用词精准 这些做起来不是一蹴而就，需要多看、多学、多训练。我在这里对一些非常重要的短句、词汇以及长句表达做一些收集，便于以后写文章参考。 重要词汇regulation ：管理，控制enduring： 持久的,不朽的exclusively： 专门地;专有地unprecedented： 前所未有的, 无前例的paradigmatic： 范例的a plethora of phenomena： 种类繁多的现象tailoring the emission： 调整发射 典型专业表述收集 Quantum optics is an exciting field, in which many fundamental experiments, rev ...

MathematicaI have been using Mathemtica to write notes. It’s really convenient to use notebook files because we can write program, formular, figures and text in one file. However, not everyone uses Mathematica, If we want to share our notes with others, exporting nb files into pdf format is necessary. In default, the generated pdf by Mathematica is relatively simple, and there is no TOC or bookmarks in the generated pdf files. It will be very inconvenient to read if the notes are very long. I ...

Meep说明回归科研以后，一些爱好需要重新拾取起来了。最近学会了用COMSOL计算光子晶体能带，现在尝试来用Meep计算光子晶体能带，希望可以以后有机会作为主力的程序来使用。 我主要参考的是官方文档：https://mpb.readthedocs.io/en/latest/Python_Tutorial/ 简单二维光子晶体能带例子能带计算本次主要想重复书籍Photonic Crystals: Molding the Flow of Light second edition 第五章的结构的能带计算：按照官方文档的介绍，直接设置结构，计算相应能带即可，具体可以看代码和注释123456789101112131415161718192021222324252627282930313233import mathimport meep as mpfrom meep import mpbnum_bands = 8# 定义好扫描的k波矢k_points = [mp.Vector3(), # Gamma mp.Vector3(0.5), # X ...

物理差分化表示天下武功出少林，物理中经典光学问题都全部用麦克斯韦方程组描述，量子力学的关键则是薛定谔方程。很多经典的量子物理模型如无限深势阱、简谐振子等模型，薛定谔方程的解都可以解析的表示，但是大部分物理系统得用数值方法来计算。 本篇笔记将会总结如何通过数值方法分析特定物理系统的本征值以及本征矢量。主要参考了该回答：Numerical solution to Schrödinger equation - eigenvalues 一维空间中一般薛定谔方程可以表示为： i \hbar \frac{\partial}{\partial t} \Psi(x, t)=\left[-\frac{\hbar^2}{2 m} \frac{\partial^2}{\partial x^2}+V(x, t)\right] \Psi(x, t)我们假设势能项不含时间，考虑稳态的薛定谔方程，并且近似认为，此时相应的方程变为 V(x) \Psi (x)- \psi''(x)=E \Psi (x)采用差分方法，格点化 为 \Psi (x_0),\Psi (x_1),\Psi (x_2),\text{.. ...

Mathematica目的：绘制给定直角坐标系下三个垂直的三维箭头，需要绘制出三维效果。 MATLAB没有天然的三维箭头绘制支持，需要用到别人写好的包。我经过多次尝试发现这些包的实现结果都不如人意。最后不得不采用Mathematica来实现，实现的过程中也发现有许多的坑，本篇笔记将会对实现方法进行记录。 角度的对比用Mathematica绘制时，12345678910111213141516171819202122232425262728293031323334Table[ viewpt = {Sin[\[Theta]view]*Cos[\[Phi]view], Sin[\[Theta]view]*Sin[\[Phi]view], Cos[\[Theta]view]}; Module[{ thickline = 0.02, thickhead = 0.05}, {d1, d2, d3} = Dipole3DfromRot[{1, 1, 0, 0, 0}]; p2 = Graphics3D[{(*Text[Style["d1",18,Red],d1] ...

绘图本页面将会将常用的MATLAB和Mathematica命令总结出来，便于快速使用 MATLABMATLAB如何绘图而不显示？123f = figure('visible','off');plot(a)saveas(f,'newout','fig') MATLAB如何设置图片宽度和大小1set(gcf,'units','normalzied','position',[x,y,width,height]) MATLAB绘图使用全屏1f=figure('visible','on','units','normalized','outerposition',[0 0 1 1]); MATLAB如何使用-tight-subplot1234567891011%% Directly show figuresfig1=figure(1);set(gcf,'units','normalized','Position', [0.1 0.1 0.6 0.8])[ha, pos] = tight_subplot(4,4,[.05 .01],[.1 .1],[.05 .05]);for l=1:16 a ...

MATLAB #DataprocessingGIF现如今广泛流行于网络平台，因为其比图片具有更多的信息，同时又比视频具有更小的体积。 在科研、工作中，经常会需要将一个过程动态的展示出来，之前都是登录一些在线网站，上传生成的图片，没有网络的话就不方便，同时也不利于信息安全。 我最近发现MATLAB直接生成GIF也是非常方便的，写一篇笔记记录。 方法1：直接绘图并生成GIF第一种方法是直接生成GIF，实现的方法可以通过直接在MATLAB搜索GIF关键词来找到。需要首先通过代码绘制一系列图片，并且将这些图存储起来，然后通过代码将图片导入到一个GIF文件中。我这里直接给出代码：12345678910111213141516171819202122232425262728293031323334% To generate a gif files from plotnumx=100;numt=20;x=linspace(0,6*pi,numx);t=linspace(0,10,numt);%% 1. Directly export gif figures needed%-------------- ...