Schmidt decomposition - example
$begingroup$
I'm trying to compute the Schmidt decomposition of $left| psi right> = (left| 00 right> + left| 01 right> + left| 10 right>)/sqrt{3}$. This should be possible by first computing the reduced density matrices $rho_A=Tr_Bleft| psi right> left< psi right|=sum_i p_i left| a_i right> left< a_i right|$ and $rho_B=Tr_Aleft| psi right> left< psi right|=sum_i p_i left| b_i right> left< b_i right|$, and then by identifying $left| psi right>=sum_i sqrt{p_i}left| a_i right> otimes left| b_i right>$.
However, somewhere in the computation I'm doing a mistake: I find that
$$rho_A=rho_B=left(2left| 0 right> left< 0 right| + left| 0 right> left< 1 right| + left| 1 right> left< 0 right|+left| 1 right> left< 1 right|right)/3,$$
which has eigenvalues $lambda_1=0.87$ and $lambda_2=0.13$ with corresponding eigenvectors
$$ left| a_1 right>=0.85 left| 0 right> + 0.53 left| 1 right>,quad left| a_2 right>=-0.53 left| 0 right> + 0.85 left| 1 right>. $$
This means that $left| psi right>$ should be given by
$$ left| psi right> = sqrt{lambda_1} left| a_1 right> otimes
left| a_1 right> + sqrt{lambda_2} left| a_2 right> otimes
left| a_2 right>.$$
This, however, is not true. If you check for example the numerical value in front of $left| 00 right>$, you find that it is not equal to $1/sqrt{3}$.
I would appreciate if someone could help me to see where I made the mistake.
Thanks in advance.
quantum-mechanics quantum-information
New contributor
$endgroup$
add a comment |
$begingroup$
I'm trying to compute the Schmidt decomposition of $left| psi right> = (left| 00 right> + left| 01 right> + left| 10 right>)/sqrt{3}$. This should be possible by first computing the reduced density matrices $rho_A=Tr_Bleft| psi right> left< psi right|=sum_i p_i left| a_i right> left< a_i right|$ and $rho_B=Tr_Aleft| psi right> left< psi right|=sum_i p_i left| b_i right> left< b_i right|$, and then by identifying $left| psi right>=sum_i sqrt{p_i}left| a_i right> otimes left| b_i right>$.
However, somewhere in the computation I'm doing a mistake: I find that
$$rho_A=rho_B=left(2left| 0 right> left< 0 right| + left| 0 right> left< 1 right| + left| 1 right> left< 0 right|+left| 1 right> left< 1 right|right)/3,$$
which has eigenvalues $lambda_1=0.87$ and $lambda_2=0.13$ with corresponding eigenvectors
$$ left| a_1 right>=0.85 left| 0 right> + 0.53 left| 1 right>,quad left| a_2 right>=-0.53 left| 0 right> + 0.85 left| 1 right>. $$
This means that $left| psi right>$ should be given by
$$ left| psi right> = sqrt{lambda_1} left| a_1 right> otimes
left| a_1 right> + sqrt{lambda_2} left| a_2 right> otimes
left| a_2 right>.$$
This, however, is not true. If you check for example the numerical value in front of $left| 00 right>$, you find that it is not equal to $1/sqrt{3}$.
I would appreciate if someone could help me to see where I made the mistake.
Thanks in advance.
quantum-mechanics quantum-information
New contributor
$endgroup$
add a comment |
$begingroup$
I'm trying to compute the Schmidt decomposition of $left| psi right> = (left| 00 right> + left| 01 right> + left| 10 right>)/sqrt{3}$. This should be possible by first computing the reduced density matrices $rho_A=Tr_Bleft| psi right> left< psi right|=sum_i p_i left| a_i right> left< a_i right|$ and $rho_B=Tr_Aleft| psi right> left< psi right|=sum_i p_i left| b_i right> left< b_i right|$, and then by identifying $left| psi right>=sum_i sqrt{p_i}left| a_i right> otimes left| b_i right>$.
However, somewhere in the computation I'm doing a mistake: I find that
$$rho_A=rho_B=left(2left| 0 right> left< 0 right| + left| 0 right> left< 1 right| + left| 1 right> left< 0 right|+left| 1 right> left< 1 right|right)/3,$$
which has eigenvalues $lambda_1=0.87$ and $lambda_2=0.13$ with corresponding eigenvectors
$$ left| a_1 right>=0.85 left| 0 right> + 0.53 left| 1 right>,quad left| a_2 right>=-0.53 left| 0 right> + 0.85 left| 1 right>. $$
This means that $left| psi right>$ should be given by
$$ left| psi right> = sqrt{lambda_1} left| a_1 right> otimes
left| a_1 right> + sqrt{lambda_2} left| a_2 right> otimes
left| a_2 right>.$$
This, however, is not true. If you check for example the numerical value in front of $left| 00 right>$, you find that it is not equal to $1/sqrt{3}$.
I would appreciate if someone could help me to see where I made the mistake.
Thanks in advance.
quantum-mechanics quantum-information
New contributor
$endgroup$
I'm trying to compute the Schmidt decomposition of $left| psi right> = (left| 00 right> + left| 01 right> + left| 10 right>)/sqrt{3}$. This should be possible by first computing the reduced density matrices $rho_A=Tr_Bleft| psi right> left< psi right|=sum_i p_i left| a_i right> left< a_i right|$ and $rho_B=Tr_Aleft| psi right> left< psi right|=sum_i p_i left| b_i right> left< b_i right|$, and then by identifying $left| psi right>=sum_i sqrt{p_i}left| a_i right> otimes left| b_i right>$.
However, somewhere in the computation I'm doing a mistake: I find that
$$rho_A=rho_B=left(2left| 0 right> left< 0 right| + left| 0 right> left< 1 right| + left| 1 right> left< 0 right|+left| 1 right> left< 1 right|right)/3,$$
which has eigenvalues $lambda_1=0.87$ and $lambda_2=0.13$ with corresponding eigenvectors
$$ left| a_1 right>=0.85 left| 0 right> + 0.53 left| 1 right>,quad left| a_2 right>=-0.53 left| 0 right> + 0.85 left| 1 right>. $$
This means that $left| psi right>$ should be given by
$$ left| psi right> = sqrt{lambda_1} left| a_1 right> otimes
left| a_1 right> + sqrt{lambda_2} left| a_2 right> otimes
left| a_2 right>.$$
This, however, is not true. If you check for example the numerical value in front of $left| 00 right>$, you find that it is not equal to $1/sqrt{3}$.
I would appreciate if someone could help me to see where I made the mistake.
Thanks in advance.
quantum-mechanics quantum-information
quantum-mechanics quantum-information
New contributor
New contributor
New contributor
asked yesterday
ScottScott
161
161
New contributor
New contributor
add a comment |
add a comment |
2 Answers
2
active
oldest
votes
$begingroup$
The phase of your eigenvectors is not correct (or rather, it is not determined, so you need to make a judicious choice). If you put a minus sign in front of the 2nd term, it works out.
To elaborate on that: If you want to find the Schmidt decomposition, you can proceed e.g. as in Preskill's lecture notes: Diagonalize the reduced state of A, which yields eigenvalues $lambda_i$ and eigenvectors $|a_irangle$. Then, rewrite
$$
|psirangle = sum_i |a_irangleotimes |b_irangle .tag{*}
$$
($|b_irangle$ can be determined e.g. as $|b_irangle = langle a_i|psirangle$.) Then, the $|b_irangle$ are orthogonal with $langle b_i|b_irangle=lambda_i$ (cf. Preskill), i.e., the form $(*)$ above is the Schmidt decomposition (upon normalizing the $|b_irangle$).
(The latter can be seen by computing the reduced density matrix of A from $(*)$, which yields
$$
sum |a_irangle langle a_j | ; langle b_j|b_jrangle = sum lambda_i |a_irangle langle a_i| ,
$$
which yields $langle b_j|b_jrangle = lambda_idelta_{ij}$ as the $|a_iranglelangle a_j|$ are linearly independent.)
$endgroup$
add a comment |
$begingroup$
Norbert Schuch's answer is correct. Just for fun, here's the exact decomposition:
$$
|psiranglepropto
(3+sqrt{5})|Arangleotimes |Arangle
-(3-sqrt{5})|Brangleotimes |Brangle
$$
with
begin{align}
|Arangle &= 2|0rangle+(sqrt{5}-1)|1rangle \
|Brangle &= 2|0rangle-(sqrt{5}+1)|1rangle.
end{align}
Using these equations, we can verify that the coefficient of $|11rangle$ is zero and that the coefficients of $|00rangle$, $|01rangle$, and $|10rangle$ are all equal to each other, and
$$
langle A|Brangle = 0.
$$
Unnormalized vectors $A,B$ are used here to simplify the coefficients in the overall expression for $|psirangle$.
$endgroup$
add a comment |
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: "151"
};
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
});
}
});
Scott is a new contributor. Be nice, and check out our Code of Conduct.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fphysics.stackexchange.com%2fquestions%2f468380%2fschmidt-decomposition-example%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
$begingroup$
The phase of your eigenvectors is not correct (or rather, it is not determined, so you need to make a judicious choice). If you put a minus sign in front of the 2nd term, it works out.
To elaborate on that: If you want to find the Schmidt decomposition, you can proceed e.g. as in Preskill's lecture notes: Diagonalize the reduced state of A, which yields eigenvalues $lambda_i$ and eigenvectors $|a_irangle$. Then, rewrite
$$
|psirangle = sum_i |a_irangleotimes |b_irangle .tag{*}
$$
($|b_irangle$ can be determined e.g. as $|b_irangle = langle a_i|psirangle$.) Then, the $|b_irangle$ are orthogonal with $langle b_i|b_irangle=lambda_i$ (cf. Preskill), i.e., the form $(*)$ above is the Schmidt decomposition (upon normalizing the $|b_irangle$).
(The latter can be seen by computing the reduced density matrix of A from $(*)$, which yields
$$
sum |a_irangle langle a_j | ; langle b_j|b_jrangle = sum lambda_i |a_irangle langle a_i| ,
$$
which yields $langle b_j|b_jrangle = lambda_idelta_{ij}$ as the $|a_iranglelangle a_j|$ are linearly independent.)
$endgroup$
add a comment |
$begingroup$
The phase of your eigenvectors is not correct (or rather, it is not determined, so you need to make a judicious choice). If you put a minus sign in front of the 2nd term, it works out.
To elaborate on that: If you want to find the Schmidt decomposition, you can proceed e.g. as in Preskill's lecture notes: Diagonalize the reduced state of A, which yields eigenvalues $lambda_i$ and eigenvectors $|a_irangle$. Then, rewrite
$$
|psirangle = sum_i |a_irangleotimes |b_irangle .tag{*}
$$
($|b_irangle$ can be determined e.g. as $|b_irangle = langle a_i|psirangle$.) Then, the $|b_irangle$ are orthogonal with $langle b_i|b_irangle=lambda_i$ (cf. Preskill), i.e., the form $(*)$ above is the Schmidt decomposition (upon normalizing the $|b_irangle$).
(The latter can be seen by computing the reduced density matrix of A from $(*)$, which yields
$$
sum |a_irangle langle a_j | ; langle b_j|b_jrangle = sum lambda_i |a_irangle langle a_i| ,
$$
which yields $langle b_j|b_jrangle = lambda_idelta_{ij}$ as the $|a_iranglelangle a_j|$ are linearly independent.)
$endgroup$
add a comment |
$begingroup$
The phase of your eigenvectors is not correct (or rather, it is not determined, so you need to make a judicious choice). If you put a minus sign in front of the 2nd term, it works out.
To elaborate on that: If you want to find the Schmidt decomposition, you can proceed e.g. as in Preskill's lecture notes: Diagonalize the reduced state of A, which yields eigenvalues $lambda_i$ and eigenvectors $|a_irangle$. Then, rewrite
$$
|psirangle = sum_i |a_irangleotimes |b_irangle .tag{*}
$$
($|b_irangle$ can be determined e.g. as $|b_irangle = langle a_i|psirangle$.) Then, the $|b_irangle$ are orthogonal with $langle b_i|b_irangle=lambda_i$ (cf. Preskill), i.e., the form $(*)$ above is the Schmidt decomposition (upon normalizing the $|b_irangle$).
(The latter can be seen by computing the reduced density matrix of A from $(*)$, which yields
$$
sum |a_irangle langle a_j | ; langle b_j|b_jrangle = sum lambda_i |a_irangle langle a_i| ,
$$
which yields $langle b_j|b_jrangle = lambda_idelta_{ij}$ as the $|a_iranglelangle a_j|$ are linearly independent.)
$endgroup$
The phase of your eigenvectors is not correct (or rather, it is not determined, so you need to make a judicious choice). If you put a minus sign in front of the 2nd term, it works out.
To elaborate on that: If you want to find the Schmidt decomposition, you can proceed e.g. as in Preskill's lecture notes: Diagonalize the reduced state of A, which yields eigenvalues $lambda_i$ and eigenvectors $|a_irangle$. Then, rewrite
$$
|psirangle = sum_i |a_irangleotimes |b_irangle .tag{*}
$$
($|b_irangle$ can be determined e.g. as $|b_irangle = langle a_i|psirangle$.) Then, the $|b_irangle$ are orthogonal with $langle b_i|b_irangle=lambda_i$ (cf. Preskill), i.e., the form $(*)$ above is the Schmidt decomposition (upon normalizing the $|b_irangle$).
(The latter can be seen by computing the reduced density matrix of A from $(*)$, which yields
$$
sum |a_irangle langle a_j | ; langle b_j|b_jrangle = sum lambda_i |a_irangle langle a_i| ,
$$
which yields $langle b_j|b_jrangle = lambda_idelta_{ij}$ as the $|a_iranglelangle a_j|$ are linearly independent.)
edited yesterday
answered yesterday
Norbert SchuchNorbert Schuch
9,30722639
9,30722639
add a comment |
add a comment |
$begingroup$
Norbert Schuch's answer is correct. Just for fun, here's the exact decomposition:
$$
|psiranglepropto
(3+sqrt{5})|Arangleotimes |Arangle
-(3-sqrt{5})|Brangleotimes |Brangle
$$
with
begin{align}
|Arangle &= 2|0rangle+(sqrt{5}-1)|1rangle \
|Brangle &= 2|0rangle-(sqrt{5}+1)|1rangle.
end{align}
Using these equations, we can verify that the coefficient of $|11rangle$ is zero and that the coefficients of $|00rangle$, $|01rangle$, and $|10rangle$ are all equal to each other, and
$$
langle A|Brangle = 0.
$$
Unnormalized vectors $A,B$ are used here to simplify the coefficients in the overall expression for $|psirangle$.
$endgroup$
add a comment |
$begingroup$
Norbert Schuch's answer is correct. Just for fun, here's the exact decomposition:
$$
|psiranglepropto
(3+sqrt{5})|Arangleotimes |Arangle
-(3-sqrt{5})|Brangleotimes |Brangle
$$
with
begin{align}
|Arangle &= 2|0rangle+(sqrt{5}-1)|1rangle \
|Brangle &= 2|0rangle-(sqrt{5}+1)|1rangle.
end{align}
Using these equations, we can verify that the coefficient of $|11rangle$ is zero and that the coefficients of $|00rangle$, $|01rangle$, and $|10rangle$ are all equal to each other, and
$$
langle A|Brangle = 0.
$$
Unnormalized vectors $A,B$ are used here to simplify the coefficients in the overall expression for $|psirangle$.
$endgroup$
add a comment |
$begingroup$
Norbert Schuch's answer is correct. Just for fun, here's the exact decomposition:
$$
|psiranglepropto
(3+sqrt{5})|Arangleotimes |Arangle
-(3-sqrt{5})|Brangleotimes |Brangle
$$
with
begin{align}
|Arangle &= 2|0rangle+(sqrt{5}-1)|1rangle \
|Brangle &= 2|0rangle-(sqrt{5}+1)|1rangle.
end{align}
Using these equations, we can verify that the coefficient of $|11rangle$ is zero and that the coefficients of $|00rangle$, $|01rangle$, and $|10rangle$ are all equal to each other, and
$$
langle A|Brangle = 0.
$$
Unnormalized vectors $A,B$ are used here to simplify the coefficients in the overall expression for $|psirangle$.
$endgroup$
Norbert Schuch's answer is correct. Just for fun, here's the exact decomposition:
$$
|psiranglepropto
(3+sqrt{5})|Arangleotimes |Arangle
-(3-sqrt{5})|Brangleotimes |Brangle
$$
with
begin{align}
|Arangle &= 2|0rangle+(sqrt{5}-1)|1rangle \
|Brangle &= 2|0rangle-(sqrt{5}+1)|1rangle.
end{align}
Using these equations, we can verify that the coefficient of $|11rangle$ is zero and that the coefficients of $|00rangle$, $|01rangle$, and $|10rangle$ are all equal to each other, and
$$
langle A|Brangle = 0.
$$
Unnormalized vectors $A,B$ are used here to simplify the coefficients in the overall expression for $|psirangle$.
edited yesterday
answered yesterday
Chiral AnomalyChiral Anomaly
12.6k21542
12.6k21542
add a comment |
add a comment |
Scott is a new contributor. Be nice, and check out our Code of Conduct.
Scott is a new contributor. Be nice, and check out our Code of Conduct.
Scott is a new contributor. Be nice, and check out our Code of Conduct.
Scott is a new contributor. Be nice, and check out our Code of Conduct.
Thanks for contributing an answer to Physics 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.
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
StackExchange.ready(
function () {
StackExchange.openid.initPostLogin('.new-post-login', 'https%3a%2f%2fphysics.stackexchange.com%2fquestions%2f468380%2fschmidt-decomposition-example%23new-answer', 'question_page');
}
);
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Post as a guest
Required, but never shown
Sign up or log in
StackExchange.ready(function () {
StackExchange.helpers.onClickDraftSave('#login-link');
});
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
Sign up using Google
Sign up using Facebook
Sign up using Email and Password
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