Formation mechanism of polycyclic aromatic hydrocarbons in benzene combustion: Quantum chemical molecular dynamics simulations

Molecular Dynamics Study of Nanoribbon Formation by Encapsulating Cyclic Hydrocarbon Molecules inside Single-Walled Carbon Nanotube

Carbon nanotubes filled with organic molecules can serve as chemical nanoreactors. Recent experimental results show that, by introducing cyclic hydrocarbon molecules inside carbon nanotubes, they can be transformed into nanoribbons or inner tubes, depending on the experimental conditions. In this paper, we present our results obtained as a continuation of our previous molecular dynamics simulation work. In our previous work, the initial geometry consisted of independent carbon atoms. Now, as an initial condition, we have placed different molecules inside a carbon nanotube (18,0): C5H5 (fragment of ferrocene), C5, C5+H2; C6H6 (benzene), C6, C6+H2; C20H12 (perylene); and C24H12 (coronene). The simulations were performed using the REBO-II potential of the LAMMPS software package, supplemented with a Lennard-Jones potential between the nanotube wall atoms and the inner atoms. The simulation proved difficult due to the slow dynamics of the H abstraction. However, with a slight modification of the parameterization, it was possible to model the formation of carbon nanoribbons inside the carbon nanotube.

Insights on the initial stages of carbonization of sub-bituminous coal

Detailed studies were carried out on the initial stages of carbonization of sub-bituminous coal using reactive force field (ReaxFF) molecular dynamics simulations. Evolution of different gaseous species during carbonization were analyzed at different temperatures and densities. Elementary reactions were identified leading to the formation of small gaseous species. Cleavage of homolytic O-H bond was found to be the first step of sub-bituminous coal carbonization. CH₄ formed mainly due to reaction of ·CH₃ with hydrogen radical/hydrogen abstraction. Among others, C₉H_mO_n radicals were the most abundant species at any simulations condition considered here, where m = 9, 8 or 7 and n = 1, 2 or 3. Bond dissociation energy (BDE) of the identified reactions decreased by ∼3.0 kcal/mol for endothermic reactions and increased by ∼3.0 kcal/mol for exothermic due to change in enthalpy at higher temperature. It was observed that the formation of 5-membered carbon only rings and sp-hybridized carbon species played crucial role in the subsequent growth. The results agree with available experimental observations and computational studies. This work provides new insights on the carbonization mechanism of sub-bituminous coal.

Theoretical investigation of the electronic relaxation in highly excited chrysene and tetracene: The effect of armchair vs zigzag edge

Non-adiabatic molecular dynamics of neutral chrysene and tetracene molecules is investigated using Tully's fewest switches surface hopping algorithm coupled to the time-dependent density functional based tight-binding (TD-DFTB) method for electronic structure calculations. We first assess the performance of two DFTB parameter sets based on the computed TD-DFTB absorption spectra. The main focus is given to the analysis of the electronic relaxation from the brightest excited state following absorption of a UV photon. We determine the dynamical relaxation times and discuss the underlying mechanisms. Our results show that the electronic population of the brightest excited singlet state in armchair-edge chrysene decays an order-of-magnitude faster than the one in zigzag-edge tetracene. This is correlated with a qualitatively similar difference of energy gaps between the brightest state and the state lying just below in energy, which is also consistent with our previous study on polyacenes.

Density-functional tight-binding: basic concepts and applications to molecules and clusters

The scope of this article is to present an overview of the Density Functional based Tight Binding (DFTB) method and its applications. The paper introduces the basics of DFTB and its standard formulation up to second order. It also addresses methodological developments such as third order expansion, inclusion of non-covalent interactions, schemes to solve the self-interaction error, implementation of long-range short-range separation, treatment of excited states via the time-dependent DFTB scheme, inclusion of DFTB in hybrid high-level/low level schemes (DFT/DFTB or DFTB/MM), fragment decomposition of large systems, large scale potential energy landscape exploration with molecular dynamics in ground or excited states, non-adiabatic dynamics. A number of applications are reviewed, focusing on -(i)- the variety of systems that have been studied such as small molecules, large molecules and biomolecules, bare orfunctionalized clusters, supported or embedded systems, and -(ii)- properties and processes, such as vibrational spectroscopy, collisions, fragmentation, thermodynamics or non-adiabatic dynamics. Finally outlines and perspectives are given.

Simulations of the synthesis of boron-nitride nanostructures in a hot, high pressure gas volume

We performed nanosecond timescale computer simulations of clusterization and agglomeration processes of boron nitride (BN) nanostructures in hot, high pressure gas, starting from eleven different atomic and molecular precursor systems containing boron, nitrogen and hydrogen at various temperatures from 1500 to 6000 K. The synthesized BN nanostructures self-assemble in the form of cages, flakes, and tubes as well as amorphous structures. The simulations facilitate the analysis of chemical dynamics and we are able to predict the optimal conditions concerning temperature and chemical precursor composition for controlling the synthesis process in a high temperature gas volume, at high pressure. We identify the optimal precursor/temperature choices that lead to the nanostructures of highest quality with the highest rate of synthesis, using a novel parameter of the quality of the synthesis (PQS). Two distinct mechanisms of BN nanotube growth were found, neither of them based on the root-growth process. The simulations were performed using quantum-classical molecular dynamics (QCMD) based on the density-functional tight-binding (DFTB) quantum mechanics in conjunction with a divide-and-conquer (DC) linear scaling algorithm, as implemented in the DC-DFTB-K code, enabling the study of systems as large as 1300 atoms in canonical NVT ensembles for 1 ns time.

Effect of ammonia on chemical vapour deposition and carbon nanotube nucleation mechanisms

Chemical vapour deposition (CVD) growth of carbon nanotubes is currently the most viable method for commercial-scale nanotube production. However, controlling the 'chirality', or helicity, of carbon nanotubes during CVD growth remains a challenge. Recent studies have shown that adding chemical 'etchants', such as ammonia and water, to the feedstock gas can alter the diameter and chirality of nanotubes produced with CVD. To date, this strategy for chirality control remains sub-optimal, since we have a poor understanding of how these etchants change the CVD and nucleation mechanisms. Here, we show how ammonia alters the mechanism of methane CVD and single-walled carbon nanotube nucleation on iron catalysts, using quantum chemical molecular dynamics simulations. Our simulations reveal that ammonia is selectively activated by the catalyst, and this enables ammonia to play a dual role during methane CVD. Following activation, ammonia nitrogen removes carbon from the catalyst surface exclusively via the production of hydrogen (iso)cyanide, thus impeding the growth of extended carbon chains. Simultaneously, ammonia hydrogen passivates carbon dangling bonds, which impedes nanotube nucleation and promotes defect healing. Combined, these effects lead to slower, more controllable nucleation and growth kinetics.

Quantum Chemical Molecular Dynamics Simulations of 1,3-Dichloropropene Combustion

Oxidative decomposition of 1,3-dichloropropene was investigated using quantum chemical molecular dynamics (QM/MD) at 1500 and 3000 K. Thermal oxidation of 1,3-dichloropropene was initiated by (1) abstraction of allylic H/Cl by O2 and (2) intra-annular C-Cl bond scission and elimination of allylic Cl. A kinetic analysis shows that (2) is the more dominant initiation pathway, in agreement with QM/MD results. These QM/MD simulations reveal new routes to the formation of major products (H2O, CO, HCl, CO2), which are propagated primarily by the chloroperoxy (ClO2), OH, and 1,3-dichloropropene derived radicals. In particular, intra-annular C-C/C-H bond dissociation reactions of intermediate aldehydes/ketones are shown to play a dominant role in the formation of CO and CO2. Our simulations demonstrate that both combustion temperature and radical concentration can influence the product yield, however not the combustion mechanism.

Fullerenes patched by flowers with octagonal core

The aromatic character of three different flowers of general formula [n:(p i,p j)n/2], namely [8:(5,7)4], [8:(6,6)4] and [8:(5,6)4], has been evaluated by means of geometric (HOMA index), energetic (heats of formation) and magnetic criteria (NICS index, exaltation of magnetic susceptibility). Also, the reactivity descriptors within the DFT approach -absolute hardness, electrophilicity, Fukui functions — have been computed. All the different methods used for estimating the aromaticity led to a unitary conclusion. Two fullerene structures, patched by the most stable 8-Sumanene flower, have been designed and evaluated by means of ab initio computations

Computational investigations on covalent dimerization/oligomerization of polyacenes: is it relevant to soot formation?

We have postulated a novel pathway that could assist in the nucleation of soot particles through covalent dimerization and oligomerizations of a variety of PAHs. DFT calculations were performed with the objective of obtaining the relative thermal stabilities and formation probabilities of oligomeric species that exploit the facile dimerization that is known to occur in linear oligoacenes. We propose that the presence of small stretches of linear oligoacence (tetracene or longer) in extended PAH, either embedded or tethered, would be adequate for enabling the formation of such dimeric and oligomeric adducts; these could then serve as nuclei for the growth of soot particles. Our studies also reveal the importance of π-stacking interactions between extended aromatic frameworks in governing the relative stabilities of the oligomeric species that are formed.