Abstract
在高反射镜的热噪声是几种类型的高精度干涉实验,旨在达到标准量子极限或机械系统冷却到其量子基态的一个主要障碍。这是,例如预计将在最敏感的频段限制的情况下,未来的引力波观测站,其引力波信号的敏感性,通过原子振动镜群众。一所追求的有前途的方法来克服这个限制是采用高阶拉盖尔 - 高斯(LG)光束代替传统使用的基本模式。由于他们的更均匀的光强度分布的光束平均更有效地在热驱动的镜子,这反过来又降低了检测出的激光在反射镜的位置的不确定性的波动。
我们展示了一个有前途的方法来生成高阶LG光束的帮助下,衍射光学元件的一个基本的高斯光束整形。我们表明,与传统的传感和控制技术是已知的用于稳定基波激光束的高阶LG模可以被纯化,并稳定在一个相对较高的水平一样好。了一套诊断工具,使我们能够控制和定制生成的LG光束的属性。这使我们生产的LG电子束报告日期纯度最高。证明高阶LG干涉测量技术与标准和使用标准球面光学模式的兼容性,使得应用在下一代的高精度干涉他们的理想人选。
Introduction
在过去的几十年里,高精度干涉实验推向最终灵敏度制度量子效应开始发挥了决定性的作用。在这些目前和未来的实验,如激光冷却机械振荡器1,光学陷阱一代纠缠测试群众镜2,3,量子非拆迁干涉的4刚性腔5,稳频激光器,引力波探测6 ,7,8,研究人员正面临着许多限制基本面和技术面的噪声源。最严重的问题之一是谐振器反射镜的干涉的设置,这是造成由热激发的原子构成的反射镜基板和反射涂层的反射镜7,8,9的热噪声。这样的效果,也被称为布朗运动 ,将导致一个不确定性的相位反射光的任何测试的群众,所以会表现为根本干涉仪的输出噪声限制。例如,在项目设计先进的引力波天线,如高级,高级LIGO的处女座,和爱因斯坦望远镜的灵敏度是有限的这种类型的噪声最敏感的区域的观测频段10,11,12。
实验物理学家在社区努力做出持续的努力,以尽量减少这些噪声贡献,并提高他们的仪器的灵敏度。在镜子布朗噪声的特定情况下,缓解的一个方法是采用更大范围内的光束点的大小,因为目前使用的标准的基本HG 00束测试质量表面的更有效的表面的随机运动的更大范围内的平均束13,14。已经显示出扩展的反射镜的热噪声的功率谱密度逆高斯光束的大小的反射镜基板和镜子9平方成反比。然而,束斑变大时,较大部分的光功率丢失的反射面的边缘。如果使用比通常使用的HG 00束(见例如图1)的径向光强分布更均匀的光束, 布朗热噪声电平可以减小而不会增加这种类型的损失。在所有的更均匀的光束类型,已经提出了新版本的高精度的干涉,例如梅萨梁或圆锥形模式13,14,最有前途的是高阶LG与目前使用的球形光束由于其潜在的兼容性反射镜表面15。例如,二进制中子星螺旋系统的检出率 - 这被认为是最有前途的天体物理源第一毛重检测-离子约为系数2或更大的16上面的最小量的成本在目前正在建设10,11的第二代干涉仪的设计修改,可以增强。除了 热噪声的好处,更广泛的高阶LG光束的强度分布(见,例如, 图2)已被证明,以减轻热像差的光学干涉仪内的幅度。这将减少热补偿系统在何种程度上依赖在未来的实验,以达到设计的敏感性19。
我们也调查,成功证明了可行性产生LG光束成功运作最好的灵敏度16,18,19,20,21,22 GW干涉所需的纯度和稳定性水平。该方法结合在不同的领域物理和光学齐全的开发技术和专业知识h为新一代的高稳定性,低噪声的单模激光束23,利用空间光调制器,衍射光学元件的操纵空间分布的光束18,22,24,25,26,和使用旨在进一步纯化和稳定的激光的光学谐振腔27的感测,控制和稳定化的先进技术。该方法已成功地证明,在实验室实验中,出口为测试在大型原型干涉仪20,用于产生在高的激光功率的LG模式高达80 W 21。在这篇文章中,我们提出高阶LG光束的方法的细节,并讨论所产生的光束的表征和验证的方法论。另外,在步骤4中概述了数值研究了一种用于与非完美反射镜19的空腔。
Protocol
序言:在本协议部分中,我们假设,一个纯粹的,噪音低,功率稳定的基模高斯光束,例如含图3中所示的标准设置,通过商业的Nd:YAG激光器产生连续波在1064纳米的波长的光;法拉第隔离器(FI),以避免对激光源的背面反射的光的一个电光调制器(EOM)的红外光的光的相位调制。由此产生的光束被注入到一个三角形的光学谐振腔,其中通过有效的控制回路27的激光频率的光功率稳定,而谐振腔提供不必要的波束形状的空间滤波。
上面描述的设置, 如图3所示,是一个传统的实验安排,科学装置中使用的精度要求低噪声激光稳定测量1-8。协议部分说明这个基本模式高斯光束可以有效地转换成高阶拉盖尔 - 高斯光束相媲美的表演,如果不相同,在纯度,噪声和稳定。这是通过所示的装置在图4中 ,其设计,建设,和操作在下面的章节描述的装置。在这项工作中提出了在这个例子中,所生成的模式将是一个所有33。然而,值得强调的是,该技术具有普遍的有效性和描述的协议,适用于任何所需的高阶LG模式。
1。光学模式转换为基本模式激光束的优化转换高阶LG梁的设计和原型
一个基本模式束转换成高阶LG光束的相位调制档案中的要求是复制的相位CROβ-边的所需的LG模式,该模式将被印上的波阵面的入射光束26通过成比例的相移。两种类型的模式转换器以这种方式工作:空间光调制器(SLM) - 计算机控制的液晶显示器,其可以控制到压印相的像素上的入射光转变 - 衍射相位板 - 蚀刻玻璃基板所需在传输中产生的相移由故意不同厚度的玻璃元件。 SLM的是灵活的,但缺乏稳定性和效率,而相位板是稳定的,有效的,但缺乏灵活性。因此,我们建议使用最初的研究和原型设计和使用的相位板的长期操作的SLM。
最佳的转换依赖于精确的选择参数(腰围的大小和位置)的形状的光束。因此,注入到模式转换器之前,最初的基本模式时,必须将其特征在于,和它的参数重新塑造,以提供最佳的转换的匹配 - 该操作被称为“模式匹配的”。
- 拿起束从图3中描述的基本模式设置。
- 使用光束分析仪配备了实时图像分析软件,可以衡量光束半径沿光路。一旦已获得足够一套半径(一般至少10个数据点都需要一个良好的质量结果),适合测量半径和提取束腰的规模和它的位置。
- 建立所需的光束半径转换点。使用大束数毫米大小的顺序才能使用相位转换区域的充分程度。
- 选择一组沿光路的透镜和它们的位置,将重新塑造成所需的入射光束参数(腰围的大小和位置)。对于对齐的目的,可以很方便地放置的模式转换器吨的入射光束的腰身。
- 重复步骤1.2和1.4,通过连续的透镜位置调整,直到已获得所期望的光束模式转换参数。
- SLM模式转换沿入射光束路径,并注入到SLM束。反射型SLM,我们建议使用一个小的入射角度,5度或以下的订单。在生成的光束入射角大,会造成散光,打破了LG模式圆柱对称。
- 应用相位分布到SLM的液晶显示装置 - ,相位将其转换为所需的高阶LG电子束的横截面。 33的LG模式,这是目前的研究在未来的毛重检测器16的应用程序,在图5中所示的例子中的相位调制模式。
- 选择适当的相位图案的大小(对应的光束的相位图案的大小)的大小的基础上的注入TED束表1到9的顺序,使用数值模拟28导出最佳的光束大小比为LG模式的列表。或者,找到最佳的光束是通过改变施加到空间光调制器的相位图案的大小和分析所得到的光束的图像的图像尺寸比实验。
- 观察反射光束从SLM使用的CCD照相机,在距离中的一个或多个瑞利范围相差从SLM。仔细对准,以优化在CCD上的光束的对称图像的SLM。
一些注入光的相位调制装置的交互期间,仍然未调制的相位调制电平由于量化。此未转化的光转换后的光束沿同一轴线上,损坏所需的相位调制效应。为了规避这个问题,可以对LG模式转换相位图像覆盖闪耀光栅轮廓。调制后的升搭载LG模式相位轮廓的洞察力将闪耀光栅偏转,而未灯光,不与基材互动的,将继续不受干扰。这将导致在两种类型之间的光束的空间分离。
- 熊熊结构重叠先前SLM产生的相位分布。 LG模与方位指数升> 0的相位图案将有一个叉形光栅的特点,如在图6的例子中看到。
- 优化燃烧的角度,使得入一阶的衍射角的光束的发散角大于。继续进行,直到较高的衍射级之间的一个合理的分离结果(使用作为外圈本身的直径一样大的连续梁的外圈之间的分离)。
- 一旦实现最佳的转换模式,进行相板的制造。这些是商业上阿瓦伊标签和制造,以满足广泛的定制要求。使用与SLM在优化过程中所获得的结果,确定最佳的相位变换的图案被蚀刻到的相位板。可选步骤:应用中的至少一个的相位板的表面上的防反射涂层,以尽量减少散射的光反射回的光功率向所述激光源和分散。
2。相板,模式转换和纯度增强操作
- 更换空间光调制器的相位板。至于在SLM,它是方便的,以将其定位在注入的基本模式要转换的光束的腰身。
- 小心地将相位板等的初始光束的相位板是垂直于梁,相对于光束的中心的相结构。
- 传播的光束透过的相位板,直到较高的差异分离raction订单发生。梁与梁卡可以很容易地可视。
- 当达到一个足够的“好”的分离(在步骤1.12中所述),遮掩较高的衍射级次光束的孔为中心的主衍射阶。
所讨论的相位板设计,不能调节幅度以及相位意味着他们将不能转换成所需的模式所有传入的根本束。其结果是一个组合梁的显性其他高阶模式轻微的强度,如在图7中示出的背景下进行所需的LG光束。为了在空间上过滤掉不想要的LG模式和增强模式的纯度,可以将转换后的光束的光学谐振腔的注入。这样的空腔,也可作为一个“模式选择”只允许特定的光学模式的传输,相对于光的波长的谐振腔长度根据。
- 设计莫清洁腔。如果在其实施简单,使用由两个反光镜的线性腔结构, 如图4所示,在这种反射镜是平的(通常是输入镜)和其他的反射镜(输出)是凹的。这提供了实施的光稳定性和简单。甲行之有效的具体设计是一个的输出反射镜的曲率半径为1米,而反射镜的反射表面之间的距离是21厘米29。在这种情况下,最佳的输入光束半径大约是365微米在腰部,位于平面镜的反射表面。
- 选择腔镜的反射率来确定的腔技巧。使用低为了几百技巧,有良好的抑制不必要的模式订单,而不会引入大的扭曲,由于耦合退化模式(见第4步)。最好是使用具有相同反射率的反射镜,以最大限度地提高模腔通过。
- 使用RIGID间隔为两腔的支持反映机械振动增强免疫力。胶的反射镜在隔板上,并插入一个压电环元件之间的两个反射镜的一个,与间隔件的纵向长度的控制和稳定的目的,以允许谐振腔长度的微观调整。
- 模式匹配的光束所产生的相位板模式清洁器腔体本征模式。的LG光束的束仿形可以不被执行使用相同的工具用于根本模式束,因此,记录的强度分布的光束与一个放置在不同的位置沿光束路径的CCD相机和分析采用量身定制合身的记录,图像脚本,它可以识别主导的LG模式和光束半径估计在给定的位置30。该光束的强度分布拟合程序的一个例子示于图8。
- 一旦有足够的光束直径已成为连接测量(通常,至少10个数据点所需要的良好的质量的结果),适合测得的半径和外推的束腰直径和它的位置。一个很好的光束轮廓看起来像一个如图9所示。在1.2和1.4,选择透镜,并重复该过程,直到最优波束的大小和位置被发现,2.7,2.8和2.9中描述。一旦实现模式匹配,注入所产生的光束进入模式吸尘器腔,确保输入(平坦)的反射面的反射镜是否正确放置在腰部注入光束。
- 到空腔优化注入的光束的对准,同时,通过移动反射镜与压电扫描腔长,监视发送波束。
- 使用由模式清洁器腔作为调查对LG光束所产生的相位板的模式内容的谐振腔长度(也称为空腔扫描)的函数的透射光的测量,前夕ntually评估的相位板本身的转换效率。
- 通过检查的CCD图像识别有关的寄生模式。评估这种模式下的功率通过在光电二极管的信号其振幅和计算准确的整体光束的模式的内容。的测量结果和确切形式的内容可以进行复制,数值模拟21。的分析结果在图10中给出一个很好的例子,模式含量结果示于表2。
一旦成束的最佳取向模式清洁器腔被实现,分析已注入的光束模式的内容,可以最终实现方式清洁和增强的复合LG束纯度。一斤德雷弗霍尔锁定方案27可用于稳定腔长到所需的谐振模式。光传输模式吸尘器CAVITY可以读取由一个光电二极管,它可以提供错误的信号所必需的控制回路,控制腔长。
- 腔长锁定到发送的诊断所产生的光束,并限定其纯度与CCD摄像机的空腔产生的光束的档案的主共振和记录图像。
3。生成的LG光束诊断和表征
在这个实验中,两个主要的属性定义一个'好'光束质量的高精度干涉测量的成功实施:束功率和光束的纯度。使用上实现的基本模式的光束中,如上所述的相同的控制技术,可以保存其他相关的属性如频率或功率稳定。
- LG电子束功率测量,通过激光功率计。注意梁剪切:LG光束具有较大EXTEN锡永相比,传统的高斯光束,它可能会超过大多数商业工具的敏感面积的尺寸。最高权力明显。
- 通过比较与理论光束剖面评估生成的LG电子束的纯度。要做到这一点,通过CCD相机探查和拍摄照片的光束强度估计它的光束半径,推导出理论束振幅剖面比较测量。通过评估的纯度平方内积理论和测量振幅分布。高纯度的建议。
两个重要的品质因数是有用的的整个模式转换过程中的质量评价:相板的转换效率和整体设置。
- 为了评估CONVERS离子相板效率,遵循腔扫描程序步骤2.11和2.12。
- 评估的生成所需的LG光束的功率与注入基模高斯光束的功率之间的比率设置为整体的转换效率。高转换效率显然是可取的。
4。注入大干涉模拟研究
这个协议的应用之一是调查他们的引力波探测器使用LG光束。这些长基线高精度干涉。基线需要比较大的反射镜和光束大小。然而,这提高了不完善的光学元件的影响,特别是当使用高次模。本节介绍了一种基于模拟的方法进行调查的行为高阶LG模式在现实探测器。
- 选择光场建模仿真工具,在interferomet的呃为了测试高阶LG模式。仿真软件应该能够在设置不完善的影响(错位,模式不匹配,镜子数字错误等 )建模的模式内容的光束。一个例子是仿真工具FINESSE 28。
- 建立一个模型,一个真正的探测器使用选定的模拟工具。在高级LIGO的情况下,这是一个双循环迈克耳孙干涉仪臂腔的Fabry-Perot。这些初始模拟的目的是验证了模型的可靠性,假设完美的光学系统。
- 测试模型与完善的基本模式横梁。为了验证该模型的可靠性,这应该允许再生的列表,在实际的检测器进行的实验程序,如:对臂腔,腔的扫描循环中的电源,如预期的数字误差信号,并检查,并角度和长度的干涉仪,其子系统的控制S通过传感和控制方案。进一步的模拟应包括干涉仪引力波信号的响应。一旦模拟按预期执行,该模型可以适于高阶LG模。
- 测试模型与完美的LG33梁:适应干涉仪的设计,使用LG模式。这就需要减少上的空腔,其中可以实现通过改变的反射镜的曲率半径的反射镜的光束的大小。一旦模型已被改编为LG模式,应重复进行的测试4.3的新的输入束。如果在完美的光学系统的情况下,结果应该是非常相似,使用HG 00(见实例19)。
使用更高阶的光束引入了简并性的光学腔,因为有几个不同的光束形状,争取主导地位。高斯模式的光学谐振腔的谐振该命令适用于所有模式。一个HG00 0阶模式是唯一的模式,因此所有其他模式抑制。例如,LG 33模式10模式9阶的一个,所有这些都将在干涉仪得到加强。镜面扭曲,总是存在于真实的干涉夫妇到其他的事件模式。如果这些新的模式是相同的顺序,作为入射光束在臂腔增强,导致高度扭曲的循环光束。这最终会恶化仪器的灵敏度。
- 设置一个现实的干涉模型:将腔反射镜的表面数字的真实数据。此数据的属性几何高度或反射率的反射镜的表面,如一个“地图”的形式,在图11中的高级LIGO镜子看到一个例子。包括这些影响后,高阶模的性能,应进行调查,特别是在在检测器的输出和误差信号中的多个过零点的可能性的对比度缺陷。在这些领域中,高阶模式预计将表现逊于HG 00。
- 模拟子系统:为了更好地了解模型中存在的退化效果,模拟子系统简,例如高级LIGO的法布里 - 珀罗臂腔。这些子系统模拟产生腔内扫描和错误信号,以识别任何的频率分裂和流通领域的检测,可以分析其模式内容。
- 镜面要求:获得更严格的要求的情况下镜面的平整度4.6步的结果显示这将使执行高阶LG光束不可能在其他模式下的频率分裂或电源不可接受的程度。对于这一点,分析耦合等造成的直接跨阶其中可以实现数字或使用的解析近似19的表面。使用这种方法来识别任何特定的镜子的形状,是造成大量相同的顺序输入光束和模式之间的耦合。通过比较这些结果与模拟,估算为这些形状的镜像要求,对于一个特定的循环束纯度。最后模拟满量程的干涉仪模型修改为新规格的镜像地图,示出了改进对比度的缺陷和频率分裂。
Materials
Name | Company | Catalog Number | Comments |
The experimental apparatus discussed in this paper requires the following types of instruments: | |||
Instrument | |||
Solid state Laser source, Nd:YAG 1064 nm CW laser | Quantity: 1 | ||
Faraday Isolator | Quantity: 1 | ||
Electro-Optic Modulator (EOM) | Quantity: 1 | ||
CCDcamera beam profiler | Quantity: 1 | ||
Lenses | Quantity: depending on apparatus design | ||
Steering Mirrors | Quantity: depending on apparatus design | ||
Aperture | Quantity: 1 | ||
High reflectivity mirrors (for normal incidence) | Quantity: 2 | ||
Piezoelectric ring | Quantity: 1 | ||
Cavity spacer | Quantity: 1 | ||
Photodiodes and related control electronics | Quantity: 1 or more, depending on apparatus design | ||
Spatial light modulator | Quantity: 1 Holoeye LCR-2500 |
||
All the above instruments are commercially available and no particular specification is required. We leave the choice of the most suitable instruments to the experimenter’s discretion. | |||
For the interest of the experimenter interested in reproducing the protocol, we recommend the following tools used in our experiment: | |||
Tools | |||
Innolight OEM 300NE, 1064 nm, 300 mW | Laser Source: | ||
SIMTOOLs | Software for data analysis, available at www.gwoptics.org/simtools/ | ||
FINESSE | Software for optical simulations, www.gwoptics.org/finesse/ | ||
Finally, the phase plate employed in the present experiment was manufactured by Jenoptik GmbH, based on a custom design provided by the Authors. |
References
- Cohadon, P. F., Heidmann, A., Pinard, M. Cooling of a Mirror by Radiation Pressure. Physical Review Letters. 83, 3174-3177 (1999).
- Corbitt, T., et al. An All-Optical Trap for a Gram-Scale Mirror. Physical Review Letters. 98, 150802 (2007).
- Müller-Ebhardt, H., Rehbein, H., Schnabel, R., Danzmann, K., Chen, Y. Entanglement of Macroscopic Test Masses and the Standard Quantum Limit in Laser Interferometry. Physical Review Letters. 100, 013601 (2008).
- Kimble, H. J., Levin, Y., Matsko, A. B., Thorne, K. S., Vyatchanin, S. P. Conversion of conventional gravitational-wave interferometers into quantum nondemolition interferometers by modifying their input and output optics. Physical Review D. 65, 022002 (2001).
- Numata, K., Kemery, A., Camp, J. Thermal-Noise Limit in the Frequency Stabilization of Lasers with Rigid Cavities. Physical Review Letters. 93, 250602 (2004).
- Aufmuth,, Danzmann, K.
Gravitational wave detectors. New Journal of Physics. 7, 202 (2005). - Harry, G. M., et al. Thermal noise in interferometric gravitational wave detectors due to dielectric optical coatings. Classical and Quantum Gravity. 19, 897-917 (2002).
- Crooks, D. R. M., et al. Excess mechanical loss associated with dielectric mirror coatings on test masses in interferometric gravitational wave detectors. Classical and Quantum Gravity. 19, 883-896 (2002).
- Yu, L. Internal thermal noise in the LIGO test masses: A direct approach. Physical Review D. 57, 659-663 (1998).
- Harry, G. M. the LIGO Collaboration. Advanced LIGO: the next generation of gravitational wave detectors. Classical and Quantum Gravity. 27, 084006 (2010).
- Accadia, T., Acernese, F., Antonucci, F., Astone, P., Ballardin, G., Barone, F., et al.
Status of the VIRGO project. Classical and Quantum Gravity. 28, 114002 (2011). - Sathyaprakash, B., Abernathy, M., Acernese, F., Ajith, P., Allen, B., Amaro-Seoane, P., et al.
Scientific objectives of Einstein Telescope. Classical and Quantum Gravity. 29 (12), 124013 (2012). - Bondarescu, M., Thorne, K. S. New family of light beams and mirror shapes for future LIGO interferometers. Physical Review D. 74, 082003 (2006).
- Bondarescu, M., Kogan, O., Chen, Y. Optimal light beams and mirror shapes for future LIGO interferometers. Physical Review D. 78, 082002 (2008).
- Tournefier Mours, E., Vinet, J. Y. Thermal noise reduction in interferometric gravitational wave antennas: using high order TEM modes. Classical and Quantum Gravity. 23, 5777 (2006).
- Chelkowski, S., Hild, S., Freise, A. Prospects of higher-order Laguerre-Gauss modes in future gravitational wave detectors. Physical Review D. 79, 122002 (2009).
- Vinet, J. Y. Reducing thermal effects in mirrors of advanced gravitational wave interferometric detectors. Classical and Quantum Gravity. 24, 3897 (2007).
- Fulda, P., Kokeyama, K., Chelkowski, S., Freise, A. Experimental demonstration of higher-order Laguerre-Gauss mode interferometry. Physical Review D. 82, 012002 (2010).
- Bond, C., Fulda, P., Carbone, L., Kokeyama, K., Freise, A. Higher order Laguerre-Gauss mode degeneracy in realistic, high finesse cavities. Physical Review D. 84, 102002 (2011).
- Sorazu, B., et al. Experimental test of higher-order Laguerre-Gauss modes in the 10 m Glasgow prototype interferometer. Class. Quantum Grav. 30, 035004 (2013).
- Carbone, L., Bogan, C., Fulda, P., Freise, A., Willke, W. Generation of High-Purity Higher-Order Laguerre-Gauss Beams at High Laser Power. Physical Review Letters. 110, 251101 (2013).
- Fulda, P., et al. Phaseplate design for Laguerre-Gauss mode conversion. , In Preparation (2013).
- Willke, B. Stabilized lasers for advanced gravitational wave detectors. Laser & Photonics Rev. , 1-15 (2010).
- Granata, M., Buy, C., Ward, R., Barsuglia, M. Higher-Order Laguerre-Gauss Mode Generation and Interferometry for Gravitational Wave Detectors. Physical Review Letters. 105, 231102 (2010).
- Matsumoto, N., Ando, T., Inoue, T., Ohtake, Y., Fukuchi, N., Hara, T. Generation of high-quality higher-order Laguerre Gaussian beams using liquid-crystal-on-silicon spatial light modulators. J. Opt. Soc. Am. 25, 1642-1651 (2008).
- Kennedy, S. A., Szabo, M. J., Teslow, H., Porterfield, J. Z., Abraham, E. R. I. Creation of Laguerre-Gaussian laser modes using diffractive optics. Physical Review A. 66, 043801 (2002).
- Black, E. D. An introduction to Pound-Drever-Hall laser frequency stabilization. American Journal of Physics. 69, 79-87 (2001).
- Freise, A., Heinzel, G., Lueck, H., Schilling, R., Willke, B., Danzmann, K. Frequency-domain interferometer simulation with higher-order spatial. Class.Quant.Grav. 21, (2004).
- Uehara, N. Mode cleaner for the Initial LIGO 10 Watt Laser. LIGO internal report. 12, (1997).
- Simtools, a collection of Matlab tools for optical simulations [Internet]. , Available from: http://www.gwoptics.org/simtools/ (2013).