Define the coupled two-level systems:
Then run the Schmidt decompotions:
We can easily print the modes:
## [[1]]
## [[1]]$eigenvalue
## [1] 0.7
##
## [[1]]$sys1vector
## [1] -1 0
##
## [[1]]$sys2vector
## [1] 0 -1
##
##
## [[2]]
## [[2]]$eigenvalue
## [1] 0.3
##
## [[2]]$sys1vector
## [1] 0 -1
##
## [[2]]$sys2vector
## [1] -1 0
And their corresponding entanglement entropy:
## [1] 0.6108643