For a system without spin–orbit coupling, the (i) nuclear plus electronic linear momentum and (ii) nuclear plus orbital electronic angular momentum are good quantum numbers. Thus, when a molecular system undergoes a nonadiabatic transition, there should be no change in the total linear or angular momentum. Now, the standard surface hopping algorithm ignores the electronic momentum and indirectly equates the momentum of the nuclear degrees of freedom to the total momentum. However, even with this simplification, the algorithm still does not conserve either the nuclear linear or the nuclear angular momenta. Here, we show that one way to address these failures is to dress the derivative couplings (i.e., the hopping directions) in two ways: (i) we disallow changes in the nuclear linear momentum by working in a translating basis (which is well known and leads to electron translation factors) and (ii) we disallow changes in the nuclear angular momentum by working in a basis that rotates around the center of mass [which is not well-known and leads to a novel, rotationally removable component of the derivative coupling that we will call electron rotation factors below, cf. Eq. (96)]. The present findings should be helpful in the short term as far as interpreting surface hopping calculations for singlet systems (without spin) and then developing the new surface hopping algorithm in the long term for systems where one cannot ignore the electronic orbital and/or spin angular momentum.
This content will become publicly available on March 28, 2025
Within the context of fewest-switch surface hopping (FSSH) dynamics, one often wishes to remove the angular component of the derivative coupling between states J and K. In a previous set of papers, Shu et al. [J. Phys. Chem. Lett. 11, 1135–1140 (2020)] posited one approach for such a removal based on direct projection, while we isolated a second approach by constructing and differentiating a rotationally invariant basis. Unfortunately, neither approach was able to demonstrate a one-electron operatorÔ whose matrix element JÔK was the angular component of the derivative coupling. Here, we show that a one-electron operator can, in fact, be constructed efficiently in a semi-local fashion. The present results yield physical insight into designing new surface hopping algorithms and are of immediate use for FSSH calculations.
more » « less- Award ID(s):
- 2102402
- NSF-PAR ID:
- 10520187
- Publisher / Repository:
- jcp
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 12
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
more » « less
-
null (Ed.)Pinhole-free ultrathin films of polyelectrolyte complex assembled using layer-by-layer deposition were used to evaluate electron transfer from a redox species in solution to an electrode over the distance range 1 to 9 nm. Over this thickness, the polyelectrolytes employed wet the surface and the polymer molecules flatten to less than their equilibrium size in 3-dimensions. A decay constant β for current as a function of distance of about 0.3 nm-1 placed this system in the regime expected for multistep hopping versus a one-step tunneling event. Discreet hopping sites within the films were identified as ferrocyanide ions with an equilibrium concentration of 0.032 M and an average separation of 3.7 nm. The Butler-Volmer (BV) expression for electron transfer as a function of overpotential was modified by distributing the applied voltage evenly amongst the hopping sites. This modified BV expression fit both the distance dependence and the applied potential dependence well, wherein the only freely adjustable parameter was the electron transfer coefficient. The finding that β is simply the inverse of the hopping range is consistent with previous conclusions that electrons within conjugated molecule sites are delocalized, or, for non-conjugated systems, spread over more than one repeat unit by lattice distortionsmore » « less
-
A generalized quantum master equation approach is introduced to describe electron transfer in molecular junctions that spans both the off-resonant (tunneling) and resonant (hopping) transport regimes. The model builds on prior insights from scattering theory but is not limited to a certain parameter range with regard to the strength of the molecule–electrode coupling. The framework is used to study the simplest case of energy and charge transfer between the molecule and the electrodes for a single site noninteracting Anderson model in the limit of symmetric and asymmetric coupling between the molecule and the electrodes. In the limit of elastic transport, the Landauer result is recovered for the current by invoking a single active electron Ansatz and a binary collision approximation for the memory kernel. Inelastic transport is considered by allowing the excitation of electron–hole pairs in the electrodes in tandem with charge transport. In the case of low bias voltages where the Fermi levels of the electrodes remain below the molecular state, it is shown that the current arises from tunneling and the molecule remains neutral. However, once the threshold is reached for aligning the fermi level of one electrode with the molecular orbital, a small amount of charge transfer occurs with a negligible amount of hopping current. While inelasticity in the current has a minimal impact on the shape of the current–voltage curve in the case of symmetric electrode coupling, the results for a slight asymmetry in coupling demonstrate complete charge transfer and a significant drop in current. These results provide encouraging confirmation that the framework can describe charge transport across a wide range of electrode–molecule coupling and provide a unique perspective for developing new master equation treatments for energy and charge transport in molecular junctions. An extension of this work to account for inelastic scattering from electron–vibrational coupling at the molecule is straightforward and will be the subject of subsequent work.more » « less
-
more » « less
Abstract A microtubule is a cylindrical biological polymer that plays key roles in cellular structure, transport, and signalling. In this work, based on studies of electronic properties of polyacetelene and mechanical properties of microtubules themselves (Spakowitz A. J.,
Phys. Rev. Lett. ,103 (2009) 248101), we explore the possibility that microtubules could act as topological insulators that are gapped to electronic excitations in the bulk but possess robust electronic bounds states at the tube ends. Through analyses of structural and electronic properties, we model the microtubule as a cylindrical stack of Su-Schrieffer-Heeger chains (originally proposed in the context of polyacetylene) describing electron hopping between the underlying dimerized tubulin lattice sites. We postulate that the microtubule is mostly uniform, dominated purely by GDP-bound dimers, and is capped by a disordered regime due to the presence of GTP-bound dimers as well. In the uniform region, we identify the electron hopping parameter regime in which the microtubule is a topological insulator. We then show the manner in which these topological features remain robust when the hopping parameters are disordered. We briefly mention possible biological implications for these microtubules to possess topologically robust electronic bound states. -
more » « less
We demonstrate that, for systems with spin–orbit coupling and an odd number of electrons, the standard fewest switches surface hopping algorithm does not conserve the total linear or angular momentum. This lack of conservation arises not so much from the hopping direction (which is easily adjusted) but more generally from propagating adiabatic dynamics along surfaces that are not time reversible. We show that one solution to this problem is to run along eigenvalues of phase-space electronic Hamiltonians H(R, P) (i.e., electronic Hamiltonians that depend on both nuclear position and momentum) with an electronic–nuclear coupling Γ · P [see Eq. (25)], and we delineate the conditions that must be satisfied by the operator Γ. The present results should be extremely useful as far as developing new semiclassical approaches that can treat systems where the nuclear, electronic orbital, and electronic spin degrees of freedom altogether are all coupled together, hopefully including systems displaying the chiral-induced spin selectivity effect.
