Search NASA⌕ Search

SEARCH · Search NASA

Results for “PAIRING INTERACTIONS”

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 1,693 records · Page 94

Phosphorus limits nitrogen dynamics in the O horizon of a forested watershed in Maine, USA

Anthropogenic emissions have altered the N status of northeastern USA forests, highlighting the importance of our imperfect understanding of the interaction between nitrogen and phosphorus in response to atmospheric N deposition. We investigated soil N and P dynamics at the Bear Brook Watershed in Maine, a long-term paired watershed experiment established to study ecosystem response to experimental N enrichment. One watershed (West Bear, WB) was treated bimonthly for 25 years with ammonium sulfate, and the reference East Bear (EB) received only ambient N deposition. We conducted field and laboratory studies to examine the interaction between N and P in O horizon soils in both EB and WB. In the “Pulse-Chase” Study, a field addition of 100 kg P ha-1 resulted in an overall 38% decrease in extractable ammonium (NH4+-N), due to increased immobilization. Noteworthy was that even under ambient N deposition rates, EB displayed parallel, although lower, evidence of P limitation compared to WB. In contrast, in a “Response Surface” laboratory study, soils incubated in the laboratory at rates of 0 to 200 kg P ha-1 showed increases in PNNM and NH4+-N concentrations in response to P additions. There was no effect of N-status in the Response Surface Study, but hardwood soils showed a faster response and greater rates of potential net N mineralization (PNNM) compared to softwood soils. Despite contrasting responses, both the field and laboratory studies suggested that N dynamics at BBWM were P limited, and that N status and forest type influenced the N and P interaction.

Salvino, Cayce↗

Centrally Concentrated X-Ray Radiation from an Extended Accreting Corona in Active Galactic Nuclei

The X-ray emission from bright active galactic nuclei (AGNs) is believed to originate in a hot corona lying above a cold, geometrically thin accretion disk. A highly concentrated corona located within ∼10 gravitational radii above the black hole is inferred from observations. Based on the accretion of interstellar medium/wind, a disk corona model has been proposed in which the corona is well coupled to the disk by radiation, thermal conduction, as well as by mass exchange. Such a model avoids artificial energy input to the corona and has been used to interpret the spectral features observed in AGN. In this work, it is shown that the bulk emission size of the corona is very small for the extended accretion flow in our model. More than 80% of the hard X-ray power is emitted from a small region confined within 10 Schwarzschild radii around a non-spinning black hole, which is expected to be even smaller accordingly for a spinning black hole. Here, the corona emission is more extended at higher Eddington ratios. The compactness parameter of the corona, l=(L/R)((σ{sub T})/(m{sub e}c{sup 3})), is shown to be in the range of 1–33 for Eddington ratios of 0.02–0.1. Combined with the electron temperature in the corona, this indicates that electron–positron pair production is not dominant in this regime. A positive relation between the compactness parameter and photon index is also predicted. By comparing the above model predictions with observational features, we find that the model is in agreement with observations.

79 ASTRONOMY AND ASTROPHYSICS↗

Graph network heterogeneity predicts interplant wake losses

Wind plants generate large-scale wakes, which can affect the performance of neighboring installations. Such wakes are challenging to estimate due to the inherent complexity in modeling wake interactions between large quantities of turbines at various distances. Weighted directed graph networks can inform complex models by linking turbine pairs into chains of upstream and downstream neighbors for a given wind direction. A novel interpretation of the graph network adjacency matrix is proposed where each element of the matrix represents the cumulative impact of upstream turbines on an individual. In this study, wake losses were estimated with an engineering wake model across a range of inflow conditions for nine parametric variations of a system containing two neighboring wind plants. The parametric nature of the study isolates turbine spacing within the plant, separation distance between plants, and wind direction as the main drivers of wake losses. Spatial heterogeneity is computed from the weighted average adjacency matrix of each plant arrangement. The proposed method is orders of magnitude faster than wake modeling and does not require detailed turbine information or atmospheric conditions. Furthermore, the weighted average adjacency matrix provides insight on the spatial organization of wake losses at various scales. Plant heterogeneity is correlated with wake losses within and among plants. Framing wind plant wake interaction in terms of graph network spatial heterogeneity provides an efficient approach for predicting wake losses within and among neighboring wind plants with applications to other complex systems where wake interactions are key factor.

17 WIND ENERGY↗

Microscopic origin of reflection-asymmetric nuclear shapes

The presence of nuclear ground states with stable reflection-asymmetric shapes is supported by rich experimental evidence. Theoretical surveys of odd-multipolarity deformations predict the existence of pear-shaped isotopes in several fairly localized regions of the nuclear landscape in the vicinity of near-lying single-particle shells with Δℓ = Δj = 3. We analyze the role of isoscalar, isovector, neutron-proton, neutron-neutron, and proton-proton multipole interaction energies in inducing the onset of reflection-asymmetric ground-state deformations. The calculations are performed in the framework of axial reflection-asymmetric Hartree-Fock-Bogoliubov theory using two Skyrme energy-density functionals and density-dependent pairing force. We show that reflection-asymmetric ground-state shapes of atomic nuclei are driven by the odd-multipolarity neutron-proton (or isoscalar) part of the nuclear interaction energy. This result is consistent with the particle-vibration picture, in which the main driver of octupole instability is the isoscalar octupole-octupole interaction giving rise to large E3 polarizability. The necessary condition for the appearance of localized regions of pear-shaped nuclei in the nuclear landscape is the presence of parity doublets involving Δℓ = Δj = 3 proton or neutron single-particle shells. This condition alone is, however, not sufficient to determine whether pear shapes actually appear, and, if so, what the corresponding reflection-asymmetric deformation energies are. The predicted small reflection-asymmetric deformation energies result from dramatic cancellations between even- and odd-multipolarity components of the nuclear binding energy.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Concentration Dependent Interfacial Chemistry of the NaOH(aq):Gibbsite Interface

Caustic conditions are often employed for dissolution of a wide variety of minerals, where ion sorption, surface diffusion, and interfacial organization impact surface reactivity. In the case of gibbsite, gamma-Al(OH)3, the chemistry at the NaOH(aq) interface is deeply intertwined with industrial processing of aluminum, including metal production and the disposition of Al-containing wastes. To date, little is known about the structure, speciation, and dynamic behavior of gibbsite interfaces (and that of many other minerals) with NaOH(aq)—particularly as a function of ionic strength. Yet concentration-dependent interfacial organization and dynamics are a critical starting point to develop a fundamental understanding of the factors that influence dissolution. This work reports equilibrium molecular dynamics simulations of the gamma-Al(OH)3:NaOH(aq) interface, revealing the sorption behavior and speciation of ions from 0.5–10 M [NaOH]. As inner-sphere complexes, Na+ primarily coordinates to the side of the gibbsite hexagonal cavities, while OH accepts hydrogen-bonding from the surface-OH groups. The mobility of inner-sphere Na+ and OH ions is significantly reduced due to a strong surface affinity in comparison to previous reports of NaCl, CaCl2, or BaCl2 electrolytes. At high [NaOH], contact ion pairing that is observed in the bulk solution is partially disrupted upon sorption to the gibbsite surface by the individual ion–surface interactions. The molecular-scale changes to surface speciation and competition between ion–surface vs. ion–ion interactions influence surface characterization of gibbsite and potential dissolution processes, providing a valuable baseline for starting conditions needed within future reactive molecular simulations.

alkaline solutions, IDREAM, gibbsite, molecular dy↗

Search for contact interactions and large extra dimensions in the dilepton mass spectra from proton-proton collisions at $ \sqrt{s}=13 $ TeV

A search for nonresonant excesses in the invariant mass spectra of electron and muon pairs is presented. The analysis is based on data from proton-proton collisions at a center-of-mass energy of 13 TeV recorded by the CMS experiment in 2016, corresponding to a total integrated luminosity of 36 fb$^{–1}$. No significant deviation from the standard model is observed. Limits are set at 95% confidence level on energy scales for two general classes of nonresonant models. For a class of fermion contact interaction models, lower limits ranging from 20 to 32 TeV are set on the characteristic compositeness scale Λ. For the Arkani-Hamed, Dimopoulos, and Dvali model of large extra dimensions, the first results in the dilepton final state at 13 TeV are reported, and values of the ultraviolet cutoff parameter Λ$_{T}$ below 6.9 TeV are excluded. A combination with recent CMS diphoton results improves this exclusion to Λ$_{T}$ below 7.7 TeV, providing the most sensitive limits to date in nonhadronic final states.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Geometrically frustrated interactions drive structural complexity in amorphous calcium carbonate

Amorphous calcium carbonate is an important precursor for biomineralization in marine organisms. Key outstanding problems include understanding the structure of amorphous calcium carbonate and rationalizing its metastability as an amorphous phase. Here we report high-quality atomistic models of amorphous calcium carbonate generated using state-of-the-art interatomic potentials to help guide fits to X-ray total scattering data. Exploiting a recently developed inversion approach, we extract from these models the effective Ca∙∙∙Ca interaction potential governing the structure. This potential contains minima at two competing distances, corresponding to the two different ways that carbonate ions bridge Ca 2+ -ion pairs. We reveal an unexpected mapping to the Lennard-Jones–Gauss model normally studied in the context of computational soft matter. The empirical model parameters for amorphous calcium carbonate take values known to promote structural complexity. We thus show that both the complex structure and its resilience to crystallization are actually encoded in the geometrically frustrated effective interactions between Ca 2+ ions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Phonon screening and dissociation of excitons at finite temperatures from first principles

The properties of excitons, or correlated electron–hole pairs, are of paramount importance to optoelectronic applications of materials. A central component of exciton physics is the electron–hole interaction, which is commonly treated as screened solely by electrons within a material. However, nuclear motion can screen this Coulomb interaction as well, with several recent studies developing model approaches for approximating the phonon screening of excitonic properties. While these model approaches tend to improve agreement with experiment, they rely on several approximations that restrict their applicability to a wide range of materials, and thus far they have neglected the effect of finite temperatures. Here, we develop a fully first-principles, parameter-free approach to compute the temperature-dependent effects of phonon screening within the ab initio GW -Bethe–Salpeter equation framework. We recover previously proposed models of phonon screening as well-defined limits of our general framework, and discuss their validity by comparing them against our first-principles results. We develop an efficient computational workflow and apply it to a diverse set of semiconductors, specifically AlN, CdS, GaN, MgO, and SrTiO 3 . We demonstrate under different physical scenarios how excitons may be screened by multiple polar optical or acoustic phonons, how their binding energies can exhibit strong temperature dependence, and the ultrafast timescales on which they dissociate into free electron–hole pairs.

Science & Technology - Other Topics↗

Room‐Temperature Formate Ester Transfer Hydrogenation Enables an Electrochemical/Thermal Organometallic Cascade for Methanol Synthesis from CO 2

Abstract The reduction of CO 2 to synthetic fuels is a valuable strategy for energy storage. However, the formation of energy‐dense liquid fuels such as methanol remains rare, particularly under low‐temperature and low‐pressure conditions that can be coupled to renewable electricity sources via electrochemistry. Here, a multicatalyst system pairing an electrocatalyst with a thermal organometallic catalyst is introduced, which enables the reduction of CO 2 to methanol at ambient temperature and pressure. The cascade methanol synthesis proceeds via CO 2 reduction to formate by electrocatalyst [Cp*Ir(bpy)Cl] + (Cp*=pentamethylcyclopentadienyl, bpy=2,2′‐bipyridine), Fischer esterification of formate to isopropyl formate catalyzed by trifluoromethanesulfonic acid (HOTf), and thermal transfer hydrogenation of isopropyl formate to methanol facilitated by the organometallic catalyst (H‐PNP)Ir(H) 3 (H‐PNP=HN(C 2 H 4 P i Pr 2 ) 2 ). The isopropanol solvent plays several crucial roles: activating formate ion as isopropyl formate, donating hydrogen for the reduction of formate ester to methanol via transfer hydrogenation, and lowering the barrier for transfer hydrogenation through hydrogen bonding interactions. In addition to reporting a method for room‐temperature reduction of challenging ester substrates, this work provides a prototype for pairing electrochemical and thermal organometallic reactions that will guide the design and development of multicatalyst cascades.

Fernández, Sergio [Department of Chemistry Univers↗

Temperature-dependent changes in gas chromatographic separation metrics for trihexyl(tetradecyl)phosphonium-based ionic liquid stationary phases and comparison to conventional polysiloxane stationary phases

Here, to assess whether the trihexyl(tetradecyl)phosphonium chloride ([P 66614 + ][Cl - ]) ionic liquid (IL), employed as a gas chromatographic stationary phase, exhibits temperature-dependent phase transitions or structural heterogeneity, retention thermodynamics for probe molecules of varied chemical structure were investigated by gas chromatography (GC) over the temperature range of 25 to 64 °C using van’t Hoff analysis. All van’t Hoff plots for n-alkanes (C8–C12) and other probe analytes were linear, indicating the absence of detectable phase-transition behavior. The [P 66614 + ][Cl - ] IL was further compared with the [P 66614 + ] bis[(trifluoromethyl)sulfonyl]imide ([P 66614 + ][NTf 2 - ]) IL and commercial polysiloxane stationary phases in terms of temperature-dependent resolution behavior and chromatographic selectivity over the temperature range of 28–52 °C. The resolution factors of analytes in mixtures obtained using gas chromatography–mass spectrometry (GC–MS) and the summation of resolution values decreased smoothly with temperature on both IL stationary phases and on the dimethylphenylcyano-substituted polymeric stationary phase (OV-1701), whereas the poly(50% diphenyl/50% dimethyl siloxane) stationary phase (SPB-50) exhibited irregularities in its temperature-dependent summation of resolution profile, showing a distinct maximum in the temperature range of 34–40 °C. Selectivity values for the IL stationary phases were largely determined by the anion, with the [P 66614 + ][Cl - ] IL providing enhanced separation between polar and aromatic analytes, while the [P 66614 + ][NTf 2 - ] IL favored dispersive and π-π interactions that resulted in a reversal of elution order for some aliphatic–aromatic and alkyne-containing analyte pairs. Capillary columns prepared using different surface activation methods and coated with identical [P 66614 + ][Cl - ] IL stationary phase showed indistinguishable temperature-dependent trends in summation of resolution value and selectivity. An evaluation of analyte-specific selectivity provides insight into the characteristic separation behavior of the stationary phases and enables identification of analyte classes that are more efficiently separated on the [P 66614 + ]-based IL stationary phases compared to conventional stationary phase materials.

Chromatographic selectivity↗

Supported bicyclic amidine ionic liquids as a potential CO 2 /N 2 separation medium

In this paper, five new room temperature ionic liquids were tested as supported ionic liquid membranes (SILMs) for CO 2 /N 2 separation capability. This series of ionic liquids (ILs) contains bicyclic amidine cations of 1,8-diazabicyclo[5.4.0]undec-7-ene (DBU) or 1,5-diazabicyclo[4.3.0]non-5-ene (DBN) with small alkyl chain substituents. These ILs have been synthesized to examine the effects of cation structure on free volume. Relevant physical and chemical properties have been examined, including CO 2 solubility and viscosity. The relationship between fractional free volume and permeability of gases is explored, with an interesting note on solubility. Each cation was paired with two different anions: bis(trifluoromethanesulfonyl)imide and tetracyanoborate, for further investigation of the cation-anion interaction on SILM permeability and selectivity. Finally, three of these new ILs have exceeded Robeson's upper bound for CO 2 /N 2 gas separations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Distinctive Trend of Metal Binding Affinity via Hydration Shell Breakage in Nanoconfined Cavity

Here, we report a distinctive trend in the binding affinity of metal cations when they interact with lacunary Keggin-type polyoxometalate K 7 [α-PW 11 O 39 ] clusters in aqueous media. The binding affinity of metal cations follows the special trends of Li + > Na + > K + for monovalent cations and Mg 2+ > Ca 2+ > Sr 2+ > Ba 2+ for divalent cations, completely reversed from the common trends of Li + < Na + < K + and Mg 2+ < Ca 2+ < Sr 2+ < Ba 2+ for macroion–counterion ion-pair formation in dilution solution. Isothermal titration calorimetry and NMR measurements confirm that the hydration shells of metal cations are significantly broken during the interaction in the nanoscale cavity, making the binding affinity stronger for the cations with smaller ionic sizes instead of hydrated sizes. Furthermore, the same phenomenon is observed for the lacunary Dawson-type clusters, which contain well-defined cavity in their structures as well, demonstrating the universality of this reverse trend and its connection to nanoconfinement effect in small submolecular cavities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Bridging Atomic Solvation Environment with Electrochemical Properties for the Bis(trifluoromethylsulfonyl)imide-Based Divalent Cation Electrolytes for the Next-Generation Energy Storage Systems

A deep molecular-level understanding of the multivalent electrolyte and its correlation with the electrochemical properties is crucial for designing optimized electrolytes for next-generation rechargeable batteries. Comprehensive knowledge of the atomic level of the solvation structure and its connection with electrochemical stability and ion transport properties is especially critical. However, the interaction of these three components coupled with clear atomistic insights is lacking in the literature. Here, our current contribution evaluates representative electrolytes with the bis(trifluoromethanesulfonyl)imide (TFSI) anions for multivalent cations of Mg, Ca, and Zn, at different ionic conditions with and without a cosolvated environment in ether-based solvent. Two critical problems are investigated: first, resolving the solvation structures in the electrolyte solutions as a function of concentrations through pair distribution function analysis and the corresponding electrochemical transport properties; second, unmasking the quantitative correlation of the atomistic environment with both electrochemical kinetics and cation dependence. We discovered that the magnesium- and calcium-based electrolytes display versatile coordination lengths but poor average anodic stability due to ion pairing with TFSI - . On the contrary, the zinc-based electrolytes show the shortest solvent coordination lengths, shielding the Zn cation from rigid solvent interactions and resulting in the highest anodic stabilities. Calcium-based electrolytes exhibit the longest and most concentration-independent coordination lengths. This work provides valuable insights into the molecular structural and electrochemical features of diverse multivalent electrolyte systems with cations in various solvation environments, emphasizing the importance of the solvation structure and construction in designing high-performance electrolytes.

cation coordination↗

Solubility Model for Ferrous Iron Hydroxide, Hibbingite, Siderite, and Chukanovite in High Saline Solutions of Sodium Chloride, Sodium Sulfate, and Sodium Carbonate

Here, a solubility model is presented for ferrous iron hydroxide (Fe(OH) 2 (s)), hibbingite (Fe 2 Cl(OH) 3 (s)), siderite (FeCO 3 (s)), and chukanovite (Fe 2 CO 3 (OH) 2 (s)). The Pitzer activity coefficient equation was utilized in developing the model to account for the excess free energies of aqueous species in the background solutions of high ionic strength. Solubility limiting minerals were analyzed before and after experiments using X-ray diffraction. Formation of Fe(OH) 2 (s) was observed in the experiments that were initiated with Fe 2 Cl(OH) 3 (s) in Na 2 SO 4 solution. Coexistence of siderite and chukanovite was observed in the experiments in Na 2 CO 3 + NaCl solutions. Two equilibrium constants that had been reported by us for the dissolution of Fe(OH) 2 (s) and Fe 2 Cl(OH) 3 (s) (Nemer et al.) were rederived in this paper, using newer thermodynamic data selected from the literature to maintain internal consistency of the series of our data analyses in preparation, including this paper. Three additional equilibrium constants for the following reactions were determined in this paper: dissolution of siderite and chukanovite and dissociation of the aqueous species Fe(CO 3 ) 2 –2 . Five Pitzer interaction parameters were derived in this paper: β (0) , β (1) , and C φ parameters for the species pair Fe +2 /SO 4 –2 ; β (0) and β (1) parameters for the species pair Na+/Fe(CO3)2–2. Our model predicts that, among the four inorganic ferrous iron minerals, siderite is the stable mineral in two WIPP-related brines (WIPP: Waste Isolation Pilot Plant), i.e., GWB and ERDA6 (Brush and Domski), and the electrochemical equilibrium between elemental iron and siderite provides a low oxygen fugacity (10 –91.2 atm) that can keep the actinides at their lowest oxidation states. (Nemer et al., Brush and Domski; references numbered 1 and 2 in the main text).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effect of Protein Incorporation on the Nanostructure of the Bicontinuous Microemulsion Phase of Winsor-III Systems: A Small-Angle Neutron Scattering Study

Small-angle neutron scattering (SANS) analysis using the Teubner₋Strey model has been employed to evaluate the effect of protein incorporation into the middle, bicontinuous microemulsion (BμE) phase of Winsor-III (W III ) systems formed by an aerosol-OT (AOT)/alkyl ethoxylate mixed surfactant system to understand better the extraction of proteins into and out of BμEs and to study the effect of proteins on a system that serves as a biomimetic analog of cell membranes. Under conditions of high salinity, the incorporation of positively charged proteins cytochrome c, lysozyme, and α-chymotrypsin, near their solubilization limit in the BμEs promoted the release of water and oil from the BμEs, a decrease in the quasi-periodic repeat distance (d), an increase in ordering (a decrease in the amphiphilicity factor, f a ) for the surfactant monolayers, and a decrease in the surface area per surfactant headgroup, suggesting that the proteins affected the self-assembly of components in the BμE phase and produced Debye shielding of AOTs sulfonate headgroup. For W III systems possessing lower salinity, cytochrome c reduced the efficiency of surfactant in the BμE phase, noted by increases in d and f a , suggesting that the enzyme and AOT underwent ion pairing. We find that the results of this study demonstrate the importance of ionic strength to modulate proteinsurfactant interactions, which in turn will control the release of proteins encapsulated in the BμEs, relevant to W III -based protein extraction and controlled release from BμE delivery systems, and demonstrate the utility of BμEs as a model system to understand the effect of proteins on biomembranes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

A compendium of multi-omic sequence information from the Saanich Inlet water column

Marine oxygen minimum zones (OMZs) are widespread regions of the ocean that are currently expanding due to global warming. While inhospitable to most metazoans, OMZs are hotspots for microbial mediated biogeochemical cycling of carbon, nitrogen and sulphur, contributing disproportionately to marine nitrogen loss and climate active trace gas production. Our current understanding of microbial community responses to OMZ expansion is limited by a lack of time-resolved data sets linking multi-omic sequence information (DNA, RNA, protein) to geochemical parameters and process rates. Here, we present six years of time-resolved multi-omic observations in Saanich Inlet, a seasonally anoxic fjord on the coast of Vancouver Island, British Columbia, Canada that undergoes recurring changes in water column oxygenation status. This compendium provides a unique multi-omic framework for studying microbial community responses to ocean deoxygenation along defined geochemical gradients in OMZ waters.

Hawley, Alyse K. (ORCID:0000000235164043)↗

Vibrational probe of the origin of singlet exciton fission in TIPS-pentacene solutions

We use native vibrational modes of the model singlet fission chromophore 6,13-bis(triisopropylsilylethynyl)pentacene (TIPS-Pn) to examine the origins of singlet fission in solution between molecules that are not tethered by a covalent linkage. We use the C—H stretch modes of TIPS side groups of TIPS-Pn to demonstrate that singlet fission does not occur by diffusive encounter of independent molecules in solution. Instead, TIPS-Pn molecules aggregate in solution through their TIPS side groups. This aggregation breaks the symmetry of the TIPS-Pn molecules and enables the formation of triplets to be probed through the formally symmetry forbidden symmetric alkyne stretch mode of the TIPS side groups. The alkyne stretch modes of TIPS-Pn are sensitive to the electronic excited states present during the singlet fission reaction and provide unique signatures of the formation of triplets following the initial separation of triplet pair intermediates. These findings highlight the opportunity to leverage structural information from vibrational modes to better understand intermolecular interactions that lead to singlet fission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effective charges and virial pressure of concentrated macroion solutions

The stability of colloidal suspensions is crucial in a wide variety of processes, including the fabrication of photonic materials and scaffolds for biological assemblies. The ionic strength of the electrolyte that suspends charged colloids is widely used to control the physical properties of colloidal suspensions. The extensively used two-body Derjaguin-Landau-Verwey-Overbeek (DLVO) approach allows for a quantitative analysis of the effective electrostatic forces between colloidal particles. DLVO relates the ionic double layers, which enclose the particles, to their effective electrostatic repulsion. Nevertheless, the double layer is distorted at high macroion volume fractions. Therefore, DLVO cannot describe the many-body effects that arise in concentrated suspensions. In this paper, we show that this problem can be largely resolved by identifying effective point charges for the macroions using cell theory. This extrapolated point charge (EPC) method assigns effective point charges in a consistent way, taking into account the excluded volume of highly charged macroions at any concentration, and thereby naturally accounting for high volume fractions in both salt-free and added-salt conditions. We provide an analytical expression for the effective pair potential and validate the EPC method by comparing molecular dynamics simulations of macroions and monovalent microions that interact via Coulombic potentials to simulations of macroions interacting via the derived EPC effective potential. The simulations reproduce the macroion-macroion spatial correlation and the virial pressure obtained with the EPC model. Finally, our findings provide a route to relate the physical properties such as pressure in systems of screened Coulomb particles to experimental measurements.

36 MATERIALS SCIENCE↗