---
title: Pressure Threshold Model (PT)
url: https://www.emergentmind.com/topics/pressure-threshold-model-pt
type: topic
---

# Pressure Threshold Model (PT)

Pressure Threshold Model (PT) designates threshold-driven dynamics in which state transitions depend on an accumulated notion of pressure, but the term is used for two technically distinct model families. In clarinet acoustics, PT denotes a dynamic-threshold framework for predicting the onset of self-sustained oscillations in an iterated-map model under time-varying mouth pressure, extending the static oscillation threshold to regimes with bifurcation delay [1207.4636]. In social-network diffusion, PT denotes a graph process that extends the Linear Threshold (LT) model by increasing a newly activated node’s outgoing influence proportionally to the pressure it received at activation [2509.12822]. The common vocabulary of pressure and threshold therefore spans two separate mathematical settings: a nonlinear delayed map for reed–bore interaction, and a progressive synchronous diffusion process on weighted directed graphs.

## 1. Terminological scope and domain-specific meanings

In the clarinet literature, the central problem is the onset of oscillation when the musician’s mouth pressure is not constant but increases through time. Simple clarinet models based on iterated maps successfully estimate the threshold of oscillation as a function of a constant blowing pressure, yet when the blowing pressure gradually increases, oscillations appear at a much higher value than in the static case; this is the dynamic oscillation threshold associated with bifurcation delay [1407.3547]. In that usage, PT is effectively a pressure-threshold framework for onset prediction in a non-autonomous nonlinear dynamical system.

In network science, PT is explicitly introduced as the “Pressure Threshold model,” a diffusion model for influence propagation on social networks. It preserves the LT activation rule but adds a feedback mechanism: outgoing influence from a newly activated node is amplified according to the influence pressure that caused its activation [2509.12822]. Here PT is not an onset threshold in an acoustic oscillator but a threshold diffusion process with adaptive edge weights.

A common misconception is that PT denotes a single standardized model across disciplines. The available literature instead uses the same label for unrelated systems with different state variables, update operators, and analytical objectives. In the clarinet case, the state is the outgoing acoustic wave \(p_n^+\) and the threshold concerns a flip bifurcation; in the network case, the state is a node activation indicator \(s_v(t)\) and the threshold concerns cumulative in-neighbor influence.

## 2. Clarinet PT as an iterated-map onset model

The clarinet PT model is built on a generator–resonator decomposition. The exciter is the reed–mouthpiece system, modeled by a nonlinear characteristic \(F\) linking mouthpiece pressure \(p\) to volume flow \(u\). The resonator is a lossless cylindrical bore with one-delay sign-inverting reflection, so that in discrete time \(p_n^-=-r p_{n-1}^+\), with \(r=1\) in the lossless case treated in the paper [1407.3547].

The normalization is centered on the reed-closing pressure \(P_M=kH\), where \(k\) is reed stiffness per unit displacement and \(H\) is the reed tip opening at rest. The dimensionless control parameter is
\[
\gamma=\frac{P_m}{P_M},
\]
with analogous normalization for mouthpiece pressure and flow. The reed opening parameter is
\[
\zeta=\frac{Z_c U_A}{P_M}
= Z_c wH \sqrt{\frac{2}{\rho P_M}},
\]
where \(Z_c\) is the characteristic impedance, \(w\) the effective reed width, and \(\rho\) air density [1407.3547].

Using the wave decomposition
\[
p_n=p_n^+ + p_n^-,
\qquad
u_n=p_n^+ - p_n^-,
\]
and the lossless reflection law, the clarinet becomes a one-step non-autonomous iterated map
\[
p_n^+ = G(p_{n-1}^+,\gamma_n).
\]
The nonlinear reed characteristic is expressed in terms of the pressure drop \(\Delta p_n=\gamma_n-p_n\). In the standard Bernoulli-based piecewise law used in the paper,
\[
\Delta p_n \ge 1 \Rightarrow u_n=0,
\]
\[
0<\Delta p_n<1 \Rightarrow u_n=\zeta (1-\Delta p_n)\sqrt{\Delta p_n},
\]
\[
\Delta p_n \le 0 \Rightarrow u_n=-\zeta \sqrt{-\Delta p_n}.
\]
The factor \((1-\Delta p)\) captures the reduction of effective opening area as the reed approaches closure, while the reverse-flow branch omits that factor in the standard formulation used by the paper [1407.3547].

For constant \(\gamma\), the non-oscillating regime is a fixed point \(p^{+*}\) satisfying
\[
p^{+*}=G(p^{+*},\gamma).
\]
Under the assumptions of a lossless resonator, ideal spring reed, and \(\zeta<1\), the fixed point loses stability through a flip bifurcation at the static threshold
\[
\gamma_{st}=\frac{1}{3},
\]
beyond which the steady regime is a 2-cycle [1207.4636].

## 3. Static and dynamic oscillation thresholds under a linear pressure ramp

For a linear mouth-pressure ramp,
\[
\gamma_n=\gamma_0+\epsilon n,
\qquad 0<\epsilon\ll 1,
\]
the clarinet map becomes non-autonomous:
\[
p_n^+=G(p_{n-1}^+,\gamma_n).
\]
The key phenomenon is bifurcation delay: the orbit continues to track the unstable fixed-point branch past \(\gamma_{st}\), so the actual onset occurs at a dynamic threshold \(\gamma_{dt}\) strictly above the static threshold [1207.4636].

The dynamic analogue of the fixed point is the invariant curve \(\phi_\epsilon(\gamma)\), defined by
\[
\phi_\epsilon(\gamma)=G(\phi_\epsilon(\gamma-\epsilon),\gamma).
\]
For small \(\epsilon\), the invariant curve admits a perturbation expansion, with \(\phi_0(\gamma)=p^{+*}(\gamma)\). Bergeot et al. derive the threshold condition by examining the growth of deviations from this invariant curve. The deterministic dynamic threshold \(\gamma_{dt}^{th}\) is given by
\[
\int_{\gamma_0+\epsilon}^{\gamma_{dt}^{th}+\epsilon}
\ln\left|
\partial_x G\!\left(\phi_\epsilon(\gamma'-\epsilon),\gamma'\right)
\right|
\, d\gamma' = 0.
\]
In the deterministic regime, this threshold is largely independent of \(\epsilon\) provided \(\epsilon\) is small enough and the initial state is close to the stable branch [1407.3547].

The same analysis also yields a first-order approximation to the invariant curve in the linear-ramp case:
\[
\phi(\gamma,\epsilon)\approx p^{+*}(\gamma)
+\epsilon \frac{d p^{+*}(\gamma)}{d\gamma}
\frac{\partial_x G(p^{+*}(\gamma),\gamma)}
{\partial_x G(p^{+*}(\gamma),\gamma)-1},
\]
with
\[
p^{+*}(\gamma)=\frac{\zeta}{2}(1-\gamma)\sqrt{\gamma}.
\]
This formulation makes explicit that the onset delay is governed by the derivative \(\partial_x G\) evaluated along the invariant curve rather than by the static fixed-point criterion alone [1207.4636].

Finite numerical precision or physical noise suppresses bifurcation delay. Modeling perturbations as additive white noise with variance \(\sigma^2\), Bergeot et al. obtain the sweep-dominant estimator
\[
\hat{\gamma}_{dt}^{th}
=
\gamma_{st}
+
\sqrt{
-\frac{2\epsilon}{K}
\ln\!\left[
\left(\frac{\pi}{K}\right)^{1/4}
\frac{\sigma}{\epsilon^{5/4}}
\right]
},
\]
where \(K\) depends on the local slope of \(\partial_x G\) near the static bifurcation [1407.3547]. In this regime, larger \(\epsilon\) produces larger overshoot above \(\gamma_{st}\), while larger noise yields earlier onset and smaller delay. Numerical experiments also show strong sensitivity to precision: at low precision the delay collapses toward \(\gamma_{st}\), whereas very high precision recovers the deterministic prediction [1207.4636].

## 4. Exponential stabilization of mouth pressure and note-attack dynamics

A linear ramp is analytically convenient but physically incomplete because a musician does not increase mouth pressure indefinitely during a note attack. The 2014 extension therefore studies an exponential approach to a target pressure,
\[
\gamma_n=a\gamma_{n-1}+\gamma_M(1-a),
\qquad
\gamma_0=0,
\qquad
a=e^{-\epsilon},
\]
equivalently
\[
\gamma_n=\gamma_M(1-e^{-n\epsilon}),
\]
with \(\gamma_M=1\) in the paper’s examples [1407.3547].

The analysis is reduced to a linear sweep by introducing the auxiliary variable
\[
\eta(\gamma)=\ln\!\left(\frac{\gamma_M}{\gamma_M-\gamma}\right),
\qquad
\gamma(\eta)=\gamma_M(1-e^{-\eta}),
\]
so that
\[
\eta_n=\eta_{n-1}+\epsilon.
\]
Defining \(H(x,\eta)=G(x,\gamma(\eta))\), the dynamics becomes
\[
p_n^+=H(p_{n-1}^+,\eta_n),
\qquad
\eta_n=\eta_{n-1}+\epsilon.
\]
The invariant-curve condition in \(\eta\)-space is
\[
\phi_\epsilon(\eta)=H(\phi_\epsilon(\eta-\epsilon),\eta),
\]
and the corresponding deterministic threshold satisfies
\[
\int_{\eta_0+\epsilon}^{\eta_{dt}^{th}+\epsilon}
\ln\left|
\partial_x H\!\left(\phi_\epsilon(\eta'-\epsilon),\eta'\right)
\right|
\, d\eta' = 0.
\]
The mouth-pressure threshold is then recovered as
\[
\Gamma_{dt}^{th}=\gamma_M(1-e^{-\eta_{dt}^{th}}).
\]
A sweep-dominant expression of the same form is obtained for \(\hat{\eta}_{dt}^{th}\), with a constant \(K_\eta\) defined by the linearization of \(H\) near the bifurcation [1407.3547].

The comparison between linear and exponential profiles is one of the main substantive results. In deterministic simulations, representative thresholds are around \(0.9\) for a linear ramp and near \(0.7\) for the exponential profile when \(\zeta=0.5\), \(r=1\), and \(\gamma_M=1\). The paper interprets this as a consequence of the slowing of \(\gamma\)’s rate as it approaches the target, which reduces the bifurcation overshoot [1407.3547]. The numerical study identifies two regimes: DReg, with weak dependence on \(\epsilon\), and SDReg, with strong dependence on \(\epsilon\). Using a common rise-time measure \(N\) to reach \(99\%\) of \(\gamma_M\), the exponential profile reaches threshold in fewer time steps even when, in SDReg, its threshold in pressure may exceed that of the linear profile.

The numerical methodology also differs across the two clarinet papers. In the exponential-rise study, the dynamic threshold is the first time the orbit’s distance to the invariant curve exceeds \(\epsilon\), with simulations run in arbitrary precision arithmetic using mpmath in Python and effective noise level \(\sigma\approx 10^{-\text{precision}}\) [1407.3547]. In the earlier linear-ramp study, \(\gamma_{dt}^{num}\) is detected when the second-order difference of \(p_n^+\) changes sign between successive samples [1207.4636]. Both analyses underscore the same point: static thresholds alone are insufficient for note-attack transients.

## 5. PT as a pressure-amplified diffusion model on social networks

In social networks, the Pressure Threshold model is defined on a directed graph \(G=(V,E)\) with edge weights \(w_{uv}\in[0,1]\) and node thresholds \(\theta_v\in(0,1]\). Under the weighted-cascade initialization used in experiments,
\[
w_{uv}=\frac{1}{\mathrm{in\_deg}(v)}
\quad\text{for all }u\in N^-(v),
\]
so that initially \(\sum_{u\in N^-(v)} w_{uv}=1\) [2509.12822].

Each node has an activation state \(s_v(t)\in\{0,1\}\), with active set
\[
A_t=\{v\in V:s_v(t)=1\}.
\]
The process is progressive: once active, a node remains active. The incoming pressure at node \(v\) is
\[
p_v(t)=\sum_{u\in A_t,\,(u,v)\in E} w_{uv},
\]
and the activation rule is LT-style:
\[
s_v(t+1)=1 \quad \text{if } p_v(t)\ge \theta_v.
\]
The PT extension is the influence adjustment performed when a node activates. If \(v\) activates at time \(t\), let
\[
I_v=p_v(t)=\sum_{u\in A_t,\,(u,v)\in E} w_{uv}.
\]
For each outgoing neighbor \(s\) with \((v,s)\in E\) and \(s\notin A_t\),
\[
w_{vs}\leftarrow w'_{vs}=\min(1,\,w_{vs}+\alpha I_v),
\]
where \(\alpha\ge 0\) is the amplification parameter [2509.12822].

The update schedule is synchronous and two-phase. First, all inactive nodes are tested for activation by comparing \(p_v(t)\) to \(\theta_v\). Second, outgoing edges of the newly activated nodes are amplified toward inactive neighbors. The process halts at a fixed point when the newly activated set is empty. There is no global re-normalization after amplification: per-edge clamping ensures \(w_{vs}\le 1\), but the sum of in-neighbor weights at a recipient may exceed \(1\) after updates [2509.12822].

A recurrent misunderstanding is to treat PT as merely LT with different notation. The reduction
\[
\alpha=0 \Rightarrow w'_{vs}=w_{vs}
\]
shows instead that LT is a special case of PT. PT adopts the LT activation criterion but adds a state-dependent feedback loop through adaptive outgoing influence.

## 6. Influence maximization, structural properties, and empirical behavior

The influence-maximization problem under PT is to choose a seed set \(S\subseteq V\) with \(|S|\le k\) maximizing expected final spread,
\[
\max_{|S|\le k}\ \mathbb{E}[|\sigma_{PT}(S)|],
\]
where \(\sigma_{PT}(S)\) is the final active set under PT diffusion [2509.12822]. Because \(\alpha=0\) recovers LT, influence maximization under PT is NP-hard via reduction from LT. The spread function is monotone: if \(A\subseteq B\), then \(\sigma_{PT}(A)\subseteq \sigma_{PT}(B)\). However, PT is not submodular in general when \(\alpha>0\).

The non-submodularity is exhibited by a four-node counterexample with \(V=\{a,b,c,d\}\), edges \(a\to c\) and \(b\to c\) each of weight \(0.4\), edge \(c\to d\) of weight \(0.3\), thresholds \(\theta_c=0.8\) and \(\theta_d=0.4\), and \(\alpha=1.0\). For \(S=\{a\}\), only \(a\) remains active and \(\sigma(S)=1\). For \(S\cup\{c\}=\{a,c\}\), the spread is \(2\). For \(T=\{a,b\}\), node \(c\) activates with pressure \(0.8\), its edge to \(d\) is amplified to \(1.0\), and the final spread is \(4\). For \(T\cup\{c\}=\{a,b,c\}\), seeded \(c\) does not adjust its outgoing weight, so \(d\) remains inactive and the spread is \(3\). Hence
\[
\sigma(S\cup\{c\})-\sigma(S)=1,
\qquad
\sigma(T\cup\{c\})-\sigma(T)=-1,
\]
violating submodularity [2509.12822]. Greedy procedures such as CELF and CELF++ remain usable as heuristics, but the classical \((1-1/e)\) guarantee does not apply.

The experiments use Monte Carlo spread estimation with \(1{,}000\) simulations per evaluation. On Facebook with \(k=20\), CELF produces distinct seed sequences under LT and PT; the two agree early but diverge later, and PT introduces vertices such as \(1215\), \(1426\), \(3234\), \(2145\), and \(2832\) that never appear under LT [2509.12822]. Average influence also increases with \(\alpha\). Selected endpoints reported in the paper include: Facebook at \(k=60\), with LT spread \(2{,}461.219\), PT spread \(2{,}810.784\) for \(\alpha=0.001\), and \(3{,}490.280\) for \(\alpha=0.005\); Bitcoin at \(k=60\), with LT spread \(4{,}565.969\) and PT spread \(4{,}673.804\) for \(\alpha=0.001\); Wikipedia, where PT with \(\alpha=0.001\) reaches full coverage \(7{,}115\) by approximately \(k\approx 30\), and PT with \(\alpha=0.005\) reaches full coverage by \(k=5\); and an Erdős–Rényi network, where PT with \(\alpha=0.005\) reaches full coverage \(5{,}000\) by \(k=24\) [2509.12822].

The density effect is explained directly in terms of aggregate amplification. If \(N_t\) is the set of nodes newly activated at round \(t\), the total added incoming weight across their outgoing edges to inactive neighbors is
\[
\Delta W_t
=
\sum_{v\in N_t}
\sum_{s\in N^+(v),\,s\notin A_t}
\alpha I_v
=
\alpha
\sum_{v\in N_t}
I_v \cdot
\bigl|\{s\in N^+(v):s\notin A_t\}\bigr|.
\]
Approximating the inactive-neighbor count by \(\deg^+(v)\) early in the process gives
\[
\Delta W_t \approx \alpha \sum_{v\in N_t} I_v\,\deg^+(v).
\]
This explains why higher edge–node ratios, such as Facebook \(21.846\) and Wikipedia \(14.573\), display stronger PT amplification than the sparser Bitcoin network with edge–node ratio \(6.052\) [2509.12822].

## 7. Assumptions, limitations, and interpretive boundaries

The clarinet PT framework is intentionally idealized. The reed is treated as an ideal spring; reed motion-induced flow is neglected in the 2012 analysis; no reed mass, inertia, or contact dynamics are included; the bore is a lossless straight cylinder with perfect reflection; and the Bernoulli nonlinearity is quasi-steady [1207.4636]. In the 2014 exponential-rise study, losses are ignored by setting \(r=1\), and threshold predictions depend on the piecewise map \(G\) and especially on \(\partial_x G\) near the fixed point or invariant curve [1407.3547]. The sweep-dominant regime further shows that the dynamic threshold depends strongly on precision or noise, so measured onset can approach the static threshold even when the deterministic theory predicts substantial delay.

The network PT model also imposes restrictive assumptions. Diffusion is progressive, thresholds are fixed per run and independent across nodes, the graph is stationary during diffusion, and amplification is controlled by a single scalar \(\alpha\). Because updates are clamped edgewise rather than renormalized, large \(\alpha\) can produce saturated edges \(w_{uv}=1\) and may violate the initial in-weight normalization at recipients [2509.12822]. The lack of submodularity for \(\alpha>0\) removes standard approximation guarantees for greedy influence maximization, although the paper notes that approximate submodularity can hold empirically for small \(\alpha\) on large networks.

These limitations clarify the scope of the term. In clarinet acoustics, PT is a predictive framework for dynamic onset under prescribed mouth-pressure profiles, with the central distinction between static and dynamic oscillation thresholds. In network science, PT is an adaptive-threshold diffusion model in which activation pressure feeds back into subsequent influence transmission. The two usages share threshold logic and path dependence, but the literature treats them as separate model classes rather than as instances of a single unified formalism [1407.3547].

Source: https://www.emergentmind.com/topics/pressure-threshold-model-pt