-
Rational Design of Molybdenum Transition-Metal subnanoclusters catalysts with Particle Swarm Optimization
Authors:
Yao Wei,
Alejandro Santana-Bonilla,
Lev Kantorovich
Abstract:
The development of novel sub-nanometer clusters (SNCs) catalysts with superior catalytic performance depends on the precise control of clusters' atomistic sizes, shapes, and accurate deposition onto surfaces. Recent advancements in manufacturing and characterization techniques have paved the way for the production and atomic resolution characterization of transition-metal SNCs catalysts, positioni…
▽ More
The development of novel sub-nanometer clusters (SNCs) catalysts with superior catalytic performance depends on the precise control of clusters' atomistic sizes, shapes, and accurate deposition onto surfaces. Recent advancements in manufacturing and characterization techniques have paved the way for the production and atomic resolution characterization of transition-metal SNCs catalysts, positioning them as a promising new class of materials for this application. Nevertheless, the intrinsic complexity of the adsorption process complicates the ability to achieve an atomistic understanding of the most relevant structure-reactivity relationships hampering the rational design of novel catalytic materials. In most cases, existing computational approaches rely on just a few structures to conclude clusters' reactivity thereby neglecting the complexity of the existing energy landscapes thus leading to insufficient sampling and, most likely, unreliable predictions. Moreover, modelling of the actual experimental procedure that is responsible for the deposition of SNCs on surfaces is often not done even though in some cases this procedure may enhance the significance of certain (e.g., metastable) adsorption geometries. This study proposes a novel approach that utilizes particle swarm optimization (PSO) method, in conjunction with ab-initio calculations, to predict the most relevant SNCs structures on a surface in beam experiments, and consequently their reactivity. To illustrate the main steps of our approach, we consider the deposition of Molybdenum SNC of 6 Mo atoms on a free-standing graphene surface, as well as their catalytic properties concerning the CO molecule dissociation reaction. This study demonstrates the feasibility of the PSO technique for studying catalyst transition-metal SNCs and establishes a reliable procedure for performing theoretical rational design predictions.
△ Less
Submitted 9 February, 2024;
originally announced February 2024.
-
A comprehensive exploration of structural and electronic properties of Molybdenum clusters
Authors:
Yao Wei,
Lev Kantorovich
Abstract:
Molybdenum clusters, characterised by their unique structure and intriguing catalytic properties, have gained significant attention in recent years. In several existing studies density functional theory (DFT) methods have been used to find the lowest energy Mo clusters and explore their electronic and magnetic structure. In all cases, with the exception of a single recent study, where a genetic al…
▽ More
Molybdenum clusters, characterised by their unique structure and intriguing catalytic properties, have gained significant attention in recent years. In several existing studies density functional theory (DFT) methods have been used to find the lowest energy Mo clusters and explore their electronic and magnetic structure. In all cases, with the exception of a single recent study, where a genetic algorithm was employed, initial geometries of the clusters, prior to geometry optimisation, were chosen using heuristic approaches based on symmetry considerations and known structures. DFT calculations were performed using different types of pseudopotentials, from hard to soft, and different types of basis sets. However, no comprehensive study has yet been done in which a DFT method with the best control on its precision would be complemented by a reliable global minimum search method to find the lowest energy Mo clusters. In this work, we employ a combination of a plane wave-based DFT method and \emph{ab initio} random structure searching (AIRSS) technique to find the lowest energy clusters of up to 10 Mo atoms. In each case, the search has been performed for clusters with different spin multiplicities, which enabled us to explore their magnetic structure. The results are compared for both hard and soft pseudopotentials stressing the importance of treating more electrons explicitly, in agreement with some of the previous studies. For most of the low-energy magnetic structures found, we investigate the distribution of their spin densities, and for all low energy clusters, we confirm their stability by calculating their phonon structure. Finally, free energies of the Mo clusters, within the quasi-harmonic approximation, are also calculated and discussed.
△ Less
Submitted 2 October, 2023;
originally announced October 2023.
-
Stochastic entropy production for continuous measurements of an open quantum system
Authors:
D. Matos,
L. Kantorovich,
I. J. Ford
Abstract:
We investigate the total stochastic entropy production of a two-level bosonic open quantum system under protocols of time dependent coupling to a harmonic environment. These processes are intended to represent the measurement of a system observable, and consequent selection of an eigenstate, whilst the system is also subjected to thermalising environmental noise. The entropy production depends on…
▽ More
We investigate the total stochastic entropy production of a two-level bosonic open quantum system under protocols of time dependent coupling to a harmonic environment. These processes are intended to represent the measurement of a system observable, and consequent selection of an eigenstate, whilst the system is also subjected to thermalising environmental noise. The entropy production depends on the evolution of the system variables and their probability density function, and is expressed through system and environmental contributions. The continuous stochastic dynamics of the open system is based on the Markovian approximation to the exact, noise-averaged stochastic Liouville-von Neumann equation, unravelled through the addition of stochastic environmental disturbance mimicking a measuring device. Under the thermalising influence of time independent coupling to the environment, the mean rate of entropy production vanishes asymptotically, indicating equilibrium. In contrast, a positive mean production of entropy as the system responds to time dependent coupling characterises the irreversibility of quantum measurement, and a comparison of its production for two coupling protocols, representing connection to and disconnection from the external measuring device, satisfies a detailed fluctuation theorem.
△ Less
Submitted 15 May, 2022;
originally announced May 2022.
-
Subsurface Carbon-Induced Local Charge of Copper for On-Surface Displacement Reaction
Authors:
Shaoshan Wang,
Pengcheng Ding,
Zhuo Li,
Cristina Mattioli,
Wenlong E,
Ye Sun,
André Gourdon,
Lev Kantorovich,
Flemming Besenbacher,
Xueming Yang,
Miao Yu
Abstract:
Transition metal carbides have sparked unprecedented enthusiasm as high-performance catalysts in recent years. Still, the catalytic properties of copper (Cu) carbide remain unexplored. By introducing subsurface carbon (C) to Cu(111), displacement reaction of proton in carboxyl acid group with single Cu atom is demonstrated at the atomic scale and room temperature. Its occurrence is attributed to t…
▽ More
Transition metal carbides have sparked unprecedented enthusiasm as high-performance catalysts in recent years. Still, the catalytic properties of copper (Cu) carbide remain unexplored. By introducing subsurface carbon (C) to Cu(111), displacement reaction of proton in carboxyl acid group with single Cu atom is demonstrated at the atomic scale and room temperature. Its occurrence is attributed to the C-doping induced local charge of surface Cu atoms (up to +0.30 e/atom), which accelerates the rate of on-surface deprotonation via reduction of the corresponding energy barrier, thus enabling the instant displacement of a proton with a Cu atom when the molecules land on the surface. Such well-defined and robust Cu$^{δ+}$ surface based on the subsurface C doping offers a novel catalytic platform for on-surface synthesis.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
On-Surface Decarboxylation Coupling Facilitated by Lock-to-Unlock Variation of Molecules upon the Reaction
Authors:
Shaoshan Wang,
Zhuo Li,
Pengcheng Ding,
Cristina Mattioli,
Wujun Huang,
Yang Wang,
André Gourdon,
Ye Sun,
Mingshu Chen,
Lev Kantorovich,
Xueming Yang,
Federico Rosei,
Miao Yu
Abstract:
On-surface synthesis (OSS) involving relatively high energy barriers remains challenging due to a typical dilemma: firm molecular anchor is required to prevent molecular desorption upon the reaction, whereas sufficient lateral mobility is crucial for subsequent coupling and assembly. By locking the molecular precursors on the substrate then unlocking them during the reaction, we present a strategy…
▽ More
On-surface synthesis (OSS) involving relatively high energy barriers remains challenging due to a typical dilemma: firm molecular anchor is required to prevent molecular desorption upon the reaction, whereas sufficient lateral mobility is crucial for subsequent coupling and assembly. By locking the molecular precursors on the substrate then unlocking them during the reaction, we present a strategy to address this challenge. High-yield synthesis based on well-defined decarboxylation, intermediate transition and hexamerization is demonstrated, resulting in an extended and ordered network exclusively composed of the newly-synthesized macrocyclic compound. Thanks to the steric hindrance of its maleimide group, we attain a preferential selection of the coupling. This work unlocks a promising path to enrich the reaction types and improve the coupling selectivity hence the structual homogeneity of the final product for OSS.
△ Less
Submitted 14 September, 2021;
originally announced September 2021.
-
Describing non-Hermitian dynamics using a Generalized Three-Time NEGF for a Partition-free Molecular Junction with Electron-Phonon Coupling
Authors:
Matthew A Lane,
Lev Kantorovich
Abstract:
In this paper we develop the Non-Equilibrium Green's Function (NEGF) formalism for a dissipative molecular junction that consists of a central molecular system with one-dimensional electronic transport coupled to a phonon environment and attached to multiple electronic leads. Our approach is partitionless - initial preparation of the system places the whole system in the correct canonical equilibr…
▽ More
In this paper we develop the Non-Equilibrium Green's Function (NEGF) formalism for a dissipative molecular junction that consists of a central molecular system with one-dimensional electronic transport coupled to a phonon environment and attached to multiple electronic leads. Our approach is partitionless - initial preparation of the system places the whole system in the correct canonical equilibrium state - and is valid for an external bias with arbitrary time dependence. Using path integrals as an intermediary tool, we apply a two-time Hubbard-Stratonovich transformation to the phonon influence functional with mixed real and imaginary times to obtain an exact expression for the electronic density matrix at the expense of introducing coloured Gaussian noises whose properties are rigorously derived from the environment action. This results in a unique stochastic Hamiltonian on each branch of the Konstantinov-Perel' contour (upper, lower, vertical) such that the time evolution operators in the Liouville equation no longer form a Hermitian conjugate pair, thus corresponding to non-Hermitian dynamics. To account for this we develop a generalized three-time NEGF which is sensitive to all branches of the contour, and relate it to the standard NEGF in the absence of phonons via a perturbative expansion of the noises. This approach is exact and fully general, describing the non-equilibrium driven dynamics from an initial thermal state while subject to inelastic scattering, and can be applied to non-Hermitian dynamics in general.
△ Less
Submitted 11 August, 2021;
originally announced August 2021.
-
Quantum interference and the time-dependent radiation of nanojunctions
Authors:
Michael Ridley,
Lev Kantorovich,
Robert van Leeuwen,
Riku Tuovinen
Abstract:
Using the recently developed time-dependent Landauer-Büttiker formalism and Jefimenko's retarded solutions to the Maxwell equations, we show how to compute the time-dependent electromagnetic field produced by the charge and current densities in nanojunctions out of equilibrium. We then apply this formalism to a benzene ring junction, and show that geometry-dependent quantum interference effects ca…
▽ More
Using the recently developed time-dependent Landauer-Büttiker formalism and Jefimenko's retarded solutions to the Maxwell equations, we show how to compute the time-dependent electromagnetic field produced by the charge and current densities in nanojunctions out of equilibrium. We then apply this formalism to a benzene ring junction, and show that geometry-dependent quantum interference effects can be used to control the magnetic field in the vicinity of the molecule. Then, treating the molecular junction as a quantum emitter, we demonstrate clear signatures of the local molecular geometry in the non-local radiated power.
△ Less
Submitted 1 February, 2021;
originally announced February 2021.
-
Efficient choice of coloured noises in stochastic dynamics of open quantum systems
Authors:
Daniel Matos,
Matthew A Lane,
Ian J Ford,
Lev Kantorovich
Abstract:
The Stochastic Liouville-von Neumann (SLN) equation describes the dynamics of an open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. The interaction with the environment is represented by complex coloured noises which drive the system, and whose correlation functions are set by the properties of the environment. We present a number of schemes capable of gene…
▽ More
The Stochastic Liouville-von Neumann (SLN) equation describes the dynamics of an open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. The interaction with the environment is represented by complex coloured noises which drive the system, and whose correlation functions are set by the properties of the environment. We present a number of schemes capable of generating coloured noises of this kind that are built on a noise amplitude reduction procedure [Imai et al, Chem. Phys. 446, 134 (2015)], including two analytically optimised schemes. In doing so, we pay close attention to the properties of the correlation functions in Fourier space, which we derive in full. For some schemes the method of Wiener filtering for deconvolutions leads to the realisation that weakening causality in one of the noise correlation functions improves numerical convergence considerably, allowing us to introduce a well controlled method for doing so. We compare the ability of these schemes, along with an alternative optimised scheme [Schmitz and Stockburger, Eur. Phys. J.: Spec. Top. 227, 1929 (2019)], to reduce the growth in the mean and variance of the trace of the reduced density matrix, and their ability to extend the region in which the dynamics is stable and well converged for a range of temperatures. By numerically optimising an additional noise scaling freedom, we identify the scheme which performs best for the parameters used, improving convergence by orders of magnitude and increasing the time accessible by simulation.
△ Less
Submitted 18 February, 2021; v1 submitted 2 June, 2020;
originally announced June 2020.
-
Exactly Thermalised Quantum Dynamics of the Spin-Boson Model coupled to a Dissipative Environment
Authors:
M. A. Lane,
D. Matos,
I. J. Ford,
L. Kantorovich
Abstract:
We present an application of the Extended Stochastic Liouville-von Neumann equations (ESLN) method introduced earlier [PRB 95, 125124 (2017); PRB 97, 224310 (2018)] which describes the dynamics of an exactly thermalised open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is t…
▽ More
We present an application of the Extended Stochastic Liouville-von Neumann equations (ESLN) method introduced earlier [PRB 95, 125124 (2017); PRB 97, 224310 (2018)] which describes the dynamics of an exactly thermalised open quantum system reduced density matrix coupled to a non-Markovian harmonic environment. Critically, the combined system of the open system fully coupled to its environment is thermalised at finite temperature using an imaginary time evolution procedure before the application of real time evolution. This initialises the combined system in the correct canonical equilibrium state rather than being initially decoupled. We apply our theory to the spin-boson Hamiltonian and develop a number of competing ESLN variants designed to reduce the numerical divergence of the trace of the open system density matrix. We find that a careful choice of the driving noises is essential for improving numerical stability. We also investigate the effect of applying higher order numerical schemes for solving stochastic differential equations, such as the Stratonovich-Heun scheme, and conclude that stochastic sampling dominates convergence with the improvement associated with the numerical scheme being less important for short times but required for late times. To verify the method and its numerical implementation, we consider evolution under a fixed Hamiltonian and show that the system either remains in, or approaches, the correct canonical equilibrium state at long times. Additionally, evolution of the open system under non-equilibrium Landau-Zener (LZ) driving is considered and the asymptotic convergence to the LZ limit was observed for vanishing system-environment coupling and temperature. When coupling and temperature are non-zero, initially thermalising the combined system at a finite time in the past was found to be a better approximation of the true LZ initial state than a pure state.
△ Less
Submitted 19 February, 2020; v1 submitted 18 February, 2020;
originally announced February 2020.
-
Non-adiabatic dynamics of electrons and atoms under non-equilibrium conditions
Authors:
L. Kantorovich
Abstract:
An approach to non-adiabatic dynamics of atoms in molecular and condensed matter systems under general non-equilibrium conditions is proposed. In this method interaction between nuclei and electrons is considered explicitly up to the second order in atomic displacements defined with respect to the mean atomic trajectory. This method enables one to consider movement of atoms beyond their simple vib…
▽ More
An approach to non-adiabatic dynamics of atoms in molecular and condensed matter systems under general non-equilibrium conditions is proposed. In this method interaction between nuclei and electrons is considered explicitly up to the second order in atomic displacements defined with respect to the mean atomic trajectory. This method enables one to consider movement of atoms beyond their simple vibrations. Both electrons and nuclei are treated fully quantum-mechanically using a combination of path integrals applied to nuclei and non-equilibrium Green's functions (NEGF) to elections. Our method is partition-less: initially, the entire system is coupled and assumed to be at thermal equilibrium. Then, the exact application of the Hubbard-Stratanovich transformation in mixed real and imaginary times enables us to obtain, without doing any additional approximations, an exact expression for the reduced density matrix for nuclei and hence an effective quantum Liouville equation for them, both containing Gaussian noises. It is shown that the time evolution of the expectation values for atomic positions is described by an infinite hierarchy of stochastic differential equations for atomic positions and momenta and their various fluctuations. The actual dynamics is obtained by sampling all stochastic trajectories. It is expected that applications of the method may include photo-induced chemical reactions (e.g. dissociation), electromigration, atomic manipulation in scanning tunneling microscopy, to name just a few.
△ Less
Submitted 22 March, 2018; v1 submitted 15 March, 2018;
originally announced March 2018.
-
Driving Spin-Boson Models From Equilibrium Using Exact Quantum Dynamics
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual a…
▽ More
We present an application of the Extended Stochastic Liouville Equation (ESLE) Phys. Rev. B 95, 125124, which gives an exact solution for the reduced density matrix of an open system surrounded by a harmonic heat bath. This method considers the extended system (the open system and the bath) being thermally equilibrated prior to the action of a time dependent perturbation, as opposed to the usual assumption that system and bath are initially partitioned. This is an exact technique capable of accounting for arbitrary parameter regimes of the model. Here we present our first numerical implementation of the method in the simplest case of a Caldeira-Leggett representation of the bath Hamiltonian, and apply it to a spin-Boson system driven from coupled equilibrium. We observe significant behaviours in both the transient dynamics and asymptotic states of the reduced density matrix not present in the usual approximation.
△ Less
Submitted 18 December, 2017;
originally announced December 2017.
-
Modelling a Bistable System Strongly Coupled to a Debye Bath: A Quasiclassical Approach Based on the Generalised Langevin Equation
Authors:
L. Stella,
H. Ness,
C. D. Lorenz,
L. Kantorovich
Abstract:
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised…
▽ More
Bistable systems present two degenerate metastable configurations separated by an energy barrier. Thermal or quantum fluctuations can promote the transition between the configurations at a rate which depends on the dynamical properties of the local environment (i.e., a thermal bath). In the case of classical systems, strong system-bath interaction has been successfully modelled by the Generalised Langevin Equation (GLE) formalism. Here we show that the efficient GLE algorithm introduced in Phys. Rev. B 89, 134303 (2014) can be extended to include some crucial aspects of the quantum fluctuations. In particular, the expected isotopic effect is observed along with the convergence of the quantum and classical transition rates in the strong coupling limit. Saturation of the transition rates at low temperature is also retrieved, in qualitative, yet not quantitative, agreement with the analytic predictions. The discrepancies in the tunnelling regime are due to an incorrect sampling close to the barrier top. The domain of applicability of the quasiclassical GLE is also discussed.
△ Less
Submitted 12 April, 2017;
originally announced April 2017.
-
Partition-free approach to open quantum systems in harmonic environments: an exact stochastic Liouville equation
Authors:
G. M. G. McCaul,
C. D. Lorenz,
L. Kantorovich
Abstract:
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation.…
▽ More
We present a partition-free approach to the evolution of density matrices for open quantum systems coupled to a harmonic environment. The influence functional formalism combined with a two-time Hubbard-Stratonovich transformation allows us to derive a set of exact differential equations for the reduced density matrix of an open system, termed the Extended Stochastic Liouville-von Neumann equation. Our approach generalises previous work based on Caldeira-Leggett models and a partitioned initial density matrix. This provides a simple, yet exact, closed-form description for the evolution of open systems from equilibriated initial conditions. The applicability of this model and the potential for numerical implementations are also discussed.
△ Less
Submitted 3 January, 2017; v1 submitted 16 December, 2016;
originally announced December 2016.
-
Nonequilibrium Generalised Langevin Equation for the calculation of heat transport properties in model 1D atomic chains coupled to two 3D thermal baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The…
▽ More
We use a Generalised Langevin Equation (GLE) scheme to study the thermal transport of low dimensional systems. In this approach, the central classical region is connected to two realistic thermal baths kept at two different temperatures [H. Ness et al., Phys. Rev. B {\bf 93}, 174303 (2016)]. We consider model Al systems, i.e. one-dimensional atomic chains connected to three-dimensional baths. The thermal transport properties are studied as a function of the chain length $N$ and the temperature difference $ΔT$ between the baths. We calculate the transport properties both in the linear response regime and in the non-linear regime. Two different laws are obtained for the linear conductance versus the length of the chains. For large temperatures ($T \gtrsim 500$ K) and temperature differences ($ΔT \gtrsim 500$ K), the chains, with $N > 18$ atoms, present a diffusive transport regime with the presence of a temperature gradient across the system. For lower temperatures($T \lesssim 500$ K) and temperature differences ($ΔT \lesssim 400$ K), a regime similar to the ballistic regime is observed. Such a ballistic-like regime is also obtained for shorter chains ($N \le 15 $). Our detailed analysis suggests that the behaviour at higher temperatures and temperature differences is mainly due to anharmonic effects within the long chains.
△ Less
Submitted 8 April, 2017; v1 submitted 3 December, 2016;
originally announced December 2016.
-
Partition-free theory of time-dependent current correlations in nanojunctions in response to an arbitrary time-dependent bias
Authors:
Michael Ridley,
Angus MacKinnon,
Lev Kantorovich
Abstract:
Working within the Nonequilibrium Green's Function (NEGF) formalism, a formula for the two-time current correlation function is derived for the case of transport through a nanojunction in response to an arbitrary time-dependent bias. The one-particle Hamiltonian and the Wide Band Limit Approximation (WBLA) are assumed, enabling us to extract all necessary Green's functions and self energies for th…
▽ More
Working within the Nonequilibrium Green's Function (NEGF) formalism, a formula for the two-time current correlation function is derived for the case of transport through a nanojunction in response to an arbitrary time-dependent bias. The one-particle Hamiltonian and the Wide Band Limit Approximation (WBLA) are assumed, enabling us to extract all necessary Green's functions and self energies for the system, extending the analytic work presented previously [Ridley et al. Phys. Rev. B (2015)]. We show that our new expression for the two-time correlation function generalises the Büttiker theory of shot and thermal noise on the current through a nanojunction to the time-dependent bias case including the transient regime following the switch-on. Transient terms in the correlation function arise from an initial state that does not assume (as is usually done) that the system is initially uncoupled, i.e. our approach is partition-free. We show that when the bias loses its time-dependence, the long time-limit of the current correlation function depends on the time difference only, as in this case an ideal steady state is reached. This enables derivation of known results for the single frequency power spectrum and for the zero frequency limit of this power spectrum. In addition, we present a technique which for the first time facilitates fast calculations of the transient quantum noise, valid for arbitrary temperature, time and voltage scales. We apply this to the quantum dot and molecular wire systems for both DC and AC biases, and find a novel signature of the traversal time for electrons crossing the wire in the time-dependent cross-lead current correlations.
△ Less
Submitted 16 November, 2016;
originally announced November 2016.
-
c-number Quantum Generalised Langevin Equation for an open system
Authors:
L. Kantorovich,
H. Ness,
L. Stella,
C. Lorenz
Abstract:
We derive a $c-$number Generalised Langevin Equation (GLE) describing the evolution of the expectation values $\left\langle x_{i}\right\rangle_{t}$ of the atomic position operators $x_{i}$ of an open system. The latter is coupled linearly to a harmonic bath kept at a fixed temperature. The equations of motion contain a non-Markovian friction term with the classical kernel {[}L. Kantorovich - PRB 7…
▽ More
We derive a $c-$number Generalised Langevin Equation (GLE) describing the evolution of the expectation values $\left\langle x_{i}\right\rangle_{t}$ of the atomic position operators $x_{i}$ of an open system. The latter is coupled linearly to a harmonic bath kept at a fixed temperature. The equations of motion contain a non-Markovian friction term with the classical kernel {[}L. Kantorovich - PRB 78, 094304 (2008){]} and a zero mean \emph{non-Gaussian} random force with correlation functions that depend on the initial preparation of the open system. We used a density operator formalism without assuming that initially the combined system was decoupled. The only approximation made in deriving quantum GLE consists in assuming that the Hamiltonian of the open system at time $t$ can be expanded up to the second order with respect to operators of atomic displacements $u_{i}=x_{i}-\left\langle x_{i}\right\rangle_{t}$ in the open system around their exact atomic positions $\left\langle x_{i}\right\rangle_{t}$ (the "harmonisation" approximation). The noise is introduced to ensure that sampling many quantum GLE trajectories yields exactly the average one. An explicit expression for the pair correlation function of the noise, consistent with the classical limit, is also proposed. Unlike the usually considered quantum operator GLE, the proposed $c-$number quantum GLE can be used in direct molecular dynamic simulations of open systems under general equilibrium or non-equilibrium conditions.
△ Less
Submitted 8 July, 2016;
originally announced July 2016.
-
Growth of Epitaxial Graphene: Theory and Experiment
Authors:
H. Tetlow,
J. Posthuma de Boer,
I. J. Ford,
D. D. Vvedensky,
J. Coraux,
L. Kantorovich
Abstract:
A detailed review of the literature for the last 5-10 years on epitaxial growth of graphene is presented. Both experimental and theoretical aspects related to growth on transition metals and on silicon carbide are thoroughly reviewed. Thermodynamic and kinetic aspects of growth on all these materials, where possible, are discussed. To make this text useful for a wider audience, a range of importan…
▽ More
A detailed review of the literature for the last 5-10 years on epitaxial growth of graphene is presented. Both experimental and theoretical aspects related to growth on transition metals and on silicon carbide are thoroughly reviewed. Thermodynamic and kinetic aspects of growth on all these materials, where possible, are discussed. To make this text useful for a wider audience, a range of important experimental techniques that have been used over the last decade to grow (e.g. CVD, TPG and segregation) and characterize (STM, LEEM, etc.) graphene are reviewed, and a critical survey of the most important theoretical techniques is given. Finally, we critically discuss various unsolved problems related to growth and its mechanism which we believe require proper attention in future research.
△ Less
Submitted 22 February, 2016;
originally announced February 2016.
-
Nonequilibrium processes from Generalised Langevin Equations: realistic nanoscale systems connected to two thermal baths
Authors:
H. Ness,
A. Genina,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the d…
▽ More
We extend the Generalised Langevin Equation (GLE) method [Phys. Rev. B 89, 134303 (2014)] to model a central classical region connected to two realistic thermal baths at two different temperatures. In such nonequilibrium conditions a heat flow is established, via the central system, in between the two baths. The GLE-2B (GLE two baths) scheme permits us to have a realistic description of both the dissipative central system and its surrounding baths. Following the original GLE approach, the extended Langevin dynamics scheme is modified to take into account two sets of auxiliary degrees of freedom corresponding to the mapping of the vibrational properties of each bath. These auxiliary variables are then used to solve the non-Markovian dissipative dynamics of the central region. The resulting algorithm is used to study a model of a short Al nanowire connected to two baths. The results of the simulations using the GLE-2B approach are compared to the results of other simulations that were carried out using standard thermostatting approaches (based on Markovian Langevin and Nose-Hoover thermostats). We concentrate on the steady state regime and study the establishment of a local temperature profile within the system. The conditions for obtaining a flat profile or a temperature gradient are examined in detail, in agreement with earlier studies. The results show that the GLE-2B approach is able to treat, within a single scheme, two widely different thermal transport regimes, i.e. ballistic systems, with no temperature gradient, and diffusive systems with a temperature gradient.
△ Less
Submitted 30 April, 2016; v1 submitted 11 December, 2015;
originally announced December 2015.
-
Calculation of the current response in a nanojunction for an arbitrary time-dependent bias: application to the molecular wire
Authors:
Michael Ridley,
Angus MacKinnon,
Lev Kantorovich
Abstract:
Recently [Phys. Rev. B 91, 125433 (2015)] we derived a general formula for the time-dependent quantum electron current through a molecular junction subject to an arbitrary time-dependent bias within the Wide Band Limit Approximation (WBLA) and assuming a single particle Hamiltonian. Here we present an efficient numerical scheme for calculating the current and particle number. Using the Padé expans…
▽ More
Recently [Phys. Rev. B 91, 125433 (2015)] we derived a general formula for the time-dependent quantum electron current through a molecular junction subject to an arbitrary time-dependent bias within the Wide Band Limit Approximation (WBLA) and assuming a single particle Hamiltonian. Here we present an efficient numerical scheme for calculating the current and particle number. Using the Padé expansion of the Fermi function, it is shown that all frequency integrals occurring in the general formula for the current can be removed analytically. Furthermore, when the bias in the reservoirs is assumed to be sinusoidal it is possible to manipulate the general formula into a form containing only summations over special functions. To illustrate the method, we consider electron transport through a one-dimensional molecular wire coupled to two leads subject to out-of-phase biases. We also investigate finite size effects in the current response and particle number that results from the switch-on of such a bias.
△ Less
Submitted 16 November, 2015;
originally announced November 2015.
-
Applications of the Generalised Langevin Equation: towards a realistic description of the baths
Authors:
H. Ness,
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding…
▽ More
The Generalised Langevin Equation (GLE) method, as developed in Ref. [Phys. Rev. B 89, 134303 (2014)], is used to calculate the dissipative dynamics of systems described at the atomic level. The GLE scheme goes beyond the commonly used bilinear coupling between the central system and the bath, and permits us to have a realistic description of both the dissipative central system and its surrounding bath. We show how to obtain the vibrational properties of a realistic bath and how to convey such properties into an extended Langevin dynamics by the use of the mapping of the bath vibrational properties onto a set of auxiliary variables. Our calculations for a model of a Lennard-Jones solid show that our GLE scheme provides a stable dynamics, with the dissipative/relaxation processes properly described. The total kinetic energy of the central system always thermalises toward the expected bath temperature, with appropriate fluctuation around the mean value. More importantly, we obtain a velocity distribution for the individual atoms in the central system which follows the expected canonical distribution at the corresponding temperature. This confirms that both our GLE scheme and our mapping procedure onto an extended Langevin dynamics provide the correct thermostat. We also examined the velocity autocorrelation functions and compare our results with more conventional Langevin dynamics.
△ Less
Submitted 18 December, 2014;
originally announced December 2014.
-
Current through a multi-lead junction caused by applied bias with arbitrary time-dependence
Authors:
Michael Ridley,
Angus MacKinnon,
Lev Kantorovich
Abstract:
We apply the Nonequilibrium Green's Function (NEGF) formalism to the problem of a multi-terminal nanojunction subject to an arbitrary time-dependent bias. In particular, we show that taking a generic one-particle system Hamiltonian within the wide band limit approximation (WBLA), it is possible to obtain a closed analytical expression for the current in each lead. Our formula reduces to the well-k…
▽ More
We apply the Nonequilibrium Green's Function (NEGF) formalism to the problem of a multi-terminal nanojunction subject to an arbitrary time-dependent bias. In particular, we show that taking a generic one-particle system Hamiltonian within the wide band limit approximation (WBLA), it is possible to obtain a closed analytical expression for the current in each lead. Our formula reduces to the well-known result of Jauho et. al. [doi:10.1103/PhysRevB.50.5528] in the limit where the switch-on time is taken to the remote past, and to the result of Tuovinen et. al. [doi:10.1088/1742-6596/427/1/012014] when the bias is maintained at a constant value after the switch-on. As we use a partition-free approach, our formula contains both the long-time current and transient effects due to the sudden switch-on of the bias. Numerical calculations performed for the simple case of a single-level quantum dot coupled to two leads are performed for a sinusoidally-varying bias. At certain frequencies of the driving bias, we observe `ringing' oscillations of the current, whose dependence on the dot level, level width, oscillation amplitude and temperature is also investigated.
△ Less
Submitted 12 December, 2014;
originally announced December 2014.
-
Pauli's Principle in Probe Microscopy
Authors:
Samuel Jarvis,
Adam Sweetman,
Lev Kantorovich,
Enda McGlynn,
Philip Moriarty
Abstract:
Exceptionally clear images of intramolecular structure can be attained in dynamic force microscopy through the combination of a passivated tip apex and operation in what has become known as the "Pauli exclusion regime" of the tip-sample interaction. We discuss, from an experimentalist's perspective, a number of aspects of the exclusion principle which underpin this ability to achieve submolecular…
▽ More
Exceptionally clear images of intramolecular structure can be attained in dynamic force microscopy through the combination of a passivated tip apex and operation in what has become known as the "Pauli exclusion regime" of the tip-sample interaction. We discuss, from an experimentalist's perspective, a number of aspects of the exclusion principle which underpin this ability to achieve submolecular resolution. Our particular focus is on the origins, history, and interpretation of Pauli's principle in the context of interatomic and intermolecular interactions.
△ Less
Submitted 5 August, 2014;
originally announced August 2014.
-
The Generalized Langevin Equation: An efficient approach to non-equilibrium molecular dynamics of open systems
Authors:
L. Stella,
C. D. Lorenz,
L. Kantorovich
Abstract:
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmo…
▽ More
The Generalized Langevin Equation (GLE) has been recently suggested to simulate the time evolution of classical solid and molecular systems when considering general non-equilibrium processes. In this approach, a part of the whole system (an open system), which interacts and exchanges energy with its dissipative environment, is studied. Because the GLE is derived by projecting out exactly the harmonic environment, the coupling to it is realistic, while the equations of motion are non-Markovian. Although the GLE formalism has already found promising applications, e.g., in nanotribology and as a powerful thermostat for equilibration in classical molecular dynamics simulations, efficient algorithms to solve the GLE for realistic memory kernels are highly non-trivial, especially if the memory kernels decay non-exponentially. This is due to the fact that one has to generate a colored noise and take account of the memory effects in a consistent manner. In this contribution, we present a simple, yet efficient, algorithm for solving the GLE for practical memory kernels and we demonstrate its capability for the exactly solvable case of a harmonic oscillator coupled to a Debye bath.
△ Less
Submitted 18 March, 2014; v1 submitted 17 December, 2013;
originally announced December 2013.
-
Non-equilibrium statistical mechanics of classical nuclei interacting with the quantum electron gas
Authors:
Yu Wang,
Lev Kantorovich
Abstract:
Kinetic equations governing time evolution of positions and momenta of atoms in extended systems are derived using quantum-classical ensembles within the Non-Equilibrium Statistical Operator Method (NESOM). Ions are treated classically, while their electrons quantum mechanically; however, the statistical operator is not factorised in any way and no simplifying assumptions are made concerning the…
▽ More
Kinetic equations governing time evolution of positions and momenta of atoms in extended systems are derived using quantum-classical ensembles within the Non-Equilibrium Statistical Operator Method (NESOM). Ions are treated classically, while their electrons quantum mechanically; however, the statistical operator is not factorised in any way and no simplifying assumptions are made concerning the electronic subsystem. Using this method, we derive kinetic equations of motion for the classical degrees of freedom (atoms) which account fully for the interaction and energy exchange with the quantum variables (electrons). Our equations, alongside the usual Newtonian-like terms normally associated with the Ehrenfest dynamics, contain additional terms, proportional to the atoms velocities, which can be associated with the electronic friction. Possible ways of calculating the friction forces which are shown to be given via complicated non-equilibrium correlation functions, are discussed. In particular, we demonstrate that the correlation functions are directly related to the thermodynamic Matsubara Green's functions, and this relationship allows for the diagrammatic methods to be used in treating electron-electron interaction perturbatively when calculating the correlation functions. This work also generalises previous attempts, mostly based on model systems, of introducing the electronic friction into Molecular Dynamics equations of atoms.
△ Less
Submitted 19 July, 2007;
originally announced July 2007.
-
Arrow diagram method based on overlapping electronic groups: corrections to the linked AD theorem
Authors:
Yu Wang,
Lev Kantorovich
Abstract:
Arrow diagram (AD) method (L. Kantorovich and B. Zapol, J. Chem. Phys. \textbf{96}, 8420 (1992); \emph{ibid}, 8427) provides a convenient means of systematic calculation of arbitrary matrix elements, $<Ψ|\hat{O}|Ψ>$, of symmetrical operators, $\hat{O}$, in quantum chemistry when the total system wavefunction $Ψ$ is represented as an antisymmetrised product of overlapping many-electron group func…
▽ More
Arrow diagram (AD) method (L. Kantorovich and B. Zapol, J. Chem. Phys. \textbf{96}, 8420 (1992); \emph{ibid}, 8427) provides a convenient means of systematic calculation of arbitrary matrix elements, $<Ψ|\hat{O}|Ψ>$, of symmetrical operators, $\hat{O}$, in quantum chemistry when the total system wavefunction $Ψ$ is represented as an antisymmetrised product of overlapping many-electron group functions, $Φ_{A}$, corresponding to each part (group) $A$ of the system: $Ψ=\hat{A}\prod_{A}Φ_{A}$. For extended (e.g. infinite) systems the calculation is somewhat difficult, however, as mean values of the operators require that each term of the diagram expansion is to be divided by the normalisation integral $S=<Ψ|Ψ>$, which is given by an AD expansion as well. A linked AD theorem suggested previously (L. Kantorovich, Int. J. Quant. Chem. \textbf{76}, 511 (2000)) to deal with this problem is reexamined in this paper using a simple Hartree-Fock problem of a one-dimensional ring of infinite size which is found to be analytically solvable. We find that corrections to the linked AD theorem are necessary in a general case of a finite overlap between different electronic groups. A general method of constructing these corrections in a form of a power series expansion with respect to overlap is suggested. It is illustrated on the ring model system.
△ Less
Submitted 16 January, 2005; v1 submitted 14 December, 2004;
originally announced December 2004.
-
Strongly localised molecular orbitals for $α$-quartz
Authors:
Oleh Danyliv,
Lev Kantorovich
Abstract:
A previously proposed computational procedure for constructing a set of nonorthogonal strongly localised one-electron molecular orbitals (O. Danyliv, L. Kantorovich - physics/0401107) is applied to a perfect $α$-quartz crystal characterised by an intermediate type of chemical bonding. The orbitals are constructed by applying various localisation methods to canonical Hartree-Fock orbitals calcula…
▽ More
A previously proposed computational procedure for constructing a set of nonorthogonal strongly localised one-electron molecular orbitals (O. Danyliv, L. Kantorovich - physics/0401107) is applied to a perfect $α$-quartz crystal characterised by an intermediate type of chemical bonding. The orbitals are constructed by applying various localisation methods to canonical Hartree-Fock orbitals calculated for a succession of finite molecular clusters of increased size with appropriate boundary conditions. The calculated orbitals span the same occupied Fock space as the canonical HF solutions, but have an advantage of reflecting the true chemical nature of the bonding in the system. The applicability of several localisation techniques as well as of a number of possible choices of localisation regions (structure elements) are discussed for this system in detail.
△ Less
Submitted 20 May, 2004;
originally announced May 2004.
-
Comparison of localization procedures for applications in crystal embedding
Authors:
Oleh Danyliv,
Lev Kantorovich
Abstract:
With the aim of future applications in quantum mechanical embedding in extended systems such as crystals, we suggest a simple and computationally efficient method which enables construction of a set of nonorthogonal highly localized one-electron orbitals for periodic nonmetallic crystals which reflect their chemical nature. The orbitals are also used to build up the Hartree-Fock (HF) electron de…
▽ More
With the aim of future applications in quantum mechanical embedding in extended systems such as crystals, we suggest a simple and computationally efficient method which enables construction of a set of nonorthogonal highly localized one-electron orbitals for periodic nonmetallic crystals which reflect their chemical nature. The orbitals are also used to build up the Hartree-Fock (HF) electron density of the entire crystals. The simplicity of the method stems from the fact that it does not require usage and/or modification of periodic electronic structure codes, and is instead based on the HF calculation of a sequence of finite clusters with subsequent application of a localization procedure to transform the HF canonical molecular orbitals. Two extreme cases of chemical bonding, ionic (MgO crystal) and covalent (Si crystal), are considered for which a number of known localization schemes are applied and compared. With some modifications our method can also be applied to nonperiodic nonmetallic systems as well.
△ Less
Submitted 21 January, 2004;
originally announced January 2004.
-
Calculation of electron density of periodic systems using non-orthogonal localised orbitals
Authors:
Lev Kantorovich,
Oleh Danyliv
Abstract:
Methods for calculating an electron density of a periodic crystal constructed using non-orthogonal localised orbitals are discussed. We demonstrate that an existing method based on the matrix expansion of the inverse of the overlap matrix into a power series can only be used when the orbitals are highly localised (e.g. ionic systems). In other cases including covalent crystals or those with an i…
▽ More
Methods for calculating an electron density of a periodic crystal constructed using non-orthogonal localised orbitals are discussed. We demonstrate that an existing method based on the matrix expansion of the inverse of the overlap matrix into a power series can only be used when the orbitals are highly localised (e.g. ionic systems). In other cases including covalent crystals or those with an intermediate type of chemical bonding this method may be either numerically inefficient or fail altogether. Instead, we suggest an exact and numerically efficient method which can be used for orbitals of practically arbitrary localisation. Theory is illustrated by numerical calculations on a model system.
△ Less
Submitted 21 January, 2004;
originally announced January 2004.
-
The Influence of Gradient Corrections on Bulk and Surface Properties of TiO2 and SnO2
Authors:
J. Goniakowski,
J. M. Holender,
L. N. Kantorovich,
M. J. Gillan,
J. A. White
Abstract:
First-principles calculations based on density functional theory and the pseudo\-potential method have been used to investigate the influence of gradient corrections to the standard LDA technique on the equilibrium structure and energetics of rutile TiO$_2$ and SnO$_2$ perfect crystals and their (110) surfaces. We find that gradient corrections increase the calculated lattice parameters by rough…
▽ More
First-principles calculations based on density functional theory and the pseudo\-potential method have been used to investigate the influence of gradient corrections to the standard LDA technique on the equilibrium structure and energetics of rutile TiO$_2$ and SnO$_2$ perfect crystals and their (110) surfaces. We find that gradient corrections increase the calculated lattice parameters by roughly 3~\%, as has been found for other types of material. Gradient corrections give only very minor changes to the equilibrium surface structure, but reduce the surface energies by about 30~\%.
△ Less
Submitted 24 August, 1995;
originally announced August 1995.
-
The Structure of the Stoichiometric and Reduced SnO2 (110) Surface
Authors:
I. Manassidis,
J. Goniakowski,
L. N. Kantorovich,
M. J. Gillan
Abstract:
First-principles calculations based on density functional theory (DFT) and the pseudopotential method have been used to study the stoichiometric and reduced SnO2 (110) surface. The ionic relaxations are found to be moderate for both the stoichiometric and reduced surfaces, and are very similar to those found in recent DFT-pseudopotential work on TiO2. Removal of neutral oxygen leaves two electro…
▽ More
First-principles calculations based on density functional theory (DFT) and the pseudopotential method have been used to study the stoichiometric and reduced SnO2 (110) surface. The ionic relaxations are found to be moderate for both the stoichiometric and reduced surfaces, and are very similar to those found in recent DFT-pseudopotential work on TiO2. Removal of neutral oxygen leaves two electrons per oxygen on the surface, which are distributed in channels passing through bridging oxygen sites. The associated electron density can be attributed to reduction of tin from Sn4+ to Sn2+, but only if the charge distribution on Sn2+ is recognized to be highly asymmetric. Reduction of the surface gives rise to a broad distribution of gap states, in qualitative agreement with spectroscopic measurements.
△ Less
Submitted 31 May, 1995;
originally announced May 1995.
-
The Energetics and Electronic Structure of Defective and Irregular Surfaces on MgO
Authors:
L. N. Kantorovich,
J. M. Holender,
M. J. Gillan
Abstract:
Ab initio calculations based on the density-functional pseudopotential approach have been used to study the fully relaxed structure, the electron distribution and the electronic density of states of (001) terraces, steps, corners and reverse corners, and of F-centers at these surface features on MgO. The calculations confirm earlier predictions of the relaxed structures of surface irregularities…
▽ More
Ab initio calculations based on the density-functional pseudopotential approach have been used to study the fully relaxed structure, the electron distribution and the electronic density of states of (001) terraces, steps, corners and reverse corners, and of F-centers at these surface features on MgO. The calculations confirm earlier predictions of the relaxed structures of surface irregularities based on simple interaction models. A substantial narrowing of the band-gap is found at the surface, which for terraces and steps is due to surface states at the bottom of the conduction band, but for the corner and reverse corner is also due to surface states at the top of the valence band. The F-center formation energy decreases steadily as the coordination of the oxygen site is reduced. The energy of the F-center level shows a tendency to approach the top of the valence band as the coordination of its site decreases.
△ Less
Submitted 1 May, 1995;
originally announced May 1995.