Influence of nuclear quantum effects on the electronic properties of amorphous carbon.

We carry out quantum simulations to study the physical properties of diamond-like amorphous carbon by coupling first-principles molecular dynamics with a quantum thermostat, and we analyze multiple samples representative of different defective sites present in the disordered network. We show that quantum vibronic coupling is critical in determining the electronic properties of the system, in particular its electronic and mobility gaps, while it has a moderate influence on the structural properties. We find that despite localized electronic states near the Fermi level, the quantum nature of the nuclear motion leads to a renormalization of the electronic gap surprisingly similar to that found in crystalline diamond. We also discuss the notable influence of nuclear quantum effects on band-like and variable-hopping mechanisms contributing to electrical conduction. Our calculations indicate that methods often used to evaluate electron-phonon coupling in ordered solids are inaccurate to study the electronic and transport properties of amorphous semiconductors composed of light atoms.

Quantum simulations of thermally activated delayed fluorescence in an all-organic emitter.

We investigate the prototypical NAI-DMAC thermally activated delayed fluorescence (TADF) emitter in the gas phase- and high-packing fraction limits at finite temperature, by combining first principles molecular dynamics with a quantum thermostat to account for nuclear quantum effects (NQE). We find a weak dependence of the singlet-triplet energy gap (ΔEST) on temperature in both the solid and the molecule, and a substantial effect of packing. While the ΔEST vanishes in the perfect crystal, it is of the order of ∼0.3 eV in the molecule, with fluctuations ranging from 0.1 to 0.4 eV at 300 K. The transition probability between the HOMOs and LUMOs has a stronger dependence on temperature than the singlet-triplet gap, with a desirable effect for thermally activated fluorescence; such temperature effect is weaker in the condensed phase than in the molecule. Our results on ΔEST and oscillator strengths, together with our estimates of direct and reverse intersystem crossing rates, show that optimization of packing and geometrical conformation is critical to increase the efficiency of TADF compounds. Our findings highlight the importance of considering thermal fluctuations and NQE to obtain robust predictions of the electronic properties of NAI-DMAC.

Ab Initio Prediction of Adsorption Isotherms for Small Molecules in Metal-Organic Frameworks.

For CO and N2 on Mg2+ sites of the metal-organic framework CPO-27-Mg (Mg-MOF-74), ab initio calculations of Gibbs free energies of adsorption have been performed. Combined with the Bragg-Williams/Langmuir model and taking into account the experimental site availability (76.5%), we obtained adsorption isotherms in close agreement with those in experiment. The remaining deviations in the Gibbs free energy (about 1 kJ/mol) are significantly smaller than the "chemical accuracy" limit of about 4 kJ/mol. The presented approach uses (i) a DFT dispersion method (PBE+D2) to optimize the structure and to calculate anharmonic frequencies for vibrational partition functions and (ii) a "hybrid MP2:(PBE+D2)+ΔCCSD(T)" method to determine electronic energies. With the achieved accuracy (estimated uncertainty ±1.4 kJ/mol), the ab initio energies become useful benchmarks for assessing different DFT + dispersion methods (PBE+D2, B3LYP+D*, and vdW-D2), whereas the ab initio heats, entropies, and Gibbs free energies of adsorption are used to assess the reliability of experimental values derived from fitting isotherms or from variable-temperature IR studies.

Electroreflectance imaging of gold-H3PO4 supercapacitors. Part I: experimental methodology.

Electroreflectance microscopy is demonstrated as a high-resolution, non-contact method to image dynamic charge distribution in integrated microsupercapacitor structures during fast voltage cycling. Electroreflectance camera images of a gold electrode H3PO4 polymer electrolyte microsupercapacitor reveal time varying charge distribution with submicron spatial resolution, millisecond time resolution, and electroreflectance resolution on the order of 500 nC cm(-2). A model describing changes in the metal electrode's optical constants as a function of free electron concentration shows good agreement with measured electroreflectance. The proposed method can be used for sensitive, non-contact measurements of charge spatial distribution, and defect and performance characterization in electrode-electrolyte microdevices.

Electroreflectance imaging of gold-H3PO4 supercapacitors. Part II: microsupercapacitor ageing characterization.

This microsupercapacitor ageing study demonstrates the usefulness of the electroreflectance technique by quantifying local charge accumulation. Two separate devices with interdigitated electrodes were evaulated over a period of 4.1 million charge/discharge cycles. The key results are spatial mapping of charge accumulation in the gold electrodes derived from variation in the observed electrode reflectance. The nominal device exhibited little change in spatial distribution throughout the ageing cycle and serves as a comparison for the test device, which exhibited some nonuniform charge accumulation behavior. Further, an accelerated ageing test was completed by applying increasing voltage pulses up to 1.46 V to the device. Visual evidence of electrode ageing emerged in the reflectance distribution. An equivalent circuit model was developed to assess the evolution of individual circuit elements that correlate to the physical causes of ageing.

Quantum Vibronic Effects on the Excitation Energies of the Nitrogen-Vacancy Center in Diamond.

We investigated the impact of quantum vibronic coupling on the electronic properties of solid-state spin defects using stochastic methods and first-principles molecular dynamics with a quantum thermostat. Focusing on the negatively charged nitrogen-vacancy center in diamond as an exemplary case, we found a significant dynamic Jahn-Teller splitting of the doubly degenerate single-particle levels within the diamond's band gap, even at 0 K, with a magnitude exceeding 180 meV. This pronounced splitting leads to substantial renormalizations of these levels and, subsequently, of the vertical excitation energies of the doubly degenerate singlet and triplet excited states. Our findings underscore the pressing need to incorporate quantum vibronic effects into first-principles calculations, particularly when comparing computed vertical excitation energies with experimental data. Our study also reveals the efficiency of stochastic thermal line sampling for studying phonon renormalizations of solid-state spin defects.

Nuclear Quantum Effects on the Electronic Structure of Water and Ice.

The electronic properties and optical response of ice and water are intricately shaped by their molecular structure, including the quantum mechanical nature of the hydrogen atoms. Despite numerous previous studies, a comprehensive understanding of the nuclear quantum effects (NQEs) on the electronic structure of water and ice at finite temperatures remains elusive. Here, we utilize molecular simulations that harness efficient machine-learning potentials and many-body perturbation theory to assess how NQEs impact the electronic bands of water and hexagonal ice. By comparing path-integral and classical simulations, we find that NQEs lead to a larger renormalization of the fundamental gap of ice, compared to that of water, ultimately yielding similar bandgaps in the two systems, consistent with experimental estimates. Our calculations suggest that the increased quantum mechanical delocalization of protons in ice, relative to water, is a key factor leading to the enhancement of NQEs on the electronic structure of ice.

Quantum Vibronic Effects on the Electronic Properties of Molecular Crystals.

We present a study of molecular crystals, focused on the effect of nuclear quantum motion and anharmonicity on their electronic properties. We consider a system composed of relatively rigid molecules, a diamondoid crystal, and one composed of floppier molecules, NAI-DMAC, a thermally activated delayed fluorescence compound. We compute fundamental electronic gaps at the density functional theory (DFT) level of theory, with the Perdew-Burke-Erzenhof (PBE) and strongly constrained and approximately normed (SCAN) functionals, by coupling first-principles molecular dynamics with a nuclear quantum thermostat. We find a sizable zero-point renormalization (ZPR) of the band gaps, which is much larger in the case of diamondoids (0.6 eV) than for NAI-DMAC (0.22 eV). We show that the frozen phonon (FP) approximation, which neglects intermolecular anharmonic effects, leads to a large error (∼50%) in the calculation of the band gap ZPR. Instead, when using a stochastic method, we obtain results in good agreement with those of our quantum simulations for the diamondoid crystal. However, the agreement is worse for NAI-DMAC where intramolecular anharmonicities contribute to the ZPR. Our results highlight the importance of accurately including nuclear and anharmonic quantum effects to predict the electronic properties of molecular crystals.

Computational Protocol to Evaluate Electron-Phonon Interactions Within Density Matrix Perturbation Theory.

We present a computational protocol, based on density matrix perturbation theory, to obtain non-adiabatic, frequency-dependent electron-phonon self-energies for molecules and solids. Our approach enables the evaluation of electron-phonon interaction using hybrid functionals, for spin-polarized systems, and the computational overhead to include dynamical and non-adiabatic terms in the evaluation of electron-phonon self-energies is negligible. We discuss results for molecules, as well as pristine and defective solids.

Combined First-Principles Calculations of Electron-Electron and Electron-Phonon Self-Energies in Condensed Systems.

We present a method to efficiently combine the computation of electron-electron and electron-phonon self-energies, which enables the evaluation of electron-phonon coupling at the G0W0 level of theory for systems with hundreds of atoms. In addition, our approach, which is a generalization of a method recently proposed for molecules [J. Chem. Theory Comput. 2018, 14, 6269-6275], enables the inclusion of nonadiabatic and temperature effects at no additional computational cost. We present results for diamond and defects in diamond and discuss the importance of numerically accurate G0W0 band structures to obtain robust predictions of zero point renormalization (ZPR) of band gaps, and of the inclusion of nonadiabatic effects to accurately compute the ZPR of defect states in the band gap.

Predicting adsorption selectivities from pure gas isotherms for gas mixtures in metal-organic frameworks.

We perform Grand Canonical Monte Carlo simulations on a lattice of Mg2+ sites (GCMC) for adsorption of four binary A/B mixtures, CH4/N2, CO/N2, CO2/N2, and CO2/CH4, in the metal-organic framework Mg2(2,5-dioxidobenzedicarboxylate), also known as CPO-27-Mg or Mg-MOF-74. We present a mean field co-adsorption isotherm model and show that its predictions agree with the GCMC results if the same quantum chemical ab initio data are used for Gibbs free energies of adsorption at the individual sites and for lateral interaction energies between the same, A⋯A and B⋯B, and unlike, A⋯B, adsorbed molecules. We use both approaches to test the assumption underlying Ideal Adsorbed Solution Theory (IAST), namely approximating A⋯B interaction energies as the arithmetic mean of A⋯A and B⋯B interaction energies. While IAST works well for mixtures with weak lateral interactions, CH4/N2 and CO/N2, the deviations are large for mixtures with stronger lateral interactions, CO2/N2 and CO2/CH4. Motivated by the theory of London dispersion forces, we propose use of the geometric mean instead of the arithmetic mean and achieve substantial improvements. For CO2/CH4, the lateral interactions become anisotropic. To include this in the geometric mean co-adsorption model, we introduce an anisotropy factor. We propose a protocol, named co-adsorption mean field theory (CAMT), for co-adsorption selectivity prediction from known (experiment or simulation) pure component isotherms which is similar to the IAST protocol but uses the geometric mean to approximate mixed pair interaction energies and yields improved results for non-ideal mixtures.

Ab Initio Prediction of Adsorption Isotherms for Gas Mixtures by Grand Canonical Monte Carlo Simulations on a Lattice of Sites.

Gibbs free energies of adsorption on individual sites and the lateral (adsorbate-adsorbate) interaction energies are obtained from quantum chemical ab initio methods and molecular statistics. They define a Grand Canonical Monte Carlo (GCMC) Hamiltonian for simulations of gas mixtures on a lattice of adsorption sites. Coadsorption of CO2 and CH4 at Mg2+ sites in the pores of the metal-organic framework CPO-27-Mg (Mg-MOF-74) is studied as an example. Simulations with different approximations as made in widely used coadsorption models such as the ideal adsorbed solution theory (IAST) show their limitations in describing adsorption selectivities for binary mixtures.