skip to main content

Title: Ab initio metadynamics calculations of dimethylamine for probing p K b variations in bulk vs. surface environments
The basicity constant, or p K b , is an intrinsic physical property of bases that gives a measure of its proton affinity in macroscopic environments. While the p K b is typically defined in reference to the bulk aqueous phase, several studies have suggested that this value can differ significantly at the air–water interface (which can have significant ramifications for particle surface chemistry and aerosol growth modeling). To provide mechanistic insight into surface proton affinity, we carried out ab initio metadynamics calculations to (1) explore the free-energy profile of dimethylamine and (2) provide reasonable estimates of the p K b value in different solvent environments. We find that the free-energy profiles obtained with our metadynamics calculations show a dramatic variation, with interfacial aqueous dimethylamine p K b values being significantly lower than in the bulk aqueous environment. Furthermore, our metadynamics calculations indicate that these variations are due to reduced hydrogen bonding at the air–water surface. Taken together, our quantum mechanical metadynamics calculations show that the reactivity of dimethylamine is surprisingly complex, leading to p K b variations that critically depend on the different atomic interactions occurring at the microscopic molecular level.
; ; ; ; ;
Award ID(s):
Publication Date:
Journal Name:
Physical Chemistry Chemical Physics
Page Range or eLocation-ID:
26265 to 26277
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 and 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.more »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 develop Mg/C/O/H ReaxFF parameter sets for two environments: an aqueous force field for magnesium ions in solution and an interfacial force field for minerals and mineral–water interfaces. Since magnesium is highly ionic, we choose to fix the magnesium charge and model its interaction with C/O/H through Coulomb, Lennard-Jones, and Buckingham potentials. We parameterize the forcefields against several crystal structures, including brucite, magnesite, magnesia, magnesium hydride, and magnesium carbide, as well as Mg 2+ water binding energies for the aqueous forcefield. Then, we test the forcefield for other magnesium-containing crystals, solvent separated and contact ion-pairs and single-molecule/multilayer water adsorption energies on mineral surfaces. We also apply the forcefield to the forsterite–water and brucite–water interface that contains a bicarbonate ion. We observe that a long-range proton transfer mechanism deprotonates the bicarbonate ion to carbonate at the interface. Free energy calculations show that carbonate can attach to the magnesium surface with an energy barrier of about 0.22 eV, consistent with the free energy required for aqueous Mg–CO 3 ion pairing. Also, the diffusion constant of the hydroxide ions in the water layers formed on the forsterite surface are shown to be anisotropic and heterogeneous. These findings can help explain the experimentally observedmore »fast nucleation and growth of magnesite at low temperature at the mineral–water–CO 2 interface in water-poor conditions.« 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 can 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 acidmore »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. Natural aquatic environments such as oceans, lakes, and rivers are home to a tremendous diversity of microorganisms. Some may cross the air-water interface within droplets and become airborne, with the potential to impact the Earth’s radiation budget, precipitation processes, and spread of disease. Larger droplets are likely to return to the water or adjacent land, but smaller droplets may be suspended in the atmosphere for transport over long distances. Here, we report on a series of controlled laboratory experiments to quantify wind-driven droplet production from a freshwater source for low wind speeds. The rate of droplet production increased quadratically with wind speed above a critical value (10-m equivalent 5.7 m/s) where droplet production initiated. Droplet diameter and ejection speeds were fit by a gamma distribution. The droplet mass flux and momentum flux increased with wind speed. Two mechanisms of droplet production, bubble bursting and fragmentation, yielded different distributions for diameter, speed, and angle. At a wind speed of about 3.5 m/s, aqueous suspensions of the ice-nucleating bacteriumPseudomonas syringaewere collected at rates of 283 cells m−2 s−1at 5 cm above the water surface, and at 14 cells m−2 s−1at 10 cm above the water surface. At a wind speed of about 4.0 m/s, aqueous suspensions ofP. syringaewere collectedmore »at rates of 509 cells m−2 s−1at 5 cm above the water surface, and at 81 cells m−2 s−1at 10 cm above the water surface. The potential for microbial flux into the atmosphere from aquatic environments was calculated using known concentrations of bacteria in natural freshwater systems. Up to 3.1 × 104 cells m−2 s−1of water surface were estimated to leave the water in potentially suspended droplets (diameters <100 µm). Understanding the sources and mechanisms for bacteria to aerosolize from freshwater aquatic sources may aid in designing management strategies for pathogenic bacteria, and could shed light on how bacteria are involved in mesoscale atmospheric processes.

    « less
  5. Abstract
    <p>This data set for the manuscript entitled &#34;Design of Peptides that Fold and Self-Assemble on Graphite&#34; includes all files needed to run and analyze the simulations described in the this manuscript in the molecular dynamics software NAMD, as well as the output of the simulations. The files are organized into directories corresponding to the figures of the main text and supporting information. They include molecular model structure files (NAMD psf or Amber prmtop format), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD configuration files, Colvars configuration files, NAMD log files, and NAMD output including restart files (in binary NAMD format) and trajectories in dcd format (downsampled to 10 ns per frame). Analysis is controlled by shell scripts (Bash-compatible) that call VMD Tcl scripts or python scripts. These scripts and their output are also included.</p> <p>Version: 2.0</p> <p>Changes versus version 1.0 are the addition of the free energy of folding, adsorption, and pairing calculations (Sim_Figure-7) and shifting of the figure numbers to accommodate this addition.</p> <p><br /> Conventions Used in These Files<br /> &#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;&#61;</p> <p>Structure Files<br /> ----------------<br /> - graph_*.psf or sol_*.psf (original NAMD (XPLOR?) format psf file including atom details (type, charge, mass),More>>