No Arabic abstract
We present an implementation in a linear-scaling density-functional theory code of an electronic enthalpy method, which has been found to be natural and efficient for the ab initio calculation of finite systems under hydrostatic pressure. Based on a definition of the system volume as that enclosed within an electronic density isosurface [Phys. Rev. Lett., 94, 145501 (2005)], it supports both geometry optimizations and molecular dynamics simulations. We introduce an approach for calibrating the parameters defining the volume in the context of geometry optimizations and discuss their significance. Results in good agreement with simulations using explicit solvents are obtained, validating our approach. Size-dependent pressure-induced structural transformations and variations in the energy gap of hydrogenated silicon nanocrystals are investigated, including one comparable in size to recent experiments. A detailed analysis of the polyamorphic transformations reveals three types of amorphous structures and their persistence on depressurization is assessed.
We demonstrate that Daubechies wavelets can be used to construct a minimal set of optimized localized contracted basis functions in which the Kohn-Sham orbitals can be represented with an arbitrarily high, controllable precision. Ground state energies and the forces acting on the ions can be calculated in this basis with the same accuracy as if they were calculated directly in a Daubechies wavelets basis, provided that the amplitude of these contracted basis functions is sufficiently small on the surface of the localization region, which is guaranteed by the optimization procedure described in this work. This approach reduces the computational costs of DFT calculations, and can be combined with sparse matrix algebra to obtain linear scaling with respect to the number of electrons in the system. Calculations on systems of 10,000 atoms or more thus become feasible in a systematic basis set with moderate computational resources. Further computational savings can be achieved by exploiting the similarity of the contracted basis functions for closely related environments, e.g. in geometry optimizations or combined calculations of neutral and charged systems.
In orbital-free density functional theory the kinetic potential (KP), the functional derivative of the kinetic energy density functional, appears in the Euler equation for the electron density and may be more amenable to simple approximations. We study properties of two solid-state systems, Al and Si, using two nonlocal KPs that gave good results for atoms. Very accurate results are found for Al, but results for Si are much less satisfactory, illustrating the general need for a better treatment of extended covalent systems. A different integration pathway in the KP formalism may prove useful in attacking this fundamental problem.
We present a quantitative exploration, combining experiment and simulation, of the mechanical and electronic properties, as well as the modifications induced by an alkylthiolated coating, at the single NP level. We determine the response of the NPs to external pressure in a controlled manner by using an atomic force microscope tip. We find a strong reduction of their Young modulus, as compared to bulk gold, and a significant influence of strain in the electronic properties of the alkylthiolated NPs. Electron transport measurements of tiny molecular junctions (NP/alkylthiol/CAFM tip) show that the effective tunnelling barrier through the adsorbed monolayer strongly decreases with increasing the applied load, which translates in a remarkable and unprecedented increase of the tunnel current. These observations are successfully explained using simulations based on finite element analysis (FEA) and first-principles calculations that permit to consider the coupling between the mechanical response of the system and the electric dipole variations at the interface.
We propose a computationally efficient approach to the nonadiabatic time-dependent density functional theory (TDDFT) which is based on a representation of the frequency-dependent exchange correlation kernel as a response of a set of damped oscillators. The requirements to computational resources needed to implement our approach do not differ from those of the standard real-time TDDFT in the adiabatic local density approximation (ALDA). Thus, our result offers an exciting opportunity to take into account temporal nonlocality and memory effects in calculations with TDDFT in quantum chemistry and solid state physics for unprecedentedly low costs.
Although poorly understood, cation-exchange reactions are increasingly used to dope or transform colloidal semiconductor nanocrystals (quantum dots). We used density-functional theory and kinetic Monte Carlo simulations to develop a microscopic theory that explains structural, optical, and electronic changes observed experimentally in Ag-cation-exchanged CdSe nanocrystals. We find that Coulomb interactions, both between ionized impurities and with the polarized nanocrystal surface, play a key role in cation exchange. Our theory also resolves several experimental puzzles related to photoluminescence and electrical behavior in CdSe nanocrystals doped with Ag.