---
title: 'FAUST: Protostellar Science, Formal Methods, DSP'
url: https://www.emergentmind.com/topics/faust
type: topic
---

# FAUST: Protostellar Science, Formal Methods, DSP

Searching arXiv for FAUST-related papers to ground the article in current research.
{"query":"FAUST arXiv ALMA Large Program FAUST² Faust compiler", "max_results": 10}
In the literature summarized here, **FAUST** denotes several distinct technical enterprises rather than a single framework. The dominant usage is the ALMA Large Program **“Fifty AU STudy of the chemistry in the disk/envelope system of Solar-like protostars”**, which was designed to resolve the chemistry, kinematics, and physical structure of young protostellar systems on scales of roughly \(50\) au. Separate usages denote **FAUST\(^2\)**, a MATLAB toolbox for formal abstractions of uncountable-state discrete-time Markov processes, and **Faust**, a purely functional, block-diagram-oriented language for real-time audio and DSP whose recent compiler work includes automatic fixed-point format determination and CLAP plugin generation [2111.14121] [1403.3286] [2403.06527] [2606.13193].

## 1. FAUST as an ALMA Large Program in protostellar science

The ALMA FAUST program was introduced to reveal and quantify the variety of chemical composition of the envelope/disk system at scales of \(\sim 50\) au in a sample of Class 0 and I protostars representative of the chemical diversity observed at larger scales. Its stated objectives were to disentangle the components of the \(50\)–\(2000\) au envelope/disk system, characterize organic complexity in each component, probe the ionization structure, and measure molecular deuteration. The original program description specifies a sample of **13 Class 0/I protostars** at distances \(\le 250\) pc and luminosities \(\le 25\,L_\odot\), observed with three frequency setups—\(214\)–\(219\) GHz and \(229\)–\(234\) GHz, \(242.5\)–\(247.5\) GHz and \(257.5\)–\(262.5\) GHz, and \(85\)–\(89\) GHz and \(97\)–\(101\) GHz—at \(0.2\) km s\(^{-1}\) spectral resolution and \(250\)–\(350\) mas angular resolution, using both the 12 m array and the 7 m ACA. One spectral window per setup was dedicated to continuum measurements from \(3\) mm to \(1\) mm [2111.14121].

The molecular strategy was explicitly partitioned into tracer groups. Envelope probes included **c-C\(_3\)H\(_2\)** and **CS**; centrifugal-barrier probes included **CH\(_3\)OH**, **SO**, and **SiO**; disk probes included **H\(_2\)CO**, **C\(^{18}\)O**, and **HC\(_3\)N**. Organic-complexity diagnostics targeted **CH\(_3\)OH**, **NH\(_2\)CHO**, **CH\(_3\)CHO**, **CH\(_3\)OCH\(_3\)**, and **HCOOCH\(_3\)**. Ionization diagnostics used **H\(^{13}\)CO\(^+\)**, **DCO\(^+\)**, and **N\(_2\)H\(^+\)**, while deuteration diagnostics included **c-C\(_3\)HD**, **N\(_2\)D\(^+\)**, **HDCO**, **D\(_2\)CO**, and **CH\(_2\)DOH** [2111.14121].

Later program-context descriptions retain the same \(\sim 50\) au emphasis while broadening the characterization of the survey. One such summary defines FAUST as a program designed to resolve and characterize young Class 0/I disks and their infall envelopes on \(\sim 50\) au scales, combining multi-frequency continuum and spectral-line mapping to study both chemistry and dust growth, and states that the sample includes **\(\gtrsim 16\)** Sun-like protostars [2606.22489]. This suggests that FAUST functions both as a tightly specified observing program and as an evolving publication series centered on uniform, high-resolution, multi-scale studies of young Solar-like systems.

## 2. Chemical regimes, molecular complexity, and deuteration in the FAUST series

A central result of the astronomy FAUST papers is that chemically rich inner regions are not restricted to the earliest Class 0 phase. **FAUST I** reported a hot corino toward the Class I protostar **L1551 IRS5**, with methanol non-LTE analysis giving a gas temperature of \(\sim 100\) K, density \(\ge 1.5\times 10^8\) cm\(^{-3}\), and an emitting radius of \(\sim 10\) au. The methyl formate and ethanol relative abundances were found to be compatible with those measured in Class 0 hot corinos, and the study concluded that little chemical evolution from Class 0 to I hot corinos occurs [2007.10275].

FAUST also identifies sources in which hot-corino chemistry coexists with carbon-chain chemistry. In **CB68**, CH\(_3\)OH, HCOOCH\(_3\), and CH\(_3\)OCH\(_3\) were detected toward the protostar, with a CH\(_3\)OH rotation temperature of \(131\pm 11\) K and an emitting region of \(\sim 10\) au, while c-C\(_3\)H\(_2\) and CCH were detected on a \(1000\) au scale. The source was therefore described as having a **hybrid chemistry** [2206.06603]. In **VLA 1623–2417 B**, compact CH\(_3\)OH emission and methyl formate were detected; LVG analysis yielded a size of \(0.11\)–\(0.34\) arcsec (\(14\)–\(45\) au), \(N({\rm CH_3OH}) = 10^{16}\)–\(10^{17}\) cm\(^{-2}\), kinetic temperature \(>170\) K, and volume density \(>10^8\) cm\(^{-3}\), while the line profiles were interpreted in terms of a chemically enriched ring of radius \(33\) au close to the centrifugal barrier [2206.13339].

Several FAUST studies connect chemical enrichment to shocks and filamentary accretion. In **[BHB2007] 11**, more than **45 CH\(_3\)OCHO lines**, **8 CH\(_3\)OCH\(_3\)** transitions, one H\(_2\)CCO transition, and four trans-HCOOH transitions were detected, and the compact iCOM emission was found to encompass both protostars while tending to align with the southern filament. The study concluded that the detected methanol and the other iCOMs are generated by shocked gas from incoming filaments streaming toward the two sources [2403.07757]. In **Elias 29**, SO rotational temperatures reached \(53^{+25}_{-15}\) K at the interaction point of the outflow and the southern ridge and \(72^{+66}_{-29}\) K within the southeastern outflow, while \(T_{\rm rot}({\rm SO})\) remained \(20\)–\(30\) K in the quiescent southern ridge; the warm condition was attributed to the nearby B-type star HD147889, whereas local hot spots were associated with outflow or jet interactions [2501.10634].

Deuteration is a second major axis of the FAUST chemical program. **FAUST X** measured H\(_2\)CO, HDCO, and D\(_2\)CO toward **[BHB2007] 11**, deriving \(N({\rm H_2CO}) \simeq (3\)–\(8)\times 10^{14}\) cm\(^{-2}\), \(N({\rm HDCO}) \simeq (0.8\)–\(2.9)\times 10^{13}\) cm\(^{-2}\), and \(N({\rm D_2CO}) \simeq (2.6\)–\(4.3)\times 10^{12}\) cm\(^{-2}\), corresponding to an average D/H ratio of \(0.01\)–\(0.04\); the same work tentatively interpreted a second large-scale H\(_2\)CO feature as an asymmetric molecular outflow launched by a wide-angle disk wind [2309.00501]. **FAUST XVII** reported the first detection of deuterated formaldehyde in a planet-forming disk, **IRS 63**, with \([{\rm HDCO}]/[{\rm H_2CO}] \sim 0.1\)–\(0.3\) and \([{\rm D_2CO}]/[{\rm H_2CO}] \sim 0.1\); notably, D\(_2\)CO was strongly asymmetric and peaked where the streamer strikes the disk, and it was also detected in two outflow spots, leading to an interpretation in which HDCO is dominated by gas-phase formation while D\(_2\)CO is mainly grain-mantle material released by shocks [2407.04813]. **FAUST-XXII** extended this approach to **VLA 1623–2417**, finding \(T \sim 125\) K toward the hot corino, \(20\)–\(40\) K in the outflow cavities, and \(\le 15\) K in the streamers, with \([{\rm HDCO}]/[{\rm H_2CO}] \sim 0.16\) in the hot corino, \(\sim 0.07\)–\(0.13\) in the outflow cavities, and \(\sim 0.3\) in a streamer. The similar deuteration values across components were interpreted as evidence that prestellar deuteration is inherited mostly unaltered into the protostellar phase [2502.18675].

Taken together, these results show that FAUST chemistry is not limited to cataloging line detections. It uses spatially resolved abundances, isotopologue ratios, and excitation diagnostics to distinguish hot corinos, warm carbon-chain chemistry, shock-enhanced chemistry, grain-surface inheritance, and gas-phase reprocessing on scales from tens to thousands of au.

## 3. Disk-envelope coupling, streamers, outflows, and multiplicity

FAUST astronomy also emphasizes kinematic continuity and discontinuity across scales. In **VLA 1623AB**, the envelope rotation axis traced by H\(^{13}\)CO\(^+\) was measured as \(127^\circ \pm 4^\circ\), whereas the circum-binary disk minor axis of VLA 1623A was \(115^\circ \pm 4^\circ\), implying a misalignment of about \(12^\circ \pm 6^\circ\). The outflow velocity gradient seen in CCH is perpendicular to the outflow axis and agrees with the envelope rotation axis rather than the disk minor axis; under the assumption of outflow rotation with constant specific angular momentum \(j \approx 100\) au km s\(^{-1}\), the launching radius was estimated as \(5\)–\(16\) au. The same study detected, for the first time, a velocity gradient associated with rotation toward the VLA 1623B disk, with a sense opposite to that of the large-scale envelope, outflow, and circum-binary disk [2201.07334].

A subsequent FAUST study of the same protocluster detected an accelerating **SO streamer** plausibly feeding VLA 1623B. The streamer extends over \(\sim 1500\)–\(2000\) au, has \(T_{\rm ex}=33\pm 9\) K, \(N_{\rm SO}\simeq 2\times 10^{14}\) cm\(^{-2}\), \(X({\rm SO})\simeq 1\times 10^{-8}\), a total mass of \(3\times 10^{-3}\,M_\odot\), and an accretion rate of \(3\)–\(5\times 10^{-7}\,M_\odot\,{\rm yr}^{-1}\). That rate is close to the inferred mass accretion rate of VLA 1623B itself, which the study used to argue that asymmetric streamer infall is an important contributor to protostellar disk growth. The same work reported the first SiO detection toward VLA 1623–2417, identifying a compact \(100\)–\(200\) au jet around source B and shock-heated high-velocity SO with \(T_{\rm ex}\simeq 92\pm 18\) K and \(102\pm 19\) K in the red and blue components [2402.10258].

FAUST has also been applied to binaries with circumbinary structure and entangled outflow morphology. In **L1551 IRS 5**, continuum and C\(^{18}\)O analysis yielded a circumbinary disk mass of \(0.018\,M_\odot\), northern and southern circumstellar disk masses of \(0.004\,M_\odot\) and \(0.002\,M_\odot\), a centrifugal barrier radius \(r_{CB}=55\) au, and a specific angular momentum \(j\sim 270\) au km s\(^{-1}\). The line-of-sight velocity field required an infalling-rotating envelope rather than a pure Keplerian curve, and an analytic jet-driven outflow model fit the broad X-shaped C\(^{18}\)O cavities with density power-law index \(\alpha \simeq 1.7\) and envelope rotation velocity \(V_c \simeq 2\) km s\(^{-1}\) [2506.11363].

A related L1551 IRS 5 study resolved two continuum over-densities at the edge of the circumbinary cavity, with the northern feature being about **20% brighter** than the southern one. Using C\(^{18}\)O kinematics and previous astrometry, it derived \(M_b \simeq 1.1\,M_\odot\), semimajor axis \(a_b = 42.2\) au, eccentricity \(e_b = 0.175\), inclination \(i_b = 122.3^\circ\), longitude of ascending node \(\Omega = 168.9^\circ\), and argument of periastron \(\omega = 176.5^\circ\). Three-dimensional gas-dust SPH simulations with **Phantom**, followed by MCFOST post-processing, reproduced the brightness contrast and the S-shaped velocity twist, supporting a binary–disc interaction origin for the dust concentration [2512.10486].

Another FAUST development is methodological rather than purely descriptive. **FAUST XXIX** proposed a luminosity estimator for embedded protobinaries based on the sublimation radius of **OCS**, aided by quantum-mechanical calculations of the OCS binding-energy distribution. Applied to **NGC 1333 IRAS 4A**, the method yielded luminosities of \(7.5\pm 2.5\,L_\odot\) for A1 and \(7\pm 1\,L_\odot\) for A2 [2512.04674]. In **GSS 30**, FAUST multi-wavelength continuum analysis showed that the spectral index in IRS3 increases radially from \(2.0\) at the center to \(2.5\) at the disk edge while dropping to \(1.6\)–\(1.8\) along the outflow direction; SED fitting gave \(a_{\max}=30^{+12}_{-18}\,\mu{\rm m}\) and a dust mass of \(650\)–\(1510\,M_\oplus\). The same paper reported millimetre variability in IRS2, with the \(1.2\) mm flux decreasing from \(8.5\pm 1.1\) mJy to \(2.7\pm 1.2\) mJy over \(\simeq 420\) s, consistent with a magnetic flare [2606.22489].

| Paper | System | Key result |
|---|---|---|
| [2201.07334] | VLA 1623AB | \(12^\circ \pm 6^\circ\) disk-envelope misalignment |
| [2402.10258] | VLA 1623–2417 | \(\sim 2000\) au streamer and compact SiO jet |
| [2506.11363] | L1551 IRS 5 | \(r_{CB}=55\) au, \(j\sim270\) au km s\(^{-1}\) |
| [2512.10486] | L1551 IRS 5 | \(20\%\) dust over-density contrast from binary–disc interaction |
| [2512.04674] | NGC 1333 IRAS 4A | OCS-based luminosities for A1 and A2 |
| [2606.22489] | GSS 30 | Radial \(\alpha\) structure and rapid mm variability |

These studies show that FAUST treats chemistry and dynamics as a coupled problem. Misalignments, counter-rotation, streamer impact zones, cavity walls, and circumbinary pressure maxima are not secondary complications; they are the structures within which the observed chemistry is organized.

## 4. FAUST\(^2\): formal abstractions of uncountable-state stochastic processes

Outside astronomy, **FAUST\(^2\)** is a formal-methods software tool for generating abstractions of possibly non-deterministic **discrete-time Markov processes** defined over uncountable state spaces. A dtMP model \(\mathfrak S=(\mathcal S,\mathcal U,T_s)\) is specified in MATLAB and abstracted as a finite-state Markov chain or Markov decision process, with the abstraction procedure formally related to the concrete model through a user-defined maximum threshold on the approximation error. The toolbox can export abstract models to **PRISM** or **MRMC**, or can compute **PCTL** properties internally and refine the results back on the concrete dtMP using quantified error bounds [1403.3286].

Its architecture is divided into four components: a **Model Definition Module**, an **Abstraction Engine**, an **Error Analysis Module**, and an **Output & Verification Interface**. State and input spaces are given as box-shaped bounds in \(\mathbb R^n\) and \(\mathbb R^m\), and the stochastic kernel can be supplied by a linear-Gaussian template, a nonlinear-Gaussian template, or a user-defined density function. The GUI exposes “Formula-free,” “PCTL Safety,” and “PCTL Reach-Avoid” modes, as well as gridding assumptions and error controls [1403.3286].

The abstraction itself is partition based. FAUST\(^2\) chooses measurable cells \(\{A_i\}_{i=1}^p\), assigns representative points \(z_i\), and defines transition probabilities by integrating the concrete kernel over target cells:
\[
T_p(z_i\to z_{i'}) = T_s(A_{i'}\mid z_i)
= \int_{A_{i'}} t_s(s'\mid z_i)\,ds' .
\]
For controlled systems, the same construction yields an MDP kernel \(T_p(u_j,z_i\to z_{i'})\). The implementation relies on MATLAB’s vectorized `integral` routine, or `quad` in older releases, and can parallelize these computations if the Parallel Computing Toolbox is available [1403.3286].

A distinguishing feature of FAUST\(^2\) is explicit error quantification. Under a global Lipschitz assumption on the transition density, the horizon-\(N\) abstraction error satisfies
\[
E_N \le N\,h_s\,\delta_s\,\mathscr L(\mathcal S)\le \epsilon ,
\]
where \(\delta_s\) is the partition diameter and \(\epsilon\) is the user-supplied tolerance. For safety and reach-avoid properties, FAUST\(^2\) also supports a local Lipschitz, formula-dependent bound that drives adaptive sequential refinement of cells with large local errors [1403.3286].

The verification interface supports two workflows. One exports the abstract MC or MDP and labeling map to PRISM or MRMC. The other evaluates bounded-until PCTL formulas internally by dynamic programming on the abstract graph and then returns a guaranteed interval enclosure for the concrete probability. In the controlled case, the toolbox additionally solves a max or min over inputs to synthesize an \(\epsilon\)-optimal policy \(\pi^\*\) [1403.3286].

The reported case studies illustrate the intended operating regime. A **2D room-temperature control** example with safe set \([19.7,20.3]\times[4.7,5.3]\), horizon \(N=3\), and \(\epsilon=0.5\) produced about \(36\) states in under \(10\) s, with maximum one-step error about \(0.16\) and total safety error about \(0.48\). A **3D two-room extension** using adaptive gridding produced about \(150\) cells in under \(30\) s [1403.3286]. In this context, FAUST\(^2\) occupies a specific position between stochastic control, probabilistic verification, and numerical abstraction of continuous-state systems.

## 5. Faust in real-time audio DSP, fixed-point inference, and CLAP compilation

In audio and DSP research, **Faust** is a purely functional, block-diagram-oriented language for real-time audio and DSP. Faust programs are parameterized over the abstract type **FAUSTFLOAT**, which typically maps to IEEE-754 32-bit float in generated C++ or LLVM back ends. Recent work has focused on extending the compiler toward hardware-efficient fixed-point realizations and modern plugin targets [2403.06527].

The fixed-point format work augments the compiler’s symbolic propagation pass with interval analysis and precision inference. For each signal \(s\), the compiler computes an interval \([a_s,b_s]\) and chooses
\[
m_s=\left\lceil \log_2\bigl(\max(|a_s|,|b_s|)\bigr)\right\rceil+1
\]
to prevent overflow. The least-significant-bit position \(\ell_s\) is then chosen using a pseudo-injectivity criterion intended to preserve distinctness through quantized operators. In the generated code, signal types become `Ap_fixed<w_s,m_s>` with \(w_s=m_s-\ell_s+1\), and casts are inserted at core operations [2403.06527].

The numerical motivation is practical. The paper reports that on a mid-range FPGA a single-precision floating-point adder uses about **313 LUTs** and **11.4 ns** of combinational delay, whereas a **24-bit** two’s-complement integer adder uses **24 LUTs** and **1.7 ns**. Preliminary audio tests yielded SNR values of **32 dB** for a sine generator with step \(1/64\), **25 dB** for a sine generator with step \(0.01\), and **33 dB** for Karplus–Strong [2403.06527]. The work therefore treats fixed-point inference not as a manual back-end tweak but as a compiler-level semantic analysis problem.

A separate compiler-development line is **faust2clap**, described as the first officially maintained compilation pathway from Faust DSP specifications to the **CLAP** format. The framework has two modes. A **static** mode performs ahead-of-time compilation to native binaries through `faust2clap.py`, `clap-arch.cpp`, and a CMake build. A **dynamic** mode uses the Faust interpreter backend and a filesystem watcher to hot-reload DSP code without restarting the host application [2606.13193].

The major algorithmic issue in faust2clap is **parameter identity** under structural mutations of the DSP graph. The framework therefore uses an **address-based identity matching algorithm** to preserve parameter values across reloads and a **stable slot allocation** scheme to preserve host automation bindings. The implementation comprises approximately **2,400 lines of C++** architecture and plugin code plus approximately **200 lines of Python**, and has been integrated into the main Faust distribution [2606.13193].

The reported performance figures place the dynamic mode well within interactive use. Interpreter overhead per **256-sample block at 48 kHz** ranges from **0.008 ms** for a resonant low-pass filter to **0.27 ms** for `dm.zita` reverb, all below the **5.33 ms** real-time deadline. Reload latency scales with DSP size but remains **\(\le 60\) ms**, while compilation during hot reload typically takes **5–60 ms** [2606.13193]. In this literature, Faust is thus simultaneously a language, a compiler infrastructure, and a target-adaptation framework for embedded, FPGA, and plugin workflows.

## 6. Disambiguation, common patterns, and scholarly significance

The term **FAUST** is therefore genuinely polysemous across current technical literature. In astronomy it refers to an ALMA large program on the chemistry and dynamics of Solar-like protostars [2111.14121]. In formal methods it denotes a MATLAB toolbox for abstraction and verification of uncountable-state stochastic processes [1403.3286]. In audio DSP, stylistically often written **Faust**, it denotes a language and compiler ecosystem for real-time signal processing, with current research on automatic fixed-point formats and CLAP plugin generation [2403.06527] [2606.13193].

The domains are unrelated in subject matter, but a common structural pattern is visible. ALMA FAUST translates multi-line observations into a coherent view from envelope to disk and outflow scales; FAUST\(^2\) translates continuous-state stochastic dynamics into finite probabilistic models with quantified error; Faust compiler research translates high-level DSP graphs into fixed-point circuits or CLAP plugins. This suggests that the shared name recurs in contexts concerned with representation change across scales or abstractions, even though the scientific objects—protostars, dtMPs, and audio signal graphs—are entirely different.

For scholarly use, explicit disambiguation is essential. Citations such as **FAUST I**, **FAUST XVII**, or **FAUST XXXI** belong to the ALMA protostellar survey; **FAUST\(^2\)** belongs to stochastic-process abstraction and probabilistic model checking; and **Faust** in compiler papers denotes the audio DSP language. Treated separately, each of these literatures is technically coherent. Treated together under a single bare acronym, they require context, arXiv identifier, or subtitle to avoid category errors.

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