Averaging over atom snapshots in linear-response TDDFT of disordered systems: A case study of warm dense hydrogen.

Linear-response time-dependent density functional theory (LR-TDDFT) simulations of disordered extended systems require averaging over different snapshots of ion configurations to minimize finite size effects due to the snapshot-dependence of the electronic density response function and related properties. We present a consistent scheme for the computation of the macroscopic Kohn-Sham (KS) density response function connecting an average over snapshot values of charge density perturbations to the averaged values of KS potential variations. This allows us to formulate the LR-TDDFT within the adiabatic (static) approximation for the exchange-correlation (XC) kernel for disordered systems, where the static XC kernel is computed using the direct perturbation method [Moldabekov et al., J. Chem. Theory Comput. 19, 1286 (2023)]. The presented approach allows one to compute the macroscopic dynamic density response function as well as the dielectric function with a static XC kernel generated for any available XC functional. The application of the developed workflow is demonstrated for the example of warm dense hydrogen. The presented approach is applicable for various types of extended disordered systems, such as warm dense matter, liquid metals, and dense plasmas.

Stretch Evolution of Electronic Coupling of the Thiophenyl Anchoring Group with Gold in Mechanically Controllable Break Junctions.

The current-voltage characteristics of a single-molecule junction are determined by the electronic coupling Γ between the electronic states of the electrodes and the dominant transport channel(s) of the molecule. Γ is profoundly affected by the choice of the anchoring groups and their binding positions on the tip facets and the tip-tip separation. In this work, mechanically controllable break junction experiments on the N,N'-bis(5-ethynylbenzenethiol-salicylidene)ethylenediamine are presented, in particular, the stretch evolution of Γ with increasing tip-tip separation. The stretch evolution of Γ is characterized by recurring local maxima and can be related to the deformation of the molecule and sliding of the anchoring groups above the tip facets and along the tip edges. A dynamic simulation approach is implemented to model the stretch evolution of Γ, which captures the experimentally observed features remarkably well and establishes a link to the microscopic structure of the single-molecule junction.

Assessing the accuracy of hybrid exchange-correlation functionals for the density response of warm dense electrons.

We assess the accuracy of common hybrid exchange-correlation (XC) functionals (PBE0, PBE0-1/3, HSE06, HSE03, and B3LYP) within the Kohn-Sham density functional theory for the harmonically perturbed electron gas at parameters relevant for the challenging conditions of the warm dense matter. Generated by laser-induced compression and heating in the laboratory, the warm dense matter is a state of matter that also occurs in white dwarfs and planetary interiors. We consider both weak and strong degrees of density inhomogeneity induced by the external field at various wavenumbers. We perform an error analysis by comparing with the exact quantum Monte Carlo results. In the case of a weak perturbation, we report the static linear density response function and the static XC kernel at a metallic density for both the degenerate ground-state limit and for partial degeneracy at the electronic Fermi temperature. Overall, we observe an improvement in the density response when the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are used, compared with the previously reported results for the PBE, PBEsol, local-density approximation, and AM05 functionals; B3LYP, on the other hand, does not perform well for the considered system. Additionally, the PBE0, PBE0-1/3, HSE06, and HSE03 functionals are more accurate for the density response properties than SCAN in the regime of partial degeneracy.

Non-empirical Mixing Coefficient for Hybrid XC Functionals from Analysis of the XC Kernel.

We present an analysis of the static exchange-correlation (XC) kernel computed from hybrid functionals with a single mixing coefficient such as PBE0 and PBE0-1/3. We break down the hybrid XC kernels into the exchange and correlation parts using the Hartree-Fock functional, the exchange-only PBE, and the correlation-only PBE. This decomposition is combined with exact data for the static XC kernel of the uniform electron gas and an Airy gas model within a subsystem functional approach. This gives us a tool for the non-empirical choice of the mixing coefficient under ambient and extreme conditions. Our analysis provides physical insights into the effect of the variation of the mixing coefficient in hybrid functionals, which is of immense practical value. The presented approach is general and can be used for other types of functionals like screened hybrids.