基于第一性原理量子输运模拟的并行计算

为了解决基于第一性原理分析计算大尺度量子输运体系时遇到的耗时长久问题,挖掘密度泛函理论与非平衡格林函数相结合方法(DFT+NEGF方法)在自洽迭代过程中的计算热点,就计算电子密度矩阵时的能量点积分和计算格林函数时的矩阵求逆/乘法运算提出MPI/Open MP并行计算方案.能量点积分采用MPI多进程并行方案,在数据初始化时需要将稀疏矩阵和积分能量点依照轮询调度算法分配给各进程.矩阵求逆/乘法的并行化既可调用ScaLAPACK子程序实现又可调用IntelMKL数学库中的OpenMP多线程加速函数实现.由于不同能量点计算的独立性,能量点积分采用的MPI并行计算获得近乎线性的加速比曲线.由于Open MP多线程并行采用的是基于共享内存的数据交换机制以及线程间切换通信开销小,矩阵求逆/乘法运算的OpenMP并行实现在计算效率上要优于而在程序的可扩展性上要劣于MPI多进程并行实现.     

团簇LaO的理论研究

采用密度泛函理论研究LaO团簇体系。中性分子LaO的基态是两重态(2Σ),阴离子LaO-和阳离子LaO 的基态都是单重态(1Σ)。使用不同的方法计算团簇LaO的电子亲和能和电离能。计算结果表明用BLYP方法和弥散极化基组计算结果和实验数据吻合较好。用含时密度泛函理论计算团簇LaO的低能激发态,从理论上归属LaO-的光电子能谱的谱峰和LaO的吸收光谱的谱峰。计算得到与实验一致的结果。     

电子激发和溶剂效应对芴酮-甲醇分子间氢键增强现象的理论研究 (cited: 1)

理论研究了电子激发和溶剂效应导致的芴酮-甲醇复合体系中分子间氢键增强现象．通过基态和激发态性质的计算,不仅展示了分子间氢键键长的变化以及变化在振动光谱中的影响,而且揭示了导致氢键变化的内在物理机制：溶质分子的电子激发及溶剂化效应引起的电子重新分布,增大了溶质和溶剂分子的偶极矩,导致了它们之间的相互作用的增大,并最终加强了分子间氢键的强度．还分别对处于液相及气相中的复合体的基态和激发态的几何结构、红外谱、复合体及构成分子的偶极矩进行了理论计算,结果阐明了电子激发与溶剂化效应对氢键变化的贡献,同时还发现只有进一步引入溶剂化效应,复合体的基态、激发态的性质才能与实验达到精确一致．所有激发态均采用所开发的基于含时密度泛函理论解析计算一阶、二阶激发态能量导数的方法．     

遗传算法结合密度泛函理论研究Cu13团簇的几何结构

基于Gupta 多体势采用遗传算法详细计算了金属团簇Cu13的几何结构,对所得到的大量低能（基态及低激发态）结构,利用密度泛函理论方法作进一步优化计算。结果表明：尽管遗传算法所得基于Gupta原子间多体势的基态结构（Ih高对称性的紧致结构）并不对应第一性原理计算结果（非紧致低对称性基态）,但遗传算法给予的大量候选结构经密度泛函理论再次计算仍然可以高效地得到真实基态结构,体现出遗传算法在计算具有非紧致低对称性基态体系时的有效性。     

遗传算法结合密度泛函理论研究Cu13团簇的几何结构

基于Gupta 多体势采用遗传算法详细计算了金属团簇Cu13的几何结构,对所得到的大量低能（基态及低激发态）结构,利用密度泛函理论方法作进一步优化计算。结果表明：尽管遗传算法所得基于Gupta原子间多体势的基态结构（Ih高对称性的紧致结构）并不对应第一性原理计算结果（非紧致低对称性基态）,但遗传算法给予的大量候选结构经密度泛函理论再次计算仍然可以高效地得到真实基态结构,体现出遗传算法在计算具有非紧致低对称性基态体系时的有效性。     

KSSOLV-GPU：一款利用GPU高效求解Kohn-Sham方程的平面波基组密度泛函理论MATLAB程序包

KSSOLV(Kohn-Sham Solver)是一款用于求解平面波基组下Kohn-Sham方程(KS-DFT)的MATLAB(Matrix Laboratory)工具箱. 在KS-DFT的基态计算中,通常自洽场迭代中Kohn-Sham哈密顿量的对角化是最昂贵的部分. 为了使得个人计算机也能够执行数百个原子的中等大小KS-DFT计算,本文提出了一种CPU-GPU的混合编程方案,通过调用MATLAB内置的并行计算工具箱来加速在KSSOLV中实现的迭代对角化算法. 比较了KSSOLV-GPU在RTX3090、V100、A100三种GPU上的性能;结果表明,对于包含128个原子的块状硅体系,与串行的CPU计算相比,混合CPU-GPU的编程可以实现约10倍的加速. 特别是其在最新的民用GPU显卡RTX3090上也具有优秀的表现,可以预想到在不远的将来,KSSOLV-GPU借助MATLAB强大的可视化能力与GPU的加速支持可以在一台配备了民用GPU显卡的个人电脑上实现常规的DFT计算分析与可视化,从而降低了材料模拟与计算领域的门槛.     

Density Functional Study on Relative Energies, Structures, and Bonding of Low-lying Electronic States of Lutetium Dimer

用密度泛函理论方法研究了镥二聚体(Lu₂)低能量电子态的性质,计算了电子态相对能量、平衡键长、振动频率以及基态解离能,考察了密度泛函性质、相对论有效势种类以及Hartree-Fock交换作用大小对计算结果的影响．结果表明,无论采用何种密度泛函和相对论有效势,体系的基态都为三重态,与其他一些基于分子轨道理论的从头计算方法得到的结论是一致的．另外,与分子轨道从头计算结果以及实验结果比较发现,采用杂化密度泛函理论和Stuttgart小核有效势计算得到的结果总体吻合最好．最后,特别分析研究了B3LYP计算中Hartree-Fock交换作用大小对基态键长和基态解离能的影响,发现随着交换作用的增大,键长增长,解离能减小,这是由于5d轨道杂化导致的共价成键作用减弱造成的．     

苋菜红分子基态和激发态结构与光谱性质的量子化学研究

应用英国Edinburgh FLS920P光谱仪对苋菜红的吸收光谱和荧光光谱进行了实验检测．同时,分别采用密度泛函理论（DFT）和含时密度泛函理论（TD-DFT）对苋菜红分子的基态和激发态构型进行优化,经振动分析合理后,比较这两种能态下分子结构的差异,并对其前线分子轨道和发光机制进行了分析．在此基础上,选用6种泛函并结合溶剂化模型（PCM）在6-311++G(d, p)水平上分别计算苋菜红的吸收光谱和荧光光谱．计算结果表明：苋菜红含有分子内氢键,基态结构非平面,两个萘环所在平面有一定的夹角,激发态时两个萘环共平面;CAM-B3LYP泛函得到光谱的理论值与实验结果基本吻合;421 nm处的荧光峰值波长对应的轨道跃迁为LUMO→HOMO-1．     

三碘甲状腺素团簇结构与光谱性质的密度泛函理论研究

我们利用密度泛函理论(DFT)，在B3LYP/Lan12mb基组水平上，得到了三碘甲状腺素团簇的几何和电子结构.在此基础上，利用含时密度泛函理论(TDDFT)，使用相同的基组和采用极化连续介质模型(PCM)，对其溶剂效应下的吸收光谱进行研究．研究结果表明，优化所得三碘甲状腺素团簇的对称性为C₁;在基态稳定结构基础上，研究了该分子的红外和拉曼分子振动谱特性，同时研究了其输运性质，即三碘甲状腺素团簇为p型输运材料;通过含时密度泛函理论，在优化好的基态结构基础上，又计算了它的溶剂效应，进一步得出该分子在水溶剂中的吸收光谱特性.     

丙酮激发态α解离的理论研究

利用密度泛函理论对丙酮分子基态和最低三重激发态的结构进行了构型优化和振动频率计算,并对丙酮三重激发态α解离反应进行了理论分析。     

高激发态原子的相干效应

本文以高激发态原子为研究对象, 由基态、激发态和高激发态能级形成阶梯型三能级系统, 理论上求解阶梯型三能级系统的密度矩阵方程, 研究了高激发态原子的相干效应, 计算获得探测光的吸收和色散曲线. 并研究了高激发态原子间相互作用以及外加电场对相干效应的影响. 结果表明, 外加场可以使吸收和色散曲线产生频移.     

杀草强分子的拉曼光谱的密度泛函理论研究

本文采用密度泛函理论(density functional theory,DFT),在B3LYP杂化泛函,6-31++g(d,p)(C,H,N)和LanL2DZ(Ag)基组下对杀草强分子及其Ag复合物的结构进行了优化;通过频率计算,获得了杀草强分子及其Ag复合物的拉曼光谱,并利用势能函数分布(PED)对拉曼光谱进行了指认,结合SERS光谱推测了杀草强和增强基底之间的吸附方式;采用含时密度泛函理论(time dependent density functional theory,TDDFT)对杀草强分子和杀草强分子-Ag复合物进行了激发态的分析计算。     

氢键效应对2,3-二氢-3-酮基-1H-吡啶并[3,2,1-kl]吩噻嗪光物理性质的影响

利用含时密度泛函理论(TDDFT)方法和飞秒时间分辨的瞬态吸收光谱技术对2,3-二氢-3-酮基-1H-吡啶并[3,2,1-kl]吩噻嗪(PTZ4)和3-酮基-1H-吡啶并[3,2,1-kl]吩噻嗪(PTZ5)这两种荧光探针分子的光物理性质进行了研究. TDDFT结果表明PTZ4和PTZ5在甲醇溶液形成了氢键络合物导致它们吸收峰的红移. PTZ4分子在基态有四种稳定构型,其在四氢呋喃溶液中的双荧光峰正是来自于四种构型下的内部电荷转移态. PTZ4分子在四氢呋喃和甲醇溶液中的瞬态吸收光谱表明,从局域态到转移态的弛豫时间常数在四氢呋喃中为16.0 ps,在甲醇中为7.5 ps;PTZ4分子在甲醇中的激发态寿命为53.8 ps,而这种超短的寿命可能是由于PTZ4分子在激发态时形成的面外型氢键络合物导致的.     

Steady state quantum discord for circularly accelerated atoms

We study, in the framework of open quantum systems, the dynamics of quantum entanglement and quantum discord of two mutually independent circularly accelerated two-level atoms in interaction with a bath of fluctuating massless scalar fields in the Minkowski vacuum. We assume that the two atoms rotate synchronically with their separation perpendicular to the rotating plane. The time evolution of the quantum entanglement and quantum discord of the two-atom system is investigated. For a maximally entangled initial state, the entanglement measured by concurrence diminishes to zero within a finite time, while the quantum discord can either decrease monotonically to an asymptotic value or diminish to zero at first and then followed by a revival depending on whether the initial state is antisymmetric or symmetric. When both of the two atoms are initially excited, the generation of quantum entanglement shows a delayed feature, while quantum discord is created immediately. Remarkably, the quantum discord for such a circularly accelerated two-atom system takes a nonvanishing value in the steady state, and this is distinct from what happens in both the linear acceleration case and the case of static atoms immersed in a thermal bath.     

Iterative submatrix diagonalisation for large configuration interaction problems

The Davidson method has been highly successful for solving for eigenpairs of the large matrices that are common in quantum chemical simulations. Electronic structure simulations, however, can still easily generate matrices that are too large for current computational resources to handle. Therefore, many strategies have arisen to obtain eigenpairs of sufficient accuracy without considering the full Hamiltonian matrix. This article introduces one such strategy by creating a systematic series of submatrix approximations to the full matrix using natural orbitals. By solving for eigenpairs in this series, the eigenvalue accuracy can be gradually increased until a convergence threshold is reached. Importantly, this allows the series to terminate without ever reaching the full matrix, resulting in lower computational costs and reduced memory demands. Application of the method to the full configuration interaction problem for ground states, excited states and potential energy scans of various systems shows that the iterative submatrix diagonalisation method can systematically control eigenvalue errors and provide substantial cost-savings. This method is therefore expected to be highly useful for large-scale diagonalisation problems in electronic structure theory.     

Results of numerically solving an integral equation for a two-fermion system

A two-particle system is described by integral equations whose kernels are dependent on the total energy of the system. Such equations can be reduced to an eigenvalue problem featuring an eigenvalue-dependent operator. This nonlinear eigenvalue problem is solved by means of an iterative procedure developed by the present authors. The energy spectra of a two-fermion system formed by particles of identical masses are obtained for two cases, that where the total spin of the system is equal to zero and that where the total spin of the system is equal to unity. The splitting of the ground-state levels of positronium and dimuonium, the frequency of the transition from the ground state of orthopositronium to its first excited state, and the probabilities of parapositronium and paradimuonium decays are computed. The results obtained in this way are found to be in good agreement with experimental data.     

Charge transfer, double and bond-breaking excitations with time-dependent density matrix functional theory

Time-dependent density functional theory (TDDFT) in its current adiabatic implementations exhibits three remarkable failures: (a) completely wrong behavior of the excited state surface along a bond-breaking coordinate; (b) lack of doubly excited configurations; (c) much too low charge transfer excitation energies. These TDDFT failure cases are all strikingly exhibited by prototype two-electron systems such as dissociating H2 and HeH+. We find for these systems with time-dependent density matrix functional theory that: (a) Within previously formulated simple adiabatic approximations, the bonding-to-antibonding excited state surface as well as charge transfer excitations are described without problems, but not the double excitations; (b) An adiabatic approximation is formulated in which also the double excitations are fully accounted for.     

Time-dependent density-functional theory in massively parallel computer architectures: the OCTOPUS project

Octopus is a general-purpose density-functional theory (DFT) code, with a particular emphasis on the time-dependent version of DFT (TDDFT). In this paper we present the ongoing efforts to achieve the parallelization of octopus. We focus on the real-time variant of TDDFT, where the time-dependent Kohn-Sham equations are directly propagated in time. This approach has great potential for execution in massively parallel systems such as modern supercomputers with thousands of processors and graphics processing units (GPUs). For harvesting the potential of conventional supercomputers, the main strategy is a multi-level parallelization scheme that combines the inherent scalability of real-time TDDFT with a real-space grid domain-partitioning approach. A scalable Poisson solver is critical for the efficiency of this scheme. For GPUs, we show how using blocks of Kohn-Sham states provides the required level of data parallelism and that this strategy is also applicable for code optimization on standard processors. Our results show that real-time TDDFT, as implemented in octopus, can be the method of choice for studying the excited states of large molecular systems in modern parallel architectures.     

The single photon superradiance from the eigenmode analysis

Using the eigenmode analysis of the scalar photon theory, I compute the probability of the atoms remaining excited and the probability for the atoms remaining in the initial quantum state of a system of two-level atoms cloud in a sphere initially prepared to radiate in the forward direction, i.e., the single photon superradiance problem. The convergence in the results obtained for increasingly larger radii for the sphere suggests that the asymptotic limits for these quantities are obtained for a sphere with a radius equal to six times the resonant wavelength. I predict the maximal value of the probability of secondary excited states from large spheres at 17.1%.