Robust non-adiabatic molecular dynamics for metals and insulators

We present a new formulation of the correlated electron-ion dynamics (CEID) scheme, which systematically improves Ehrenfest dynamics by including quantum fluctuations around the mean-field atomic trajectories. We show that the method can simulate models of non-adiabatic electronic transitions, and test it against exact integration of the time-dependent Schroedinger equation. Unlike previous formulations of CEID, the accuracy of this scheme depends on a single tunable parameter which sets the level of atomic fluctuations included. The convergence to the exact dynamics by increasing the tunable parameter is demonstrated for a model two level system. This algorithm provides a smooth description of the non-adiabatic electronic transitions which satisfies the kinematic constraints (energy and momentum conservation) and preserves quantum coherence. The applicability of this algorithm to more complex atomic systems is discussed.Comment: 36 pages, 5 figures. Accepted for publication in Journal of Chemical Physic

Power dissipation in nanoscale conductors: classical, semi-classical and quantum dynamics

Modelling Joule heating is a difficult problem because of the need to introduce correct correlations between the motions of the ions and the electrons. In this paper we analyse three different models of current induced heating (a purely classical model, a fully quantum model and a hybrid model in which the electrons are treated quantum mechanically and the atoms are treated classically). We find that all three models allow for both heating and cooling processes in the presence of a current, and furthermore the purely classical and purely quantum models show remarkable agreement in the limit of high biases. However, the hybrid model in the Ehrenfest approximation tends to suppress heating. Analysis of the equations of motion reveals that this is a consequence of two things: the electrons are being treated as a continuous fluid and the atoms cannot undergo quantum fluctuations. A means for correcting this is suggested

Efficient simulations with electronic open boundaries

We present a reformulation of the Hairy Probe method for introducing electronic open boundaries that is appropriate for steady state calculations involving non-orthogonal atomic basis sets. As a check on the correctness of the method we investigate a perfect atomic wire of Cu atoms, and a perfect non-orthogonal chain of H atoms. For both atom chains we find that the conductance has a value of exactly one quantum unit, and that this is rather insensitive to the strength of coupling of the probes to the system, provided values of the coupling are of the same order as the mean inter-level spacing of the system without probes. For the Cu atom chain we find in addition that away from the regions with probes attached, the potential in the wire is uniform, while within them it follows a predicted exponential variation with position. We then apply the method to an initial investigation of the suitability of graphene as a contact material for molecular electronics. We perform calculations on a carbon nanoribbon to determine the correct coupling strength of the probes to the graphene, and obtain a conductance of about two quantum units corresponding to two bands crossing the Fermi surface. We then compute the current through a benzene molecule attached to two graphene contacts and find only a very weak current because of the disruption of the π-conjugation by the covalent bond between the benzene and the graphene. In all cases we find that very strong or weak probe couplings suppress the current

Adulteration and Poor Quality of Ginkgo biloba Supplements

Adulteration of Ginkgo products sold as unregistered supplements within the very large market of Ginkgo products (reputedly £650 million annually) through the post-extraction addition of cheaper (e.g. buckwheat derived) rutin is suspected to allow sub-standard products to appear satisfactory to third parties, e.g. secondary buyers along the value chain or any regulatory authorities. This study was therefore carried out to identify products that did not conform to their label specification and may have been actively adulterated to enable access to the global markets. 500 MHz Bruker NMR spectroscopy instrumentation combined with Topspin version 3.2 and a CAMAG HPTLC system (HPTLC Association for the analysis of Ginkgo biloba leaf) were used to generate NMR spectra (focusing on the 6–8 ppm region for analysis) and chromatograms, respectively. Out of the 35 samples of Ginkgo biloba analysed, 33 were found to contain elevated levels of rutin and/or quercetin, or low levels of Ginkgo metabolites when compared with the reference samples. Samples with disproportional levels of rutin or quercetin compared with other gingko metabolites are likely to be adulterated, either by accident or intentionally, and those samples with low or non-existent gingko metabolite content may have been produced using poor extraction techniques. Only two of the investigated samples were found to match with the High-Performance Thin-Layer Chromatography (HPTLC) fingerprint of the selected reference material. All others deviated significantly. One product contained a 5-hydroxytryptophan derivative, which is not a natural constituent of Ginkgo biloba. Overall, these examples either suggest a poor extraction technique or deliberate adulteration along the value chain. Investigating the ratio of different flavonoids e.g. quercetin and kaempferol using NMR spectroscopy and HPTLC will provide further evidence as to the degree and kind of adulteration of Gingko supplements. From a consumer perspective the equivalence in identity and overall quality of the products needs to be guaranteed for supplements too and not only for products produced according to a quality standard or pharmacopoeial monograph

Block bond-order potential as a convergent moments-based method

The theory of a novel bond-order potential, which is based on the block Lanczos algorithm, is presented within an orthogonal tight-binding representation. The block scheme handles automatically the very different character of sigma and pi bonds by introducing block elements, which produces rapid convergence of the energies and forces within insulators, semiconductors, metals, and molecules. The method gives the first convergent results for vacancies in semiconductors using a moments-based method with a low number of moments. Our use of the Lanczos basis simplifies the calculations of the band energy and forces, which allows the application of the method to the molecular dynamics simulations of large systems. As an illustration of this convergent O(N) method we apply the block bond-order potential to the large scale simulation of the deformation of a carbon nanotube.Comment: revtex, 43 pages, 11 figures, submitted to Phys. Rev.

Inelastic quantum transport: the self-consistent Born approximation and correlated electron-ion dynamics

A dynamical method for inelastic transport simulations in nanostructures is compared with a steady-state method based on non-equilibrium Green's functions. A simplified form of the dynamical method produces, in the steady state in the weak-coupling limit, effective self-energies analogous to those in the Born Approximation due to electron-phonon coupling. The two methods are then compared numerically on a resonant system consisting of a linear trimer weakly embedded between metal electrodes. This system exhibits enhanced heating at high biases and long phonon equilibration times. Despite the differences in their formulation, the static and dynamical methods capture local current-induced heating and inelastic corrections to the current with good agreement over a wide range of conditions, except in the limit of very high vibrational excitations, where differences begin to emerge.Comment: 12 pages, 7 figure

Could humans recognize odor by phonon assisted tunneling?

Our sense of smell relies on sensitive, selective atomic-scale processes that are initiated when a scent molecule meets specific receptors in the nose. However, the physical mechanisms of detection are not clear. While odorant shape and size are important, experiment indicates these are insufficient. One novel proposal suggests inelastic electron tunneling from a donor to an acceptor mediated by the odorant actuates a receptor, and provides critical discrimination. We test the physical viability of this mechanism using a simple but general model. Using values of key parameters in line with those for other biomolecular systems, we find the proposed mechanism is consistent both with the underlying physics and with observed features of smell, provided the receptor has certain general properties. This mechanism suggests a distinct paradigm for selective molecular interactions at receptors (the swipe card model): recognition and actuation involve size and shape, but also exploit other processes.Comment: 10 pages, 1 figur

A comprehensive computational model of sound transmission through the porcine lung

A comprehensive computational simulation model of sound transmission through the porcine lung is introduced and experimentally evaluated. This subject-specific model utilizes parenchymal and major airway geometry derived from x-ray CT images. The lung parenchyma is modeled as a poroviscoelastic material using Biot theory. A finite element (FE) mesh of the lung that includes airway detail is created and used in COMSOL FE software to simulate the vibroacoustic response of the lung to sound input at the trachea. The FE simulation model is validated by comparing simulation results to experimental measurements using scanning laser Doppler vibrometry on the surface of an excised, preserved lung. The FE model can also be used to calculate and visualize vibroacoustic pressure and motion inside the lung and its airways caused by the acoustic input. The effect of diffuse lung fibrosis and of a local tumor on the lung acoustic response is simulated and visualized using the FE model. In the future, this type of visualization can be compared and matched with experimentally obtained elastographic images to better quantify regional lung material properties to noninvasively diagnose and stage disease and response to treatment