Triangular Lattice Filling Fraction Calculator
Compute area filling fraction for disks on a 2D triangular lattice with occupancy and geometry controls.
Formula used: η = f × πr² / ((√3/2)a²), where f = occupancy fraction.
How to Calculate Filling Fraction in a Triangular Lattice: Complete Expert Guide
If you are studying condensed matter physics, materials science, soft matter, photonic crystals, or nanostructured surfaces, you will repeatedly encounter one geometric quantity: the filling fraction. In a two-dimensional triangular lattice, filling fraction tells you what percentage of the lattice area is occupied by particles, pores, inclusions, cylinders, or disks. It sounds simple, but accurate filling-fraction calculations are foundational for predicting transport properties, optical response, mechanical stiffness, and thermodynamic phase behavior.
This guide shows exactly how to calculate filling fraction in a triangular lattice, step by step, from first principles. You will also learn how occupancy defects change the result, why the triangular lattice gives the highest monodisperse disk packing in 2D, and how to avoid common mistakes that produce unphysical values above 1.0. Use the calculator above for fast estimates, then use the derivation below when you need publishable rigor.
1) What Filling Fraction Means in a Triangular Lattice
In 2D, filling fraction (often written as η, ϕ, or f depending on field convention) is:
- Filled area divided by total reference area.
- Dimensionless and typically reported from 0 to 1 or 0% to 100%.
- Dependent on geometry (radius, spacing) and occupancy (how many lattice sites are actually filled).
For a triangular lattice of identical circular particles, one natural reference area is the primitive unit cell: Acell = (√3/2)a², where a is the nearest-neighbor lattice constant. If one disk of radius r is associated with each occupied lattice site, area contributed per site is Adisk = πr².
If only a fraction of sites are occupied, define occupancy fraction as f = occupancy% / 100. Then: η = f × (πr²) / ((√3/2)a²).
2) Core Derivation from Lattice Geometry
A triangular Bravais lattice is generated by two equal-magnitude basis vectors separated by 60°. The area of the primitive parallelogram is the magnitude of the cross product:
Acell = |a1 × a2| = a² sin(60°) = (√3/2)a².
For monodisperse disks centered on lattice points:
- Compute particle area per occupied site: Adisk = πr².
- Compute cell area: Acell = (√3/2)a².
- Apply occupancy factor f if vacancies are present.
- Evaluate η = f × Adisk / Acell.
In the special touching limit, nearest neighbors just touch, so a = 2r. Substituting gives:
ηmax,2D triangular = πr² / ((√3/2)(2r)²) = π/(2√3) ≈ 0.9069.
That value is the well-known densest packing fraction for identical circles in two dimensions under periodic triangular ordering.
3) Worked Numerical Example
Suppose your lattice constant is a = 500 nm, particle radius is r = 210 nm, and 92% of lattice sites are occupied.
- Adisk = π(210²) = 138,544 nm² (approx).
- Acell = (√3/2)(500²) = 216,506 nm² (approx).
- Raw geometric ratio = 138,544 / 216,506 = 0.6399.
- Apply occupancy f = 0.92: η = 0.92 × 0.6399 = 0.5887.
Final filling fraction is 58.87%. Void fraction is 41.13%. This direct calculation is exactly what the calculator performs.
4) Comparison Table: Common 2D Packing Statistics
The table below compares widely used 2D packing references. These are standard geometric or simulation-backed benchmarks frequently used in materials and statistical mechanics discussions.
| Packing or State | Expression / Typical Value | Decimal Value | Interpretation |
|---|---|---|---|
| Triangular lattice (monodisperse, touching) | π/(2√3) | 0.9069 | Maximum ordered 2D circle packing fraction |
| Square lattice (monodisperse, touching) | π/4 | 0.7854 | Lower packing efficiency than triangular lattice |
| Random loose hard-disk regime | Typical simulation range | ~0.77 to 0.80 | Disordered low-density jammed-like states |
| Random close-like disordered hard-disk regime | Typical simulation range | ~0.82 to 0.84 | Dense disordered assemblies without long-range crystalline order |
5) Occupancy and Defects: Why Real Samples Differ from Ideal Geometry
In laboratory systems, perfect site occupation is rare. Missing particles, substitutional defects, polydispersity, and domain boundaries all reduce effective filling fraction. If your imaging analysis identifies occupied lattice sites, you should include occupancy explicitly rather than forcing an idealized 100% value.
For example, if your ideal touching triangular geometry gives η = 0.9069 and your occupancy is 95%, then effective area fraction becomes 0.8616. This is a meaningful shift that can alter optical bandgaps, permeability, and percolation thresholds in patterned media.
- Use occupancy correction when analyzing microscopy-derived lattice maps.
- Report both geometric packing fraction and defect-corrected effective fraction.
- Track uncertainty in radius and lattice constant; η is sensitive to both.
6) Sensitivity Analysis and Error Propagation
Because η scales as r²/a², small measurement errors can produce amplified uncertainty. A practical differential estimate:
Δη/η ≈ 2(Δr/r) + 2(Δa/a) + (Δf/f), using conservative linearized propagation.
If r has 2% uncertainty, a has 1.5% uncertainty, and occupancy has 1% uncertainty, relative error in η can be around 8% using worst-case addition. If you use independent random-error propagation (root-sum-square), uncertainty is smaller, but still nontrivial. This is why calibrated imaging and robust segmentation matter in quantitative materials characterization.
7) Comparison Table: Example Triangular-Lattice Scenarios
| Case | a | r | Occupancy | Computed η | Void Fraction (1 – η) |
|---|---|---|---|---|---|
| Ideal touching crystal | 1.00 | 0.50 | 100% | 0.9069 | 0.0931 |
| Touching with 90% occupancy | 1.00 | 0.50 | 90% | 0.8162 | 0.1838 |
| Non-touching sparse geometry | 1.00 | 0.40 | 100% | 0.5804 | 0.4196 |
| Dense but defective film | 1.00 | 0.48 | 93% | 0.7732 | 0.2268 |
8) Practical Use Cases Across Disciplines
In photonics, filling fraction directly controls refractive index contrast in 2D periodic media. In catalysis and porous membranes, it influences active area and flow resistance. In colloidal science, it helps determine phase state and ordering transitions. In battery electrodes and printed electronics, it contributes to conductivity pathways and mechanical integrity.
Even when your system is not perfectly circular or monodisperse, triangular-lattice filling-fraction calculations provide a useful baseline model. You can extend the same framework by replacing πr² with your measured average inclusion area and adjusting occupancy to account for defects.
9) Common Mistakes and How to Avoid Them
- Using the wrong unit cell area: triangular lattice is (√3/2)a², not a².
- Mixing diameter and radius: verify whether your measured quantity is r or 2r.
- Ignoring occupancy: vacancy-rich samples can be significantly overestimated otherwise.
- Inconsistent units: keep a and r in the same unit before computing.
- Accepting η > 1 without checking: this indicates inconsistent geometry or input error.
10) Authority References for Further Study (.gov and .edu)
For deeper context on crystal geometry, close packing, and lattice physics, consult these authoritative educational and government sources:
- Georgia State University HyperPhysics: Close Packing in Solids (.edu)
- Georgia State University HyperPhysics: Bravais Lattices (.edu)
- National Institute of Standards and Technology (NIST) (.gov)
11) Quick Recap
To calculate filling fraction in a triangular lattice, use one clean equation: η = f × πr² / ((√3/2)a²). For touching monodisperse disks at full occupancy, η reaches 0.9069. Real systems usually fall below that because of vacancies, spacing variation, and finite polydispersity. If you consistently define radius, lattice constant, and occupancy, you can produce accurate, reproducible filling-fraction values for both research and engineering applications.
Use the calculator above to run fast scenarios, compare defect levels, and visualize occupied versus void area immediately with a chart. For reporting, include assumptions, measurement uncertainty, and whether your value is geometric ideal or defect-corrected effective filling fraction.