# ⚛  L1 Principle — Photoacoustic-Ultrasound Dual-mode Imaging (PAUS)

**ID:** `L1-509` · **Status:** ⊙ Testnet (genesis catalog)

> **🌐 Domain:** Medical Imaging — *Concurrent functional + anatomical imaging via dual-contrast acquisition (multi-physics joint inverse)*
> **🎯 Problem class:** nonlinear inverse · **🧮 Solution space:** 4D dual contrast functional anatomical
> **📡 Carrier:** photon_acoustic_dual · **🌫 Noise:** gaussian
> **⚖ Difficulty (δ):** 5 · **⛓ Block:** 41553387

---

## 🧠 1. Introduction

**Photoacoustic-Ultrasound Dual-mode Imaging (PAUS)** is a **nonlinear inverse problem** whose unknown lives in **4D dual contrast functional anatomical** space, within the **Concurrent functional + anatomical imaging via dual-contrast acquisition (multi-physics joint inverse)** sub-domain of **Medical Imaging**.

Measurements consist of photon acoustic dual via a **joint thermoelastic and pulse echo** sensing mechanism.

The forward operator applies, in order: L · optical source operator; L · optical diffusion operator; exponential attenuation along the propagation path; L · thermoelastic grueneisen operator; L · acoustic wave propagation operator; L · us transmit beam operator; L · us pulse echo operator; L · transducer detection dual operator; int · spectral temporal spatial operator.

Observations are corrupted by additive Gaussian noise. Existence of recovered joint state (mu_a, Z, c, rho) is guaranteed within the declared Omega bounds. Uniqueness holds under multi-illumination + multi-wavelength PA acquisition combined with US transmit-receive at adequate SNR; single-wavelength single-illumination PAUS suffers the same multiplicative non-uniqueness as single-wavelength qPAT (excluded from the spec range). Stability is moderately conditioned (kappa_eff ~ 40 after joint model-based reconstruction) — pa_us_synchronization_jitter dominates registration error; fluence_inhomogeneity dominates qPAT-class chromophore bias; acoustic_attenuation dominates US-class deep-tissue contrast; transducer_bandwidth_difference_pa_vs_us contributes a frequency-dependent reconstruction artifact (PA typically lower-frequency band 0.5-10 MHz, US higher 5-20 MHz). Joint Hadamard well-posedness for the coupled qPAT + US Born forward is established by the qPAT references (Bal-Uhlmann 2010, Bal-Ren 2011) plus the dual-mode literature: Niederhauser et al. 2005 (combined OAT/US), Beard 2011 (review), Wang LV 2008 (multiscale photoacoustic), Mehrmohammadi et al. 2013 (PA-US clinical), Park et al. 2017 (PAUS guided interventions).

## ⚙ 2. Forward Model

Physical chain: **x** → L · optical source → L · optical diffusion → Beer-Lambert attenuation → L · thermoelastic grueneisen → L · acoustic wave propagation → L · us transmit beam → L · us pulse echo → L · transducer detection dual → int · spectral temporal spatial → **y** (detector).

```
y = `int.spectral_temporal_spatial` `L.transducer_detection_dual` `L.us_pulse_echo` `L.us_transmit_beam` `L.acoustic_wave_propagation` `L.thermoelastic_grueneisen` exp(-∫µ dl) `L.optical_diffusion` `L.optical_source` x + n,    n ~ 𝒩(0, σ²)
```

**Measurement DAG:**

| Primitive | What it does |
|---|---|
| `L.optical_source` | L · optical source operator |
| `L.optical_diffusion` | L · optical diffusion operator |
| `L.beer_lambert` | Exponential attenuation along the propagation path |
| `L.thermoelastic_grueneisen` | L · thermoelastic grueneisen operator |
| `L.acoustic_wave_propagation` | L · acoustic wave propagation operator |
| `L.us_transmit_beam` | L · us transmit beam operator |
| `L.us_pulse_echo` | L · us pulse echo operator |
| `L.transducer_detection_dual` | L · transducer detection dual operator |
| `int.spectral_temporal_spatial` | Int · spectral temporal spatial operator |

## 🔬 3. Physics Fingerprint

| Property | Value |
|---|---|
| Domain | Medical Imaging |
| Sub domain | Concurrent functional + anatomical imaging via dual-contrast acquisition (multi-physics joint inverse) |
| Carrier | photon_acoustic_dual |
| Problem class | nonlinear_inverse |
| Solution space | 4D_dual_contrast_functional_anatomical |
| Noise model | gaussian |
| Integration axis | spectral_temporal_spatial |
| Difficulty delta | 5 |
| L dag | 9.5 |

## 📡 4. Measurement Model

Existence of recovered joint state (mu_a, Z, c, rho) is guaranteed within the declared Omega bounds. Uniqueness holds under multi-illumination + multi-wavelength PA acquisition combined with US transmit-receive at adequate SNR; single-wavelength single-illumination PAUS suffers the same multiplicative non-uniqueness as single-wavelength qPAT (excluded from the spec range). Stability is moderately conditioned (kappa_eff ~ 40 after joint model-based reconstruction) — pa_us_synchronization_jitter dominates registration error; fluence_inhomogeneity dominates qPAT-class chromophore bias; acoustic_attenuation dominates US-class deep-tissue contrast; transducer_bandwidth_difference_pa_vs_us contributes a frequency-dependent reconstruction artifact (PA typically lower-frequency band 0.5-10 MHz, US higher 5-20 MHz). Joint Hadamard well-posedness for the coupled qPAT + US Born forward is established by the qPAT references (Bal-Uhlmann 2010, Bal-Ren 2011) plus the dual-mode literature: Niederhauser et al. 2005 (combined OAT/US), Beard 2011 (review), Wang LV 2008 (multiscale photoacoustic), Mehrmohammadi et al. 2013 (PA-US clinical), Park et al. 2017 (PAUS guided interventions).

| Metric | Value |
|---|---|
| Metric | PSNR_dB |
| Secondary | RMSE_per_modality |

## 📏 5. Operating Range (Ω)

**Center problem class:** `paus_dual_contrast` · **Forward operator:** `paus_joint_forward`

**Center point:**

| Parameter | Unit | Value |
|---|---|---|
| H | px | 256 |
| W | px | 256 |
| Z | — | 64 |
| Snr pa db | dB | 22 |
| Snr us db | dB | 28 |
| N us modes | — | 2 |
| Voxel size mm | mm | 0.2 |
| Lambda range nm | nm | 700 – 1064 |
| N pa wavelengths | — | 5 |
| Us bandwidth mhz | MHz | 5 |
| Speed of sound mps | — | 1500 |
| Acoustic attenuation | — | 0 |
| Fluence inhomogeneity | — | 0 |
| Us center frequency mhz | MHz | 7.5 |
| Motion during acquisition | — | 0 |
| Acoustic speed heterogeneity | — | 0 |
| Pa us synchronization jitter | — | 0 |
| Transducer bandwidth difference pa vs us | µs | 0 |

**Allowed bounds:**

| Parameter | Unit | Range |
|---|---|---|
| H | px | 128 – 512 |
| W | px | 128 – 512 |
| Z | — | 32 – 256 |
| Snr pa db | dB | 5 – 40 |
| Snr us db | dB | 5 – 45 |
| N us modes | — | 1 – 5 |
| Voxel size mm | mm | 0.05 – 1.0 |
| N pa wavelengths | — | 2 – 16 |
| Us bandwidth mhz | MHz | 1 – 20 |
| Speed of sound mps | — | 1400 – 1600 |
| Acoustic attenuation | — | 0.0 – 0.4 |
| Fluence inhomogeneity | — | 0.0 – 0.4 |
| Us center frequency mhz | MHz | 1 – 30 |
| Motion during acquisition | — | 0.0 – 0.3 |
| Acoustic speed heterogeneity | — | 0.0 – 0.1 |
| Pa us synchronization jitter | — | 0.0 – 0.3 |
| Transducer bandwidth difference pa vs us | µs | 0.0 – 0.5 |

## 🎯 6. Tolerance (ε)

**Center tolerance:** 26.0

| Metric | Range |
|---|---|
| Psnr db | 8.0 – 42.0 |

## ⚖ 7. Hardness Function

Hardness scales as **`epsilon_fn`** on **PSNR_dB**, with κ = `280` and δ = `5`.

## 💾 8. Reference Dataset

- **primary** · weight 1.0 · IPFS _(not pinned yet)_

## 9. On-chain Registration

- **Chain hash:** `0xedc6787d7a57c3ffbe0fcc816d4c1b47dd38ddfa2914b2392b57b48d19033e1d`
- **Chain tx hash:** `0x6d485b3261bc469ae58032e3f62431adb715c465c744d975ddc8ec63b9b20fc0`
- **Chain block:** `41553387`

---

## File Mapping

This bundle consists of: `L1-509.md`, `L1-509.json`.

| File | Role | How to regenerate |
|------|------|-------------------|
| `L1-509.md` | Source of truth — edit this | Human or LLM |
| `L1-509.json` | Structured metadata for the registry | LLM regenerates from the sections above |

**Prompt for your LLM after editing this Markdown:**

> Read the attached Markdown. Regenerate the sibling `.json` so every field matches.
> Preserve the schema documented in the rows above.
> Output each file in its own fenced code block tagged with the filename.
> Output only the JSON object.

_This Markdown was auto-synthesized from the catalog row for `L1-509`._
_Edit it, regenerate the JSON, and submit at [/submit](/submit) to claim the artifact._