Modeling the electroluminescence of atomic wires from quantum dynamics simulations.

Static and time-dependent quantum-mechanical approaches have been employed in the literature to characterize the physics of light-emitting molecules and nanostructures. However, the electromagnetic emission induced by an input current has remained beyond the realm of molecular simulations. This is the challenge addressed here with the help of an equation of motion for the density matrix coupled to a photon bath based on a Redfield formulation. This equation is evolved within the framework of the driven-Liouville von Neumann approach, which incorporates open boundaries by introducing an applied bias and a circulating current. The dissipated electromagnetic power can be computed in this context from the time derivative of the energy. This scheme is applied in combination with a self-consistent tight-binding Hamiltonian to investigate the effects of bias and molecular size on the electroluminescence of metallic and semiconducting chains. For the latter, a complex interplay between bias and molecular length is observed: there is an optimal number of atoms that maximizes the emitted power at high voltages but not at low ones. This unanticipated behavior can be understood in terms of the band bending produced along the semiconducting chain, a phenomenon that is captured by the self-consistency of the method. A simple analytical model is proposed that explains the main features revealed by the simulations. The methodology, applied here at a self-consistent tight-binding level but extendable to more sophisticated Hamiltonians such as density functional tight binding and time dependent density functional theory, promises to be helpful for quantifying the power and quantum efficiency of nanoscale electroluminescent devices.

Coordinate-Free and Low-Order Scaling Machine Learning Model for Atomic Partial Charge Prediction for Any Size of Molecules.

The atomic partial charge is of great importance in many fields, such as chemistry and drug-target recognition. However, conventional quantum-based computing of atomic charges is relatively slow, limiting further applications of atomic charge analysis. With the help of machine learning methods, various kinds of models appear to speed up atomic charge calculations. However, there are still some concerning problems. Some models based on geometric coordinates require high-accuracy geometry optimization as a preprocess, while other models have a limitation on the size of input molecules that narrow the applications of the model. Here, we propose a machine learning atomic charge model based on a message-passing featurizer. This preprocessing featurizer can quickly extract atomic environment information from a molecule according to the connectivity inside the molecule. The resulting descriptor can be used with a neural network to quickly predict the atomic partial charge. The model is able to automatically adapt to any size of molecule while remaining efficient and achieves a root-mean-square error in the Hirshfeld charge prediction of 0.018e, with an overall time complexity of O(n2). Thus, this model could enlarge the range of applications of atomic partial charge to more fields and cases.

Revealing Interface Polarization Effects on the Electrical Double Layer with Efficient Open Boundary Simulations under Potential Control.

A major challenge in modeling interfacial processes in electrochemical (EC) devices is performing simulations at constant potential. This requires an open-boundary description of the electrons, so that they can enter and leave the computational cell. To enable realistic modeling of EC processes under potential control we have interfaced density functional theory with the hairy probe method in the weak coupling limit (Zauchner et al. Phys. Rev. B 2018, 97, 045116). Our implementation was systematically tested using simple parallel-plate capacitor models with pristine surfaces and a single layer of adsorbed water molecules. Remarkably, our code's efficiency is comparable with a standard DFT calculation. We reveal that local field effects at the electrical double layer induced by the change of applied potential can significantly affect the energies of chemical steps in heterogeneous electrocatalysis. Our results demonstrate the importance of an explicit modeling of the applied potential in a simulation and provide an efficient tool to control this critical parameter.

A computational model for a molecular chemical sensor.

In this study, we propose that a molecular junction with a sharp Negative Differential Resistance (NDR) current peak could improve the selectivity, thereby functioning as a potential molecular sensor for molecule recognition. Using DFT-NEGF simulations, we investigate the connection between molecule-molecule coupling, molecule-electrode coupling and the corresponding NDR peak shape. Based on this analysis we propose three design rules to control the sensitivity of a sensor and determine that one mechanism for NDR is for a localised molecular orbital involved in resonant tunneling to enter and leave the bias window. Our findings provide useful insight into the development of single molecule sensors for molecule recognition.

Fluorescence in quantum dynamics: Accurate spectra require post-mean-field approaches.

Real time modeling of fluorescence with vibronic resolution entails the representation of the light-matter interaction coupled to a quantum-mechanical description of the phonons and is therefore a challenging problem. In this work, taking advantage of the difference in timescales characterizing internal conversion and radiative relaxation-which allows us to decouple these two phenomena by sequentially modeling one after the other-we simulate the electron dynamics of fluorescence through a master equation derived from the Redfield formalism. Moreover, we explore the use of a recent semiclassical dissipative equation of motion [C. M. Bustamante et al., Phys. Rev. Lett. 126, 087401 (2021)], termed coherent electron electric-field dynamics (CEED), to describe the radiative stage. By comparing the results with those from the full quantum-electrodynamics treatment, we find that the semiclassical model does not reproduce the right amplitudes in the emission spectra when the radiative process involves the de-excitation to a manifold of closely lying states. We argue that this flaw is inherent to any mean-field approach and is the case with CEED. This effect is critical for the study of light-matter interaction, and this work is, to our knowledge, the first one to report this problem. We note that CEED reproduces the correct frequencies in agreement with quantum electrodynamics. This is a major asset of the semiclassical model, since the emission peak positions will be predicted correctly without any prior assumption about the nature of the molecular Hamiltonian. This is not so for the quantum electrodynamics approach, where access to the spectral information relies on knowledge of the Hamiltonian eigenvalues.

Role of electron localisation in H adsorption and hydride formation in the Mg basal plane under aqueous corrosion: a first-principles study.

Understanding hydrogen-metal interactions is important in various fields of surface science, including the aqueous corrosion of metals. The interaction between atomic H and a Mg surface is a key process for the formation of sub-surface Mg hydride, which may play an important role in Mg aqueous corrosion. In the present work, we performed first-principles Density Functional Theory (DFT) calculations to study the mechanisms for hydrogen adsorption and crystalline Mg hydride formation under aqueous conditions. The Electron Localisation Function (ELF) is found to be a promising indicator for predicting stable H adsorption in the Mg surface. It is found that H adsorption and hydride layer formation is dominated by high ELF adsorption sites. Our calculations suggest that the on-surface adsorption of atomic H, OH radicals and atomic O could enhance the electron localisation at specific sites in the sub-surface region, thus forming effective H traps locally. This is predicted to result in the formation of a thermodynamically stable sub-surface hydride layer, which is a potential precursor of the crucial hydride corrosion product of magnesium.

QM/MM optimization with quantum coupling: Host-guest interactions in a pentacene-doped p-terphenyl crystal.

In this work, we present a novel force-based scheme to perform hybrid quantum mechanics/molecular mechanics (QM/MM) computations. The proposed scheme becomes especially relevant for the simulation of host-guest molecular systems, where the description of the explicit electronic interactions between a guest molecule and a classically described host is of key importance. To illustrate its advantages, we utilize the presented scheme in the geometry optimization of a technologically important host-guest molecular system: a pentacene-doped p-terphenyl crystal, a core component of a room-temperature MASER device. We show that, as opposed to the simpler and widely used hybrid scheme ONIOM, our Quantum-Coupling QM/MM scheme was able to reproduce explicit interactions in the minimum energy configuration for the host-guest complex. We also show that, as a result of these explicit interactions, the host-guest complex exhibits an oriented net electric dipole moment that is responsible for red-shifting the energy of the first singlet-singlet electronic excitation of pentacene.

Dissipative Equation of Motion for Electromagnetic Radiation in Quantum Dynamics.

The dynamical description of the radiative decay of an electronically excited state in realistic many-particle systems is an unresolved challenge. In the present investigation electromagnetic radiation of the charge density is approximated as the power dissipated by a classical dipole, to cast the emission in closed form as a unitary single-electron theory. This results in a formalism of unprecedented efficiency, critical for ab initio modeling, which exhibits at the same time remarkable properties: it quantitatively predicts decay rates, natural broadening, and absorption intensities. Exquisitely accurate excitation lifetimes are obtained from time-dependent DFT simulations for C^{2+}, B^{+}, and Be, of 0.565, 0.831, and 1.97 ns, respectively, in accord with experimental values of 0.57±0.02, 0.86±0.07, and 1.77-2.5 ns. Hence, the present development expands the frontiers of quantum dynamics, bringing within reach first-principles simulations of a wealth of photophysical phenomena, from fluorescence to time-resolved spectroscopies.

A simple approximation to the electron-phonon interaction in population dynamics.

The modeling of coupled electron-ion dynamics including a quantum description of the nuclear degrees of freedom has remained a costly and technically difficult practice. The kinetic model for electron-phonon interaction provides an efficient approach to this problem, for systems evolving with low amplitude fluctuations, in a quasi-stationary state. In this work, we propose an extension of the kinetic model to include the effect of coherences, which are absent in the original approach. The new scheme, referred to as Liouville-von Neumann + Kinetic Equation (or LvN + KE), is implemented here in the context of a tight-binding Hamiltonian and employed to model the broadening, caused by the nuclear vibrations, of the electronic absorption bands of an atomic wire. The results, which show close agreement with the predictions given by Fermi's golden rule (FGR), serve as a validation of the methodology. Thereafter, the method is applied to the electron-phonon interaction in transport simulations, adopting to this end the driven Liouville-von Neumann equation to model open quantum boundaries. In this case, the LvN + KE model qualitatively captures the Joule heating effect and Ohm's law. It, however, exhibits numerical discrepancies with respect to the results based on FGR, attributable to the fact that the quasi-stationary state is defined taking into consideration the eigenstates of the closed system rather than those of the open boundary system. The simplicity and numerical efficiency of this approach and its ability to capture the essential physics of the electron-phonon coupling make it an attractive route to first-principles electron-ion dynamics.

Rectification and negative differential resistance via orbital level pinning.

A donor-acceptor system, 4-thiophenyl-azafulleroid (4TPA-C60), is investigated at the point of HOMO/LUMO resonance and beyond to understand how negative differential resistance (NDR) features may be observed in such systems. Our previous investigation showed that charge transfer between the occupied and unoccupied states at resonance hindered crossing of the HOMO and LUMO levels, thus preventing the formation of an NDR feature. In this work, it is shown that the negative differential resistance feature of 4TPA-C60 can be tailored based on the couplings at the metal/molecule interface. Ab initio calculations show that limited charge extraction from atomically sharp contacts results in a HOMO-LUMO pinning effect which delays the onset of the NDR feature. Subsequent unpinning of the states can only occur when additional charge extraction channels enter the bias window, highlighting an important role which non-frontier states play in charge transport. The proposed charge transfer mechanism is then exploited by introducing a fluorine atom into the C60 cage to tune the energies of the acceptor, and narrow the width of the current peak. These findings not only demonstrate the importance of the metal/molecule interface in the design of molecular electronic architectures but also serve to inform future design of molecular diodes and RTDs.

Efficient local-orbitals based method for ultrafast dynamics.

Computer simulations are invaluable for the study of ultrafast phenomena, as they allow us to directly access the electron dynamics. We present an efficient method for simulating the evolution of electrons in molecules under the influence of time-dependent electric fields, based on the Gaussian tight binding model. This model improves upon standard self-charge-consistent tight binding by the inclusion of polarizable orbitals and a self-consistent description of charge multipoles. Using the examples of bithiophene, terthiophene, and tetrathiophene, we show that this model produces electrostatic, electrodynamic, and explicitly time-dependent properties in strong agreement with density-functional theory, but at a small fraction of the cost.

HOMO-LUMO coupling: the fourth rule for highly effective molecular rectifiers.

Three rules for creating highly effective unimolecular rectifiers that utilize asymmetric anchoring groups have been proposed by Van Dyck and Ratner [Ratner et al., Nano Lett., 2015, 15, 1577-1584]. This study investigates their proposed rectification mechanism in a functionalised azafullerene system (4TPA-C60) and identifies a fourth rule. NEGF-DFT shows that 4TPA-C60 fulfills the three design rules and finds that a saturated bridge is not required to fulfil the third rule, contrary to previous belief. Instead a twisted-π bridge decouples the donor and acceptor states whilst still providing a high conductance pathway. The molecular junction has a calculated rectification ratio of 145 at a bias of ±1 V and the U-type rectification mechanism is driven by the pinning of the HOMO to the LUMO when the device is forward biased, but not when reverse biased. The switching behaviour is a result of a charge dipole forming at different interfaces for different bias directions. An additional design rule is thus proposed: charge transport should allow bias dependent coupling of filled to unfilled states. The findings in this work not only help in understanding charge transport in molecular rectifiers, but also have wider implications for the design of molecular resonant tunneling devices.

Gaussian polarizable-ion tight binding.

To interpret ultrafast dynamics experiments on large molecules, computer simulation is required due to the complex response to the laser field. We present a method capable of efficiently computing the static electronic response of large systems to external electric fields. This is achieved by extending the density-functional tight binding method to include larger basis sets and by multipole expansion of the charge density into electrostatically interacting Gaussian distributions. Polarizabilities for a range of hydrocarbon molecules are computed for a multipole expansion up to quadrupole order, giving excellent agreement with experimental values, with average errors similar to those from density functional theory, but at a small fraction of the cost. We apply the model in conjunction with the polarizable-point-dipoles model to estimate the internal fields in amorphous poly(3-hexylthiophene-2,5-diyl).

Correlation of Local Structure and Diffusion Pathways in the Modulated Anisotropic Oxide Ion Conductor CeNbO(4.25).

CeNbO4.25 is reported to exhibit fast oxygen ion diffusion at moderate temperatures, making this the prototype of a new class of ion conductor with applications in a range of energy generation and storage devices. To date, the mechanism by which this ion transport is achieved has remained obscure, in part due to the long-range commensurately modulated structural motif. Here we show that CeNbO4.25 forms with a unit cell ∼12 times larger than the stoichiometric tetragonal parent phase of CeNbO4 as a result of the helical ordering of Ce(3+) and Ce(4+) ions along z. Interstitial oxygen ion incorporation leads to a cooperative displacement of the surrounding oxygen species, creating interlayer "NbO6" connectivity by extending the oxygen coordination number to 7 and 8. Molecular dynamic simulations suggest that fast ion migration occurs predominantly within the xz plane. It is concluded that the oxide ion diffuses anisotropically, with the major migration mechanism being intralayer; however, when obstructed, oxygen can readily move to an adjacent layer along y via alternate lower energy barrier pathways.

Electron spin changes during general anesthesia in Drosophila.

We show that the general anesthetics xenon, sulfur hexafluoride, nitrous oxide, and chloroform cause rapid increases of different magnitude and time course in the electron spin content of Drosophila. With the exception of CHCl3, these changes are reversible. Anesthetic-resistant mutant strains of Drosophila exhibit a different pattern of spin responses to anesthetic. In two such mutants, the spin response to CHCl3 is absent. We propose that these spin changes are caused by perturbation of the electronic structure of proteins by general anesthetics. Using density functional theory, we show that general anesthetics perturb and extend the highest occupied molecular orbital of a nine-residue α-helix. The calculated perturbations are qualitatively in accord with the Meyer-Overton relationship and some of its exceptions. We conclude that there may be a connection between spin, electron currents in cells, and the functioning of the nervous system.

Nonlocal propagation and tunnelling of surface plasmons in metallic hourglass waveguides.

The nanofocusing performance of hourglass plasmonic waveguides is studied analytically and numerically. Nonlocal effects in the linearly tapered metal-air-metal stack that makes up the device are taken into account within a hydrodynamical approach. Using this hourglass waveguide as a model structure, we show that spatial dispersion drastically modifies the propagation of surface plasmons in metal voids, such as those generated between touching particles. Specifically, we investigate how nonlocal corrections limit the enormous field enhancements predicted by local electromagnetic treatments of geometric singularities. Finally, our results also indicate the emergence of nonlocality assisted tunnelling of plasmonic modes across hourglass contacts as thick as 0.5 nm.

Electron-energy loss study of nonlocal effects in connected plasmonic nanoprisms.

We investigate the emergence of nonlocal effects in plasmonic nanostructures through electron-energy loss spectroscopy. To theoretically describe the spatial dispersion in the metal permittivity, we develop a full three-dimensional nonlocal hydrodynamic solution of Maxwell's equations in frequency domain that implements the electron beam as a line current source. We use our numerical approach to perform an exhaustive analysis of the impact of nonlocality in the plasmonic response of single triangular prisms and connected bowtie dimers. Our results demonstrate the complexity of the interplay between nonlocal and geometric effects taking place in these structures. We show the different sensitivities to both effects of the various plasmonic modes supported by these systems. Finally, we present an experimental electron-energy loss study on gold nanoprisms connected by bridges as narrow as 1.6 nm. The comparison with our theoretical predictions enables us to reveal in a phenomenological fashion the enhancement of absorption damping that occurs in these atomistic junctions due to quantum confinement and grain boundary electron scattering.

The swipe card model of odorant recognition.

Just how we discriminate between the different odours we encounter is not completely understood yet. While obviously a matter involving biology, the core issue isa matter for physics: what microscopic interactions enable the receptors in our noses-small protein switches­to distinguish scent molecules? We survey what is and is not known about the physical processes that take place when we smell things, highlighting the difficulties in developing a full understanding of the mechanics of odorant recognition. The main current theories, discussed here, fall into two major groups. One class emphasises the scent molecule's shape, and is described informally as a "lock and key" mechanism. But there is another category, which we focus on and which we call "swipe card" theories:the molecular shape must be good enough, but the information that identifies the smell involves other factors. One clearly-defined "swipe card" mechanism that we discuss here is Turin's theory, in which inelastic electron tunnelling is used to discern olfactant vibration frequencies. This theory is explicitly quantal, since it requires the molecular vibrations to take in or give out energy only in discrete quanta. These ideas lead to obvious experimental tests and challenges. We describe the current theory in a form that takes into account molecular shape as well as olfactant vibrations. It emerges that this theory can explain many observations hard to reconcile in other ways. There are still some important gaps in a comprehensive physics-based description of the central steps in odorant recognition. We also discuss how far these ideas carry over to analogous processes involving other small biomolecules, like hormones, steroids and neurotransmitters. We conclude with a discussion of possible quantum behaviours in biology more generally, the case of olfaction being just one example. This paper is presented in honour of Prof. Marshall Stoneham who passed away unexpectedly during its writing.

Nonlocal effects in the nanofocusing performance of plasmonic tips.

The nanofocusing performance of plasmonic tips is studied analytically and numerically. The effects of electron-electron interactions in the dielectric response of the metal are taken into account through the implementation of a nonlocal, spatially dispersive, hydrodynamic permittivity. We demonstrate that spatial dispersion only slightly modifies the device parameters which maximize its field enhancement capabilities. The interplay between nonlocality, tip bluntness, and surface roughness is explored. We show that, although spatial dispersion reduces the field enhancement taking place at the structure apex, it also diminishes the impact that geometric imperfections have on its performance.

Plasmonic sinks for the selective removal of long-lived states.

The use of plasmonic nanostructures for the removal of unwanted long-lived states is investigated. We show that the total decay rate of such a state can be increased by up to 4 orders of magnitude, as compared to its intrinsic radiative decay rate, while leaving other neighboring optical transitions unaffected. For the specific case of molecular triplet excited states, we show that the use of a "plasmonic sink" has the potential to reduce photobleaching and ground-state depletion by at least 2 orders of magnitude. We consider, in addition, the impact of such structures on the performance of organic semiconductor lasers and show that, under realistic device conditions, plasmonic sinks have the capacity to increase the achievable laser repetition rate by a factor equal to the triplet decay rate enhancement. We conclude by studying the effect of exciton diffusion on the triplet density in the presence of metallic nanoparticles.