Performance of different density functionals for the calculation of vibrational frequencies with vibrational coupled cluster method in bosonic representation

Photophysical properties of donor (D)-acceptor (A)-donor (D) diketopyrrolopyrrole (A) systems as donors for applications to organic electronic devices

Fourteen substituted diketopyrrolopyrrole (DPP) molecules in a donor (D)-acceptor (DPP)-donor (D) arrangement were designed. We employed density functional theory, time-dependent DFT, DFT-MRCI and the ab initio wave function second-order algebraic diagrammatic construction (ADC(2)) methods to investigate theoretically these systems. The examined aromatic substituents have one, two, or three hetero- and non-hetero rings. We comprehensively investigated their optical, electronic, and charge transport properties to evaluate potential applications in organic electronic devices. We found that the donor substituents based on one, two, or three aromatic rings bonded to the DPP core can improve the efficiency of an organic solar cell by fine-tuning the highest occupied molecular orbital/lowest unoccupied molecular orbital levels to match acceptors in typical bulk heterojunctions acceptors. Several properties of interest for organic photovoltaic devices were computed. We show that the investigated molecules are promising for applications as donor materials when combined with typical acceptors in bulk heterojunctions because they have appreciable energy conversion efficiencies resulting from their low ionization potentials and high electron affinities. This scenario allows a more effective charge separation and reduces the recombination rates. A comprehensive charge transfer analysis shows that D-A (DDP)-D systems have significant intramolecular charge transfer, further confirming their promise as candidates for donor materials in solar cells. The significant photophysical properties of DPP derivatives, including the high fluorescence emission, also allow these materials to be used in organic light-emitting diodes.

Harmonic and anharmonic vibrational computations for biomolecular building blocks: Benchmarking DFT and basis sets by theoretical and experimental IR spectrum of glycine conformers

Advanced vibrational spectroscopic experiments have reached a level of sophistication that can only be matched by numerical simulations in order to provide an unequivocal analysis, a crucial step to understand the structure-function relationship of biomolecules. While density functional theory (DFT) has become the standard method when targeting medium-size or larger systems, the problem of its reliability and accuracy are well-known and have been abundantly documented. To establish a reliable computational protocol, especially when accuracy is critical, a tailored benchmark is usually required. This is generally done over a short list of known candidates, with the basis set often fixed a priori. In this work, we present a systematic study of the performance of DFT-based hybrid and double-hybrid functionals in the prediction of vibrational energies and infrared intensities at the harmonic level and beyond, considering anharmonic effects through vibrational perturbation theory at the second order. The study is performed for the six-lowest energy glycine conformers, utilizing available "state-of-the-art" accurate theoretical and experimental data as reference. Focusing on the most intense fundamental vibrations in the mid-infrared range of glycine conformers, the role of the basis sets is also investigated considering the balance between computational cost and accuracy. Targeting larger systems, a broad range of hybrid schemes with different computational costs is also tested.

CO2 electroreduction on single atom catalysts: the role of the DFT functional

One key process involving single atom catalysts (SACs) is the electroreduction of CO2 to fuels. The chemistry of SACs differs largely from that of extended catalytic surfaces, presenting an opportunity to improve the ability to activate very stable molecules, such as CO2. In this work, we performed a density functional theory (DFT) study of CO2 activation on a series of SACs, focusing on the role played by the adopted functional in activity predictions. The role of the exchange-correlation functional has been widely investigated in heterogenous catalysts, but it is less explored in SACs. We tested the widely used PBE and the PBE+U corrected functionals against the more robust hybrid PBE0 functional. The results show that PBE is reliable if one is interested in qualitative predictions, but it leads to some inaccuracies in other cases. A possible way to attenuate this effect is by adopting the PBE+U framework, as it gives results that are very similar to PBE0 at an acceptable computational cost. The results of this study further underline the importance of the computational framework adopted in predicting the activity of SACs. The work suggests that one needs to go beyond PBE for quantitative estimates, an important consideration when performing screening and high-throughput calculations.

An approximation to the vibrational coupled-cluster method for CH-stretching of large molecules: application to naphthalene and anthracene

We propose an approximation to the vibrational coupled-cluster method (VCCM) to describe the CH-stretching region of the vibrational spectrum of large molecules. The vibrational modes of a molecule are divided into two sets: the target set and the bath set. The target set includes the CH stretches and the modes that are strongly coupled with the CH stretches and/or involve strong Fermi resonances with a CH stretch fundamental. The rest of the modes are in the bath set. First, the effective harmonic oscillator (EHO) approximation is invoked for the whole system to obtain the zeroth-order frequencies and modified potentials. The effects of interaction between the bath set and the target sets are included in the modified potential from the EHO calculation. The VCCM equations are constructed with the modified potential from the EHO calculations and for the target set only. The transition energies and intensities are calculated using such a truncated VCCM approximation. The proposed method is applied to calculate the IR spectra of naphthalene and anthracene. The results with three different criteria for selecting the modes in the target set are compared with the experimental IR spectra.

Analytical harmonic vibrational frequencies with VV10-containing density functionals: Theory, efficient implementation, and benchmark assessments

VV10 is a powerful nonlocal density functional for long-range correlation that is used to include dispersion effects in many modern density functionals, such as the meta-generalized gradient approximation (mGGA), B97M-V, the hybrid GGA, ωB97X-V, and the hybrid mGGA, ωB97M-V. While energies and analytical gradients for VV10 are already widely available, this study reports the first derivation and efficient implementation of the analytical second derivatives of the VV10 energy. The additional compute cost of the VV10 contributions to analytical frequencies is shown to be small in all but the smallest basis sets for recommended grid sizes. This study also reports the assessment of VV10-containing functionals for predicting harmonic frequencies using the analytical second derivative code. The contribution of VV10 to simulating harmonic frequencies is shown to be small for small molecules but important for systems where weak interactions are important, such as water clusters. In the latter cases, B97M-V, ωB97M-V, and ωB97X-V perform very well. The convergence of frequencies with respect to the grid size and atomic orbital basis set size is studied, and recommendations are reported. Finally, scaling factors to allow comparison of scaled harmonic frequencies with experimental fundamental frequencies and to predict zero-point vibrational energy are presented for some recently developed functionals (including r2SCAN, B97M-V, ωB97X-V, M06-SX, and ωB97M-V).

Insights on α-Glucose Biosensors/Carriers Based on Boron-Nitride Nanomaterials from an Atomistic and Electronic Point of View

The interaction of α-glucose with a BN-nanosheet, BN-nanotube, and BN-fullerene, was analyzed from an atomistic and electronic point of view, to evaluate such nanostructures as possible carriers and/or biosensors of the α-glucose molecule. Adsorption energies are in the range of physisorption (-0.79 eV to -0.91 eV) for the BN-nanosheet and -nanotube, and chemisorption (-2.24 eV to -2.35 eV), for the BN-fullerene. All systems, exhibit semiconductor-like behavior and great stability according to |LUMO-HOMO| energy gap [Gap_LH ] and chemical potential values, respectively. For the BN-nanosheet and -nanotube, the stabilization of the complexes is through hydrogen bonds, while for BN-fullerene is through a covalent bond and charge transfer. Furthermore, the BN-fullerene is able to dissociate the α-glucose molecule, which could help to decomposer such a compound, and be used for biological applications. The data taking into consideration solvent effects have no significant impact with respect to gas phase, except in the dipole moment (M_d ) where we noticed an increase up to ∼45 %. Our results suggest that BN-nanosheet and -nanotube, may act as biosensors, while BN-fullerene, may serve as a carrier or degrader of the α-glucose molecule.

Density inversion method for local basis sets without potential auxiliary functions: inverting densities from RDMFT

A density inversion method is presented, to obtain the constrained, optimal, local potential that has a prescribed asymptotic behaviour and reproduces optimally any given ground-state electronic density. This work builds...

Understanding of the C-H stretch region of infra-red spectroscopy: an analysis of the final state wavefunctions

The nature of the wavefunctions associated with the final states in the CH stretching region of several medium sized molecules is analysed. The number of optically bright transitions is much larger than the number of CH oscillators present in the molecule, and they are spread over a range of about 300 cm^-1. Several of them are clustered together within about 5 cm^-1 with near equal intensities. The final states of all these transitions are superpositions of multiple zeroth order states. In almost all of such superpositions, no single zeroth order state has more than 50% weight. Several multiquantum states, with three to four quanta of excitation dominate the final states, with the CH chromophore contributing only a small weightage. Thus the band structure of the CH stretch region is due to several optically bright transitions whose final states are superpositions of low frequency multiquantum states with the CH chromophore contributing only a small weight to make them spectroscopically active.

On the role of steric and exchange–correlation effects in halogenated complexes

Density functional theory formalisms of energy partitioning schemes are utilized to find out what energetic components govern interactions in halogenated complexes.

Can Density Functional Theory Be Trusted for High-Order Electric Properties? The Case of Hydrogen-Bonded Complexes

This work reports on an extensive assessment of the performance of a wide palette of density functional approximations in predicting the (high-order) electric properties of hydrogen-bonded complexes. To this end, we compute the electronic and vibrational contributions to the electric polarizability and the first and second hyperpolarizabilities, using the CCSD(T)/aug-cc-pVTZ level of theory as reference. For all the studied properties, the average absolute errors below 20% can only be obtained using the CAM-B3LYP functional, while LC-BLYP and MN15 are shown to be only slightly less accurate (average absolute errors not exceeding 30%). Among Minnesota density functionals, i.e., M06, M06-2X, and MN15, we only recommend the latter one, which quite accurately predicts the electronic and vibrational (hyper)polarizabilities. We also analyze the optimal tuning of the range-separation parameter μ for the LC-BLYP functional, finding that this approach does not bring any systematic improvement in the predictions of electronic and vibrational (hyper)polarizabilities and the accuracy of computed properties is largely system-dependent. Finally, we report huge errors in predicting the vibrational second hyperpolarizability by ωB97X, M06, and M06-2X functionals. Based on the explicit evaluation of anharmonic terms contributing to the second hyperpolarizability, this failure is traced down to a poor determination of third- and fourth-order energy derivatives with respect to normal modes. These results reveal serious flaws of some density functional approximations and suggest caution in selecting the appropriate functional to calculate not only electronic and vibrational (hyper)polarizabilities but also other molecular properties that contain vibrational anharmonic contributions.