Using Density Functional Theory to Model Realistic TiO2 Nanoparticles, Their Photoactivation and Interaction with Water

Computational modeling of titanium dioxide nanoparticles of realistic size is extremely relevant for the direct comparison with experiments but it is also a rather demanding task. We have recently worked on a multistep/scale procedure to obtain global optimized minimum structures for chemically stable spherical titania nanoparticles of increasing size, with diameter from 1.5 nm (~300 atoms) to 4.4 nm (~4000 atoms). We use first self-consistent-charge density functional tight-binding (SCC-DFTB) methodology to perform thermal annealing simulations to obtain globally optimized structures and then hybrid density functional theory (DFT) to refine them and to achieve high accuracy in the description of structural and electronic properties. This allows also to assess SCC-DFTB performance in comparison with DFT(B3LYP) results. As a further step, we investigate photoexcitation and photoemission processes involving electron/hole pair formation, separation, trapping and recombination in the nanosphere of medium size by hybrid DFT. Finally, we show how a recently defined new set of parameters for SCC-DFTB allows for a proper description of titania/water multilayers interface, which paves the way for modeling large realistic nanoparticles in aqueous environment

Failure of hydrogenation in protecting polycyclic aromatic hydrocarbons from fragmentation

A recent study of soft X-ray absorption in native and hydrogenated coronene cations, C$_{24}$H$_{12+m}^+$ $m=0-7$, led to the conclusion that additional hydrogen atoms protect (interstellar) Polycyclic Aromatic Hydrocarbon (PAH) molecules from fragmentation [Reitsma et al., Phys. Rev. Lett. 113, 053002 (2014)]. The present experiment with collisions between fast (30-200 eV) He atoms and pyrene (C$_{16}$H$_{10+m}^+$, $m=0$, 6, and 16) and simulations without reference to the excitation method suggests the opposite. We find that the absolute carbon-backbone fragmentation cross section does not decrease but increases with the degree of hydrogenation for pyrene molecules.Comment: 10 pages, 5 figure

Path integral evaluation of equilibrium isotope effects

A general and rigorous methodology to compute the quantum equilibrium isotope effect is described. Unlike standard approaches, ours does not assume separability of rotational and vibrational motions and does not make the harmonic approximation for vibrations or rigid rotor approximation for the rotations. In particular, zero point energy and anharmonicity effects are described correctly quantum mechanically. The approach is based on the thermodynamic integration with respect to the mass of isotopes and on the Feynman path integral representation of the partition function. An efficient estimator for the derivative of free energy is used whose statistical error is independent of the number of imaginary time slices in the path integral, speeding up calculations by a factor of 60 at 500 K. We describe the implementation of the methodology in the molecular dynamics package Amber 10. The method is tested on three [1,5] sigmatropic hydrogen shift reactions. Because of the computational expense, we use ab initio potentials to evaluate the equilibrium isotope effects within the harmonic approximation, and then the path integral method together with semiempirical potentials to evaluate the anharmonicity corrections. Our calculations show that the anharmonicity effects amount up to 30% of the symmetry reduced reaction free energy. The numerical results are compared with recent experiments of Doering and coworkers, confirming the accuracy of the most recent measurement on 2,4,6,7,9-pentamethyl-5-(5,5-$^2$H$_2$)methylene-11,11a-dihydro-12H-naphthacene as well as concerns about compromised accuracy, due to side reactions, of another measurement on 2-methyl-10-(10,10-$^2$H$_2$)methylenebicyclo[4.4.0]dec-1-ene.Comment: 14 pages, 8 figures, 6 table

Ewald summation on a helix : a route to self-consistent charge density-functional based tight-binding objective molecular dynamics

We explore the generalization to the helical case of the classical Ewald method, the harbinger of all modern self-consistent treatments of waves in crystals, including ab initio electronic structure methods. Ewald-like formulas that do not rely on a unit cell with translational symmetry prove to be numerically tractable and able to provide the crucial component needed for coupling objective molecular dynamics with the self-consistent charge density-functional based tight-binding treatment of the inter-atomic interactions. The robustness of the method in addressing complex hetero-nuclear nano- and bio-systems is demonstrated with illustrative simulations on a helical boron nitride nanotube, a screw dislocated zinc oxide nanowire, and an ideal DNA molecule

Mapping the structural diversity of C60 carbon clusters and their infrared spectra

The current debate about the nature of the carbonaceous material carrying the infrared (IR) emission spectra of planetary and proto-planetary nebulae, including the broad plateaus, calls for further studies on the interplay between structure and spectroscopy of carbon-based compounds of astrophysical interest. The recent observation of C60 buckminsterfullerene in space suggests that carbon clusters of similar size may also be relevant. In the present work, broad statistical samples of C60 isomers were computationally determined without any bias using a reactive force field, their IR spectra being subsequently obtained following local optimization with the density-functional-based tight-binding theory. Structural analysis reveals four main structural families identified as cages, planar polycyclic aromatics, pretzels, and branched. Comparison with available astronomical spectra indicates that only the cage family could contribute to the plateau observed in the 6-9 micron region. The present framework shows great promise to explore and relate structural and spectroscopic features in more diverse and possibly hydrogenated carbonaceous compounds, in relation with astronomical observations

Structure and electronic structure of Metal-Organic Frameworks within the Density-Functional based Tight-Binding method

Density-functional based tight-binding is a powerful method to describe large molecules and materials. Metal-Organic Frameworks (MOFs), materials with interesting catalytic properties and with very large surface areas have been developed and have become commercially available. Unit cells of MOFs typically include hundreds of atoms, which make the application of standard Density-Functional methods computationally very expensive, sometimes even unfeasible. The aim of this paper is to prepare and to validate the Self-Consistent Charge Density-Functional based Tight Binding (SCC-DFTB) method for MOFs containing Cu, Zn and Al metal centers. The method has been validated against full hybrid density-functional calculations for model clusters, against gradient corrected density-functional calculations for supercells, and against experiment. Moreover, the modular concept of MOF chemistry has been discussed on the basis of their electronic properties. We concentrate on MOFs comprising three common connector units: copper paddlewheels (HKUST-1), zinc oxide Zn4O tetrahedron (MOF-5, MOF-177, DUT-6 (MOF-205)) and aluminium oxide AlO4(OH)2 octahedron (MIL-53). We show that SCC-DFTB predicts structural parameters with a very good accuracy (with less than 5% deviation, even for adsorbed CO and H2O on HKUST-1), while adsorption energies differ by 12 kJ mol-1 or less for CO and water compared to DFT benchmark calculations.Comment: Submitted to Phys. Status Solidi

Mechanical Properties and Fracture Dynamics of Silicene Membranes

As graphene became one of the most important materials today, there is a renewed interest on others similar structures. One example is silicene, the silicon analogue of graphene. It share some the remarkable graphene properties, such as the Dirac cone, but presents some distinct ones, such as a pronounced structural buckling. We have investigated, through density functional based tight-binding (DFTB), as well as reactive molecular dynamics (using ReaxFF), the mechanical properties of suspended single-layer silicene. We calculated the elastic constants, analyzed the fracture patterns and edge reconstructions. We also addressed the stress distributions, unbuckling mechanisms and the fracture dependence on the temperature. We analysed the differences due to distinct edge morphologies, namely zigzag and armchair