We study current-induced vibrational cooling, heating, and instability in a donor-acceptor rectifying molecular junction using a full counting statistics approach. In our model, electron-hole pair excitations are coupled to a given molecular vibrational mode which is either harmonic or highly anharmonic. This mode may be further coupled to a dissipative thermal environment. Adopting a master equation approach, we confirm the charge and heat exchange fluctuation theorem in the steady-state limit, for both harmonic and anharmonic models. Using simple analytical expressions, we calculate the charge current and several measures for the mode effective temperature. At low bias, we observe the effect of bias-induced cooling of the vibrational mode. At higher bias, the mode effective temperature is higher than the environmental temperature, yet the junction is stable. Beyond that, once the vibrational mode (bias-induced) excitation rate overcomes its relaxation rate, instability occurs. We identify regimes of instability as a function of voltage bias and coupling to an additional phononic thermal bath. Interestingly, we observe a reentrant behavior where an unstable junction can properly behave at a high enough bias. The mechanism for this behavior is discussed.

No abstract

We expand iterative numerically-exact influence functional path-integral tools and present a method capable of following the nonequilibrium time evolution of subsystems coupled to multiple bosonic and fermionic reservoirs simultaneously. Using this method, we study the real-time dynamics of charge transfer and vibrational mode excitation in an electron conducting molecular junction. We focus on nonequilibrium vibrational effects, particularly, the development of vibrational instability in a current-rectifying junction. Our simulations are performed by assuming large molecular vibrational anharmonicity (or low temperature). This allows us to truncate the molecular vibrational mode to include only a two-state system. Exact numerical results are compared to perturbative Master equation calculations demonstrating an excellent agreement in the weak electron-phonon coupling regime. Significant deviations take place only at strong coupling. Our simulations allow us to quantify the contribution of different transport mechanisms, coherent dynamics and inelastic transport, in the overall charge current. This is done by studying two model variants: The first admits inelastic electron transmission only, while the second one allows for both coherent and incoherent pathways.

Key words non-Markovian decay, power law decay, bound stateIt is known that quantum systems yield non-exponential (power law) decay on long time scales, associated with continuum threshold effects contributing to the survival probability for a prepared initial state. For an open quantum system consisting of a discrete state coupled to continuum, we study the case in which a discrete bound state of the full Hamiltonian approaches the energy continuum as the system parameters are varied. We find in this case that at least two regions exist yielding qualitatively different power law decay behaviors; we term these the long time 'near zone' and long time 'far zone.' In the near zone the survival probability falls off according to a t −1 power law, and in the far zone it falls off as t −3 . We show that the timescale TQ separating these two regions is inversely related to the gap between the discrete bound state energy and the continuum threshold. In the case that the bound state is absorbed into the continuum and vanishes, then the time scale TQ diverges and the survival probability follows the t −1 power law even on asymptotic scales. Conversely, one could study the case of an anti-bound state approaching the threshold before being ejected from the continuum to form a bound state. Again the t −1 power law dominates precisely at the point of ejection.

The backbone conformation of conjugated polymers affects, to a large extent, their optical and electronic properties. The usually flexible substituents provide solubility and influence the packing behavior of conjugated polymers in films or in bad solvents. However, the role of the side chains in determining and potentially controlling the backbone conformation, and thus the optical and electronic properties on the single polymer level, is currently under debate. Here, we investigate directly the impact of the side chains by studying the bulky-substituted poly(3-(2,5-dioctylphenyl)thiophene) (PDOPT) and the common poly(3-hexylthiophene) (P3HT), both with a defined molecular weight and high regioregularity, using low-temperature single-chain photoluminescence (PL) spectroscopy and quantum-classical simulations. Surprisingly, the optical transition energy of PDOPT is significantly (∼2,000 cm or 0.25 eV) red-shifted relative to P3HT despite a higher static and dynamic disorder in the former. We ascribe this red shift to a side-chain induced backbone planarization in PDOPT, supported by temperature-dependent ensemble PL spectroscopy. Our atomistic simulations reveal that the bulkier 2,5-dioctylphenyl side chains of PDOPT adopt a clear secondary helical structural motif and thus protect conjugation, i.e., enforce backbone planarity, whereas, for P3HT, this is not the case. These different degrees of planarity in both thiophenes do not result in different conjugation lengths, which we found to be similar. It is rather the stronger electronic coupling between the repeating units in the more planar PDOPT which gives rise to the observed spectral red shift as well as to a reduced calculated electron-hole polarization.

The optoelectronic properties of amorphous conjugated polymers are sensitive to the details of the conformational disorder, and spectroscopy provides the means for structural characterization of the fragments of the chain that interact with light-"chromophores". A faithful interpretation of spectroscopic conformational signatures, however, presents a theoretical challenge. Here we investigate the relationship between the ground-state optical gaps, the properties of the excited states, and the structural features of chromophores of a single molecule poly(3-hexyl)-thiophene (P3HT) using quantum-classical atomistic simulations. Our results demonstrate that chromophoric disorder arises through the interplay between excited-state delocalization and electron-hole polarization, controlled by the torsional disorder introduced by side chains. Within this conceptual framework, we predict and explain the counterintuitive spectral behavior of P3HT, a red-shifted absorption, despite shortening of chromophores, with increasing temperature. This discussion introduces the concept of disorder-induced separation of charges in amorphous conjugated polymers.

Significance It is of fundamental importance to establish the microscopic picture of crystal nucleation for both academic and industrial research. Here, we report an in situ time-resolved Raman spectroscopy to study crystallization from solution, one crystal nucleation at a time. The observed dynamics of α-glycine crystallization from water without any additive supports the nonclassical nucleation pathway, where prenucleation aggregates form and convert to a crystal. By the direct comparison of Raman spectrum of the aggregates between the experiment and simulation, we propose the structure of these aggregates as linear hydrogen-bonded networks. We demonstrate the power of studying one nucleation event at a time, which can accelerate the investigation of crystal nucleation by optical spectroscopy.

We look for manifestations of quantum interference effects in the Seebeck coefficient of a molecular junction, when the electronic conductance exhibits pronounced destructive interference features due to the presence of quasi-degenerate electronic states which differ in their spatial symmetry. We perform our analysis by considering three separate limits for electron transport: coherent, fully dephased, and suffering inelastic scattering with molecular vibrations. We find that while the conductance displays strong signatures of the underling transport mechanisms: destructive quantum interference features in the coherent case and thermal activation characteristics in the inelastic limit, the Seebeck coefficient conceals details of electron dynamics while it robustly reveals information about the energy characteristics of the junction. We provide closed-form expressions for the electronic conductance and the thermopower of our system as a function of temperature, gate voltage and hybridization energy in different transport limits, then exemplify our analysis on a specific conjugated molecule with quasi-degenerate orbitals of different spatial symmetry.