Search NASA⌕ Search

SEARCH · Search NASA

Results for “electrostatic interaction”

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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

An Analysis of the Novel Fluorocycline TP-6076 Bound to Both the Ribosome and Multidrug Efflux Pump AdeJ from Acinetobacter baumannii

Antibiotic resistance among bacterial pathogens continues to pose a serious global health threat. Multidrug-resistant (MDR) strains of the Gram-negative organism Acinetobacter baumannii utilize a number of resistance determinants to evade current antibiotics. One of the major resistance mechanisms employed by these pathogens is the use of multidrug efflux pumps. These pumps extrude xenobiotics directly out of bacterial cells, resulting in treatment failures when common antibiotics are administered. Here, the structure of the novel tetracycline antibiotic TP-6076, bound to both the Acinetobacter drug efflux pump AdeJ and the ribosome from Acinetobacter baumannii, using single-particle cryo-electron microscopy (cryo-EM), is elucidated. In this work, the structure of the AdeJ–TP-6076 complex is solved, and we show that AdeJ utilizes a network of hydrophobic interactions to recognize this fluorocycline. Concomitant with this, we elucidate three structures of TP-6076 bound to the A. baumannii ribosome and determine that its binding is stabilized largely by electrostatic interactions. We then compare the differences in binding modes between TP-6076 and the related tetracycline antibiotic eravacycline in both targets. These differences suggest that modifications to the tetracycline core may be able to alter AdeJ binding while maintaining interactions with the ribosome. Together, this work highlights how different mechanisms are used to stabilize the binding of tetracycline-based compounds to unique bacterial targets and provides guidance for the future clinical development of tetracycline antibiotics.

36 MATERIALS SCIENCE↗
NASA NTRS1981

The plasma wake of mesosonic conducting bodies. I - An experimental parametric study of ion focusing by the plasma sheath

The experimental investigation considered is concerned with the deflection of ion streams resulting from the interaction of conducting test bodies with an unmagnetized, mesosonic (supersonic with respect to ions but subsonic with respect to electrons) plasma stream. The investigation is, therefore, limited to plasma-electrostatic interactions. The experimental conditions are similar to those of the spacecraft-ionospheric interaction in that the ionic mass and number density, the electron temperature, and the plasma drift velocity ranges include values appropriate for small satellites or diagnostic probes at 200 to 400 km altitude. The study provides direct observations of deflected ion streams for cylindrical test bodies and gives a detailed description of the effects of the governing, dimensionless parameter ratios obtained from the steady-state, nonmagnetic Maxwell-Vlasov system of equations.

Stone, N. H.↗
DOE OSTI2020

Experimental and DFT Calculated IR Spectra of Guests in Zeolites: Acyclic Olefins and Host-Guest Interactions

We performed experimental and periodic density functional theory (DFT) IR spectroscopy to investigate the adsorption of acyclic olefins over both acidic and non-acidic zeolites. Two conjugated polyenes, 2,4-dimethyl-1,3-pentadiene (I) and 2,6-dimethyl-2,4,6-octatriene (II) were studied to probe organic intermediates that can form during methanol conversion, and that can lead to deactivating species known collectively as “coke.” We computed vibrational spectra using zeolite-adsorbed and gas-phase models for both neutral and protonated forms of I and II, and compared these DFT results to diffuse reflectance IR Fourier transform (DRIFT) spectra of zeolite-guest systems. Our experimental and computational results are precise enough to pinpoint that the gauche s-cis conformation of species I predominates during adsorption over de-aluminated zeolite beta. Computed zeolite-adsorbed spectra of the protonated species I and II best represent the DRIFT spectra obtained after the adsorption of the olefins on HMOR at 20 ?, with computed bands at 1543 and 1562 cm–1 for molecules I+ and II+, respectively, attributed to the allylic stretching mode, ?(C=C–C+). These computed band frequencies are within 6 cm–1 of experimental data and confirm that the interaction between neutral acyclic olefins and acidic zeolites leads to protonation of the olefin. Comparison of computed spectra of the protonated species in the gas phase to those in the zeolite indicates that the electrostatic interaction between alkenyl and alkadienyl cations and negative zeolite framework does not significantly impact the position of the allylic stretching bands. These results highlight that computed spectroscopy and thermodynamics coupled with experimental spectra can be used for elucidating complex mixtures in zeolites, and certain spectral features of adsorbed olefins can be accurately modeled by gas-phase calculations. ACKNOWLEDGMENTS BM and SMA were supported by the National Science Foundation under Award # CBET-1512442. Acknowledgment is made to the Donors of the American Chemical Society Petroleum Research Fund for support (or partial support) of this research (#589978-ND5, EDH and FCJ). MDB was supported by the US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences (BES), Division of Materials Science and Engineering. CJM was supported by the DOE, Office of Science, BES, Division of Chemical Sciences, Geosciences, and Biosciences. Pacific Northwest National Laboratory (PNNL) is a multiprogram national laboratory operated by Battelle for the U.S. DOE. Computing resources were generously provided through PNNL’s institutional computing (PIC).

Manookian, Babgen↗
DOE OSTI2018

Interaction of dequalinium chloride with phosphatidylcholine bilayers: A biophysical study with consequences on the development of lipid-based mitochondrial nanomedicines

Dequalinium (DQ) has been proposed as a mitochondrial targeting ligand for nanomedicines, including liposomes, given the implication of these organelles in many diseases. This original study focuses on the interactions of DQ with phosphatidylcholine bilayers during the formation of liposomes. Firstly, PEGylated liposomes suitable for drug delivery were studied and were found to be more stable when made in water than in phosphate-buffered saline, emphasizing the role of electrostatic interactions between positive charges on DQ and the polar head groups of the lipids. To gain more information, differential scanning calorimetry, small- and wide-angle X-ray scattering and diffraction, 31 P and 2 H NMR spectroscopy and freeze-fracture electron microscopy were performed on dimyristoylphosphatidylcholine (DMPC) model membranes in the presence of DQ. This molecule was shown to be located at the level of polar head groups and to induce electrostatic repulsions between adjacent lipid bilayers leading to membrane budding in water. These findings indicate that DQ is not completely inert towards lipid membranes and therefore is not an ideal candidate for encapsulation in liposomes. As a whole, our work stresses the necessity for thorough physico-chemical characterization to better understand the mechanisms underlying the development of nanomedicines.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗
DOE OSTI2021

Polymorphism and Optoelectronic Properties in Crystalline Supramolecular Polymers

Supramolecular polymers can emulate some of the physical properties of covalent polymers but offer new opportunities given the possibility of designing monomers that will form highly ordered assemblies with defined shapes. Internally ordered supramolecular polymers formed through nucleation–elongation self-assembly are well-known but highly crystalline examples which exhibit important properties such as light harvesting, charge transport, and ferroelectricity are not common. We report here on a detailed study of supramolecular polymers formed in water by carboxylated naphtho-p-quinodimethane amphiphiles. We found that supramolecular polymerization of these amphiphiles in aqueous media yields crystalline assemblies with morphologies that included ribbons, helically rolled ribbons, and twisted filaments. This polymorphism was found to be controlled exclusively by repulsive electrostatic interactions controlled by the degree of protonation of the carboxylic head groups which also dictates the nature of supramolecular packing. Substoichiometric amounts of base lead to highly crystalline ribbons due to a decreased surface charge density and less electrostatic repulsion. Increasing deprotonation results in helically rolled ribbons with a different polymorph crystal lattice, whereas excessive deprotonation leads to twisted filaments with maximum surface charge density. Ribbons, helical rolled ribbons, and twisted filaments revealed an increasing red shift in their visible absorption maxima. These crystalline assemblies could be potential candidates for solar energy materials and photocatalytic systems.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Harnessing Intermolecular Interactions to Promote Long-Lived Photoinduced Charge Separation from Copper Phenanthroline Chromophores

In this study, facilitating photoinduced electron transfer (PET) while minimizing rapid charge-recombination processes to produce a long-lived charge-separated (CS) state represents a primary challenge associated with achieving efficient solar fuel production. Natural photosynthetic systems employ intermolecular interactions to arrange the electron-transfer relay in reaction centers and promote a directional flow of electrons. This work explores a similar tactic through the synthesis and ground- and excited-state characterization of two Cu(I)bis(phenanthroline) chromophores with homoleptic and heteroleptic coordination geometries and which are functionalized with negatively charged sulfonate groups. The addition of sulfonate groups enables solubility in pure water, and it also induces assembly with the dicationic electron acceptor methyl viologen (MV 2+ ) via bimolecular, dynamic electrostatic interactions. The effect of the sulfonate groups on the ground- and excited-state properties was evaluated by comparison with the unsulfonated analogues in 1:1 acetonitrile/water. The excited-state lifetimes for all sulfonated complexes are similar to what we expect from previous literature, with the exception of the sulfonated heteroleptic complex whose metal-to-ligand charge-transfer (MLCT) lifetime in water has two components that are fit to 10 and 77 ns. For the sulfonated complexes, we detected reduced MV +$\bullet$ in both solvent environments following MLCT excitation, but control measurements in 1:1 acetonitrile/water with the unsulfonated analogues showed no PET to MV 2+ , indicating that electrostatically driven supramolecular assemblies of the sulfonated complexes with MV 2+ facilitate the observed PET. Additionally, the strength of the intermolecular interactions driving the formation of these assemblies changes drastically with the solvent environment. In 1:1 acetonitrile/water, PET occurred from both sulfonated complexes with quantum yields (Φ ET ) of 2-3% but increased to a remarkable 98% for the sulfonated heteroleptic complex with a 3 μs CS-state lifetime in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Interacting dust grains in complex plasmas: Ion wake formation and the electric potential

Dust grains have been used as minimally invasive probes to determine plasma parameters including the plasma density, temperature, and electric field in a plasma discharge. However, the dust grains in a plasma generate local potential disturbances due to the collection of charge and the subsequent electrostatic interactions between the dust and charged plasma particles. Dust grains in close proximity to one another exhibit interesting non-reciprocal interactions and self-organize into structures such as one-dimensional filamentary chains, two-dimensional “zigzags,” and three-dimensional helices, among others. The formation of these structures suggests that although the dust grains may be less invasive than traditional plasma probes, the disturbance to the local plasma environment introduced by dust grains is non-trivial. Commonly used analytic forms of the electric potential describing complex plasmas have failed to resolve the near-dust region, and as a result are insufficient to provide insight about the formation of complex dust structures. Here, we use an N-body simulation to compute the electric potential from ion densities near various dust grain configurations. We provide an alternative description to the standard analytic model for the electric potential of dust and ion wakes based on a Gaussian shaped cloud of ions. The electric potential obtained from simulations is used to identify minimum energy configurations for two and three dust grains. It is further demonstrated that the minimum potential region identified for N dust grains and their associated ion wakes does not predict the minimum-energy configuration of N + 1 dust grains.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Inhomogeneous Distribution of Cationic Surfactants around Anionic Molecular Clusters

An interesting phenomenon is reported when uranyl peroxide nanoclusters U 60 (Li 48+ m K 12 (OH) m [UO 2 (O 2 )(OH)] 60 (H 2 O) n , m ≈20 and n ≈310) interact with a small number of cationic surfactant molecules. Cationic surfactant molecules do not distribute evenly around the U 60 clusters during the interaction as expected. Instead, a small fraction of U 60 clusters attract almost all the surfactant molecules, leading to the self-assembly into supramolecular structures by using surfactant–U 60 complexes as building locks, and later further aggregate and precipitate based on hydrophobic interaction, whereas the rest of the clusters remained unbounded soluble macroions in bulk dispersion. This phenomenon nicely demonstrates a unique feature of macroion solutions. Furthermore, considering that Debye–Hückel approximation is no longer valid in such solutions, the competition between the local electrostatic interaction and hydrophobic interaction becomes important to regulate the solution behaviors of macroions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Nonlinear Poisson–Boltzmann solutions for charged parallel plates: When opposite charges repel

I present an exact solution of the Poisson–Boltzmann equation for two parallel plates and discuss the solution properties. I discuss in more detail plates with opposite charges: In this case, there are two critical separations, L c,1 < L c,2 . For separations less than L c,1 , the force between plates is repulsive. It switches to attractive at L c,1 , but with the electric potential having the same sign on both plates. For L > L c,2 , the force remains attractive, and the potential at the plates has the same sign as the charge on each plate. I also describe charge regulation, determined by pK a , and provide formulas for both the critical distance where oppositely charged plates repel and their charging process. Finally, the implications of these results for the nanoparticle assembly, as driven by electrostatic interactions, are also discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2022

Computational Approach for Probing Redox Potential for Iron-Sulfur Clusters in Photosystem I

Photosystem I is a light-driven electron transfer device. Available X-ray crystal structure from Thermosynechococcus elongatus showed that electron transfer pathways consist of two nearly symmetric branches of cofactors converging at the first iron–sulfur cluster F X , which is followed by two terminal iron–sulfur clusters F A and F B . Experiments have shown that F X has lower oxidation potential than F A and F B , which facilitates the electron transfer reaction. Here, we use density functional theory and Multi-Conformer Continuum Electrostatics to explain the differences in the midpoint Em potentials of the F X , F A and F B clusters. Our calculations show that F X has the lowest oxidation potential compared to F A and F B due to strong pairwise electrostatic interactions with surrounding residues. These interactions are shown to be dominated by the bridging sulfurs and cysteine ligands, which may be attributed to the shorter average bond distances between the oxidized Fe ion and ligating sulfurs for F X compared to F A and F B . Moreover, the electrostatic repulsion between the 4Fe-4S clusters and the positive potential of the backbone atoms is lowest for F X compared to both F A and F B . These results agree with the experimental measurements from the redox titrations of low-temperature EPR signals and of room temperature recombination kinetics.

59 BASIC BIOLOGICAL SCIENCES↗
NASA NTRS2016

Evaluation of Methods for Molecular Dynamics Simulation of Ionic Liquid Electric Double Layers

We investigate how systematically increasing the accuracy of various molecular dynamics modeling techniques influences the structure and capacitance of ionic liquid electric double layers (EDLs). The techniques probed concern long-range electrostatic interactions, electrode charging (constant charge versus constant potential conditions), and electrolyte polarizability. Our simulations are performed on a quasi-two-dimensional, or slab-like, model capacitor, which is composed of a polarizable ionic liquid electrolyte, [EMIM][BF4], interfaced between two graphite electrodes. To ensure an accurate representation of EDL differential capacitance, we derive new fluctuation formulas that resolve the differential capacitance as a function of electrode charge or electrode potential. The magnitude of differential capacitance shows sensitivity to different long-range electrostatic summation techniques, while the shape of differential capacitance is affected by charging technique and the polarizability of the electrolyte. For long-range summation techniques, errors in magnitude can be mitigated by employing two-dimensional or corrected three dimensional electrostatic summations, which lead to electric fields that conform to those of a classical electrostatic parallel plate capacitor. With respect to charging, the changes in shape are a result of ions in the Stern layer (i.e. ions at the electrode surface) having a higher electrostatic affinity to constant potential electrodes than to constant charge electrodes. For electrolyte polarizability, shape changes originate from induced dipoles that soften the interaction of Stern layer ions with the electrode. The softening is traced to ion correlations vertical to the electrode surface that induce dipoles that oppose double layer formation. In general, our analysis indicates an accuracy dependent differential capacitance profile that transitions from the characteristic camel shape with coarser representations to a more diffuse profile with finer representations.

Haskins, Justin B.↗
DOE OSTI2023

Scattering evidence of positional charge correlations in polyelectrolyte complexes

Polyelectrolyte complexation plays an important role in materials science and biology. The internal structure of the resultant polyelectrolyte complex (PEC) phase dictates properties such as physical state, response to external stimuli, and dynamics. Small-angle scattering experiments with X-rays and neutrons have revealed structural similarities between PECs and semidilute solutions of neutral polymers, where the total scattering function exhibits an Ornstein–Zernike form. In spite of consensus among different theoretical predictions, the existence of positional correlations between polyanion and polycation charges has not been confirmed experimentally. Here, we present small-angle neutron scattering profiles where the polycation scattering length density is matched to that of the solvent to extract positional correlations among anionic monomers. The polyanion scattering functions exhibit a peak at the inverse polymer screening radius of Coulomb interactions, q* ≈ 0.2 Å –1 . This peak, attributed to Coulomb repulsions between the fragments of polyanions and their attractions to polycations, is even more pronounced in the calculated charge scattering function that quantifies positional correlations of all polymer charges within the PEC. Screening of electrostatic interactions by adding salt leads to the gradual disappearance of this correlation peak, and the scattering functions regain an Ornstein–Zernike form. Experimental scattering results are consistent with those calculated from the random phase approximation, a scaling analysis, and molecular simulations.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

pH‐ and Redox‐Responsive Pickering Emulsions Based on Cellulose Nanocrystal Surfactants

Abstract Nanoparticle surfactants (NPSs), formed by using dynamic interactions between nanoparticles and complementary ligands at the liquid‐liquid interface, have emerged as “smart emulsifiers” with attributes of high emulsification efficiency, long‐term stability, and on‐demand emulsification/demulsification capabilities. However, only pH‐responsiveness can be adopted for the assembly of reported NPSs formed by electrostatic interactions. Here, we propose an alternative design strategy, by taking advantage of the ferrocenium (Fc + ) sulfate ion pair, to develop a new type of cellulose nanocrystal (CNC) surfactant. The Fc + groups are sensitive to pH, redox reagents and voltage, imparting the CNC surfactants and derived Pickering emulsions with multi‐stimuli‐responsiveness, and showing promising applications in controllable delivery, release, and biphasic biocatalysis.

Yang, Yang↗
NASA NTRS2001

Carbon Nanotubes in Water: MD Simulations of Internal and External Flow, Self Organization

We have developed computational tools, based on particle codes, for molecular dynamics (MD) simulation of carbon nanotubes (CNT) in aqueous environments. The interaction of CNTs with water is envisioned as a prototype for the design of engineering nano-devices, such as artificial sterocillia and molecular biosensors. Large scale simulations involving thousands of water molecules are possible due to our efficient parallel MD code that takes long range electrostatic interactions into account. Since CNTs can be considered as rolled up sheets of graphite, we expect the CNT-water interaction to be similar to the interaction of graphite with water. However, there are fundamental differences between considering graphite and CNTs, since the curvature of CNTs affects their chemical activity and also since capillary effects play an important role for both dynamic and static behaviour of materials inside CNTs. In recent studies Gordillo and Marti described the hydrogen bond structure as well as time dependent properties of water confined in CNTs. We are presenting results from the development of force fields describing the interaction of CNTs and water based on ab-initio quantum mechanical calculations. Furthermore, our results include both water flows external to CNTs and the behaviour of water nanodroplets inside heated CNTs. In the first case (external flows) the hydrophobic behaviour of CNTs is quantified and we analyze structural properties of water in the vicinity of CNTs with diagnostics such as hydrogen bond distribution, water dipole orientation and radial distribution functions. The presence of water leads to attractive forces between CNTs as a result of their hydrophobicity. Through extensive simulations we quantify these attractive forces in terms of the number and separation of the CNT. Results of our simulations involving arrays of CNTs indicate that these exhibit a hydrophobic behaviour that leads to self-organising structures capable of trapping water clusters. In the second case (internal flows) we study the behaviour of water droplets confined inside CNTs. Constant temperature simulations allow us to capture structural properties such as the contact angles and density profiles of the equilibrated drops. By heating and subsequently cooling of the CNT, we are able to measure the evaporation and the condensation rate of the entrapped water.

Jaffe, Richard L.↗
DOE OSTI2019

Human Glycerol 3-Phosphate Dehydrogenase: X-ray Crystal Structures That Guide the Interpretation of Mutagenesis Studies

Human liver glycerol 3-phosphate dehydrogenase (hlGPDH) catalyzes the reduction of dihydroxyacetone phosphate (DHAP) to form glycerol 3-phosphate, using the binding energy associated with the nonreacting phosphodianion of the substrate to properly orient the enzyme–substrate complex within the active site. Herein, we report the crystal structures for unliganded, binary E·NAD, and ternary E·NAD·DHAP complexes of wild type hlGPDH, illustrating a new position of DHAP, and probe the kinetics of multiple mutant enzymes with natural and truncated substrates. Mutation of Lys120, which is positioned to donate a proton to the carbonyl of DHAP, results in similar increases in the activation barrier to hlGPDH-catlyzed reduction of DHAP and to phosphite dianion-activated reduction of glycolaldehyde, illustrating that these transition states show similar interactions with the cationic K120 side chain. The K120A mutation results in a 5.3 kcal/mol transition state destabilization, and 3.0 kcal/mol of the lost transition state stabilization is rescued by 1.0 M ethylammonium cation. The 6.5 kcal/mol increase in the activation barrier observed for the D260G mutant hlGPDH-catalyzed reaction represents a 3.5 kcal/mol weakening of transition state stabilization by the K120A side chain and a 3.0 kcal/mol weakening of the interactions with other residues. The interactions, at the enzyme active site, between the K120 side chain and the Q295 and R269 side chains were likewise examined by double-mutant analyses. Here, these results provide strong evidence that the enzyme rate acceleration is due mainly or exclusively to transition state stabilization by electrostatic interactions with polar amino acid side chains.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Energy-efficient, Large-scale Molecular Dynamics Simulations via Hardware- and Algorithm-level Optimization

This work aims to develop a framework for energy-efficient computing that will enable molecular dynamics (MD) simulations of large-scale phenomena with atomic precision and simultaneously remove computational bottlenecks limiting the speed of MD simulations. We seek to implement such an approach through the development of surrogate models for the interatomic force calculation combined with the use of mixed numerical precision formats. For a model system of neutral atoms (only pairwise interactions), significant force calculation efficiency improvements were achieved, without detrimental effects on atomic structures or average energies, using single precision, by developing a surrogate model (deep neural network), and by quantizing this surrogate model. For a model system of charged atoms, the reciprocal-space calculation of electrostatic interactions was identified as the main bottleneck, and the development of a surrogate model should be pursued to achieve an estimated one-order-of-magnitude additional speedup.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2018

Charged-cell periodic DFT simulations via an impurity model based on density embedding: Application to the ionization potential of liquid water

Calculations of charged systems in periodic boundary conditions (PBC) are problematic because there are spurious interactions between the charges in different periodic images that can affect the physical picture. In addition, the intuitive limit of Coulomb interactions decaying to zero as the interacting charges are placed at infinite separation no longer applies, and for example total energies become undefined. Leveraging subsystem density functional theory (also known as density embedding) we define an impurity model that embeds a finite neutral or charged subsystem within an extended (infinite) surrounding subsystem. The combination of the impurity model and a consistent choice of the Coulomb reference provides us with an algorithm for evaluating the ionization potential (IP) in extended systems. Additionally, we demonstrate our approach in a pilot calculation of the IP of liquid water, based on a configuration from a prior ab initio molecular dynamics (AIMD) simulation of liquid water (Genova et al., J. Chem. Phys. 2016, 144, 234105). The calculations with the impurity model capture the broadening on the ionization energies introduced by the interactions between the water molecules. Furthermore, the calculated average IP value (10.5 eV) compare favorably to experiments (9.9-10.06 eV) and very recent simulations based on the GW approximation (10.55 eV), while at the same time outperforming density embedding calculations carried out with a naïve handling of the electrostatic interactions (about 7 eV).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗