Effects of Molecular Dynamics on Electrical Conductance of Single Molecular Junction in Aqueous Solution: First Principles Calculations

Total Page:16

File Type:pdf, Size:1020Kb

Effects of Molecular Dynamics on Electrical Conductance of Single Molecular Junction in Aqueous Solution: First Principles Calculations e-Journal of Surface Science and Nanotechnology 23 January 2010 e-J. Surf. Sci. Nanotech. Vol. 8 (2010) 38-43 Conference - ACSIN-10 - Effects of Molecular Dynamics on Electrical Conductance of Single Molecular Junction in Aqueous Solution: First Principles Calculations¤ Arihiro Tawara,y Tomofumi Tada, and Satoshi Watanabe Department of Materials Engineering, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-8656, Japan (Received 9 October 2009; Accepted 4 December 2009; Published 23 January 2010) The electronic transport in benzene-1,4-dithiolate molecule in aqueous solution sandwiched between gold elec- trodes have been investigated by the ab initio nonequilibrium Green’s function method combined with Car- Parrinello molecular dynamics. We have found that the C–S bond length shows clear negative correlation with the conductance both in aqueous solution and vacuum, whereas the Au–S and C–C bond lengths have little correlation. This originates from large local density of states around C–S bonds at the Fermi level. [DOI: 10.1380/ejssnt.2010.38] Keywords: Electrical transport; Benzene dithiolate (BDT); Water; Density functional calculations; Nonequilibrium Green’s function method; Molecular dynamics I. INTRODUCTION neither been understood enough nor explored in theoreti- cal fashion so far except for our previous calculation [32]. Recently, electronic transport properties of single In our previous study [32], we calculated a large number molecular junctions have been increasingly important and of conductances of a system consisting of a BDT molecule actively investigated in experimental and theoretical stud- and water molecules sandwiched between Au(100) sur- ies to deepen our understanding of nanoscale devices on faces in order to clarify bare effects of aqueous solution single molecule level. Benzene-1,4-dithiolate (BDT) be- on conductance at room temperature, using the ab ini- tween gold electrodes has been studied extensively as a tio NEGF-DFT method and Car-Parrinello molecular dy- benchmark molecule both in experimental [1–4] and the- namics (CPMD). Analyses of conductance histograms re- oretical [5–13] studies. In particular, theoretical inves- vealed the effects of the aqueous solution on the conduc- tigations on various effects such as adsorption configura- tance of the BDT: the peak of the conductance histogram tion [14–18], inelastic current [19–23] and limitation of the shifts downward by 0.01-0.02 G0 (5-10 %), which is at- density functional theory (DFT) [24–27] deepened our un- tributed to the electrostatic effects of the aqueous solu- derstanding on the transport properties of molecular junc- tion. We also found that subsidiary peaks in histograms tions. In experimental studies, conductance is mostly ob- emerge or disappear due to dynamical effects of the aque- served at room temperature and in solution, and the fluc- ous solution, which is correlated with the C–S stretching tuation of conductance is usually analyzed through his- mode of the BDT molecule. However, the analyses were tograms. Thus theoretical analyses based on conductance focused only on the C–S stretching, and the dynamical ef- histograms considering the fluctuations of the molecular fects induced by other modes have not been investigated junction are expected to provide a better understanding in the previous work [32]. In this paper, we assess the of the electronic transport through a single molecule. A correlations of conductance with Au–S and C–C lengths few such studies have already been reported for single as well as C–S length to understand the effects on con- molecular junctions in vacuum [28–30]. ductance in terms of the dynamical behavior of the BDT Several interesting studies on solution effects have also molecular junction further. been reported. In a mechanically controllable break junction experiment, Taniguchi et al. observed that the ¡ Au¡[Ni(dmit)2] ¡Au molecular junction shows almost II. COMPUTATIONAL METHOD the same conductance but different stability between the cases in a trichlorobenzene solution and in vacuum at Figure 1 shows the schematic of our computational room temperature [31]. Cao et al. studied the calculated model for the molecular junction adopted in our study, conductance histograms of a perylene tetracarboxylic di- where a BDT molecule is sandwiched between Au(100) imides molecule between gold electrodes in aqueous solu- surfaces. Fig. 1(a) and 1(b) show the unit cells of the sys- tion at different temperatures and showed a temperature tems with and without the aqueous solution (1.0 g/cm3, dependence of conductance due to water molecules [30]. 22 H2O), respectively. We abbreviate the systems in the However, much room for discussion is left and further in- aqueous solution and in vacuum as “Aq” and “Vac”. The vestigations are needed in the electronic transport proper- BDT is assumed to be adsorbed at the hollow site of ties affected by solution. Even for the simple BDT molec- Au(100) surface shown in Fig. 1(c) with the distance be- ular junction, the solution effects on conductance have tween electrodes of 9.8 A.˚ The computational model and procedures are identical to those given in our previous study [32]. Here, the computational procedures are de- scribed briefly. The procedure consists of three steps as ¤This paper was presented at 10th International Conference on follows. Atomically Controlled Surfaces, Interfaces and Nanostructures (ACSIN-10), Granada Conference Centre, Spain, 21-25 September, Step I: We perform CPMD calculations using the 2009. CPMD code [33, 34] to obtain a large number of configu- yCorresponding author: [email protected] rations in the Aq/Vac systems. Note that water and BDT ISSN 1348-0391 °c 2010 The Surface Science Society of Japan (http://www.sssj.org/ejssnt) 38 e-Journal of Surface Science and Nanotechnology Volume 8 (2010) (a) Aq (b) Vac 10 0.22 Conductance (G 9 z 0.20 8 0.18 7 0.16 6 0.14 5 (Å) 0.12 z 4 0.10 3 0.08 2 0 ) 0.06 1 0 0.04 0.0 0.5 1.0 1.5 Time (ps) (c) FIG. 2: Calculated time evolution of oxygen positions (left vertical axis, black lines) and conductance (right vertical axis, 1st layer red points) of BDT in the aqueous solution. Here, the z-axis on top bridge is perpendicular to the Au(100) surface, and z = 0 and 9.8 A˚ + + j 2nd layer z correspond to the topmost layers of the Au(100) surfaces. + hollow i (a) (b) 60 60 FIG. 1: (a) Snapshot of the system with the aqueous solu- Aq Vac tion (Aq), (b) the system without the aqueous solution (Vac) 50 50 and (c) the top view of a Au(100) surface together with the 40 40 adsorption sites (hollow, bridge, and on-top). Au, S, C, O, 30 30 Counts Counts and H atoms are depicted with gold, yellow, gray, red, and 20 20 white, respectively, in the snapshots. The z direction is per- 10 10 pendicular to the Au(100) surface. The white squares define 0 0 the scattering region in conductance calculations with ATK. 0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 Conductance (G ) Conductance (G ) 0 0 (c) (d) molecules are relaxed while Au atoms are fixed during the 60 60 Aq Vac dynamics. 50 50 Step II: We perform NEGF-DFT calculations using the 40 40 code Atomistix ToolKit (ATK) [35–38] to obtain zero-bias 30 30 Counts Counts conductances of BDT in Aq/Vac for hundreds of configu- 20 20 rations selected from the CPMD calculations in Step I. 10 10 Step III: We make conductance histograms using the 0 0 calculated conductance and analyze the Aq/Vac his- 0.0 0.1 0.2 0.3 0.0 0.1 0.2 0.3 Conductance (G ) Conductance (G ) tograms to characterize the effects of the aqueous solution 0 0 on conductance. Note that in the present study we use the same compu- FIG. 3: Calculated conductance histograms for the cases (a) tational data obtained in our previous work [32] to analyze in the aqueous solution (Aq) and (b) in vacuum (Vac). Fig- histogram in more detail. ure 3(a) is obtained from the calculated conductance shown as red points in Fig 2. The total number of samples is 400 in the Aq histogram (Fig. 3(a)) and 448 in the Vac one ( Fig. 3(b)). Figures 3(c) and 3(d) show conductance histograms obtained III. RESULTS AND DISCUSSION from 300 samples randomly selected from the histograms of Figs. 3(a) and 3(b), respectively. The bin sizes are 0.006 G0 in Figure 2 shows the time evolution of oxygen z-positions all the histograms. The light blue and black lines denote the of water molecules and conductance of BDT in water. Wa- respective Gaussian functions obtained by the fitting and the ter molecules confined in the nano-gapped space between sum of the Gaussians, respectively. the electrodes form three layers with each interlayer dis- tance of about 2.5 A,˚ and the uppermost and lowermost layers appear at about 2.5 A˚ apart from the gold surfaces. lines, we found two peaks in the histogram of Fig. 3(a) In Step II of our method, we picked up the CPMD con- and three in Fig. 3(b) as listed in Table I. The main figurations every 20 time-steps (1.94 fs) and calculated peaks are the 1st peak in the Aq histogram and the 3rd the zero-bias conductances for the selected configurations peak in the Vac, respectively. The downward shift of the using ATK to make conductance histograms. Figures 3(a) main peak in Aq is caused by the electrostatic effect of and (b) show the conductance histograms of the Aq and water molecules [32], i.e.
Recommended publications
  • Density Functional Theory
    Density Functional Approach Francesco Sottile Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) 22 October 2010 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Density Functional Theory 1. Any observable of a quantum system can be obtained from the density of the system alone. < O >= O[n] 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Hohenberg, P. and W. Kohn, 1964, Phys. Rev. 136, B864 Kohn, W. and L. Sham, 1965, Phys. Rev. 140, A1133 Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Nitrogen (7 electron) 1. Any observable of a quantum Ψ(r1; ::; r7) 21 coordinates system can be obtained from 10 entries/coordinate ) 1021 entries the density of the system alone. 8 bytes/entry ) 8 · 1021 bytes 4:7 × 109 bytes/DVD ) 2 × 1012 DVDs 2. The density of an interacting-particles system can be calculated as the density of an auxiliary system of non-interacting particles. Density Functional ... Why ? Density Functional ... Why ? Density Functional ... Why ? Basic ideas of DFT Importance of the density Example: atom of Oxygen (8 electron) 1. Any (ground-state) observable Ψ(r1; ::; r8) 24 coordinates of a quantum system can be 24 obtained from the density of the 10 entries/coordinate ) 10 entries 8 bytes/entry ) 8 · 1024 bytes system alone. 5 · 109 bytes/DVD ) 1015 DVDs 2.
    [Show full text]
  • Simulation Studies of the Mechanisms of Interaction Between Carbon Nanotubes and Amino Acids
    Simulation Studies of the Mechanisms of Interaction between Carbon Nanotubes and Amino Acids by George Bassem Botros Abadir B.Sc., Ain Shams University, Cairo, Egypt, 2000 M.Sc., Ain Shams University, Cairo, Egypt, 2005 A THESIS SUBMITTED IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY in The Faculty of Graduate Studies (Electrical and Computer Engineering) THE UNIVERSITY OF BRITISH COLUMBIA (Vancouver) August 2010 c George Bassem Botros Abadir 2010 Abstract In this thesis, molecular dynamics and ab initio density functional theory/non- equilibrium Green’s function simulations are used to study the interaction between carbon nanotubes and amino acids. Firstly, rules for the proper choice of the parameters used in these simulations are established. It is demonstrated how the improper choice of these parameters (particularly the basis set used in ab initio simulations) can lead to quantitatively and qualitatively erroneous conclusions regarding the bandgap of the nanotubes. It is then shown that the major forces responsible for amino-acid adsorp- tion on carbon nanotubes are van der Waals forces, and that hydrophobic interactions may accelerate the adsorption process, but are not necessary for it to occur. The mechanisms of interaction between carbon nanotubes and amino acids are elucidated. It is found that geometrical deformations do not play a major role in the sensing process, and that electrostatic inter- actions represent the major interaction mechanism between the tubes and amino acids. Fully metallic armchair tubes are found to be insensitive to various amino acids, while small-radius nanotubes are shown to be inade- quate for sensing in aqueous media, as their response to the motion of the atoms resulting from the immersion in water is comparable to that of an- alyte adsorption.
    [Show full text]
  • Newsletter 85
    Ψk Newsletter AB INITIO (FROM ELECTRONIC STRUCTURE) CALCULATION OF COMPLEX PROCESSES IN MATERIALS Number 85 February 2008 Editor: Z. (Dzidka) Szotek Sponsored by: UK’s CCP9 E-mail: [email protected] and ESF Psi-k Programme 1 Contents 1 Editorial 4 2 General News 5 2.1 Psi-k Portal and Psi-k Mailing List . ........ 5 3 Psi-k School/Workshop/Conference Announcements 6 3.1 THIRD Psi-k/NANOQUANTA SCHOOL& WORKSHOP . 6 3.2 Psi-kTraining-MarieCurieSummerSchool . ........ 11 4 General School/Workshop/Conference Announcements 12 4.1 2nd International Workshop on “Ab initio Description of Iron and Steel (ADIS2008): MagnetismandPhaseDiagrams” . 12 4.2 TheCAMDSummerSchool............................. 13 4.3 International Center for Materials Research (ICMR) Summer School on Multi- ferroicMaterialsandBeyond . 16 4.4 15th WIEN2khands-onworkshop . 18 4.5 Meeting on Molecular Dynamics for Non-Adiabatic Processes ........... 19 4.6 12th European Conference ”Physics of Magnetism” . ........... 20 4.7 SymposiumG,EMRSFallMeeting. 21 4.8 Summer School on Ab Initio Modelling in Solid State Chemistry - MSSC2008 . 23 4.9 Computational Molecular Science 2008 . ......... 25 5 General Job Announcements 26 6 Abstracts 55 7 New Book Announcements 73 8 SCIENTIFIC HIGHLIGHT OF THE MONTH: Local Self-Interaction Cor- rection of a Scattering Resonance: The Atom in Jellium Model 74 1 Introduction 74 2 Local Self-interaction Correction (LSIC) 76 2 3 LSDA solution of an atom in jellium 77 4 LSIC solution for an atom in jellium 78 5 Self-interaction correction of scattering resonances for Ce in jellium 79 6 Effective Medium Theory of Ce metal 82 7 Conclusions 84 3 1 Editorial As we are currently trying to change the format of the future Psi-k newsletters, the present newsletter in a way already reflects this transition period.
    [Show full text]
  • Quantum Chemical Calculations of NMR Parameters
    Quantum Chemical Calculations of NMR Parameters Tatyana Polenova University of Delaware Newark, DE Winter School on Biomolecular NMR January 20-25, 2008 Stowe, Vermont OUTLINE INTRODUCTION Relating NMR parameters to geometric and electronic structure Classical calculations of EFG tensors Molecular properties from quantum chemical calculations Quantum chemistry methods DENSITY FUNCTIONAL THEORY FOR CALCULATIONS OF NMR PARAMETERS Introduction to DFT Software Practical examples Tutorial RELATING NMR OBSERVABLES TO MOLECULAR STRUCTURE NMR Spectrum NMR Parameters Local geometry Chemical structure (reactivity) I. Calculation of experimental NMR parameters Find unique solution to CQ, Q, , , , , II. Theoretical prediction of fine structure constants from molecular geometry Classical electrostatic model (EFG)- only in simple ionic compounds Quantum mechanical calculations (Density Functional Theory) (EFG, CSA) ELECTRIC FIELD GRADIENT (EFG) TENSOR: POINT CHARGE MODEL EFG TENSOR IS DETERMINED BY THE COMBINED ELECTRONIC AND NUCLEAR WAVEFUNCTION, NO ANALYTICAL EXPRESSION IN THE GENERAL CASE THE SIMPLEST APPROXIMATION: CLASSICAL POINT CHARGE MODEL n Zie 4 V2,k = 3 Y2,k ()i,i i=1 di 5 ATOMS CONTRIBUTING TO THE EFG TENSOR ARE TREATED AS POINT CHARGES, THE RESULTING EFG TENSOR IS THE SUM WITH RESPECT TO ALL ATOMS VERY CRUDE MODEL, WORKS QUANTITATIVELY ONLY IN SIMPLEST IONIC SYSTEMS, BUT YIELDS QUALITATIVE TRENDS AND GENERAL UNDERSTANDING OF THE SYMMETRY AND MAGNITUDE OF THE EXPECTED TENSOR ELECTRIC FIELD GRADIENT (EFG) TENSOR: POINT CHARGE MODEL n Zie 4 V2,k = 3 Y2,k ()i,i i=1 di 5 Ze V = ; V = 0; V = 0 2,0 d 3 2,±1 2,±2 2Ze V = ; V = 0; V = 0 2,0 d 3 2,±1 2,±2 3 Ze V = ; V = 0; V = 0 2,0 2 d 3 2,±1 2,±2 V2,0 = 0; V2,±1 = 0; V2,±2 = 0 MOLECULAR PROPERTIES FROM QUANTUM CHEMICAL CALCULATIONS H = E See for example M.
    [Show full text]
  • UOW High Performance Computing Cluster User's Guide
    UOW High Performance Computing Cluster User’s Guide Information Management & Technology Services University of Wollongong ( Last updated on February 2, 2015) Contents 1. Overview6 1.1. Specification................................6 1.2. Access....................................7 1.3. File System.................................7 2. Quick Start9 2.1. Access the HPC Cluster...........................9 2.1.1. From within the UOW campus...................9 2.1.2. From the outside of UOW campus.................9 2.2. Work at the HPC Cluster.......................... 10 2.2.1. Being familiar with the environment................ 10 2.2.2. Setup the working space...................... 11 2.2.3. Initialize the computational task.................. 11 2.2.4. Submit your job and check the results............... 14 3. Software 17 3.1. Software Installation............................ 17 3.2. Software Environment........................... 17 3.3. Software List................................ 19 4. Queue System 22 4.1. Queue Structure............................... 22 4.1.1. Normal queue............................ 22 4.1.2. Special queues........................... 23 4.1.3. Schedule policy........................... 23 4.2. Job Management.............................. 23 4.2.1. PBS options............................. 23 4.2.2. Submit a batch job......................... 25 4.2.3. Check the job/queue status..................... 27 4.2.4. Submit an interactive job...................... 29 4.2.5. Submit workflow jobs....................... 30 4.2.6. Delete jobs............................. 31 5. Utilization Agreement 32 5.1. Policy.................................... 32 5.2. Acknowledgements............................. 32 5.3. Contact Information............................. 33 Appendices 35 A. Access the HPC cluster from Windows clients 35 A.1. Putty..................................... 35 A.2. Configure ‘Putty’ with UOW proxy.................... 36 A.3. SSH Secure Shell Client.......................... 38 B. Enable Linux GUI applications using Xming 41 B.1.
    [Show full text]
  • LOUISIANA SCIENTIST Vol. 5A No. 1
    Louisiana Scientist Bulletin of the Louisiana Academy of Sciences Volume 5A Number 1 (2015 Annual Meeting Abstracts) Published by THE LOUISIANA ACADEMY OF SCIENCES 05 July 2016 1 Louisiana Scientist Bulletin of the Louisiana Academy of Sciences IN THIS ISSUE Louisiana Academy of Sciences Abstracts of Presentations 2015 Annual Meeting Nicholls State University Thibodaux, Louisiana 14 March 2015 Division/Section Page Division of Agriculture, Forestry, and Wildlife . 4 Division of Biological Sciences . 7 Environmental Sciences Section . 7 Microbiology Section . 10 Molecular and Biomedical Biology Section . 11 Zoology Section . 16 Division of Physical Sciences . 23 Chemistry Section . 23 Computer Science Section . 26 Materials Science and Engineering Section . 31 Math and Statistics Section . 35 Physics Section . 36 Division of Science Education . .. 41 Higher Education Section . 41 Division of Sciences and Humanities . 43 Division of Social Sciences . 46 Acknowledgement . 50 2 The following abstracts of oral and poster presentations represent those received by the Abstract Editor. Authors’ affiliations are abbreviated as follows: CPMC Calcasieu Parish Mosquito Control FHS Franklinton High School GSU Grambling State University HSCA Harvard-Smithsonian Center for Astrophysics LSU-A Louisiana State University at Alexandria LSU-BR Louisiana State University, Baton Rouge LSU-E Louisiana State University, Eunice LTU Louisiana Tech University LU-NO Loyola University, New Orleans McSU McNeese State University NCTR National Centre for Toxicological
    [Show full text]
  • Virt&L-Comm.3.2012.1
    Virt&l-Comm.3.2012.1 A MODERN APPROACH TO AB INITIO COMPUTING IN CHEMISTRY, MOLECULAR AND MATERIALS SCIENCE AND TECHNOLOGIES ANTONIO LAGANA’, DEPARTMENT OF CHEMISTRY, UNIVERSITY OF PERUGIA, PERUGIA (IT)* ABSTRACT In this document we examine the present situation of Ab initio computing in Chemistry and Molecular and Materials Science and Technologies applications. To this end we give a short survey of the most popular quantum chemistry and quantum (as well as classical and semiclassical) molecular dynamics programs and packages. We then examine the need to move to higher complexity multiscale computational applications and the related need to adopt for them on the platform side cloud and grid computing. On this ground we examine also the need for reorganizing. The design of a possible roadmap to establishing a Chemistry Virtual Research Community is then sketched and some examples of Chemistry and Molecular and Materials Science and Technologies prototype applications exploiting the synergy between competences and distributed platforms are illustrated for these applications the middleware and work habits into cooperative schemes and virtual research communities (part of the first draft of this paper has been incorporated in the white paper issued by the Computational Chemistry Division of EUCHEMS in August 2012) INTRODUCTION The computational chemistry (CC) community is made of individuals (academics, affiliated to research institutions and operators of chemistry related companies) carrying out computational research in Chemistry, Molecular and Materials Science and Technology (CMMST). It is to a large extent registered into the CC division (DCC) of the European Chemistry and Molecular Science (EUCHEMS) Society and is connected to other chemistry related organizations operating in Chemical Engineering, Biochemistry, Chemometrics, Omics-sciences, Medicinal chemistry, Forensic chemistry, Food chemistry, etc.
    [Show full text]
  • Author's Personal Copy
    Author's personal copy Computer Physics Communications 183 (2012) 2272–2281 Contents lists available at SciVerse ScienceDirect Computer Physics Communications journal homepage: www.elsevier.com/locate/cpc Libxc: A library of exchange and correlation functionals for density functional theory✩ Miguel A.L. Marques a,b,∗, Micael J.T. Oliveira c, Tobias Burnus d a Université de Lyon, F-69000 Lyon, France b LPMCN, CNRS, UMR 5586, Université Lyon 1, F-69622 Villeurbanne, France c Center for Computational Physics, University of Coimbra, Rua Larga, 3004-516 Coimbra, Portugal d Peter Grünberg Institut and Institute for Advanced Simulation, Forschungszentrum Jülich, and Jülich Aachen Research Alliance, 52425 Jülich, Germany a r t i c l e i n f o a b s t r a c t Article history: The central quantity of density functional theory is the so-called exchange–correlation functional. This Received 8 March 2012 quantity encompasses all non-trivial many-body effects of the ground-state and has to be approximated Received in revised form in any practical application of the theory. For the past 50 years, hundreds of such approximations have 7 May 2012 appeared, with many successfully persisting in the electronic structure community and literature. Here, Accepted 8 May 2012 we present a library that contains routines to evaluate many of these functionals (around 180) and their Available online 19 May 2012 derivatives. Keywords: Program summary Density functional theory Program title: LIBXC Density functionals Catalogue identifier: AEMU_v1_0 Local density approximation Generalized gradient approximation Program summary URL: http://cpc.cs.qub.ac.uk/summaries/AEMU_v1_0.html Hybrid functionals Program obtainable from: CPC Program Library, Queen’s University, Belfast, N.
    [Show full text]
  • Jacob Davis Gayles
    Jacob Davis Gayles Graduate Research Assistant CONTACT Johannes Gutenberg-Universität Mainz Work: +49 6131 39 27424 INFORMATION INSPIRE group Mobile: +49 173 622 74 14 Staudingerweg 7 E-mail: [email protected] Mainz, DE 55128 WWW: INSPIRE/jacob-gayles PhD Candidate Texas A&M University Department of Physics & Astronomy College Station, TX 77843-4242 RESEARCH Spin-Orbit Effects, First Principle Calculations, Spin-Transport, Non-Trivial Magnetism, INTERESTS Hall Effects, Berry Phase Physics, Skyrmions, Magnetic Exchange Interactions, Weyl Semimetals, Transitions in Magnetic Systems, Kubo Formalism, Boltzmann Tranport, Non-equilibrium Green Functions, Tight-Binding models ACADEMIC Graduate Research Assistant, Johannes Gutenberg Universität Mainz Jan 2015 to present APPOINTMENTS INSPIRE GROUP • Supervisor: Jairo Sinova − Spintronics − Spin orbit torques using Scattering Approach − Anomalous and Topological Hall effect in non-trivial magnetic systems − First Principle calculations of the Dzyaloshinskii-Moriya Interaction − First Principle calculations of B20 compounds, Half-Heuslers, Graphene, and Stanene Graduate Research Assistant, Forschungszentrum Jülich Summer 2013 Prof. Mokoursov Topoligical Group − Helicity and size dependence of skyrmions in B20 compounds Graduate Research Assistant, Institute of Physics, Academy of Sciences Czech Republic Summer 2012 Prof. Tomas Jungwirth − Anomalous Hall effect in 2D multiband systems with strong disorder Graduate Student Researcher, Texas A&M University August 2011 to present Department of Physics & Astronomy • Supervisor: Jairo Sinova • Supervisor: Artem Abanov Undergraduate Student Researcher, California State University Northridge January 2007 to August 2011 Department of Physics & Astronomy • National Science Foundation PREM Center • W.M. Keck Computational Materials Theory Center • Supervisor (PI): Professor Nicholas Kioussis • co-Supervisor (co-PI): Professor Donna Sheng − Transferable Tight-Binding parameters for CdTe (2011) − Transport Properties of DNA Base Pairs-Prof.
    [Show full text]
  • Nemoviz: a Visual Interactive System for Atomistic Simulations Design Daniel Mejia* , Tillmann Kubis and Gerhard Klimeck
    Mejia et al. Visualization in Engineering (2018) 6:6 https://doi.org/10.1186/s40327-018-0067-4 CASESTUDY Open Access NemoViz: a visual interactive system for atomistic simulations design Daniel Mejia* , Tillmann Kubis and Gerhard Klimeck Abstract Introduction: Given that the nanoscale regime has been reached, atomistic simulations are being used as predictive tools on a nanoscopic scale in nanoelectronics, materials science, and computational fluid dynamics. These simulations are generally supported by complex and specific simulation engines that require a deep knowledge of the engine as well as the field. The inputs required by these engines are usually described in long-text files, and their composition is error prone. These files contain details about the corresponding materials, geometries, algorithms, initial values, etc. The introduction of visualizations to the simulation creation process could alleviate some common user issues with simulation creation, in particular regarding the input-text files. In other fields, visual analytics have been used to facilitate user interaction with complex data. This work features a case study on the creation of a nanoelectronic device simulation and provides evidence of how the use of visual analytics reduces the cognitive complexity of defining an atomistic simulation. Case description: NEMO5 is a tool designed to simulate the electronic properties of nanoelectronic devices on an atomistic level. In this work we introduce NemoViz, an interactive visualization tool that enables users to define the simulation inputs required by NEMO5. Regular NEMO5 users were exposed to NemoViz, and users’ effectiveness and efficiency was measured while debugging the input for a simulation. Discussion and Evaluation: The results of this work show that NemoViz reduced the time that users spent defining the inputs of a NEMO5 simulation.
    [Show full text]
  • Quantum Many-Particle Electron Transport in Time-Dependent Systems with Bohmian Trajectories by Alfonso Alarcón Pardo B.S. in P
    Quantum many-particle electron transport in time-dependent systems with Bohmian trajectories by Alfonso Alarc´onPardo B.S. in Physics, 2007 M.S. in Nanoelectronics 2009 Submitted to the Departament d'Enginyeria Electr`onica in partial fulfillment of the requirements for the degree of PhD in Electronic Engineering at the Universitat Aut`onomade Barcelona Thesis supervisor: Xavier Oriols Pladevall Date of defense: April 2011 i Summary It is known that at nanoscale regime we must deal with the many-particle problem in order to study electronic devices. In this scenario, the time-dependent many-particle Schr¨odinger equation is only directly solvable for very few degrees of freedom. However, there are many electrons (degrees of freedom) in any electron device. In this sense, many-particle quantum electron formalisms (such as time-dependent Density Functional Theory, Green's functions techniques or Quantum Monte Carlo techniques) have been developed in the literature to provide reasonable approximations to model many-particle electron transport. An alterna- tive proposal has been developed by Dr. Oriols to decompose the N-particle Schr¨odinger equation into a N-single particle Schr¨odingerequation using Bohmian trajectories. Based on this proposal a general, versatile and time-dependent 3D electron transport simulator for na- noelectronic devices, named BITLLES (Bohmian Interacting Transport for non-equiLibrium eLEctronic Structures) is presented. The novelty of the BITLLES simulator is based on two points. First, it presents a many- particle quantum electron transport model taking into account explicitly the Coulomb and exchange correlations among electrons using Bohmian trajectories. Second, it provides a full information of the all current distribution moments (i.e.
    [Show full text]
  • Photoinduced Charge Transfer in Quasi-One-Dimensional Polymers in Two-Photon Absorption
    RSC Advances View Article Online PAPER View Journal | View Issue Photoinduced charge transfer in quasi-one- dimensional polymers in two-photon absorption† Cite this: RSC Adv., 2020, 10,33288 Pen-ji Yan,‡a Xijiao Mu,‡b Jun Dong *c and Mengtao Sun *b In this work, we theoretically investigate the structure and the transition characteristics of one- (OPA) and two-photon absorption (TPA) spectra of different length neutral and charged thiophene polymers. The effects and regulation of different charges on photoinduced charge transfer are discovered and their physical mechanisms are explained. We find that both the OPA and TPA spectra undergo a sizeable redshift after the charge is injected into the polymer, and the redshift after the positive charge injection is excellent. The alternating charge transfer that occurs in a two-photon transition of a charged system is derived from the alternating distribution of charge (dipole moment) in the dynamics of the system. To Received 24th July 2020 study the gradual behavior of infinite polymers, we also theoretically calculated the optical properties Accepted 2nd September 2020 and electronic structures of infinitely long thiophene polymers under different electrical charge DOI: 10.1039/d0ra06436j injections by a one-dimensional periodic model. The redshift of the OPA and TPA spectra is found to be Creative Commons Attribution 3.0 Unported Licence. rsc.li/rsc-advances due to orbital energy level movement. 1. Introduction charge injection, so its optical properties also change greatly. This change cannot be ignored when designing and In recent years, with the development of organic electronics,1,2 manufacturing devices.
    [Show full text]