Asymmetric Kedem–Katchalsky Conditions

• Asymmetric Kedem–Katchalsky conditions are interface laws that govern transport across membranes by incorporating unequal side contributions via parameters like permeability and partition coefficients.
• They are implemented in various settings including two-layer diffusion, interface logistic systems, and thin-membrane limits, capturing direction-dependent flux and concentration discontinuities.
• Applications span drug-eluting stents, ecological models, and reaction–diffusion systems, where the asymmetry leads to non-self-adjoint operators and unique bifurcation phenomena.

Asymmetric Kedem–Katchalsky boundary conditions are interface laws for transport across a membrane or sharp internal boundary in which the two sides do not enter on an equivalent footing. In their linear diffusion form, they couple flux continuity to a membrane law such as $J=P(c_1-\sigma c_2)$, so that the driving jump is not the symmetric difference $c_1-c_2$ but a weighted difference determined by the partition coefficient $\sigma$; in other settings, asymmetry is expressed by distinct interface coefficients on the two sides, as in $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$ and $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$, or by nonlinear jumps such as $\llbracket \Pi(u)\rrbracket$ in thin-membrane limits. Recent work places these conditions at the intersection of membrane transport, stochastic simulation, elliptic and parabolic transmission problems, spectral theory, and singular perturbation analysis (Farago et al., 2020, Álvarez-Caudevilla et al., 24 Jul 2025, Ciavolella et al., 2021).

1. Canonical formulations

The most common linear setting is a multilayer diffusion problem in which each layer satisfies

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$

while the interface enforces flux continuity together with a Kedem–Katchalsky law,

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$

Here $P$ is the permeability coefficient and $\sigma$ is the partition coefficient. This pair of relations is the complete interface law in the two-layer diffusion model developed for Langevin simulation, and it is asymmetric whenever $c_1-c_2$0; in practice it is also asymmetric when $c_1-c_2$1 (Farago et al., 2020).

A second, explicitly directional formulation appears in interface logistic systems. There the two populations live in adjacent subdomains and interact only through the common interface, where the membrane law is written

$c_1-c_2$2

The coefficients $c_1-c_2$3 and $c_1-c_2$4 are generally different, so the law is asymmetric in a stronger sense than the partition-weighted diffusion law: the two sides are assigned distinct transfer coefficients (Álvarez-Caudevilla et al., 2024).

A third formulation arises in thin-membrane limits for porous-medium or Darcy systems. After collapse of a membrane of thickness $c_1-c_2$5 to a zero-thickness interface, the effective transmission condition becomes

$c_1-c_2$6

with $c_1-c_2$7. For the porous-medium pressure law $c_1-c_2$8,

$c_1-c_2$9

This law is compatible with nonlinear generalized Kedem–Katchalsky conditions, although it does not introduce separate one-sided membrane permeabilities (Ciavolella et al., 2021).

Setting	Interface law	Asymmetry mechanism
Two-layer diffusion	$\sigma$0	$\sigma$1, $\sigma$2
Interface logistic system	$\sigma$3, $\sigma$4	$\sigma$5
Thin-membrane porous-medium limit	$\sigma$6 common normal flux	nonlinear jump, heterogeneous one-sided traces

2. Mechanisms of asymmetry

The partition coefficient is the simplest asymmetry mechanism. In the diffusion law $\sigma$7, the interface flux is driven by $\sigma$8 rather than by a plain concentration difference. If $\sigma$9, the two sides are weighted differently; physically this corresponds to a jump in chemical potential across the interface. The same framework makes clear that finite permeability and partitioning are distinct effects: when $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$0, the law reduces to the partition condition $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$1, while for finite $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$2 and $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$3 one obtains the pure membrane law $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$4, which still permits a concentration jump because the membrane has finite resistance (Farago et al., 2020).

In Brownian-diffusion formulations, the same asymmetry is expressed through the relation

$\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$5

where $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$6 is the chemical-potential jump. The boundary condition

$\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$7

is therefore direction-dependent in effect whenever $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$8, even though the membrane permeability parameter $\partial_{\nu_1}u_1=\gamma_1(u_2-u_1)$9 itself is taken symmetric (Farago, 2020).

A stronger notion of asymmetry appears when the interface law itself uses different coefficients on the two sides. In the two-species stationary interface logistic system, adding the two interface relations yields

$\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$0

so continuity of total normal flux fails unless $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$1. In that setting the operator is non-self-adjoint, and the asymmetry is structural rather than merely a consequence of unequal bulk coefficients (Álvarez-Caudevilla et al., 24 Jul 2025).

Not all generalized Kedem–Katchalsky laws are asymmetric in this explicit directional sense. In the porous-medium thin-membrane limit, the interface constitutive law uses a single effective permeability $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$2, and the surviving asymmetry comes from unequal bulk mobilities $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$3 and from the one-sided traces at the limit interface. The paper is explicit that this is not a law with distinct directional membrane coefficients (Ciavolella et al., 2021).

A more mechanistic perspective comes from the kinetic-theory derivation of osmotic transport in porous media. There, $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$4 is interpreted as the fraction of solute molecules reflected at pore entrances and exits, and a further parameter $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$5 measures momentum transfer from rebounding solute molecules to nearby solvent molecules. The paper does not derive asymmetric left/right boundary conditions directly, but this suggests that genuinely directional Kedem–Katchalsky laws could be built by allowing $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$6 or $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$7 to depend on the side of approach (Cardoso et al., 2014).

3. Stochastic realization and computational implementation

A major development is the translation of asymmetric Kedem–Katchalsky interface laws into particle-based Langevin dynamics. In the two-layer diffusion model, particle motion is generated by the underdamped Langevin equation

$\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$8

with Einstein relation $\partial_{\nu_2}u_2=\gamma_2(u_1-u_2)$9, and numerically integrated by the GJF scheme. The two Kedem–Katchalsky parameters are represented separately. First, finite permeability is converted into a transmission probability

$\llbracket \Pi(u)\rrbracket$0

where

$\llbracket \Pi(u)\rrbracket$1

Second, the partition coefficient is represented as a chemical-potential step

$\llbracket \Pi(u)\rrbracket$2

approximated by a thin force layer together with the trajectory weight

$\llbracket \Pi(u)\rrbracket$3

When a trajectory attempts to cross the interface, a uniform random number determines reflection or transmission; after transmission, the step is recomputed with an effective friction

$\llbracket \Pi(u)\rrbracket$4

The method reproduces the continuum law

$\llbracket \Pi(u)\rrbracket$5

shows excellent agreement with the continuum diffusion problem, captures the concentration jump at the interface, and converges as the interface-layer thickness parameter $\llbracket \Pi(u)\rrbracket$6 is reduced. In the drug-eluting-stent example, the error at $\llbracket \Pi(u)\rrbracket$7 on the right side at $\llbracket \Pi(u)\rrbracket$8 decreases from $\llbracket \Pi(u)\rrbracket$9 for $$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$0 to $$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$1 for $$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$2; the timestep must remain in the ballistic regime,

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$3

so that particles resolve the interface layer (Farago et al., 2020).

A closely related Brownian-dynamics framework decomposes the general interface into three primitive discontinuities: a friction jump, a semi-permeable membrane, and a step chemical potential. In that formulation the membrane law is again

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$4

with transmission probability

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$5

while the partition step is realized by a narrow interfacial potential layer and Boltzmann reweighting. The general Kedem–Katchalsky boundary is then implemented by combining ballistic averaging across the diffusivity jump, probabilistic transmission/reflection for permeability, and interface-layer forcing for the chemical-potential jump (Farago, 2020).

4. Variational structure, operator theory, and nonlinear steady states

A useful baseline is provided by symmetric membrane problems. In a reaction–diffusion system on two subdomains separated by a membrane, the condition

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$6

is conservative and dissipative: flux continuity follows from $$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$7, and the weak formulation contains the nonnegative jump penalty

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$8

The associated bilinear form has the structure

$$\frac{\partial c_i}{\partial t}=D_i\frac{\partial^2 c_i}{\partial x^2},$$9

which underlies the global weak-solution theory for $$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$0 data. This is Kedem–Katchalsky in a conservative symmetric form, not yet an asymmetric law (Ciavolella et al., 2020).

Membrane spectral theory extends naturally to reaction–diffusion systems used for Turing analysis. For each species, the membrane conditions

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$1

lead to symmetric bilinear forms with positive jump terms on $$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$2, compact self-adjoint membrane operators, and permeability-dependent eigenvalue sequences. The resulting patterns may be discontinuous at the membrane because the law permits nonzero jumps in concentration for finite $$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$3 (Ciavolella, 2021).

The fully asymmetric case changes the analytical picture. In the stationary two-species/interface logistic system with

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$4

the interface term in the bilinear form becomes

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$5

which is symmetric only when $$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$6. The paper establishes a unique positive principal eigenvalue

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$7

an upper threshold

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$8

and the sharp existence criterion

$$-D_1\frac{\partial c_1}{\partial x}=-D_2\frac{\partial c_2}{\partial x}=P(c_1-\sigma c_2).$$9

It also proves that $P$0 is the unique bifurcation point from extinction and that, when $P$1, the positive steady state exhibits non-simultaneous blow-up: $P$2 uniformly on $P$3 while $P$4 remains bounded in $P$5 as $P$6 (Álvarez-Caudevilla et al., 24 Jul 2025).

5. Thin-membrane limits and singular-perturbation regimes

The most rigorous derivation of a generalized Kedem–Katchalsky interface law from a resolved membrane model is obtained for a porous-medium equation under Darcy’s law. The prelimit problem is posed on three subdomains, with a membrane of thickness $P$7 and mobility scaling

$P$8

Before passing to the limit, the model imposes continuity of density and continuity of normal flux across both membrane faces. As $P$9, the membrane collapses to a zero-thickness interface and the effective law becomes

$\sigma$0

where $\sigma$1. For $\sigma$2,

$\sigma$3

When $\sigma$4, one recovers the classical linear jump in $\sigma$5, so the result can be read as a nonlinear generalization of the classical Kedem–Katchalsky interface law (Ciavolella et al., 2021).

Singular-perturbation analysis in elliptic interface logistic problems shows how membrane asymmetry persists in effective limits. For the asymmetric interface law

$\sigma$6

the principal eigenvalue satisfies

$\sigma$7

while the nonlinear positive solution obeys

$\sigma$8

By contrast, as $\sigma$9,

$c_1-c_2$00

The small-diffusion regime is therefore dominated by local carrying capacities, whereas the large-diffusion regime produces a membrane-weighted homogenization law (Álvarez-Caudevilla et al., 2024).

6. Applications, misconceptions, and scope

The interface law is not merely formal. In the drug-eluting-stent model, the coating occupies $c_1-c_2$01, the arterial wall $c_1-c_2$02, and the membrane at $c_1-c_2$03 is modeled by

$c_1-c_2$04

with

$c_1-c_2$05

$c_1-c_2$06

Because $c_1-c_2$07, the interface energetically disfavors the right side and produces a substantial concentration jump at the membrane (Farago et al., 2020).

In ecological interface problems, asymmetric Kedem–Katchalsky conditions determine existence thresholds, bifurcation from extinction, and singular behavior near refuge zones. In particular, the directional coefficients $c_1-c_2$08 enter the principal eigenvalue, alter the effective weighted averages seen in large-diffusion limits, and permit non-self-adjoint coupling strong enough to support non-simultaneous blow-up (Álvarez-Caudevilla et al., 24 Jul 2025).

In reaction–diffusion pattern formation, membrane conditions with species-dependent permeabilities $c_1-c_2$09 and $c_1-c_2$10 modify the membrane Laplacian spectrum and admit discontinuous stationary patterns. This shows that Kedem–Katchalsky conditions should not be conflated with continuity transmission: for finite permeability, the solution itself may jump across the membrane, while the interface law ties that jump to the normal flux (Ciavolella, 2021).

A common misconception is that “asymmetric” necessarily means separate left-to-right and right-to-left permeability coefficients. The recent literature shows several inequivalent meanings. In the two-layer diffusion and Brownian formulations, a single permeability $c_1-c_2$11 is retained and asymmetry is encoded by $c_1-c_2$12 and by heterogeneous bulk transport coefficients. In the porous-medium thin-membrane limit, the generalized law is nonlinear and uses one effective membrane permeability, so asymmetry is implicit through traces and bulk mobilities. Only some interface logistic models are asymmetric in the strongest sense, with distinct side coefficients $c_1-c_2$13 and $c_1-c_2$14 (Farago, 2020, Ciavolella et al., 2021, Álvarez-Caudevilla et al., 2024).

Taken together, these formulations define a broad research program. Asymmetric Kedem–Katchalsky boundary conditions may describe weighted concentration jumps, side-dependent membrane exchange, or nonlinear effective transmission laws derived from thin layers. Their mathematical signatures include concentration discontinuities, non-self-adjoint interface operators, membrane-weighted homogenized limits, and interface spectra that differ fundamentally from those of standard continuity transmission problems.