密度矩阵重整化群及其在凝聚态物理中的应用

Slides:



Advertisements
Similar presentations
定 格 入 格 破 格 —— 新诗仿写复习训练 仿照下列句子,再把 “ 人生 ” 比喻成 “ 大海 ”“ 天空 ” , 造两个句子。 如果说人生是一首优美的乐曲,那么痛苦则 是其中一个不可或缺的音符。 参考答案: 1 、如果说人生是一望无际的大海,那么挫折则 是其中一个骤然翻起的浪花。 2 、如果说人生是一片湛蓝的天空,那么失意则.
Advertisements

第五节 函数的微分 一、微分的定义 二、微分的几何意义 三、基本初等函数的微分公式与微分运算 法则 四、微分形式不变性 五、微分在近似计算中的应用 六、小结.
2.8 函数的微分 1 微分的定义 2 微分的几何意义 3 微分公式与微分运算法则 4 微分在近似计算中的应用.
第七节 函数的微分 一 、微分 概念 二、微分的几何意义 三、 基本初等函数的微分公 式与 微分运算法则 四 、小结.
2.5 函数的微分 一、问题的提出 二、微分的定义 三、可微的条件 四、微分的几何意义 五、微分的求法 六、小结.
第二章 导数与微分. 二、 微分的几何意义 三、微分在近似计算中的应用 一、 微分的定义 2.3 微 分.
全微分 教学目的:全微分的有关概念和意义 教学重点:全微分的计算和应用 教学难点:全微分应用于近似计算.
第三节 微分 3.1 、微分的概念 3.2 、微分的计算 3.3 、微分的应用. 一、问题的提出 实例 : 正方形金属薄片受热后面积的改变量.
传媒学生应该如何度 过四年大学生活?. 进入大学一个多月了,用一个词形容大 学生活 自卑感 不适应 空虚感 被动感 孤独感 失望感 一、大学新生不适应大学生活的表现:
信号与系统 第三章 傅里叶变换 东北大学 2017/2/27.
王同学的苦恼﹗ MC 4.1 诚可贵﹗.
*第 7 章 量子力学中的矩阵形式 与表象变换.
学党章党规、学系列讲话,做合格党员 学习教育
政府採購法規概要 報告人:杜國正 行政院公共工程委員會企劃處.
第五章 二次型. 第五章 二次型 知识点1---二次型及其矩阵表示 二次型的基本概念 1. 线性变换与合同矩阵 2.
第三章 二次量子化之基礎理論.
之 魔 析 妖 鬼 解 怪 大 沈家仪小组出品.
第1章第3节 量化研究与质化研究 案例1:关于中学思想政治教师专业发展现状和需求的调查研究
团队介绍 (1)西湖区社区街道挂职社会实践基地 (2)武义、缙云、双浦乡镇挂职社会实践基地 (3)BOX企业实习社会实践基地
第九章 长期资产及摊销 2017/3/21.
第5章 定积分及其应用 基本要求 5.1 定积分的概念与性质 5.2 微积分基本公式 5.3 定积分的换元积分法与分部积分法
项目申报及投资推进工作实务 更多模板、视频教程: 兰溪市发展和改革局 2013年9月 1.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
2-7、函数的微分 教学要求 教学要点.
§5 微分及其应用 一、微分的概念 实例:正方形金属薄片受热后面积的改变量..
第三讲 矩阵特征值计算及其应用 — 正交变换与QR方法.
电子教案 量子力学教程(第二版) 湖州师范学院 编 主 编 于少英 沈彩万 参 编 刘艳鑫 董永胜 董国香 邱为钢 李艳霞
組員: 陳曉東(2), 張梓錕(3), 蔡浩維(5), 鍾智灝(7), 馮浩然(8)
§3.7 热力学基本方程及麦克斯韦关系式 热力学状态函数 H, A, G 组合辅助函数 U, H → 能量计算
盒子中之自由電子 Particle in a box
第六章 自旋和角动量 复旦大学 苏汝铿.
Raman Spectra 姚思嘉 合作者:蔺楠、尹伊伦.
Oct 杨 金 民 中国科学院理论物理研究所 怀柔暑期学校( )
动态规划(Dynamic Programming)
排列组合 1. 两个基本原理 分类加法计数原理 分步乘法计数原理.
第三章 矩阵力学基础 ——力学量和算符 复旦大学 苏汝铿.
若2002年我国国民生产总值为 亿元,如果 ,那么经过多少年国民生产总值 每年平均增长 是2002年时的2倍? 解:设经过 年国民生产总值为2002年时的2倍, 根据题意有 , 即.
第4章 非线性规划 4.5 约束最优化方法 2019/4/6 山东大学 软件学院.
从物理角度浅谈 集成电路 中的几个最小尺寸 赖凯 电子科学与技术系 本科2001级.
第7讲 自旋与泡利原理.
光子能量线性_不同灵敏层厚度 photon,Cell Size 5x5mm
1.5 函数y=Asin(ωx+φ)的图象.
线 性 代 数 厦门大学线性代数教学组 2019年4月24日6时8分 / 45.
第三节 常见天气系统.
第9讲 原子光谱项.
WPT MRC. WPT MRC 由题目引出的几个问题 1.做MRC-WPT的多了,与其他文章的区别是什么? 2.Charging Control的手段是什么? 3.Power Reigon是什么东西?
3. 分子动力学 (Molecular Dynamics,MD) 算法
北师大版五年级数学下册 分数乘法(一).
第十章重整化群理论 10.1 引言 相变理论的核心问题:求临界指数和标度律,阐明相变普适类的根源。
第16讲 相似矩阵与方阵的对角化 主要内容: 1.相似矩阵 2. 方阵的对角化.
红利、年金、满期金自动转入聚宝盆,收益有保底,升值空间更大
§5.3 泡利原理和同科电子 一、确定电子状态的量子数 标志电子态的量子数有五个:n,l,s,ml,ms。
第18 讲 配合物:晶体场理论.
第七、八次实验要求.
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
第15讲 特征值与特征向量的性质 主要内容:特征值与特征向量的性质.
线 性 代 数 厦门大学线性代数教学组 2019年5月12日4时19分 / 45.
GIS基本功能 数据存储 与管理 数据采集 数据处理 与编辑 空间查询 空间查询 GIS能做什么? 与分析 叠加分析 缓冲区分析 网络分析
第二节 函数的极限 一、函数极限的定义 二、函数极限的性质 三、小结 思考题.
热力学第一定律的应用 --理想气体等容过程、定容摩尔热容 --理想气体等压过程 、定压摩尔热容.
§2 方阵的特征值与特征向量.
衰变中的末态相互作用研究 梁伟红 宫昊 (广西师范大学) 研究目的 衰变中的末态相互作用 工作的不足与展望.
百艳图.
第三章 从概率分布函数的抽样 (Sampling from Probability Distribution Functions)
本底对汞原子第一激发能测量的影响 钱振宇
第十七讲 密码执行(1).
第十二讲 密码执行(上).
§4.5 最大公因式的矩阵求法( Ⅱ ).
§1 向量的内积、长度及正交性 1. 内积的定义及性质 2. 向量的长度及性质 3. 正交向量组的定义及求解 4. 正交矩阵与正交变换.
海报题目 简介: 介绍此项仿真工作的目标和需要解决的问题。 可以添加合适的图片。
Presentation transcript:

密度矩阵重整化群及其在凝聚态物理中的应用 向 涛 中科院理论物理所

凝聚态物理多体理论需要解决的问题是什么? 一个多体相互作用系统在某个特定状态 (例如基态) 下的物理性质 困难点:不可微扰

优化处理多粒子相互作用体系的一种数值重整化群方法 密度矩阵重整化群 系统的总自由度随粒子数呈指数增长: mN (m = 2, 3, …, N ~1023) 优化处理多粒子相互作用体系的一种数值重整化群方法 用密度矩阵挑选所要保留的基矢 用有限的几个基矢来近似表示一个无穷维空间中的一些状态

S = 1/2 Heisenberg 模型 Total degrees of freedom: 2N 量子效应:

Heisenberg 相互作用: H2 分子 能量 三重态 单态 J

Particle in a box

所研究的矩阵的特点 维数高: mN 稀疏:90%或更多矩阵元为零 有一定的对称性(或守恒量〕:矩阵可分块对角化

重整化群思想 标度变换: 作用量A与A’具有相同的泛函形式(称之为可重整性),这也是量子场论方法的基础 重正化群:只是一个半群

经典重整化群方法:按能量保留状态 保留 H2 的 p 最小本征态 保留 H4 的 p 最小本征态

两个开边界子系统合在一起其衔接部分的状态与实际差的很远 经典重整化群方法失败的原因 边界误差太大 切断误差太大 共 p2个状态 仅 p 个被保留 按能量取舍状态 有可能丢掉了一些 有用的状态而保留 了一些无用的状态 两个开边界子系统合在一起其衔接部分的状态与实际差的很远

改进的重整化群方法 边界误差减小 切断误差减小 2p个状态,保留p个

密度矩阵重整化群 Superblock 按系统的约化密度矩阵的本征值保留状态 系统 环境

约化密度矩阵 约化密度矩阵的本征值  等于其对应的本征态 |> 在基态上的投影振幅

DMRG 迭代过程 系统和环境中各加进一个点并初始化或更新 H=Hsys+Henv+Hsys,env 用Lanczos或其它稀疏矩阵对角化方法对角化H 求出基态波函数 做基矢切断并求出变换矩阵 Unp 构造并对角化约化密度矩阵

Lanczos方法

DMRG与其它方法比较 Monte Carlo或其它近似方法 误差 ~ 1% 总自由度数:2L 1D量子系统DMRG的误差远小于其它近似方法

密度矩阵重整化群方法发展的主要进展 零温,实空间:1992 热力学计算(TMRG):经典系统 1995,1D量子系统 1996 高维空间:动量空间1995,分子第一性原理计算1998 ,待进一步发展 动力学关联函数计算:零温及1D有限温度 1999 非平衡态 (含时演化) 问题:2001,待进一步发展 与Monte Carlo方法的结合:1999,有很大的发展空间

计算量 主要CPU时间用于矩阵的对角化 实际计算的矩阵的维数:104 - 106 稀疏程度:10-30% 需要对角化的矩阵的个数:103 - 105 矩阵与矢量相乘的总次数:105 - 107 硬盘:10G - 200G

转移矩阵重整化群:有限温度DMRG方法 转移矩阵 空间 时间

转移矩阵重整化群与DMRG的比较 T=0 DMRG TMRG Target Matrix Hamiltonian H Symmetric Transfer Matrix T Non-symmetric Target State Ground state max | max> Density matrix Lattice size Finite Infinity (Finite time slices)

S=1/2 Heiserberg 模型的磁化率

S=1/2 Heiserberg 模型的关联长度

二维密度矩阵重整化群方法 核心问题:2D格子如何向1D格子映射? 多链方法 2D方法 Condmat/0102200

Heisenberg 模型的基态性质 Square Lattice Triangle Lattice Square Triangle DMRG -0.3346 -0.1814 MC -0.334719 -0.1819 SW -0.33475 -0.1822

小 结 密度矩阵重整化群是目前研究一维量子多体系统最为精确的数值计算方法 但在研究高维或非平衡态系统的物理性质方面还有许多需要解决的数学问题