Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userStatic dielectric response of icosahedral fullerenes from C60 to C2160 by an all electron density functional theory
122
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
The static dielectric response of C60, C180, C240, C540, C720, C960, C1500, and C2160 fullerenes is characterized by an all-electron density-functional method. First, the screened polarizabilities of C60, C180, C240, and C540, are determined by the finite-field method using Gaussian basis set containing 35 basis functions per atom. In the second set of calculations, the unscreened polarizabilities are calculated for fullerenes C60 through C2160 from the self-consistent Kohn-Sham orbitals and eigen-values using the sum-over-states method. The approximate screened polarizabilities, obtained by applying a correction determined within linear response theory show excellent agreement with the finite-field polarizabilities. The static dipole polarizability per atom in C2160 is (4 Angstrom^3) three times larger than that in C60 (1.344 Angstrom^3). Our results reduce the uncertainty in various theoretical models used previously to describe the dielectric response of fullerenes and show that quantum size effects in polarizability are significantly smaller than previously thought.
rate research
Read More
We present a benchmark of the density functional linear response calculation of NMR shieldings within the Gauge-Including Projector-Augmented-Wave method against all-electron Augmented-Plane-Wave$+$local-orbital and uncontracted Gaussian basis set results for NMR shieldings in molecular and solid state systems. In general, excellent agreement between the aforementioned methods is obtained. Scalar relativistic effects are shown to be quite large for nuclei in molecules in the deshielded limit. The small component makes up a substantial part of the relativistic corrections.
We implement and benchmark the frozen core approximation, a technique commonly adopted in electronic structure theory to reduce the computational cost by means of mathematically fixing the chemically inactive core electron states. The accuracy and efficiency of this approach are well controlled by a single parameter, the number of frozen orbitals. Explicit corrections for the frozen core orbitals and the unfrozen valence orbitals are introduced, safeguarding against seemingly minor numerical deviations from the assumed orthonormality conditions of the basis functions. A speedup of over two-fold can be achieved for the diagonalization step in all-electron density-functional theory simulations containing heavy elements, without any accuracy degradation in terms of the electron density, total energy, and atomic forces. This is demonstrated in a benchmark study covering 103 materials across the periodic table, and a large-scale simulation of CsPbBr3 with 2,560 atoms. Our study provides a rigorous benchmark of the precision of the frozen core approximation (sub-meV per atom for frozen core orbitals below -200 eV) for a wide range of test cases and for chemical elements ranging from Li to Po. The algorithms discussed here are implemented in the open-source Electronic Structure Infrastructure software package.
The non-local van der Waals density functional (vdW-DF) has had tremendous success since its inception in 2004 due to its constraint-based formalism that is rigorously derived from a many-body starting point. However, while vdW-DF can describe binding energies and structures for van der Waals complexes and mixed systems with good accuracy, one long-standing criticism---also since its inception---has been that the $C_6$ coefficients that derive from the vdW-DF framework are largely inaccurate and can be wrong by more than a factor of two. It has long been thought that this failure to describe the $C_6$ coefficients is a conceptual flaw of the underlying plasmon framework used to derive vdW-DF. We prove here that this is not the case and that accurate $C_6$ coefficient can be obtained without sacrificing the accuracy at binding separations from a modified framework that is fully consistent with the constraints and design philosophy of the original vdW-DF formulation. Our design exploits a degree of freedom in the plasmon-dispersion model $omega_{mathbf{q}}$, modifying the strength of the long-range van der Waals interaction and the cross-over from long to short separations, with additional parameters tuned_ to reference systems. Testing the new formulation for a range of different systems, we not only confirm the greatly improved description of $C_6$ coefficients, but we also find excellent performance for molecular dimers and other systems. The importance of this development is not necessarily that particular aspects such as $C_6$ coefficients or binding energies are improved, but rather that our finding opens the door for further conceptual developments of an entirely unexplored direction within the exact same constrained-based non-local framework that made vdW-DF so successful in the first place.
The treatment of atomic anions with Kohn-Sham density functional theory (DFT) has long been controversial since the highest occupied molecular orbital (HOMO) energy, $E_{HOMO}$, is often calculated to be positive with most approximate density functionals. We assess the accuracy of orbital energies and electron affinities for all three rows of elements in the periodic table (H-Ar) using a variety of theoretical approaches and customized basis sets. Among all of the theoretical methods studied here, we find that a non-empirically tuned range-separated approach (constructed to satisfy DFT-Koopmans theorem for the anionic electron system) provides the best accuracy for a variety of basis sets - even for small basis sets where most functionals typically fail. Previous approaches to solve this conundrum of positive $E_{HOMO}$ values have utilized non-self-consistent methods; however electronic properties, such as electronic couplings/gradients (which require a self-consistent potential and energy), become ill-defined with these approaches. In contrast, the non-empirically tuned range-separated procedure used here yields well-defined electronic couplings/gradients and correct $E_{HOMO}$ values since both the potential and resulting electronic energy are computed self-consistently. Orbital energies and electron affinities are further analyzed in the context of the electronic energy as a function of electronic number (including fractional numbers of electrons) to provide a stringent assessment of self-interaction errors for these complex anion systems.
We present a study of the optical response of compact and hollow icosahedral clusters containing up to 868 silver atoms by means of time-dependent density functional theory. We have studied the dependence on size and morphology of both the sharp plasmonic resonance at 3-4 eV (originated mainly from $sp$-electrons), and the less studied broader feature appearing in the 6-7 eV range (interband transitions). An analysis of the effect of structural relaxations, as well as the choice of exchange correlation functional (local density versus generalized gradient approximations) both in the ground state and optical response calculations is also presented. We have further analysed the role of the different atom layers (surface versus inner layers) and the different orbital symmetries on the absorption cross-section for energies up to 8 eV. We have also studied the dependence on the number of atom layers in hollow structures. Shells formed by a single layer of atoms show a pronounced red shift of the main plasmon resonances that, however, rapidly converge to those of the compact structures as the number of layers is increased. The methods used to obtain these results are also carefully discussed. Our methodology is based on the use of localized basis (atomic orbitals, and atom-centered- and dominant- product functions), which bring several computational advantages related to their relatively small size and the sparsity of the resulting matrices. Furthermore, the use of basis sets of atomic orbitals also brings the possibility to extend some of the standard population analysis tools (e.g., Mulliken population analysis) to the realm of optical excitations. Some examples of these analyses are described in the present work.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
