CN115698897A - 具有多体相互作用和全对全连接的超级伊辛仿真器 - Google Patents

具有多体相互作用和全对全连接的超级伊辛仿真器 Download PDF

Info

Publication number
CN115698897A
CN115698897A CN202180040984.4A CN202180040984A CN115698897A CN 115698897 A CN115698897 A CN 115698897A CN 202180040984 A CN202180040984 A CN 202180040984A CN 115698897 A CN115698897 A CN 115698897A
Authority
CN
China
Prior art keywords
optical
phase
computing system
generate
pump beam
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202180040984.4A
Other languages
English (en)
Inventor
S·库马尔
H·张
黄瑜平
T·卜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Stevens Polytechnic Board Of Directors
Original Assignee
Stevens Polytechnic Board Of Directors
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Stevens Polytechnic Board Of Directors filed Critical Stevens Polytechnic Board Of Directors
Publication of CN115698897A publication Critical patent/CN115698897A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/35—Non-linear optics
    • G02F1/37—Non-linear optics for second-harmonic generation
    • G—PHYSICS
    • G02—OPTICS
    • G02B—OPTICAL ELEMENTS, SYSTEMS OR APPARATUS
    • G02B27/00—Optical systems or apparatus not provided for by any of the groups G02B1/00 - G02B26/00, G02B30/00
    • G02B27/10—Beam splitting or combining systems
    • G02B27/14—Beam splitting or combining systems operating by reflection only
    • G02B27/141—Beam splitting or combining systems operating by reflection only using dichroic mirrors
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/01—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics for the control of the intensity, phase, polarisation or colour 
    • G02F1/0102—Constructional details, not otherwise provided for in this subclass
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/35—Non-linear optics
    • G02F1/355—Non-linear optics characterised by the materials used
    • G02F1/3551—Crystals
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/35—Non-linear optics
    • G02F1/355—Non-linear optics characterised by the materials used
    • G02F1/3558—Poled materials, e.g. with periodic poling; Fabrication of domain inverted structures, e.g. for quasi-phase-matching [QPM]
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F1/00—Devices or arrangements for the control of the intensity, colour, phase, polarisation or direction of light arriving from an independent light source, e.g. switching, gating or modulating; Non-linear optics
    • G02F1/35—Non-linear optics
    • G02F1/39—Non-linear optics for parametric generation or amplification of light, infrared or ultraviolet waves
    • G02F1/392—Parametric amplification
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06E—OPTICAL COMPUTING DEVICES; COMPUTING DEVICES USING OTHER RADIATIONS WITH SIMILAR PROPERTIES
    • G06E1/00—Devices for processing exclusively digital data
    • G06E1/02—Devices for processing exclusively digital data operating upon the order or content of the data handled
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N5/00—Computing arrangements using knowledge-based models
    • G06N5/01—Dynamic search techniques; Heuristics; Dynamic trees; Branch-and-bound
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F2202/00—Materials and properties
    • G02F2202/20—LiNbO3, LiTaO3
    • G—PHYSICS
    • G02—OPTICS
    • G02F—OPTICAL DEVICES OR ARRANGEMENTS FOR THE CONTROL OF LIGHT BY MODIFICATION OF THE OPTICAL PROPERTIES OF THE MEDIA OF THE ELEMENTS INVOLVED THEREIN; NON-LINEAR OPTICS; FREQUENCY-CHANGING OF LIGHT; OPTICAL LOGIC ELEMENTS; OPTICAL ANALOGUE/DIGITAL CONVERTERS
    • G02F2203/00—Function characteristic
    • G02F2203/12—Function characteristic spatial light modulator

Landscapes

  • Physics & Mathematics (AREA)
  • Nonlinear Science (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Optics & Photonics (AREA)
  • Computing Systems (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Artificial Intelligence (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computational Linguistics (AREA)
  • Probability & Statistics with Applications (AREA)
  • Optical Modulation, Optical Deflection, Nonlinear Optics, Optical Demodulation, Optical Logic Elements (AREA)

Abstract

一种光学计算系统包括：光源，被配置为产生泵浦光束；光学调制器，被配置为基于调制掩模调制泵浦光束，以生成调制光束；非线性介质，被配置为将调制光束的一部分转换为二次谐波(SH)光束并产生包括SH光束和泵浦光束的未转换部分的输出；以及分色镜，被配置为接收非线性介质的输出并将SH光束和泵浦光束的未转换部分去耦合；检测器，被配置为检测泵浦光束的未转换部分的第一光功率，并检测SH光束的第二光功率；以及控制器，被配置为基于第一和第二光功率生成更新调制掩模，以传输到光学调制器。

Description

具有多体相互作用和全对全连接的超级伊辛仿真器
交叉引用相关申请
本申请要求2020年5月13日在美国专利商标局提交的美国临时专利申请号63/024,257的优先权和权益，其全部内容通过引用结合于此。
技术领域
本公开的实施例的各方面通常涉及用于解决组合优化问题的处理系统及其使用方法。
发明背景
光子系统目前正用于寻找复杂系统的近似基态，例如，具有时间复用脉冲的光纤环形光学参量振荡器腔中的快速相干伊辛机，尽管其具有有限的自旋数或依赖于光电检测和电子反馈来模拟自旋间相互作用。此外，基于空间光调制的线性光学伊辛机通过将其编码为空间光调制器(SLM)上的像素的二进制相位，被显示为对着大约80,000个自旋。但是，该机器仅限于两体相互作用。然而，有一些物理系统和数值模型的动力学不能被两体相互作用完全捕获，需要对多体相互作用进行适当的描述，例如，k-SAT(k-可满足性)问题。这提出了一个重大的计算挑战，其复杂性和体积远远超过只有两体相互作用的伊辛问题，即使对于中等数量的自旋也是如此。
背景技术部分中公开的上述信息仅用于增强对本公开的理解，因此其可能包含不构成本领域普通技术人员已知的现有技术的信息。
发明内容
本公开的实施例的各方面涉及用于仿真伊辛机并且能够同时支持高连通性、多体相互作用和大量自旋的系统和方法。
根据本发明的一些实施例，提供了一种光学计算系统，包括：光源，被配置为产生泵浦光束；光学调制器，被配置为接收调制掩模并基于调制掩模调制泵浦光束，以生成调制光束；非线性介质，被配置为将调制光束的一部分转换为二次谐波(SH)光束并产生包括SH光束和泵浦光束的未转换部分的输出；以及分色镜，被配置为接收非线性介质的输出并将SH光束和泵浦光束的未转换部分去耦合；检测器，被配置为检测泵浦光束的未转换部分的第一光功率，并检测SH光束的第二光功率；以及控制器，被配置为基于第一和第二光功率生成更新调制掩模，以传输到光学调制器。
在一些实施例中，光学调制器包括空间光调制器，空间光调制器被配置为调制泵浦光束的小波的相位，以生成调制光束。
在一些实施例中，调制掩模是随机二进制相位图案。
在一些实施例中，空间光调制器被配置为通过基于调制掩模将小波中的每一个的相位编码为第一相位或第二相位，来调制泵浦光束的小波的相位，并且第一和第二相位相差180度。
在一些实施例中，光学调制器包括数字微镜器件，数字微镜器件被配置为调制泵浦光束的小波的强度，以生成调制光束。
在一些实施例中，光学计算系统进一步包括：第一透镜，被配置为将调制光束聚焦在非线性介质内部；以及第二透镜，被配置为将非线性介质的输出准直到分色镜上。
在一些实施例中，检测器包括：第一光电检测器，被配置为检测泵浦光束的未转换部分的第一光功率；以及第二光电检测器，被配置为检测SH光束的第二光功率，并且其中，非线性介质包括周期极化铌酸锂晶体。
在一些实施例中，调制掩模包括多个像素，并且控制器被配置为基于第一和第二光功率来识别像素簇，并且通过反转像素簇的像素来生成更新调制掩模。
根据本发明的一些实施例，提供了一种控制光学计算系统以确定系统能量的近似基态的方法，该方法包括：驱动光源，以生成泵浦光束；生成相位掩模；将相位掩模应用于空间光调制器，以调制泵浦光束的小波的相位，从而基于相位掩模生成调制光束；以及从检测器接收泵浦光束的第一光功率和基于泵浦光束生成的二次谐波(SH)光束的第二光功率，其中，生成相位掩模包括：识别相位掩模的簇；基于第一光功率确定当前系统能量；至少基于系统能量和先前最小系统能量来确定是否更新相位掩模；以及响应于确定更新相位掩模，基于相位掩模的簇更新相位掩模。
在一些实施例中，生成相位掩模进一步包括：响应于确定不更新相位掩模，将当前系统能量和先前最小系统能量中的较小者识别为近似基态。
在一些实施例中，生成相位掩模进一步包括：识别热能；基于当前系统能量和先前最小系统能量计算系统能量的变化；以及基于系统能量的变化和热能确定玻尔兹曼概率，确定是否更新相位掩模进一步基于玻尔兹曼概率。
在一些实施例中，确定是否更新相位掩模包括：确定反馈迭代计数小于阈值；确定更新相位掩模。
根据本发明的一些实施例，提供了一种光学计算系统，包括：光源，被配置为产生泵浦光束；光学调制器，被配置为接收调制掩模并基于调制掩模调制泵浦光束，以生成调制光束；非线性介质，被配置为将调制光束的一部分转换为二次谐波(SH)光束并产生包括SH光束和泵浦光束的未转换部分的输出；以及分色镜，被配置为接收非线性介质的输出并将SH光束和泵浦光束的未转换部分去耦合；第一相机，被配置为捕获泵浦光束的未转换部分并生成第一强度矩阵；第二相机，被配置为捕获泵浦光束的未转换部分并生成第二强度矩阵；以及控制器，被配置为基于第一和第二强度矩阵中的至少一个生成更新调制掩模，以传输到光学调制器。
在一些实施例中，光学调制器包括空间光调制器，空间光调制器被配置为调制泵浦光束的小波的相位，以生成调制光束，并且调制掩模是二进制相位图案。
在一些实施例中，空间光调制器被配置为通过基于调制掩模将小波中的每一个的相位编码为第一相位或第二相位，来调制泵浦光束的小波的相位，并且第一和第二相位相差180°。
在一些实施例中，光学调制器包括数字微镜器件，数字微镜器件被配置为调制泵浦光束的小波的强度，以生成调制光束。
在一些实施例中，控制器被配置为通过以下方式生成更新调制掩模：将相机状态确定为第一强度矩阵、第二强度矩阵或者第一和第二强度矩阵的加权和；以及基于调制掩模和相机状态来确定更新调制掩模。
在一些实施例中，基于调制掩模和相机状态确定更新调制掩模包括：将更新调制掩模确定为调制掩模与相机状态和反馈步长的乘积之间的差。
在一些实施例中，泵浦光束包括由光源生成的高斯脉冲序列中的高斯脉冲，并且反馈步长对应于高斯脉冲序列的周期性。
在一些实施例中，控制器被配置为在生成更新调制掩模之前适应性地调整反馈步长。
根据本发明的一些实施例，提供了一种光学计算系统，包括：光学放大器，被配置为生成具有平坦相位波前的泵浦光束；简并光学参量放大器，被配置为基于泵浦光束在小波中生成双波长信号；滤光器，被配置为从信号中过滤泵浦光束；光学分束器，被配置为读出信号的一部分；非线性光学晶体，被配置为基于信号生成二次谐波光；第一双波段反射器，被配置为分离信号和二次谐波光；光学接收器，被配置为将二次谐波光耦合到光学放大器中；以及第二双波段反射器，被配置为将信号和二次谐波光组合到简并光学参量放大器中。
在一些实施例中，光学计算系统进一步包括：第一傅立叶透镜，被配置为将信号聚焦到非线性光学晶体中，以生成二次谐波光；以及第二傅立叶透镜，其位于第一和第二双波段反射器之间，并被配置为使信号散焦。
在一些实施例中，第一和第二双波段反射器中的每一个都包括分色镜或滤光器。
在一些实施例中，随着信号连续穿过光学计算系统，光学放大器的光学增益增加，直到光学参量振荡发生。
在一些实施例中，光学计算系统还包括位于简并光学参量放大器和非线性光学晶体之间的第一光学调制器，第一光学调制器被配置为在信号到达非线性光学晶体之前调制信号，以改变光学计算系统的有效四体相互作用哈密顿函数。
在一些实施例中，光学计算系统还包括位于非线性光学晶体和光学接收器之间的第二光学调制器，第二光学调制器被配置为调制二次谐波光，以改变光学计算系统的有效四体相互作用哈密顿函数。
根据本发明的一些实施例，提供了一种光学计算系统，包括：光学放大器，被配置为生成具有平坦相位波前的放大信号；非线性光学器件，被配置为接收放大信号并基于放大信号通过二次谐波生成来生成泵浦光束；简并光学参量放大器，被配置为基于泵浦光束在小波中生成双波长信号；滤光器，被配置为从信号中过滤泵浦光束；第一光学分束器，被配置为读出信号的一部分；第二光学分束器，被配置为将信号的第一部分引导至用于调制信号的该部分的光学调制器，并将信号的第二部分引导至分色镜；以及光学接收器，被配置为将信号的调制的第一部分耦合到光学放大器中，以生成放大信号，其中，分色镜被配置为将信号的第二部分与泵浦光束组合到简并光学参量放大器中。
在一些实施例中，光学调制器包括空间光调制器或数字微镜器件，并且非线性光学器件包括非线性光学晶体或光波导。
附图说明
附图与说明书一起示出了本公开的示例实施例，并且与说明书一起用于解释本公开的原理。
图1示出了根据本公开的一些实施例的光学计算系统。
图2示出了根据本公开的一些实施例的控制光学计算系统以识别系统能量的近似基态的过程。
图3-图6是示出根据本公开的一些示例实施例的光学计算系统将其输出演化到其能量基态附近的能力的示图。
图7示出了根据本公开的一些实施例的利用自松弛反馈机制的光学计算系统。
图8示出了根据本公开的一些实施例的由控制器进行的状态更新的示例。
图9A-图9B和图10示出了根据本公开的一些实施例的利用光学反馈的光学计算系统。
图11(a)-图11(c)示出了根据本公开的一些实施例的在伊辛模型中分割的示图。
具体实施方式
以下阐述的详细描述旨在作为根据本公开提供的用于仿真伊辛机以解决不可能或难以解决的组合优化问题的系统和方法的示例实施例的描述，并且不旨在表示可以构建或利用本公开的唯一形式。该描述结合所示的实施例阐述了本公开的特征。然而，应当理解，相同或等效的功能和结构可以通过不同的实施例来实现，这些实施例也旨在包含在本公开的范围内。如本文别处所示，相似的元件编号旨在指示相似的元件或特征。
在整个说明书中，术语在上下文中可能具有暗示或暗指的微妙含义，而不仅仅是明确陈述的含义。
本发明涉及一种能够解决许多组合优化问题的光学计算系统/伊辛机，这些问题使用现有的现代计算机是不可能或难以解决的。更具体地说，本发明涉及可以解决组合优化问题的机器，这些问题涉及诸如蛋白质折叠和药物发现的生命科学、研究社交网络影响的社会学、寻找具有最大边数的切割(MAX-CUT)的图论以及确定布尔可满足性(k-SAT)的计算机科学等领域。更一般地，本发明涉及非确定性的、多项式的、时间困难的问题，例如，以上那些问题，这些问题的根源在于多体相互作用。
根据一些实施例，本公开的系统和方法涉及一种光学计算系统(在本文中也称为伊辛机/仿真器)，其能够进行可调节的两体相互作用、四体相互作用，在大量自旋上具有全对全连接，以托管伊辛机。使用空间光调制器(SLM)可以灵活地控制每个自旋的能量(例如，化学势)及其对总能量的共同贡献。泵浦功率及其二次谐波(SH)功率的光学测量可以分别模拟两体和四体相互作用。基于蒙特卡罗的自适应反馈控制器用于演化有效自旋系统，以找到给定伊辛问题的近似基态解。由本公开的系统模拟的高阶、多体相互作用也可以用作光学机器学习的强大激活功能。诸如矩阵乘法、傅立叶变换和四体相互作用等主要计算任务经由非线性光学系统来执行或由非线性光学系统来模拟。目前的伊辛模拟器可能会为大数据分析和量子模拟的其他不可及的领域开辟一条道路。
根据一些实施例，光学计算系统包括非线性光学系统，该非线性光学系统能够仿真包含非常大量的自旋(例如，一百万个自旋或更多)的伊辛机。光学计算系统能够操纵具有全对全连接的两体和四体相互作用。利用空间光调制器以相干激光束中的小波的二进制相位值的形式编码和控制自旋，并在傅立叶平面上模拟非线性晶体中频率转换的高阶相互作用。通过自适应反馈控制，系统可以演化成有效的自旋配置，这些配置很好地逼近具有全对全连接的多体相互作用的伊辛哈密顿函数的基态。
实际上，光学计算系统使用激光器和空间光调制器对给定的多体问题进行编码，并通过光学元件和由其执行的非线性处理来执行计算操作，例如，傅立叶变换、矩阵乘法和非线性相互作用。光学计算系统使用与光功率计或相机连接的单模光纤来执行测量，并使用检测到的信号来提供自适应反馈(例如，经由蒙特卡罗自旋翻转)，以演进系统。光学计算系统包括激光器、空间光调制器、光学元件、非线性晶体、单模光纤以及一个或多个相机和/或功率计。
图1示出了根据本公开的一些实施例的光学计算系统100。
光学计算系统100(也称为基于光学的计算系统)可以是能够在大量自旋上模拟化学势、两体相互作用和四体相互作用的非线性光学伊辛机。因此，光学计算系统是在执行计算时操纵光子束的系统。
根据一些实施例，光学计算系统(例如，伊辛仿真器)100包括光源(例如，激光器)102、空间光调制器(SLM)104、非线性介质(例如，周期极化铌酸锂(PPLN)晶体)106、分色镜108、检测器110和控制器112。
光源102被配置为产生输入到SLM 104的高斯泵浦光束。在一些示例中，光源102生成光脉冲序列。入射到SLM 104上的泵浦光束的横向FWHM可以是大约2.6mm或者不同的值，这取决于SLM104的尺寸。
SLM 104被配置为接收调制掩模(例如，相位掩模)并调制泵浦光束的小波的相位，以基于相位掩模生成调制的(例如，空间调制的)小波。在一些示例中，SLM 104可以具有足够高的分辨率(例如，1440×1050像素)，以将超过一百万像素/自旋的相位掩模编码到高斯光束上。在一些实施例中，相位掩模是二进制相位掩模，每个像素对应于0或π的相位值。在一些实施例中，光学计算系统100还包括第一透镜105(例如，具有大约200mm焦距的聚焦透镜)，被配置为将调制光束聚焦在非线性介质106内部(例如，靠近中心)。
在一些实施例中，可以是温度稳定的PPLN晶体的非线性介质106被配置为执行非线性过程，以基于调制的泵浦光束生成二次谐波(SH)光束。SH生成实现了自旋之间的有效相互作用。在一些示例中，对于SH生成，PPLN晶体可以具有大约19.36μm的极化周期和大约1cm的长度。晶体内的泵浦束腰大约为45μm。
光学计算系统利用第二透镜(例如，准直透镜)107将由非线性介质106输出的SH光束和剩余的泵浦光束(即，未转换的泵浦光束)准直(例如，对准)到分色镜108上。在一些实施例中，分色镜108(或双波长反射器)被配置为过滤非线性介质106的准直输出，以去耦合(例如，空间分离)SH光束和剩余的/未转换的泵浦光束，该光束可以使用光纤耦合器/准直器(例如，非球面透镜)111a和111b耦合到单模光纤(SMF)109a和109b中。在一些实施例中，检测器110包括分别耦合到SMF 109a和109b的两个光功率计(例如，光电二极管)110a和110b，这两个光功率计分别测量SH光束和剩余泵浦光束的光功率。
在一些实施例中，光学计算系统使用测量的泵浦功率来模拟与两体(即自旋间)相互作用相关联的能量，并且测量的SH功率用于捕获自旋之间的四体相互作用。总势能(例如，总化学势能)是在SLM104处编码的所有自旋的加权和。控制器112然后基于测量的泵浦和SH功率更新SLM 104处的相位掩模，以迭代地找到由光学计算系统100模拟的伊辛问题的最佳基态解或其近似值。
根据一些实施例，检测器110使用具有许多像素的两个CCD相机来测量泵浦和SH光束的空间光功率。控制器112然后使用CCD相机结果作为SLM 104上的下一个迭代步骤的输入。
在其他实施例中，使用光学环路来实现反馈控制，并且空间光调制器由非线性光学器件代替。
光学计算系统的有效哈密顿函数表示为：
H＝αH1+βH2+γH4， 等式(1)
其中，α、β和γ是哈密顿函数的可调参数，H1、H2和H4分别表示化学势、两体相互作用和四体相互作用。这些可以表示为
Figure BDA0003985256310000101
Figure BDA0003985256310000102
以及 等式(3)
Figure BDA0003985256310000103
在此处，光源102产生具有峰值振幅E0和束腰wp的波长为λp的高斯泵浦光束。泵浦光束入射到SLM 104上，SLM 104的相位掩模由以(x′m，y′m)为中心的m×n个像素(m，n)组成。紧接在SLM 104之后的横向电场是
Figure BDA0003985256310000104
其中，
Figure BDA0003985256310000111
是像素(m，n)处的幅度，Π是宽度a的矩形函数，σmn＝±1(二进制相位调制)。
焦距为F的傅立叶透镜可以转换长度为L的周期极化铌酸锂(PPLN)晶体中心的电场。给出
Figure BDA0003985256310000112
其中，
Figure BDA0003985256310000113
κp＝2πnp/λp，np是PPLN晶体中泵浦的折射率。由于只有近轴光被耦合到SMF 109a和109b中，这使得sinc(axπ/λpf)＝sinc(ayπ/λpf)≈1，并且将上述等式简化为
Figure BDA0003985256310000114
通过使用简化的符号
Figure BDA0003985256310000115
和
Figure BDA0003985256310000116
来简化上述等式中的二重求和，其中，i＝1，2，...N，具有N＝N1×N1个自旋，等式7可以表示为
Figure BDA0003985256310000117
非线性介质106内部的泵浦光束和SH光束的动态特性可以通过以下等式来评估
Figure BDA0003985256310000121
Figure BDA0003985256310000122
其中，κSH＝2πnSH/λSH＝nSHωSH/c是SH光在折射率为nSH的晶体中的波数，ωSH是SH光的频率，Δκ＝2κp-κSH-2π/Λ是相位失配，极化周期为Λ。假设具有未耗尽泵浦和可忽略不计的小衍射的相位匹配条件，
Figure BDA0003985256310000123
其中，
Figure BDA0003985256310000124
在非线性介质106的输出端，泵浦波和SH波均耦合到单模光纤中，并被光功率计(例如，光电二极管)110a和110b检测到。泵浦波和SH波的光功率分别是
Figure BDA0003985256310000125
和
Figure BDA0003985256310000126
其中，
Figure BDA0003985256310000127
是束腰分别为
Figure BDA0003985256310000128
和
Figure BDA0003985256310000129
的归一化反向传播光纤模式。
因此，检测到的泵浦波和SH波的光功率是
Figure BDA00039852563100001210
和
Figure BDA0003985256310000131
其中，
Figure BDA0003985256310000132
以及
Figure BDA0003985256310000133
因此：
Figure BDA0003985256310000134
和
Figure BDA0003985256310000135
其中，
Figure BDA0003985256310000136
和
Figure BDA0003985256310000137
分别为两体和四体相互作用项。
在此处，系统的总能量(包括所有自旋及其两体和四体相互作用)可以用定义为下式的单个参数来表征
E＝αC+βPp+γPh， 等式(18)
其中，α、β和γ是自由参数，分别定义化学势、两体和四体相互作用能的贡献。通过适当地设置自由参数，光学计算系统100可以建模伊辛问题。
Figure BDA0003985256310000138
是表示化学能的自旋的加权和，局部化学势为μi∈[-1，1]。等式(18)等效于等式(1)中定义的有效哈密顿函数，以求系统的基态解。系统的磁化可以定义为
Figure BDA0003985256310000139
为了求出系统总能量的最小值或最小值的近似值，控制器112通过自适应反馈迭代地更新(例如，优化)SLM 104处的相位掩模。
图2示出了根据本公开的一些实施例的控制光学计算系统100以识别系统能量的近似基态的过程200。应当理解，该过程的步骤顺序不是固定的，并且可以修改、改变顺序、不同地执行、依次执行、并发地或同时地执行，或者改变成任何期望的顺序，如本领域技术人员所认识到的。
在框202，控制器112驱动光源102生成主泵浦光束，生成可以是随机二进制相位掩模的初始相位掩模，并将初始相位掩模应用于SLM 104，以调制泵浦光束的小波的相位，从而基于初始相位掩模生成调制光束。
在框204，控制器112经由等式18定义由光学计算系统100模拟的系统的能量。在这样，控制器112确定自由参数α、β和γ，以分别定义化学势、两体和四体相互作用能的贡献，以便建模由光学计算系统100模拟的伊辛问题。
在框206，控制器112在具有系统的尺寸s和热能τ的相位掩模中定义(例如，随机选择)像素簇。在一些实施例中，相位掩模可以具有带有s个像素的正方形或矩形形状。然而，本公开的实施例不限于此，簇可以具有任何合适的形状。在此处，热能τ可以简单地是该优化算法的参数，而不是物理量。
在框208，控制器112从检测器110接收主泵浦光束的第一光功率和SH光束的第二光功率。
在框210，控制器112使用等式18并基于第一和第二光焦度以及在框204确定的自由参数α、β和γ来确定当前系统能量E新。控制器112将系统能量的变化U计算为当前系统能量和先前最小系统能量E旧之间的差E旧(U＝E新-E旧))。控制器112进一步确定玻尔兹曼概率P＝exp(-U/τ)。
在框212和214，控制器112确定是否在SLM更新相位掩模。这样，在框212，控制器检查是否满足基于当前系统能量和玻尔兹曼概率的优化标准。当系统能量U的变化为负时(即，当前系统能量E新是迄今为止计算的最低系统能量)，或者如果玻尔兹曼概率P大于随机变量F，其为0和1之间的真实值，则满足该标准。当这些条件都不满足时，在框216，控制器112确定更新新相位掩模。这样，控制器112通过翻转先前相位掩模的随机选择的大小为s的簇中的自旋来生成新的相位掩模(框206)。当满足优化标准时，控制器还确定迭代次数(例如，反馈迭代计数)是否小于阈值(例如，1200)。如果是，则控制器112继续在框216更新相位掩模。
当控制器112确定不再更新相位掩模时，在框218，控制器112收集结果，以将当前系统能量和先前最小系统能量中的较小者识别为近似基态。
因此，为了减少(例如，最小化)总能量，控制器112遵循蒙特卡罗方法，适应性地翻转随机选择的大小为s的簇内的所有自旋。在每次迭代中，根据能量变化函数U和玻尔兹曼概率P＝exp(-U/τ)接受或拒绝翻转的自旋，其中，U＝E新-E旧是由于自旋翻转引起的能量变化，τ是热能。为了避免陷入局部最小值，在迭代过程中簇的大小s和热能都是变化的。
图3-图6是示出根据本公开的一些示例实施例的光学计算系统100将其输出演化到其能量基态附近的能力的示图。
图3示出了纯四体相互作用和800×800自旋的全对全连接的结果，即α＝0和β＝0(见等式18)。图3(a)绘出了γ＝-1的SH光束功率演变，其中，系统的基态对应于最大SH功率。在这种情况下，系统表现出类似铁磁性的行为。相反，图3(b)绘出了γ＝1时SH功率的演变，以通过最小化SH功率来求出基态。最小可检测SH功率接近光学传感器的极限(～5nW)。在这个示例中，系统演变成一种类似顺磁性的状态。在图3(a)和图3(b)中，即使图中显示的数据以不同的初始相位掩模(对应于不同的随机自旋配置)开始，系统也演变，以求出基态解。这表明根据本发明的实施例制造的机器可以可靠且有效地进化到其基态附近。
图4示出了800×800自旋的两体相互作用和全对全连接，即α＝0、β＝-1和γ＝0。图4(a)示出了泵浦光束功率的演变。随着泵浦功率的增加，使系统演化，以使总能量E最小化。图4(b)示出了磁化M的演化。当α＝0时，存在自发的对称性破缺，因为如果所有自旋翻转，系统能量保持不变。这样，反馈控制将以相等的概率优化朝向正或负磁状态的自旋。这在图4(b)中很清楚，其中，磁化趋势是双向的。
图5示出了对于800×800自旋的优化，α＝-1、β＝-1和γ＝-1的测量结果。以幂对迭代次数的对数标度示出结果。图5(a)和图5(b)绘出了四种不同初始相位掩模的SH和泵浦功率的演变。图5(c)描绘了系统的相应磁化。图5(d)示出了初始和最终相位掩模，其中，黑色和白色像素分别表示正和负取向自旋。图5(d)(1)-(4)表示对具有不同输入二进制自旋配置的SH和泵浦功率的评估。在所有情况下，这些结果表明，泵浦的光功率及其SH收敛到最佳值，该最佳值属于给定伊辛问题的基态解。
图6示出了α＝1、β＝-0.5和γ＝-1时的结果。图6(a)和图6(b)分别示出了泵浦和SH功率的演变，图6(c)绘出了系统的相应磁化。当泵浦和SH功率都随着迭代次数增加时，磁化强度在负方向上获取最佳值。图6(d)示出了初始和最终相位掩模图案，其中，相位掩模图案的黑色和白色区域分别是正和负取向自旋。图6(d)(1)-(3)表示对具有不同输入二进制相位掩模的SH和泵浦功率的评估。这些结果表明光学伊辛机可以求出给定伊辛问题的近似基态解。
图7示出了根据本公开的一些实施例的利用自松弛反馈机制的光学计算系统700。除了调制器704、相机710a和710b以及控制器712之外，图7的系统700基本上类似于图1的系统100。这样，为了便于描述，系统100和700共有的部件在此不再描述。
参考图7，根据一些实施例，高斯泵浦激光束(例如，波长为1551.5nm)入射到调制器704上，该调制器可以是用于调制入射泵浦光束的相位的SLM 104，或者可以是用于调制入射泵浦光束的强度的数字微镜器件(DMD)。分色镜108被安装成将SH光(例如，波长为775.75nm)从未转换的/剩余的泵浦光束(可以是ab红外光束)中分离出来。在一些实施例中，第一相机710a(例如，红外相机)捕获来自分色镜108的未转换/剩余的泵浦光束的强度，第二相机710b(例如，可见光相机)捕获来自分色镜108的SH光束的强度。由第二相机710b捕获的SH光可以模拟各种四体相互作用哈密顿函数。由第二相机710b处的像素矩阵捕获的强度矩阵/图案可以被写为：
Figure BDA0003985256310000171
其中，J′ijsr是四体相互作用强度。同时，由第一相机710a捕获的未转换的光可以模拟两体相互作用，其中，第一相机处的强度矩阵/图案是：
Figure BDA0003985256310000172
这两个相机可以单独工作，也可以计算其加权和αI2(x，y)+βI4(x，y)，来实现自放松反馈。相机710a和710b的每个像素的强度读数是量化的灰度值。控制器712处理来自相机710a和710b的图像(或强度矩阵)，并控制调制器(例如，SLM或DMD)704，以用新的相位或强度信息更新调制器704处的相位或强度掩模。
根据一些实施例，光学计算系统700将四体相互作用结合到光学环路中，以形成自松弛伊辛机，该自松弛伊辛机可以自动演化成期望的自旋配置。这样，光学计算系统700利用量子系统的虚数时间演化来实现自松弛反馈。为了说明这一点，假设在时间t＝0时调制器704处的相位/强度信息/状态是|ψ(0)>，并且整个系统700的哈密顿函数是H(包括相位/强度调制、通过光学器件的光束传播、二次谐波生成、相机上的强度测量以及像素值的任何数字处理)，则在时间t时的状态是
|ψ(t)>＝e-iHt|ψ(0)>。 等式(21)
从数学的角度来看，实值时间t可以用虚构值-it来代替。然后，该状态将是|ψ(-it)>＝e-iH(-it)|ψ(0)>＝e-Ht|ψ(0)>。初始状态|ψ(0)>是随机态，该状态是不同本征态的线性组合，|ψ(0)>＝∑ncn|n>。因此，时间t的状态可以表示为
Figure BDA0003985256310000181
众所周知，较低的本征态具有较低的本征值，即E0＜E1＜E2＜…因此，随着时间的推移，基态|0>将比其他状态更突出，因为其衰变常数E0导致最慢的衰变速度。在一些实施例中，量子状态在进化的每个步骤之后被重新归一化。
图8示出了根据本公开的一些实施例的由控制器712进行的状态更新的示例。
在图8的示例中，在时间t的调制器状态是|ψ(t)>，反馈环路之间的时间步长是dt。因此，相机状态I4(x，y，t)可以写成H|ψ(t)>。因为在时间t+dt的新状态是t+dt是|ψ(t+dt)>＝e-Hdt|ψ(t)>，并且时间步长dt足够小，以至于可以忽略指数展开的高阶：
Figure BDA0003985256310000182
然后得出：
Figure BDA0003985256310000183
因此，新的调制器状态(|ψ(t+dt)>)可以表示为先前调制器状态|ψ(t)和dt乘以相机状态(H|ψ(t)>)之间的差。
在启动时，控制器712选择随机状态|ψ0>，作为初始调制器相位状态。在穿过非线性晶体之后，上转换光(即SH光束)和未转换光(即剩余的泵浦光束)被相机710a和710b检测到。根据光学计算系统700所模拟的问题，控制器712计算总的相机状态|φ0>，作为来自第一和第二相机710a和710b之一的强度矩阵(即，I2(x，y)or I4(x，y)))或者强度矩阵的加权和(即，αI2(x，y)+βI4(x，y))。
根据一些实施例，控制器712通过计算新状态为|ψ1>＝|ψ0>-Δt|φ0>，并通过初始状态
Figure BDA0003985256310000191
将其归一化，来确定更新调制器相位状态。因此，控制器712确定下一循环迭代的新调制器状态相位为：
Figure BDA0003985256310000192
在此处，反馈环路的时间步长Δt足够小，可以减少或最小化数值误差。在一些示例中，反馈步长Δt表示光源102的脉冲序列的周期性。在一些实施例中，控制器712自适应地调整Δt，以实现自松弛，因为光学计算系统700的动态可以加速或变慢。可以选择步长，使得等式(22)中进行的近似是有效的。
根据一些实施例，光学计算系统可以自动对准其自旋，以最大化或最小化其相互作用能量。
图9A示出了根据本公开的一些实施例的利用光学反馈的光学计算系统900。
在图9的示例中，使用光学环路来提供反馈，其中，只有具有最高四体相互作用能量的自旋配置将满足光学振荡条件并积累。光学环路从光学放大器(OA)902开始，以生成具有平坦相位波前的二次谐波波长的弱泵浦光束，从而驱动简并光学参量放大器(DOPA)904。DOPA 904可以是二阶非线性的相位匹配或准相位匹配光学介质，例如，β硼酸钡(BBO)晶体、磷酸钛氧钾(KTP)晶体或周期极化铌酸锂(PPLN)晶体。在DOPA 904中，泵浦光束生成并放大两倍波长的信号波。泵浦光束是宽的，使得在DOPA 904内部几乎没有衍射，并且在泵浦光束的横截面上独立地生成/放大许多信号小波，每个信号小波表示一个自旋。
在DOPA 904之后，通过使用滤波器906过滤/去除泵浦光束，滤波器906可以是带通滤波器、边缘滤波器、分色镜等，并且使用第一傅立叶透镜908将信号聚焦到二次谐波生成(SHG)装置910中。SHG装置910可以是二阶非线性的相位匹配或准相位匹配光学介质，例如，BBO晶体、KTP晶体或PPLN晶体。在SH生成期间，实现自旋之间的有效相互作用，以产生二次谐波光。在输出端，剩余的泵浦信号和生成的SH光大致在第一双波段反射器(例如，第一分色镜(DM)或滤光器)912或任何其他光学器件处分离。剩余的泵浦信号然后通过第二傅立叶透镜914，以取消聚焦。然后经由第二双波段反射器(例如，第二分色镜或滤光器)916将其导向DOPA 904。另一方面，SH光由光学接收器(OR)918收集，该光学接收器可以是光纤或自由空间光学器件的准直光束。然后，SH光被OA 902放大，并与剩余的泵浦信号组合，并被引导到DOPA 904，作为新的泵浦。该环路将运行多次，直到产生稳定的结果，在此基础上，将经由分束器(BS)920读取信号，以测量其在横向剖面上的自旋。在此处，随着信号连续穿过光学计算系统，光学放大器的光学增益增加，直到发生光学参量振荡。
在这种光学配置中，DOPA 904通过简并光学参量放大生成并放大泵浦光束横截面上的信号小波。因为泵浦光束具有平坦的相位波前，所以在每个空间网格(即像素)处生成的信号小波可以与泵浦光束具有相同的相位，或者相反。这种机制定义了信号小波的离散相位，使得信号小波中的每一个都可以表示自旋向上(相对于泵浦的0相位)或自旋向下(相对于泵浦的pi相位)。最初，没有信号，DOPA 904可以被理解为放大真空噪音的量子过程。真空噪音包含了所有可能的自旋配置。将过滤放大的真空噪音，然后，排空泵浦。然后，噪音将聚焦到二阶非线性晶体中，用于生成SH。在SH生成期间可以产生的二次谐波光的量取决于自旋配置。事实上，其强度可以被写成自旋的四体相互作用的能量。所得到的二次谐波信号然后被OA 902放大，以产生DOPA 904的泵浦，从而在下一轮放大信号。环路将持续运行。
以这种方式，光学计算系统900有效地实现了光学参量振荡器(OPO)，区别在于，在此处，通过SH生成和光学放大的中间步骤来实现光学增益。当操作时，OA 902的光学增益可以缓慢增加，直到出现第一光学参量振荡信号，即，直到可以读出强且稳定的信号。当发生这种情况时，信号将以理想的自旋配置准备，以使四体相互作用最大化。这是因为其他配置将产生较少的二次谐波光，因此有效OPO增益较小。由于OPO要求每个光学环路期间的光学增益大于净损耗，通过缓慢增加OA增益，当期望的自旋配置的光学增益刚好超过净损耗时，只出现第一OPO信号。所有其它的自旋配置将导致更多的损耗而不是增益，因此对于这些其它的自旋配置不可能有光学振荡。在此处，在OPO发生之前，在所有可能的自旋配置中，DOPA过程可以被理解为真空噪音的放大。
为了进一步确保只存在一种自旋配置，可以将OA设置为在光功率饱和极限下操作，即，如果输出功率达到极限，OA将停止放大。以这种方式，在为期望的自旋配置建立了OPO之后，信号和泵浦功率都增加，以使OA饱和，从而降低了其他自旋配置的有效增益。
在一些实施例中，可以通过首先用具有平坦相位波前的本地振荡器对其进行干涉，然后使用CCD相机对其进行测量，来读出信号的自旋配置。在一些实施例中，信号可以通过线性光学元件，例如，傅立叶透镜，从而可以在不同模式的基础上测量其自旋配置。
在一些实施例中，光学环路相对于任何热、机械或其他漂移是稳定的。可以通过反馈控制机制来实现稳定。
根据一些实施例，修改有效相互作用哈密顿函数，以模制感兴趣的实际问题。这可以通过在SHG之前对信号应用光学调制，和/或在被OR收集之前对生成的二次谐波光进行调制来实现。在图9B中描述这个选项。在光学计算系统900-1中，用于泵浦(MP)的光学调制器922和用于信号(MS)的光学调制器924的调制可以具有相位调制、振幅调制或这两者的形式。可以使用空间光调制器、数字微镜器件或这两者的组合来实现调制。对于泵浦光，可以通过其他方式有效地实现调制，例如，将其耦合到单模光纤中。
图9A和图9B具体描述了如何通过优化自旋配置使四体相互作用最大化。通过使DOPA的泵浦功率与二次谐波光的检测功率成反比，可以开发类似的方案，来最小化四体相互作用。
除了最大化或最小化四体相互作用之外，同样的原理可以用来寻找最大化或最小化两体相互作用的自旋配置。
图10给出了这样一个示例，该示例遵循图9A的基本结构，但是部分分离信号，并且不通过SHG 910而直接检测该信号。为了调制有效的两体相互作用哈密顿函数，在检测之前由光学调制器(MS)924调制信号。在一些实施例中，这种调制不是必需的，并且可以被规避。检测到的信号产生具有平坦相位波前的泵浦来驱动DOPA 904。在一些实施例中，这是通过首先放大信号功率，然后将其倍频来实现的，如图9B所示。在一些实施例中，使用可调增益激光源或固定增益激光源(但是具有可调分束器或可变功率衰减器)对信号进行光电检测，并且所得的电子信号用于控制泵浦功率，以驱动DOPA。
如图10所示，根据一些实施例，光学计算系统1000包括光学放大器(OA)902，被配置为生成具有平坦相位波前的放大信号；非线性光学器件(SHG)910，被配置为接收放大信号并基于放大信号通过二次谐波生成来生成泵浦光束；简并光学参量放大器(DOPA)904，被配置为基于泵浦光束在小波中生成双波长信号；滤光器(FT)906，被配置为从信号中过滤泵浦光束；第一光学分束器(BS)920，被配置为读出信号的一部分；第二光学分束器921，被配置为将信号的第一部分引导至光学调制器(MS)924，以调制信号的该部分，并将信号的第二部分引导至分色镜(DM)916；以及光学接收器(OR)918，被配置为将信号的调制的第一部分耦合到光学放大器902中，以生成放大信号。在此处，分色镜916被配置为将信号的第二部分与泵浦光束组合到简并光学参量放大器中。根据一些示例，光学调制器924可以是空间光调制器或数字微镜器件，非线性光学器件910可以是非线性光学晶体或光波导。
在一些实施例中，可以组合两体和四体相互作用的方案，以找到最大化或最小化一个相互作用同时也最大化或最小化另一个相互作用或者最大化或最小化这两者的加权和的自旋配置。
虽然本文公开的光学计算系统仅支持两体和四体相互作用，但是可以采用非线性光学器件的其他过程，例如，和频、四波混频和高阶谐波生成，来实现更高阶的相互作用，其计算复杂度和强度可能会快速增长到超出现有计算平台的能力。
将NP-难度问题映射到伊辛模型的示例：
许多NP-难度问题可以被绘制成图，并且这些图可以被映射成伊辛模型。下面是将划分图问题映射成伊辛问题的示例。
考虑一个简单的无向图G＝(V，E)，其具有N个V顶点的节点和E条边，问题可能是通过划分其节点集来将图分成“m”个子图。有m(N-1)-1种方法将图分成具有N个节点的子图。为了最佳地计算顶点的划分，V＝V0UV1UV2UVm-1，需要满足以下条件：
1){Vi}是不相交的→Vi∩Vj＝Φ
2)所有{Vi}大致平衡→|Vi|～|Vj|
3)切割边的数量应为最小值/最大值，即，|E切割|最小/最大＝{(u，v)|uεVi，vεVj，i≠j}。
求出大量节点的图划分的精确解是NP组合优化问题。为了简单起见，可以考虑一个图，该图具有5个与6条边连接的节点，权重为wij，其中，i≠j是节点，wij＝wji＞0，如图11(a)所示。有许多方法可以将这个图切割成两个子图，以找到通过最大数量的边的切割。图11(b)和11(c)示出了其中的两个。当分割图时，值Si＝{+，-}可以被分配给分离的子图的节点。切割的成本函数可以表示为
Figure BDA0003985256310000241
Figure BDA0003985256310000242
其中，如果在同一个子图中的i，j是1，则
Figure BDA0003985256310000243
的值为零，否则为1。成本函数可以重写为
C＝∑i＜jwij-∑i＜jSiSjwij。 等式(24)
可以通过最小化成本函数的第二项来找到最大切割解。这表明成本函数的最小化等价于伊辛哈密顿函数的能量的最小化。
具有成对和高阶相互作用的伊辛模型是统计物理和机器学习的活跃研究领域。根据本公开的一些实施例，具有大量节点和全对全连接的非线性光学系统是更快解决这些问题的可行解决方案。
由本公开的光学计算系统的组成部件(例如，控制器)执行的操作可以由“处理电路”或“处理器”来执行，该“处理电路”或“处理器”可以包括用于处理数据或数字信号的硬件、固件和软件的任意组合。处理电路硬件可以包括例如专用集成电路(ASIC)、通用或专用中央处理单元(CPU)、数字信号处理器(DSP)、图形处理单元(GPU)以及诸如现场可编程门阵列(FPGA)等可编程逻辑器件。如本文所用，在处理电路中，每个功能由被配置(即硬连线)为执行该功能的硬件来执行，或者由被配置为执行存储在非暂时性存储介质中的指令的更通用的硬件(例如，CPU)来执行。可以在单个印刷线路板(PWB)上制造处理电路，或者处理电路可以分布在几个互连的PWB上。处理电路可以包含其他处理电路；例如，处理电路可以包括在PWB上互连的两个处理电路，即FPGA和CPU。
应当理解，尽管术语“第一”、“第二”、“第三”等可以在此处用来描述各种元件、部件、区域、层和/或部分，但是这些元件、部件、区域、层和/或部分不应该受到这些术语的限制。这些术语用于将一个元件、部件、区域、层或部分与另一个元件、部件、区域、层或部分区分开来。因此，下面讨论的第一元件、部件、区域、层或部分可以被称为第二元件、部件、区域、层或部分，而不脱离本发明概念的范围。
在此处使用的术语是为了描述特定的实施例，而不是为了限制本发明的概念。如本文所用，单数形式“一种”(a)和“一种”(an)也旨在包括复数形式，除非上下文另有明确指示。还应当理解，当在本说明书中使用时，术语“包括”(include)、“包括”(including)、“包含”(comprises)和/或“包含”(comprising)指定存在所述特征、整数、步骤、操作、元件和/或部件，但是不排除存在或添加一个或多个其他特征、整数、步骤、操作、元件、部件和/或其组。如本文所用，术语“和/或”包括一个或多个相关列出项目的任何和所有组合。此外，当描述本发明构思的实施例时，使用“可以”是指“本发明构思的一个或多个实施例”。此外，术语“示例性”旨在表示示例或说明。
如本文所用，术语“使用”(use)、“使用”(using)和“使用的”(used)可以被认为分别与术语“利用”(utilize)、“利用”(utilizing)和“利用的”(utilized)同义。
此外，当描述本发明构思的实施例时，使用“可以”是指“本发明构思的一个或多个实施例”。此外，术语“示例性”旨在表示示例或说明。
虽然已经结合某些示例性实施例描述了本发明，但是应当理解，本发明不限于所公开的实施例，而是相反，旨在覆盖包括在所附权利要求及其等同物的精神和范围内的各种修改和等同布置。

Claims (28)

1.一种光学计算系统，包括：
光源，其配置为产生泵浦光束；
光学调制器，其配置为接收调制掩模并基于所述调制掩模调制所述泵浦光束，以生成调制光束；
非线性介质，其配置为将所述调制光束的一部分转换为二次谐波(SH)光束并产生包括所述SH光束和所述泵浦光束的未转换部分的输出；以及
分色镜，其配置为接收所述非线性介质的所述输出并将所述SH光束和所述泵浦光束的所述未转换部分去耦合；
检测器，其配置为检测所述泵浦光束的所述未转换部分的第一光功率，并检测所述SH光束的第二光功率；和
控制器，其配置为基于所述第一光功率和所述第二光功率生成更新调制掩模，以传输到所述光学调制器。
2.根据权利要求1所述的光学计算系统，其中，所述光学调制器包括空间光调制器，所述空间光调制器配置为调制所述泵浦光束的小波的相位，以生成所述调制光束。
3.根据权利要求2所述的光学计算系统，其中，所述调制掩模是随机二进制相位图案。
4.根据权利要求2所述的光学计算系统，其中，所述空间光调制器配置为通过基于所述调制掩模将小波中的每个小波的相位编码为第一相位或第二相位，来调制所述泵浦光束的所述小波的相位，并且
其中，所述第一相位和所述第二相位相差180度。
5.根据权利要求1所述的光学计算系统，其中，所述光学调制器包括数字微镜器件，所述数字微镜器件配置为调制所述泵浦光束的小波的强度，以生成所述调制光束。
6.根据权利要求1所述的光学计算系统，进一步包括：
第一透镜，其配置为将所述调制光束聚焦在所述非线性介质内部；以及
第二透镜，其配置为将所述非线性介质的所述输出准直到所述分色镜上。
7.根据权利要求1所述的光学计算系统，其中，所述检测器包括：
第一光电检测器，其配置为检测所述泵浦光束的所述未转换部分的所述第一光功率；以及
第二光电检测器，其配置为检测所述SH光束的所述第二光功率，并且
其中，所述非线性介质包括周期极化铌酸锂晶体。
8.根据权利要求1所述的光学计算系统，其中，所述调制掩模包括多个像素，并且
其中，所述控制器配置为基于所述第一光功率和所述第二光功率来识别像素簇，并且通过反转所述像素簇的像素来生成所述更新调制掩模。
9.一种控制光学计算系统以确定系统能量的近似基态的方法，所述方法包括：
驱动光源，以生成泵浦光束；
生成相位掩模；
将所述相位掩模应用于空间光调制器，以调制所述泵浦光束的小波的相位，从而基于所述相位掩模生成调制光束；以及
从检测器接收所述泵浦光束的第一光功率和基于所述泵浦光束生成的二次谐波(SH)光束的第二光功率，
其中，生成所述相位掩模包括：
识别所述相位掩模的簇；
基于所述第一光功率确定当前系统能量；
至少基于所述系统能量和先前最小系统能量来确定是否更新所述相位掩模；以及
响应于确定更新所述相位掩模，
基于所述相位掩模的簇更新所述相位掩模。
10.根据权利要求9所述的方法，其中，生成所述相位掩模进一步包括：
响应于确定不更新所述相位掩模，
将所述当前系统能量和所述先前最小系统能量中的较小者识别为所述近似基态。
11.根据权利要求9所述的方法，其中，生成所述相位掩模进一步包括：
识别热能；
基于所述当前系统能量和所述先前最小系统能量计算系统能量的变化；以及
基于所述系统能量的变化和所述热能确定玻尔兹曼概率，并且
其中，确定是否更新所述相位掩模进一步基于所述玻尔兹曼概率。
12.根据权利要求9所述的方法，其中，确定是否更新所述相位掩模包括：
确定反馈迭代计数小于阈值；
确定更新所述相位掩模。
13.一种光学计算系统，包括：
光源，其配置为产生泵浦光束；
光学调制器，其配置为接收调制掩模并基于所述调制掩模调制所述泵浦光束，以生成调制光束；
非线性介质，其配置为将所述调制光束的一部分转换为二次谐波(SH)光束并产生包括所述SH光束和所述泵浦光束的未转换部分的输出；以及
分色镜，其配置为接收所述非线性介质的所述输出并将所述SH光束和所述泵浦光束的所述未转换部分去耦合；
第一相机，其配置为捕获所述泵浦光束的所述未转换部分并生成第一强度矩阵；
第二相机，其配置为捕获所述泵浦光束的所述未转换部分并生成第二强度矩阵；以及
控制器，其配置为基于所述第一强度矩阵和所述第二强度矩阵中的至少一者生成更新调制掩模，以传输到所述光学调制器。
14.根据权利要求13所述的光学计算系统，其中，所述光学调制器包括空间光调制器，所述空间光调制器配置为调制所述泵浦光束的小波的相位，以生成所述调制光束，并且
其中，所述调制掩模是二进制相位图案。
15.根据权利要求14所述的光学计算系统，其中，所述空间光调制器配置为通过基于所述调制掩模将小波中的每个小波的相位编码为第一相位或第二相位，来调制所述泵浦光束的所述小波的相位，并且
其中，所述第一相位和所述第二相位相差180度。
16.根据权利要求13所述的光学计算系统，其中，所述光学调制器包括数字微镜器件，所述数字微镜器件配置为调制所述泵浦光束的小波的强度，以生成所述调制光束。
17.根据权利要求13所述的光学计算系统，其中，所述控制器配置为通过以下方式生成所述更新调制掩模：
将相机状态确定为所述第一强度矩阵、所述第二强度矩阵或者所述第一强度矩阵和所述第二强度矩阵的加权和；以及
基于所述调制掩模和所述相机状态来确定所述更新调制掩模。
18.根据权利要求17所述的光学计算系统，其中，基于所述调制掩模和所述相机状态确定所述更新调制掩模包括：
将所述更新调制掩模确定为所述调制掩模与所述相机状态和反馈步长的乘积之间的差。
19.根据权利要求18所述的光学计算系统，其中，所述泵浦光束包括由所述光源生成的高斯脉冲序列中的高斯脉冲，并且
其中，所述反馈步长对应于所述高斯脉冲序列的周期性。
20.根据权利要求18所述的光学计算系统，其中，所述控制器配置为在生成所述更新调制掩模之前适应性地调整所述反馈步长。
21.一种光学计算系统，包括：
光学放大器，其配置为生成具有平坦相位波前的泵浦光束；
简并光学参量放大器，其配置为基于所述泵浦光束在小波中生成双波长信号；
滤光器，其配置为从所述信号中过滤所述泵浦光束；
光学分束器，其配置为读出所述信号的一部分；
非线性光学晶体，其配置为基于所述信号生成二次谐波光；
第一双波段反射器，其配置为分离所述信号和所述二次谐波光；
光学接收器，其配置为将所述二次谐波光耦合到所述光学放大器中；和
第二双波段反射器，其配置为将所述信号和所述二次谐波光组合到所述简并光学参量放大器中。
22.根据权利要求21所述的光学计算系统，进一步包括：
第一傅立叶透镜，其配置为将所述信号聚焦到所述非线性光学晶体中，以生成所述二次谐波光；和
第二傅立叶透镜，其位于所述第一双波段反射器和所述第二双波段反射器之间，并配置为使所述信号散焦。
23.根据权利要求21所述的光学计算系统，其中，所述第一双波段反射器和所述第二双波段反射器中的每一者都包括分色镜或滤光器。
24.根据权利要求21所述的光学计算系统，其中，随着所述信号连续穿过所述光学计算系统，所述光学放大器的光学增益增加，直到光学参量振荡发生。
25.根据权利要求21所述的光学计算系统，还包括位于所述简并光学参量放大器和所述非线性光学晶体之间的第一光学调制器，所述第一光学调制器配置为在所述信号到达所述非线性光学晶体之前调制所述信号，以改变所述光学计算系统的有效四体相互作用哈密顿函数。
26.根据权利要求24所述的光学计算系统，还包括位于所述非线性光学晶体和所述光学接收器之间的第二光学调制器，所述第二光学调制器配置为调制所述二次谐波光，以改变所述光学计算系统的有效四体相互作用哈密顿函数。
27.一种光学计算系统，包括：
光学放大器，其配置为生成具有平坦相位波前的放大信号；
非线性光学器件，其配置为接收所述放大信号并基于所述放大信号通过二次谐波生成来生成泵浦光束；
简并光学参量放大器，其配置为基于所述泵浦光束在小波中生成双波长信号；
滤光器，其配置为从所述信号中过滤所述泵浦光束；
第一光学分束器，其配置为读出所述信号的一部分；
第二光学分束器，其配置为将所述信号的第一部分引导至用于调制所述信号的该部分的光学调制器，并将所述信号的第二部分引导至分色镜；以及
光学接收器，其配置为将所述信号的调制第一部分耦合到所述光学放大器中，以生成所述放大信号，
其中，所述分色镜配置为将所述信号的第二部分与所述泵浦光束组合到所述简并光学参量放大器中。
28.根据权利要求27所述的光学计算系统，其中，所述光学调制器包括空间光调制器或数字微镜器件，并且
其中，所述非线性光学器件包括非线性光学晶体或光波导。
CN202180040984.4A 2020-05-13 2021-05-13 具有多体相互作用和全对全连接的超级伊辛仿真器 Pending CN115698897A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
US202063024257P 2020-05-13 2020-05-13
US63/024,257 2020-05-13
PCT/US2021/032338 WO2021231794A1 (en) 2020-05-13 2021-05-13 Super ising emulator with multi-body interactions and all-to-all connections

Publications (1)

Publication Number Publication Date
CN115698897A true CN115698897A (zh) 2023-02-03

Family

ID=76284206

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202180040984.4A Pending CN115698897A (zh) 2020-05-13 2021-05-13 具有多体相互作用和全对全连接的超级伊辛仿真器

Country Status (7)

Country Link
US (1) US20230185160A1 (zh)
EP (1) EP4150430A1 (zh)
JP (1) JP2023527135A (zh)
KR (1) KR20230011342A (zh)
CN (1) CN115698897A (zh)
CA (1) CA3178484A1 (zh)
WO (1) WO2021231794A1 (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2021077021A1 (en) 2019-10-18 2021-04-22 The Trustees Of The Stevens Institute Of Technology Systems and methods for quantum-secured, private-preserving computations

Also Published As

Publication number Publication date
KR20230011342A (ko) 2023-01-20
JP2023527135A (ja) 2023-06-27
EP4150430A1 (en) 2023-03-22
CA3178484A1 (en) 2021-11-18
US20230185160A1 (en) 2023-06-15
WO2021231794A1 (en) 2021-11-18

Similar Documents

Publication Publication Date Title
US20200327403A1 (en) All optical neural network
Rafayelyan et al. Large-scale optical reservoir computing for spatiotemporal chaotic systems prediction
Dong et al. Optical reservoir computing using multiple light scattering for chaotic systems prediction
Goy et al. Low photon count phase retrieval using deep learning
Pierangeli et al. Scalable spin-glass optical simulator
Kumar et al. Large-scale Ising emulation with four body interaction and all-to-all connections
Ong et al. Photonic convolutional neural networks using integrated diffractive optics
Zuo et al. Scalability of all-optical neural networks based on spatial light modulators
Huisman et al. Programming balanced optical beam splitters in white paint
Calvanese Strinati et al. All-optical scalable spatial coherent Ising machine
Tezak et al. A coherent perceptron for all-optical learning
US12020150B2 (en) Optical neural network
Zhou et al. Ultrafast dynamic machine vision with spatiotemporal photonic computing
Suprano et al. Dynamical learning of a photonics quantum-state engineering process
CN115698897A (zh) 具有多体相互作用和全对全连接的超级伊辛仿真器
Hazineh et al. D-flat: A differentiable flat-optics framework for end-to-end metasurface visual sensor design
Ouyang et al. On-demand photonic Ising machine with simplified Hamiltonian calculation by phase encoding and intensity detection
Vorontsov Problems of large neurodynamics system modeling: Optical synergetics and neural networks
Lu et al. Enhanced phase retrieval using nonlinear dynamics
Nakajima et al. Physical deep learning with biologically plausible training method
Chen et al. Physics-aware complex-valued adversarial machine learning in reconfigurable diffractive all-optical neural network
Hermans et al. Towards trainable media: Using waves for neural network-style training
Bacot et al. Revisiting time reversal and holography with spacetime transformations
Rozenberg et al. Designing nonlinear photonic crystals for high-dimensional quantum state engineering
Ramesh et al. A spatial-photonic Ising machine to solve the two-way number-partitioning problem

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination