QETLAB njohnsto_qetlab https://qetlab.com/Main_Page MediaWiki 1.45.3 first-letter Media Special Talk User User talk QETLAB QETLAB talk File File talk MediaWiki MediaWiki talk Template Template talk Help Help talk Category Category talk 0 1 35997 35532 2023-08-01T01:54:15Z Nathaniel 1 LaTeX fix 35997 wikitext text/x-wiki __NOTOC____NOEDITSECTION__== QETLAB: A MATLAB Toolbox for Quantum Entanglement == '''QETLAB''' ('''Q'''uantum '''E'''ntanglement '''T'''heory '''LAB'''oratory) is a MATLAB toolbox for exploring quantum entanglement theory. While there are many quantum information theory toolboxes that allow the user to perform basic operations such as the [[PartialTranspose|partial transposition]], new tests are constantly discovered. The goal of QETLAB is to remain up-to-date and contain an ever-growing catalogue of separability criteria, positive maps, and related functions of interest. Furthermore, QETLAB is designed to work well both with full matrices and with large sparse matrices, and makes use of many advanced techniques based on semidefinite programming. {{DownloadBox}} ==Getting started: some introductory examples== ===Pure state entanglement and the Schmidt decomposition=== The heart of QETLAB is a suite of functions that make working with quantum entanglement quick and easy. Working with pure state entanglement is a breeze, as demonstrated by the following code snippet, which generates a random pure state (according to uniform spherical measure/Haar measure) and computes its Schmidt decomposition: <code><nowiki>>> d = 3; % this example will be in 3 \otimes 3 dimensions >> phi = </nowiki>[[RandomStateVector]]<nowiki>(d^2); >> [s,a,b] = </nowiki>[[SchmidtDecomposition]]<nowiki>(phi) s = 0.9159 0.3951 0.0709 a = -0.8062 -0.5411 0.2392 -0.2685 - 0.0526i 0.0305 + 0.2555i -0.8357 + 0.4006i 0.1580 + 0.5002i -0.3599 - 0.7152i -0.2814 + 0.0681i b = -0.1401 - 0.3521i 0.0829 - 0.7178i -0.3997 - 0.4176i -0.2513 + 0.1492i -0.3837 + 0.3479i 0.0874 - 0.7992i 0.5520 - 0.6828i -0.3670 + 0.2737i -0.1378 - 0.0208i </nowiki></code> We could verify that this Schmidt decomposition is correct as follows: <code><nowiki>>> phi2 = s(1)*kron(a(:,1),b(:,1)) + s(2)*kron(a(:,2),b(:,2)) + s(3)*kron(a(:,3),b(:,3)); >> norm(phi - phi2) ans = 5.8378e-016 </nowiki></code> However, QETLAB has a [[TensorSum|built-in function]] that acts as the inverse of the Schmidt decomposition, to make this process easier: <code><nowiki>>> phi3 = </nowiki>[[TensorSum]]<nowiki>(s,a,b); >> norm(phi - phi3) ans = 5.8378e-016 </nowiki></code> ===Mixed state entanglement and detecting bound entanglement=== QETLAB can also handle mixed state entanglement without trouble. To demonstrate these features, let's start by creating a bound entangled state and then detecting its entanglement. One well-known way to create bound entangled states is via unextendible product bases (UPBs)<ref name="BDM99">C.H. Bennett, D.P. DiVincenzo, T. Mor, P.W. Shor, J.A. Smolin, and B.M. Terhal. Unextendible product bases and bound entanglement. ''Phys. Rev. Lett.'' 82, 5385&ndash;5388, 1999. E-print: [http://arxiv.org/abs/quant-ph/9808030 arXiv:quant-ph/9808030]</ref><ref name="DMS03">D.P. DiVincenzo, T. Mor, P.W. Shor, J.A. Smolin, and B.M. Terhal. Unextendible product bases, uncompletable product bases and bound entanglement. ''Commun. Math. Phys.'' 238, 379&ndash;410, 2003. E-print: [http://arxiv.org/abs/quant-ph/9908070 arXiv:quant-ph/9908070]</ref>. Let's first generate the "Tiles" UPB in <math>\mathbb{C}^3 \otimes \mathbb{C}^3</math> and construct a bound entangled state from it: <code><nowiki>>> u = </nowiki>[[UPB]]<nowiki>('Tiles') % generate the "Tiles" unextendible product basis u = 0.7071 0 0 0 0.3333 -0.7071 0 0 0 0.3333 0 0.7071 0 0 0.3333 0 0 0 0.7071 0.3333 0 0 0 0 0.3333 0 -0.7071 0 0 0.3333 0 0 0 -0.7071 0.3333 0 0 0.7071 0 0.3333 0 0 -0.7071 0 0.3333 >> rho = eye(9) - u*u'; % rho is the projection onto the orthogonal complement of the UPB vectors >> rho = rho/trace(rho); % scale rho to have trace 1 </nowiki></code> Now <tt>rho</tt> is a bound entangled state. We can verify this in QETLAB as follows: <code><nowiki>>> </nowiki>[[IsPPT]]<nowiki>(rho) % this will show that rho has positive partial transpose ans = 1 >> </nowiki>[[IsSeparable]]<nowiki>(rho) % this will show that rho is entangled Determined to be entangled via the realignment criterion. Reference: K. Chen and L.-A. Wu. A matrix realignment method for recognizing entanglement. Quantum Inf. Comput., 3:193-202, 2003. ans = 0 </nowiki></code> QETLAB can also carry out many related functions, such as [[SymmetricExtension|finding entanglement witnesses that detect entanglement in a given state]], [[IsBlockPositive|determining whether or not a given operator is an entanglement witness]], and [[IsUPB|deciding whether or not a given set of vectors form an unextendible product basis]]. Furthermore, there are numerous examples scattered throughout the documentation that make learning how to carry out these operations a breeze. ==References== <references /> hvljigul1oxjikuv0r63ovxar6iexff