Current-induced bond rupture in single-molecule junctions: Effects of multiple electronic states and vibrational modes

Yaling Ke,Jan Dvořák,Martin Čížek,Raffaele Borrelli,Michael Thoss

DOI: https://doi.org/10.1063/5.0155290

2023-04-19

Abstract:Current-induced bond rupture is a fundamental process in nanoelectronic architectures such as molecular junctions and in scanning tunneling microscopy measurements of molecules at surfaces. The understanding of the underlying mechanisms is important for the design of molecular junctions that are stable at higher bias voltages and is a prerequisite for further developments in the field of current-induced chemistry. In this work, we analyse the mechanisms of current-induced bond rupture employing a recently developed method, which combines the hierarchical equations of motion approach in twin space with the matrix product state formalism, and allows accurate, fully quantum mechanical simulations of the complex bond rupture dynamics. Extending previous work [J. Chem. Phys. 154, 234702 (2021)], we consider specifically the effect of multiple electronic states and multiple vibrational modes. The results obtained for a series of models of increasing complexity show the importance of vibronic coupling between different electronic states of the charged molecule, which can enhance the dissociation rate at low bias voltages profoundly.

Mesoscale and Nanoscale Physics,Chemical Physics

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to understand the mechanism of current - induced bond breaking in single - molecule junctions, especially how these mechanisms change when there are multiple electronic states and multiple vibrational modes. Specifically:

1. **Basic mechanism of current - induced bond breaking**: The author has previously studied the mechanism of current - induced bond breaking in the minimal model containing one electronic state and a single vibrational reaction mode. However, in polyatomic molecules, multiple electronic states and multiple vibrational modes are usually involved in the reaction mechanism, which makes the situation more complicated.

2. **Considering the influence of multiple electronic states and multiple vibrational modes**: In this work, the author extends the previous research and considers the influence of multiple electronic states and multiple vibrational modes on current - induced bond breaking. The research shows that the vibrational coupling between the electronic states of different charged molecules can significantly enhance the dissociation rate at lower bias voltages.

3. **Developing new simulation methods**: In order to accurately simulate the complex bond - breaking dynamics, the author adopts a new method that combines the Hierarchical Equations of Motion (HEOM) method with the Matrix Product State (MPS) form. This method allows for fully quantum - mechanical simulations and can handle complex systems containing multiple electronic states and vibrational modes.

### Formula summary

- **System Hamiltonian**:
\[
H = H_S + H_B + H_{SB}
\]
where,
- \(H_S\) is the system Hamiltonian that describes the molecule.
- \(H_B\) is the bath Hamiltonian that describes the macroscopic electrode.
- \(H_{SB}\) is the coupling Hamiltonian between the molecule and the electrode.

- **Specific expression of the system Hamiltonian**:
\[
H_S = T_{\text{nuc}} + V_0(x_1, x_2) + \epsilon_1(x_1)d_1^{\dagger}d_1 + \epsilon_2(x_1)d_2^{\dagger}d_2 + U d_1^{\dagger}d_1d_2^{\dagger}d_2 + \Delta(x_2)(d_1^{\dagger}d_2 + d_2^{\dagger}d_1)
\]

- **Potential energy surface**:
\[
V_0(x_1, x_2) = D_e(1 - e^{-a(x_1 - x_0)})^2+\frac{1}{2}\omega_bx_2^2
\]

- **Coupling parameter**:
\[
\Delta(x_2)=\Delta_0+\Delta_1x_2e^{-\lambda x_2^2}
\]

Through these formulas and methods, the author aims to reveal the design principles for the molecular junction to remain stable at higher bias voltages and provide theoretical support for the development of nanoscale chemical catalysis.