---
title: 'Peeps: Gravitational Bursts and Peer Effects'
url: https://www.emergentmind.com/topics/peeps
type: topic
---

# Peeps: Gravitational Bursts and Peer Effects

“Peeps” is a polysemous technical term in recent research usage. In gravitational-wave astrophysics, a peep is the recurring, short-duration periapsis burst emitted during the early, very highly eccentric phase of an extreme mass-ratio inspiral (EMRI), when a stellar-mass compact object repeatedly sweeps through periapsis around a massive black hole yet remains silent for most of each long orbit [2507.19704]. In microeconomic theory, “peeps” is shorthand for peer effects that operate through either consideration sets, preference rankings, or both in dynamic discrete-choice environments, with identification based on how social interactions distort observed conditional choice probabilities [2310.12272]. The two usages are unrelated substantively, but both concern mechanisms whose individual realizations can be weak or intermittent while their aggregate consequences are analytically decisive [2305.05793].

## 1. Definitions and terminological scope

In the gravitational-wave literature, a peep is defined by recurrence. The underlying system is a stellar-mass compact object of mass \(\mu\sim 10\!-\!30\,M_\odot\) orbiting a massive black hole of mass \(M\sim 10^{5.2}-10^8\,M_\odot\), with mass ratio \(q\equiv \mu/M \sim 10^{-7}-10^{-4}\). After capture into a near-parabolic bound orbit, the object spends most of its orbital period far from the massive black hole and radiates little in the LISA band, then emits a brief millihertz burst at periapsis; because the orbit is bound with \(e\lesssim 1\) rather than exactly parabolic, the burst repeats once per orbit [2507.19704]. This distinguishes a peep both from an extreme-mass-ratio burst treated as a one-off parabolic flyby and from the later, quasi-continuous EMRI waveform near merger.

In the discrete-choice literature, “peeps” denotes peer effects in a model where social interactions can change either what an agent considers or how she ranks what she considers. The core objects are two directed edge sets, \(\Gamma=(\Gamma_C,\Gamma_R)\), where \(\Gamma_C\) affects consideration and \(\Gamma_R\) affects preferences. A peer can therefore affect only consideration, only preferences, or both [2310.12272]. This terminology is narrower than generic “social interactions,” because the paper ties it to a specific identification program based on behavioral signatures in observed choices.

A common misconception in both literatures is to collapse the term into a nearby but distinct category. In the astrophysical case, peeps are not merely parabolic bursts with slightly different notation; the bound recurrence is the central physical feature. In the economic case, peer effects are not restricted to preference spillovers; limited consideration is a distinct margin with different empirical implications [2507.19704].

## 2. Gravitational-wave peeps as an EMRI regime

The gravitational-wave peep arises in the standard EMRI capture channel through stellar dynamical interactions in galactic nuclei. Two-body or multi-body scattering occasionally places a compact object on an extremely eccentric bound orbit around a central massive black hole, after which gravitational radiation reaction removes energy and angular momentum and converts the capture into an inspiral [2507.19704]. The key observational fact is a strong separation of timescales: the orbital period is set mainly by the very large apoapsis of the nearly parabolic orbit, whereas the LISA-band emission is concentrated into a short periapsis passage.

This regime can persist for most of the inspiral lifetime. One explicit capture-stage example has \(T \approx 185000\ \text{years}\) for an orbit with \(a=0.9M\), \(p=120M\), and \(e=0.999999\), while a full inspiral example lasts about \(t \approx 6\times 10^8\ \text{years}\); early in that evolution, a 7-day zoom-in contains one peep, whereas near plunge the same interval contains about 115 peaks and the signal becomes nearly continuous [2305.05793]. The early peep phase is therefore not a different source class from EMRIs, but an early dynamical regime of the same system.

The morphology is governed primarily by periapsis conditions. The emitted waveform can exhibit zoom-whirl or homoclinic behavior, with a long-period structure plus a short, rapid periapsis feature, and the amplitude evolves only weakly over much of the inspiral because the periapsis remains roughly fixed while the apoapsis shrinks [2507.19704]. This is why peeps can remain astrophysically long-lived yet individually faint.

The distinction from parabolic-burst modeling is consequential. For fixed parameters, a bound source with \(e=0.99\) can complete 3 orbits in one year and emit 3 peeps, while an almost parabolic source emits only one burst; in the example reported, the peep spectrum has a peak amplitude about 4 times larger and extends to higher frequencies, and at \(p=15M\) the peep source gives about 50 complete orbits in one year [2305.05793]. Approximating peeps as single-pass parabolic bursts therefore undercounts the finite-time frequency-domain content relevant to space-based detection.

## 3. Modeling peeps for LISA and confusion-noise estimation

The waveform model used for peeps is the Numerical Kludge (NK) model of Babak et al. (2007). The parameter space is reduced to the 12 quantities
\[
\{M,\mu,a,p_0,e_0,\iota_0,d_l,\theta_s,\phi_s,\theta_k,\phi_k,\chi\},
\]
with redshift \(z\) as input and conversion to luminosity distance \(d_l\). The workflow is to define the orbit through the constants of motion \((E,L_z,Q)\), integrate the Kerr geodesic equations to obtain the trajectory in Boyer–Lindquist coordinates \(\{t,r,\theta,\phi\}\), map that trajectory into flat-space Cartesian coordinates, and then apply the quadrupole formula to compute the metric perturbation [2507.19704]. The model is taken to retain sufficient strong-field structure for the peep regime, especially for periapses \(r_p\ge 5M\), corresponding roughly to \(p\gtrsim 10M\) at peep eccentricities.

The source population is built from the Illustris-1 cosmological hydrodynamic simulation out to \(z=3.01\), with the EMRI-relevant massive-black-hole mass range restricted to
\[
\log_{10}(M/M_\odot)=5.2-8.0,
\]
mass functions tabulated in \(0.1\)-dex bins, and shell-by-shell integration in \(dz=0.1\) increments beyond the local \(z=0\) volume [2507.19704]. Baseline EMRI rates are taken from Babak et al. (2017), using the median rate from their Eqs. 23–31 and plunge-to-EMRI ratio \(N_p=10\). For each massive black hole in the synthesized population, a compact-object mass is drawn, a 4-year EMRI occurrence probability is assigned, orbital and angular parameters are sampled, an NK waveform is generated, the signal is cosmologically redshifted via
\[
t_{\rm det} = (1+z)\, t_{\rm src}, \qquad h_{\rm det} = \frac{h_{\rm src}}{1+z},
\]
and the result is passed through the LISA response using `fastlisaresponse` with second-generation time-delay interferometry and ESA trailing orbits, retaining only the \(A\) and \(E\) channels because \(T\) is approximately null [2507.19704].

The study does not construct the background via the standard \(\Omega_{\rm GW}(f)\) or \(dE_s/df_s\)-based stochastic-background formalism. Instead it performs direct time-domain Monte Carlo superposition of many individually unresolvable detector responses. Characteristic strain is taken as
\[
h_c(f) = 2 f\left(\left|\tilde{X}(f)\right|^2\right)^{1 / 2},
\]
and single-channel SNR is evaluated with the standard matched-filter-like expression
\[
\mathrm{SNR}^2 = 4\int_0^{f_{\max}} df\, \frac{|\tilde X(f)|^2}{S_n(f)},
\]
with the \(A\) and \(E\) channel SNRs then combined in quadrature [2507.19704].

Three background scenarios summarize the astrophysical uncertainty. Background 1 assumes \(15M \le p_0 \le 120M\), \(0.99 \le e_0 \le 0.99999\), and at most one highly eccentric EMRI per massive black hole during the mission. Background 2 keeps the one-or-fewer-per-massive-black-hole assumption but broadens the capture distribution to \(8M \le p_0 \le 120M\) and \(0.9 \le e_0 \le 0.999999\). Background 3 adopts Background 2 and then multiplies it by \(\sqrt{1000}\) to approximate an incoherent sum of 1000 signals per host, motivated by the much more abundant early-EMRI picture discussed by Amaro-Seoane et al. (2024) [2507.19704].

The reported Monte Carlo realization generates 1470 unique peep waveforms. In Background 1, the per-channel SNRs are of order \(0.2\) and the combined SNR is \(0.33\). In Background 2, the per-channel SNRs are of order \(1.7\) and the combined SNR is \(2.4\). In Background 3, the \(A\) and \(E\) channels each have SNR \(\sim 55\), with combined SNR \(\sim 77\) [2507.19704]. Under the first two, peeps produce only a slight rise in the LISA noise floor; under the third, they would constitute a directly detectable foreground capable of obscuring other sources.

## 4. “Peeps” in dynamic discrete choice

In the economic model, there is a finite set of agents
\[
\mathcal A=\{1,2,\dots,A\},
\]
a finite menu
\[
\mathcal Y=\{0,1,\dots,Y\},
\]
and a choice configuration \(\mathbf y=(y_a)_{a\in\mathcal A}\in\mathcal Y^A\). Time is continuous, and each agent \(a\) receives revision opportunities according to an independent Poisson clock with rate \(\lambda_a\) [2310.12272]. When a clock rings, the decision problem is two-stage: the agent first forms a consideration set \(\mathcal C\subseteq\mathcal Y\), then chooses an option from \(\mathcal C\).

The model’s distinctive feature is the separation of social influence into two directed networks. For agent \(a\), \(\mathcal{NC}_a\) denotes peers affecting consideration and \(\mathcal{NR}_a\) denotes peers affecting preferences. Consideration for alternative \(v\) is summarized by \(Q_a(v\mid \mathbf y,\mathcal{NC}_a)\), with the baseline restriction that inclusion is independent across alternatives:
\[
C_a(\mathcal C\mid \mathbf y,\mathcal{NC}_a,\mathcal Y) = \prod_{v\in\mathcal C} Q_a(v\mid \mathbf y,\mathcal{NC}_a) \prod_{v\in\mathcal Y\setminus\mathcal C}\bigl(1-Q_a(v\mid \mathbf y,\mathcal{NC}_a)\bigr),
\]
and \(Q_a(0\mid \mathbf y,\mathcal{NC}_a)=1\). For \(v\neq 0\), the probability of considering \(v\) depends on the number, not the identity, of consideration peers choosing \(v\):
\[
Q_a(v\mid \mathbf y,\mathcal{NC}_a)\equiv Q_a\bigl(v\mid NC_a^v(\mathbf y)\bigr).
\]
This makes consideration effects alternative-specific and non-crossing.

Conditional on \(\mathcal C\), the agent chooses according to
\[
R_a(\cdot\mid \mathbf y,\mathcal{NR}_a,\mathcal C),
\]
again with an anonymity restriction: the rule depends on counts of peers choosing each alternative in \(\mathcal C\), not their identities. Unlike consideration, preference effects are intrinsically cross-alternative because peer shifts among alternatives change relative desirability. The paper gives a logit example,
\[
R_a(v\mid \mathbf y,\mathcal{NR}_a,\mathcal C) = \frac{\exp\bigl(u_{a,v,\mathcal C}(\mathbf y,\mathcal{NR}_a)\bigr)} {\sum_{v'\in\mathcal C}\exp\bigl(u_{a,v',\mathcal C}(\mathbf y,\mathcal{NR}_a)\bigr)},
\]
and explicitly allows \(\mathcal C\) itself to affect utilities, so IIA can fail [2310.12272].

The observed conditional choice probability is
\[
P_a(v\mid \mathbf y) = \sum_{\mathcal C\subseteq\mathcal Y} R_a(v\mid \mathbf y,\mathcal{NR}_a,\mathcal C) \prod_{v'\in\mathcal C}Q_a(v'\mid \mathbf y,\mathcal{NC}_a) \prod_{v'\in\mathcal Y\setminus\mathcal C}\bigl(1-Q_a(v'\mid \mathbf y,\mathcal{NC}_a)\bigr),
\]
which factorizes into the probability that \(v\) is considered times the probability that \(v\) is chosen conditional on being considered [2310.12272]. This decomposition is the basis for identification.

## 5. Identification, recovery, and empirical implementation

Because only one agent updates at a time almost surely, the model induces a continuous-time Markov chain over configurations. Under positive consideration and positive choice probabilities, the paper establishes existence of a unique full-support invariant distribution \(\mu\) [2310.12272]. More important than equilibrium existence, however, is the identification logic: consideration effects are separable and alternative-specific, whereas preference effects generate cross-order effects across alternatives.

The first step is to detect whether another agent \(a'\) is any kind of peer by examining
\[
\Delta_{a'}^v \ln P_a(v\mid \mathbf 0),
\]
where \(\mathbf 0=(0,\dots,0)\). If this is nonzero, \(a'\) belongs to the relevant peer set. The second step uses the cross-alternative double difference
\[
\Delta_{a''}^w \Delta_{a'}^v \ln P_a(v\mid \mathbf 0), \qquad w\neq v.
\]
Because the consideration term for \(v\) does not depend on peers choosing \(w\neq v\), this object vanishes for consideration-only links and is nonzero when the link affects preferences. A further same-alternative double-difference argument, together with the exclusion condition
\[
|\mathcal{NC}_a\setminus\mathcal{NR}_a| + |\mathcal{NR}_a\setminus\mathcal{NC}_a| \ge 1,
\]
separates preference-only links from links affecting both margins [2310.12272].

Once the network structure is known, the paper shows how to recover ratios of consideration probabilities, synthetic reduced-menu conditional choice probabilities, and—under sufficient support—both the full consideration mechanism \(Q_a\) and the choice rule \(R_a\). The key constructive device is that variation in consideration-only peers mimics menu variation. If a consideration-only peer switches to \(v'\), only the chance that \(v'\) enters the consideration set changes; preferences do not. This makes it possible to identify counterfactual reduced-menu probabilities without actual menu variation [2310.12272].

The paper also distinguishes two data environments. With continuous-time observations, the transition-rate matrix is directly observable, so \(\lambda_a P_a(v\mid \mathbf y)\) is identified and summing over \(v\) recovers \(\lambda_a\). With discrete-time panel observations at fixed interval \(\Delta\), the transition matrix satisfies \(\mathcal P(\Delta)=e^{\Delta \mathcal M}\), and under a generic distinct-eigenvalue condition from Blevins (2018), \(\mathcal M\), \(P\), and \(\lambda\) are generically identified [2310.12272].

The empirical application studies Starbucks and Luckin expansion across 152 Chinese markets, treating each firm-market pair \(a=(f,m)\) as an agent in a binary-choice environment. Preference peers are assumed known from industrial-organization structure and are same-market only, whereas consideration peers can come from neighboring markets and are estimated. The recovered consideration network has 960 directed links, average in-degree and out-degree about 6.3, and about 30% of markets have no incoming consideration effects from other markets, though almost all affect some other market’s consideration [2310.12272]. The main findings are that Starbucks is close to full consideration, Luckin exhibits substantial limited consideration, ignoring limited consideration misreads profitability, and forcing full consideration causes duopolies to emerge about 17 months sooner.

## 6. Caveats, misconceptions, and broader analytical significance

In gravitational-wave analysis, the central caveat is that the peep “background” is not a standard isotropic Gaussian stochastic gravitational-wave background derived from \(\Omega_{\rm GW}(f)\)-type integrals. It is a direct waveform-summation estimate of a confusion foreground generated by many individually unresolvable deterministic sources [2507.19704]. Over the affected band, the frequency range is \(\Omega \sim 4\times 10^{-3}\ {\rm Hz}\), the 4-year bin width is \(\Delta f \sim 8\times 10^{-9}\ {\rm Hz}\), and the number of bins is therefore \(\Omega/\Delta f \sim 0.5\times 10^6\). In the first two scenarios there are only \(\mathcal O(10^3)\) peep signals, far fewer than the number of bins, so the foreground is better described as sparse, burst-train, or “popcorn-like” than as naturally Gaussianized noise [2507.19704]. Additional uncertainties arise from EMRI rate prescriptions, assumed \(p_0\) and \(e_0\) distributions, Illustris massive-black-hole demographics, the neglect of orbital evolution over 4 years, and the approximation of the most abundant case by multiplying Background 2 by \(\sqrt{1000}\).

In the discrete-choice setting, the main limitation is that the clean cross-alternative logic requires \(Y\ge 2\). In the binary case \(Y=1\), one can identify the total peer set but cannot separate consideration and preference effects without additional structure, such as one known network component, stronger network restrictions, or a non-exponential curvature condition for preference ratios [2310.12272]. The baseline recovery of consideration probabilities also relies on the independent-inclusion consideration structure, although the paper provides a separate large-support identification result for a fully general consideration mechanism using excluded covariates that affect only preferences.

A plausible commonality across the two literatures is methodological rather than substantive. In both, the relevant micro-events are individually easy to misclassify as negligible: a single peep waveform is usually unresolvable, and a single peer perturbation can be confounded with standard utility variation. Yet in both cases the aggregate object—confusion noise in one literature, nonparametrically identified network structure in the other—depends precisely on taking those weak or intermittent effects seriously [2507.19704].

Source: https://www.emergentmind.com/topics/peeps