No abstract

DFTB+ is a versatile community developed open source software package offering fast and efficient methods for carrying out atomistic quantum mechanical simulations. By implementing various methods approximating density functional theory (DFT), such as the density functional based tight binding (DFTB) and the extended tight binding method, it enables simulations of large systems and long timescales with reasonable accuracy while being considerably faster for typical simulations than the respective ab initio methods. Based on the DFTB framework, it additionally offers approximated versions of various DFT extensions including hybrid functionals, time dependent formalism for treating excited systems, electron transport using non-equilibrium Green’s functions, and many more. DFTB+ can be used as a user-friendly standalone application in addition to being embedded into other software packages as a library or acting as a calculation-server accessed by socket communication. We give an overview of the recently developed capabilities of the DFTB+ code, demonstrating with a few use case examples, discuss the strengths and weaknesses of the various features, and also discuss on-going developments and possible future perspectives.

Deep eutectic solvents (DESs) are a mixture of a salt and a molecular hydrogen bond donor, which form a eutectic liquid with a depressed melting point. Quantum mechanical molecular dynamics (QM/MD) simulations have been used to probe the 1 : 2 choline chloride-urea (ChCl : U), choline chloride-ethylene glycol (ChCl : EG) and choline chloride-glycerol (ChCl : Gly) DESs. DES nanostructure and interactions between the ions is used to rationalise differences in DES eutectic point temperatures and viscosity. Simulations show that the structure of the bulk hydrogen bond donor is largely preserved for hydroxyl based hydrogen bond donors (ChCl:Gly and ChCl:EG), resulting in a smaller melting point depression. By contrast, ChCl:U exhibits a well-established hydrogen bond network between the salt and hydrogen bond donor, leading to a larger melting point depression. This extensive hydrogen bond network in ChCl:U also leads to substantially higher viscosity, compared to ChCl:EG and ChCl:Gly. Of the two hydroxyl based DESs, ChCl:Gly also exhibits a higher viscosity than ChCl:EG. This is attributed to the over-saturation of hydrogen bond donor groups in the ChCl:Gly bulk, which leads to more extensive hydrogen bond donor self-interaction and hence higher cohesive forces within the bulk liquid.

The atomic scale details of single-walled carbon nanotube (SWNT) nucleation on metal catalyst particles are elusive to experimental observations. Computer simulation of metal-catalyzed SWNT nucleation is a challenging topic but potentially of great importance to understand the factors affecting SWNT diameters, chirality, and growth efficiency. In this work, we use nonequilibrium density functional tight-binding molecular dynamics simulations and report nucleation of sp(2)-carbon cap structures on an iron particle consisting of 38 atoms. One C(2) molecule was placed every 1.0 ps around an Fe(38) cluster for 30 ps, after which a further 410 ps of annealing simulation without carbon supply was performed. We find that sp(2)-carbon network nucleation and annealing processes occur in three sequential and repetitive stages: (A) polyyne chains on the metal surface react with each other to evolve into a Y-shaped polyyne junction, which preferentially form a five-membered ring as a nucleus; (B) polyyne chains on the first five-membered ring form an additional fused five- or six-membered ring; and (C) pentagon-to-hexagon self-healing rearrangement takes place with the help of short-lived polyyne chains, stabilized by the mobile metal atoms. The observed nucleation process resembles the formation of a fullerene cage. However, the metal particle plays a key role in differentiating the nucleation process from fullerene cage formation, most importantly by keeping the growing cap structure from closing into a fullerene cage and by keeping the carbon edge "alive" for the addition of new carbon material.

Since their discovery in the early 1990s, single-walled carbon nanotubes (SWNTs) have spawned previously unimaginable commercial and industrial technologies. Their versatility stems from their unique electronic, physical/chemical, and mechanical properties, which set them apart from traditional materials. Many researchers have investigated SWNT growth mechanisms in the years since their discovery. The most prevalent of these is the vapor-liquid-solid (VLS) mechanism, which is based on experimental observations. Within the VLS mechanism, researchers assume that the formation of a SWNT starts with co-condensation of carbon and metal atoms from vapor to form liquid metal carbide. Once the liquid reaches supersaturation, the solid phase nanotubes begin to grow. The growth process is partitioned into three distinct stages: nucleation of a carbon "cap-precursor," "cap-to-tube" transformation, and continued SWNT growth. In recent years, molecular dynamics (MD) simulations have come to the fore with respect to SWNT growth. MD simulations lead to spatial and temporal resolutions of these processes that are superior to those possible using current experimental techniques, and so provide valuable information regarding the growth process that researchers cannot obtain experimentally. In this Account, we review our own recent efforts to simulate SWNT nucleation, growth, and healing phenomena on transition-metal catalysts using quantum mechanical molecular dynamics (QM/MD) methods. In particular, we have validated each stage of the SWNT condensation mechanism using a self-consistent-charge density-functional tight-binding (SCC-DFTB) methodology. With respect to the nucleation of a SWNT cap-precursor (stage 1), we have shown that the presence of a transition-metal carbide particle is not a necessary prerequisite for SWNT nucleation, contrary to conventional experimental presumptions. The formation and coalescence of polyyne chains on the metal surface occur first, followed by the formation of the SWNT cap-precursor, "ring condensation", and the creation of an sp(2)-hybridized carbon structure. In our simulations, the nucleation process takes approximately 400 ps. This first step occurs over a much longer time scale than the second stage of SWNT condensation (approximately 50 ps). We therefore observe SWNT nucleation to be akin to the rate-limiting step of the SWNT formation process. In addition to the QM/MD simulation of various stages of SWNT nucleation, growth, and healing processes, we have determined the effects of temperature, catalyst composition, and catalyst size on the kinetics and mechanism of SWNT growth. With respect to temperature dependence, we observe a "sweet-spot" with respect to the efficiency of SWNT growth. In addition, Ni-catalyzed SWNT growth is observed to be 70-100% faster compared to Fe-catalyzed SWNT growth, depending on the catalyst particle size. We also observe a noticeable increase in SWNT growth rates using smaller catalyst particles. Finally, we review our recent QM/MD investigation of SWNT heali...