Complexation Scenario: Mechanisms & Models
- Complexation scenario is a set of processes in which distinct chemical species form bound or partially bound states driven by coupled electrostatic, entropic, energetic, interfacial, and conformational effects.
- Mechanisms include counterion release, sliding-rod pathways, and multivalent ion screening, which dictate whether binding is entropy- or enthalpy-driven while influencing intermediate states.
- Structural outcomes range from finite complexes and coacervates to interfacial membranes and nanoparticle aggregates, with formulation variables and mixing order tuning the final morphology.
Searching arXiv for the cited complexation papers to ground the article in current records. arxiv_search({"query":"(Singh et al., 2019) complexation charged polypeptides (Whittaker et al., 2020) ion complexation waves (Baubigny et al., 2017) interfacial H-bond polymer complexation (1410.06547) PDADMAC PANa complexation (Baulin et al., 2012) interpolyelectrolyte complexes (Xu et al., 2016) potential of mean force transient states in polyelectrolyte complexation (Mitra et al., 2021) minimal theory (Vitorazi et al., 2024) mixing order asymmetry complexation precipitation (Buyukdagli et al., 2019) like-charge polymer-membrane complexation multivalent cations (Buyukdagli et al., 2019) DNA anionic membranes weak intermediate electrostatic coupling", "max_results": 10}) Complexation scenario denotes a class of association processes in which distinct chemical species, macromolecules, ions, interfaces, or colloids form bound or partially bound states whose stability, pathway, and morphology are controlled by coupled electrostatic, entropic, energetic, interfacial, and conformational effects. Across the literature considered here, the term encompasses oppositely charged polypeptides, polyelectrolytes, polymer–surfactant aggregates, interfacial hydrogen-bonded membranes, like-charge adsorption mediated by multivalent ions, transient dimer-driven mesoscale clustering, molten-salt impurity association, and surface complexation on porous sorbents. A unifying feature is that complexation is rarely a single elementary event: it is often pathway-dependent, can switch between entropy- and enthalpy-dominated regimes, and may terminate in finite complexes, coacervates, precipitates, adsorbed interfacial states, or metastable dense phases rather than a unique equilibrium aggregate (Singh et al., 2019, 1410.06547, Baulin et al., 2012, Davydov et al., 9 Sep 2025).
1. Thermodynamic driving forces
The thermodynamic character of complexation depends strongly on charge density, ion content, and the existence of bound or redistributed small ions. In molecular dynamics simulations of poly(lysine) and poly(glutamate), association of fully charged peptides is thermodynamically favorable and predominantly entropy-driven when neutralizing counterions are present or when excess salt is added, whereas removing all small ions makes complexation energetic rather than entropic (Singh et al., 2019). The same study reports that lowering the peptide charge density also shifts the driving force toward energetics, even though the partially charged peptides remain above the nominal Manning threshold (Singh et al., 2019). This suggests that small-ion content and charge density are control parameters for whether complexation is entropy- or enthalpy-driven.
A related but more explicitly analytical result appears in the minimal two-chain theory for symmetric oppositely charged polyelectrolytes. There, the free energy of complexation is decomposed into the entropy of condensed counterions, the entropy of free counterions and salt ions, Debye–Hückel ionic fluctuations, the Coulomb energy of bound ion pairs, conformational free energy, and screened Coulomb attraction between dangling charged segments (Mitra et al., 2021). That theory concludes that the entropy of free counterions and salt ions and the Coulomb enthalpy of bound ion-pairs dictate the equilibrium of PE complexation, with three regimes: enthalpy-driven complexation at low Coulomb strengths, entropy-driven but enthalpy-opposed complexation at moderate Coulomb strengths, and complexation prohibited by enthalpy loss at higher Coulomb strengths (Mitra et al., 2021). This suggests that electrostatics can either favor or suppress association, depending on whether it stabilizes the complexed state or over-stabilizes the separated counterion-condensed state.
The PMF study of oppositely charged coarse-grained polyelectrolytes reaches a complementary conclusion. Above the counterion-condensation threshold, the effective attraction is well described by the sum of screened rod–rod electrostatics with renormalized charges and an explicit entropy contribution from condensed-counterion release (Xu et al., 2016). In that framework, complexation of highly charged chains is driven primarily by counterion release, while direct electrostatics alone describe the weaker-charge regime less completely (Xu et al., 2016). The counterion theme recurs in calorimetric work on PDADMAC/PANa and related systems, where small positive enthalpies and large positive entropies are reported for both the primary and secondary association steps, consistent with a counterion-release scenario (1410.06547, Vitorazi et al., 2024).
Not all complexation scenarios are dominated by counterion entropy. In Ni–Cr association in eutectic FLiBe, the key result is a weak but thermodynamically meaningful binding free energy between dissolved and of approximately , interpreted as a fluorine-mediated second-shell cation association rather than a direct metal–metal bond (Attarian et al., 2024). In concentrated aqueous $\ce{FeCl2}$ and $\ce{MgCl2}$, by contrast, the main inference is that viscosity is highly sensitive to speciation, and that long-lived monochloro and dichloro complexes reduce viscosity by reducing the net amount of encumbered water molecules in contact with the ions (Goswami et al., 5 Mar 2026). In these cases, thermodynamics and transport are linked through speciation rather than through canonical counterion release.
2. Pathways, intermediate states, and kinetic structure
Complexation often proceeds through distinct intermediate states rather than direct collapse into the final morphology. In the PMF analysis of oppositely charged polyelectrolytes, the onset and intermediate stages are characterized by a sliding-rod-like pathway: the chains first stretch out at a critical distance close to their contour length, then “shake hand” and slide along each other in a parallel fashion before eventually folding into a neutral complex (Xu et al., 2016). The mean force is discontinuous at the onset of complex formation, reflecting coexistence between isolated-coil and interconnected-rod states separated by a free-energy barrier due to unfavorable stretching entropy (Xu et al., 2016). This is a mechanistic example in which complexation requires conformational reorganization before the main free-energy gain can be realized.
A two-step sequence also appears in isothermal titration calorimetry of PDADMAC/PANa. The first process is the formation of highly charged primary polyelectrolyte complexes of size about $100$–, and the second is the transition toward a coacervate phase made of rich and poor polymer droplets (1410.06547). The same two-step logic extends to nanoparticle–polymer systems, where titration first forms finite charged PECs and then induces phase separation or microphase separation depending on the chemistry (Vitorazi et al., 2024). In PDADMAC/CeO@PAA, the second step yields micron-scale precipitates, whereas in PTEA-b-PAm/-FeO0@PAA it yields finite neutral aggregates stabilized by the neutral block (Vitorazi et al., 2024).
Mixing order can alter these pathways. In PDADMAC/PANa, the charge of the primary PECs, the transition point to coacervation, and the sign of the secondary enthalpy all depend on whether PDADMAC is added into PANa or PANa into PDADMAC (1410.06547). A broader comparison across polymer/polymer and polymer/nanoparticle systems reaches the same conclusion: cationic species added to anionic ones produce exothermic secondary profiles, whereas the reverse order produces endothermic ones (Vitorazi et al., 2024). This suggests that primary complexes retain structural memory of the order in which local charge imbalance was established.
Kinetic structure is even more central in the dimer-based complexation scenario for mesoscopic clusters. In that framework, transient dimerization allows a steady-state ensemble of finite-sized droplets of a metastable dense liquid, with the characteristic droplet radius estimated as 1 in the earlier phenomenology that the later thermodynamic model seeks to justify (Davydov et al., 9 Sep 2025). The explicit dimer model then shows that the cluster phase is predicted to be a metastable monomer-rich dense liquid, whereas the macroscopically phase-separated dense liquid should be dimer-rich (Davydov et al., 9 Sep 2025). This suggests that long-lived finite complexes can arise from coupling between reaction kinetics and hidden metastable phases rather than from equilibrium finite aggregation.
In clay suspensions, time-dependent restructuring is expressed as waves of ion complexation and repeated exfoliation/restacking events. The paper does not provide a predictive continuum wave equation, but combines cryoET, cryoEM, and XPCS to argue that curvature-biased ion binding, elasticity, and hydration feed back on one another, producing abrupt and repetitive mesoscale rearrangements whose timescale depends strongly on counterion identity (Whittaker et al., 2020). Here the pathway itself is collective and non-equilibrium.
3. Structural motifs and morphologies
The structural outcome of complexation is highly system-specific. In the mean-field theory of spherical interpolyelectrolyte complexes, oppositely charged polymers condense into a dense spherical core permeable to micro-ions, while neutral blocks or loop-and-tail coronas arrest further growth and stabilize finite aggregates (Baulin et al., 2012). The equilibrium aggregation numbers 2 are selected near the electroneutrality line by balancing Poisson–Boltzmann electrostatics, corona entropy, and translational entropy (Baulin et al., 2012). This provides a canonical finite-complex scenario distinct from macroscopic phase separation.
In the thermoresponsive PECop/DTAB system, the complexes are spherical, water-rich, and often surprisingly monodisperse despite the broad polydispersity of the parent copolymer (Ritacco et al., 4 Aug 2025). At 3 and 4, static light scattering gives 5, 6, and Porod exponent 7, while 8, all interpreted as spherical particles (Ritacco et al., 4 Aug 2025). At 9, the same analysis gives 0 and 1, described as compatible with ellipsoidal particles with fractal surfaces (Ritacco et al., 4 Aug 2025). This suggests concentration-dependent structural reorganization within a single polymer–surfactant chemistry.
Interfacial complexation generates yet another morphology. At oil/water interfaces, protonated PMAA in water and PPO in oil form a membrane localized at the interface by hydrogen-bond complexation, without layer-by-layer deposition, templates, or covalent crosslinking (Baubigny et al., 2017). In planar geometry the membrane thickens over time according to 2, reaching about 3 after 100 days, with an effective diffusion coefficient of order 4, indicating diffusion-limited growth through the already formed membrane (Baubigny et al., 2017). The same chemistry yields centimeter-scale capsules by dripping, polydisperse microcapsules of 5 by rotor-stator emulsification, and monodisperse micron-scale capsules by microfluidics (Baubigny et al., 2017).
In thermo-enhanced surfactant–polymer interfacial complexation, a germinal interfacial polymer layer forms first by electrostatic association between PLL-g-PNIPAM and Krytox at the water/fluorocarbon-oil interface, then above 6 the collapse of PNIPAM enriches and densifies that shell (Sixdenier et al., 2021). Surface excess almost doubles on heating, and for droplets smaller than 7 essentially all polymer can be segregated into the interfacial layer (Sixdenier et al., 2021). This yields a persistent, particle-capturing mixed shell rather than a bulk aggregate.
At the interface between DNA and cationic gold nanoparticles, morphology is controlled by particle curvature and loading stoichiometry. A single 30 bp dsDNA wraps more tightly around smaller cationic AuNPs because of higher curvature, but when two DNAs bind the same small particle, electrostatic repulsion between exposed DNA ends prevents efficient wrapping and induces twisting from the original orientation (Monisha et al., 8 Sep 2025). Larger particles bind more strongly and reduce this twisting frustration (Monisha et al., 8 Sep 2025). This suggests that the structure of a complex can invert the naive expectation that stronger curvature always improves loading.
4. Ions, screening, and interfacial mediation
A recurrent mechanism across complexation scenarios is the action of interfacial or local ionic atmospheres rather than direct bare electrostatics alone. In the charged-polypeptide study, the favorable entropy of complexation is attributed not to release from a static condensed state, but to redistribution of ion probability density away from isolated peptide neighborhoods when oppositely charged chains neutralize one another (Singh et al., 2019). The direct evidence is that radial distribution functions and coordination numbers change gradually with peptide separation, with no sudden jump at the thermodynamic crossover distance (Singh et al., 2019).
In like-charge adsorption of a rodlike anionic polymer onto a like-charged membrane in the presence of tri- and tetravalent cations, the polymer grand potential is lowered because multivalent ions condense strongly at the membrane and maximize local screening ability, producing attraction even though the direct polymer–membrane coupling is repulsive (Buyukdagli et al., 2019). The theory separates a direct polymer–membrane term, a polymer self-energy term, and explicit polymer–counterion screening terms, showing that interfacial counterion excess acts as the adhesive force (Buyukdagli et al., 2019). Higher counterion valency lowers the critical concentration for adsorption, whereas added monovalent salt weakens interfacial multivalent condensation and causes desorption (Buyukdagli et al., 2019).
A related but distinct one-loop analysis for DNA near anionic membranes identifies two correlation-driven mechanisms. At intermediate membrane charge, an attractive salt-induced image force generated by enhanced screening near the membrane produces short-ranged like-charge adsorption and an orientational transition from parallel to perpendicular (Buyukdagli et al., 2019). At higher membrane charge, surface-induced ionic correlations can generate charge inversion, producing a separate mechanism for adsorption over an extended distance interval (Buyukdagli et al., 2019). This suggests that “image” forces in electrolyte complexation need not be dielectric in origin; they can arise from spatially varying screening alone.
In clay suspensions, curvature delocalizes electrostatics and makes convex and concave sides of the same layer inequivalent. The resulting asymmetry in inner-sphere and outer-sphere complexation is interpreted as a “complexation dipole,” not necessarily aligned normal to the layer, whose collective action produces an emergent long-range force and drives exfoliation and restacking (Whittaker et al., 2020). The measured layer thickness 8 decreases from 9 Å to $\ce{FeCl2}$0 Å as LiCl concentration rises from $\ce{FeCl2}$1 M to $\ce{FeCl2}$2 M, interpreted as a shift from more outer-sphere to more inner-sphere Li$\ce{FeCl2}$3 complexation (Whittaker et al., 2020). Curvature then breaks the symmetry of those equilibria.
In surface complexation of Mn(II) by biochar, the ionic mechanism is a sequence of cation release, pH increase, deprotonation of surface groups, and then inner- or outer-sphere binding (Ngambia et al., 23 Mar 2026). Molecular simulations show that neutral surfaces bind Mn weakly, whereas deprotonated low-temperature biochar models exhibit inner-sphere uptake of $\ce{FeCl2}$4 and $\ce{FeCl2}$5 for W400-DP and S400-DP, respectively, plus additional outer-sphere association (Ngambia et al., 23 Mar 2026). This suggests that interfacial charge generation by deprotonation is the chemical prerequisite for true surface-complexation-led sequestration.
5. Charge density, composition, and formulation variables
Charge density is one of the most consistent organizing variables across the literature. In realistic peptide systems, reducing each 11-mer to 6 charges per chain shifts the driving force from entropy to energetics under practical salt-free conditions, even though the estimated Manning-like parameters remain above unity (Singh et al., 2019). In the minimal theory of symmetric PE complexation, increasing the effective Coulomb strength $\ce{FeCl2}$6 first changes complexation from enthalpy-driven to entropy-driven and then suppresses complexation entirely (Mitra et al., 2021). In both cases, “more charge” does not simply mean “more binding”; it changes which free-energy term dominates.
Surfactant concentration can play an analogous role. In PECop/DTAB, low $\ce{FeCl2}$7 preserves thermoresponsive collapse, intermediate $\ce{FeCl2}$8 around $\ce{FeCl2}$9 maximizes compaction and eliminates thermal responsiveness, and high $\ce{MgCl2}$0 inverts the thermal response so that complexes grow rather than shrink above LCST (Ritacco et al., 4 Aug 2025). The binding isotherm is described as two consecutive and superimposed sigmoidal curves, with Region 1 associated mainly with backbone charge pairing and Region 2 with stronger hydrophobic contribution and probable structural rearrangement (Ritacco et al., 4 Aug 2025). This suggests that composition variables can change both the amount and the mechanism of complexation.
Mixing order is a formulation variable with similarly strong consequences. In PDADMAC/PANa, the primary complexes formed by PDADMAC into PANa are negatively charged, whereas those formed by PANa into PDADMAC are positively charged, and the secondary coacervation step is exothermic in the former case but endothermic in the latter (1410.06547). The 2024 nanoparticle/polymer calorimetry generalizes this asymmetry across four systems and identifies it as a consistent feature of the secondary process (Vitorazi et al., 2024). A plausible implication is that complexation pathways should be treated as state variables in formulation studies, not as experimental details.
Particle or substrate size can also determine which stoichiometries are feasible. In DNA-functionalized gold nanoparticles, a single DNA covers 67.50% of the circumference of the 2 nm particle, 52.85% of the 3 nm particle, and 41.99% of the 4 nm particle, helping explain why two DNAs twist and wrap incompletely on the smaller particles but remain roughly vertical on the 4 nm particle (Monisha et al., 8 Sep 2025). In the same paper, surface-area-based loading-capacity estimates rise from $\ce{MgCl2}$1 to $\ce{MgCl2}$2 DNAs as particle size increases from 2 to 4 nm, though the authors note that geometry alone overestimates practical loading on the smaller particles (Monisha et al., 8 Sep 2025).
6. Misconceptions, limits, and broader significance
A common misconception is that complexation can always be reduced to direct Coulomb attraction between opposite charges. Several papers explicitly reject that simplification. Realistic peptide association is not explained by bare peptide–peptide contact alone; the presence or absence of small ions changes the sign of the driving mechanism (Singh et al., 2019). Like-charge polymer–membrane adsorption is impossible at mean-field level and instead requires interfacial ionic correlations (Buyukdagli et al., 2019, Buyukdagli et al., 2019). In concentrated $\ce{MgCl2}$3 and $\ce{MgCl2}$4, the key indicator of speciation is not a direct structural count from unbiased simulation but the viscosity response required to reconcile experiment and simulation (Goswami et al., 5 Mar 2026).
A second misconception is that condensed or complexed ions must be static. The charged-polypeptide work argues against a literal picture of immobilized counterions released into translational freedom; all counterions retain nonzero self-diffusion coefficients and short lifetimes near the chains (Singh et al., 2019). The relevant entropy gain arises from removing positional bias rather than restoring motion from zero (Singh et al., 2019). The same caution appears in the FLiBe impurity study, where Ni–Cr association is explicitly not presented as a long-lived molecular adduct but as weak, transient short-range order with a negative PMF minimum (Attarian et al., 2024).
Several limitations recur across the corpus. Short oligomers are not long polymers (Singh et al., 2019). Mean-field theories neglect internal core structure, image chemistry, or strong-ion correlations outside their formal regime (Baulin et al., 2012, Mitra et al., 2021, Buyukdagli et al., 2019). Binding analyses may assume temperature-independent enthalpy and entropy over a finite interval (Singh et al., 2019). Mobility-based binding isotherms can oversimplify the notion of a binding site in graft copolymers with hydrophobic side chains (Ritacco et al., 4 Aug 2025). CryoET- and XPCS-based “complexation waves” remain mechanistic interpretations rather than a closed dynamical field theory (Whittaker et al., 2020). These caveats do not negate the underlying scenarios, but they delimit which claims are firmly established.
The broader significance of complexation scenarios lies in their role as precursors or regulators of larger-scale phenomena: complex coacervation (Singh et al., 2019, 1410.06547), finite PIC micelles (Baulin et al., 2012), responsive drug-delivery vectors (Ritacco et al., 4 Aug 2025, Sixdenier et al., 2021, Monisha et al., 8 Sep 2025), corrosion chemistry in molten salts (Attarian et al., 2024), concentrated-electrolyte transport (Goswami et al., 5 Mar 2026), clay swelling and rheology (Whittaker et al., 2020), and contaminant sequestration on engineered sorbents (Ngambia et al., 23 Mar 2026). Across these systems, the central encyclopedic conclusion is that complexation is best understood as a coupled free-energy problem in which local structure, ion distributions, pathway history, and molecular architecture determine not only whether association occurs, but which associated state is selected and whether it is stable, metastable, finite, interfacial, or phase-separating.