Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interactions and forces”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Impact of Ligands on the Ion–Ion Selectivity of Ligand-Appended-Pillar[ n ]arene Channels

Ligand-appended pillar[5]arene (LAP) channels are emerging as a promising platform for ion-selective membranes. In this study, we investigated the ion selectivities of three LAP channels functionalized with distinct ligands: diglycolamine, carbamoylmethyl phosphine oxide, and propionamide phosphonic acid. These ligands were chosen based on their demonstrated affinities for lanthanides in solvent extraction studies. We examined the selectivity of each channel toward a series of monovalent, divalent, and trivalent ions: Li + , Na + , K + , Mg 2+ , Ca 2+ , La 3+ , Eu 3+ , and Yb 3+ . To quantify ion-channel interactions and the energetics of translocation, we computed the potential of mean force (PMF) profiles for each ion. These PMFs were subsequently recast into permeability coefficients, enabling a direct evaluation of ion permselectivity. Our results reveal distinct selectivity patterns governed by the chemical nature of the appended ligands. The diglycolamine-functionalized channel exhibits strong interactions with monovalent ions, leading to replication of the Li + bulk hydration shell within the channel. This promotes preferential transport Li + ions while effectively excluding divalent and trivalent species. In contrast, the carbamoylmethyl phosphine oxide-functionalized channel displays a reduced selectivity between monovalent ions while similarly exhibiting an effective rejection of divalent and trivalent ions. This arises from steric hindrance around channel-lining oxygens imposed by the ligand's bulky phenyl groups which limits heavy ion coordination and promotes retention of hydration shells effectively increasing effective ion size. The propionamide phosphonic acid-functionalized channel exhibits a different trend, favoring ions with lower hydration free energies, with K + ions showing enhanced permeation. These findings highlight the critical role of ligand-ion interactions in modulating ion selectivity within LAP channels. In particular, the strong coordination exhibited by diglycolamine and propionamide phosphonic acid ligands significantly influences the ion transport energetics and selectivity profiles.

hydration↗
DOE OSTI2024

Understanding aerosol–cloud interactions using a single-column model for a cold-air outbreak case during the ACTIVATE campaign

Marine boundary layer clouds play a critical role in Earth's energy balance. Their microphysical and radiative properties are highly impacted by ambient aerosols and dynamic forcings. In this study, we evaluate the representation of these clouds and related aerosol–cloud interaction processes in the single-column version of the E3SM climate model (SCM) against field measurements collected during the NASA ACTIVATE campaign over the western North Atlantic, as well as intercompare results with high-resolution process level models. We show that E3SM SCM reproduces the macrophysical properties of post-frontal boundary layer clouds in a cold-air outbreak (CAO) case well. However, it generates fewer but larger cloud droplets compared to aircraft measurements. Further sensitivity tests show that the underestimation of both aerosol number concentration and vertical velocity variance contributes to this bias. Aerosol–cloud interactions are examined by perturbing prescribed aerosol properties in E3SM SCM with fixed dynamics. Higher aerosol number concentration or hygroscopicity leads to more numerous but smaller cloud droplets, resulting in a stronger cooling via shortwave cloud forcing. This apparent Twomey effect is consistent with prior climate model studies. The cloud liquid water path shows a weakly positive relation with cloud droplet number concentration due to precipitation suppression. This weak aerosol effect on cloud macrophysics may be attributed to the dominant impact of strong dynamical forcing associated with the CAO. Our findings indicate that the SCM framework is a key tool to bridge the gap between climate models, process level models, and field observations to facilitate process level understanding.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

HIFiRE-1 Turbulent Shock Boundary Layer Interaction - Flight Data and Computations

The Hypersonic International Flight Research Experimentation (HIFiRE) program is a hypersonic flight test program executed by the Air Force Research Laboratory (AFRL) and Australian Defence Science and Technology Organisation (DSTO). This flight contained a cylinder-flare induced shock boundary layer interaction (SBLI). Computations of the interaction were conducted for a number of times during the ascent. The DPLR code used for predictions was calibrated against ground test data prior to exercising the code at flight conditions. Generally, the computations predicted the upstream influence and interaction pressures very well. Plateau pressures on the cylinder were predicted well at all conditions. Although the experimental heat transfer showed a large amount of scatter, especially at low heating levels, the measured heat transfer agreed well with computations. The primary discrepancy between the experiment and computation occurred in the pressures measured on the flare during second stage burn. Measured pressures exhibited large overshoots late in the second stage burn, the mechanism of which is unknown. The good agreement between flight measurements and CFD helps validate the philosophy of calibrating CFD against ground test, prior to exercising it at flight conditions.

Turbulence↗
DOE OSTI2022

Using aircraft measurements to characterize subgrid-scale variability of aerosol properties near the Atmospheric Radiation Measurement Southern Great Plains site

Abstract. Complex distributions of aerosol properties evolve in space and time as a function of emissions, new particle formation, coagulation, condensational growth, chemical transformation, phase changes, turbulent mixing and transport, removal processes, and ambient meteorological conditions. The ability of chemical transport models to represent the multi-scale processes affecting the life cycle of aerosols depends on their spatial resolution since aerosol properties are assumed to be constant within a grid cell. Subgrid-scale-dependent processes that affect aerosol populations could have a significant impact on the formation of particles, their growth to cloud condensation nuclei (CCN) sizes, aerosol–cloud interactions, dry deposition and rainout and hence their burdens, lifetimes, and radiative forcing. To address this issue, we characterize subgrid-scale variability in terms of measured aerosol number, size, composition, hygroscopicity, and CCN concentrations made by repeated aircraft flight paths over the Atmospheric Radiation Measurement (ARM) program's Southern Great Plains (SGP) site during the Holistic Interactions of Shallow Clouds, Aerosols and Land Ecosystem (HI-SCALE) campaign. Subgrid variability is quantified in terms of both normalized frequency distributions and percentage difference percentiles using grid spacings of 3, 9, 27, and 81 km that represent those typically used by cloud-system-resolving models as well as the current and next-generation climate models. Even though the SGP site is a rural location, surprisingly large horizontal gradients in aerosol properties were frequently observed. For example, 90 % of the 3, 9, and 27 km cell mean organic matter concentrations differed from the 81 km cell around the SGP site by as much as ∼ 46 %, large spatial variability in aerosol number concentrations and size distributions were found during new particle formation events, and consequently 90 % of the 3, 9, and 27 km cell mean CCN number concentrations differed from the 81 km cell mean by as much as ∼ 38 %. The spatial variability varied seasonally for some aerosol properties, with some having larger spatial variability during the spring and others having larger variability during the late summer. While measurements at a single surface site cannot reflect the surrounding variability of aerosol properties at a given time, aircraft measurements that are averaged within an 81 km cell were found to be similar to many, but not all, aerosol properties measured at the ground SGP site. This analysis suggests that it is reasonable to directly compare most ground SGP site aerosol measurements with coarse global climate model predictions. In addition, the variability quantified by the aircraft can be used as an uncertainty range when comparing the surface point measurements with model predictions that use coarse grid spacings.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Scale translation yields insights into gas adsorption under nanoconfinement

This work describes a scale-translating simulation framework to investigate gas adsorption behavior in nanoconfined pores. The framework combines molecular simulations (MSs), equation of state (EoS), and lattice Boltzmann (LB) simulations. MSs reveal the physics of methane adsorption in nano-sized pores, where input values of fugacity coefficients are optimized based on EoS predictions. Then, an LB free-energy model, which incorporates a viral EoS, upscales intermolecular forces and estimates adsorption behavior via a proposed fluid–wall interaction model. Armed with the values of the LB interaction parameter as a function of pressure, the LB model is used to predict fluid behavior in irregular nanopores, and the results are validated against reference MS data. The LB model is then used to study adsorption behavior at a continuum scale in representative organic shale nanopores based on finely characterized Vaca Muerta shale samples. Furthermore, the results show that methane adsorption could significantly increase contained fluids by 10%–25% in pores smaller than 20 nm. However, in larger pores (40 nm to 90 nm), adsorption's impact diminishes to 2%–3%, suggesting sorption's negligible role beyond a 40 nm pore size.

74 ATOMIC AND MOLECULAR PHYSICS↗
NASA NTRS1993

Observations of nonlinear interactions in directionally spread shoaling surface gravity waves

Shoaling wave fields generated in laboratory experiments were analyzed to determine the sensitivity of nonlinear interactions to the directional distributions of incident waves. Peaks in the directional spectra observed in shallow water were consistent with near resonating, quadratic interactions between two primary waves transferring energy to a third wave with the sum frequency and vector sum wavenumber of the primary waves. Directionally colinear waves forced a higher-frequency wave propagating in the same direction as the primary waves, while directionally spread (i.e., noncolinear) primary waves forced a higher-frequency wave that propagated in a direction between those of the interacting primary waves. Deepwater wave fields with similar frequency spectra but different directional spectra evolved to different shallow-water directional spectra, yet their shallow-water frequency spectra were remarkably similar. This result suggests that the shape of the directional spectrum of the incident wave field has only a small effect on the magnitudes of nonlinear energy transfer during shoaling. The principal effect of directionality in the incident wave field is on the directions, not the amplitudes, of the nonlinearly generated waves. The laboratory data demonstrate clearly the importance of triad interactions between noncolinear and colinear shoaling waves.

Elgar, Steve↗
DOE OSTI2025

The Juno mission as a probe of long-range new physics

Orbits of celestial objects, especially the geocentric and heliocentric ones, have been well explored to constrain new long-range forces beyond the Standard Model (SM), often referred to as fifth forces. In this paper, for the first time, we apply the motion of a spacecraft around Jupiter to probe fifth forces that don’t violate the equivalence principle. The spacecraft is the Juno orbiter, and ten of its early orbits already allow a precise determination of the Jovian gravitational field. We use the shift in the precession angle as a proxy to test non-gravitational interactions between Juno and Jupiter. Requiring that the contribution from the fifth force does not exceed the uncertainty of the precession shift inferred from data, we find that a new parameter space with the mass of the fifth-force mediator around 10 −14 eV is excluded at 95% C.L.

new light particles↗
NASA NTRS1994

Climatology and natural variability of the global hydrologic cycle in the GLA atmospheric general circulation model

Time average climatology and low-frequency variabilities of the global hydrologic cycle (GHC) in the Goddard Laboratory for Atmospheres (GLA) general circulation model (GCM) were investigated in the present work. A 730-day experiment was conducted with the GLA GCM forced by insolation, sea surface temperature, and ice-snow undergoing climatological annual cycles. Ifluences of interactive soil moisture on time average climatology and natural variability of the GHC were also investigated by conducting 365-day experiments with and without interactive soil moisture. Insolation, sea surface temperature, and ice-snow were fixed at their July levels in the latter two experiments. Results show that the model's time average hydrologic cycle variables for July in all three experiments agree reasonably well with observations. Except in the case of precipitable water, the zonal average climates of the annual cycle experiment and the two perpetual July experiments are alike, i.e., their differences are within limits of the natural variability of the model's climate. Statistics of various components of the GHC, i.e., water vapor, evaporation, and precipitation, are significantly affected by the presence of interactive soil moisture. A long-term trend is found in the principal empirical modes of variability of ground wetness, evaporation, and sensible heat. Dominant modes of variability of these quantities over land are physically consistent with one another and with land surface energy balance requirements. The dominant mode of precipitation variability is found to be closely related to organized convection over the tropical western Pacific Ocean. The precipitation variability has timescales in the range of 2 to 3 months and can be identified with the stationary component of the Madden-Julian Oscillation. The precipitation mode is not sensitive to the presence of interactive soil moisture but is closely linked to both the rotational and divergent components of atmospheric moisture transport. The present results indicate that globally coherent natural variability of the GHC in the GLA GCM has two basic timescales in the absence of annual cycles of external forcings: a long-term trend associated with atmosphere-soil moisture interaction which affects the model atmosphere mostly over midlatitude continental regions and a large-scale 2- to 3-month variability associated with atmospheric moist processes over the western Pacific Ocean.

Lau, K.-M.↗
DOE OSTI2022

Dispersive interactions between standard and Dirac materials and the role of dimensionality

The van der Waals (vdW) interaction plays a prominent role between neutral objects at separations where short ranged chemical forces are negligible. This type of dispersive coupling is determined by the interplay between geometry and response properties of the materials making up the objects. Here, we investigate the vdW interaction between 1D, 2D, and 3D standard and Dirac materials within the Random Phase Approximation, which takes into account collective excitations originating from the electronic Coulomb potential. A comprehensive understanding of characteristic functionalities and scaling laws are obtained for systems with parabolic energy dispersion (standard materials) and crossing linear bands (Dirac materials). By comparing the quantum mechanical and thermal limits the onset of thermal fluctuations in the vdW interaction is discussed showing that thermal effects are significantly pronounced at smaller scales in reduced dimensions.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Toward a First-Principles Framework for Predicting Collective Properties of Electrolytes

Conspectus Given the universal importance of electrolyte solutions, it is natural to expect that we have a nearly complete understanding of the fundamental properties of these solutions (e.g., the chemical potential) and that we can therefore explain, predict, and control the phenomena occurring in them. In fact, reality falls short of these expectations. But, recent advances in the simulation and modeling of electrolyte solutions indicate that it should soon be possible to make progress toward these goals. In this Account, we will discuss the use of first-principles interaction potentials based in quantum mechanics (QM) to enhance our understanding of electrolyte solutions. Specifically, we will focus on the use of quantum density functional theory (DFT) combined with molecular dynamics simulation (DFT-MD) as the foundation for our approach. The overarching concept is to understand and accurately reproduce the balance between local or short-ranged (SR) structural details and long-range (LR) correlations, allowing the prediction of the thermodynamics of both single ions in solution as well as the collective interactions characterized by activity/osmotic coefficients. In doing so, relevant collective motions and driving forces characterized by chemical potentials can be determined. Here, in this Account, we will make the case that understanding electrolyte solutions requires a faithful QM representation of the SR nature of the ion–ion, ion–water, and water–water interactions. However, the number of molecules that is required for collective behavior makes the direct application of high-level QM methods that contain the best SR physics untenable, making methods that balance accuracy and efficiency a practical goal. Alternatives such as continuum solvent models (CSMs) and empirically based classical molecular dynamics have been extensively employed to resolve this problem but without yet overcoming the fundamental issue of SR accuracy. We will demonstrate that accurately describing the SR interaction is imperative for predicting both intrinsic properties, namely, at infinite dilution, and collective properties of electrolyte solutions. DFT has played an important role in our understanding of condensed phase systems, e.g., bulk liquid water, the air–water interface, ions in bulk, and at the air–water interface. This approach holds huge promise to provide benchmark calculations of electrolyte solution properties that will allow for the development and improvement of more efficient methods, as well as an enhanced understanding of fundamental phenomena. However, the standard protocol using the generalized gradient approximation with van der Waals (vdW) correction requires improvement in order to achieve a high level of quantitative accuracy. Simply simulating with higher level DFT functionals may not be the best route considering the significant computational cost. Alternative methods of incorporating information from higher levels of QM should be explored; e.g., using force matching techniques on small clusters, where high level benchmark calculations are possible, to develop ideal correction terms to the DFT functional is a promising possibility. We argue that DFT with statistical mechanics is becoming an increasingly useful framework enabling the prediction of collective electrolyte properties.

Duignan, Timothy T.↗
DOE OSTI2023

Multi-level Force-dependent Allosteric Enhancement of αE-catenin Binding to F-actin by Vinculin

Classical cadherins are transmembrane proteins whose extracellular domains link neighboring cells, and whose intracellular domains connect to the actin cytoskeleton via β-catenin and α-catenin. The cadherin-catenin complex transmits forces that drive tissue morphogenesis and wound healing. In addition, tension-dependent changes in αE-catenin conformation enables it to recruit the actin-binding protein vinculin to cell-cell junctions, which contributes to junctional strengthening. How and whether multiple cadherin-complexes cooperate to reinforce cell-cell junctions in response to load remains poorly understood. Here, we used single-molecule optical trap measurements to examine how multiple cadherin-catenin complexes interact with F-actin under load, and how this interaction is influenced by the presence of vinculin. We show that force oriented toward the (-) end of the actin filament results in mean lifetimes 3-fold longer than when force was applied towards the barbed (+) end. We also measured force-dependent actin binding by a quaternary complex comprising the cadherin-catenin complex and the vinculin head region, which cannot itself bind actin. Binding lifetimes of this quaternary complex increased as additional complexes bound F-actin, but only when load was oriented toward the (-) end. In contrast, the cadherin-catenin complex alone did not show this form of cooperativity. Furthermore, these findings reveal multi-level, force-dependent regulation that enhances the strength of the association of multiple cadherin/catenin complexes with F-actin, conferring positive feedback that may strengthen the junction and polarize F-actin to facilitate the emergence of higher-order cytoskeletal organization.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Fluctuations and Pinning for Individually Manipulated Skyrmions

We numerically examine the dynamics of individually dragged skyrmions interacting simultaneously with an array of other skyrmions and quenched disorder. For drives just above depinning, we observe a broadband noise signal with a 1/ f characteristic, while at higher drives, narrowband or white noise appears. Even in the absence of quenched disorder, the threshold force that must be applied to translate the driven skyrmion is finite due to elastic interactions with other skyrmions. The depinning threshold increases as the strength of the quenched disorder is raised. Above the depinning force, the skyrmion moves faster in the presence of quenched disorder than in a disorder-free system since the pinning sites prevent other skyrmions from being dragged along with the driven skyrmion. For strong pinning, we find a stick-slip motion of the driven skyrmion which produces a telegraph noise signature. The depinning threshold increases monotonically with skyrmion density in the absence of quenched disorder, but when pinning is present, the depinning threshold changes nonmonotonically with skyrmion density, and there are reentrant pinned phases due to a competition between pinning induced by the quenched disorder and that produced by the elastic interactions of the skyrmion lattice.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Impact of three-body forces on elastic nucleon-nucleus scattering observables

In a previous series of papers we investigated the domain of applicability of chiral potentials to the construction of a microscopic optical potential (OP) for elastic nucleon-nucleus scattering. The OP was derived at the first order of the spectator expansion of the Watson multiple scattering theory and its final expression was a folding integral between the nucleon-nucleon ( N N ) t matrix and the nuclear density of the target. In the calculations N N and three-nucleon ( 3 N ) chiral interactions were used for the target density and only the N N interaction for the N N t matrix. The purpose of this study is to achieve another step towards the calculation of a more consistent OP introducing the 3 N force also in the dynamic part of the OP. The full treatment of the 3 N interaction is beyond our present capabilities. Thus, here it is approximated with a density dependent N N interaction obtained after the averaging over the Fermi sphere. In practice, in our model the 3 N force acts as a medium correction of the bare N N interaction used to calculate the t matrix. Even if the 3 N force is treated in an approximate way, this method naturally extends our previous model of the OP and allows a direct comparison of our present and previous results. We consider as case studies the elastic scattering of nucleons off C 12 and O 16 . We present results for the differential cross section and the spin observables for different values of the projectile energy. From the comparison with the experimental data and with the results of our previous model we assess the importance of the 3 N interaction in the dynamic part of the OP. Conclusions: Our analysis indicates that the contribution of the 3 N force in the t matrix is small for the differential cross section and it is sizable for the spin observables, in particular, for the analyzing power. We find that the two-pion exchange term is the major contributor to the 3 N force. A chiral expansion order-by-order analysis of the scattering observables confirms the convergence of our results at the next-to-next-to-next-to-leading-order, as already established in our previous work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

System and method for voltage-based real-time micro-object position control with the aid of a digital computer

The system and method described allow for real-time control over positioning of a micro-object. A movement of at least one micro-object suspended in a medium can be induced by a generation of one or more forces by electrodes proximate to the micro-object. Prior to inducing the movement, a simulation is used to develop a model describing a parameter of an interaction between each of the electrodes and the micro-object. A function describing the forces generated by an electrode and an extent of the movement induced due to the forces is generated using the model. The function is used to design closed loop policy control scheme for moving the micro-object towards a desired position. The position of the micro-object is tracked and taken into account when generating voltage patterns in the scheme.

Matei, Ion↗
DOE OSTI2026

Substantial Diel Changes of Cloud Adjustments to Aerosols in Ship‐Tracks

Clouds adjust their albedo, amount and liquid water path to anthropogenic aerosols, and produce a climate forcing. Such cloud adjustments are complex functions of both environmental conditions and interaction time. However, direct observations of temporal changes in cloud adjustments are limited. Here, we develop a method to show observational evidence of substantial diel changes in cloud adjustments within ship-tracks. Polluted clouds experiencing extended nighttime exposure exhibit higher increases in cloud fraction while liquid-water path adjustments follow a non-monotonic evolution. The higher cloud fraction increase persists through daytime. Such diel variations likely stem from diel changes in solar heating and background cloud precipitation, and are important because nighttime changes persist and affect solar radiation during the day. Accounting for diel variations produces a stronger aerosol forcing, by −0.1 ∼−0.4 W⁢m −2 . Our approach reveals how the temporal evolution of cloud adjustments can have important implications for aerosol-cloud interactions.

Yuan, Tianle [NASA Goddard Space Flight Center (GS↗
DOE OSTI2020

Multi-Scale Simulation of Wind Farm Performance during a Frontal Passage

Predicting the response of wind farms to changing flow conditions is necessary for optimal design and operation. Here in this work, simulation and analysis of a frontal passage through a utility scale wind farm is achieved for the first time using a seamless multi-scale modeling approach. A generalized actuator disk (GAD) wind turbine model is used to represent turbine–flow interaction, and results are compared to novel radar observations during the frontal passage. The Weather Research and Forecasting (WRF) model is employed with a nested grid setup that allows for coupling between multi-scale atmospheric conditions and turbine response. Starting with mesoscale forcing, the atmosphere is dynamically downscaled to the region of interest, where the interaction between turbulent flows and individual wind turbines is simulated with 10 m grid spacing. Several improvements are made to the GAD model to mimic realistic turbine operation, including a yawing capability and a power output calculation. Ultimately, the model is able to capture both the dynamics of the frontal passage and the turbine response; predictions show good agreement with observed background velocity, turbine wake structure, and power output after accounting for a phase shift in the mesoscale forcing. This study demonstrates the utility of the WRF-GAD model framework for simulating wind farm performance under complex atmospheric conditions.

17 WIND ENERGY↗
DOE OSTI2020

Finite-size effects on the self-force

Electromagnetic and linear gravitational radiation do not solely propagate on the null cone in 3 + 1 dimensions in curved spacetimes, contrary to their well-known behavior in flat spacetime. Their additional propagation inside the null cone is known as the tail effect. A compact body will produce a signal whose tail will interact with its future worldline, thus producing a tail-induced self-force. In this work, we present new results for the tail-induced scalar, electromagnetic, and gravitational self-force for a test mass in orbit around a central mass, including effects from the internal structure of that body.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1999

Effective Forces Between Colloidal Particles

Colloidal suspensions have proven to be excellent model systems for the study of condensed matter and its phase behavior. Many of the properties of colloidal suspensions can be investigated with a systematic variation of the characteristics of the systems and, in addition, the energy, length and time scales associated with them allow for experimental probing of otherwise inaccessible regimes. The latter property also makes colloidal systems vulnerable to external influences such as gravity. Experiments performed in micro-ravity by Chaikin and Russell have been invaluable in extracting the true behavior of the systems without an external field. Weitz and Pusey intend to use mixtures of colloidal particles with additives such as polymers to induce aggregation and form weak, tenuous, highly disordered fractal structures that would be stable in the absence of gravitational forces. When dispersed in a polarizable medium, colloidal particles can ionize, emitting counterions into the solution. The standard interaction potential in these charged colloidal suspensions was first obtained by Derjaguin, Landau, Verwey and Overbeek. The DLVO potential is obtained in the mean-field linearized Poisson-Boltzmann approximation and thus has limited applicability. For more precise calculations, we have used ab initio density functional theory. In our model, colloidal particles are charged hard spheres, the counterions are described by a continuum density field and the solvent is treated as a homogeneous medium with a specified dielectric constant. We calculate the effective forces between charged colloidal particles by integrating over the solvent and counterion degrees of freedom, taking into account the direct interactions between the particles as well as particle-counterion, counterion-counterion Coulomb, counterion entropic and correlation contributions. We obtain the effective interaction potential between charged colloidal particles in different configurations. We evaluate two- and three-body forces in the bulk as well as study the influence of soft walls. We qualitatively explain the effects of the walls on the forces and demonstrate that many-body effects are negligible in our system. With adjustments in the parameters, the DLVO pair-potential can describe the results quantitatively. Besides electrostatic interactions, entropic depletion effects that arise from (hard-core) exclusion play an important role in determining the behavior of multi-component colloidal suspensions. A standard theory for depletion forces is due to Asakura and Oosawa and is based on the ideal gas approximation. To go beyond this approximation, we have studied entropic forces in molecular dynamics simulations of systems of hard spheres (the effects of the solvent have been ignored). The effective depletion forces for these systems can be found either from equilibrium distribution functions or from direct momentum transfer calculations. Our results obtained by either method show qualitative differences from the Asakura-Oosawa forces, indicating a longer range, higher value at contact and most importantly a more complicated structure, comprising of several maxima and minima. Our calculations include the determination of effective forces between two spheres, a hard sphere and a wall, and the behavior of a hard sphere near a step-edge and a corner. We also demonstrate that such entropic forces do not necessarily satisfy pairwise additivity.

Tehver, Riina↗