---
title: 'OffsetCrust: Variable-Radius Offset Surfaces'
url: https://www.emergentmind.com/topics/offsetcrust
type: topic
---

# OffsetCrust: Variable-Radius Offset Surfaces

Searching arXiv for the exact term to ground the article in the relevant paper and adjacent uses.
OffsetCrust is a framework for approximating **variable-radius offset surfaces** by recasting offset computation as a **power-diagram separation problem**. In its explicit, formal usage, the term refers to the geometry-processing method introduced in "OffsetCrust: Variable-Radius Offset Approximation with Power Diagrams" [2507.10924], where a base surface is sampled, paired with off-surface weighted sites, and the desired offset is extracted from power-diagram facets separating the two site classes. In a distinct and nonstandard secondary usage, closely related neutron-star papers use “offset crust” only as an interpretive description for crust-induced shifts, matching effects, or low-density systematic biases, rather than as a named framework [2003.03330], [1902.04616], [2406.14906].

## 1. Formal definition in geometry processing

In the geometry-processing literature, OffsetCrust addresses the problem of computing **variable-radius offset surfaces**. Given a base surface \(\mathcal S\) and a positive radius function \(\mathcal R:\mathcal S\to \mathbb R_{>0}\), the offset is formulated through the generalized distance field
\[
\phi(x)\coloneqq \min_{p\in\mathcal S}\big(\|x-p\|-\mathcal R(p)\big),
\]
with the variable-radius offset surface defined as the zero level set
\[
\mathcal S_{\mathcal R}^{\text{off}=\{x\mid \phi(x)=0\}.
\]
For constant radius \(r\), this reduces to the usual Minkowski-sum view \(\mathcal S \oplus B_r\); for variable radius, each point \(p\in\mathcal S\) contributes its own ball \(B_{\mathcal R(p)}(p)\), and the offset is the envelope of these balls [2507.10924].

The radius field is defined directly on the base surface,
\[
\mathcal R:\mathcal S\to \mathbb R_{>0},
\]
and for mesh input it is typically given at vertices and assumed linear over each triangle. If sparse values are given, missing values are interpolated by solving the biharmonic equation
\[
\Delta^2 \mathcal R = 0.
\]
A key assumption used by the method is
\[
\|\nabla \mathcal R(p)\|\le 1,
\]
which ensures that the displacement-direction construction is valid and is also consistent with medial axis transforms, where balls are maximal and non-nested [2507.10924].

This formulation distinguishes OffsetCrust from constant-radius offsetting methods. The central difficulty is that, in the variable-radius setting, the envelope geometry is no longer captured by a uniform normal displacement. OffsetCrust therefore replaces the standard normal-offset construction with a weighted arrangement in which the target surface appears as a bisector structure in power distance [2507.10924].

## 2. Power-diagram construction and geometric principle

OffsetCrust uses a **power diagram** of weighted sites. For weighted sites \(\{(p_i,w_i)\}_{i=1}^n\), the power distance is
\[
d^{\mathrm{pow}_{p_i}(x)=\|x-p_i\|^2-w_i.
\]
The method constructs two classes of sites: **base sites** on \(\mathcal S\), weighted by \(\mathcal R^2(p)\), and **displaced sites** off the surface, weighted by \((\mathcal R(p)-\epsilon)^2\). The offset surface is then approximated by the power-diagram facets separating these two groups [2507.10924].

The continuous geometric argument is based on a contributing point \(p\in\mathcal S\) and its offset point \(p_{\mathcal R^{\text{off}\). For small \(\epsilon>0\), the paper defines
\[
p_{\mathcal R^\epsilon} := p + \epsilon\cdot \boldsymbol n_{\mathcal R^{\text{off}(p)},
\]
with \(p\), \(p_{\mathcal R^\epsilon}\), and \(p_{\mathcal R^{\text{off}\) colinear. The offset point satisfies
\[
\|p_{\mathcal R^{\text{off}-p\|-\mathcal R(p) = \|p_{\mathcal R^{\text{off}-p_{\mathcal R^\epsilon}\|-(\mathcal R(p)-\epsilon) =0,
\]
and in squared form
\[
\|p_{\mathcal R^{\text{off}-p\|^2-\mathcal R^2(p) = \|p_{\mathcal R^{\text{off}-p_{\mathcal R^\epsilon}\|^2-(\mathcal R(p)-\epsilon)^2 =0.
\]
This is the reason the envelope can be represented through weighted bisectors, and hence through a power-diagram construction [2507.10924].

The site construction follows this derivation directly. For each contributing base point \(p\), OffsetCrust creates the base site \((p,\mathcal R^2(p))\) and the displaced site \((p_{\mathcal R^\epsilon},(\mathcal R(p)-\epsilon)^2)\). The extracted surface is the subset of the power diagram that separates base from displaced sites [2507.10924].

The framework also explains why non-contributing points do not affect the envelope. If an offset point \(p_{\mathcal R^{\text{off}\) lies inside the envelope, then some other site \(q\) satisfies
\[
\|q-p_{\mathcal R^{\text{off}\|<\mathcal R(q),
\]
equivalently
\[
\|q-p_{\mathcal R^{\text{off}\|^2-\mathcal R^2(q)<0.
\]
Since \(\|p-p_{\mathcal R^{\text{off}\|^2-\mathcal R^2(p)=0\), it follows that \(p_{\mathcal R^{\text{off}\) belongs to another power cell. This means only contributors to the envelope matter in the final separation structure [2507.10924].

## 3. Displacement directions and the variable-radius deviation from normal offsets

A defining feature of OffsetCrust is that, in the variable-radius case, the displacement direction is generally **not** the surface normal. For constant \(\mathcal R\), the offset direction aligns with the surface normal; this is the classical normal-offset regime. For spatially varying radii, the paper proves that if \(\|\nabla \mathcal R(p)\|\le 1\), the offset direction \(\boldsymbol n_{\mathcal R^{\text{off}(p)\) is obtained by rotating the unit surface normal \(\boldsymbol n_p\) by the angle
\[
\alpha=\arcsin(\|\nabla \mathcal R(p)\|),
\]
around the axis
\[
\frac{\nabla \mathcal R(p)}{\|\nabla \mathcal R(p)\|}\times \boldsymbol n_p.
\]
The geometric interpretation given is that spatial variation in \(\mathcal R\) tilts the tangent of the envelope away from the base normal; only when \(\nabla \mathcal R=0\) do the directions coincide with the normals [2507.10924].

This displacement-direction theorem is the principal reason a direct transplant of constant-radius crust methods is insufficient. The method must account for a direction field determined jointly by the base geometry and the radius gradient. In the language of the paper, the off-surface points are displaced along **\(R\)-dependent directions**, and only in the constant-radius case do these directions align exactly with the surface normals of \(S\) [2507.10924].

A practical implication is that variable-radius offsetting is not merely a weighted version of standard offsetting. The site construction, power weights, and later refinement stages are all designed to manage the mismatch between base normals and true envelope directions. This suggests that the principal technical novelty of OffsetCrust lies not only in adopting power diagrams, but in making them compatible with non-normal displacement geometry [2507.10924].

## 4. Misalignment, fine-tuning, and sampling pipeline

OffsetCrust inherits a known weakness of crust-based approaches: **misaligned facets**. A facet in the separating structure may be formed either by a base point \(p_i\) and its own displaced point or by a base point \(p_i\) and the displaced point of another site. These “cross-pair” facets generate misalignment, and the paper notes that the problem is more likely in the variable-radius setting because displacement directions differ from the normal and vary spatially [2507.10924].

To reduce this artifact, the method uses a **local least-squares vertex refinement**. If \(\{p_i\}_{i=1}^k\) are the base points contributing to a power-diagram vertex \(v\), with
\[
d_i:=\mathcal R(p_i),\qquad \boldsymbol n_i^{\text{off}:=\boldsymbol n_{\mathcal R^{\text{off}(p_i)},
\]
then the ideal constraint is
\[
(v-p_i)\cdot \boldsymbol n_i^{\text{off}-d_i=0.
\]
OffsetCrust therefore minimizes
\[
\min_v \sum_{i=1}^k \left((v-p_i)\cdot \boldsymbol n_i^{\text{off}-d_i\right)^2,
\]
and stabilizes it with a small regularizer,
\[
\min_v \sum_{i=1}^k \left((v-p_i)\cdot \boldsymbol n_i^{\text{off}-d_i\right)^2 +\lambda\|v-v_0\|^2,
\]
where \(v_0\) is the original vertex and \(\lambda\) is small, typically \(10^{-4}\). The closed-form solution is
\[
v^*=\boldsymbol H^{-1}(\lambda v_0+\boldsymbol b),
\]
with
\[
\boldsymbol H=\sum_{i=1}^k \boldsymbol n_i^{\text{off}(\boldsymbol n_i^{\text{off})^T+\lambda I,
\]
\[
\boldsymbol b=\sum_{i=1}^k (p_i\cdot \boldsymbol n_i^{\text{off}+d_i)\boldsymbol n_i^{\text{off}.
\]
The paper characterizes this step as “lightweight” because it is per-vertex, local, quadratic, and solved by a small linear system; no global optimization or explicit self-intersection repair is needed [2507.10924].

The sampling pipeline is correspondingly specialized. Triangle-interior points are sampled by blue noise, each generating one displaced point. Edge-type points and vertex-type points use a **one base point, multiple displaced points (1vN)** strategy. For a manifold edge shared by faces \(f_1,f_2\), multiple directions are generated by spherical linear interpolation between face-based directions. Around each vertex \(v\), the method places a small sphere of radius \(\rho l\), intersects it with the surface, samples the resulting closed trajectory, and again uses spherical linear interpolation to generate displaced points. Interior triangle points inside this protected sphere are removed so that vertex-generated displaced points have priority. An optional dihedral-angle threshold can pre-detect sharp feature lines so that 1vN is only applied where needed [2507.10924].

These design choices indicate that OffsetCrust is not only a theoretical reformulation. It is also a sampling and reconstruction pipeline adapted to sharp features, anisotropic offset directions, and local ambiguity in the crust-style separating structure [2507.10924].

## 5. Implementation, evaluation, and application domain

The implementation reported for OffsetCrust is in C++, using **CGAL** exact predicates/exact constructions, **TBB** for parallelization, **Eigen LDLT** for the refinement solve, and **AABB trees**, **PQP**, and **libigl** for distance and inside/outside queries. Experiments were run on a machine with Intel i9-13900K and 64 GB RAM [2507.10924].

Typical settings reported are **70K blue-noise samples**, \(\lambda=0.01\), \(\epsilon=10^{-6}\), \(\rho=5\%\), and a discrete spherical surface with **642 vertices**. For evaluation, the paper uses **CD** (Chamfer Distance), **HD** (Hausdorff Distance), and **NC** (Normal Consistency). For constant-radius accuracy it also uses a one-sided normalized distance metric,
\[
|\mathbf D(x)-d|,\qquad d=\delta\cdot l_{\text{diag}.
\]
The reported average runtime is about **100 seconds** with roughly \(10^6\) samples. Power-diagram computation is the main bottleneck, facet and adjacency extraction is also costly, and the refinement step is relatively cheap [2507.10924].

The empirical comparison is made against dual contouring at \(300^3\). OffsetCrust is reported to achieve comparable or better CD/HD stability, fewer outliers, and strongly improved NC after refinement. For inward offsets it often performs better than dual contouring in NC, while outward offsets are competitive though not always superior. The paper also reports good performance on **2K Thingi10K models**, strong robustness for \(\delta=\pm 2\%\), and a roughly order-of-magnitude improvement in normal consistency after refinement in many cases. Ablation studies indicate that more blue-noise samples improve accuracy but increase runtime, smaller \(\rho\) helps small offsets preserve sharpness, finer spherical interpolation improves rounded regions and accuracy, and overly aggressive sharp-feature filtering hurts quality [2507.10924].

A prominent application is **medial axis transform (MAT) reconstruction**. In this setting, the MAT surface is the base surface, radii are already defined at MAT elements, and the recovery of the original boundary becomes a variable-radius offset problem. The method is particularly natural here because MAT balls satisfy the condition
\[
\|\nabla \mathcal R\|<1.
\]
The paper reports faithful reconstructions of original surfaces with good CD, HD, and NC, while preserving geometry and topology without requiring explicit envelope handling [2507.10924].

## 6. Secondary usages in neutron-star research

Outside geometry processing, “offset crust” is not a standardized term but appears as an interpretive label for several distinct crust-related effects in neutron-star research. These usages should be distinguished from the formal framework named OffsetCrust [2003.03330], [1902.04616], [2406.14906].

One usage concerns a **radius posterior shift** caused by the low-density crust equation of state. In "The impact of the crust equation of state on the analysis of GW170817" [1902.04616], the central result is that different crust models do not strongly impact the mass or tidal deformability of a neutron star, but they do affect the inferred radius. Earlier GW170817 analyses fixed densities below \(\rho \approx 10^{14}\,\mathrm{g/cm^3}\) to the SLy description, yielding
\[
R_1 = 11.9_{-1.4}^{+1.4}\,\mathrm{km}, \qquad R_2 = 11.9_{-1.4}^{+1.4}\,\mathrm{km}.
\]
Reanalysis with a changed crust gave
\[
R_1 = 11.7_{-1.4}^{+1.4}\,\mathrm{km}, \qquad R_2 = 11.7_{-1.4}^{+1.3}\,\mathrm{km},
\]
and the paper estimates a crust-induced systematic radius shift of about \(\sim 0.3~\mathrm{km}\), approximately \(\sim 3\%\) of the neutron-star radius. In that context, “offset crust” refers to the fact that the credible region in \(\Lambda\) is nearly unchanged while the radius posterior is shifted [1902.04616].

A related but distinct usage concerns **non-unified crust-core matching** and inference bias. "Inference of neutron-star properties with unified crust-core equations of state for parameter estimation" [2406.14906] develops **CUTER** to consistently match a nuclear-physics-informed crust to an arbitrary high-density EoS. The paper argues that a fixed, realistic-but-inconsistent crust causes small but avoidable errors in the estimation of global neutron-star properties and leads to an underestimation of uncertainties. Quantitatively, for tested models the relative error in radius after CUTER reconstruction is \(\lesssim 0.5\%\), while fixed-crust treatments can slightly shift medians and shrink posterior widths. For example, under LD+HD filters only with \(n_{\rm match}=0.16\) fm\(^{-3}\), unified treatment gives
\[
R_{1.4}=13.57^{+0.54}_{-0.81}\ \mathrm{km},\qquad \Lambda_{1.4}=925^{+274}_{-332},
\]
whereas a unique SLy4 crust gives
\[
R_{1.4}=13.83^{+0.38}_{-0.75}\ \mathrm{km},\qquad \Lambda_{1.4}=1025^{+226}_{-349}.
\]
Here the “offset” is an inference-level bias induced by inconsistent low-density modeling [2406.14906].

A third usage is found in "Probing Crust Meltdown in Inspiraling Binary Neutron Stars" [2003.03330], where the relevant effect is a **crust-meltdown-induced phase offset** in the gravitational-wave signal. There the inspiraling companion excites crust-core interface modes, the crust yields plastically once the local elastic strain exceeds a breaking strain of order
\[
\epsilon_b\sim 0.1,
\]
and dissipative heating eventually melts the crust. The resulting waveform modification is written as
\[
\delta\Psi(f)=\sum_{i=1,2}\delta\phi_i\left(1-\frac{f}{f_i}\right)\Theta(f-f_i) \approx \delta\phi_a\left(1-\frac{f}{f_a}\right)\Theta(f-f_a),
\]
with the paper quoting a \(\sim\mathcal O(0.1)\) phase shift and noting that, for an equal-mass \(1.3M_\odot+1.3M_\odot\) binary, \(\delta\phi_a\) varies from about \(0.03\) to \(0.6\) depending on the EoS and \(n_{\rm b,cc}\). In that setting, “offset crust” denotes a localized, resonance-triggered phase and time jump in the inspiral waveform, not a surface-reconstruction method [2003.03330].

A further neutron-star usage concerns **crust-core interpolation**. "Consistent crust-core interpolation and its effect on non-radial neutron star oscillations" [2502.02373] studies a thermodynamically and causally consistent interpolation in the pressure–chemical-potential plane, contrasting it with classical concatenation. The interpolation satisfies
\[
P(\mu)=P_1+\Delta P\,\frac{\mu^{b}-\mu_1^{b}}{\mu_2^{b}-\mu_1^{b}}, \qquad
P=c_s^2(\varepsilon-\tilde{\varepsilon}),
\]
introduces an energy-density gap corresponding to a first-order phase transition, and affects low-mass-star radii, crust thickness, and especially the \(p_1\)-mode frequency, while leaving the \(f\)-mode essentially unchanged. Here again the phrase is descriptive rather than terminologically fixed [2502.02373].

## 7. Conceptual scope and distinctions

The term OffsetCrust therefore has a sharply defined meaning in one domain and only a loose descriptive role in another. In geometry processing, it denotes a named framework whose central ingredients are a variable-radius offset formulation, weighted base and displaced sites, power-diagram extraction, and lightweight local refinement [2507.10924]. In neutron-star studies, by contrast, “offset crust” has been used only to describe crust-induced shifts in radii, waveform phase, or crust-core matching, and the underlying phenomena are physically unrelated to the geometric method [1902.04616], [2003.03330], [2406.14906].

This distinction matters because the two literatures use similar words for entirely different technical objects. In geometry processing, the offset is a surface envelope generated by balls of spatially varying radius, and the “crust” refers to a class of reconstruction methods inspired by crust-based surface extraction [2507.10924]. In neutron-star astrophysics, the offset is a shift in an inferred or observable quantity produced by crust microphysics or EOS treatment, while “crust” refers to the star’s low-density solid outer layers [1902.04616], [2003.03330].

A plausible implication is that the geometry-processing usage is likely to dominate the exact proper noun **OffsetCrust**, because it appears as the title of a dedicated framework [2507.10924]. The neutron-star usages remain valuable for query interpretation, but they are better understood as context-dependent descriptive phrases rather than as a unified concept.

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