Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userCharge transfer excitations with range separated functionals using improved virtual orbitals
70
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
We present an implementation of range separated functionals utilizing the Slater-function on grids in real space in the projector augmented waves method. The screened Poisson equation is solved to evaluate the necessary screened exchange integrals on Cartesian grids. The implementation is verified against existing literature and applied to the description of charge transfer excitations. We find very slow convergence for calculations within linear response time-dependent density functional theory and unoccupied orbitals of the canonical Fock operator. Convergence can be severely improved by using Huzinagas virtual orbitals instead. This combination furthermore enables an accurate determination of long-range charge transfer excitations by means of ground-state calculations.
rate research
Read More
We present a solution of the full TDDFT eigenvalue equation in the linear response formalism exhibiting a linear-scaling computational complexity with system size, without relying on the simplifying Tamm-Dancoff approximation (TDA). The implementation relies on representing the occupied and unoccupied subspace with two different sets of in situ optimised localised functions, yielding a very compact and efficient representation of the transition density matrix of the excitation with the accuracy associated with a systematic basis set. The TDDFT eigenvalue equation is solved using a preconditioned conjugate-gradients algorithm that is very memory-efficient. The algorithm is validated on a test molecule and a good agreement with results obtained from standard quantum chemistry packages is found, with the preconditioner yielding a significant improvement in convergence rates. The method developed in this work is then used to reproduce experimental results of the absorption spectrum of bacteriochlorophyll (BChl) in an organic solvent, where it is demonstrated that the TDA fails to reproduce the main features of the low energy spectrum, while the full TDDFT equation yields results in good qualitative agreement with experimental data. Furthermore, the need for explicitly including parts of the solvent into the TDDFT calculations is highlighted, making the treatment of large system sizes necessary that are well within reach of the capabilities of the algorithm introduced here. Finally, the linear-scaling properties of the algorithm are demonstrated by computing the lowest excitation energy of BChl in solution. The largest systems considered in this work are of the same order of magnitude as a variety of pigment-protein complexes, opening up the possibility of studying their properties without having to resort to any semiclassical approximations to parts of the protein environment.
Koopmans-compliant (KC) functionals have been shown to provide accurate spectral properties through a generalized condition of piece-wise linearity of the total energy as a function of the fractional addition/removal of an electron to/from any orbital. We analyze the performance of different KC functionals on the GW100 test-set, comparing the ionization potentials (as opposite of the energy of the highest occupied orbital) of these 100 molecules to those obtained from CCSD(T) total energy differences, and experimental results, finding excellent agreement with a mean absolute error of 0.20 eV for the KIPZ functional, that is state-of-the-art for both DFT-based calculations and many-body perturbation theory. We highlight similarities and differences between KC functionals and other electronic-structure approaches, such as dielectric-dependent hybrid functionals and G$_0$W$_0$, both from a theoretical and from a practical point of view, arguing that Koopmans-compliant potentials can be considered as a local and orbital-dependent counterpart to the electronic GW self-energy, albeit already including approximate vertex corrections.
Motivated by an experiment in which the singlet-triplet gap in triphenylene based copolymers was effectively tuned, we used time dependent density functional theory (TDDFT) to reproduce the main results. By means of conventional and long-range corrected exchange correlation functionals, the luminescence energies and the exciton localization were calculated for a triphenylene homopolymer and several different copolymers. The phosphorescence energy of the pure triphenylene chain is predicted accurately by means of the optimally tuned long-range corrected LC-PBE functional and slightly less accurate by the global hybrid B3LYP. However, the experimentally observed fixed phosphorescence energy could not be reproduced because the localization pattern is different to the expectations: Instead of localizing on the triphenylene moiety - which is present in all types of polymers - the triplet state localizes on the different bridging units in the TDDFT calculations. This leads to different triplet emission energies for each type of polymer. Yet, there are clear indications that long-range corrected TDDFT has the potential to predict the triplet emission energies as well as the localization behavior more accurate than conventional local or semi-local functionals.
The traditional fluid perturbation theory is improved by taking electronic excitations and ionizations into account, in the framework of average ion spheres. It is applied to calculate the equation of state for fluid Xenon, which turns out in good agreement with the available shock data.
We propose a fast method for the calculation of short-range interactions in molecular dynamics simulations. The so-called random-batch list method is a stochastic version of the classical neighbor-list method to avoid the construction of a full Verlet list, which introduces two-level neighbor lists for each particle such that the neighboring particles are located in core and shell regions, respectively. Direct interactions are performed in the core region. For the shell zone, we employ a random batch of interacting particles to reduce the number of interaction pairs. The error estimate of the algorithm is provided. We investigate the Lennard-Jones fluid by molecular dynamics simulations, and show that this novel method can significantly accelerate the simulations with a factor of several fold without loss of the accuracy. This method is simple to implement, can be well combined with any linked cell methods to further speed up and scale up the simulation systems, and can be straightforwardly extended to other interactions such as Ewald short-range part, and thus it is promising for large-scale molecular dynamics simulations.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
