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.

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.

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.

Analog of Rabi oscillations in resonant electron-ion systems.

Quantum coherence between electron and ion dynamics, observed in organic semiconductors by means of ultrafast spectroscopy, is the object of recent theoretical and computational studies. To simulate this kind of quantum coherent dynamics, we have introduced in a previous article [L. Stella, M. Meister, A. J. Fisher, and A. P. Horsfield, J. Chem. Phys. 127, 214104 (2007)] an improved computational scheme based on Correlated Electron-Ion Dynamics (CEID). In this article, we provide a generalization of that scheme to model several ionic degrees of freedom and many-body electronic states. To illustrate the capability of this extended CEID, we study a model system which displays the electron-ion analog of the Rabi oscillations. Finally, we discuss convergence and scaling properties of the extended CEID along with its applicability to more realistic problems.

Nonconservative current-induced forces: A physical interpretation.

We give a physical interpretation of the recently demonstrated nonconservative nature of interatomic forces in current-carrying nanostructures. We start from the analytical expression for the curl of these forces, and evaluate it for a point defect in a current-carrying system. We obtain a general definition of the capacity of electrical current flow to exert a nonconservative force, and thus do net work around closed paths, by a formal noninvasive test procedure. Second, we show that the gain in atomic kinetic energy over time, generated by nonconservative current-induced forces, is equivalent to the uncompensated stimulated emission of directional phonons. This connection with electron-phonon interactions quantifies explicitly the intuitive notion that nonconservative forces work by angular momentum transfer.

Molecular conduction: do time-dependent simulations tell you more than the Landauer approach?

A dynamical method for simulating steady-state conduction in atomic and molecular wires is presented which is both computationally and conceptually simple. The method is tested by calculating the current-voltage spectrum of a simple diatomic molecular junction, for which the static Landauer approach produces multiple steady-state solutions. The dynamical method quantitatively reproduces the static results and provides information on the stability of the different solutions.

Correlated electron-ion dynamics with open boundaries: formalism.

We extend a new formalism, which allows correlated electron-ion dynamics to be applied to the problem of open boundary conditions. We implement this at the first moment level (allowing heating of ions by electrons) and observe the expected cooling in the classical part of the ionic kinetic energy and current-induced heating in the quantum contribution. The formalism for open boundaries should be easily extended to higher moments of the correlated electron-ion fluctuations.