Search NASA⌕ Search

SEARCH · Search NASA

Results for “exploitations”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 559 records · Page 31

A Synchrotron-Based Vacuum Ultraviolet Photoionization Mass Spectrometer-Coupled Microreactor To Probe Thermocatalysis

Vacuum ultraviolet photoionization (VUV-PI) mass spectrometry offers an isomer-selective and universal ionization with minimal fragmentation detection of organics in complex chemical systems such as pyrolysis and combustion. Here, in this study, we report a state-of-the-art experimental setup of a universal catalytic microreactor combined with a molecular beam to investigate the thermocatalytic oxidation of a heterogeneous system relevant for probing reactions at gas–solid interfaces. In strong contrast to traditional off-line analytical methods, this technique is capable of identifying and quantifying short-lived species (radicals) as well as stable products to decipher initial reaction steps via the detection of nascent products. The thermocatalytic oxidative degradation of exo-tetrahydrodicyclopentadiene (JP-10), a high energy-density hydrocarbon fuel, over solid titanium–aluminum–boron reactive mixed metal nanopowder (Ti-Al-B RMNP) is exploited to showcase potential applications. Overall, some 59 nascent gas-phase products are identified via photoionization efficiency (PIE) curves, including oxygenated species and hydrocarbons comprising closed-shell molecules and radicals. The critical temperature for complete oxidative decomposition of JP-10 was lowered by 450 K from 1400 K to 950 K, indicating an efficient thermocatalytic action of Ti-Al-B nanoparticles on JP-10. The enabling of a universal chemical microreactor along with VUV-PI mass spectrometry broadens the applicability of this technique to hydrocarbon fuel oxidation and pyrolysis characterization. This isomer-selective sensitive probing along with the detection of radical transients makes the aforementioned technique superior to other conventional analytical techniques such as microflow tube and pyrolysis-gas chromatography coupled with mass spectrometry for investigating similar pyrolysis reactions and comprehensive quantification.

Paul, Dababrata [Univ. of Hawaii at Manoa, Honolul↗

Mechanisms of Three-Dimensional Solid-Phase Epitaxial Crystallization of Strontium Titanate

Strontium titanate (SrTiO 3 , STO) is a complex metal oxide with a cubic perovskite crystal structure. Due to its easily described and understood crystal structure in the cubic phase, STO is an ideal model system for exploring the mechanistic details of solid-phase epitaxy (SPE) in complex oxides. SPE is a crystallization approach that aims to guide crystal growth at low homologous temperatures to achieve targeted microstructures. Beyond planar thin films, SPE can also exploit the addition of a chemically inert, noncrystallizing, amorphous obstacle in the path of crystallization to generate complex three-dimensional structures. The introduction of this mask fundamentally alters the SPE process, inducing a transition from two- to three-dimensional geometries and from vertical to lateral crystal growth under the influence of the crystal/mask/amorphous boundary. Using a combination of molecular dynamics simulations and experiments, we identify several unique phenomena in the nanoscale growth behaviors in both conventional (unmasked) and masked SPE. Examining conventional SPE of STO, we find that crystallization at the interface is strongly correlated to, and potentially driven by, density fluctuations in the region of the amorphous STO near the crystalline/amorphous interface with a strong facet dependence. In the masked case, we find that the crystalline growth front becomes nonplanar near contact with the mask. We also observe a minimum vertical growth requirement prior to lateral crystallization. Both phenomena depend on the relative bulk and interfacial free energies of the three-phase (crystal/mask/amorphous) system.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

On the Structural Origin of Fast Li-Ion Cycling in Tetragonal Bronze-Type Nb 8 W 9 O 47

Bronze and bronze-derived, mixed metal oxides are materials of growing interest for lithium-ion battery anodes due to their high-rate capabilities and the potential for high energy densities via multielectron redox. We report on the synthesis and electrochemical properties of the tetragonal tungsten bronzetype phase Nb 8 W 9 O 47 and our investigation into the structural evolution of this phase upon lithium (de)intercalation, with particular attention to how it relates to the fast-cycling capabilities. Electrochemical cycling shows Nb 8 W 9 O 47 can achieve greater than one Li + per transition metal at rates of C/2 or slower and maintains a capacity equivalent to 0.65 Li + per transition metal at a rate of 20C. Sequential Rietveld analysis of operando X-ray diffraction data reveals anisotropic changes to the unit cell parameters during discharge, accompanied by an off-centering of the transition metals within their octahedral environments. While higher capacities can be accessed below 1.2 V, cells discharged to 1.0 V exhibit a significant expansion of the unit cell volume and reduced Coulombic efficiencies. This report suggests that rigid crystal structures incapable of undergoing polyhedral rotational deformations can instead exploit displacive distortions of the transition metal redox center to enable fast Li-ion cycling with minimal atomic rearrangements.

25 ENERGY STORAGE↗

Single Component Dye-Sensitized Solar Cells Enabled by Copper Chemistry: Introduction of the Retro Cell

The possibility of utilizing a single molecule to act as both a chromophore and a redox shuttle in a new configuration of a dye-sensitized solar cell (DSSC) is investigated. This design, termed a retro cell, exploits the copper chromophore, [Cu(dsbtmp) 2 ] + (dsbtmp = bis(2,9-di(sec-butyl)-3,4,7,8- tetramethyl-1,10-phenanthroline), which has been shown to have excited state lifetimes in excess of a microsecond, enabling sensitization of TiO 2 while dissolved in solution. The oxidized chromophore can then diffuse to the counter electrode to be regenerated. This concept simplifies the device components and fabrication and eliminates a charge transfer step compared to that of traditional DSSCs. Initial investigations show the concept is viable; however, the performance is suboptimal. We have found the addition of 4-tertbutylpyridine (TBP) to the electrolyte plays a crucial role in enabling solar energy conversion. Evidence of TBP displacing one of the dsbtmp ligands in the oxidized [Cu(dsbtmp) 2 ] 2+ complex has been presented, which likely plays an important role in reducing recombination and enabling charge collection. The performance-limiting steps and routes to improved performance and viability of a retro cell are further discussed.

14 SOLAR ENERGY↗

Mechanisms of Polyethylene Terephthalate Pellet Fragmentation into Nanoplastics and Assimilable Carbons by Wastewater Comamonas

Comamonadaceae bacteria are enriched on poly(ethylene terephthalate) (PET) microplastics in wastewaters and urban rivers, but the PET-degrading mechanisms remain unclear. Here, we investigated these mechanisms with Comamonas testosteroniKF-1, a wastewater isolate, by combining microscopy, spectroscopy, proteomics, protein modeling, and genetic engineering. Compared to minor dents on PET films, scanning electron microscopy revealed significant fragmentation of PET pellets, resulting in a 3.5-fold increase in the abundance of small nanoparticles (<100 nm) during 30-day cultivation. Infrared spectroscopy captured primarily hydrolytic cleavage in the fragmented pellet particles. Solution analysis further demonstrated double hydrolysis of a PET oligomer, bis(2-hydroxyethyl) terephthalate, to the bioavailable monomer terephthalate. Supplementation with acetate, a common wastewater co-substrate, promoted cell growth and PET fragmentation. Of the multiple hydrolases encoded in the genome, intracellular proteomics detected only one, which was found in both acetate-only and PET-only conditions. Homology modeling of this hydrolase structure illustrated substrate binding analogous to reported PET hydrolases, despite dissimilar sequences. Mutants lacking this hydrolase gene were incapable of PET oligomer hydrolysis and had a 21% decrease in PET fragmentation; re-insertion of the gene restored both functions. Thus, we have identified constitutive production of a key PET-degrading hydrolase in wastewater Comamonas, which could be exploited for plastic bioconversion.

54 ENVIRONMENTAL SCIENCES↗

Adsorption and Chromatographic Behavior of Dispersed Sodium Bismuthate Systems for the Separation of Americium from Curium

The selective partitioning of americium (Am) and curium (Cm) is integral for nuclear science areas such as the nuclear fuel cycle, stockpile stewardship, and isotope production but remains a long-standing radioanalytical challenge due to nearly identical chemical properties. Differences in redox chemistry can be exploited since Am can be oxidized to the hexavalent oxidation state. Sodium bismuthate (NaBiO 3 ) exhibits favorable oxidation and ion exchange characteristics conducive to a rapid and efficient Am/Cm chromatographic separation. Contact of nitric acid with NaBiO 3 was shown to significantly decrease the nitric acid concentration and solution volume. The adsorption, kinetic, and chromatographic behavior of Am and Cm in systems that disperse NaBiO3 in filter aids was evaluated. Further, dispersion of NaBiO 3 increased separation factors to over 100, with rapid kinetics achieved within 1 min of contact and sustained for at least 2 h. The adsorption capacity was more than doubled from 0.066 mequiv g –1 for solid NaBiO 3 alone to 0.149(6) and 0.156(9) meq g –1 when dispersed in Celite 545 and silica gel, respectively. Complete separation was achieved in all systems with ~100% recovery in the respective eluted fractions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Benders Decomposition Using Graph Modeling and Multi-Parametric Programming

Benders decomposition is a widely used method for solving large and structured optimization problems, but its performance is affected by the repeated solution of subproblems. We propose a flexible and modular algorithmic framework for accelerating Benders decomposition. Specifically, we express the problem structure by using a graph-theoretic modeling abstraction in which nodes represent optimization subproblems and edges represent connectivity between subproblems. A key innovation of our approach is that we embed multiparametric programming (mp) surrogates for node subproblems, which maps the exact analytical map of the subproblem solution space. The use of mp surrogates allows us to replace subproblem solves with fast look-ups and function evaluations for primal and dual variables during the iterative Benders process. We formally show the equivalence between classical Benders cuts and those derived from the mp solution. We implement our framework in the open-source PlasmoBenders.jl software package. To demonstrate the capabilities of the proposed framework, we apply it to a two-stage stochastic programming problem, which aims to make optimal capacity expansion decisions under market uncertainty. We evaluate both single-cut and multicut variants of Benders decomposition and show that the use of mp surrogates achieves substantial speedups in subproblem solve time, while preserving the convergence guarantees of Benders decomposition. We highlight advantages in solution analysis and interpretability that is enabled by mp critical region tracking; specifically, we show that these reveal how decisions evolve geometrically across the Benders search. Our results aim to demonstrate that combining surrogate modeling with graph modeling offers a promising and extensible foundation for structure-exploiting decomposition. In addition, by decomposing the problem into more tractable subproblems, the proposed approach also aims to overcome scalability issues of mp. Finally, the use of mp surrogates provides a unifying and modular optimization framework that enables the representation of heterogeneous node subproblems as modeling objects with a homogeneous structure.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The Curious Case of Pu + : Insight on 5f Orbital Activity from Inductively Coupled Plasma Tandem Mass Spectrometry (ICP-MS/MS) Reactions

A comprehensive understanding of when and how 5f orbitals participate in complex chemical bonding is important for a variety of applications. The actinides are unique in that they possess 5f orbitals and can access high oxidation states, which make them attractive for use in catalysis. Fundamental studies of actinide–ligand interactions offer a mechanism to examine the activation of the 5f orbitals so that the selectivity of 5f orbitals can be assessed. A previous study examined the reaction of Pu + + CO 2 and determined that the reaction efficiency is restricted by a barrier, namely, promotion from the Pu + ground-state configuration, 5f 6 7s, to a reactive-state configuration, 5f 5 6d 2 . Here, the present study illustrates the benefit of activation of Pu’s 5f orbitals when studying the reaction of Pu + + NO. In this reaction, PuO + forms in an exothermic, barrierless process. The 5f orbitals can and do participate in forming a linear intermediate, [N-Pu-O] + , and this drives the exothermic reaction. Understanding the conditions under which 5f orbitals are active in chemical bonding is the key to exploiting the actinides’ selective catalytic capabilities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Impact of Pendant Amine Basicity on Electrochemically-Promoted Cobalt Hydride Formation: Kinetic and Mechanistic Analysis

Here, we report the role of pendant amine basicity on the proton-coupled electron transfer (PCET) reactivity for the conversion of [Co III Cp(P Ph 2 N R 2 )(CH 3 CN)] 2+ complexes to [HCo III Cp(P Ph 2 N R 2 )] + , which is a key transformation involved in catalytic CO 2 conversion to formate and in H 2 evolution. Three complexes were studied, where the amine substituent (R) varies from benzyl, methoxyphenyl, or phenyl. In previous work on the benzyl system, we showed that the amine on the P Ph 2 N Bn 2 ligand serves as a kinetically accessible protonation site and enables three participating hydride formation mechanisms. In this work, a combination of electrochemical measurements and theoretical calculations were used to show that the electronic donation at the pendant amine influences the accessible PCET mechanism and proton transfer kinetics related to cobalt hydride formation under analogous reaction conditions. Notably, the amine with the most electron-donating substituent correlates to the lowest barrier for amine protonation, and specific cobalt hydride formation mechanisms can be shut off for the amine with the least electron-donating substituent. The mechanistic and kinetic changes upon modulation of the amine substituent have great implications for overall catalytic efficiency and selectivity, especially to generate the cobalt hydride intermediate involved in selective CO 2 reduction to formate. This work shows how to exploit kinetic basicity using ligand-cooperative design to facilitate PCET reactions involved in energy related transformations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Multi-Spectroscopic Determination of Exchange Coupling, Zero-Field Splitting, and g-Matrices in Radical-Bridged Dinuclear Fe(III) Complexes

When the energy gap, Δ, between the lowest-lying spin manifolds within a spin-exchange coupled molecule approaches Δ/k B ≈ 300 K, the traditional temperature-dependence (T < 400 K) of the molar magnetic susceptibility is not always a reliable way to obtain a good estimate of intramolecular exchange couplings. We develop a spectroscopic approach capable of accurately parametrizing complex magnetic Hamiltonians by exploiting the separation of the anisotropy and exchange energy scales in strongly coupled magnetic molecules. Specifically, we combine inelastic neutron scattering, high-frequency electron paramagnetic resonance, far-infrared magneto-spectroscopy and magnetometry, and obtain detailed information about the magnetic properties of a series of diiron complexes derived from [[Fe(cth)] 2 (dxbq)] 3+ (H 2 dxbq: 2,5-dihydroxy-1,4-benzoquinone (x = h) or 3,6-dichloro-2,5-dihydroxy-1,4-benzoquinone (x = c), cth: 5,5,7,12,12,14-hexamethyl-1,4,8,11-tetraazacyclotetradecane). Well-isolated S = 9/2 ground states emerge due to strong direct antiferromagnetic exchange between the Fe 3+ centers (S = 5/2) and the radical bridging benzoquinone ligand (S = 1/2). The specific sensitivities and transition selection rules of the applied methods allow us to determine the parameters of the microscopic Hamiltonian including exchange coupling, fourth-order Stevens operators and g-factors. Our methodology is directly portable to other strongly coupled molecular compounds.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Efficient Mixed-Precision Matrix Factorization of the Inverse Overlap Matrix in Electronic Structure Calculations with AI-Hardware and GPUs

In recent years, a new kind of accelerated hardware has gained popularity in the artificial intelligence (AI) community which enables extremely high-performance tensor contractions in reduced precision for deep neural network calculations. In this article, we exploit Nvidia Tensor cores, a prototypical example of such AI-hardware, to develop a mixed precision approach for computing a dense matrix factorization of the inverse overlap matrix in electronic structure theory, S –1 . This factorization of S –1 , written as ZZT = S –1 , is used to transform the general matrix eigenvalue problem into a standard matrix eigenvalue problem. Here we present a mixed precision iterative refinement algorithm where Z is given recursively using matrix–matrix multiplications and can be computed with high performance on Tensor cores. To understand the performance and accuracy of Tensor cores, comparisons are made to GPU-only implementations in single and double precision. Additionally, we propose a nonparametric stopping criteria which is robust in the face of lower precision floating point operations. The algorithm is particularly useful when we have a good initial guess to Z, for example, from previous time steps in quantum-mechanical molecular dynamics simulations or from a previous iteration in a geometry optimization.

36 MATERIALS SCIENCE↗

An Orbital Basis Set for Double Photoionization of Atoms and Molecules

The ab initio theoretical treatment of one-photon double photoionization processes has been limited to atoms and diatomic molecules by the challenges posed by large grid-based representations of the double ionized continuum wave function. To provide a path for extensions to polyatomics, an energy-adapted orbital basis approach is demonstrated that reduces the dimensions of such representations and simultaneously allows larger time steps in time-dependent computational descriptions of double ionization. Additionally, an algorithm that exploits the diagonal nature of the two-electron integrals in the grid basis and dramatically accelerates the transformation between grid and orbital representations is presented. Excellent agreement between the present results and benchmark theoretical calculations is found for H – and Be atoms, as well as the hydrogen molecule, including for the triply differential cross sections that relate the angular distribution and energy sharing of all of the particles in the molecular frame.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

GPU-Accelerated Solution of the Bethe–Salpeter Equation for Large and Heterogeneous Systems

We present a massively parallel GPU-accelerated implementation of the Bethe–Salpeter equation (BSE) for the calculation of the vertical excitation energies (VEEs) and optical absorption spectra of condensed and molecular systems, starting from single-particle eigenvalues and eigenvectors obtained with density functional theory. The algorithms adopted here circumvent the slowly converging sums over empty and occupied states and the inversion of large dielectric matrices through a density matrix perturbation theory approach and a low-rank decomposition of the screened Coulomb interaction, respectively. Further computational savings are achieved by exploiting the nearsightedness of the density matrix of semiconductors and insulators to reduce the number of screened Coulomb integrals. We scale our calculations to thousands of GPUs with a hierarchical loop and data distribution strategy. The efficacy of our method is demonstrated by computing the VEEs of several spin defects in wide-band-gap materials, showing that supercells with up to 1000 atoms are necessary to obtain converged results. We discuss the validity of the common approximation that solves the BSE with truncated sums over empty and occupied states. In conclusion, we then apply our GW-BSE implementation to a diamond lattice with 1727 atoms to study the symmetry breaking of triplet states caused by the interaction of a point defect with an extended line defect.

Absorption spectra↗

A Combined Crossed Molecular Beam and Theoretical Investigation of the Elementary Reaction of Tricarbon (C 3 (X 1 Σ g + )) with Diacetylene (C 4 H 2 (X 1 Σ g + )): Gas Phase Formation of the Heptatriynylidyne Radical ( l -C 7 H(X 2 Π))

An elucidation of the underlying formation pathways to acyclic hydrocarbons such as polyynes (C n H 2 ), cumulenes (C n H 2 ), and linear resonantly stabilized linear radicals (l-C n H) is indispensable to understand the hydrocarbon chemistry in extreme low and high temperature environments. In this study, we exploited the crossed molecular beam technique to investigate the reaction of tricarbon C 3 (X 1 Σ g + ) with diacetylene (butadiyne; HCCCCH; X 1 Σ g + ) at a collision energy of 47 ± 1 kJ mol ⁻1 . The experimental data were merged with ab initio calculations of the singlet C 7 H 2 potential energy surface (PES) revealing that the reaction is initiated via the formation of an initial van der Waals reactant complex in the entrance channel. Subsequent rearrangements lead to various carbene-type and cyclic intermediates via ring-opening, ring-closure, and hydrogen migration processes eventually forming acyclic C 7 H 2 isomers prior to their barrierless unimolecular decomposition to the most stable linear isomer, heptatriynylidyne (C 7 H, X 2 Π) in an overall endoergic reaction (+57 kJ mol ⁻1 ). The reaction exhibits strong similarities to the tricarbon – acetylene (C 3 – C 2 H 2 ). Furthermore, the significant energy threshold suggests that the tricarbon reaction with (poly)acetylenes forming resonantly stabilized linear radicals are open in high-temperature environments such as combustion flames and circumstellar envelopes of carbon stars and planetary nebulae as their descendants; however, these reactions are closed in low-temperature environments as in cold molecular clouds and hydrocarbon-rich atmospheres of planets and their moons such as in Titan.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Collisional Excitation of HCN by CO to Refine the Modeling of Cometary Comae

Here, we present the first dataset of collisional (de)-excitation rate coefficients of HCN induced by CO, one of the main perturbing gases in cometary atmospheres. The dataset spans the temperature range of 5–50 K. It includes both state-to-state rate coefficients involving the lowest ten and nine rotational levels of HCN and CO, respectively, and the so-called “thermalized” rate coefficients over the rotational population of CO at each kinetic temperature. The derivation of these coefficients exploited the good performance of the statistical adiabatic channel model (SACM) on top of an accurate interaction potential computed at the CCSD(T)-F12b/CBS level of theory. The reliability of the SACM approach was validated by comparison with full quantum calculations restricted at the lowest total angular momentum of the system. These results provide essential input to accurately model the distribution among the rotational energy levels and the abundance of HCN in cometary atmospheres, accounting for deviations from local thermodynamic equilibrium that typically occurs in such environments.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Reconciling the Different Apparent Decay Pathways for Iron(III) Tetraphenylporphyrin Chloride: Evidence for Excited State Branching

Detailed understanding of the photophysical properties of metalloporphyrins is key to rationally exploiting them in a variety of applications ranging from photocatalysis to opto-magnetics. Previous studies of the ferric-tetraphenylporphyrin chloride have provided contradictory descriptions of the excited state evolution. Optical transient-absorption suggested initial formation of a ππ* excited state, followed by ligand-to-metal charge-transfer (LMCT) from the porphyrin ring and then decay on the ∼2 ps time scale to a metal-centered excited state that had a lifetime of ∼15 ps. In contrast, femtosecond extreme ultraviolet transient-absorption at the Fe M 2,3 -edge, while agreeing on the initial formation of an LMCT state, found evidence that this decayed in ∼2 ps to the ground state. Here, we have used K-edge transient X-ray absorption and X-ray emission, together with time-resolved X-ray solution scattering to explore this system. Based on these data, we propose a new model, consistent with both the earlier and the current data, in which photoexcited FeTPPCl evolves through three different states on the LMCT manifold, with the ∼2 ps decay now seen to involve a branching between return to the ground state (∼70%) and formation of a long-lived LMCT state (∼30%).

Excited states↗

Cryogenic Vibrational Spectroscopy of the Deprotonated Dimer of Phosphoric Acid

Phosphate-containing molecules are ubiquitous in nature, where they play crucial roles in biochemical processes. Further, they are of technical importance, for example, in certain batteries and in fuel cells, where a unique property of phosphoric acid is exploited its exceptionally high proton conductivity. Proton transport in phosphoric acid is known to involve proton shuttling; however, the elementary steps involved are not clear. To elucidate the hydrogen bonding preferences of phosphoric acid, we investigate the dihydrogen phosphate anion as well as the deprotonated dimer of phosphoric acid (H 3 PO 4 ·H 2 PO 4 − ) in the gas phase using infrared action spectroscopy in helium nanodroplets and infrared D 2 -tagging photodissociation spectroscopy, and the experimental spectra are compared to theoretical ones. Theory finds for H 3 PO 4 ·H 2 PO 4 − two different structures that are predicted to be nearly isoenergetic. The comparison to the experimental spectra, however, allows for a clear assignment and structure identification. The resulting structure has an interesting binding motif, which might be of relevance to interactions of phosphoric acid in the condensed phase and which can serve as a benchmark for quantum chemical calculations.

Chemical structure↗

Accurate Prediction of pKb in Amines: Validation of the CAM-B3LYP/6-311+G(d,p)/SMD Model

Amines play several key roles in chemistry and biology and are involved in numerous industrial processes, often with significant economic impacts. Recently, amines are also garnering interest as catalysts for polymer synthesis and for CO 2 fixation, incentivizing the need to rapidly design and screen new amino compounds. Hence, developing reliable methods to predict their physicochemical properties, e.g., the base dissociation constant (pKb), is pivotal. Here, a density functional theory (DFT)-based approach was employed to compute the pKb of substituted amines, exploring the impact of several key parameters, including (i) the number of explicit water molecules at the reaction center, (ii) the van der Waals (vdW) surface, and (iii) solvent polarizability. In previous work, it was determined that including two explicit water molecules at the reaction center resulted in highly accurate pKb estimates for primary amines. Here, we find that including a third water molecule at the reaction center is essential for accurate pKb for secondary and tertiary amines. The revised methodology was then applied to a wider selection of amines, obtaining a minimum average error (MAE) < 0.4. In conclusion, this result represents an extension of our “easy-to-use method,” a simple and direct DFT approach exploiting CAM-B3LYP/SMD/6-311G+(d,p) to compute pKb without post facto modifications.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗