diff --git a/honeycomb_strip_crossing_mass_v2.txt b/honeycomb_strip_crossing_mass_v2.txt
index 86882ee..1c5fc38 100644
--- a/honeycomb_strip_crossing_mass_v2.txt
+++ b/honeycomb_strip_crossing_mass_v2.txt
@@ -20,7 +20,7 @@ The reciprocal is the honeycomb connective constant, shown equal to
 sqrt(2 + sqrt 2) by Duminil-Copin and Smirnov. Thus rho is the largest
 weight for which the path sums below converge.
 
-https://arxiv.org/pdf/1007.0575
+./pdf/unknown-nd-1007.0575.pdf
 
 Let S_N be the infinite horizontal strip of N triangular bands.
 Consecutive band boundaries are sqrt(3)/2 apart. Index the bottom ports by
@@ -60,7 +60,7 @@ Lawler, Schramm and Werner supplied the conformally invariant scaling
 picture, with boundary exponent 5/8, predicting the -5/4 boundary mass and
 hence the -1/4 strip power.
 
-https://arxiv.org/pdf/math/0204277
+./pdf/unknown-nd-math0204277.pdf
 
 Duminil-Copin and Smirnov proved the critical point by a parafermionic
 observable whose local cancellation yields a positive boundary identity,
@@ -77,7 +77,7 @@ quantitative sub-ballisticity.
 
 https://doi.org/10.1214/19-AIHP1024
 
-https://arxiv.org/pdf/2310.17299
+./pdf/unknown-nd-2310.17299.pdf
 
 The route computes m_N first and recovers B_N from the increments
 m_{N+1} - m_N. It passes from finite-width strip transfers to a polynomial
@@ -177,7 +177,7 @@ Yang-Baxter operators are algebraic, and the gauge-invariant periodic
 partition functions must be distinguished from other gauge-dependent
 quantities, as emphasised by Grimm and Pearce.
 
-https://arxiv.org/pdf/hep-th/9303161
+./pdf/unknown-nd-hep-th9303161.pdf
 
 The two-slot operator R(v) has inputs (south, west) and outputs (east,
 north), with coefficients 1, A(v), B(v), U(v), E(v), F(v) on the empty
@@ -366,11 +366,11 @@ points forces equality by interpolation. Exchange relations, degree bounds
 and reductions at special parameters have this role in earlier dense and
 dilute loop calculations.
 
-https://arxiv.org/pdf/math-ph/0410061
+./pdf/unknown-nd-math-ph0410061.pdf
 
-https://arxiv.org/pdf/1411.7020
+./pdf/unknown-nd-1411.7020.pdf
 
-https://arxiv.org/pdf/math-ph/0504032
+./pdf/unknown-nd-math-ph0504032.pdf
 
 Write X = (u_1, ..., u_N) and a_i = a(u_i). Define
 
@@ -653,9 +653,9 @@ Related finite-strip current calculations use bilinear link-state pairings,
 fusion reductions and interpolation, though those O(1) current formulas
 assume a cancellation of a normalising factor or an extra symmetry.
 
-https://arxiv.org/pdf/1004.4037
+./pdf/unknown-nd-1004.4037.pdf
 
-https://arxiv.org/pdf/1510.02721
+./pdf/unknown-nd-1510.02721.pdf
 
 Arch formula. Let eta = 2d - 1 = sqrt 2 - 1, so that c = eta C and
 C = eta(2C + c), and let
@@ -962,7 +962,7 @@ positive Vandermonde determinants, prod x_i, positive weights, and
 
 which is positive by clearing denominators and checking the constant one.
 
-https://arxiv.org/pdf/0904.2602
+./pdf/unknown-nd-0904.2602.pdf
 
 A solution of the homogeneous problem for g, of degree at most N - 1, has
 [b^(4N+1)] Q = 0, but [b^(4N+1)] Q = (-8)^N D_0, so D_0 = 0. Positivity of M
@@ -1057,7 +1057,7 @@ one-variable weight exponent is a - 1, so this a exceeds the parameter of
 Forrester and Kieburg by one; the ordered normaliser is the unordered
 integral divided by n!.
 
-https://arxiv.org/pdf/1410.6883
+./pdf/unknown-nd-1410.6883.pdf
 
 Association and monotone comparison. Let I be an interval in (0, infinity),
 n >= 1, and w_1, ..., w_n positive smooth on I. If
@@ -1080,7 +1080,7 @@ increasing functions. Conditional on the largest coordinate t, the rest have
 a law of the same form on the shorter interval with modified one-variable
 factors, so they are associated by induction. For s < t the ratio of the
 conditional density at s to that at t, on the former support, is
-coordinatewise decreasing by the mixed-derivative inequality, and extended
+coordinatewise decreasing by the mixed-derivative ineuality, and extended
 by zero where a coordinate exceeds s it stays decreasing. Association at t
 makes a decreasing likelihood tilt decrease every increasing expectation;
 an unbounded ratio is truncated and handled by monotone convergence. Hence
@@ -1095,7 +1095,7 @@ expectations, nonnegative by
 
 T' an independent copy of T. Truncation and integrable approximation extend
 the result to the stated class. Finally, for an increasing likelihood ratio
-L, E[F L]/E L >= E F for increasing F, with the reversed inequality for
+L, E[F L]/E L >= E F for increasing F, with the reversed ineuality for
 decreasing L, which is the stochastic comparison. End of proof.
 
 This is the continuous MTP_2 principle of Sarkar, stated for log-supermodular
@@ -1142,7 +1142,7 @@ pairing two opposite arguments gives that constant recursively as one; for
 odd size add an argument tending to infinity, whose augmenting column
 consists of ones.
 
-https://arxiv.org/pdf/math/0411280
+./pdf/unknown-nd-math0411280.pdf
 
 Expand the determinant and the Pfaffian and integrate over the unordered
 region, divided by n!. This is de Bruijn's determinant-Pfaffian integration
@@ -1150,9 +1150,9 @@ formula, with the general antisymmetric-kernel form including its odd-size
 border; absolute integrability follows from |(s - r)/(s + r)| <= 1 and finite
 gamma moments for a > 0.
 
-https://pure.tue.nl/ws/files/1920642/597510.pdf
+./pdf/unknown-nd.pdf
 
-https://arxiv.org/pdf/math/9905083
+./pdf/unknown-nd-math9905083.pdf
 
 The double moment with indices i, j is