|
23 | 23 | "Block encoding is a well known technique in quantum computing used to embed non-unitary matrices on a quantum computer that only allows for unitary evolution. \n", |
24 | 24 | "\n", |
25 | 25 | "***Definition:*** \n", |
26 | | - " Let $ a, n, m \\in \\mathbb{N} $ such that $ m = a + n $. A $ m $ -qubit unitary $ U $ is said to be an $(\\alpha, a)$ -block-encoding of an $ n $ -qubit operator $ A $ if \n", |
| 26 | + " Let $ a, n, m \\in \\mathbb{N} $ such that $ m = a + n $. An $ m $ -qubit unitary $ U $ is said to be an $(\\alpha, a)$ -block-encoding of an $ n $ -qubit operator $ A $ if \n", |
27 | 27 | "\n", |
28 | 28 | "$$\n", |
29 | 29 | "\\tilde{A} = \\left( \\langle 0 |^{\\otimes a} \\otimes I_n \\right) U \\left( |0 \\rangle^{\\otimes a} \\otimes I_n \\right)\n", |
|
32 | 32 | "\n", |
33 | 33 | "where, $A = \\alpha \\tilde{A}$ . The parameters $(\\alpha, a)$ represent the *subnormalization factor* (which adjusts for encoding matrices of any norm), and the *number of ancilla qubits* used in the block-encoding scheme respectively. \n", |
34 | 34 | "\n", |
35 | | - "Efficiently block encoding arbitrary matrices is a very difficult problem and this task is not trivial even for well structured and sparse matrices. The paper by Sturm et al. 2025 (https://arxiv.org/abs/2509.02429), provides for efficient quantum circuits for block encoding $N-$ dimensional Laplacians, with periodic boundary conditions along all dimensions. This is useful in many applications, especially ones involving problems of linear algebra. Moreover, given an efficient block encoding of a matrix $\\tilde{A}$, its possible to efficiently construct a block encoding of certain polynomials of $\\tilde{A}$ through quantum singular value transformations (QSVT). \n", |
| 35 | + "Efficiently block encoding arbitrary matrices is a very difficult problem and this task is not trivial even for well structured and sparse matrices. The paper by Sturm et al. 2025 (https://arxiv.org/abs/2509.02429), provides efficient quantum circuits for block encoding $N-$ dimensional Laplacians, with periodic boundary conditions along all dimensions. This is useful in many applications, especially ones involving linear algebra problems. Moreover, given an efficient block encoding of a matrix $\\tilde{A}$, it's possible to efficiently construct a block encoding of certain polynomials of $\\tilde{A}$ through quantum singular value transformations (QSVT). \n", |
36 | 36 | "\n", |
37 | 37 | "### **Notebook contents**\n", |
38 | | - "- ##### Block encoding circuits for D- dimensional Laplacian matrix with periodic boundary conditions along all dimnesions.\n", |
| 38 | + "- ##### Block encoding circuits for D-dimensional Laplacian matrix with periodic boundary conditions along all dimensions.\n", |
39 | 39 | "___" |
40 | 40 | ] |
41 | 41 | }, |
|
62 | 62 | "\\tag{1}\n", |
63 | 63 | "\\end{equation}\n", |
64 | 64 | "\n", |
65 | | - "We impose **periodic boundary conditions**, meaning the function repeats itself across opposite sides of the domain:\n", |
| 65 | + "Consider a function $u$ on which the operator $L_D$ acts. We impose **periodic boundary conditions**, meaning the function repeats itself across opposite sides of the domain:\n", |
66 | 66 | "$$\n", |
67 | | - "v(x^{(0)},\\ldots,0,\\ldots,x^{(D-1)}) \n", |
| 67 | + "u(x^{(0)},\\ldots,0,\\ldots,x^{(D-1)}) \n", |
68 | 68 | "= \n", |
69 | | - "v(x^{(0)},\\ldots,1,\\ldots,x^{(D-1)}).\n", |
| 69 | + "u(x^{(0)},\\ldots,1,\\ldots,x^{(D-1)}).\n", |
70 | 70 | "$$\n", |
71 | 71 | "\n", |
72 | 72 | "This ensures the domain “wraps around” in **every dimension**.\n", |
73 | 73 | "\n", |
74 | 74 | "\n", |
| 75 | + "\n", |
75 | 76 | "To approximate the operator on a computer, we replace each interval with a **uniform grid** (*we assume equidistant grid points as per the paper*). \n", |
76 | 77 | "\n", |
77 | 78 | "For one dimension, choose\n", |
|
109 | 110 | "j = \\ket{j^{(D-1)}}\\otimes \\cdots \\otimes\\ket{j^{(1)}}\\otimes\\ket{j^{(0)}},\n", |
110 | 111 | "$$\n", |
111 | 112 | "\n", |
| 113 | + "- ***Example***\n", |
| 114 | + "\n", |
| 115 | + " If $n = 3$ (so $N = 2^n = 8$), the index $j^{(d)} = 5$ is represented in the computational basis as\n", |
| 116 | + " $\n", |
| 117 | + " \\ket{5} = \\ket{101},\n", |
| 118 | + " $\n", |
| 119 | + " corresponding to the binary representation of the integer $5$.\n", |
| 120 | + " Explicitly, the basis state $\\lvert 5 \\rangle$ corresponds to the $2^n$-dimensional column vector\n", |
| 121 | + " $$\n", |
| 122 | + " \\lvert 5 \\rangle\n", |
| 123 | + " =\n", |
| 124 | + " \\begin{pmatrix}\n", |
| 125 | + " 0 \\\\\n", |
| 126 | + " 0 \\\\\n", |
| 127 | + " 0 \\\\\n", |
| 128 | + " 0 \\\\\n", |
| 129 | + " 0 \\\\\n", |
| 130 | + " 1 \\\\\n", |
| 131 | + " 0 \\\\\n", |
| 132 | + " 0\n", |
| 133 | + " \\end{pmatrix},\n", |
| 134 | + " $$\n", |
| 135 | + " where the single nonzero entry appears at position $5$, following zero-based indexing.\n", |
| 136 | + "\n", |
| 137 | + "\n", |
| 138 | + "\n", |
112 | 139 | "This ordering turns the entire grid $\\Omega_{D,h}$ into a vector of $\\mathbf{\\Omega_{D,h}}$ size $N^D= 2^{nD}$ where,\n", |
113 | 140 | "\n", |
114 | 141 | "$$\n", |
|
583 | 610 | }, |
584 | 611 | { |
585 | 612 | "cell_type": "code", |
586 | | - "execution_count": 5, |
| 613 | + "execution_count": 3, |
587 | 614 | "id": "11", |
588 | 615 | "metadata": {}, |
589 | 616 | "outputs": [], |
|
0 commit comments