A prominent tool to study the dynamics of open quantum systems is the reduced density matrix. Yet, approaching open quantum systems by means of state vectors has well known computational advantages. In this respect, the physical meaning of the so-called conditional states in Markovian and non-Markovian scenarios has been a topic of recent debate in the construction of stochastic Schrödinger equations. We shed light on this discussion by acknowledging the Bohmian conditional wavefunction (linked to the corresponding Bohmian trajectory) as the proper mathematical object to represent, in terms of state vectors, an arbitrary subset of degrees of freedom. As an example of the practical utility of these states, we present a time-dependent quantum Monte Carlo algorithm to describe electron transport in open quantum systems under general (Markovian or non-Markovian) conditions. By making the most of trajectory-based and wavefunction methods, the resulting simulation technique extends to the quantum regime, the computational capabilities that the Monte Carlo solution of the Boltzmann transport equation offers for semi-classical electron devices.
A prominent tool to study the dynamics of open quantum systems is the reduced density matrix. Yet, approaching open quantum systems by means of state vectors has well known computational advantages. In this respect, the physical meaning of the so-called conditional states in Markovian and non-Markovian scenarios has been a topic of recent debate in the construction of stochastic Schrödinger equations. We shed light on this discussion by acknowledging the Bohmian conditional wavefunction (linked to the corresponding Bohmian trajectory) as the proper mathematical object to represent, in terms of state vectors, an arbitrary subset of degrees of freedom. As an example of the practical utility of these states, we present a time-dependent quantum Monte Carlo algorithm to describe electron transport in open quantum systems under general (Markovian or non-Markovian) conditions. By making the most of trajectory-based and wavefunction methods, the resulting simulation technique extends to the quantum regime, the computational capabilities that the Monte Carlo solution of the Boltzmann transport equation offers for semi-classical electron devices.
A prominent tool to study the dynamics of open quantum systems is the reduced density matrix. Yet, approaching open quantum systems by means of state vectors has well known computational advantages. In this respect, the physical meaning of the so-called conditional states in Markovian and non-Markovian scenarios has been a topic of recent debate in the construction of stochastic Schrödinger equations. We shed light on this discussion by acknowledging the Bohmian conditional wavefunction (linked to the corresponding Bohmian trajectory) as the proper mathematical object to represent, in terms of state vectors, an arbitrary subset of degrees of freedom. As an example of the practical utility of these states, we present a time-dependent quantum Monte Carlo algorithm to describe electron transport in open quantum systems under general (Markovian or non-Markovian) conditions. By making the most of trajectory-based and wavefunction methods, the resulting simulation technique extends to the quantum regime, the computational capabilities that the Monte Carlo solution of the Boltzmann transport equation offers for semi-classical electron devices.
Without access to the full quantum state, modeling dissipation in an open system requires approximations. The physical soundness of such approximations relies on using realistic microscopic models of dissipation that satisfy completely positive dynamical maps. Here we present an approach based on the use of the Bohmian conditional wave function that, by construction, ensures a completely positive dynamical map for either Markovian or non-Markovian scenarios while allowing the implementation of realistic dissipation sources. Our approach is applied to compute the current-voltage characteristic of a resonant tunneling device with a parabolic-band structure, including electron-lattice interactions. A stochastic Schr\"odinger equation is solved for the conditional wave function of each simulated electron. We also extend our approach to (graphenelike) materials with a linear band structure using Bohmian conditional spinors for a stochastic Dirac equation.
After the seminal work of R. Landauer in 1957 relating the electrical resistance of a conductor to its scattering properties, much progress has been made in our ability to predict the performance of electron devices in the DC (stationary) regime. Computational tools to describe their dynamical behavior (including the AC, transient and noise performance), however, are far from being as trustworthy as would be desired by the electronic industry. While there is no fundamental limitation to correctly modeling the high-frequency quantum transport and its fluctuations, certainly more careful attention must be paid to delicate issues such as overall charge neutrality, total current conservation, or the back action of the measuring apparatus. In this review, we will show how the core ideas behind the Bohmian formulation of quantum mechanics can be exploited to design an efficient Monte Carlo algorithm that provides a quantitative description of electron transport in open quantum systems. By making the most of trajectory-based and wave function methods, the BITLLES simulator, a free software developed by the authors, extends the capabilities that the semi-classical Monte Carlo simulation technique has offered for decades (DC, AC, noise, transients) to the quantum regime.
We present a stochastic projection formalism for the description of quantum dynamics in bosonic or spin environments. The Schrödinger equation in the coherent state representation with respect to the environmental degrees of freedom can be reformulated by employing the Feshbach partitioning technique for open quantum systems based on the introduction of suitable non-Hermitian projection operators. In this picture the reduced state of the system can be obtained as a stochastic average over pure state trajectories, for any temperature of the bath. The corresponding non-Markovian stochastic Schrödinger equations include a memory integral over the past states. In the case of harmonic environments and linear coupling the approach gives a new form of the established non-Markovian quantum state diffusion stochastic Schrödinger equation without functional derivatives. Utilizing spin coherent states, the evolution equation for spin environments resembles the bosonic case with, however, a non-Gaussian average for the reduced density operator.
Moore’s law has been the milestone for electronics improvement, and the cause for the exponential growth in our computational abilities and for reaching nanoscale electron devices. At such dimensions, classical simulation tools must be substituted by quantum ones. Such substitution implies to tackle new fundamental problems. First, the many-body problem makes (almost) impossible the exact simulation of many particles scenarios. Second, the measurement problem is especially relevant in high-frequency scenarios where multi-time measurements are mandatory. Alternatively to the orthodox theory, Bohmian mechanics emerges as a quantum theory which is specially well-equipped for the simulation of high-frequency characteristic of quantum electron devices. The Bohmian theory provides the conditional wave function that guides particles with well-defined positions. In this thesis, I explored quantum transport using Bohmian mechanics, putting special emphasis to graphene, a 2D material with linear bandstructure, which is expected to play an important role in the next future electronics. Differently from other materials, graphene obeys the bispinor Dirac equation, electrons behave as massless particles, exhibiting exotic behaviors, such as the Klein Tunneling effect. During the thesis, the quantum BITLLES simulator has been improved to correctly model electron nanodevices with either linear (by the inclusion of the complex bispinor Dirac equation) or parabolic bandstructures in either ballistic or dissipative (by the inclusion of the complete positive Bohmian scattering approach) systems. Thus, the BITLLES has become the candidate for substituting the versatile semiclassical Monte Carlo approach in the quantum regime, while keeping the versatility to predict DC, AC and noise performances. As practical applications, during the thesis, unexpected scattering probabilities in a Hong-Ou-Mandel experiment were predicted. In addition, different current-voltage characteristics were analyzed. Finally, a new limit (due to the discrete nature of charge) for ultra-small logic applications working at THz frequencies was predicted.
Evaluation of the charge transport property of organic semiconductors requires exact quantum dynamics simulation of large systems. We present a numerically nearly exact approach to investigate carrier transport dynamics in organic semiconductors by extending the non-Markovian stochastic Schrödinger equation with complex frequency modes to a forward-backward scheme and by solving it using the matrix product state (MPS) approach. By utilizing the forward-backward formalism for noise generation, the bath correlation function can be effectively treated as a temperature-independent imaginary part, enabling a more accurate decomposition with fewer complex frequency modes. Using this approach, we study the carrier transport and mobility in the one-dimensional Peierls model, where the nonlocal electron-phonon interaction is taken into account. The reliability of this approach was validated by comparing carrier diffusion motion with those obtained from the hierarchical equations of motion method across various parameter regimes of the phonon bath. The efficiency was demonstrated by the modest virtual bond dimensions of MPS and the low scaling of the computational time with the system size.