SUPPLEMENTARY INFORMATION
|
|
- Myles Elliott
- 5 years ago
- Views:
Transcription
1 In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION DOI: /NCHEM.2824 Supramolecular heterostructures formed by sequential epitaxial deposition of two-dimensional hydrogen-bonded arrays Vladimir V. Korolkov 1, Matteo Baldoni 2, Kenji Watanabe 3, Takashi Taniguchi 3, Elena Besley 2 and Peter H. Beton 1 * 1 School of Physics & Astronomy, University of Nottingham, Nottingham NG7 2RD, UK 2 School of Chemistry, University of Nottingham, Nottingham NG7 2RD, UK 3 National Institute for Materials Science, 1-1 Namiki, Tsukuba, Ibaraki , Japan *elena.besley@nottingham.ac.uk; peter.beton@nottingham.ac.uk Supplementary Information 1. Additional AFM Images 2 2. Simulations : Methods and Additional Results 4 NATURE CHEMISTRY Macmillan Publishers Limited, part of Springer Nature. All rights reserved.
2 1. Additional AFM images hbn 1.2nm 2.0nm Supplementary Figure 1 Contact mode image of CA.M network and underlying hbn lattice showing a misalignment of principal axes by 3 ± 1 o. Both image were acquired in contact mode from the same area of sample. The scan of hbn lattice was taken at a slightly elevated setpoint (typical value 2-3 nn) allowing a punch-through of the CA.M layer to reveal the hbn lattice. 200 Z[nm] µm X[µm] Supplementary Figure 2 An AFM image showing multilayers of CA.M network on a 102nm thick hbn flake. The corresponding cross-section provides height information for the thicknesses of these layers which ranges between 16 nm to 50 nm. Multilayers were grown from 240 μm aqueous solution with a depositon time of 10 minute. 2
3 a 400nm b 400nm Supplementary Figure 3 AFM scans showing surface coverage of TMA on CA.M. Samples were prepared via 100 sec (a) and 300 sec (b-d) immersion of a CA.M/BN sample into a TMA solution in EtOH (0.04mg/ml). Scans were acquired in AC mode. 2.4nm Supplementary Figure 4 AFM showing lattice resolution of a TMA overlayer and CA.M underlayer. The dotted lines highlight the misalignment of pores of TMA and pores CA.M; note that the pores are not coincident. 3
4 2. Simulations: methods and additional results 2.1 General Description of the Methods Molecular mechanics (MM) calculations have been performed using the LAMMPS simulation package 1. The OPLS potential 2 has been employed, which was previously used to model the self-assembly of CA molecules 3 and molecular adsorption on 2D materials 4. The OPLS parameters for TPA and TMA can be found in the general OPLS database, however for CA and M molecules these parameters are not included. To describe CA and M molecules, the parameters for the Lennard-Jones (LJ) potential and bonding terms (bonds, angles dihedrals and impropers) were inferred from uracil and 1,3,5-triazine/aniline molecules, respectively. The atomic charges were calculated using the restrained electrostatic potential (RESP) scheme 5 based on density functional theory (DFT) as implemented in the CP2K program package 6. Identical types of atoms were constrained to have the same charge. DFT calculations were performed at the gradient-corrected level by applying the BLYP 7,8 exchange-correlation functional within the GPW approach 9,10. Electronic states were expanded by a double-ζ plus polarization basis set, DZVP, 11 with norm-conserving pseudopotentials for the description of core levels and a plane-wave representation of the charge density with a cut-off of 500 Ry. The Martyna-Tuckerman Poisson solver 15 was used in combination with a relatively large simulation box (up to 3.0 nm in the main dimension) to avoid computational artefacts. Long-range electrostatic interactions were treated by the particle-particle particle-mesh method 16 with the relative error in electrostatic forces of A cut-off distance of 1.4 nm was used for both Coulomb and van der Waals (vdw) interactions. Energy minimisation, using the conjugate gradient algorithm with tight convergence criteria, was used to locate the minimum potential energy structures. During energy minimisation, the LJ potential was shifted at the cut-off value to 0. Molecular dynamics (MD) simulations have been performed with a time step of 1 fs. The Nosé Hoover thermostat was used for simulations in the canonical NVT ensemble with a time constant of 1.0 ps. Periodic boundary conditions (PBC) in 3 dimensions were used for the simulation of the 2D systems, by inserting a vacuum region of about 4.0 nm along the z-direction. The structures and potential were built using VMD 20 and Moltemplate 21. 4
5 Supplementary Figure 5 (a) CA.M (in red and blue respectively) on hbn (in white) model shown as a snapshot from an MD at 300 K;(b) 2D radial distribution function and coordination number of CA+M (bin size 0.01 nm) 5
6 2.2 Interaction between CA.M and hbn The hexagonal BN layer (hbn) was treated as rigid with atoms fixed at their crystallographic positions. The LJ parameters were taken from Ref. 22. The values for atomic charges of B and N, ranging from 0 to ±0.6 e, were initially tested leading to the same qualitative results; hence B and N atoms were kept neutral in all simulations. A square hbn model with dimensions of 61 nm x 61 nm with PBC was used to represent an infinite hbn plane. A layer of alternating CA and M molecules with a hexagonal motif was placed on top of hbn (Supplementary Figure 5a). The CA.M layer was made a bit smaller in dimensions than hbn in order to remove any bias due to the registry mismatch at the edges of the simulation cell. In MD simulations, the CA.M layer has been equilibrated in an NVT ensemble at 300 K for 5 ns and then propagated for a further 2 ns in an NVE ensemble. The effect of relative orientation of the CA.M and hbn lattices has been extensively studied by performing a lateral slide analysis of the potential energy surface (PES) (followed by relaxation of the CA.M layer) for different orientations over a range of 30 o, i.e. ranging from armchair to zigzag direction in the hbn lattice. The minimum value for the interaction energy between the layers was found to be around ±1-3, in good agreement with the experimental results, although the difference in the interaction energy remains below kt/molecule at all angles. For the CA.M layer, the optimised value of the average lattice vector at 0K is found to be nm and the angle between lattice vectors is 120. The stability of the CA.M layer has also been probed by MD simulations at 300K, which confirmed that the adsorbed molecular layer remains stable at room temperature. The lattice constant of CA.M layer at 300K has been estimated through the 2D radial distribution function between equivalent atoms in molecules, which peaks at nm (Supplementary Figure 5b), in excellent agreement with experiment. 2.3 Interaction between TMA, TPA and CA.M/hBN The interaction between TPA and TMA arrangements and CA.M has been investigated by the lateral slide analysis of PES. A CA.M hexagonal arrangement consisting of 8024 molecules has been prepared with a periodicity of 9.75 Å (the lattice constant predicted by 6
7 MD simulations at 300K) and minimised in vacuum keeping the central ring of each molecule fixed at its original position. PBCs applied to the simulation box were chosen to be an exact multiple of the CA.M lattice vector in the xy-directions, namely nm x nm. A row of TPA containing 40 molecules (Supplementary Figure 6a) was optimised in vacuum producing an intermolecular distance of nm. The energy of a single hydrogen bond, defined as ( ) ( ), where n is the number of TPA molecules in a row, is ( ) calculated to be 0.27 ev. The optimised TPA row has been placed onto the CA.M layer, separated by 0.35 nm. The TPA row is then moved as a rigid body along x and y with a step size of 0.01 nm. At each step, the adsorption energy was calculated as ) ( ) (. ). This gives a pure adsorption contribution to the total energy excluding the energy required to form hydrogen bonds. The adsorption energy as a function of the displacement along the XY plane for different angles between the principal axis of the row and the CA.M lattice vector has been calculated (see Supplementary Fig 6). PES profiles for the rotation angles of 0, 10, 20, and 30 are presented in Supplementary Figures 6b-e. The lowest energy corresponds to orientations where a TPA row is aligned with the CA.M lattice vector, i.e. at the rotation angles of 0, 60 and 120 ; the adsorption process at these angles is favoured by about 0.09 ev per molecule (see Supplementary Figure 10). After full relaxation of TPA molecules on the CA.M layer, the adsorption energy at 0 is ev/molecule. The energy decomposition, shown in Supplementary Figure 6b, indicates that although the van der Waals interactions account for the largest contribution to the total energy of adsorption, the exact position of the minimum is strongly influenced by the Coulombic energy. Except for preferred orientations, the PES is fairly flat representing an average of inregistry and out-of-registry positions of the interacting layers. The low energy regions near the minima extend to 0.2 to 0.4 nm in each direction allowing for some flexibility in the positioning of TPA molecules on CA.M layer, which is important in the multi-row systems. 7
8 Supplementary Figure 6: a) set-up of the model and the minimum energy configuration. CA, M and TPA molecules are colored in red, blue and green, repsectively; Lateral sliding analysis of PES at 0 rotation, including Lennard-Jones and Coulombic contributions (b); at 10 (c); at 20 (d); at 30 (e) 8
9 Supplementary Figure 7: a) the interaction potential energy for 14 TPA rows as function of MD simulation time calculated for different initial rotation angles between the TPA bundle and CA.M lattice. b) a snapshot of MD simulation showing that for the initial rotation angle of 10 the ends of TPA bundle spontaneously align along the preferred orientation of 0 during equilibration time. CA, M and TPA molecules are colored in red, blue and green, respectively. The underlying hbn layer was removed for sake of clarity. In MD simulations, the three-layered system (as shown in Figure 5a of the main text) was equilibrated for 4 ns in an NVT ensemble and then propagated for a further 2 ns with an NVE integrator. For the initial rotation angle of 0 o, 10 o, 20 o and 30 o between the TPA row bundle and CA.M lattice vector, the interaction potential energy of TPA rows was recorded every 10 fs and averaged every 100 fs (Supplementary Figure 7a). MD simulations show that the orientation at 0 is the most stable, by about 0.04 ev/molecule. This energy value is somewhat smaller than that found for the rigid body scan because the increased flexibility allows the assembly in the non-optimal orientations accommodated through local distortions. It is also interesting to note that in simulations with the initial angle of 10, both ends of the TPA bundle rotated to the preferred orientation of 0 during the equilibration time, thus lowering the potential energy of interaction (Supplementary Figure 7b) and demonstrating a tendency for the rows to align along the CA.M lattice vector. A complete rotation of the whole assembly is difficult to observe in a short time scale of MD simulation as it implies simultaneous rotations of all the rows, which are strongly interacting with one another. These simulations merely provide a comparison between different orientations of the TPA bundle on CA.M/hBN support, and are not an attempt to model a direct assembly of TPA rows. 9
10 In Supplementary Fig. 8a and b two example of snapshots from MD runs with different relative displacement in neighboring rows are shown. The configuration in Supplementary Fig. 8b has the correct inter-row separation, but not the right displacement. The configuration in Supplementary Fig. 8a shows a displacement closer to the experimental results, but the inter-row separation is slightly broader. In both cases we can observe defect formation in terms of row dislocations. As mentioned in the main text, this could be due to a an intrinsic limitation of the classical potential. We tried to improve it by recalculating the TPA charges from DFT-based RESP calculations. In the Standard OPLS set of charges all the atoms in the phenyl ring have the same negative charge ( e). In our new set those attached to the carboxylic group are slightly positive (0.095 e). This make sense from a chemical perspective given the electron-withdrawing nature of the COOH group. The C atom of the COOH group becomes less positive (0.532 e instead of e) while all the other atom types remain very close to the original OPLS values (the largest difference is e). With this new set of charges we can obtain a perfect match with the experimental results for both inter-row separation and displacement as shown in Supplementary Fig 8c and d. 10
11 Supplementary Figure 8: a-b) two examples of different configurations of TPA rows taken from different MD simulations of the three-layered system at 300 K using the standard OPLS set of charges for TPA. CA, M and TPA molecules are colored in red, blue and green, respectively. The underlying hbn layer was removed for sake of clarity. c) and d) show contact mode images of the same area showing rows of TPA molecules (d). At higher deflection setpoint (c) the underlying CA.M is revealed and shows that every sixth TPA row is positioned over every fourth CA.M row. The angle between the CA.M and TPA lattice direction is highlited. e) a snaphost taken from MD simulations with new set of RESP charges for TPA. Highlighted angle show the perfect match with experimental results. Similarly, a TMA molecular network was studied. A hexagonal arrangement of 486 TMA molecules was formed and optimised in vacuum. The gain in energy corresponding to hydrogen bond formation is the same as in the case of TPA (0.27 ev). The TMA molecular arrangement was placed 0.35 nm above the CA.M layer (Supplementary Figure 9a). Lateral 11
12 sliding analysis of the PES clearly identified the preferred angle of 30 between the CA.M and TMA lattices, for which the total energy and different contributions to it are shown in Supplementary Figure 9c. The energy profiles at 0, 10 and 20 are shown in Supplementary Figure 9d-f. The adsorption energy profile (see below) indicates a minor secondary minimum at 0. In all other orientations, the PES is essentially flat. As in the case of TPA, most of the adsorption energy in the preferred orientation comes from the vdw interaction (Supplementary Figure 9c). However, for TMA two distinct minima in the total energy profile are present: one is due to the vdw interaction and corresponds to a TMA molecule positioned on-top of M, while the other is driven by Coulombic interaction and corresponds to TMA being on-top of CA. Geometry optimisation in the two positions gives an almost identical adsorption energy of about ev/molecule in both cases. To investigate further these two positions, the three-layer system with TMA on-top of CA and on top of M was studied using MD simulations. The potential energy of TMA molecules in both cases was recorded, as described above for the case of TPA, and is shown in Supplementary Figure 10a. The configuration with TMA on-top M is slightly more stable, but the difference in energy is comparable to kt. 12
13 Supplementary Figure 9 a-b) set-up of the model and the minimum energy configuration. CA, M and TPA molecules are colored in red, blue and green, repsectively; Lateral sliding analysis of PES at 30 rotation, including Lennard-Jones and coulombic contributions (c); at 0 (d); at 10 (e); at 20 (f) 13
14 Supplementary Figure 10: a) the interaction potential energy as function of MD simulation time, calculated for TMA on top of CA and on top of M; (b) snapshot from the equilibrated MD run of the full three-layered system with TMA on-top of M. CA, M and TPA molecules are colored in red, blue and green, respectively. Underlaying hbn layer was removed for sake of clarity. Supplementary Figure 11: the adsorption energy as a function of the angle between the the CA.M lattice vector and principal axis of the TPA row (a) or the TMA lattice vector (b). 14
15 Supplementary References 1. Plimpton, S. Fast parallel algorithms for short-range molecular dynamics. J. Comput. Phys. 117, 1 19 (1995). 2. Jorgensen, W. L., Maxwell, D. S. & Tirado-Rives, J. Development and Testing of the OLPS All-Atom Force Field on Conformational Energetics and Properties of Organic Liquids. J. Am. Chem. Soc. 118, (1996). 3. Kannappan, K., Werblowsky, T. L., Rim, K. T., Berne, B. J. & Flynn, G. W. An experimental and theoretical study of the formation of nanostructures of selfassembled cyanuric acid through hydrogen bond networks on graphite. J. Phys. Chem. B 111, (2007). 4. Conti, S. & Cecchini, M. Accurate and Efficient Calculation of the Desorption Energy of Small Molecules from Graphene. J. Phys. Chem. C 119, (2015). 5. Bayly, C. C. I., Cieplak, P., Cornell, W. D. & Kollman, P. a. A well-behaved electrostatic potential based method using charge restraints for deriving atomic charges: the RESP model. J. Phys. 97, (1993). 6. Hutter, J., Iannuzzi, M., Schiffmann, F. & VandeVondele, J. CP2K: atomistic simulations of condensed matter systems. WIREs Comput Mol Sci 4, (2014). 7. Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 38, (1988). 8. Lee, C., Yang, W. & Parr, R. G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B 37, (1988). 9. Lippert, G., Hutter, J. & Parrinello, M. A hybrid Gaussian and plane wave density functional scheme. Mol. Phys. 92, (1997). 10. Vandevondele, J. et al. Quickstep: Fast and accurate density functional calculations using a mixed Gaussian and plane waves approach. Comput. Phys. Commun. 167, (2005). 11. VandeVondele, J. & Hutter, J. Gaussian basis sets for accurate calculations on molecular systems in gas and condensed phases. J. Chem. Phys. 127, (2007). 12. Goedecker, S., Teter, M. & Hutter, J. Separable Dual-Space Gaussian Pseudopotentials. Phys. Rev. B 54, (1996). 13. Hartwigsen, C., Goedecker, S. & Hutter, J. Relativistic separable dual-space Gaussian pseudopotentials from H to Rn. Phys. Rev. B 58, (1998). 14. Krack, M. Pseudopotentials for H to Kr optimized for gradient-corrected exchangecorrelation functionals. Theor. Chem. Acc. 114, (2005). 15. Martyna, G. J. & Tuckerman, M. E. A reciprocal space based method for treating long range interactions in ab initio and force-field-based calculations in clusters. J. Chem. Phys. 110, (1999). 16. Hockney, R. W. & Eastwood, J. W. Computer simulation using particles. (Adam Hilger, 1989). 17. Tuckerman, M. E., Alejandre, J., López-Rendón, R., Jochim, A. L. & Martyna, G. J. A Liouville-operator derived measure-preserving integrator for molecular dynamics simulations in the isothermal isobaric ensemble. J. Phys. A. Math. Gen. 39, (2006). 18. Shinoda, W., Shiga, M. & Mikami, M. Rapid estimation of elastic constants by molecular dynamics simulation under constant stress. Phys. Rev. B - Condens. Matter Mater. Phys. 69, (2004). 19. Martyna, G. J., Tobias, D. J. & Klein, M. L. Constant pressure molecular dynamics algorithms. J. Chem. Phys. 101, 4177 (1994). 20. Humphrey, W., Dalke, A. & Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. 14, 33 8, NaN, (1996). 21. Jewett, A. I., Zhuang, Z. & Shea, J.-E. Moltemplate a Coarse-Grained Model Assembly Tool. Biophys. J. 104, 169 (2013). 22. Lee, J. H. A Study on a Boron-Nitride Nanotube as a Gigahertz Oscillator. J. Korean 15
16 Phys. Soc. 49, (2006). 23. Leroy, F., Liu, S. & Zhang, J. Parametrizing Nonbonded Interactions from Wetting Experiments via the Work of Adhesion: Example of Water on Graphene Surfaces. J. Phys. Chem. C 119, (2015). 24. Nasrabadi, A. T. & Foroutan, M. Interactions between polymers and single-walled boron nitride nanotubes: A molecular dynamics simulation approach. J. Phys. Chem. B 114, (2010). 25. Won, C. Y. & Aluru, N. R. Structure and dynamics of water confined in a boron nitride nanotube. J. Phys. Chem. B 112, (2008). 26. Golze, D., Hutter, J. & Iannuzzi, M. Wetting of water on hexagonal boron nitride@rh(111): a QM/MM model based on atomic charges derived for nanostructured substrates. Phys. Chem. Chem. Phys. 17, (2015). 16
Universal Repulsive Contribution to the. Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information
Universal Repulsive Contribution to the Solvent-Induced Interaction Between Sizable, Curved Hydrophobes: Supporting Information B. Shadrack Jabes, Dusan Bratko, and Alenka Luzar Department of Chemistry,
More informationIon-Gated Gas Separation through Porous Graphene
Online Supporting Information for: Ion-Gated Gas Separation through Porous Graphene Ziqi Tian, Shannon M. Mahurin, Sheng Dai,*,, and De-en Jiang *, Department of Chemistry, University of California, Riverside,
More informationSupplementary Figure 1 Irregular arrangement of E,E-8-mer on TMA. STM height images formed when
Supplementary Figure 1 Irregular arrangement of E,E-8-mer on TMA. STM height images formed when a 5 µl heptanoic acid solution of E,E-8-mer is applied on: (a) a TMA templated HOPG substrate, U s = +1.00
More informationSupporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure
Supporting Information for Solid-liquid Thermal Transport and its Relationship with Wettability and the Interfacial Liquid Structure Bladimir Ramos-Alvarado, Satish Kumar, and G. P. Peterson The George
More informationShear Properties and Wrinkling Behaviors of Finite Sized Graphene
Shear Properties and Wrinkling Behaviors of Finite Sized Graphene Kyoungmin Min, Namjung Kim and Ravi Bhadauria May 10, 2010 Abstract In this project, we investigate the shear properties of finite sized
More informationSupporting Information: Surface Polarons Reducing Overpotentials in. the Oxygen Evolution Reaction
Supporting Information: Surface Polarons Reducing Overpotentials in the Oxygen Evolution Reaction Patrick Gono Julia Wiktor Francesco Ambrosio and Alfredo Pasquarello Chaire de Simulation à l Echelle Atomique
More informationSupporting Information
alladium monophosphine d(h 3 ): does it really exist in solution? ietro Vidossich*, Gregori Ujaque,* and Agustí Lledós* Supporting Information Contents: 1. Schemes S1 and S2 2. Computational details 3.
More informationSupplementary Figure 1 Experimental setup for crystal growth. Schematic drawing of the experimental setup for C 8 -BTBT crystal growth.
Supplementary Figure 1 Experimental setup for crystal growth. Schematic drawing of the experimental setup for C 8 -BTBT crystal growth. Supplementary Figure 2 AFM study of the C 8 -BTBT crystal growth
More informationSupporting Information
Projection of atomistic simulation data for the dynamics of entangled polymers onto the tube theory: Calculation of the segment survival probability function and comparison with modern tube models Pavlos
More informationSupporting information for: Charge Distribution and Fermi Level in. Bimetallic Nanoparticles
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supporting information for: Charge Distribution and Fermi Level in Bimetallic Nanoparticles
More informationSupplementary Figure 1. Schematic of rapid thermal annealing process: (a) indicates schematics and SEM cross-section of the initial layer-by-layer
Supplementary Figure 1. Schematic of rapid thermal annealing process: (a) indicates schematics and SEM cross-section of the initial layer-by-layer film configuration, (b) demonstrates schematic and cross-section
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationMechanical and Thermal Stability of Graphyne and Graphdiyne Nanoscrolls
Mechanical and Thermal Stability of Graphyne and Graphdiyne Nanoscrolls Daniel Solis 1, Cristiano F. Woellner 1,2, Daiane D. Borges 1, and Douglas S. Galvao 1 1 Applied Physics Department, University of
More informationMolecular Dynamics Simulations. Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia
Molecular Dynamics Simulations Dr. Noelia Faginas Lago Dipartimento di Chimica,Biologia e Biotecnologie Università di Perugia 1 An Introduction to Molecular Dynamics Simulations Macroscopic properties
More informationSupporting information
Electronic Supplementary Material ESI for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2018 Supporting information Understanding three-body contributions to coarse-grained force
More informationSupporting Information for. Dynamics Study"
Supporting Information for "CO 2 Adsorption and Reactivity on Rutile TiO 2 (110) in Water: An Ab Initio Molecular Dynamics Study" Konstantin Klyukin and Vitaly Alexandrov,, Department of Chemical and Biomolecular
More informationSupporting Information
Supporting Information Interface-Induced Affinity Sieving in Nanoporous Graphenes for Liquid-Phase Mixtures Yanan Hou, Zhijun Xu, Xiaoning Yang * State Key Laboratory of Material-Orientated Chemical Engineering,
More informationMolecular Dynamics Simulation of a Nanoconfined Water Film
Molecular Dynamics Simulation of a Nanoconfined Water Film Kyle Lindquist, Shu-Han Chao May 7, 2013 1 Introduction The behavior of water confined in nano-scale environment is of interest in many applications.
More informationElectronic Supporting Information Topological design of porous organic molecules
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2017 Electronic Supporting Information Topological design of porous organic molecules Valentina Santolini,
More informationIntroduction to molecular dynamics
1 Introduction to molecular dynamics Yves Lansac Université François Rabelais, Tours, France Visiting MSE, GIST for the summer Molecular Simulation 2 Molecular simulation is a computational experiment.
More informationk θ (θ θ 0 ) 2 angles r i j r i j
1 Force fields 1.1 Introduction The term force field is slightly misleading, since it refers to the parameters of the potential used to calculate the forces (via gradient) in molecular dynamics simulations.
More informationSupplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water
Supplementary Information for Electronic signature of the instantaneous asymmetry in the first coordination shell in liquid water Thomas D. Kühne 1, 2 and Rustam Z. Khaliullin 1, 1 Institute of Physical
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationSupporting Information
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supporting Information Order in the chaos: The secret of the large negative entropy
More informationDiffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study
Diffusion of Water and Diatomic Oxygen in Poly(3-hexylthiophene) Melt: A Molecular Dynamics Simulation Study Julia Deitz, Yeneneh Yimer, and Mesfin Tsige Department of Polymer Science University of Akron
More informationSupporting Information for. Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite. Edge Water Interface
Supporting Information for Ab Initio Metadynamics Study of VO + 2 /VO2+ Redox Reaction Mechanism at the Graphite Edge Water Interface Zhen Jiang, Konstantin Klyukin, and Vitaly Alexandrov,, Department
More informationWater structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
Bull. Mater. Sci., Vol. 31, No. 3, June 2008, pp. 525 532. Indian Academy of Sciences. Water structure near single and multi-layer nanoscopic hydrophobic plates of varying separation and interaction potentials
More informationLarge scale growth and characterization of atomic hexagonal boron. nitride layers
Supporting on-line material Large scale growth and characterization of atomic hexagonal boron nitride layers Li Song, Lijie Ci, Hao Lu, Pavel B. Sorokin, Chuanhong Jin, Jie Ni, Alexander G. Kvashnin, Dmitry
More informationExperiment Section Fig. S1 Fig. S2
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Materials Experiment Section The STM experiments were carried out in an ultrahigh
More informationExploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials
Exploring the Changes in the Structure of α-helical Peptides Adsorbed onto Carbon and Boron Nitride based Nanomaterials Dr. V. Subramanian Chemical Laboratory, IPC Division CSIR-Central Leather Research
More informationSUPPLEMENTARY INFORMATION
In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION DOI: 10.1038/NCHEM.2720 1 2 3 Tuning underwater adhesion with cation-π interactions Matthew A. Gebbie, Wei Wei, Alex M. Schrader,
More informationSUPPLEMENTARY INFORMATION
SUPPLEMENTARY INFORMATION Supplementary Information Figure S1: (a) Initial configuration of hydroxyl and epoxy groups used in the MD calculations based on the observations of Cai et al. [Ref 27 in the
More informationSupporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes
Supporting information of Self-Assembly, Structure and Pi- Conjugation at the Interface of Poly-3-Hexylthiophene and Carbon Nanotubes Details of the molecular dynamics simulations and computational procedures
More informationSupporting Information
Supporting Information Ionic Liquid Designed for PEDOT:PSS Conductivity Enhancement Ambroise de Izarra, a,b,1 Seongjin Park, a,1 Jinhee Lee, a Yves Lansac, b,c, * Yun Hee Jang a, * a Department of Energy
More informationMultiwalled nanotube faceting unravelled
Multiwalled nanotube faceting unravelled Itai Leven, Roberto Guerra, Andrea Vanossi, Erio Tosatti, and Oded Hod NATURE NANOTECHNOLOGY www.nature.com/naturenanotechnology 1 1. Achiral DWCNTs optimized using
More informationDetermination of Kamlet-Taft parameters for selected solvate ionic liquids.
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Determination of Kamlet-Taft parameters for selected solvate ionic liquids. Daniel
More informationIntermolecular Forces in Density Functional Theory
Intermolecular Forces in Density Functional Theory Problems of DFT Peter Pulay at WATOC2005: There are 3 problems with DFT 1. Accuracy does not converge 2. Spin states of open shell systems often incorrect
More informationInter-Layer Potential for Graphene/h-BN Heterostructures. Supplementary Information
Inter-Layer Potential for Graphene/h-BN Heterostructures Supplementary Information Itai Leven, 1 Tal Maaravi, 1 Ido Azuri, 2 Leeor Kronik, 2 and Oded Hod 1 1 Department of Physical Chemistry, School of
More informationSupporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 2016 Supporting Information Soft Nanoparticles: Nano Ionic Networks of Associated Ionic Polymers Dipak
More informationAu-C Au-Au. g(r) r/a. Supplementary Figures
g(r) Supplementary Figures 60 50 40 30 20 10 0 Au-C Au-Au 2 4 r/a 6 8 Supplementary Figure 1 Radial bond distributions for Au-C and Au-Au bond. The zero density regime between the first two peaks in g
More informationMOLECULAR DYNAMICS SIMULATION OF HYDROGEN STORAGE IN SINGLE-WALLED CARBON NANOTUBES
MOLECULAR DYNAMICS SIMULATION OF HYDROGEN STORAGE IN SINGLE-WALLED CARBON NANOTUBES Shigeo MARUYAMA Engineering Research Institute The University of Tokyo 2-11-16 Yayoi, Bunkyo-ku, Tokyo 113-8656, Japan
More informationMechanical Properties of Phagraphene Membranes: A Fully Atomistic Molecular Dynamics Investigation
Mechanical Properties of Phagraphene Membranes: A Fully Atomistic Molecular Dynamics Investigation J. M. de Sousa 1,2, A. L. Aguiar 2, E. C. Girão 2, Alexandre F. Fonseca 1, A. G. Sousa Filho 3, and Douglas
More informationElectronic excitations in conjugated. Many-Body Green's Functions. Behnaz Bagheri Varnousfaderani. Behnaz Bagheri Varnousfaderani
Technische Universiteit Eindhoven University of Technology Electronic excitations in conjugated polymers Electronic from excitations QM/MM simulations in conjugated pol based from on Many-Body QM/MM simulations
More informationNew Six-site Acetonitrile Model for Simulations of Liquid Acetonitrile and its Aqueous Mixtures
New Six-site Acetonitrile Model for Simulations of Liquid Acetonitrile and its Aqueous Mixtures ALEXEI M. NIKITIN, 1,2 ALEXANDER P. LYUBARTSEV 2 1 Engelhardt Institute of Molecular Biology Russian Academy
More informationCP2K: Past, Present, Future. Jürg Hutter Department of Chemistry, University of Zurich
CP2K: Past, Present, Future Jürg Hutter Department of Chemistry, University of Zurich Outline Past History of CP2K Development of features Present Quickstep DFT code Post-HF methods (RPA, MP2) Libraries
More informationState Key Laboratory of Molecular Reaction Dynamics, Dalian Institute of Chemical Physics, Chinese Academic of Sciences, Dalian , P. R. China.
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2015 Supplementary Information for: Theoretical exploration of MgH 2 and graphene nano-flake in cyclohexane:
More informationWater-methanol separation with carbon nanotubes and electric fields
Electronic Supplementary Material (ESI) for Nanoscale. This journal is The Royal Society of Chemistry 215 Supplementary Information: Water-methanol separation with carbon nanotubes and electric fields
More informationClar Sextet Theory for low-dimensional carbon nanostructures: an efficient approach based on chemical criteria
Clar Sextet Theory for low-dimensional carbon nanostructures: an efficient approach based on chemical criteria Matteo Baldoni Fachbereich Chemie, Technische Universität Dresden, Germany Department of Chemistry
More informationHands-on : Model Potential Molecular Dynamics
Hands-on : Model Potential Molecular Dynamics OUTLINE 0. DL_POLY code introduction 0.a Input files 1. THF solvent molecule 1.a Geometry optimization 1.b NVE/NVT dynamics 2. Liquid THF 2.a Equilibration
More informationSelf assembly of Carbon Atoms in Interstellar Space and Formation Mechanism of Naphthalene from Small Precursors: A Molecular Dynamics Study
Self assembly of Carbon Atoms in Interstellar Space and Formation Mechanism of Naphthalene from Small Precursors: A Molecular Dynamics Study Niladri Patra UIC Advisor: Dr. H. R. Sadeghpour ITAMP Large
More informationTheoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes
J. At. Mol. Sci. doi: 10.4208/jams.121011.011412a Vol. 3, No. 4, pp. 367-374 November 2012 Theoretical comparative study on hydrogen storage of BC 3 and carbon nanotubes Xiu-Ying Liu a,, Li-Ying Zhang
More informationSite dependent hydrogenation in Graphynes: A Fully Atomistic Molecular Dynamics Investigation
Site dependent hydrogenation in Graphynes: A Fully Atomistic Molecular Dynamics Investigation Pedro A. S. Autreto and Douglas S. Galvao Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas,
More informationSupplementary Information
Supplementary Information a b Supplementary Figure 1. Morphological characterization of synthesized graphene. (a) Optical microscopy image of graphene after transfer on Si/SiO 2 substrate showing the array
More informationSupplementary Information. Surface Microstructure Engenders Unusual Hydrophobicity in. Phyllosilicates
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Supplementary Information Surface Microstructure Engenders Unusual Hydrophobicity in Phyllosilicates
More informationSupporting Information
Supporting Information Structure and Dynamics of Uranyl(VI) and Plutonyl(VI) Cations in Ionic Liquid/Water Mixtures via Molecular Dynamics Simulations Katie A. Maerzke, George S. Goff, Wolfgang H. Runde,
More informationSupporting Information for: Physics Behind the Water Transport through. Nanoporous Graphene and Boron Nitride
Supporting Information for: Physics Behind the Water Transport through Nanoporous Graphene and Boron Nitride Ludovic Garnier, Anthony Szymczyk, Patrice Malfreyt, and Aziz Ghoufi, Institut de Physique de
More informationSupporting information
Supporting information Toward a Janus Cluster: Regiospecific Decarboxylation of Ag 44 (4- MBA) 30 @Ag Nanoparticles Indranath Chakraborty, Anirban Som, Tuhina Adit Maark, Biswajit Mondal, Depanjan Sarkar
More informationScientific Computing II
Scientific Computing II Molecular Dynamics Simulation Michael Bader SCCS Summer Term 2015 Molecular Dynamics Simulation, Summer Term 2015 1 Continuum Mechanics for Fluid Mechanics? Molecular Dynamics the
More informationMotion of Light Adatoms and Molecules on the Surface of Few-Layer Graphene
Supporting information Motion of Light Adatoms and Molecules on the Surface of Few-Layer Graphene Franziska Schäffel 1,*, Mark Wilson 2, Jamie H. Warner 1 1 Department of Materials, University of Oxford,
More informationMolecular Dynamics Simulation of Chemical Sputtering of Hydrogen Atom on Layer Structured Graphite
1 TH/7-1 Molecular Dynamics Simulation of Chemical Sputtering of Hydrogen Atom on Layer Structured Graphite A. Ito 1,2), Y. Wang 1), S. Irle 1), K. Morokuma 3), and H. Nakamura 2) 1) Nagoya University,
More informationSupporting Information
Supporting Information Controlled Growth of Ceria Nanoarrays on Anatase Titania Powder: A Bottom-up Physical Picture Hyun You Kim 1, Mark S. Hybertsen 2*, and Ping Liu 2* 1 Department of Materials Science
More informationExplanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph.
Supplementary Materials Explanation of Dramatic ph-dependence of Hydrogen Binding on Noble Metal Electrode: Greatly Weakened Water Adsorption at High ph. Tao Cheng,, Lu Wang, Boris V Merinov, and William
More informationInitial Stages of Growth of Organic Semiconductors on Graphene
Initial Stages of Growth of Organic Semiconductors on Graphene Presented by: Manisha Chhikara Supervisor: Prof. Dr. Gvido Bratina University of Nova Gorica Outline Introduction to Graphene Fabrication
More informationTemperature induced structural phase transitions in a twodimensional
Temperature induced structural phase transitions in a twodimensional self-assembled network Matthew O. Blunt*, Jinne Adisoejoso, Kazukuni Tahara, Keisuke Katayama, Mark Van der Auweraer*, Yoshito Tobe*,
More informationSupporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene
Supporting information for Polymer interactions with Reduced Graphene Oxide: Van der Waals binding energies of Benzene on defected Graphene Mohamed Hassan, Michael Walter *,,, and Michael Moseler, Freiburg
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.2491 Experimental Realization of Two-dimensional Boron Sheets Baojie Feng 1, Jin Zhang 1, Qing Zhong 1, Wenbin Li 1, Shuai Li 1, Hui Li 1, Peng Cheng 1, Sheng Meng 1,2, Lan Chen 1 and
More informationComparative Study of the Water Response to External Force at Nanoscale and Mesoscale
Copyright 2013 Tech Science Press CMES, vol.95, no.4, pp.303-315, 2013 Comparative Study of the Water Response to External Force at Nanoscale and Mesoscale H.T. Liu 1,2, Z. Chen 2, S. Jiang 2, Y. Gan 3,
More informationReport on Atomistic Modeling of Bonding in Carbon-Based Nanostructures
Report on Atomistic Modeling of Bonding in Carbon-Based Nanostructures Timothy Stillings Department of Physics, Astronomy and Materials Science Missouri State University Advisor: Ridwan Sakidja Abstract
More informationSupporting Information. by Hexagonal Boron Nitride
Supporting Information High Velocity Saturation in Graphene Encapsulated by Hexagonal Boron Nitride Megan A. Yamoah 1,2,, Wenmin Yang 1,3, Eric Pop 4,5,6, David Goldhaber-Gordon 1 * 1 Department of Physics,
More informationThe experimental work seems to be well carried out and the DFT calculations carefully crafted.
Reviewers' comments: Reviewer #1 (Remarks to the Author): The manuscript by Maier et al. reports on the on-surface synthesis of 1D and 2D polymers in ultra-high vacuum. A halogenated triphenylamine precursor
More informationSUPPLEMENTARY INFORMATION
DOI: 10.1038/NCHEM.2161 Coulomb explosion during the early stages of the reaction of alkali metals with water Philip E. Mason, 1 Frank Uhlig, 1 Václav Vaněk, 1 Tillmann Buttersack, 2 Sigurd Bauerecker,
More informationWATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION
WATER PERMEATION THROUGH GRAPHENE NANOSLIT BY MOLECULAR DYNAMICS SIMULATION Taro Yamada 1 and Ryosuke Matsuzaki 2 1 Department of Mechanical Engineering, Tokyo University of Science, 2641 Yamazaki, Noda,
More informationSupplementary Figure 1. Electron micrographs of graphene and converted h-bn. (a) Low magnification STEM-ADF images of the graphene sample before
Supplementary Figure 1. Electron micrographs of graphene and converted h-bn. (a) Low magnification STEM-ADF images of the graphene sample before conversion. Most of the graphene sample was folded after
More informationDislocation network structures in 2D bilayer system
Dislocation network structures in 2D bilayer system Shuyang DAI School of Mathematics and Statistics Wuhan University Joint work with: Prof. Yang XIANG, HKUST Prof. David SROLOVITZ, UPENN S. Dai IMS Workshop,
More informationComputer simulation methods (2) Dr. Vania Calandrini
Computer simulation methods (2) Dr. Vania Calandrini in the previous lecture: time average versus ensemble average MC versus MD simulations equipartition theorem (=> computing T) virial theorem (=> computing
More informationUnit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement
Unit Cell-Level Thickness Control of Single-Crystalline Zinc Oxide Nanosheets Enabled by Electrical Double Layer Confinement Xin Yin, Yeqi Shi, Yanbing Wei, Yongho Joo, Padma Gopalan, Izabela Szlufarska,
More informationA MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES
A MOLECULAR DYNAMICS STUDY OF POLYMER/GRAPHENE NANOCOMPOSITES Anastassia N. Rissanou b,c*, Vagelis Harmandaris a,b,c* a Department of Applied Mathematics, University of Crete, GR-79, Heraklion, Crete,
More informationModeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic Dynamics
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2015 Supplementary Data for Modeling Ultrafast Deactivation in Oligothiophenes via Nonadiabatic
More informationCARBON 2004 Providence, Rhode Island. Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores
CARBON Providence, Rhode Island Adsorption of Flexible n-butane and n-hexane on Graphitized Thermal Carbon Black and in Slit Pores D. D. Do* and H. D. Do, University of Queensland, St. Lucia, Qld 7, Australia
More informationThe samples used in these calculations were arranged as perfect diamond crystal of
Chapter 5 Results 5.1 Hydrogen Diffusion 5.1.1 Computational Details The samples used in these calculations were arranged as perfect diamond crystal of a2 2 2 unit cells, i.e. 64 carbon atoms. The effect
More informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
More informationAdsorption of Iodine on Pt(111) surface. Alexandre Tkachenko Marcelo Galván Nikola Batina
Adsorption of Iodine on Pt(111) surface Alexandre Tkachenko Marcelo Galván Nikola Batina Outline Motivation Experimental results Geometry Ab initio study Conclusions Motivation Unusual structural richness
More information(3) E LJ (r ab ) = 4ɛ r ab
Hybrid Ensembles for Improved Force Matching Lee-Ping Wang and Troy Van Voorhis Department of Chemistry, Massachusetts Institute of Technology 77 Massachusetts Ave. Cambridge, MA 2139 USA Dated: 24 April
More informationMolecular Dynamics Study of Carbon Dioxide Storage in Carbon-Based Organic Nanopores
Molecular Dynamics Study of Carbon Dioxide Storage in Carbon-Based Organic Nanopores Mohammad Kazemi and Ali Takbiri-Borujeni, West Virginia University Copyright 2016, Society of Petroleum Engineers This
More informationDepartment of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an
The statistical characteristics of static friction J. Wang, G. F. Wang*, and W. K. Yuan Department of Engineering Mechanics, SVL, Xi an Jiaotong University, Xi an 710049, China * E-mail: wanggf@mail.xjtu.edu.cn
More informationSupporting Material for. Microscopic origin of gating current fluctuations in a potassium channel voltage sensor
Supporting Material for Microscopic origin of gating current fluctuations in a potassium channel voltage sensor J. Alfredo Freites, * Eric V. Schow, * Stephen H. White, and Douglas J. Tobias * * Department
More informationDehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways from ab initio simulations
2nd CP2K Tutorial, ETHZ 7-11 February 2011 Materials Science & Technology Nanostructures and adsorption on metallic surfaces: Dehydrogenation of cyclohexaphenylene on Cu(111): possible reaction pathways
More informationNon-equilibrium molecular dynamics simulations of organic friction modifiers
Non-equilibrium molecular dynamics simulations of organic friction modifiers James Ewen PhD Student Tribology Group (Shell UTC) Imperial College London Session 8B Bronze 2 Lubrication Fundamentals STLE
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for ChemComm. This journal is The Royal Society of Chemistry 2018 Electronic Supplementary Information Self-assembly of conformationaly flexible molecules under
More informationAmino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information
Amino Acids and Proteins at ZnO-water Interfaces in Molecular Dynamics Simulations: Electronic Supplementary Information Grzegorz Nawrocki and Marek Cieplak Institute of Physics, Polish Academy of Sciences,
More information1. Robust hexagonal rings on Cu(111) Figure S1 2. Details of Monte Carlo simulations
Supporting Information for Influence of Relativistic Effects on Assembled Structures of V-Shaped Bispyridine Molecules on M(111) Surfaces where M = Cu, Ag, Au Xue Zhang, 1,ǁ Na Li, 1, Hao Wang, 1 Chenyang
More informationJASS Modeling and visualization of molecular dynamic processes
JASS 2009 Konstantin Shefov Modeling and visualization of molecular dynamic processes St Petersburg State University, Physics faculty, Department of Computational Physics Supervisor PhD Stepanova Margarita
More informationAri P Seitsonen CNRS & Université Pierre et Marie Curie, Paris
Self-organisation on noble metal surfaces Ari P Seitsonen CNRS & Université Pierre et Marie Curie, Paris Collaborations Alexandre Dmitriev, Nian Lin, Johannes Barth, Klaus Kern,... Thomas Greber, Jürg
More information2D MBE Activities in Sheffield. I. Farrer, J. Heffernan Electronic and Electrical Engineering The University of Sheffield
2D MBE Activities in Sheffield I. Farrer, J. Heffernan Electronic and Electrical Engineering The University of Sheffield Outline Motivation Van der Waals crystals The Transition Metal Di-Chalcogenides
More informationThe Molecular Dynamics Method
The Molecular Dynamics Method Thermal motion of a lipid bilayer Water permeation through channels Selective sugar transport Potential Energy (hyper)surface What is Force? Energy U(x) F = d dx U(x) Conformation
More informationSupplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers
Supplementary Information for: Controlling Cellular Uptake of Nanoparticles with ph-sensitive Polymers Hong-ming Ding 1 & Yu-qiang Ma 1,2, 1 National Laboratory of Solid State Microstructures and Department
More informationBiochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015,
Biochemistry,530:,, Introduc5on,to,Structural,Biology, Autumn,Quarter,2015, Course,Informa5on, BIOC%530% GraduateAlevel,discussion,of,the,structure,,func5on,,and,chemistry,of,proteins,and, nucleic,acids,,control,of,enzyma5c,reac5ons.,please,see,the,course,syllabus,and,
More informationheight trace of a 2L BN mechanically exfoliated on SiO 2 /Si with pre-fabricated micro-wells. Scale bar 2 µm.
Supplementary Figure 1. Few-layer BN nanosheets. AFM image and the corresponding height trace of a 2L BN mechanically exfoliated on SiO 2 /Si with pre-fabricated micro-wells. Scale bar 2 µm. Supplementary
More informationAB INITIO MOLECULAR-DYNAMICS SIMULATIONS OF ADSORPTION OF DYE MOLECULES AT SURFACES
AB INITIO MOLECULAR-DYNAMICS SIMULATIONS OF ADSORPTION OF DYE MOLECULES AT SURFACES M. SUGIHARA, H. MEYER, AND P. ENTEL Theoretische Tieftemperaturphysik, Universität Duisburg, 47048 Duisburg, Germany
More informationSupporting Information for. Structural and Chemical Dynamics of Pyridinic Nitrogen. Defects in Graphene
Supporting Information for Structural and Chemical Dynamics of Pyridinic Nitrogen Defects in Graphene Yung-Chang Lin, 1* Po-Yuan Teng, 2 Chao-Hui Yeh, 2 Masanori Koshino, 1 Po-Wen Chiu, 2 Kazu Suenaga
More information