Chirality-induced spin selectivity (CISS) is a general term denoting the interplay between the chiral structure of molecules and the electron spin. CISS has been studied for over two decades, leading to a consistent picture of experimental phenomena. In this thesis, we focus on transport experiments exhibiting CISS. In spite of the efforts of many scientists, there is no theoretical explanation for the high degrees of CISS that have been measured in electron transport experiments. In general, it is agreed upon by theorists that the effect is caused by the interplay between the electronic spin-orbit interaction and the helical molecule geometry, in combination with phase-breaking effects such as electron-electron or electron-phonon interactions. As of yet, no theoretical model has achieved realistic degrees of SOC without drastic inflation of the spin-orbit interaction strength.
In this thesis we explore the possibility of using matrix product states (MPS) to study spin-selectivity in boundary-driven electron transport through tight-binding models of chiral molecules, a novel approach in the field of CISS. To this end, we use a model proposed by Fransson in 2019, which considers interacting electrons in a Hubbard model with a spin-orbit interaction adapted from the Kane-Mele model. In this thesis work, the fermionic Hubbard model is mapped to a double spin chain using the Jordan-Wigner transformation. The state of the system is described by a matrix product density operator (MPDO) which is vectorised to a matrix product state (MPS). The system dynamics are described by a vectorised Lindblad equation. The advantage of this approach lies in the fact that it does not require the use of any systematic approximations to the Hamiltonian, in contrast to previous studies. The developed method is validated against the results of previous works studying the boundary-driven Heisenberg-XXZ model and the boundary-driven Hubbard model.
The method is shown to be capable of reproducing chirality-induced spin-selective effects for short chains. The results of this study show a finite magnetocurrent that is odd in bias voltage with an associated magnetoresistance of less than 1%. This is in line with previous studies of this model, but two orders of magnitude lower than experimentally measured values. However, these results are obtained using highly inflated values for the spin-orbit interaction strength. In multiple cases, the results do not satisfy the Onsager-Casimir and Büttiker reciprocity principles, which state that the magnetocurrent should vanish in the low-driving and in the non-interacting regimes. Moreover, the continuity of the current in the steady state was not fully satisfied.
We provide evidence that indicates these problems result from the time-integration error introduced by the Suzuki-Trotter decomposition. We expect that these can be mitigated using higher order time-integration schemes. From the results of this study we can conclude that matrix product states are a viable tool to study CISS in bound-electron transport. However, the method presented in this thesis suffer from numerical errors. We present several suggestions for improvement which address these shortcomings.

In this thesis we implement the Lindblad equation in the matrix product state (MPS) formalism using an operator splitting method. We developed a second-order method based on a Trotter approximation and a third-order high-dimensional midpoint method and we proposed a new fourth-order method based on Duhamel’s principle and a nested RK4 method, all of which preserve positivity and Hermiticity of the density operator. We simulated spin transport through an XXZ-Hamiltonian Heisenberg chain, for which we found the magnetisation profile and measured a spin current of 0.04-0.05. The results obtained are consistent with the existing literature. The extensive error analysis shows that the time step ∆t is the main contributor to the error, if the bond dimension χ is set to at least 15. The third-order method is in general preferred to the second-order method, as only this method preserves trace. We also analysed the Hubbard model, including a spin orbit coupling, in order to propose a method for simulating the chiral induced spin selectivity (CISS) effect.

The behaviour of the dilute Curie-Weiss model has been analysed for Gaussian, bounded and subGaussian couplings. This model is an abstraction of the Curie-Weiss model which is a model of a spin configuration where instead of a constant coupling between spins the spin coupling is a realisation of a random variable. To analyse this behaviour first the behaviour of the standard Curie-Weiss model is depicted. A notable part from this behaviour is the presence of a phase transition., indicating that this model can be used for a study of phase transitions. After having analysed the behaviour of the standard Curie-Weiss model theorems, stating the closeness between the standard Curie-Weiss model and the randomly dilute Curie-Weiss model for Gaussian, bounded and sub-Gaussian couplings, will be proven. Indicating that the models behave approximately the same way. These theorems will be in the form of a bound for the on randomly dilute Curie-Weiss model as an exponential multiplied by the standard Curie-Weiss model with the probability of this bound being a sub-Gaussian distribution.