quantum factor graphs closing the box operation and
play

Quantum Factor Graphs: Closing-the-Box Operation and Variational - PowerPoint PPT Presentation

Dec 07, 2022 •50 likes •1.84k views

Quantum Factor Graphs: Closing-the-Box Operation and Variational Approaches End-of-Second-Year Oral Exam Michael X. CAO Department of Information Engineering, CUHK August 29, 2016 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29,

  1. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In general, we associate the factorization � g ( x ) � f a ( x a ) a ∈F to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F given by E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  2. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  3. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . F We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  4. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . F V We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  5. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . F V We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  6. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . F V We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  7. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs , or classical factor graph (CFG), describe factorizations. Definition 1 In Example 1, In general, we associate the factorization � g ( x ) � f a ( x a ) f A X 1 a ∈F f B X 2 to a factor graph with variable node set V , function node set F , and edge set E ⊆ V × F X 3 given by f C X 4 E = { ( i , a ) ∈ V × F : i ∈ ∂ a } . F V We call g the global function , f the local function/factors . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 2

  8. ✶ Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs are popular in describing (large-scale) probability systems. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 3

  9. ✶ Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs are popular in describing (large-scale) probability systems. Example 2 (LDPC Code described by CFG) y 1 x 1 p Y 1 | X 1 + y 2 x 2 p Y 2 | X 2 + y 3 x 3 p Y 3 | X 3 . . . . . . . . . . . . + y n x n p Y n | X n Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 3

  10. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs are popular in describing (large-scale) probability systems. x i , y i ∈ F 2 ∀ i ∈ { 1 , · · · , n } ; Example 2 (LDPC Code described by CFG)     � f + ( x ) � ✶ x i = 0  . y 1 x 1  i incoming p Y 1 | X 1 + y 2 x 2 p Y 2 | X 2 + y 3 x 3 p Y 3 | X 3 . . . . . . . . . . . . + y n x n p Y n | X n Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 3

  11. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs are popular in describing (large-scale) probability systems. x i , y i ∈ F 2 ∀ i ∈ { 1 , · · · , n } ; Example 2 (LDPC Code described by CFG)     � f + ( x ) � ✶ x i = 0  . y 1 x 1  i incoming p Y 1 | X 1 + A problem of interest: Calculate the marginal distribution y 2 x 2 p Y 2 | X 2 + of X i given fixed { y i } n i =1 . y 3 x 3 p Y 3 | X 3 . . . . . . . . . . . . + y n x n p Y n | X n Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 3

  12. Factor Graphs/Preliminaries Factor Graphs and Normal Factor Graphs Factor graphs are popular in describing (large-scale) probability systems. x i , y i ∈ F 2 ∀ i ∈ { 1 , · · · , n } ; Example 2 (LDPC Code described by CFG)     � f + ( x ) � ✶ x i = 0  . y 1 x 1  i incoming p Y 1 | X 1 + A problem of interest: Calculate the marginal distribution y 2 x 2 p Y 2 | X 2 + of X i given fixed { y i } n i =1 . � � g y ( x ) = p Y i = y i | X i ( x i ) · f k ( x ) y 3 x 3 p Y 3 | X 3 . . i k . . . . ∝ p Y = y | X ( x ) . . . . . . + � g y ( x ) ∝ p Y = y | X i ( x i ) x j , j � = i y n x n p Y n | X n Symbol-wise ML decoding Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 3

  13. Factor Graphs/Preliminaries Partition Function/Sum Definition 3 (Partition Function/Sum) In many applications, we are interested in calculating summations like (or similar to) � � � Z � g ( x ) = f a ( x ∂ a ) , x x a ∈F which is defined to be the partition function/sum of a CFG. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 4

  14. Factor Graphs/Preliminaries Partition Function/Sum Definition 3 (Partition Function/Sum) In many applications, we are interested in calculating summations like (or similar to) � � � Z � g ( x ) = f a ( x ∂ a ) , x x a ∈F which is defined to be the partition function/sum of a CFG. Example 1 (continue) f A X 1 g = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) . f B X 2 � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) . X 3 x 1 , x 2 , x 3 , x 4 f C X 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 4

  15. Factor Graphs/Preliminaries Partition Function/Sum Definition 3 (Partition Function/Sum) In many applications, we are interested in calculating summations like (or similar to) � � � Z � g ( x ) = f a ( x ∂ a ) , x x a ∈F which is defined to be the partition function/sum of a CFG. Example 1 (continue) f A X 1 g = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) . f B X 2 � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) . X 3 x 1 , x 2 , x 3 , x 4 f C X 4 In general, calculation of Z is NP hard. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 4

  16. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A X 1 � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 f B X 2 X 3 f C X 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  17. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A X 1 � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 f B X 2 X 3 f C X 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  18. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A X 1 � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f B � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 f C X 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  19. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� � � � f A X 1 = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f B � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 f C X 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  20. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f AB � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 f C X 4 � � � ˆ � = f AB ( x 3 ) · f C ( x 3 , x 4 ) x 3 x 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  21. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f AB � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 f C X 4 � � � ˆ � = f AB ( x 3 ) · f C ( x 3 , x 4 ) x 3 x 4 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  22. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f AB � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 ˆ f C � � � ˆ � = f AB ( x 3 ) · f C ( x 3 , x 4 ) x 3 x 4 � � � ˆ · ˆ = f AB ( x 3 ) f C ( x 3 ) x 3 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  23. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 ˆ f AB � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 ˆ f C � � � ˆ � = f AB ( x 3 ) · f C ( x 3 , x 4 ) x 3 x 4 � � � ˆ · ˆ = f AB ( x 3 ) f C ( x 3 ) x 3 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  24. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Example 1 Sum over local variables first... � Z = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 1 , x 2 , x 3 , x 4 �� � �� f A � � � = f A ( x 1 ) · f B ( x 1 , x 2 , x 3 ) · f C ( x 3 , x 4 ) x 3 x 1 x 2 x 4 Z � � �� � � f A ( x 1 ) · ˆ � = f B ( x 1 , x 3 ) · f C ( x 3 , x 4 ) X 3 x 3 x 1 x 4 � � � ˆ � = f AB ( x 3 ) · f C ( x 3 , x 4 ) x 3 x 4 � � � ˆ · ˆ = f AB ( x 3 ) f C ( x 3 ) = Z x 3 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 5

  25. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Closing-the-box in general In general, “closing-the-box” means to replace the boxes with the result of the summing over the interior variable(s). z 1 w 1 f 1 w 3 x f 3 y z 3 f 2 z 2 w 2 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 6

  26. Factor Graphs/Preliminaries Partition Sum and the Closing-the-box Operations Closing-the-box in general In general, “closing-the-box” means to replace the boxes with the result of the summing over the interior variable(s). z 1 z 1 w 1 w 1 f 1 w 3 w 3 � f 1 · f 2 · f 3 x f 3 x , y y z 3 z 3 f 2 z 2 z 2 w 2 w 2 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 6

  27. Factor Graphs/Preliminaries Sum-Product Algorithm on a tree Example 1 f A X 1 f B X 2 X 3 f C X 4 Closing the boxes from the inner ones to outer ones will yield the partition sum Z. Distributive Law on R Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 7

  28. Factor Graphs/Preliminaries Sum-Product Algorithm on a tree Example 1 Sum-Product Algorithm for Trees Require: Acyclic factor graph G = ( F , V , E ); root r ∈ V ; height of the tree h � 0. f A Ensure: Partition sum Z X 1 1: for d = h − 1 , · · · , 0 do for all i ∈ V d-step reachable a from r f B X 2 2: do X 3 Let f ( i ) be the parent factor b of i ; 3: f ( i ) ← � � a ∈ ∂ i f a ( x i ); f C X 4 4: x i end for 5: 6: end for Closing the boxes from the inner ones to 7: Z ← f ( r ) . outer ones will yield the partition sum Z. Distributive Law on R a i.e., there exists a path connecting r and i passing through d factors. b i.e. the unique factor node that is both on the path from r to i and also adjacent to i . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 7

  29. Factor Graphs/Preliminaries Sum-Product Algorithm on a tree Example 1 Sum-Product Algorithm for Trees Require: Acyclic factor graph G = ( F , V , E ); root r ∈ V ; height of the tree h � 0. f A Ensure: Partition sum Z X 1 1: for d = h − 1 , · · · , 0 do for all i ∈ V d-step reachable a from r f B X 2 2: do X 3 Let f ( i ) be the parent factor b of i ; 3: f ( i ) ← � � a ∈ ∂ i f a ( x i ); f C X 4 4: x i end for 5: 6: end for Closing the boxes from the inner ones to 7: Z ← f ( r ) . outer ones will yield the partition sum Z. Distributive Law on R a i.e., there exists a path connecting r and i passing through d factors. Message-Passing Algorithm b i.e. the unique factor node that is both on the path from r to i and also adjacent to i . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 7

  30. Factor Graphs/Preliminaries Sum-Product Algorithm on a tree Example 1 Sum-Product Algorithm for Trees Require: Acyclic factor graph G = ( F , V , E ); root r ∈ V ; height of the tree h � 0. f A Ensure: Partition sum Z X 1 ← 1: for d = h − 1 , · · · , 0 do for all i ∈ V d-step reachable a from r ← f B X 2 2: → do X 3 Let f ( i ) be the parent factor b of i ; 3: f ( i ) ← � ← � a ∈ ∂ i f a ( x i ); f C X 4 4: x i end for 5: 6: end for Closing the boxes from the inner ones to 7: Z ← f ( r ) . outer ones will yield the partition sum Z. Distributive Law on R a i.e., there exists a path connecting r and i passing through d factors. Message-Passing Algorithm b i.e. the unique factor node that is both on the path from r to i and also adjacent to i . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 7

  31. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ← � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ← � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  32. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ← � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ← � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  33. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ∝ � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ∝ � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  34. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ∝ � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ∝ � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  35. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ∝ � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ∝ � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until convergence Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  36. Factor Graphs/Preliminaries Sum-Product Algorithm as a Message Passing Algorithm Sum-Product Algorithm Require: Factor graph G = ( F , V , E ); Ensure: ??? 1: for all ( i , a ) ∈ E do m i → a ← ✶ ; 2: m a → i ← ✶ ; 3: 4: end for 5: repeat for all ( i , a ) ∈ E do 6: m i → a ( x i ) ∝ � c ∈ ∂ i \ a m a → i ( x i ); 7: end for 8: for all ( i , a ) ∈ E do 9: m a → i ( x i ) ∝ � x ∂ a \ i f a ( x ∂ a ) · � j ∈ ∂ a \ i m j → a ( x j ); 10: end for 11: 12: until convergence Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 8

  37. Factor Graphs/Preliminaries Sum-Product Algorithm and the Variational Approach In acyclic case, it will always converge. b i ( x i ) � � a ∈ ∂ i m a → i ( x i ) ∝ � x V\ i g ( x ); b a ( x ∂ a ) � f a ( x ∂ a ) · � i ∈ ∂ a m i → a ( x j ) ∝ � x V\ ∂ a g ( x ). Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 9

  38. Factor Graphs/Preliminaries Sum-Product Algorithm and the Variational Approach In acyclic case, it will always converge. b i ( x i ) � � a ∈ ∂ i m a → i ( x i ) ∝ � x V\ i g ( x ); b a ( x ∂ a ) � f a ( x ∂ a ) · � i ∈ ∂ a m i → a ( x j ) ∝ � x V\ ∂ a g ( x ). In general case, if it converges, then: [Yedidia et al., 2005] The above { b i } i ∈V and { b a } a ∈F correspond to the interior stationary points of the constrained Bethe approximation problem : Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 9

  39. Factor Graphs/Preliminaries Sum-Product Algorithm and the Variational Approach In acyclic case, it will always converge. b i ( x i ) � � a ∈ ∂ i m a → i ( x i ) ∝ � x V\ i g ( x ); b a ( x ∂ a ) � f a ( x ∂ a ) · � i ∈ ∂ a m i → a ( x j ) ∝ � x V\ ∂ a g ( x ). In general case, if it converges, then: [Yedidia et al., 2005] The above { b i } i ∈V and { b a } a ∈F correspond to the interior stationary points of the constrained Bethe approximation problem : � � min F Bethe ( b a ) a ∈F , ( b i ) i ∈V � � � � � − b a ( x ∂ a ) log f a ( x ∂ a ) − H ( b a ) + ( d i − 1) · H ( b i ) x ∂ a a ∈F a ∈F i ∈V s . t . b a probability on x ∂ a , b i probability on x i , ∀ a ∈ F , ∀ i ∈ V � b i ( x i ) = b a ( x ∂ a ) ∀ x i , ∀ ( i , a ) ∈ E x ∂ a \ i Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 9

  40. Quantum Factor Graphs (QFGs) Outline Factor Graphs/Preliminaries 1 Quantum Factor Graphs (QFGs) 2 Closing-the-box Operations on QFGs 3 Variational Approach on QFGs 4 Numerical Result of QSPA 5 Conclusion & Outlook 6 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 10

  41. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) Definition 4 ([Leifer and Poulin, 2008]) A quantum Factor graph (QFG) ( V , F , E ) with local factors { ρ a } describes the “factorization” �� � � ρ � ρ a = exp log( ρ a ) (1) , a ∈F a ∈F where, for each a ∈ F , positive definite operator ρ a is an operator on � i ∈ ∂ a H i . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 10

  42. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) Definition 4 ([Leifer and Poulin, 2008]) A quantum Factor graph (QFG) ( V , F , E ) with local factors { ρ a } describes the “factorization” �� � � ρ � ρ a = exp log( ρ a ) (1) , a ∈F a ∈F where, for each a ∈ F , positive definite operator ρ a is an operator on � i ∈ ∂ a H i . Example 5 ρ A ∈ L + ( H 1 ) H 1 H 2 ρ A ρ B ρ C ρ B ∈ L + ( H 1 ⊗ H 2 ) ρ C ∈ L + ( H 2 ) A QFG describing ρ = ρ A ⊙ ρ B ⊙ ρ C . Here, L + ( H ) stands for the set of all positive semi-definite operators on the Hilbert space H . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 10

  43. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  44. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. By the Lie Product formula, we have � n 1 1 � ρ A ⊙ ρ B = lim ρ A ρ n n . (3) B n →∞ Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  45. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. By the Lie Product formula, we have � n 1 1 � ρ A ⊙ ρ B = lim ρ A ρ n n . (3) B n →∞ Equation (3) can be used to generalize the ⊙ product to PSD operators. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  46. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. By the Lie Product formula, we have � n 1 1 � ρ A ⊙ ρ B = lim ρ A ρ n n . (3) B n →∞ Equation (3) can be used to generalize the ⊙ product to PSD operators. Properties of ⊙ Associativity: ( ρ A ⊙ ρ B ) ⊙ ρ C = ρ A ⊙ ( ρ B ⊙ ρ C ); Commutativity: ρ A ⊙ ρ B = ρ B ⊙ ρ A ; Closeness: ρ A ⊙ ρ B is positive (semi) definite if ρ A , ρ B are positive (semi) definite. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  47. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. By the Lie Product formula, we have � n 1 1 � ρ A ⊙ ρ B = lim ρ A ρ n n . (3) B n →∞ Equation (3) can be used to generalize the ⊙ product to PSD operators. Properties of ⊙ Associativity: ( ρ A ⊙ ρ B ) ⊙ ρ C = ρ A ⊙ ( ρ B ⊙ ρ C ); Commutativity: ρ A ⊙ ρ B = ρ B ⊙ ρ A ; Closeness: ρ A ⊙ ρ B is positive (semi) definite if ρ A , ρ B are positive (semi) definite. �L + ( H ) , ⊙� (or �L ++ ( H ) , ⊙� ) is an Abelian group. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  48. Quantum Factor Graphs (QFGs) Quantum Factor Graphs (QFGs) For ρ A , ρ B ∈ L ++ ( H ), define [Warmuth, 2005] ρ A ⊙ ρ B � exp � � log( ρ A ) + log( ρ B ) , (2) where exp and log denote the operator exponential and the operator natural logarithm, respectively. By the Lie Product formula, we have � n 1 1 � ρ A ⊙ ρ B = lim ρ A ρ n n . (3) B n →∞ Equation (3) can be used to generalize the ⊙ product to PSD operators. ρ A ∈ L + ( H 1 ) Properties of ⊙ � H 1 � = H 2 ρ B ∈ L + ( H 2 ) Associativity: ( ρ A ⊙ ρ B ) ⊙ ρ C = ρ A ⊙ ( ρ B ⊙ ρ C ); Commutativity: ρ A ⊙ ρ B = ρ B ⊙ ρ A ; H 1 → H 3 ← H 2 Closeness: ρ A ⊙ ρ B is positive (semi) definite if ρ A , ρ B are positive (semi) definite. ρ A ⊙ ρ B � ˜ ρ A ⊙ ˜ ρ B . �L + ( H ) , ⊙� (or �L ++ ( H ) , ⊙� ) is an Abelian group. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 11

  49. Quantum Factor Graphs (QFGs) Quantum Factor Graphs Example 5: continue H 1 = C 2 H 2 = C 2 ρ A ρ B ρ C A QFG describing ρ = ρ A ⊙ ρ B ⊙ ρ C . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 12

  50. Quantum Factor Graphs (QFGs) Quantum Factor Graphs Example 5: continue H 1 = C 2 H 2 = C 2 ρ A ρ B ρ C A QFG describing ρ = ρ A ⊙ ρ B ⊙ ρ C . Suppose H 1 = H 2 = C 2 , and   0 � � � � +3 − 1 1 +3 − 1   ρ A = , ρ B =  , ρ C = .   − 1 +3 1 − 1 +3  0 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 12

  51. Quantum Factor Graphs (QFGs) Quantum Factor Graphs Example 5: continue H 1 = C 2 H 2 = C 2 ρ A ρ B ρ C A QFG describing ρ = ρ A ⊙ ρ B ⊙ ρ C . Suppose H 1 = H 2 = C 2 , and   0 � � � � +3 − 1 1 +3 − 1   ρ A = , ρ B =  , ρ C = .   − 1 +3 1 − 1 +3  0 We have,   9 − 3 − 3 1 − 3 9 1 − 3   ρ = ρ A ⊙ ρ B ⊙ ρ C =  .   − 3 1 9 − 3  1 − 3 − 3 9 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 12

  52. Quantum Factor Graphs (QFGs) Quantum Partition-Sum Problem Definition 6 (Partition Function/Sum) In a number of applications, we are interested in calculating �� � � �� �� Z � Tr ( ρ ) = Tr ρ a = Tr exp log( ρ a ) , a ∈F a ∈F which is defined to be the partition function/sum of a QFG. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 13

  53. Quantum Factor Graphs (QFGs) Quantum Partition-Sum Problem Definition 6 (Partition Function/Sum) In a number of applications, we are interested in calculating �� � � �� �� Z � Tr ( ρ ) = Tr ρ a = Tr exp log( ρ a ) , a ∈F a ∈F which is defined to be the partition function/sum of a QFG. In general calculation of the partition function/sum of a QFG is NP hard. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 13

  54. Closing-the-box Operations on QFGs Outline Factor Graphs/Preliminaries 1 Quantum Factor Graphs (QFGs) 2 Closing-the-box Operations on QFGs 3 Variational Approach on QFGs 4 Numerical Result of QSPA 5 Conclusion & Outlook 6 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  55. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in CFG X 1 X 2 f A f B f C Z = � x 1 , x 2 f A ( x 1 ) f B ( x 1 , x 2 ) f C ( x 2 ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  56. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in CFG X 1 X 2 f A f B f C Z = � x 1 , x 2 f A ( x 1 ) f B ( x 1 , x 2 ) f C ( x 2 ) = � x 1 f A ( x 1 ) � x 2 f B ( x 1 , x 2 ) f C ( x 2 ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  57. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in CFG = ⇒ Applying distributive law X 1 X 2 f A f B f C Z = � x 1 , x 2 f A ( x 1 ) f B ( x 1 , x 2 ) f C ( x 2 ) = � x 1 f A ( x 1 ) � x 2 f B ( x 1 , x 2 ) f C ( x 2 ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  58. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in QFG H 1 H 2 ρ A ρ B ρ C Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  59. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in QFG H 1 H 2 ρ A ρ B ρ C ? � � Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  60. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in QFG = ⇒ Distributive law over (partial) trace H 1 H 2 ρ A ρ B ρ C ? � � Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  61. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in QFG = ⇒ Distributive law over (partial) trace H 1 H 2 ρ A ρ B ρ C ? � � Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) However, in general, � � � � Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 Tr 2 ( ρ A ⊙ ρ B ⊙ ρ C ) � = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  62. Closing-the-box Operations on QFGs Closing-the-box operations and partial trace functions Closing the box in QFG = ⇒ Distributive law over (partial) trace H 1 H 2 ρ A ρ B ρ C ? � � Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) However, in general, � � � � Tr( ρ A ⊙ ρ B ⊙ ρ C ) = Tr 1 Tr 2 ( ρ A ⊙ ρ B ⊙ ρ C ) � = Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) . Example 7 � +1 − 1 Let H 1 = H 2 � C 2 . Suppose ρ A � 1 � , ρ B ⊙ ρ C � diag (0 , 1 , 1 , 0) . In 2 · − 1 +1 � � this case, Tr( ρ A ⊙ ρ B ⊙ ρ C ) = 0 and Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) = 1. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 14

  63. Closing-the-box Operations on QFGs Cases when factors are close to identity matrix Oftentimes, we can still have an approximate closing-the-box rule. Lemma 8 We have bounds Given ρ A ∈ L ++ ( H 1 ); � − 1 � � S κ ( ρ A ) ρ B ∈ L ++ ( H 1 ⊗ H 2 ); Tr( ρ A ⊙ ρ B ⊙ ρ C ) ρ C ∈ L ++ ( H 2 ); � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) κ ( · ) � 1 is the condition number function; � S � � κ ( ρ A ) . S ( · ) is the Specht ratio function defined as 1 S ( r ) � ( r − 1) · r r − 1 . e · log r Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 15

  64. Closing-the-box Operations on QFGs Cases when factors are close to identity matrix Oftentimes, we can still have an approximate closing-the-box rule. Lemma 8 We have bounds Given ρ A ∈ L ++ ( H 1 ); � − 1 � � S κ ( ρ A ) ρ B ∈ L ++ ( H 1 ⊗ H 2 ); Tr( ρ A ⊙ ρ B ⊙ ρ C ) ρ C ∈ L ++ ( H 2 ); � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) κ ( · ) � 1 is the condition number function; � S � � κ ( ρ A ) . S ( · ) is the Specht ratio function defined as 1 S ( r ) � ( r − 1) · r r − 1 . e · log r The proof utilizes the Golden–Thompson inequality [Bourin and Seo, 2007]. Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 15

  65. Closing-the-box Operations on QFGs Cases when factors are close to identity matrix Oftentimes, we can still have an approximate closing-the-box rule. Lemma 8 We have bounds Given Specht Ratio Function � − 1 � 8 � S κ ( ρ A ) 6 Tr( ρ A ⊙ ρ B ⊙ ρ C ) � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) 4 � S � � κ ( ρ A ) . 2 0 0 20 40 60 80 100 condition number r The proof utilizes the Golden–Thompson inequality [Bourin and Seo, 2007]. Considering ρ A ≈ I , Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 15

  66. Closing-the-box Operations on QFGs Cases when factors are close to identity matrix Oftentimes, we can still have an approximate closing-the-box rule. Lemma 8 We have bounds Given Specht Ratio Function � − 1 � 8 � S κ ( ρ A ) 6 Tr( ρ A ⊙ ρ B ⊙ ρ C ) � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) 4 � S � � κ ( ρ A ) . 2 0 0 20 40 60 80 100 condition number r The proof utilizes the Golden–Thompson inequality [Bourin and Seo, 2007]. Considering ρ A ≈ I , we expect to have � � Tr( ρ A ⊙ ρ B ⊙ ρ C ) ≈ Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 15

  67. Closing-the-box Operations on QFGs Type-1 Approximation when ρ A is “close” to identity matrix I H 1 H 2 H 1 H 2 ρ A ρ B ρ C ρ A ρ B ρ C Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  68. Closing-the-box Operations on QFGs Type-1 Approximation when ρ A is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ A ρ B ρ C ρ A ρ B ρ C Z = Tr( ρ A ⊙ ρ B ⊙ ρ C ) � � Tr 1 ρ A ⊙ Tr 2 ( ρ B ⊙ ρ C ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  69. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 ρ 1 , 2 ρ 1 ρ 1 , 2 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) � � Tr 1 ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  70. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  71. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  72. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � Lemma (Type-1 Approximation) Given X ∈ L H ( H 1 ), and Y ∈ L H ( H 1 ⊗ H 2 ), for t close to 0, we have = ( I + tX ) ⊙ Tr 2 ( I + tY ) + O ( t 3 ) . � � Tr 2 ( I + tX ) ⊙ ( I + tY ) (4) Theorem (Type-1 Approximation) Following the same setup, we have + O ( t 4 ) . � � � � Tr ( I + tX ) ⊙ ( I + tY ) = Tr 1 ( I + tX ) ⊙ Tr 2 ( I + tY ) (5) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  73. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � “Linear” close to I : i.e., ρ 1 = I + tX and ρ 1 , 2 = I + tY . Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  74. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � “Linear” close to I : i.e., ρ 1 = I + tX and ρ 1 , 2 = I + tY . Taylor Series Expansion: Tr( ρ 1 ⊙ ρ 1 , 2 ) = � � Tr 1 ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) = Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  75. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � “Linear” close to I : i.e., ρ 1 = I + tX and ρ 1 , 2 = I + tY . Taylor Series Expansion: X + Y ) + t 2 · Tr( ˜ XY + Y ˜ X ) + t 3 · 0 Tr( ρ 1 ⊙ ρ 1 , 2 ) = 1 + t · Tr( ˜ 2 � X 2 Y 2 � XY ˜ ˜ XY − ˜ Tr + t 4 · + · · · 12 = 1 + t · Tr 1 ( X + Tr 2 ( Y )) + t 2 · Tr 1 � X Tr 2 ( Y ) + Tr 2 ( Y ) X � � � Tr 1 ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) 2 X Tr 2 ( Y ) X Tr 2 ( Y ) − X 2 Tr 2 ( Y ) 2 � + t 4 · Tr 1 � + · · · 12 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  76. Closing-the-box Operations on QFGs Type-1 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � “Linear” close to I : i.e., ρ 1 = I + tX and ρ 1 , 2 = I + tY . Taylor Series Expansion: X + Y ) + t 2 · Tr( ˜ XY + Y ˜ X ) + t 3 · 0 Tr( ρ 1 ⊙ ρ 1 , 2 ) = 1 + t · Tr( ˜ 2 � X 2 Y 2 � XY ˜ ˜ XY − ˜ Tr + t 4 · + · · · 12 = 1 + t · Tr 1 ( X + Tr 2 ( Y )) + t 2 · Tr 1 � X Tr 2 ( Y ) + Tr 2 ( Y ) X � � � Tr 1 ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) 2 X Tr 2 ( Y ) X Tr 2 ( Y ) − X 2 Tr 2 ( Y ) 2 � + t 4 · Tr 1 � + · · · 12 Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 16

  77. Closing-the-box Operations on QFGs Type-2 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 17

  78. Closing-the-box Operations on QFGs Type-2 Approximation when ρ 1 or ρ 1 , 2 is “close” to identity matrix I ≈ H 1 H 2 H 1 H 2 ρ 1 , 2 ρ 1 , 2 ρ 1 ρ 1 Z = Tr( ρ 1 ⊙ ρ 1 , 2 ) Tr 1 � ρ 1 ⊙ Tr 2 ( ρ 1 , 2 ) � Lemma (Type-2 Approximation) Given X ∈ L H ( H 1 ), and Y ∈ L H ( H 1 ⊗ H 2 ), for t close to 0, we have e tX ⊙ e tY � = e tX ⊙ Tr 2 ( e tY ) + O ( t 3 ) . � Tr 2 (4) Theorem (Type-2 Approximation) Following the same setup, we have e tX ⊙ e tY � e tX ⊙ Tr 2 ( e tY ) + O ( t 4 ) . � � � Tr = Tr 1 (5) Michael X. CAO (IE@CUHK) Quantum Factor Graphs August 29, 2016 17

Download Presentation
Download Policy: The content available on the website is offered to you 'AS IS' for your personal information and use only. It cannot be commercialized, licensed, or distributed on other websites without prior consent from the author. To download a presentation, simply click this link. If you encounter any difficulties during the download process, it's possible that the publisher has removed the file from their server.

Recommend

F ACTOR graph [1], [2], or more often referred to as the classical factor graph (CFG) in this

F ACTOR graph [1], [2], or more often referred to as the classical factor graph (CFG) in this

SECOND YEAR REPORT 1 Quantum Factor Graphs: Closing-the-Box Operation and Variational Approaches Michael X. CAO, PhD Pre-Candidacy Student Abstract Factor graph model is a popular statistical graphical model, where a number of practical

624 views • 18 slides
Generalization of Factor Graphs and Belief Propagation for Quantum Information Science

Generalization of Factor Graphs and Belief Propagation for Quantum Information Science

Generalization of Factor Graphs and Belief Propagation for Quantum Information Science End-of-First-Year Oral Exam Michael X. CAO Department of Information Engineering, CUHK September 25, 2015 Michael X. CAO (IE@CUHK) Quantum Factor Graph

1.45k views • 121 slides
F ACTOR graph is a popular graphical model to represent where the inner product between A , B

F ACTOR graph is a popular graphical model to represent where the inner product between A , B

FIRST YEAR REPORT 1 Generalization of Factor Graphs and Belief Propagation for Quantum Information Science Michael X. CAO Abstract Factor graph is a useful tool to represent probability p u A u B q B b B b A systems. In this article, we

406 views • 8 slides
Post-quantum Security of the CBC, CFB, OFB, Modes of Operation. CTR, and XTS Modes of Operation.

Post-quantum Security of the CBC, CFB, OFB, Modes of Operation. CTR, and XTS Modes of Operation.

Post-quantum Security of the CBC, CFB, OFB, CTR, and XTS Post-quantum Security of the CBC, CFB, OFB, Modes of Operation. CTR, and XTS Modes of Operation. Mayuresh Anand Motivation Mayuresh Anand, Ehsan Ebrahimi Targhi, Gelo Noel Tabia,

235 views • 22 slides
Colourings of (0 , m )-graphs and the switching operation Gary MacGillivray University of

Colourings of (0 , m )-graphs and the switching operation Gary MacGillivray University of

Colourings of (0 , m )-graphs and the switching operation Gary MacGillivray University of Victoria Joint work with Chris Duffy, Ben Tremblay and J. Maria Warren GT2015, Nyborg, DK, August 26, 2015 Prelude Harary (1953) used graphs with 2

881 views • 55 slides
Leveraging Operation-Aware UREQA: Error Rates for Effective Quantum Circuit Mapping on NISQ-Era

Leveraging Operation-Aware UREQA: Error Rates for Effective Quantum Circuit Mapping on NISQ-Era

Leveraging Operation-Aware UREQA: Error Rates for Effective Quantum Circuit Mapping on NISQ-Era Quantum Computers Tirthak Patel Baolin Li Rohan Basu Roy Devesh Tiwari Quantum Computing is Coming! What is a Qubit (Quantum Bit)? A classical

347 views • 30 slides
Approximate Knowledge Compilation by Online Collapsed Importance Sampling Tal Friedman and Guy

Approximate Knowledge Compilation by Online Collapsed Importance Sampling Tal Friedman and Guy

Approximate Knowledge Compilation by Online Collapsed Importance Sampling Tal Friedman and Guy Van den Broeck Motivation Factor Graphs: 1 Motivation Factor Graphs: Great! But asking queries is hard 2 Motivation Factor Graphs: 3

629 views • 34 slides
Periodic quantum graphs in magnetic fields Differential Operators on Graphs and Waveguides

Periodic quantum graphs in magnetic fields Differential Operators on Graphs and Waveguides

Periodic quantum graphs in magnetic fields Differential Operators on Graphs and Waveguides Graz 2019 Simon Becker Cambridge Consider a lattice Z 2 or the honeycomb lattice. Consider a lattice Z 2 or the honeycomb lattice. We study quantum

1.03k views • 69 slides
Quantum Graphs! Priyanga Ganesan November 11, 2020 ASU C*-Seminar Priyanga Ganesan

Quantum Graphs! Priyanga Ganesan November 11, 2020 ASU C*-Seminar Priyanga Ganesan

Quantum Graphs! Priyanga Ganesan November 11, 2020 ASU C*-Seminar Priyanga Ganesan Quantum Graphs Classical Graphs G = ( V , E , A G ) Vertex set: V Edge set: E V V Adjacency matrix: A G = [ a ij ], where a ij = 1 if ( i , j )

808 views • 24 slides
n -particle quantum statistics on graphs Jon Harrison 1 , J.P. Keating 2 , J.M. Robbins 2 and A.

n -particle quantum statistics on graphs Jon Harrison 1 , J.P. Keating 2 , J.M. Robbins 2 and A.

Quantum statistics Statistics on graphs 3 -connected graphs n -particle quantum statistics on graphs Jon Harrison 1 , J.P. Keating 2 , J.M. Robbins 2 and A. Sawicki 2 1 Baylor University, 2 University of Bristol QMath13 10/16 Jon Harrison

745 views • 46 slides
WELCOME! Agenda Recap Task Force Report Factor 3: Growth Next Steps Closing

WELCOME! Agenda Recap Task Force Report Factor 3: Growth Next Steps Closing

FACILITIES PLANNING COMMITTEE Marble Falls ISD Meeting 4 May 23, 2018 Love and Inspire WELCOME! Agenda Recap Task Force Report Factor 3: Growth Next Steps Closing Recap The Long Range Planning Process Who we want to be

703 views • 31 slides
Experimental quantum fast Carlo Di Franco hitting on hexagonal graphs 9th International

Experimental quantum fast Carlo Di Franco hitting on hexagonal graphs 9th International

Experimental quantum fast Carlo Di Franco hitting on hexagonal graphs 9th International Conference on Quantum Simulation and Quantum Walks CIRM, Marseille, France 21st January 2020 Outline A feasible experimental platform for implementing

627 views • 39 slides
Form factor approach to the correlation functions of quantum integrable models N. Kitanine IMB,

Form factor approach to the correlation functions of quantum integrable models N. Kitanine IMB,

Form factor approach to the correlation functions of quantum integrable models N. Kitanine IMB, Universit e de Bourgogne Quantum Integrable Systems and Geometry Olh ao, September 2012 In collaboration with : K. K. Kozlowski, J.M.

299 views • 27 slides
The effect of spin in the spectral statistics of Pauli operator Spin contribution quantum graphs

The effect of spin in the spectral statistics of Pauli operator Spin contribution quantum graphs

RMT Pauli operator Spin contribution to form factor Effect of spin in spectral statistics Jon Harrison RMT The effect of spin in the spectral statistics of Pauli operator Spin contribution quantum graphs to form factor J. M. Harrison

403 views • 21 slides
Quantum graphs where back-scattering is prohibited Brian Winn School of Mathematics

Quantum graphs where back-scattering is prohibited Brian Winn School of Mathematics

Quantum graphs where back-scattering is prohibited Brian Winn School of Mathematics Loughborough University 26th February 2008 Brian Winn Quantum graphs without back-scattering Credits Joint work with Jon Harrison and Uzy Smilansky

883 views • 63 slides
Quantum Mechanics for Graphs and CW-Complexes Michael Toriyama, Zhe Hu, Boyan Xu, Chengzheng Yu

Quantum Mechanics for Graphs and CW-Complexes Michael Toriyama, Zhe Hu, Boyan Xu, Chengzheng Yu

Quantum Mechanics for Graphs and CW-Complexes Michael Toriyama, Zhe Hu, Boyan Xu, Chengzheng Yu Sarah Loeb Ivan Contreras University of Illinois at Urbana-Champaign Illinois Geometry Lab Midterm Presentation October 25, 2016 Quantum

373 views • 9 slides
New Food and Land Use Economy F O L U C O L O M B I A The FOLU coalition is a global

New Food and Land Use Economy F O L U C O L O M B I A The FOLU coalition is a global

New Food and Land Use Economy F O L U C O L O M B I A The FOLU coalition is a global initiative Works to transform food and Government land use systems across the Businesses Scientific world by identifying solutions community that can

587 views • 57 slides
Crossed products of crossed modules of Hopf monoids in a braided setting Ramn Gonzlez

Crossed products of crossed modules of Hopf monoids in a braided setting Ramn Gonzlez

The setting Some definitions of crossed modules of Hopf monoids A new definition Crossed products of crossed modules of Hopf monoids Projections Crossed products of crossed modules of Hopf monoids in a braided setting Ramn Gonzlez

1.27k views • 110 slides
A Cache-Oblivious Heap Introduced by Arge et al. [1]. Based on distribution of elements

A Cache-Oblivious Heap Introduced by Arge et al. [1]. Based on distribution of elements

A Cache-Oblivious Heap Introduced by Arge et al. [1]. Based on distribution of elements References [1] L. Arge, M. A. Bender, E. D. Demaine, B. Holland-Minkley, and J. I. Munro, Cache-oblivious priority queue and graph algorithm

520 views • 15 slides
HELIX: A Case Study of a Formal Verification of High Performance Program Generation Vadim Zaliva

HELIX: A Case Study of a Formal Verification of High Performance Program Generation Vadim Zaliva

HELIX: A Case Study of a Formal Verification of High Performance Program Generation Vadim Zaliva Franz Franchetti Department of Electrical and Computer Engineering Carnegie Mellon University FHPC18 Vadim Zaliva, Franz Franchetti (CMU)

799 views • 61 slides
On the Subadditivity of Tail Value at Risk g An Investigation with Copulas g g CAS Spring

On the Subadditivity of Tail Value at Risk g An Investigation with Copulas g g CAS Spring

On the Subadditivity of Tail Value at Risk g An Investigation with Copulas g g CAS Spring meeting, May 2009 S. Desmedt and J.F . Walhin S. Desmedt and J.F . Walhin, CAS Spring meeting, New Orleans, May 2009 p. 1/45 Outline g

714 views • 48 slides
The future of operating systems on RISC-V Alex Bradbury asb@lowrisc.org @asbradbury 4th

The future of operating systems on RISC-V Alex Bradbury asb@lowrisc.org @asbradbury 4th

The future of operating systems on RISC-V Alex Bradbury asb@lowrisc.org @asbradbury 4th March 2019 Structure of this talk Introduction to RISC-V RISC-V status Selected RISC-V topics RISC-V and open hardware: the future

733 views • 28 slides
Mid-Term Progress Report Chile National Strategy on Climate Change and Vegetation Resources San

Mid-Term Progress Report Chile National Strategy on Climate Change and Vegetation Resources San

Mid-Term Progress Report Chile National Strategy on Climate Change and Vegetation Resources San Jos, Costa Rica, November 2015 Content I. Progress in the REDD+ Readiness Process. II. Progress achieved in the REDD+ readiness activities

834 views • 38 slides
Integer Invariants of Abelian Cayley Graphs Deelan Jalil James Madison University July 26, 2013

Integer Invariants of Abelian Cayley Graphs Deelan Jalil James Madison University July 26, 2013

Integer Invariants of Abelian Cayley Graphs Deelan Jalil James Madison University July 26, 2013 Deelan Jalil (2013) Integer Invariants of Abelian Cayley Graphs July 26, 2013 1 / 29 Cube example. . Deelan Jalil (2013) Integer Invariants of

1.11k views • 71 slides

More recommend