skip to main content

Title: Absorption spectra of benzoic acid in water at different pH and in the presence of salts: insights from the integration of experimental data and theoretical cluster models
The absorption spectra of molecular organic chromophores in aqueous media are of considerable importance in environmental chemistry. In this work, the UV-vis spectra of benzoic acid (BA), the simplest aromatic carboxylic acid, in aqueous solutions at varying pH and in the presence of salts are measured experimentally. The solutions of different pH provide insights into the contributions from both the non-dissociated acid molecule and the deprotonated anionic species. The microscopic interpretation of these spectra is then provided by quantum chemical calculations for small cluster models of benzoic species (benzoic acid and benzoate anion) with water molecules. Calculations of the UV-vis absorbance spectra are then carried out for different clusters such as C 6 H 5 COOH·(H 2 O) n and C 6 H 5 COO − ·(H 2 O) n , where n = 0–8. The following main conclusions from these calculations and the comparison to experimental results can be made: (i) the small water cluster yields good quantitative agreement with observed solution experiments; (ii) the main peak position is found to be very similar at different levels of theory and is in excellent agreement with the experimental value, however, a weaker feature about 1 eV to lower energy (red more » shift) of the main peak is correctly reproduced only by using high level of theory, such as Algebraic Diagrammatic Construction (ADC); (iii) dissociation of the BA into ions is found to occur with a minimum of water molecules of n = 8; (iv) the deprotonation of BA has an influence on the computed spectrum and the energetics of the lowest energy electronic transitions; (v) the effect of the water on the spectra is much larger for the deprotonated species than for the non-dissociated acid. It was found that to reproduce experimental spectrum at pH 8.0, additional continuum representation for the extended solvent environment must be included in combination with explicit solvent molecules ( n ≥ 3); (vi) salts (NaCl and CaCl 2 ) have minimal effect on the absorption spectrum and; (vii) experimental results showed that B-band of neutral BA is not sensitive to the solvent effects whereas the effect of the water on the C-band is significant. The water effects blue-shift this band up to ∼0.2 eV. Overall, the results demonstrate the ability to further our understanding of the microscopic interpretation of the electronic structure and absorption spectra of BA in aqueous media through calculations restricted to small cluster models. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
5046 to 5056
Sponsoring Org:
National Science Foundation
More Like this
  1. Resonant tunneling diodes (RTDs) have come full-circle in the past 10 years after their demonstration in the early 1990s as the fastest room-temperature semiconductor oscillator, displaying experimental results up to 712 GHz and fmax values exceeding 1.0 THz [1]. Now the RTD is once again the preeminent electronic oscillator above 1.0 THz and is being implemented as a coherent source [2] and a self-oscillating mixer [3], amongst other applications. This paper concerns RTD electroluminescence – an effect that has been studied very little in the past 30+ years of RTD development, and not at room temperature. We present experiments andmore »modeling of an n-type In0.53Ga0.47As/AlAs double-barrier RTD operating as a cross-gap light emitter at ~300K. The MBE-growth stack is shown in Fig. 1(a). A 15-μm-diam-mesa device was defined by standard planar processing including a top annular ohmic contact with a 5-μm-diam pinhole in the center to couple out enough of the internal emission for accurate free-space power measurements [4]. The emission spectra have the behavior displayed in Fig. 1(b), parameterized by bias voltage (VB). The long wavelength emission edge is at  = 1684 nm - close to the In0.53Ga0.47As bandgap energy of Ug ≈ 0.75 eV at 300 K. The spectral peaks for VB = 2.8 and 3.0 V both occur around  = 1550 nm (h = 0.75 eV), so blue-shifted relative to the peak of the “ideal”, bulk InGaAs emission spectrum shown in Fig. 1(b) [5]. These results are consistent with the model displayed in Fig. 1(c), whereby the broad emission peak is attributed to the radiative recombination between electrons accumulated on the emitter side, and holes generated on the emitter side by interband tunneling with current density Jinter. The blue-shifted main peak is attributed to the quantum-size effect on the emitter side, which creates a radiative recombination rate RN,2 comparable to the band-edge cross-gap rate RN,1. Further support for this model is provided by the shorter wavelength and weaker emission peak shown in Fig. 1(b) around = 1148 nm. Our quantum mechanical calculations attribute this to radiative recombination RR,3 in the RTD quantum well between the electron ground-state level E1,e, and the hole level E1,h. To further test the model and estimate quantum efficiencies, we conducted optical power measurements using a large-area Ge photodiode located ≈3 mm away from the RTD pinhole, and having spectral response between 800 and 1800 nm with a peak responsivity of ≈0.85 A/W at  =1550 nm. Simultaneous I-V and L-V plots were obtained and are plotted in Fig. 2(a) with positive bias on the top contact (emitter on the bottom). The I-V curve displays a pronounced NDR region having a current peak-to-valley current ratio of 10.7 (typical for In0.53Ga0.47As RTDs). The external quantum efficiency (EQE) was calculated from EQE = e∙IP/(∙IE∙h) where IP is the photodiode dc current and IE the RTD current. The plot of EQE is shown in Fig. 2(b) where we see a very rapid rise with VB, but a maximum value (at VB= 3.0 V) of only ≈2×10-5. To extract the internal quantum efficiency (IQE), we use the expression EQE= c ∙i ∙r ≡ c∙IQE where ci, and r are the optical-coupling, electrical-injection, and radiative recombination efficiencies, respectively [6]. Our separate optical calculations yield c≈3.4×10-4 (limited primarily by the small pinhole) from which we obtain the curve of IQE plotted in Fig. 2(b) (right-hand scale). The maximum value of IQE (again at VB = 3.0 V) is 6.0%. From the implicit definition of IQE in terms of i and r given above, and the fact that the recombination efficiency in In0.53Ga0.47As is likely limited by Auger scattering, this result for IQE suggests that i might be significantly high. To estimate i, we have used the experimental total current of Fig. 2(a), the Kane two-band model of interband tunneling [7] computed in conjunction with a solution to Poisson’s equation across the entire structure, and a rate-equation model of Auger recombination on the emitter side [6] assuming a free-electron density of 2×1018 cm3. We focus on the high-bias regime above VB = 2.5 V of Fig. 2(a) where most of the interband tunneling should occur in the depletion region on the collector side [Jinter,2 in Fig. 1(c)]. And because of the high-quality of the InGaAs/AlAs heterostructure (very few traps or deep levels), most of the holes should reach the emitter side by some combination of drift, diffusion, and tunneling through the valence-band double barriers (Type-I offset) between InGaAs and AlAs. The computed interband current density Jinter is shown in Fig. 3(a) along with the total current density Jtot. At the maximum Jinter (at VB=3.0 V) of 7.4×102 A/cm2, we get i = Jinter/Jtot = 0.18, which is surprisingly high considering there is no p-type doping in the device. When combined with the Auger-limited r of 0.41 and c ≈ 3.4×10-4, we find a model value of IQE = 7.4% in good agreement with experiment. This leads to the model values for EQE plotted in Fig. 2(b) - also in good agreement with experiment. Finally, we address the high Jinter and consider a possible universal nature of the light-emission mechanism. Fig. 3(b) shows the tunneling probability T according to the Kane two-band model in the three materials, In0.53Ga0.47As, GaAs, and GaN, following our observation of a similar electroluminescence mechanism in GaN/AlN RTDs (due to strong polarization field of wurtzite structures) [8]. The expression is Tinter = (2/9)∙exp[(-2 ∙Ug 2 ∙me)/(2h∙P∙E)], where Ug is the bandgap energy, P is the valence-to-conduction-band momentum matrix element, and E is the electric field. Values for the highest calculated internal E fields for the InGaAs and GaN are also shown, indicating that Tinter in those structures approaches values of ~10-5. As shown, a GaAs RTD would require an internal field of ~6×105 V/cm, which is rarely realized in standard GaAs RTDs, perhaps explaining why there have been few if any reports of room-temperature electroluminescence in the GaAs devices. [1] E.R. Brown,et al., Appl. Phys. Lett., vol. 58, 2291, 1991. [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [2] M. Feiginov et al., Appl. Phys. Lett., 99, 233506, 2011. [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [3] Y. Nishida et al., Nature Sci. Reports, 9, 18125, 2019. [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [4] P. Fakhimi, et al., 2019 DRC Conference Digest. [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018). [5] S. Sze, Physics of Semiconductor Devices, 2nd Ed. 12.2.1 (Wiley, 1981). [6] L. Coldren, Diode Lasers and Photonic Integrated Circuits, (Wiley, 1995). [7] E.O. Kane, J. of Appl. Phy 32, 83 (1961). [8] T. Growden, et al., Nature Light: Science & Applications 7, 17150 (2018).« less
  2. We exploit gas-phase cluster ion techniques to provide insight into the local interactions underlying divalent metal ion-driven changes in the spectra of carboxylic acids at the air–water interface. This information clarifies the experimental findings that the CO stretching bands of long-chain acids appear at very similar energies when the head group is deprotonated by high subphase pH or exposed to relatively high concentrations of Ca 2+ metal ions. To this end, we report the evolution of the vibrational spectra of size-selected [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n =0 to 12 and RCO 2 − ·(Hmore »2 O) n =0 to 14 cluster ions toward the features observed at the air–water interface. Surprisingly, not only does stepwise hydration of the RCO 2 − anion and the [Ca 2+ ·RCO 2 − ] + contact ion pair yield solvatochromic responses in opposite directions, but in both cases, the responses of the 2 (symmetric and asymmetric stretching) CO bands to hydration are opposite to each other. The result is that both CO bands evolve toward their interfacial asymptotes from opposite directions. Simulations of the [Ca 2+ ·RCO 2 − ] + ·(H 2 O) n clusters indicate that the metal ion remains directly bound to the head group in a contact ion pair motif as the asymmetric CO stretch converges at the interfacial value by n = 12. This establishes that direct metal complexation or deprotonation can account for the interfacial behavior. We discuss these effects in the context of a model that invokes the water network-dependent local electric field along the C–C bond that connects the head group to the hydrocarbon tail as the key microscopic parameter that is correlated with the observed trends.« less
  3. The partitioning of medium-chain fatty acid surfactants such as nonanoic acid (NA) between the bulk phase and the air/water interface is of interest to a number of fields including marine and atmospheric chemistry. However, questions remain about the behavior of these molecules, the contributions of various relevant chemical equilibria, and the impact of pH, salt and bulk surfactant concentrations. In this study, the surface adsorption of nonanoic acid and its conjugate base is quantitatively investigated at various pH values, surfactant concentrations and the presence of salts. Surface concentrations of protonated and deprotonated species are dictated by surface-bulk equilibria which canmore »be calculated from thermodynamic considerations. Notably we conclude that the surface dissociation constant of soluble surfactants cannot be directly obtained from these experimental measurements, however, we show that molecular dynamics (MD) simulation methods, such as free energy perturbation (FEP), can be used to calculate the surface acid dissociation constant relative to that in the bulk. These simulations show that nonanoic acid is less acidic at the surface compared to in the bulk solution with a p K a shift of 1.1 ± 0.6, yielding a predicted surface p K a of 5.9 ± 0.6. A thermodynamic cycle for nonanoic acid and its conjugate base between the air/water interface and the bulk phase can therefore be established. Furthermore, the effect of salts, namely NaCl, on the surface activity of protonated and deprotonated forms of nonanoic acid is also examined. Interestingly, salts cause both a decrease in the bulk p K a of nonanoic acid and a stabilization of both the protonated and deprotonated forms at the surface. Overall, these results suggest that the deprotonated medium-chain fatty acids under ocean conditions can also be present within the sea surface microlayer (SSML) present at the ocean/atmosphere interface due to the stabilization effect of the salts in the ocean. This allows the transfer of these species into sea spray aerosols (SSAs). More generally, we present a framework with which the behavior of partially soluble species at the air/water interface can be predicted from surface adsorption models and the surface p K a can be predicted from MD simulations.« less
  4. Cyanuric acid is a triazine derivative that has been identified from reactions performed under prebiotic conditions and has been proposed as a prospective precursor of ancestral RNA. For cyanuric acid to have played a key role during the prebiotic era, it would have needed to survive the harsh electromagnetic radiation conditions reaching the Earth’s surface during prebiotic times (≥200 nm). Therefore, the photostability of cyanuric acid would have been crucial for its accumulation during the prebiotic era. To evaluate the putative photostability of cyanuric acid in water, in this contribution, we employed density functional theory (DFT) and its time-dependent variantmore »(TD-DFT) including implicit and explicit solvent effects. The calculations predict that cyanuric acid has an absorption maximum at ca. 160 nm (7.73 eV), with the lowest-energy absorption band extending to ca. 200 nm in an aqueous solution and exhibiting negligible absorption at longer wavelengths. Excitation of cyanuric acid at 160 nm or longer wavelengths leads to the population of S5,6 singlet states, which have ππ* character and large oscillator strengths (0.8). The population reaching the S5,6 states is expected to internally convert to the S1,2 states in an ultrafast time scale. The S1,2 states, which have nπ* character, are predicted to access a conical intersection with the ground state in a nearly barrierless fashion (ca. ≤ 0.13 eV), thus efficiently returning the population to the ground state. Furthermore, based on calculated spin–orbit coupling elements of ca. 6 to 8 cm−1, the calculations predict that intersystem crossing to the triplet manifold should play a minor role in the electronic relaxation of cyanuric acid. We have also calculated the vertical ionization energy of cyanuric acid at 8.2 eV, which predicts that direct one-photon ionization of cyanuric acid should occur at ca. 150 nm. Collectively, the quantum-chemical calculations predict that cyanuric acid would have been highly photostable under the solar radiation conditions reaching the Earth’s surface during the prebiotic era in an aqueous solution. Of relevance to the chemical origin of life and RNA-first theories, these observations lend support to the idea that cyanuric acid could have accumulated in large quantities during the prebiotic era and thus strengthens its candidature as a relevant prebiotic nucleobase.« less
  5. A series of complexes with low-energy Fe II to Ti IV metal-to-metal charge-transfer (MMCT) transitions, Cp 2 Ti(C 2 Fc) 2 , Cp* 2 Ti(C 2 Fc) 2 , and MeOOC Cp 2 Ti(C 2 Fc) 2 , was investigated using solvatochromism and resonance Raman spectroscopy (RRS) augmented with time-dependent density functional theory (TDDFT) calculations in order to interrogate the nature of the CT transitions. Computational models were benchmarked against the experimental UV-Vis spectra and B3LYP/6-31G(d) was found to most faithfully represent the spectra. The energy of the MMCT transition was measured in 15 different solvents and a multivariate fitmore »to the Catalán solvent parameters – solvent polarizability (SP), solvent dipolarity (SdP), solvent basicity (SB), and solvent acidity (SA) – was performed. The effect of SP indicates a greater degree of electron delocalization in the excited state (ES) than the ground state (GS). The small negative solvatochromism with respect to SdP indicates a smaller dipole moment in the ES than the GS. The effect of SB is consistent with charge-transfer to Ti. Upon excitation into the MMCT absorption band, the RRS data show enhancement of the alkyne stretching modes and of the out-of-plane bending modes of the cyclopentadienyl ring connected to Fe and the alkyne bridge. This is consistent with changes in the oxidation states of Ti and Fe, respectively. The higher-energy transitions (350–450 nm) show enhancement of vibrational modes consistent with ethnylcyclopentadienyl to Ti ligand-to-metal charge transfer (LMCT). The RRS data is consistent with the TDDFT predicted character of these transitions. TDDFT suggests that the lowest-energy transition in Cp 2 Ti(C 2 Fc) 2 CuI, where CuI is coordinated between the alkynes, retains its Fe II to Ti IV MMCT character, in agreement with the RRS data, but that the lowest-energy transitions have significant CuI to Ti character. For Cp 2 Ti(C 2 Fc) 2 CuI, excitation into the low-energy MMCT absorption band results in selective enhancement of the symmetric alkynyl stretching mode.« less