skip to main content

Title: Self-interaction error overbinds water clusters but cancels in structural energy differences

We gauge the importance of self-interaction errors in density functional approximations (DFAs) for the case of water clusters. To this end, we used the Fermi–Löwdin orbital self-interaction correction method (FLOSIC) to calculate the binding energy of clusters of up to eight water molecules. Three representative DFAs of the local, generalized gradient, and metageneralized gradient families [i.e., local density approximation (LDA), Perdew–Burke–Ernzerhof (PBE), and strongly constrained and appropriately normed (SCAN)] were used. We find that the overbinding of the water clusters in these approximations is not a density-driven error. We show that, while removing self-interaction error does not alter the energetic ordering of the different water isomers with respect to the uncorrected DFAs, the resulting binding energies are corrected toward accurate reference values from higher-level calculations. In particular, self-interaction–corrected SCAN not only retains the correct energetic ordering for water hexamers but also reduces the mean error in the hexamer binding energies to less than 14 meV/H2Ofrom about 42 meV/H2Ofor SCAN. By decomposing the total binding energy into many-body components, we find that large errors in the two-body interaction in SCAN are significantly reduced by self-interaction corrections. Higher-order many-body errors are small in both SCAN and more » self-interaction–corrected SCAN. These results indicate that orbital-by-orbital removal of self-interaction combined with a proper DFA can lead to improved descriptions of water complexes.

« less
; ; ; ; ; ; ; ;
Publication Date:
Journal Name:
Proceedings of the National Academy of Sciences
Page Range or eLocation-ID:
p. 11283-11288
Proceedings of the National Academy of Sciences
Sponsoring Org:
National Science Foundation
More Like this
  1. The product selectivity of many heterogeneous electrocatalytic processes is profoundly affected by the liquid side of the electrocatalytic interface. The electrocatalytic reduction of CO to hydrocarbons on Cu electrodes is a prototypical example of such a process. However, probing the interactions of surface-bound intermediates with their liquid reaction environment poses a formidable experimental challenge. As a result, the molecular origins of the dependence of the product selectivity on the characteristics of the electrolyte are still poorly understood. Herein, we examined the chemical and electrostatic interactions of surface-adsorbed CO with its liquid reaction environment. Using a series of quaternary alkyl ammonium cations (methyl4N+,ethyl4N+,propyl4N+, andbutyl4N+), we systematically tuned the properties of this environment. With differential electrochemical mass spectrometry (DEMS), we show that ethylene is produced in the presence ofmethyl4N+andethyl4N+cations, whereas this product is not synthesized inpropyl4N+- andbutyl4N+-containing electrolytes. Surface-enhanced infrared absorption spectroscopymore »(SEIRAS) reveals that the cations do not block CO adsorption sites and that the cation-dependent interfacial electric field is too small to account for the observed changes in selectivity. However, SEIRAS shows that an intermolecular interaction between surface-adsorbed CO and interfacial water is disrupted in the presence of the two larger cations. This observation suggests that this interaction promotes the hydrogenation of surface-bound CO to ethylene. Our study provides a critical molecular-level insight into how interactions of surface species with the liquid reaction environment control the selectivity of this complex electrocatalytic process.

    « less
  2. Abstract

    We present a Keck/MOSFIRE rest-optical composite spectrum of 16 typical gravitationally lensed star-forming dwarf galaxies at 1.7 ≲z≲ 2.6 (zmean= 2.30), all chosen independent of emission-line strength. These galaxies have a median stellar mass oflog(M*/M)med=8.290.43+0.51and a median star formation rate ofSFRHαmed=2.251.26+2.15Myr1. We measure the faint electron-temperature-sensitive [Oiii]λ4363 emission line at 2.5σ(4.1σ) significance when considering a bootstrapped (statistical-only) uncertainty spectrum. This yields a direct-method oxygen abundance of12+log(O/H)direct=7.880.22+0.25(0.150.06+0.12Z). We investigate the applicability at highzof locally calibrated oxygen-based strong-line metallicity relations, finding that the local reference calibrations of Bian et al. best reproduce (≲0.12 dex) our composite metallicity at fixed strong-line ratio. At fixedM*, our composite is well represented by thez∼ 2.3 direct-method stellar mass—gas-phase metallicity relation (MZR) of Sanders et al. When comparing to predicted MZRs from the IllustrisTNG and FIRE simulations, having recalculated our stellar masses with more realistic nonparametric star formation histories(log(M*/M)med=8.920.22+0.31), we find excellent agreement with the FIRE MZR. Our composite is consistent with no metallicity evolution, atmore »fixedM*and SFR, of the locally defined fundamental metallicity relation. We measure the doublet ratio [Oii]λ3729/[Oii]λ3726 = 1.56 ± 0.32 (1.51 ± 0.12) and a corresponding electron density ofne=10+215cm3(ne=10+74cm3) when considering the bootstrapped (statistical-only) error spectrum. This result suggests that lower-mass galaxies have lower densities than higher-mass galaxies atz∼ 2.

    « less
  3. Sub-Neptunes are common among the discovered exoplanets. However, lack of knowledge on the state of matter inH2O-rich setting at high pressures and temperatures (PT) places important limitations on our understanding of this planet type. We have conducted experiments for reactions betweenSiO2andH2O as archetypal materials for rock and ice, respectively, at highPT. We found anomalously expanded volumes of dense silica (up to 4%) recovered from hydrothermal synthesis above ∼24 GPa where theCaCl2-type (Ct) structure appears at lower pressures than in the anhydrous system. Infrared spectroscopy identified strong OH modes from the dense silica samples. Both previous experiments and our density functional theory calculations support up to 0.48 hydrogen atoms per formula unit of (Si1xH4x)O2(x=0.12). At pressures above 60 GPa,H2O further changes the structural behavior of silica, stabilizing a niccolite-type structure, which is unquenchable. From unit-cell volume and phase equilibrium considerations, we infer that the niccolite-type phase may contain H with an amount at least comparable with or higher than that of the Ct phase. Our results suggest that the phases containing both hydrogen and lithophile elements could bemore »the dominant materials in the interiors of water-rich planets. Even for fully layered cases, the large mutual solubility could make the boundary between rock and ice layers fuzzy. Therefore, the physical properties of the new phases that we report here would be important for understanding dynamics, geochemical cycle, and dynamo generation in water-rich planets.

    « less
  4. Abstract

    We report the temperature dependence of the Yb valence in the geometrically frustrated compoundYbB4from 12 to 300 K using resonant x-ray emission spectroscopy at the YbLα1transition. We find that the Yb valence,v, is hybridized between thev = 2 andv = 3 valence states, increasing fromv=2.61±0.01at 12 K tov=2.67±0.01at 300 K, confirming thatYbB4is a Kondo system in the intermediate valence regime. This result indicates that the Kondo interaction inYbB4is substantial, and is likely to be the reason whyYbB4does not order magnetically at low temperature, rather than this being an effect of geometric frustration. Furthermore, the zero-point valence of the system is extracted from our data and compared with other Kondo lattice systems. The zero-point valence seems to be weakly dependent on the Kondo temperature scale, but not on the valence change temperature scaleTv.

  5. Abstract

    We present a multiwavelength analysis of the galaxy cluster SPT-CL J0607-4448 (SPT0607), which is one of the most distant clusters discovered by the South Pole Telescope atz= 1.4010 ± 0.0028. The high-redshift cluster shows clear signs of being relaxed with well-regulated feedback from the active galactic nucleus (AGN) in the brightest cluster galaxy (BCG). Using Chandra X-ray data, we construct thermodynamic profiles and determine the properties of the intracluster medium. The cool-core nature of the cluster is supported by a centrally peaked density profile and low central entropy (K0=189+11keV cm2), which we estimate assuming an isothermal temperature profile due to the limited spectral information given the distance to the cluster. Using the density profile and gas cooling time inferred from the X-ray data, we find a mass-cooling rateṀcool=10060+90Myr−1. From optical spectroscopy and photometry around the [Oii] emission line, we estimate that the BCG star formation rate isSFR[OII]=1.70.6+1.0Myr−1, roughly two orders of magnitude lower than the predicted mass-cooling rate. In addition, using ATCA radio data at 2.1 GHz, we measure a radio jet powerPcav=3.21.3+2.1×1044erg s−1, which is consistent withmore »the X-ray cooling luminosity (Lcool=1.90.5+0.2×1044erg s−1withinrcool= 43 kpc). These findings suggest that SPT0607 is a relaxed, cool-core cluster with AGN-regulated cooling at an epoch shortly after cluster formation, implying that the balance between cooling and feedback can be reached quickly. We discuss the implications for these findings on the evolution of AGN feedback in galaxy clusters.

    « less