Partial trace over a product of matrices - one factor is in tensor product form

The name of the pictureThe name of the pictureThe name of the pictureClash Royale CLAN TAG#URR8PPP












6














$$Tr(rho^AB (sigma^A otimes I/d)) = Tr(rho^A sigma^A)$$



I came across the above, but I'm not sure how it's true. I figured they first partial traced out the B subsystem, and then trace A, but I don't see how you are allowed to partial trace out B from both the factors in the arguments. A proof or any intuition on this would be appreciated.



Edit 1:



The notation



$rho^AB$ is a state in Hilbert space $H_A otimes H_B$



$sigma^A$ is a state in Hilbert space $H_A$



$rho^A$ is $rho^AB$ with $B$ subsystem traced out.



$I/d$ is the maximally mixed state in Hilbert space $B$.



I saw this being used in Nielsen and Chuang, section 11.3.4, in the proof of subadditivity of entropy.



Edit 2:



So, I tried to write an answer based on DaftWullie's comment and Алексей Уваров's answer, but I am stuck again.



So, $$rho^AB = sum_mnop rho_mnop |morangle langle np|$$



Then $$rho^A = sum_mno rho_mnoo |mrangle langle n|$$



Let $$sigma^A = sum_ij sigma_ij |irangle langle j|$$



And $$I/d = sum_xy [I/d]_xy |xrangle langle y|$$



RHS



$$Tr(rho^A sigma^A)\
= Tr(sum_mno rho_mnoo |mrangle langle n|sum_ij sigma_ij |irangle langle j|)\
= Tr(sum_mnoj rho_mnoo sigma_nj | m rangle langle j|)\
= sum_mno rho_mnoo sigma_nm$$



LHS



$$Tr(rho^AB (sigma^A otimes I/d)\
= Tr(sum_mnop rho_mnop |morangle langle np| sum_ijxy sigma_ij [I/d]_xy |ixrangle langle jy|)\
= Tr(sum_mnoxjyrho_mnox sigma_nj [I/d]_xy | mo rangle langle jy |)\
= sum_mnyx rho_nm[I/d]_xy\
= (1/d)sum_mny rho_mnyy sigma_nm$$



Which is the same as the RHS, but there's an extra $1/d$ factor?



Also, am I thinking about this the wrong way? Is there a simpler way to look at this?










share|improve this question



















  • 2




    How do you define $rho_a$ in terms of $rho_ab$?
    – DaftWullie
    Dec 25 '18 at 8:34










  • Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
    – Blue
    Dec 25 '18 at 9:55










  • @Blue, edited title, check now?
    – Mahathi Vempati
    Dec 25 '18 at 11:05











  • @DaftWullie, Edited the question
    – Mahathi Vempati
    Dec 25 '18 at 11:06






  • 1




    @MahathiVempati Thanks. Looks good!
    – Blue
    Dec 25 '18 at 13:09















6














$$Tr(rho^AB (sigma^A otimes I/d)) = Tr(rho^A sigma^A)$$



I came across the above, but I'm not sure how it's true. I figured they first partial traced out the B subsystem, and then trace A, but I don't see how you are allowed to partial trace out B from both the factors in the arguments. A proof or any intuition on this would be appreciated.



Edit 1:



The notation



$rho^AB$ is a state in Hilbert space $H_A otimes H_B$



$sigma^A$ is a state in Hilbert space $H_A$



$rho^A$ is $rho^AB$ with $B$ subsystem traced out.



$I/d$ is the maximally mixed state in Hilbert space $B$.



I saw this being used in Nielsen and Chuang, section 11.3.4, in the proof of subadditivity of entropy.



Edit 2:



So, I tried to write an answer based on DaftWullie's comment and Алексей Уваров's answer, but I am stuck again.



So, $$rho^AB = sum_mnop rho_mnop |morangle langle np|$$



Then $$rho^A = sum_mno rho_mnoo |mrangle langle n|$$



Let $$sigma^A = sum_ij sigma_ij |irangle langle j|$$



And $$I/d = sum_xy [I/d]_xy |xrangle langle y|$$



RHS



$$Tr(rho^A sigma^A)\
= Tr(sum_mno rho_mnoo |mrangle langle n|sum_ij sigma_ij |irangle langle j|)\
= Tr(sum_mnoj rho_mnoo sigma_nj | m rangle langle j|)\
= sum_mno rho_mnoo sigma_nm$$



LHS



$$Tr(rho^AB (sigma^A otimes I/d)\
= Tr(sum_mnop rho_mnop |morangle langle np| sum_ijxy sigma_ij [I/d]_xy |ixrangle langle jy|)\
= Tr(sum_mnoxjyrho_mnox sigma_nj [I/d]_xy | mo rangle langle jy |)\
= sum_mnyx rho_nm[I/d]_xy\
= (1/d)sum_mny rho_mnyy sigma_nm$$



Which is the same as the RHS, but there's an extra $1/d$ factor?



Also, am I thinking about this the wrong way? Is there a simpler way to look at this?










share|improve this question



















  • 2




    How do you define $rho_a$ in terms of $rho_ab$?
    – DaftWullie
    Dec 25 '18 at 8:34










  • Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
    – Blue
    Dec 25 '18 at 9:55










  • @Blue, edited title, check now?
    – Mahathi Vempati
    Dec 25 '18 at 11:05











  • @DaftWullie, Edited the question
    – Mahathi Vempati
    Dec 25 '18 at 11:06






  • 1




    @MahathiVempati Thanks. Looks good!
    – Blue
    Dec 25 '18 at 13:09













6












6








6


1





$$Tr(rho^AB (sigma^A otimes I/d)) = Tr(rho^A sigma^A)$$



I came across the above, but I'm not sure how it's true. I figured they first partial traced out the B subsystem, and then trace A, but I don't see how you are allowed to partial trace out B from both the factors in the arguments. A proof or any intuition on this would be appreciated.



Edit 1:



The notation



$rho^AB$ is a state in Hilbert space $H_A otimes H_B$



$sigma^A$ is a state in Hilbert space $H_A$



$rho^A$ is $rho^AB$ with $B$ subsystem traced out.



$I/d$ is the maximally mixed state in Hilbert space $B$.



I saw this being used in Nielsen and Chuang, section 11.3.4, in the proof of subadditivity of entropy.



Edit 2:



So, I tried to write an answer based on DaftWullie's comment and Алексей Уваров's answer, but I am stuck again.



So, $$rho^AB = sum_mnop rho_mnop |morangle langle np|$$



Then $$rho^A = sum_mno rho_mnoo |mrangle langle n|$$



Let $$sigma^A = sum_ij sigma_ij |irangle langle j|$$



And $$I/d = sum_xy [I/d]_xy |xrangle langle y|$$



RHS



$$Tr(rho^A sigma^A)\
= Tr(sum_mno rho_mnoo |mrangle langle n|sum_ij sigma_ij |irangle langle j|)\
= Tr(sum_mnoj rho_mnoo sigma_nj | m rangle langle j|)\
= sum_mno rho_mnoo sigma_nm$$



LHS



$$Tr(rho^AB (sigma^A otimes I/d)\
= Tr(sum_mnop rho_mnop |morangle langle np| sum_ijxy sigma_ij [I/d]_xy |ixrangle langle jy|)\
= Tr(sum_mnoxjyrho_mnox sigma_nj [I/d]_xy | mo rangle langle jy |)\
= sum_mnyx rho_nm[I/d]_xy\
= (1/d)sum_mny rho_mnyy sigma_nm$$



Which is the same as the RHS, but there's an extra $1/d$ factor?



Also, am I thinking about this the wrong way? Is there a simpler way to look at this?










share|improve this question















$$Tr(rho^AB (sigma^A otimes I/d)) = Tr(rho^A sigma^A)$$



I came across the above, but I'm not sure how it's true. I figured they first partial traced out the B subsystem, and then trace A, but I don't see how you are allowed to partial trace out B from both the factors in the arguments. A proof or any intuition on this would be appreciated.



Edit 1:



The notation



$rho^AB$ is a state in Hilbert space $H_A otimes H_B$



$sigma^A$ is a state in Hilbert space $H_A$



$rho^A$ is $rho^AB$ with $B$ subsystem traced out.



$I/d$ is the maximally mixed state in Hilbert space $B$.



I saw this being used in Nielsen and Chuang, section 11.3.4, in the proof of subadditivity of entropy.



Edit 2:



So, I tried to write an answer based on DaftWullie's comment and Алексей Уваров's answer, but I am stuck again.



So, $$rho^AB = sum_mnop rho_mnop |morangle langle np|$$



Then $$rho^A = sum_mno rho_mnoo |mrangle langle n|$$



Let $$sigma^A = sum_ij sigma_ij |irangle langle j|$$



And $$I/d = sum_xy [I/d]_xy |xrangle langle y|$$



RHS



$$Tr(rho^A sigma^A)\
= Tr(sum_mno rho_mnoo |mrangle langle n|sum_ij sigma_ij |irangle langle j|)\
= Tr(sum_mnoj rho_mnoo sigma_nj | m rangle langle j|)\
= sum_mno rho_mnoo sigma_nm$$



LHS



$$Tr(rho^AB (sigma^A otimes I/d)\
= Tr(sum_mnop rho_mnop |morangle langle np| sum_ijxy sigma_ij [I/d]_xy |ixrangle langle jy|)\
= Tr(sum_mnoxjyrho_mnox sigma_nj [I/d]_xy | mo rangle langle jy |)\
= sum_mnyx rho_nm[I/d]_xy\
= (1/d)sum_mny rho_mnyy sigma_nm$$



Which is the same as the RHS, but there's an extra $1/d$ factor?



Also, am I thinking about this the wrong way? Is there a simpler way to look at this?







quantum-information density-matrix tensor-product nielsen-and-chuang partial-trace






share|improve this question















share|improve this question













share|improve this question




share|improve this question








edited Dec 31 '18 at 23:43









Blue

5,67221354




5,67221354










asked Dec 25 '18 at 8:15









Mahathi Vempati

3848




3848







  • 2




    How do you define $rho_a$ in terms of $rho_ab$?
    – DaftWullie
    Dec 25 '18 at 8:34










  • Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
    – Blue
    Dec 25 '18 at 9:55










  • @Blue, edited title, check now?
    – Mahathi Vempati
    Dec 25 '18 at 11:05











  • @DaftWullie, Edited the question
    – Mahathi Vempati
    Dec 25 '18 at 11:06






  • 1




    @MahathiVempati Thanks. Looks good!
    – Blue
    Dec 25 '18 at 13:09












  • 2




    How do you define $rho_a$ in terms of $rho_ab$?
    – DaftWullie
    Dec 25 '18 at 8:34










  • Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
    – Blue
    Dec 25 '18 at 9:55










  • @Blue, edited title, check now?
    – Mahathi Vempati
    Dec 25 '18 at 11:05











  • @DaftWullie, Edited the question
    – Mahathi Vempati
    Dec 25 '18 at 11:06






  • 1




    @MahathiVempati Thanks. Looks good!
    – Blue
    Dec 25 '18 at 13:09







2




2




How do you define $rho_a$ in terms of $rho_ab$?
– DaftWullie
Dec 25 '18 at 8:34




How do you define $rho_a$ in terms of $rho_ab$?
– DaftWullie
Dec 25 '18 at 8:34












Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
– Blue
Dec 25 '18 at 9:55




Hi Maharthi. It would be better if you can edit to make the title more descriptive. Questions with math-only titles are not easily searchable.
– Blue
Dec 25 '18 at 9:55












@Blue, edited title, check now?
– Mahathi Vempati
Dec 25 '18 at 11:05





@Blue, edited title, check now?
– Mahathi Vempati
Dec 25 '18 at 11:05













@DaftWullie, Edited the question
– Mahathi Vempati
Dec 25 '18 at 11:06




@DaftWullie, Edited the question
– Mahathi Vempati
Dec 25 '18 at 11:06




1




1




@MahathiVempati Thanks. Looks good!
– Blue
Dec 25 '18 at 13:09




@MahathiVempati Thanks. Looks good!
– Blue
Dec 25 '18 at 13:09










2 Answers
2






active

oldest

votes


















4














The equation at the top of the question is not correct: there is a missing factor of $1/d$ on the right-hand side. Let's eliminate this factor from the left-hand side to make it simpler, so that the equation we want is this:
$$
textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr).
$$



To see why this is true, it helps to start with an easy special case, which is that $rho^AB$ is a product state:
$$
rho^AB = rho^A otimes rho^B.
$$

In this case we have
$$
textTrbigl(bigl(rho^A otimesrho^Bbigr) bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr)textTrbigl(rho^Bbigr) = textTrbigl(rho^A sigma^Abigr),
$$

using just elementary properties of tensor products and their traces.



Now, given that the equation is true in the special case, it has to be true in general because the expressions
$$
textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr);;textand;;textTrbigl(rho^A sigma^Abigr)
$$

depend linearly on $rho^AB$, and the set of all product states $rho^Aotimesrho^B$ spans the vector space of all operators acting on $H_Aotimes H_B$.



Alternatively, we have
$$
textTr((Xotimes Y)(Zotimes I)) = textTr(XZ)textTr(Y) = textTrbigl(textTr_B(Xotimes Y), Zbigr)
$$

for all operators $X$ and $Z$ acting on $H_A$ and all $Y$ acting on $H_B$, irrespective of their traces, and therefore
$$
textTr(W(Zotimes I)) = textTrbigl(textTr_B(W) Zbigr)
$$

for all operators $W$ acting on $H_Aotimes H_B$ by linearity.






share|improve this answer




























    6














    Here the important fact is that the maximally mixed state is in fact an identity matrix.



    Let me rewrite the expression on the left in index notation (the summation sign is omitted according to the Einstein convention):



    $$
    Tr(rho^AB (sigma^A otimes I/d)) = [rho^AB]_ijkl [sigma^A]_ji [I/d]_lk
    $$



    But $[I/d]_lk = frac1d delta_lk$, therefore $[rho^AB]_ijkl [sigma^A]_ji [I/d]_lk = frac1d [rho^AB]_ijkk [sigma^A]_ji$, which is exactly what happens if you first trace out the subsystem $B$ (UPD: up to the prefactor of $1/d$ apparently).



    The physical intuition would be as follows. This expression is basically an expected value of a Hermitian operator $frac1d sigma^A otimes I$ over a state $rho$. This operator only acts nontrivially on the first subsystem, thus we can safely trace out the rest.



    EDIT: Also, this contraction problem can be understood better if you use tensor network notation. Learning it requires some time, but if you do, I suggest starting here and here.






    share|improve this answer






















    • Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
      – Mahathi Vempati
      Dec 25 '18 at 12:46










    • @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
      – Алексей Уваров
      Dec 25 '18 at 13:27











    • Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
      – Алексей Уваров
      Dec 25 '18 at 13:33






    • 1




      Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
      – Алексей Уваров
      Dec 26 '18 at 8:18







    • 1




      @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
      – Алексей Уваров
      Dec 28 '18 at 8:53










    Your Answer





    StackExchange.ifUsing("editor", function ()
    return StackExchange.using("mathjaxEditing", function ()
    StackExchange.MarkdownEditor.creationCallbacks.add(function (editor, postfix)
    StackExchange.mathjaxEditing.prepareWmdForMathJax(editor, postfix, [["$", "$"], ["\\(","\\)"]]);
    );
    );
    , "mathjax-editing");

    StackExchange.ready(function()
    var channelOptions =
    tags: "".split(" "),
    id: "694"
    ;
    initTagRenderer("".split(" "), "".split(" "), channelOptions);

    StackExchange.using("externalEditor", function()
    // Have to fire editor after snippets, if snippets enabled
    if (StackExchange.settings.snippets.snippetsEnabled)
    StackExchange.using("snippets", function()
    createEditor();
    );

    else
    createEditor();

    );

    function createEditor()
    StackExchange.prepareEditor(
    heartbeatType: 'answer',
    autoActivateHeartbeat: false,
    convertImagesToLinks: false,
    noModals: true,
    showLowRepImageUploadWarning: true,
    reputationToPostImages: null,
    bindNavPrevention: true,
    postfix: "",
    imageUploader:
    brandingHtml: "Powered by u003ca class="icon-imgur-white" href="https://imgur.com/"u003eu003c/au003e",
    contentPolicyHtml: "User contributions licensed under u003ca href="https://creativecommons.org/licenses/by-sa/3.0/"u003ecc by-sa 3.0 with attribution requiredu003c/au003e u003ca href="https://stackoverflow.com/legal/content-policy"u003e(content policy)u003c/au003e",
    allowUrls: true
    ,
    noCode: true, onDemand: true,
    discardSelector: ".discard-answer"
    ,immediatelyShowMarkdownHelp:true
    );



    );













    draft saved

    draft discarded


















    StackExchange.ready(
    function ()
    StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fquantumcomputing.stackexchange.com%2fquestions%2f5045%2fpartial-trace-over-a-product-of-matrices-one-factor-is-in-tensor-product-form%23new-answer', 'question_page');

    );

    Post as a guest















    Required, but never shown

























    2 Answers
    2






    active

    oldest

    votes








    2 Answers
    2






    active

    oldest

    votes









    active

    oldest

    votes






    active

    oldest

    votes









    4














    The equation at the top of the question is not correct: there is a missing factor of $1/d$ on the right-hand side. Let's eliminate this factor from the left-hand side to make it simpler, so that the equation we want is this:
    $$
    textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr).
    $$



    To see why this is true, it helps to start with an easy special case, which is that $rho^AB$ is a product state:
    $$
    rho^AB = rho^A otimes rho^B.
    $$

    In this case we have
    $$
    textTrbigl(bigl(rho^A otimesrho^Bbigr) bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr)textTrbigl(rho^Bbigr) = textTrbigl(rho^A sigma^Abigr),
    $$

    using just elementary properties of tensor products and their traces.



    Now, given that the equation is true in the special case, it has to be true in general because the expressions
    $$
    textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr);;textand;;textTrbigl(rho^A sigma^Abigr)
    $$

    depend linearly on $rho^AB$, and the set of all product states $rho^Aotimesrho^B$ spans the vector space of all operators acting on $H_Aotimes H_B$.



    Alternatively, we have
    $$
    textTr((Xotimes Y)(Zotimes I)) = textTr(XZ)textTr(Y) = textTrbigl(textTr_B(Xotimes Y), Zbigr)
    $$

    for all operators $X$ and $Z$ acting on $H_A$ and all $Y$ acting on $H_B$, irrespective of their traces, and therefore
    $$
    textTr(W(Zotimes I)) = textTrbigl(textTr_B(W) Zbigr)
    $$

    for all operators $W$ acting on $H_Aotimes H_B$ by linearity.






    share|improve this answer

























      4














      The equation at the top of the question is not correct: there is a missing factor of $1/d$ on the right-hand side. Let's eliminate this factor from the left-hand side to make it simpler, so that the equation we want is this:
      $$
      textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr).
      $$



      To see why this is true, it helps to start with an easy special case, which is that $rho^AB$ is a product state:
      $$
      rho^AB = rho^A otimes rho^B.
      $$

      In this case we have
      $$
      textTrbigl(bigl(rho^A otimesrho^Bbigr) bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr)textTrbigl(rho^Bbigr) = textTrbigl(rho^A sigma^Abigr),
      $$

      using just elementary properties of tensor products and their traces.



      Now, given that the equation is true in the special case, it has to be true in general because the expressions
      $$
      textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr);;textand;;textTrbigl(rho^A sigma^Abigr)
      $$

      depend linearly on $rho^AB$, and the set of all product states $rho^Aotimesrho^B$ spans the vector space of all operators acting on $H_Aotimes H_B$.



      Alternatively, we have
      $$
      textTr((Xotimes Y)(Zotimes I)) = textTr(XZ)textTr(Y) = textTrbigl(textTr_B(Xotimes Y), Zbigr)
      $$

      for all operators $X$ and $Z$ acting on $H_A$ and all $Y$ acting on $H_B$, irrespective of their traces, and therefore
      $$
      textTr(W(Zotimes I)) = textTrbigl(textTr_B(W) Zbigr)
      $$

      for all operators $W$ acting on $H_Aotimes H_B$ by linearity.






      share|improve this answer























        4












        4








        4






        The equation at the top of the question is not correct: there is a missing factor of $1/d$ on the right-hand side. Let's eliminate this factor from the left-hand side to make it simpler, so that the equation we want is this:
        $$
        textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr).
        $$



        To see why this is true, it helps to start with an easy special case, which is that $rho^AB$ is a product state:
        $$
        rho^AB = rho^A otimes rho^B.
        $$

        In this case we have
        $$
        textTrbigl(bigl(rho^A otimesrho^Bbigr) bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr)textTrbigl(rho^Bbigr) = textTrbigl(rho^A sigma^Abigr),
        $$

        using just elementary properties of tensor products and their traces.



        Now, given that the equation is true in the special case, it has to be true in general because the expressions
        $$
        textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr);;textand;;textTrbigl(rho^A sigma^Abigr)
        $$

        depend linearly on $rho^AB$, and the set of all product states $rho^Aotimesrho^B$ spans the vector space of all operators acting on $H_Aotimes H_B$.



        Alternatively, we have
        $$
        textTr((Xotimes Y)(Zotimes I)) = textTr(XZ)textTr(Y) = textTrbigl(textTr_B(Xotimes Y), Zbigr)
        $$

        for all operators $X$ and $Z$ acting on $H_A$ and all $Y$ acting on $H_B$, irrespective of their traces, and therefore
        $$
        textTr(W(Zotimes I)) = textTrbigl(textTr_B(W) Zbigr)
        $$

        for all operators $W$ acting on $H_Aotimes H_B$ by linearity.






        share|improve this answer












        The equation at the top of the question is not correct: there is a missing factor of $1/d$ on the right-hand side. Let's eliminate this factor from the left-hand side to make it simpler, so that the equation we want is this:
        $$
        textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr).
        $$



        To see why this is true, it helps to start with an easy special case, which is that $rho^AB$ is a product state:
        $$
        rho^AB = rho^A otimes rho^B.
        $$

        In this case we have
        $$
        textTrbigl(bigl(rho^A otimesrho^Bbigr) bigl(sigma^A otimes Ibigr)bigr) = textTrbigl(rho^A sigma^Abigr)textTrbigl(rho^Bbigr) = textTrbigl(rho^A sigma^Abigr),
        $$

        using just elementary properties of tensor products and their traces.



        Now, given that the equation is true in the special case, it has to be true in general because the expressions
        $$
        textTrbigl(rho^AB bigl(sigma^A otimes Ibigr)bigr);;textand;;textTrbigl(rho^A sigma^Abigr)
        $$

        depend linearly on $rho^AB$, and the set of all product states $rho^Aotimesrho^B$ spans the vector space of all operators acting on $H_Aotimes H_B$.



        Alternatively, we have
        $$
        textTr((Xotimes Y)(Zotimes I)) = textTr(XZ)textTr(Y) = textTrbigl(textTr_B(Xotimes Y), Zbigr)
        $$

        for all operators $X$ and $Z$ acting on $H_A$ and all $Y$ acting on $H_B$, irrespective of their traces, and therefore
        $$
        textTr(W(Zotimes I)) = textTrbigl(textTr_B(W) Zbigr)
        $$

        for all operators $W$ acting on $H_Aotimes H_B$ by linearity.







        share|improve this answer












        share|improve this answer



        share|improve this answer










        answered Dec 28 '18 at 13:36









        John Watrous

        1,41818




        1,41818























            6














            Here the important fact is that the maximally mixed state is in fact an identity matrix.



            Let me rewrite the expression on the left in index notation (the summation sign is omitted according to the Einstein convention):



            $$
            Tr(rho^AB (sigma^A otimes I/d)) = [rho^AB]_ijkl [sigma^A]_ji [I/d]_lk
            $$



            But $[I/d]_lk = frac1d delta_lk$, therefore $[rho^AB]_ijkl [sigma^A]_ji [I/d]_lk = frac1d [rho^AB]_ijkk [sigma^A]_ji$, which is exactly what happens if you first trace out the subsystem $B$ (UPD: up to the prefactor of $1/d$ apparently).



            The physical intuition would be as follows. This expression is basically an expected value of a Hermitian operator $frac1d sigma^A otimes I$ over a state $rho$. This operator only acts nontrivially on the first subsystem, thus we can safely trace out the rest.



            EDIT: Also, this contraction problem can be understood better if you use tensor network notation. Learning it requires some time, but if you do, I suggest starting here and here.






            share|improve this answer






















            • Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
              – Mahathi Vempati
              Dec 25 '18 at 12:46










            • @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
              – Алексей Уваров
              Dec 25 '18 at 13:27











            • Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
              – Алексей Уваров
              Dec 25 '18 at 13:33






            • 1




              Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
              – Алексей Уваров
              Dec 26 '18 at 8:18







            • 1




              @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
              – Алексей Уваров
              Dec 28 '18 at 8:53















            6














            Here the important fact is that the maximally mixed state is in fact an identity matrix.



            Let me rewrite the expression on the left in index notation (the summation sign is omitted according to the Einstein convention):



            $$
            Tr(rho^AB (sigma^A otimes I/d)) = [rho^AB]_ijkl [sigma^A]_ji [I/d]_lk
            $$



            But $[I/d]_lk = frac1d delta_lk$, therefore $[rho^AB]_ijkl [sigma^A]_ji [I/d]_lk = frac1d [rho^AB]_ijkk [sigma^A]_ji$, which is exactly what happens if you first trace out the subsystem $B$ (UPD: up to the prefactor of $1/d$ apparently).



            The physical intuition would be as follows. This expression is basically an expected value of a Hermitian operator $frac1d sigma^A otimes I$ over a state $rho$. This operator only acts nontrivially on the first subsystem, thus we can safely trace out the rest.



            EDIT: Also, this contraction problem can be understood better if you use tensor network notation. Learning it requires some time, but if you do, I suggest starting here and here.






            share|improve this answer






















            • Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
              – Mahathi Vempati
              Dec 25 '18 at 12:46










            • @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
              – Алексей Уваров
              Dec 25 '18 at 13:27











            • Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
              – Алексей Уваров
              Dec 25 '18 at 13:33






            • 1




              Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
              – Алексей Уваров
              Dec 26 '18 at 8:18







            • 1




              @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
              – Алексей Уваров
              Dec 28 '18 at 8:53













            6












            6








            6






            Here the important fact is that the maximally mixed state is in fact an identity matrix.



            Let me rewrite the expression on the left in index notation (the summation sign is omitted according to the Einstein convention):



            $$
            Tr(rho^AB (sigma^A otimes I/d)) = [rho^AB]_ijkl [sigma^A]_ji [I/d]_lk
            $$



            But $[I/d]_lk = frac1d delta_lk$, therefore $[rho^AB]_ijkl [sigma^A]_ji [I/d]_lk = frac1d [rho^AB]_ijkk [sigma^A]_ji$, which is exactly what happens if you first trace out the subsystem $B$ (UPD: up to the prefactor of $1/d$ apparently).



            The physical intuition would be as follows. This expression is basically an expected value of a Hermitian operator $frac1d sigma^A otimes I$ over a state $rho$. This operator only acts nontrivially on the first subsystem, thus we can safely trace out the rest.



            EDIT: Also, this contraction problem can be understood better if you use tensor network notation. Learning it requires some time, but if you do, I suggest starting here and here.






            share|improve this answer














            Here the important fact is that the maximally mixed state is in fact an identity matrix.



            Let me rewrite the expression on the left in index notation (the summation sign is omitted according to the Einstein convention):



            $$
            Tr(rho^AB (sigma^A otimes I/d)) = [rho^AB]_ijkl [sigma^A]_ji [I/d]_lk
            $$



            But $[I/d]_lk = frac1d delta_lk$, therefore $[rho^AB]_ijkl [sigma^A]_ji [I/d]_lk = frac1d [rho^AB]_ijkk [sigma^A]_ji$, which is exactly what happens if you first trace out the subsystem $B$ (UPD: up to the prefactor of $1/d$ apparently).



            The physical intuition would be as follows. This expression is basically an expected value of a Hermitian operator $frac1d sigma^A otimes I$ over a state $rho$. This operator only acts nontrivially on the first subsystem, thus we can safely trace out the rest.



            EDIT: Also, this contraction problem can be understood better if you use tensor network notation. Learning it requires some time, but if you do, I suggest starting here and here.







            share|improve this answer














            share|improve this answer



            share|improve this answer








            edited Dec 28 '18 at 9:06

























            answered Dec 25 '18 at 12:42









            Алексей Уваров

            1216




            1216











            • Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
              – Mahathi Vempati
              Dec 25 '18 at 12:46










            • @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
              – Алексей Уваров
              Dec 25 '18 at 13:27











            • Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
              – Алексей Уваров
              Dec 25 '18 at 13:33






            • 1




              Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
              – Алексей Уваров
              Dec 26 '18 at 8:18







            • 1




              @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
              – Алексей Уваров
              Dec 28 '18 at 8:53
















            • Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
              – Mahathi Vempati
              Dec 25 '18 at 12:46










            • @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
              – Алексей Уваров
              Dec 25 '18 at 13:27











            • Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
              – Алексей Уваров
              Dec 25 '18 at 13:33






            • 1




              Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
              – Алексей Уваров
              Dec 26 '18 at 8:18







            • 1




              @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
              – Алексей Уваров
              Dec 28 '18 at 8:53















            Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
            – Mahathi Vempati
            Dec 25 '18 at 12:46




            Hey, thank you very much for answering! But I don't understand the Einstein notation very well. Can you explain what the rhs means?
            – Mahathi Vempati
            Dec 25 '18 at 12:46












            @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
            – Алексей Уваров
            Dec 25 '18 at 13:27





            @MahathiVempati By $[M]_ij$ I mean the element of a matrix M with indices i, j. The square brackets are unnecessary in general, but I used them to separate the subsystem indices like A, B from the tensor indices. Now, whenever there are two coinciding indices, there is a summation over this index, e.g. an ordinary matrix-vector product would look like $A_ij x_j$
            – Алексей Уваров
            Dec 25 '18 at 13:27













            Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
            – Алексей Уваров
            Dec 25 '18 at 13:33




            Now, intuitively what's going on in the rhs is that the density matrix $rho^AB$ is contracted with the density matrix $sigma^A otimes I/d$, and then all of this is traced. If they both were presented as regular matrices, it would just be a product, but in tensor form it might be somewhat more confusing.
            – Алексей Уваров
            Dec 25 '18 at 13:33




            1




            1




            Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
            – Алексей Уваров
            Dec 26 '18 at 8:18





            Hmm, when I was answering the question, I thought about ij as the row and column for the first subsystem and kl as that for the second. I think it can be written as follows: $rho = sum_i,j,k,l rho_ijkl (| irangle langle j | otimes | krangle langle l |)$.
            – Алексей Уваров
            Dec 26 '18 at 8:18





            1




            1




            @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
            – Алексей Уваров
            Dec 28 '18 at 8:53




            @MahathiVempati Hmm, looks like you're right about the $1/d$ prefactor. However, in my copy of Nielsen and Chuang I can't find this exact expression. They seem to prove subadditivity somewhat differently
            – Алексей Уваров
            Dec 28 '18 at 8:53

















            draft saved

            draft discarded
















































            Thanks for contributing an answer to Quantum Computing Stack Exchange!


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            Use MathJax to format equations. MathJax reference.


            To learn more, see our tips on writing great answers.





            Some of your past answers have not been well-received, and you're in danger of being blocked from answering.


            Please pay close attention to the following guidance:


            • Please be sure to answer the question. Provide details and share your research!

            But avoid


            • Asking for help, clarification, or responding to other answers.

            • Making statements based on opinion; back them up with references or personal experience.

            To learn more, see our tips on writing great answers.




            draft saved


            draft discarded














            StackExchange.ready(
            function ()
            StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fquantumcomputing.stackexchange.com%2fquestions%2f5045%2fpartial-trace-over-a-product-of-matrices-one-factor-is-in-tensor-product-form%23new-answer', 'question_page');

            );

            Post as a guest















            Required, but never shown





















































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown

































            Required, but never shown














            Required, but never shown












            Required, but never shown







            Required, but never shown






            Popular posts from this blog

            How to check contact read email or not when send email to Individual?

            Displaying single band from multi-band raster using QGIS

            How many registers does an x86_64 CPU actually have?