Debmalya Das,George Thomas,Andrew N. Jordan
Abstract:In a quantum Stirling heat engine, the heat exchanged with two thermal baths is partly utilized for performing work by redistributing the energy levels of the working substance. We analyze the thermodynamics of a quantum Stirling engine operating in finite time. We develop a model in which a time-dependent potential barrier changes the energy-level structure of the working substance. The process takes place under a constant interaction with the thermal bath. We further show that in the limit of slow operation of the cycle and low temperature, the efficiency of such an engine approaches Carnot efficiency. We also show that the maximum output power , for the strokes that affect the energy levels, is obtained at an intermediate operating speed, demonstrating the importance of a finite-time analysis.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is the thermodynamic performance analysis of a quantum Stirling heat engine operating within a limited time. Specifically, the researchers developed a model in which part of the thermal energy is converted into work by changing the energy - level structure of the working substance. They introduced a time - dependent barrier, which changes the energy - level structure of the working substance, and the whole process is carried out while maintaining a constant interaction with the heat bath. The paper pays special attention to the situation where the efficiency of this engine approaches the Carnot efficiency in the low - temperature and slow - operation limits, and shows that in the stroke affecting the energy level, the maximum output power is obtained at the intermediate operation speed, emphasizing the importance of the finite - time analysis.
### Core problems of the paper:
1. **Influence of finite - time operation**: Explore the influence of operation within a limited time on the performance of the quantum Stirling heat engine, especially how it affects the efficiency of converting thermal energy into work and the output power.
2. **Approximation of Carnot efficiency**: Analyze whether the efficiency of the quantum Stirling heat engine can approach the Carnot efficiency in the low - temperature and slow - operation limits.
3. **Optimization of maximum output power**: Determine under which operating conditions the quantum Stirling heat engine can reach the maximum output power, especially its performance at different operating speeds.
### Main contributions:
- **Model development**: Propose a new model to change the energy - level structure of the working substance through a time - dependent barrier, thereby achieving the conversion from thermal energy to work.
- **Thermodynamic analysis**: Analyze in detail the thermodynamic performance during operation within a limited time, including the changes in efficiency and output power.
- **Numerical simulation**: Verify the validity of the model through numerical simulation and show the performance changes at different operating speeds.
### Key formulas:
- **Hamiltonian**:
\[
H(x, t)=-\frac{\hbar^{2}}{2m}\frac{\partial^{2}}{\partial x^{2}}+\alpha(t)\delta(x)
\]
where \(\alpha(t)\) is the time - dependent barrier height.
- **Time evolution of the density matrix**:
\[
\dot{\rho}=-\frac{i}{\hbar}[H, \rho]+\sum_{k}\gamma_{k}(N_{\omega_{k}} + 1)\left(L_{k}\rho L_{k}^{\dagger}-\frac{1}{2}(L_{k}^{\dagger}L_{k}\rho+\rho L_{k}^{\dagger}L_{k})\right)+\sum_{k}\gamma_{k}N_{\omega_{k}}\left(L_{k}^{\dagger}\rho L_{k}-\frac{1}{2}(L_{k}L_{k}^{\dagger}\rho+\rho L_{k}L_{k}^{\dagger})\right)
\]
where \(L_{k}\) is the jump operator, \(\gamma_{k}\) is the decay rate, and \(N_{\omega_{k}}\) is the Bose - Einstein distribution function.
- **Heat exchange and work done**:
\[
\Delta Q(t, t+\Delta t)=\text{Tr}[H(t)(\rho(t+\Delta t)-\rho(t))]
\]
\[
\Delta W(t, t+\Delta t')=\text{Tr}[(H(t+\Delta t') - H(t))\rho(t)]
\]
- **Total heat exchange and total work done**:
\[
Q_{\text{ins}}=\sum_{i = 1}^{n\Delta t}\Delta Q(t_{i}, t_{i}+\Delta t)
\]
\[
W_{\text{ins}}=\sum_{i = 1}^{n\Delta t}\Delta W(t_{i}, t_{i}+\Delta t)
\]