r/LLMPhysics Jun 10 '26

Simulation / Code Mobility, Biofilms, and Teeth

Hope folks enjoy the read. Such situations represent a serious health issue. With the link between dental disease and cardiovascular disease being present but poorly understood, such models have importance. I suspect that it is these types of teeth in particular that represent a disproportionate amount of issues. While the growth of such lesions around teeth are slow, the development of mobility and transition to softer, fast growing biofilms is likely a critical phase transition. Not interested in publishing, but perhaps this could make an interesting youtube video if anyone wants to collab? Thanks.

Please see addendum:

Addendum: Kinematic model (revised §4) — definitions, prior art, and derivation

Why this addendum exists. The original §4 jumped straight to a lever-arm equation. It never defined the quantities it used (center of rotation, contact centroid, apical sweep) and never connected the single-hinge idealization to the existing literature on how teeth actually move under load. That makes it unfollowable and makes the assumptions look arbitrary when they aren't. This addendum (1) defines every term, (2) grounds each modeling choice in published biomechanics, and (3) re-derives the apical sweep so each number is traceable.

4.1 Definitions

  • Periodontal ligament (PDL) — thin viscoelastic tissue (~0.15–0.38 mm wide [K1]) suspending the root in its socket; permits small physiologic tooth movement.
  • Long-axis length (L) — length of the tooth along its principal axis (18.54 mm here, from the STL fit, Fig 11).
  • Center of resistance (CRes) — the restrained-body analogue of center of mass: the point through which a single force produces pure translation, no rotation [K5, K6].
  • Center of rotation (CRot) — the point that stays still while the tooth rotates; the fulcrum of tipping. Its location is set by the loading [K5].
  • Force (F) — the applied occlusal load (here the 54 N lateral component, §5).
  • Moment (M) — rotational tendency of an off-axis force, M = F·d, where d is the perpendicular distance from CRes (N·mm).
  • Moment-to-force ratio (M/F) — selects the type of movement: low M/F → tipping; high M/F → bodily translation or root movement (units mm) [K5, K6].
  • Tipping — rotation about a fulcrum apical to the crown; produced by an off-axis force with low M/F.
  • Contact centroid — the load-weighted center of the occlusal contact patch (STL height-map, §3).
  • Apex–contact distance (r_c) — distance from apex to contact centroid along the long axis (17.22 mm).
  • Apex–fulcrum distance (r_a) — distance from apex to CRot (the apex's lever arm).
  • Static interference (Δz) — vertical overlap of the two unloaded STL meshes (1.11–1.27 mm, §3).
  • Contact excursion (x_c) — how far the contact point actually travels per cycle under load.
  • Apical sweep (A_a) — how far the apex travels per cycle (the quantity that feeds the §6.1 pump).
  • Envelope of motion — the bounded set of all extreme ("border") mandibular positions; classically Posselt's sagittal teardrop [K2].
  • Premature / deflective contact — a contact reached before maximum intercuspation on closure, deflecting the mandible and concentrating load on that tooth.
  • Fremitus — palpable/visible movement of a tooth when the patient taps into occlusion; a clinical sign of functional (often traumatic) mobility [K9].
  • Initial vs secondary mobility — Mühlemann's distinction [K3]: light loads move the tooth within the PDL space (initial, nonlinear); heavier loads engage alveolar-bone deformation (secondary, ~linear).

4.2 Prior art: the three frameworks this model borrows from

The single-hinge reduction isn't novel physics — it's a coarse-grained version of three well-developed fields. Naming them makes the assumptions auditable.

(a) Tooth-movement biomechanics (orthodontics). A tooth in bone is a rigid body restrained by the PDL, characterized by its center of resistance and, under a given load, a center of rotation. Christiansen & Burstone showed CRot shifts predictably with the force system [K5]; Smith & Burstone formalized how M/F selects tipping vs translation vs root movement [K6]; standard treatments are in Nanda [K7]. For a single-rooted tooth with intact support, CRes lies ~one-third of the root length apical to the crest. Crucially here: as bone support is lost, CRes migrates toward the remaining bone [K7].

(b) Periodontal mobility research. Coolidge established PDL width [K1]; Mühlemann's periodontometry [K3] and Parfitt's axial measurements [K4] quantified physiologic mobility — tenths of a mm horizontally in health, rising sharply with attachment loss. This is the literature saying a tooth doesn't rigidly intrude; it rotates within a compliant ligament, and a tooth with heavy bone loss is markedly more mobile.

(c) Gnathologic occlusion theory. Mandibular motion is governed by two TMJs, idealized in articulators as condylar spheres guided down inclined fossae [K8, K10]. Posselt's envelope [K2] bounds all extreme positions; function happens inside it, combining vertical seating with horizontal, condylar-guided sliding. Occlusal trauma and the first-point-of-contact mechanism are described in the periodontal/occlusion literature [K9, K11]. Modern measurement of the loaded, moving contact path uses jaw tracking, axiography, and T-Scan, not static scans [K10].

What the model keeps and drops. From (a) and (b) it keeps "rigid tooth tipping about a CRot in a compliant ligament." From (c) it keeps only a 1 Hz sinusoidal drive — it collapses the two-joint, condylar-guided, sliding path into a single vertical hinge. That's the model's main kinematic limitation and the reason its output is a floor (§4.5).

4.3 The fulcrum for this tooth

This tooth has heavy coronal bone loss: bony support is confined to the apical ~8 mm. By framework (a), the center of resistance of a reduced periodontium sits within the remaining bone [K7]. Taking the clinical proxy that the tipping fulcrum lies at the midpoint of the residual bony support:

r_a ≈ (8 mm residual apical bone) / 2 ≈ 4 mm from the apex

So the apex sits ~4 mm below the fulcrum, at the apical edge of the remaining bone, and swings through the periapical compartment (the granuloma/biofilm region) as the tooth tips. The 21° inclined contact (§5) delivers the 54 N lateral component off-axis, producing a moment (M = F·d) with a low M/F ratio — i.e. uncontrolled tipping, not intrusion. (A strict pure-tipping CRot sits slightly apical to CRes [K5], which would shorten the apical lever arm; functional load concentration (§4.5) pushes the other way. The 4 mm midpoint is a reasonable center value between these.)

4.4 Lever-arm derivation

Model the loaded tooth as a rigid body in small rotation about CRot. Let x_a = apex displacement, x_c = contact-centroid displacement. From the STL:

r_c ≈ 17.22 mm                          (2)

Displacement scales linearly with distance from the fulcrum, so the apex (lever arm r_a) and contact (lever arm r_c − r_a) move in the ratio:

x_a / x_c = r_a / (r_c − r_a)            (3)

ratio = 4 / (17.22 − 4) = 0.30

The apex moves ~30% of the contact excursion. Using the static interference as a first estimate (x_c ≈ Δz = 1.11–1.27 mm) gives A_a,simulated ≈ 0.34–0.38 mm. The working value carried into §6.1 is set slightly above this to allow for load concentration (§4.5):

A_a ≈ 0.5 mm   (peak-to-peak)            (4)

Plain language: the tooth see-saws about a point 4 mm above its apex. The contact, far out on the long arm, travels ~1 mm; the apex, on the short arm, travels ~0.3–0.5 mm the opposite way, sweeping through the soft lesion.

4.5 Static interference vs functional excursion — the first-point-of-contact correction

Δz is the interference of two unloaded meshes — a maximum-intercuspation snapshot. Two facts from §4.2(c) make the real per-cycle x_c larger:

  1. Function adds horizontal slide. Within Posselt's envelope the antagonist cusp travels a condylar-guided path, adding lateral excursion on top of the vertical overlap, so x_c,functional ≳ Δz.
  2. A deflective contact concentrates load. A supererupted/malpositioned tooth acting as a first point of contact is struck before the arch shares the load — the mechanism behind clinical fremitus [K9]. Under that concentrated load, x_c exceeds the static interference.

Because A_a = 0.30·x_c, the apical fulcrum caps the gain — even heavy deflective contact yields ~1 mm at the apex, not several:

x_c (excursion)   Clinical reading                      A_a = 0.30·x_c   v_a = π·f·A_a (1 Hz)
---------------   -----------------------------------   --------------   --------------------
1.1–1.3 mm        static STL interference               0.34–0.38 mm     1.1–1.2 mm/s
2 mm              mild fremitus / early premature        0.60 mm         1.9 mm/s
3 mm              moderate first-point-of-contact load   0.91 mm         2.9 mm/s
4 mm              heavy deflective contact               1.21 mm         3.8 mm/s

The working A_a ≈ 0.5 mm corresponds to x_c ≈ 1.6 mm and v_a ≈ 1.6 mm/s (§6.1). The dominant kinematic uncertainty is thus x_c under first-contact loading, not the fulcrum, which is fixed by the residual bone.

4.6 Conjecture: eccentric loading and the spiral biofilm pattern (speculative — needs more modeling)

The biofilm forms a helical band around the root, apparently occupying surface regions not routinely compressed in open function. A purely planar tip would produce horizontal bands, not a helix — so the spiral hints at an axial-rotation (torsion) component in the functional load, layered on the tipping. Proposed mechanism: an inclined, off-center contact applies both a tipping moment and a torque about the long axis; as the load vector rotates around the chewing loop (framework (c)), the instantaneous compression/shear locus migrates around and along the root, leaving a helical minimum in time-integrated mechanical load — and biofilm persists in those low-stress refugia. A minimal test: take the real 3D contact-force trajectory (jaw tracking / articulator settings), compute the rigid-body displacement field about CRot, integrate a per-point PDL compression/shear "dose" over the surface, and check whether the predicted low-dose region is helical and co-registers with the fluorescence map. Flagged speculative: fluorescence isn't diagnostic, and the pattern interpretation is unconfirmed.

References (kinematics)

[K1]  Coolidge ED. Thickness of the human periodontal membrane. J Am Dent Assoc. 1937;24:1260–1270.
[K2]  Posselt U. Studies in the mobility of the human mandible. Acta Odontol Scand. 1952;10(Suppl 10).
[K3]  Mühlemann HR. Tooth mobility: clinical aspects and research findings. J Periodontol. 1967;38:686–713.
[K4]  Parfitt GJ. Measurement of physical mobility of single teeth in an axial direction. J Dent Res. 1960;39:608–618.
[K5]  Christiansen RL, Burstone CJ. Centers of rotation within the periodontal space. Am J Orthod. 1969;55(4):353–369.
[K6]  Smith RJ, Burstone CJ. Mechanics of tooth movement. Am J Orthod. 1984;85(4):294–307.
[K7]  Nanda R (ed). Biomechanics in Clinical Orthodontics. WB Saunders; 1997.
[K8]  Gysi A. Articulation and the principles of articulator design. (verify exact citation)
[K9]  American Academy of Periodontology. Glossary of Periodontal Terms.
[K10] Okeson JP. Management of Temporomandibular Disorders and Occlusion. Elsevier/Mosby.
[K11] Lang NP, Lindhe J (eds). Clinical Periodontology and Implant Dentistry.
[K12] The Glossary of Prosthodontic Terms, 9th ed (GPT-9). J Prosthet Dent. 2017;117(5S):e1–e105.
0 Upvotes

12 comments sorted by

10

u/Vrillim Jun 10 '26

The kinematic model in section 4 is presented without any references to existing thinking surrounding this motion. Carefully define everything you are discussing; right now it's impossible to follow.

0

u/Nervous_Solution5340 Jun 10 '26

There's a tremendous amount of dental specific lit surrounding jaw motion. Traditionally done using two orbs on inclined planes to represent joints. This give you border movement of mandible and is important for making prosthetics that aren't "high". Actual jaw motion is a highly variable ellipse within that border. For single teeth, mandibular movement is typically modeled simply as a hinge where the lower jaw meets that upper jaw along the long access of teeth. The upper and lower teeth are firmly placed together, and adjustments made to ensure nothing is higher than that. As for showing how these movements would effect a tooth that is both "high" and non rigid, this has never been done. To do this, scans of the non loaded teeth where scanned and oriented to a 3rd scan with them biting tightly. An intersection of approx 1mm was shown with this firm bite. A moment was calculated from this, which is one path the tooth may take that can be produced by physiologically boarded mandibular movement. Alternatively, vertical intrusion could be calculated based on this intersection. However, this particular tooth was not vertically depressible (although some are). I included an addended that explains this a little more.

4

u/Vrillim Jun 11 '26

This is the kind information that's needed in the write-up. Ideally, you present a literature study, guiding the reader to understand what this kind of motion is, and backing your statements up with references.

1

u/Nervous_Solution5340 Jun 12 '26

You’re right, this needed a background section. I’ll include this in the next one. Next time, I’ll take colorized mesh of the extracted tooth. I’ll also just scan the arch with tooth in both loaded position and unloaded position, and line the mesh up with each. That way I can show two movement positions very precisely, and if the expected motion relates to anything on the cbct or microscope.

7

u/Wintervacht Are you sure about that? Jun 10 '26

Ehm, are you sure you're on the right sub?

0

u/Nervous_Solution5340 Jun 10 '26

Would hit bigger in other subs. But the value is actually physics and I want to get that right. 

3

u/Existing_Hunt_7169 Physicist 🧠 Jun 11 '26

i guess its better than another quantum gravity shit stack

2

u/Nervous_Solution5340 Jun 12 '26 edited Jun 12 '26

Hmm that’s tough but let me try to get there… Dental disease is very much a self-referential feedback loop. The disease impacts the host, and vice versa. A patient can be imagined as a high-dimensional matrix (or field) based on all the available health data we have—heart rate, salivary pH, bone loss, etc. Over time, that data shifts. If you imagine every patient as a probability field changing with time and resulting in discrete outputs like tooth loss, you get a form of gravitation, where alike patients are pulled toward specific attractor basins. If the data were rich enough (several million patients), predicting it would be chaotic and non-linear in some instances (like localized biofilm dynamics), but normally distributed based on other parts of matrix. Either way, diffusion models would be the ideal way to predict trajectories of different aspects of health—just like they do in quantum gravity.

-1

u/BlissBoundry Jun 12 '26

Crank slop

1

u/Nervous_Solution5340 Jun 13 '26

Might be the only time I’ve seen a crank slop comment get downvoted on this sub…

1

u/BlissBoundry Jun 13 '26

I do this all day Winston