---
title: Vertical-Slice Growth Protocol
url: https://www.emergentmind.com/topics/vertical-slice-growth-protocol
type: topic
---

# Vertical-Slice Growth Protocol

Vertical-Slice Growth Protocol denotes, across several technical literatures, a class of end-to-end procedures organized around a vertically resolved slice, a vertically stacked layer system, or a per-vertical service bundle. The surveyed usages indicate that the label does not name a single canonical method; rather, it refers to domain-specific workflows that connect formal models to operational decisions in indoor crop control, chemical vapor deposition of van der Waals materials, atmospheric slice dynamics, 5G network slicing for vertical industries, and height-aware 3D semantic occupancy prediction [2309.07540], [1911.11135], [1211.2067], [2104.04743], [2509.03999].

## 1. Terminological scope and core meanings

The collected literature exhibits two distinct senses of *vertical*. In indoor farming, van der Waals growth, atmospheric slice modeling, and semantic occupancy, the term refers to a geometric height or stacking direction. In 5G network slicing, it refers to a *vertical industry* such as automobile, manufacturing, or power grid. In all cases, the protocol is an implementable synthesis that begins from a structured model and ends in staged execution, monitoring, and adaptation.

| Domain | Vertical-slice object | Primary objective |
|---|---|---|
| Indoor vertical farming | Daily controls and annual crop schedule | Resource–yield–duration trade-off |
| van der Waals materials | Stacked layers and step boundaries | Vertical versus lateral growth control |
| Atmospheric dynamics | \(x\)–\(z\) slice with three velocity components | Frontogenesis and baroclinic-instability analysis |
| 5G network slicing | Per-vertical NSIs, US-NSIs, GN-NSIs, S-NSIs | Provisioning and life-cycle management |
| 3D semantic occupancy | Height-axis slices of voxel features | Height-aware multimodal representation |

A plausible synthesis is that a Vertical-Slice Growth Protocol is characterized by five recurrent elements: a reduced but explicit state representation; a small set of control or design variables; stage-wise execution; constraints tied to physics, service guarantees, or geometry; and a translation layer that maps the formal model to operational practice.

## 2. Indoor vertical farms: optimal-control scheduling for wheat

In "Optimal Control for Indoor Vertical Farms Based on Crop Growth" [2309.07540], the protocol operationalizes an optimal-control formulation for indoor wheat cultivation, using cultivar Batten as the reference crop. The daily control vector is
\[
u_i = [T_i, D_i, R_i]^\mathsf{T},
\]
where \(T_i\) is the mean air/leaf temperature setpoint, \(D_i\) is the drought or soil-moisture stress index ARID in \([0,1]\), and \(R_i\) is artificial plant-available radiation expressed as daily SRAD in MJ m\(^{-2}\) d\(^{-1}\). The state vector is
\[
x_i = [m_{B,i}, \tau_i, I_{50B,i}]^\mathsf{T},
\]
with biomass, cumulative temperature, and cumulative temperature to accelerate senescence. The paper fixes \(\mathrm{CO}_2\) enrichment at \(700\) ppm rather than optimizing it.

The crop dynamics are cast as a nonlinear, discrete-time hybrid model derived from SIMPLE, with a smooth approximation SIMPLE-SC used in optimization:
\[
x_{i+1} = x_i + f_s(x_i,u_i).
\]
Non-smooth \(\max\) operators are replaced by
\[
\operatorname{smax}_\epsilon(a,b)=\frac{a+b+\sqrt{(a-b)^2+\epsilon}}{2},
\]
with \(\epsilon = 10^{-4}\). Maturity is defined phenologically through the canopy radiation interception fraction \(f_\text{Solar}(\tau, I_{50B})\), and harvest is enforced by the terminal equality
\[
g(x_N)=f_\text{Solar}(\tau_N, I_{50B,N})-0.005=0.
\]

Two optimal-control problems are posed. For fixed final time,
\[
J(U,x_0)=\sum_{i=0}^{N-1} l(u_i)-V(x_N),
\]
where \(l(u_i)\) encodes energy and water costs and \(V(x_N)\) encodes crop revenue. For free final time, the dynamics become
\[
x_{i+1}=x_i+T f_s(x_i,u_i), \qquad T\in[0.5,1.5],
\]
and the annualized objective is
\[
J_t(U_T,x_0)=\sum_{i=0}^{N-1}\frac{365}{N}l(u_i)-\frac{365}{NT}V(x_N).
\]
The stage-cost weights used are \(C_\vartheta=1.8\times 10^{-6}\;€/\!^\circ\mathrm{C}^2\), \(C_D=0.02\;€\), \(C_R=0.038\;€/(\mathrm{MJ}\,\mathrm{m}^{-2}\mathrm{d}^{-1})\), \(\vartheta_0=10^\circ\mathrm{C}\), \(HI=0.3\), and \(c_\text{crop}=132.9\;€/\mathrm{kg}\).

The optimization is implemented in Python with automatic differentiation through CasADi; standard gradient-based NLP solvers such as IPOPT can be used within CasADi. The free-final-time routine initializes \(N=N_0\), solves the OCP, updates \(T\leftarrow T^*\), sets \(N\leftarrow \lfloor TN\rfloor\), and repeats until \(|T-1|<\delta\) with \(\delta=0.01\). The protocol uses a daily time step and modest computational resources, which the paper states are suitable for iterative annual scheduling updates or periodic re-optimization.

The reported optimum converges to \(N^*=102\) days with \(T\approx 1\) after nine iterations, corresponding to about \(365/102\approx 3.58\) cycles per year. The stage-wise daily schedule has three phases. During germination and early canopy closure, approximately days \(0\)–\(12\), temperature is near the upper bound at about \(35^\circ\mathrm{C}\), drought is near \(0\%\), and radiation rises to \(35\) MJ m\(^{-2}\) d\(^{-1}\) by about day \(13\). During rapid biomass accumulation, approximately days \(13\)–\(88\), temperature is held near \(T_\text{opt}=15^\circ\mathrm{C}\), drought remains at \(0\%\), and radiation remains at the upper bound \(35\) MJ m\(^{-2}\) d\(^{-1}\). During maturation and drying, approximately days \(89\)–\(102\), temperature rises gradually to about \(29^\circ\mathrm{C}\), drought ramps to about \(100\%\), and radiation decreases sharply toward zero.

The resulting biomass per cycle is \(m_B\approx 3.081\) kg m\(^{-2}\), equal to the constant-input case \([T=23^\circ\mathrm{C}, D=0, R=35\ \mathrm{MJ}]\) but achieved in fewer days. Annual biomass becomes about \(11.02\) kg m\(^{-2}\) y\(^{-1}\), compared with about \(10.22\) kg m\(^{-2}\) y\(^{-1}\) for constant inputs, an increase of about \(8\%\). The optimized cost per cycle is about \(104.55\;€/\mathrm{kg}\), versus \(149.21\;€/\mathrm{kg}\) for constant inputs, a reduction of about \(30\%\), and annual energy cost is reduced by about \(25\%\) while biomass increases. A hypothetical zero-energy-cost case yields \(N=119\) days and \(m_B\approx 4.09\) kg m\(^{-2}\), more than \(30\%\) above the cost-constrained optimum, illustrating the cost–yield tension.

The implementation layer is explicit. Because the model optimizes SRAD rather than photoperiod, PPFD, or spectrum, the protocol provides a practical conversion through the daily light integral:
\[
\mathrm{DLI}\approx \frac{R}{0.218}.
\]
At \(R=35\) MJ m\(^{-2}\) d\(^{-1}\), this implies \(\mathrm{DLI}\approx 160.6\) mol d\(^{-1}\) and, for a 16 h photoperiod, an average PPFD of about \(2785\ \mu\mathrm{mol}\ \mathrm{m}^{-2}\ \mathrm{s}^{-1}\), which the protocol identifies as very high. It therefore recommends, in practice, capping PPFD and extending photoperiod where necessary. The paper also states strong limitations: nutrients, vernalization, photoperiod, PPFD caps, spectrum, and explicit HVAC energy are omitted; the SIMPLE parameterization is field-calibrated rather than VF-calibrated; and the notional wheat price is set at \(400\times\) market to illustrate profitability, while realistic wheat pricing leaves indoor wheat currently unprofitable.

## 3. van der Waals stacked 2D materials: thermodynamic and kinetic control of vertical growth

In "Modeling the vertical growth of van der Waals stacked 2D materials using the diffuse domain method" [1911.11135], the protocol concerns CVD growth of vertically stacked monolayers such as graphene and transition-metal dichalcogenides. The physical picture is a Burton–Cabrera–Frank-type model in which adatoms are deposited with flux \(F\), diffuse with coefficient \(D\), may desorb at rate \(\tau_d^{-1}\), and attach to steps with kinetic rates \(k^\pm\). Local thermodynamics are governed by binding energies, edge energies, and curvature through Gibbs–Thomson terms. Lateral growth adds atoms to the perimeter of the first layer \(\Omega_1\), whereas vertical stacking adds atoms to the second layer \(\Omega_2\) above \(\Omega_1\), incurring a new perimeter penalty.

The energetic bias between the two pathways is expressed by
\[
\Delta \mathcal{E} = \mathcal{E}_2 - 2\mathcal{E}_1.
\]
The full free energy includes binding terms, edge terms, and regular-solution entropy. Adatom concentrations satisfy
\[
\partial_t \rho_i = D_i \Delta \rho_i + F_i - \tau_{d,i}^{-1}\rho_i
\quad \text{in }\Omega_i,\quad i=0,1,2,
\]
and the step velocities obey
\[
v_i = q_i^+ + q_i^- + \beta\, \partial_s^2 \kappa_i.
\]
Boundary conditions at \(\Gamma_1\) and \(\Gamma_2\) couple adjacent layers through attachment kinetics and local equilibrium densities.

For quasi-steady, isotropic, circular layers, the paper derives an analytic criterion for vertical growth:
\[
\Delta \mathcal{E} > \frac{\gamma_2}{R_2} - \frac{\gamma_1}{R_1}.
\]
This states that the interlayer binding advantage must exceed the edge-energy penalty per atom associated with creating the upper-layer perimeter. The corresponding critical radii are
\[
R_{2,c}=\frac{\gamma_2}{\Delta\mathcal{E}}, \qquad
R_{1,c}=\frac{\gamma_1}{\frac{\gamma_2}{R_2}-\Delta\mathcal{E}}.
\]
A kinetic critical radius \(R_{2,k}\) is also identified: when \(R_{2,c}>R_2(0)>R_{2,k}\), transient second-layer growth can occur before the first layer reaches the suppressing threshold \(R_{1,c}\). Increasing \(D_1\), \(k_2^-\), or \(k_1^+\) and decreasing \(F_0\) enlarge this transient window.

To handle arbitrary layer shapes, the model is reformulated by the diffuse domain method. The phase variable \(\varphi\) marks layers with \(\varphi\approx 0\) in \(\Omega_0\), \(\varphi\approx 1\) in \(\Omega_1\), and \(\varphi\approx 2\) in \(\Omega_2\). Adatom fields are solved on a fixed rectangular domain \(\tilde{\Omega}\), and interface motion is governed by a Cahn–Hilliard-type equation with kinetic source terms. The numerical method is mass-conservative, semi-implicit, and second-order accurate in time and space, using Crank–Nicholson time stepping, centered finite differences, a nonlinear multigrid FAS solver, and block-structured adaptive mesh refinement with four levels and finest mesh \(h_3=h_0/8\). Interface thickness \(\epsilon\) is chosen proportional to mesh size, for example \(h_3=\epsilon/6.4\), yielding overall \(O(\epsilon^2)\) accuracy in \(L^2\) and \(L^\infty\) norms.

Parameter sensitivity is central to the protocol. Larger \(\Delta\mathcal{E}\) increases vertical growth. Increasing \(\gamma\) reduces \(v_2\) and can suppress \(\Omega_2\) growth, with simulations noting suppression around \(\gamma\approx 0.16\). Larger \(\beta\) smooths corners and reduces growth rates. Decreasing \(F_0\) favors vertical growth kinetically, while increasing \(k_2^-\) generally feeds \(\Omega_2\). Kinetic and edge anisotropy are modeled through
\[
\xi_k(\theta)=1-\epsilon_{k,n}\cos\big(n(\theta-\theta_0)\big),
\]
and
\[
\gamma(\theta)=\gamma\,\big(\xi_s(\theta)+\xi_s''(\theta)\big),\qquad
\xi_s(\theta)=1-\epsilon_{s,n}\cos(n\theta),
\]
with \(n=6\) for graphene-like and \(n=3\) for MoS\(_2\)-like systems. The paper reports that, in 6-fold anisotropy, increasing \(k_1^+\) can trigger \(\Omega_1\) morphological instability, while in 3-fold anisotropy, increasing \(k_1^+\) favors \(\Omega_2\) growth without instability.

The protocol translates directly into CVD guidance: estimate \(R_1\), \(R_2\), \(\Delta\mathcal{E}\), \(\gamma_1\), and \(\gamma_2\); evaluate the growth criterion; and, if it fails, reduce substrate flux \(F_0\), increase temperature to raise \(D_1\) and \(k_2^-\), modify edge passivation to lower \(\gamma_2\) or raise \(\gamma_1\), or seed a larger \(R_2\). The model reproduces star-shaped bilayer graphene, twisted bilayer graphene, vertically stacked MoS\(_2\) triangles with negative curvature, and nearly overlapped bilayers. Its principal limitations are the neglect of strain and lattice mismatch, frequent omission of desorption, neglect of anisotropy in \(\beta\), and a primary focus on two-layer systems despite an \(n\)-layer extension.

## 4. Atmospheric science: variational vertical-slice models and growth experiments

In "A variational formulation of vertical slice models" [1211.2067], a vertical slice is a three-dimensional velocity field depending only on \(x\), \(z\), and \(t\). The slice velocity is \(u_S=(u,w)\) in the \(x\)–\(z\) plane and the transverse velocity is \(u_T=v\) in the \(y\)-direction. The specialized Lagrangian map is
\[
\phi(X,Y,Z,t)=\big(x(X,Z,t),\, y(X,Z,t)+Y,\, z(X,Z,t)\big),
\]
and the configuration space is the semidirect product \(\mathrm{Diff}(\Omega)\circledS \mathcal{F}(\Omega)\). Areal density \(D\) is advected within the slice, while potential temperature is decomposed as
\[
\theta(x,y,z,t)=\theta_S(x,z,t)+(y-y_0)s,
\]
which yields the slice evolution law
\[
\partial_t \theta_S + u_S\cdot \nabla \theta_S + u_T s = 0.
\]

For a slice Lagrangian \(l[(u_S,u_T),(\theta_S,s),D]\), the Euler–Poincaré equations lead to a Hamiltonian structure, a Lie–Poisson bracket, and conservation laws. In geometric variables,
\[
v_S=\frac{1}{D}\frac{\delta l}{\delta u_S}, \qquad
v_T=\frac{1}{D}\frac{\delta l}{\delta u_T},
\]
the Kelvin–Noether circulation theorem is
\[
\frac{d}{dt}\oint_{c(u_S)} \left(sv_S - v_T\nabla\theta_S\right)\cdot dx = 0,
\]
and the corresponding slice potential vorticity is materially conserved:
\[
\partial_t q + u_S\cdot\nabla q = 0,\qquad
q=\frac{1}{D}\left(s\,\mathrm{curl}\,v_S + \nabla \theta_S \times \nabla v_T\right)\cdot \hat y.
\]

The incompressible Eady–Boussinesq slice model is obtained from the Lagrangian
\[
l = \int_\Omega \frac{D}{2}\left(|u_S|^2+u_T^2\right) + D f u_T x
+\frac{g}{\theta_0}D\left(z-\frac{H}{2}\right)\theta_S
+ p(1-D)\, dV,
\]
which, under \(D=1\) and \(\nabla\cdot u_S=0\), yields
\[
\begin{aligned}
\partial_t u_S + u_S\cdot\nabla u_S - f u_T \hat x &= -\nabla p + \frac{g}{\theta_0}\theta_S \hat z,\\
\partial_t u_T + u_S\cdot\nabla u_T + f u_S\cdot \hat x &= -\frac{g}{\theta_0}\left(z-\frac{H}{2}\right)s,\\
\nabla\cdot u_S &= 0,\\
\partial_t \theta_S + u_S\cdot\nabla \theta_S + u_T s &= 0.
\end{aligned}
\]
The conserved energy is
\[
h=\int_\Omega \frac{1}{2}\left(|u_S|^2+u_T^2\right)
-\frac{g}{\theta_0}\left(z-\frac{H}{2}\right)\theta_S\, dV.
\]

The paper also introduces a sliced compressible model with Lagrangian
\[
l=\int_\Omega \frac{D}{2}\left(|u_S|^2+u_T^2\right)+fDu_Tx+gDz-Dc_v\theta_S\Pi\, dV,
\]
which preserves its own energy, Kelvin–Noether circulation, and potential vorticity. It is explicitly stated that this compressible model does not produce solutions that are also solutions of the three-dimensional compressible equations, but it reduces to the Eady–Boussinesq model in the low-Mach-number limit, making it useful for asymptotic limit error testing.

The protocol built on this framework prescribes numerical frontogenesis and baroclinic-instability experiments in a slice domain \(\Omega=[0,L_x]\times[0,H]\), typically periodic in \(x\) and rigid free-slip at top and bottom in \(z\). The summary specifies a uniform \(N_x\times N_z\) grid, for example \(256\times 128\) for linear stages, and a time step chosen by stability or CFL constraints. The base Eady-type incompressible state uses
\[
U(z) = -\frac{g s}{f\theta_0}\left(z-\frac{H}{2}\right), \qquad
u_S^0=(U(z),0), \qquad u_T^0=0,
\]
with small perturbations added through \(\theta_S\), a streamfunction-generated incompressible perturbation \(u_S'=\nabla^\perp \psi\), and a compatible \(u_T'\).

Growth is diagnosed numerically rather than through an explicit linear theory supplied by the paper. Perturbation kinetic energy,
\[
E'_K(t)=\frac12\int_\Omega \left(|u_S'|^2+|u_T'|^2\right)\, dV,
\]
supports an estimated exponential growth rate
\[
\sigma \approx \frac{d}{dt}\left(\frac12\ln E'_K(t)\right).
\]
Frontogenesis is tracked through \(\|\nabla\theta_S\|_\infty\), \(\|\nabla\theta_S\|_2\), slice PV, Kelvin–Noether circulation integrals, enstrophy, and energy conservation. The protocol places particular emphasis on structure-preserving discretization so that PV and circulation remain numerically compatible with the geometric form of the continuous model.

## 5. 5G network slicing: per-vertical, per-use-case life-cycle growth

In "How Should Network Slice Instances be Provided to Multiple Use Cases of a Single Vertical Industry?" [2104.04743], a Vertical-Slice Growth Protocol refers to the provisioning and evolution of per-vertical, per-use-case network slices in a 5G system. A vertical industry consists of multiple use cases with differing service, network, and connectivity requirements. The central objects are the end-to-end Network Slice Instance (NSI), the Use-case Specific NSI (US-NSI), and the Generic NSI (GN-NSI) built from Sub-network Slice Instances (S-NSIs).

A US-NSI is one end-to-end NSI per use case. The number of US-NSIs equals the number of use cases, each with its own SLA, templates, life-cycle, and isolation. Sub Network Slicing introduces S-NSIs dedicated to single use cases; a GN-NSI is a per-vertical cluster of S-NSIs offered as a single business product. GN-NSIs are either standard, with predefined templates registered in the network slice catalog, or non-standard, co-designed by operator and tenant and then inserted into the catalog. The protocol distinguishes the two models operationally: US-NSI maximizes isolation and independent life-cycle control, whereas GN-NSI supports modular reusability, coordinated governance, and selected sharing of common functions.

The architecture combines 3GPP and ETSI components. Tenant-facing exposure occurs through SCEF/NEF northbound APIs over T8, using SCS/AS; GN-NSI introduces a Master-SCS/AS coordinating multiple Slave-SCS/ASs. Admission and dynamic resource leasing are performed by a Network Slice Broker in the operator’s NMS. Management functions include CSMF for service-to-slice translation, NSMF for NSI life-cycle orchestration, NSSMF for subnetwork management, NFMF for FCAPS at network-function level, and ETSI NFV-MANO components NFVO, VNFM, and VIM. A delegated tenant-MANO is exposed through a Management Level Agreement to give tenants autonomy over selected management functions while central administrative control remains with the operator.

Each end-to-end slice is composed of 5GC, NG-RAN, and TN subnetworks. The paper explicitly notes that UEs can be served by up to eight NSIs concurrently over a single NG-RAN, and that a common AMF may be shared while SMFs remain separate per slice. In GN-NSI clustering, this concurrency bound constrains the number of S-NSIs per UE, and the protocol therefore recommends defining the cluster with up to eight S-NSIs.

The life-cycle is organized into ten steps. Step 0 captures the vertical’s portfolio of use cases, business and technical requirements, regulatory constraints, and existing catalog entries, then designs and catalogs US-NSITs or GN-NSIT plus S-NSITs. Step 1 translates tenant service requests and SLA intents into slice requirements through CSMF and, for GN-NSI, establishes Master/Slave control hierarchy. Step 2 decides among US-NSI, GN-NSI, or a hybrid. US-NSI is preferred when stringent independent isolation, different life-cycles, or different administrative ownership dominate; GN-NSI is preferred when commonalities, shared resources, and single-product governance dominate. Step 3 derives 5GC, NG-RAN, and TN subnetworks and determines what is shared or dedicated, such as shared AMF, dedicated SMF, shared gNB CU, DU partitions, and dedicated TN QoS pipes. Step 4 performs admission and instantiation through NSB, NSMF, NSSMF, MANO, and NFMF. Step 5 handles vertical and horizontal scaling. Step 6 coordinates sharing and reuse within GN-NSI. Step 7 extends to multi-domain federation where necessary. Step 8 enforces SLA assurance and run-time admission control. Step 9 implements monitoring and closed-loop automation. Step 10 adds new use cases either by instantiating new US-NSIs or extending a GN-NSI cluster with new S-NSIs.

The protocol includes explicit run-time checks. A per-use-case SLA indicator is defined as
\[
\mathcal{S}_u(x)=\mathbf{1}\big[L_u(x)\le L_u^{\max} \land R_u(x)\ge R_u^{\min} \land T_u(x)\ge T_u^{\min} \land J_u(x)\le J_u^{\max}\big],
\]
and availability for a path set \(\mathcal{P}_u\) is
\[
A_u = 1 - \prod_{i\in \mathcal{P}_u}(1-a_i).
\]
The resource-allocation model minimizes domain costs and penalties subject to latency, reliability, throughput, jitter, availability, mobility, bandwidth, compute, isolation, cluster assignment, UE concurrency, and template-compliance constraints.

Domain examples illustrate how the framework is meant to grow in practice. For automobile/V2X, autonomous driving is mapped to uRLLC with latency \(1\)–\(10\) ms, reliability \(99.999\%\), and availability \(99.9999\%\); infotainment is mapped to eMBB with latency \(<20\) ms and up to \(100\) Mbps; remote diagnostics is mapped to mMTC with latency \(<100\) ms. Manufacturing and power-grid examples similarly combine uRLLC, eMBB, and mMTC use cases under either US-NSI or GN-NSI. The principal limitations stated are the single-operator assumption, the challenge of clustered NG-RAN slicing, central MANO bottlenecks, the need for robust MLA governance in t-MANO, the eight-slice concurrency bound, and the lack of a detailed quantitative comparison between US-NSI and GN-NSI.

## 6. 3D semantic occupancy: height-aware multimodal representation and staged training

In "SliceSemOcc: Vertical Slice Based Multimodal 3D Semantic Occupancy Representation" [2509.03999], the protocol defines a vertical-slice-based multimodal framework for 3D semantic occupancy. The task is to predict a semantic label for each voxel in a fixed 3D grid from synchronized multi-view RGB images and LiDAR. The voxel feature tensor is
\[
F\in \mathbb{R}^{B\times C\times X\times Y\times Z},
\]
and semantic logits are
\[
O\in \mathbb{R}^{B\times K\times X\times Y\times Z}.
\]
The benchmark settings are explicit: on nuScenes-SurroundOcc the spatial range is \([-50,50]\) m in \(X\) and \(Y\), \([-5,3]\) m in \(Z\), with grid \(200\times 200\times 16\); on nuScenes-OpenOccupancy the range is \([-51.2,51.2]\) m in \(X\) and \(Y\), \([-5,3]\) m in \(Z\), with grid \(512\times 512\times 40\).

Image features are extracted by ResNet-50 or ResNet-101 with FPN. Visual features are lifted to the voxel grid by projecting each voxel center into each camera with camera intrinsics and extrinsics and bilinearly sampling the corresponding image feature. LiDAR uses 10 sweeps voxelized into the same grid and encoded by a VoxelNet or 3D sparse-convolution backbone to obtain a geometric feature volume \(F_L\).

The protocol’s central innovation is explicit extraction of height-axis slices. Global vertical slices are obtained by pooling over \(X\) and \(Y\) at each height layer:
\[
s_{b,z,c}=\frac{1}{XY}\sum_{x=1}^X\sum_{y=1}^Y F_{b,c,x,y,z},
\]
forming \(S\in\mathbb{R}^{B\times Z\times C}\), followed by a height-preserving 1D convolution along \(z\), implemented as a \(3\)D convolution with kernel \((1,1,3)\). Local vertical slices are extracted per \((x,y)\) column and convolved along \(Z\), then restricted to predefined height bands. For nuScenes-SurroundOcc, the data-driven bands are
\[
[-5,-3],\ [-3,-2],\ [-2,-1],\ [-1,0],\ [0,1],\ [1,3]\ \mathrm{m}.
\]
All slice operations preserve height indexing.

Before slice extraction, the feature volume is reweighted by SEAttention3D, a height-aware alternative to SENet. Instead of pooling over all spatial dimensions, SEAttention3D squeezes over \(X,Y\) only:
\[
s_{b,z,c}=\frac{1}{XY}\sum_{x=1}^X\sum_{y=1}^Y V_{b,c,x,y,z},
\]
then applies a per-height MLP
\[
a_{b,z,:}=\sigma(W_2 \delta(W_1 s_{b,z,:})),
\]
and broadcasts back:
\[
\tilde V_{b,c,x,y,z}=a_{b,z,c}\, V_{b,c,x,y,z}.
\]
The parameter count is \(2C^2/r\) per block, and the FLOPs are \(O(B Z C^2/r)\) for the MLP plus \(O(B C X Y Z)\) for pooling and scaling.

Global and local slice features are aligned by compact \(3\)D convolutions, converted to saliency maps,
\[
A_\text{global}=\sigma(\mathrm{BN}(\mathrm{Conv3D}(F_\text{global}; C\to 1, k=1))),
\]
\[
A_\text{local}=\sigma(\mathrm{BN}(\mathrm{Conv3D}(F_\text{local}; C\to 1, k=1))),
\]
and fused bidirectionally through
\[
F'_\text{global}=A_\text{local}\odot F_\text{global},\qquad
F'_\text{local}=A_\text{global}\odot F_\text{local}.
\]
The concatenated result is projected back to \(C\) channels and passed again through SEAttention3D. The full multimodal pipeline processes camera and LiDAR branches separately through this Vertical Slice Fusion, concatenates the outputs, and decodes them to per-voxel class logits.

Training is staged. Stage 0 warms up the image and LiDAR encoders for 5–10 epochs without VSF. Stage 1 enables SEAttention3D and the global slice path for 5 epochs while freezing the local branch. Stage 2 adds local slices for 10 epochs without cross-attention. Stage 3 enables bidirectional cross-attention for 10–20 epochs and trains the full model. The loss is
\[
L_\text{total}=L_\text{focal}+L_\text{lovasz}+\lambda_\text{geo}L_\text{scal}^\text{geo}+\lambda_\text{sem}L_\text{scal}^\text{sem},
\]
with softmax focal loss, Lovász-Softmax loss, and MonoScene-style affinity losses. The optimizer is Adam or AdamW with learning rate \(2\times 10^{-4}\), cosine decay, weight decay \(10^{-4}\), and batch size typically \(2\)–\(8\). Modality dropout with probability \(0.1\) is used for robustness.

On nuScenes-SurroundOcc validation, the protocol reports \(28.2\%\) mIoU. On nuScenes-OpenOccupancy validation, it reports \(22.9\%\) mIoU. The ablation trajectory on nuScenes-SurroundOcc is explicit: baseline without slices \(26.98\%\); local slices only \(27.24\%\); global slice only \(27.41\%\); both without cross-attention \(27.97\%\); both with global–local cross-attention \(28.21\%\); and SENet versus SEAttention3D improves from \(28.08\%\) to \(28.21\%\). Uniform 8 bands of 1 m yield \(27.67\%\), coarser 4 bands yield \(27.75\%\), and the data-driven 6-band strategy is best at \(28.21\%\). The paper attributes gains especially to small-object categories such as barrier, bicycle, motorcycle, pedestrian, and traffic cone.

Efficiency costs are also explicit. Memory scales as \(O(B C X Y Z)\); for SurroundOcc with \(C=128\), \(X=Y=200\), \(Z=16\), and \(B=1\), a single tensor is about \(81.9\) million floats, or about \(327\) MB in fp32. VSF adds about \(+0.6\) GB GPU memory. The paper recommends mixed precision, BatchNorm3d or GroupNorm with synchronized statistics, reflective padding along \(Z\), and careful extrinsics synchronization when temporal fusion is used. The main stated limitations are confusion among flat-ground classes and the linear growth of memory and runtime with increasing \(Z\).

## 7. Comparative structure, misconceptions, and limitations

Taken together, these sources suggest that Vertical-Slice Growth Protocol is not a discipline-specific theorem but a recurring methodological pattern [2309.07540], [1911.11135], [1211.2067], [2104.04743], [2509.03999]. In every case, the protocol serves as a translation layer between a compact formal representation and an executable workflow. In indoor farming, the compact representation is a three-state crop model and three daily controls. In CVD growth, it is a free-boundary layer system coupled through attachment kinetics. In atmospheric science, it is an \(x\)–\(z\) slice with semidirect-product structure and conserved PV. In 5G slicing, it is a catalog-and-life-cycle abstraction for US-NSIs, GN-NSIs, and S-NSIs. In semantic occupancy, it is a height-resolved voxel tensor with explicit vertical feature extraction.

A common misconception would be to treat all uses of *vertical slice* as geometrically identical. The surveyed literature does not support that. In four of the five cases, vertical refers to a spatial height or stacking coordinate; in the 5G case, it refers to a vertical industry. Another misconception would be to interpret the protocols as purely theoretical. Each source instead pushes toward implementation: indoor farming maps SRAD to DLI and PPFD practice; the vdW model maps \(\Delta\mathcal{E}\), \(F_0\), and \(k_i^\pm\) to CVD knobs; atmospheric slice theory is turned into a numerical experiment design; 5G slicing is mapped to NSMF, NSSMF, MANO, and NEF flows; and SliceSemOcc specifies pseudocode, training stages, and memory costs.

Their limitations are strongly domain-specific. The indoor farming protocol omits nutrients, vernalization, explicit HVAC, and realistic wheat economics. The vdW growth model neglects strain and mismatch, often sets desorption to zero, and primarily studies two layers. The compressible atmospheric slice model preserves energy and PV but is not a restriction of the full three-dimensional compressible equations. The 5G slicing framework assumes a single operator and leaves clustered NG-RAN slicing and comparative business quantification for future work. SliceSemOcc remains computationally heavy and still struggles on flat-ground classes. A plausible implication is that the term is most useful when understood operationally: it identifies a staged, vertically organized control or representation scheme whose value depends less on nomenclature than on how successfully it links model structure to measurable execution.

Source: https://www.emergentmind.com/topics/vertical-slice-growth-protocol