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 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.

 
more » « less
NSF-PAR ID:
10149661
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Publisher / Repository:
Proceedings of the National Academy of Sciences
Date Published:
Journal Name:
Proceedings of the National Academy of Sciences
Volume:
117
Issue:
21
ISSN:
0027-8424
Page Range / eLocation ID:
p. 11283-11288
Format(s):
Medium: X
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 spectroscopy (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.

     
    more » « less
  2. Abstract

    I employ the Lucy rectification algorithm to recover the inclination-corrected distribution of local disk galaxies in the plane of absolute magnitude (Mi) and Hivelocity width (W20). By considering the inclination angle as a random variable with a known probability distribution, the novel approach eliminates one major source of uncertainty in studies of the Tully–Fisher relation: inclination angle estimation from axial ratio. Leveraging the statistical strength derived from the entire sample of 28,264 Hi-selected disk galaxies atz< 0.06 from the Arecibo Legacy Fast ALFA survey, I show that the restored distribution follows a sharp correlation that is approximately a power law between −16 >Mi> −22:Mi=M02.5β[log(W20/250km/s)], withM0= −19.77± 0.04 andβ= 4.39 ± 0.06. At the brighter end (Mi< −22), the slope of the correlation decreases toβ≈ 3.3, confirming previous results. Because the method accounts for measurement errors, the intrinsic dispersion of the correlation is directly measured:σ(logW20)0.06dex between −17 >Mi> −23, whileσ(Mi) decreases from ∼0.8 in slow rotators to ∼0.4 in fast rotators. The statistical rectification method holds significant potential, especially in the studies of intermediate-to-high-redshift samples, where limited spatial resolution hinders precise measurements of inclination angles.

     
    more » « less
  3. The experimental sum frequency generation (SFG) spectrum is the response to an infrared pulse and a visible pulse and is a highly surface-sensitive technique. We treat the surface dangling OH bonds at the air/water interface and focus on the absolute SFG intensities for the resonant terms, a focus that permits insight into the consequences of some approximations. For the polarization combinations, the calculated linewidths for the water interface dangling OH SFG band at 3,700cm1are, as usual, too large, because of the customary neglect of motional narrowing. The integrated spectrum is used to circumvent this problem and justified here using a Kubo-like formalism and theoretical integrated band intensities rather than peak intensities. Only relative SFG intensities are usually reported. The absolute integrated SFG intensities for three polarization combinations for sum frequency, visible, and infrared beams are computed. We use molecular dynamics and the dipole and the polarizability matrix elements obtained from infrared and Raman studies ofH2O vapor. The theoretical expressions for two of the absolute susceptibilities contain only a single term and agree with experiment to about a factor of 1.3, with no adjustable parameters. The Fresnel factors are included in that comparison. One of the susceptibilities contains instead four positive and negative terms and agrees less well. The expression for the SFG correlation function is normally derived from a statistical mechanical formulation using a time-evolving density matrix. We show how a derivation based on a two-field relaxation leads to the same final result.

     
    more » « less
  4. Abstract

    M dwarfs are common host stars to exoplanets but often lack atmospheric abundance measurements. Late-M dwarfs are also good analogs to the youngest substellar companions, which share similarTeff∼ 2300–2800 K. We present atmospheric analyses for the M7.5 companion HIP 55507 B and its K6V primary star with Keck/KPIC high-resolution (R∼ 35,000)K-band spectroscopy. First, by including KPIC relative radial velocities between the primary and secondary in the orbit fit, we improve the dynamical mass precision by 60% and findMB=88.03.2+3.4MJup, putting HIP 55507 B above the stellar–substellar boundary. We also find that HIP 55507 B orbits its K6V primary star witha=383+4au ande= 0.40 ± 0.04. From atmospheric retrievals of HIP 55507 B, we measure [C/H] = 0.24 ± 0.13, [O/H] = 0.15 ± 0.13, and C/O = 0.67 ± 0.04. Moreover, we strongly detect13CO (7.8σsignificance) and tentatively detectH218O(3.7σsignificance) in the companion’s atmosphere and measure12CO/13CO=9822+28andH216O/H218O=24080+145after accounting for systematic errors. From a simplified retrieval analysis of HIP 55507 A, we measure12CO/13CO=7916+21andC16O/C18O=28870+125for the primary star. These results demonstrate that HIP 55507 A and B have consistent12C/13C and16O/18O to the <1σlevel, as expected for a chemically homogeneous binary system. Given the similar flux ratios and separations between HIP 55507 AB and systems with young substellar companions, our results open the door to systematically measuring13CO andH218Oabundances in the atmospheres of substellar or even planetary-mass companions with similar spectral types.

     
    more » « less
  5. 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, at 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.

     
    more » « less