An ion model is any simplified representation of how charged atoms or molecules behave, whether inside a crystal, dissolved in water, passing through a cell membrane, or trapped in a vacuum chamber for quantum computing. These models range from a century-old picture of point charges sitting on a lattice to modern molecular dynamics simulations tracking millions of interactions per nanosecond. What unites them is a shared challenge: capturing the physics of electrostatic attraction and repulsion accurately enough to predict real-world behavior, without making the math impossibly complex. The diversity of ion models across chemistry, physics, biology, and engineering reflects just how many problems hinge on getting the behavior of charged particles right.
Point Charges on a Lattice
The oldest and most intuitive ion model treats ions as point-like charges arranged on a regular lattice. Max Born proposed this picture in the early twentieth century to calculate the energy holding ionic crystals together. The idea is straightforward: each ion sits at a fixed position, surrounded preferentially by ions of the opposite charge. The sum of all those attractive and repulsive electrostatic interactions gives what is called the Madelung energy. In a typical ionic crystal like sodium chloride, the Madelung energy turns out to be remarkably close to the experimentally observed cohesive energy, which is the total energy required to pull the crystal apart into isolated ions.
1Solid State Physics. Cohesion of Ionic Solids in the Born ModelThe Born model was originally designed just for lattice energy calculations, but its success and simplicity led researchers to extend it to other crystal properties: elastic constants, thermal expansion, defect formation energies, and more. It works best for highly ionic compounds where the electron clouds of anions and cations barely overlap. The further a real material strays from that ideal picture, the less reliable the Born model becomes. But for something like table salt or calcium fluoride, it remains a surprisingly powerful starting point even today.
A practical difficulty with lattice energy is that you cannot measure it directly by experiment. Instead, scientists use a thermochemical cycle, combining measurable quantities like ionization energies and enthalpies of formation, to arrive at a value that serves as the experimental benchmark. Alternative computational methods have been developed to cross-check these benchmarks, and the agreement between thermochemical and computational values helps confirm the validity of point-charge models for simple ionic solids.
2Physica B: Condensed Matter. A novel lattice energy calculation technique for simple inorganic crystalsBeyond Point Charges
Treating ions as featureless point charges works fine until you need to explain what happens when two ions get very close. At short distances, their electron clouds overlap, and the quantum-mechanical requirement that no two electrons occupy the same state creates a strong repulsive force. The Born model handles this with a simple mathematical correction, typically an exponential or power-law term that rises steeply at close range. But the real shape of this repulsion depends on the geometry of the electron clouds, not just the distance between nuclei.
Modern force-field models capture this more realistically. One approach uses anisotropic, multipole-based repulsion terms, fitting them to high-level quantum mechanical calculations for pairs of molecules. A recent model of this kind derived parameters for 26 different atom classes covering most organic molecules, using data from over a hundred molecular pairs.
3The Journal of Chemical Physics. Classical Pauli repulsion: An anisotropic, atomic multipole modelThese refined repulsion models matter most in simulations of biological molecules and drug design, where the precise shape of how one molecule pushes another away can determine whether a drug binds to its target or slides off.
Ions Dissolved in Water
Modeling ions in solution brings a fundamentally different challenge. In a crystal, each ion has fixed neighbors. In a liquid, ions are surrounded by constantly shifting shells of water molecules and other ions, and the net effect of all those interactions determines everything from the solution’s electrical conductivity to whether a protein folds correctly.
The classical model for dilute electrolyte solutions comes from Peter Debye and Erich Hückel, who in the 1920s worked out how a cloud of oppositely charged ions screens the electric field around any given ion. Their theory predicts that electrostatic interactions decay exponentially with distance, with a characteristic screening length that shrinks as ion concentration increases. This makes intuitive sense: more ions means more effective shielding of any one ion’s charge.
The Debye-Hückel model works well for dilute salt solutions, but it breaks down at higher concentrations where ions are packed closely enough that their finite size, the structure of surrounding water, and correlations between ion positions all matter. Generalized versions of the theory attempt to account for these effects by treating ions and solvent molecules as objects with real volumes and shapes, incorporating how the solvent’s ability to screen charge changes with ionic strength, temperature, and pressure.
4ACS Publications. Generalized Debye–Hückel Theory for Ion Activities in Mixed-Salt and Mixed-Solvent Electrolyte SolutionsWhen Concentrated Electrolytes Break the Rules
One of the more surprising recent findings in electrolyte science is that the standard story about screening lengths breaks down dramatically in concentrated solutions. Classical theory predicts that as you add more salt, the screening length keeps getting shorter. But experiments measuring forces between charged surfaces across a range of electrolytes show the opposite: beyond a certain concentration, the screening length starts increasing again.
This counterintuitive behavior has been observed across aqueous salt solutions, ionic liquids diluted with organic solvents, and pure ionic liquids. When the screening lengths for all these different systems are scaled by the dielectric constant, they collapse onto a single curve, suggesting a universal phenomenon rather than a quirk of any particular ion pair.
5PubMed. The Electrostatic Screening Length in Concentrated Electrolytes Increases with ConcentrationThe practical significance is substantial: long-range electrostatic forces in concentrated electrolytes and ionic liquids are much stronger than classical models predict. This matters for batteries, capacitors, and biological systems where ion concentrations are high.
One proposed explanation involves the idea that in very concentrated solutions, most ions are bound up in neutral clusters or correlated domains, and only a small fraction exist as free charge carriers at any moment. These thermally activated charge fluctuations then act like the dilute ions in classical theory, but there are far fewer of them than the total ion concentration would suggest, so the effective screening length is longer.
6PubMed Central. Long-range electrostatic screening in ionic liquidsIon Models in Cell Biology
Every living cell maintains carefully controlled differences in ion concentrations across its membrane. Potassium is concentrated inside cells, sodium outside, and the resulting voltage difference across the membrane drives nerve impulses, muscle contractions, and nutrient transport. Modeling this requires treating the membrane as a barrier with selective permeability to different ions.
The foundational equation for this was published by David Goldman in 1943. It treats ions as moving independently through a membrane where the electric field is constant, and it produces a concise expression for the membrane current as a function of ion concentrations and voltage. This constant-field equation, later extended by Hodgkin and Katz, remains the theoretical pillar for interpreting how a cell’s membrane potential relates to its internal and external ionic composition.
7PubMed Central. The enduring legacy of the “constant-field equation” in membrane ion transportThe assumptions behind this equation are admittedly idealized: ions move independently, the electric field is uniform across the membrane, and there is no interaction between different ion species. Various generalizations have extended the model to handle mixtures of ions with different charges and the effects of electrogenic pumps that actively move ions against their concentration gradients.
8Mathematical Biosciences. Generalizations of the Goldman-Hodgkin-Katz equationTesting whether the Goldman equation actually holds in real tissue has produced mixed results. In gallbladder epithelium, researchers compared permeability ratios calculated from the Goldman equation against ratios measured directly using radioactive tracers. The equation turned out to be empirically valid for potassium chloride and rubidium chloride across a range of concentrations, even though none of its theoretical assumptions strictly held for the paracellular pathways being measured.
9Biochimica et Biophysica Acta (BBA) – Biomembranes. Validity of the Goldman-Hodgkin-Katz equation in paracellular ionic pathways of gallbladder epitheliumThis is a common pattern in ion modeling: simplified models give the right answer for reasons that are not entirely the ones the theory assumes.
How Ion Channels Pick Favorites
Ion channels are protein structures embedded in cell membranes that allow specific ions to pass through while blocking others. Potassium channels, for instance, can select potassium over sodium with ratios as high as 100 to 1000 to 1, despite sodium being a smaller ion.
10PubMed Central. From Selective Permeation to Physiology in Potassium ChannelsThe selectivity filter of a potassium channel contains a series of binding sites that coordinate potassium ions by mimicking the arrangement of water molecules that normally surround potassium in solution. When a potassium ion enters the filter, it trades its hydration shell for the oxygen atoms lining the channel, and the energy cost is approximately matched by the energy gain. Sodium, with its different size and hydration properties, faces a much less favorable trade.
Structural studies have added nuance to this picture. Small cations like sodium and lithium do bind within the potassium selectivity filter, but at a separate site from where potassium binds. The selectivity appears to result primarily from a large energy barrier that blocks sodium and lithium from entering the filter when potassium is present, rather than from a simple difference in how strongly each ion binds once inside.
11PubMed Central. Mechanism of potassium-channel selectivity revealed by Na(+) and Li(+) binding sites within the KcsA poreSelective conduction then proceeds through a “knock-on” mechanism: an entering ion destabilizes the next ion ahead of it in the single-file queue through the channel, pushing it forward and maintaining rapid throughput.
The Hofmeister Effect and Hydration Shells
Since 1888, scientists have known that different ions have dramatically different effects on the behavior of proteins, colloids, and other macromolecules in solution. This ordering, called the Hofmeister series, ranks ions by their ability to salt out proteins, stabilize or destabilize foams, and influence countless other aqueous processes. For over a century, the standard explanation was that some ions are “structure makers” that strengthen the hydrogen bond network of surrounding water, while others are “structure breakers” that disrupt it.
Recent experimental work has largely overturned this explanation. Time-resolved spectroscopy and thermodynamic measurements show that ions have surprisingly little effect on the bulk structure of water beyond their immediate hydration shell. Instead, the Hofmeister effect appears to arise from direct interactions between ions and macromolecular surfaces, along with interactions with water molecules in the first hydration shell of the macromolecule itself.
12PubMed Central. Interactions between macromolecules and ions: The Hofmeister seriesThis shift matters for ion models because it means you cannot simply parameterize an ion’s effect on a biological system through its bulk-water properties. You need to account for what happens at the interface between ion, water, and macromolecule, which is a much harder modeling problem.
Molecular Dynamics and the Polarization Problem
Molecular dynamics simulations track the motion of atoms and molecules by computing forces at each time step and updating positions accordingly. For ion-containing systems, the dominant force is electrostatic, so getting the charge model right is critical. The simplest approach assigns a fixed partial charge to each atom and uses Coulomb’s law. These nonpolarizable models are computationally cheap and often adequate for trends, but they miss something important: in reality, the electron clouds of atoms shift in response to nearby charges.
This omission can lead to serious quantitative errors. In simulations of a concentrated lithium battery electrolyte, for example, nonpolarizable models predicted ion diffusion coefficients more than 100 times slower than experimental measurements. The underlying problem is that nonpolarizable models overestimate the strength of electrostatic interactions at long range. Polarizable models partially fix this by allowing atomic charges to respond to their environment, producing a high-frequency dielectric constant around 2, while nonpolarizable models are stuck at 1.
13ACS Publications. Molecular Dynamics Simulations of Ionic Liquids and Electrolytes Using Polarizable Force FieldsA widely used workaround for nonpolarizable models is to scale ionic charges, typically multiplying them by about 0.7 to 0.8, which mimics some of the screening that polarization would provide. This trick improves dynamics and diffusion coefficients, but it distorts other properties like the dielectric response. The field has not settled on a single solution; the choice between polarizable and scaled-charge nonpolarizable models remains one of the most consequential decisions in any ionic simulation.
Modeling Ions in Energy Storage
Battery performance depends heavily on how ions move through electrolytes, electrode materials, and the thin films that form at electrode surfaces. In lithium-ion batteries, a layer called the solid electrolyte interphase forms on the anode during the first few charging cycles. This layer is crucial: it must allow lithium ions through while blocking electrons and solvent molecules. If it grows too thick or too resistive, the battery loses capacity.
Molecular dynamics simulations of a model solid electrolyte interphase compound, lithium ethylene dicarbonate, have identified three distinct transport regimes for lithium ions: an initial ballistic phase where ions move freely, a trapping phase where ions are temporarily stuck in local energy minima, and a longer-term diffusive phase. At moderate temperatures, the lithium diffusivity in this material is extremely low, on the order of 10⁻¹² cm²/s, and the resulting ionic conductivity is roughly 4.5 × 10⁻⁹ S/cm, consistent with both previous simulations using polarizable force fields and with experimental data.
14PubMed Central. Molecular Dynamics of Lithium Ion Transport in a Model Solid Electrolyte InterphaseMore detailed models of this interphase layer incorporate chemical reactions and multicomponent transport, tracking how different species form, migrate, and react within the film as the battery operates.
15Electrochimica Acta. Modeling detailed chemistry and transport for solid-electrolyte-interface (SEI) films in Li–ion batteriesGetting these models right has direct commercial implications: a better understanding of ion transport in the interphase could lead to batteries that charge faster, last longer, and are less prone to degradation.
Ions and Cloud Formation
Ion models also play a role far from the laboratory. In the atmosphere, cosmic rays and radioactive decay produce ions that can serve as seeds for the formation of tiny aerosol particles, a process called ion-induced nucleation. Whether these ions contribute meaningfully to cloud formation, and therefore to climate, has been debated for decades.
Models of ion-induced nucleation treat the clustering of sulfuric acid and water molecules around an atmospheric ion, calculating the free energy barrier that must be overcome for a stable particle to form. The presence of an ion significantly lowers this barrier by stabilizing the growing cluster through electrostatic attraction. One model found that the negative ion mechanism involving sulfuric acid and water is an efficient source of new particles in the middle and upper troposphere, where preexisting particle concentrations are low enough that new clusters can survive rather than being scavenged.
16Journal of Geophysical Research: Atmospheres. Atmospheric ion‐induced nucleation of sulfuric acid and waterNot all atmospheric ions are equally effective at promoting nucleation. Molecular dynamics combined with classical nucleation theory has ranked common atmospheric ions by their ability to enhance nucleation rates: sulfate is the strongest, followed by hydronium, ammonium, and nitrate. This ranking coincides with the order of their solvation free energies, suggesting that an ion’s affinity for surrounding water molecules is a rough predictor of its nucleation-promoting ability.
17The Journal of Physical Chemistry A. Free Energy Prediction of Ion-Induced Nucleation of Aqueous AerosolsThe atmospheric modeling community has also examined whether ion-induced nucleation can explain particle formation events observed near the Earth’s surface. The conclusion so far is mixed: the ion mechanism appears important at high altitudes but does not generally account for the nucleation bursts seen in the boundary layer, where other chemical pathways involving amines and organic vapors are more significant.
18Journal of Geophysical Research: Atmospheres. Model studies on ion‐induced nucleation in the atmosphereIons Under Extreme Pressure
Ion models for crystals face a severe test at high pressures, where the distances between ions shrink enough to dramatically alter the balance between attraction and repulsion. Calcium fluoride, a classic ionic crystal, undergoes a structural phase transition under pressure, shifting from the fluorite structure to the cotunnite structure. Researchers have observed this transition in real time using synchrotron X-ray diffraction during both rapid dynamic loading to about 28 GPa and slow static compression.
19PubMed. Direct Observations of a Dynamically Driven Phase Transition with in situ X-Ray Diffraction in a Simple Ionic CrystalThe ability of ion models to predict the pressure at which such transitions occur, and the structures that result, is one of the most demanding benchmarks in solid-state physics. The fact that a crystal as “simple” as calcium fluoride can adopt entirely different arrangements depending on pressure illustrates how sensitive the energy landscape is to the details of the ion-ion interaction model.
Trapped Ions as Quantum Hardware
In a very different setting, individual ions suspended in electromagnetic traps have become one of the leading platforms for quantum computing. A single trapped ion, cooled to near absolute zero, can encode quantum information in both its internal electronic states and its vibrational motion within the trap. The vibrational modes of a trapped ion behave as quantum harmonic oscillators, offering a theoretically infinite-dimensional space for encoding and processing information.
20Chinese Physics B. Quantum computation and simulation with vibrational modes of trapped ionsRecent work has explored using these vibrational modes not just as a communication bus between ion qubits, but as the primary computational resource. The long coherence times available in trapped-ion systems, combined with precise control over quantum operations, make this an attractive approach for certain types of problems, particularly those naturally suited to continuous-variable encoding rather than the discrete qubit model used in most quantum computing architectures.
21arXiv. Quantum computation over the vibrational modes of a single trapped ionThe “ion model” here is not a representation of an ion’s internal structure, but a physical ion used as a controllable quantum system, arguably the most literal possible ion model, since the model and the thing being modeled are the same object.