Virtual Winter School on Computational Chemistry

Cecam Logo

Speaker: Dr Stephan Irle
Institute: Oak Ridge National Laboratory
Country: USA
Speaker Link: https://www.ornl.gov/staff-profile/stephan-irle
Time: 17:00 CET 21-Feb-22

Dr Stephan Irle

Computational Chemistry and Nanomaterials Sciences Group, Oak Ridge National Laboratory, USA

The density-functional tight-binding (DFTB) method [1] is an approximation to density functional theory (DFT) and allows a speedup of first principles electronic structure calculations by two to three orders of magnitude.  In this talk, I will discuss DFTB-based simulations of nanoscale materials self-assembly in nonequilibrium on large length scales [2].  Fullerene, carbon nanotube, and graphene formation were simulated on the nanosecond time scale, considering experimental conditions as closely as possible.  An approximate density functional method was employed to compute energies and gradients on-the-fly in direct MD simulations, while the simulated systems were continually pushed away from equilibrium via carbon concentration or temperature gradients.  We find that carbon nanostructure formation from feedstock particles involves a phase transition of sp to sp2 carbon phases, which begins with the formation of Y-junctions, followed by a nucleus consisting of pentagons, hexagons, and heptagons.  The dominance of hexagons in the synthesized products is explained via annealing processes that occur during the cooling of the grown carbon structure, accelerated by transition metal catalysts when present.  The dimensional structures of the final synthesis products (0D spheres – fullerenes, 1D tubes – nanotubes, 2D sheets – graphenes) are induced by the shapes of the substrates/catalysts, and their interaction strength with carbon.  Our work prompts a paradigm shift away from traditional anthropomorphic formation mechanisms solely based on thermodynamic stability. Instead, we conclude that nascent carbon nanostructures at high temperatures are dissipative structures described by nonequilibrium dynamics in the manner proposed by Prigogine, Whitesides, and others.  As such, the fledgling carbon nanostructures consume energy while increasing the entropy of the environment, and only gradually anneal to achieve their familiar, final structure, maximizing hexagon formation wherever possible [2,3].

Recording:

Video is available only for registered users.

References:

[1] a) Christensen, A. S.; Kubar, T.; Cui, Q.; Elstner, M. Semiempirical Quantum Mechanical Methods for Noncovalent Interactions for Chemical and Biochemical Applications, Chem. Rev. 2016, 116, 5301-5337; b) http://www.dftbplus.org 

[2] Irle, S; Page, A. J.; Saha, B.; Wang, Y.; Chandrakumar, K. R. S.; Nishimoto, Y.; Qian, H.-J.; Morokuma, K. Atomistic mechanism of carbon nanostructure self-assembly as predicted by nonequilibrium QM/MD simulations, in: J. Leszczynski, M. K. Shukla, Eds. “Practical Aspects of Computational Chemistry II: An Overview of the Last Two Decades and Current Trends”, Springer-European Academy of Sciences, Chapter 5, pp. 105-172 (April 2, 2012).  ISBN 978-94-007-0922-5. DOI: 10.1007/978-94-007-0923-2_5 Preprint: https://www.dropbox.com/s/n2o3sjnb0t1z6mr/5_Online%20PDF.pdf?dl=0 

[3] Page, A. J.; Ding, F.; Irle, S.; Morokuma, K. Insights into carbon nanotube and graphene formation mechanisms from molecular simulations: a review, Rep. Prog. Phys. 2015, 78, 036501/1-38. 

 

One thought on “Quantum Chemical Simulations of Carbon Nanostructure Self-Assembly in Nonequilibrium”

  1. Monday, 21 February 2022 20:04

    From Stephan: Thank you for a lively and interesting discussion, I was amazed that the non-fashionable topic of nanocarbons could still spark such interest.  Here are the two papers I referred to when I cited new developments of DFTB for transition metal-catalyzed reactions:

    -Vujović M, Huynh M, Steiner S, Garcia-Fernandez P, Elstner M, Cui Q, Gruden M. Exploring the applicability of density functional tight binding to transition metal ions. Parameterization for nickel with the spin-polarized DFTB3 model. J Comput Chem. 2019 Jan 15;40(2):400-413. doi: 10.1002/jcc.25614.

    -Stepanovic S , Lai R , Elstner M , Gruden M , Garcia-Fernandez P , Cui Q . Improvement of d-d interactions in density functional tight binding for transition metal ions with a ligand field model: assessment of a DFTB3+U model on nickel coordination compounds. Phys Chem Chem Phys. 2020 Dec 7;22(46):27084-27095. doi: 10.1039/d0cp04694a.

    Here is the github for the DFTB+ code: https://github.com/dftbplus/dftbplus or simply download the executable/source from: https://dftbplus.org (Tutorial: https://dftbplus-recipes.readthedocs.io/en/latest/#)

    Regarding our parameterization tool, please send me an email at:  - As the code is still not well documented, we prefer to release it by email for now although the public github release is planned shortly.