We study how a bound notebook page can be folded once so that the page is visible when the notebook is closed and becomes a bookmark. An additional practical requirement is that the folded page stays within the same side of the bind, so that it does not get creased further when one closes the notebook. This simple problem displays a remarkably rich bifurcation behavior, which can be analyzed using undergraduate calculus. The main finding is that the optimal way of folding the page depends on the page's aspect ratio and changes abruptly when the latter exceeds the critical value of 1.207 11. This behavior also resembles a first-order phase transition in thermodynamics despite its geometric origin.

Folding, as a general method of transforming geometries, is found everywhere. Clothes, bi-fold doors, umbrellas, airbags, solar panels, and wings of insects and birds1,2 are just some common examples. Folding allows a thin object to be stored compactly while remaining easily deployable, or transformed into structures with new capabilities (e.g., origami).

In mathematics, origami, a slightly constrained form of folding, has been formalized into an axiomatic system.3 Analysis and computation aid in the creation of novel origami designs.4 Interestingly, folding has a mathematical capability that outperforms compass-and-straightedge construction and can be used to solve classical problems such as “angle trisection” or “doubling the cube.”5,6

Folding is also common in physical sciences. A prominent example is the protein folding problem, where one would like to establish the 3D protein structure from an amino acid sequence alone. Recently, deep learning allowed for a major breakthrough in answering this question.7 Interestingly, folding (together with stretching) is also a useful concept in mixing of highly viscous fluids. The diverse folding mechanisms, including mechanical, capillary, thermal, and chemical folding, can benefit a wide range of micro-scale applications.8,9 Back at the office, a violently folded (or crumpled) piece of paper has remarkable physical insight buried in its crease network.10,11

Admittedly, the mathematical, physical, or biological study of folds is full of surprises. This work presents one simple problem, in which folding a notebook page exhibits an interesting bifurcation behavior. Given a pocket notebook of length 1 in the x direction and l in the y direction (Fig. 1), a folded page can extend beyond the side bounds of its unfolded state, remaining visible when the notebook is closed and thus serving as a bookmark. Then, what is the optimal way to fold the page? Although preferences vary, we assume in this paper that the optimal fold maximizes the bookmark length e (Fig. 1), which is the length by which the page extends beyond its unfolded width along the x axis. In addition, the top edge of the page should not be torn off—a constraint similar to the one experienced by birds and insects when folding their wings due to the joint between the wings and body. Aside from this first constraint, the folded page should not extend past the top edge of the notebook, or it will be creased further when closing the notebook. This is the second constraint.

Fig. 1.

The notebook and its folded page that motivated this work. The notebook has a width 1 and a height l. The excess length e is maximized by the optimal fold (subject to constraints) so that the folded page serves as an effective bookmark.

Fig. 1.

The notebook and its folded page that motivated this work. The notebook has a width 1 and a height l. The excess length e is maximized by the optimal fold (subject to constraints) so that the folded page serves as an effective bookmark.

Close modal

We will begin by considering the first constraint alone and then solve the complete problem, using only undergraduate calculus. This paper concludes by proposing several possible directions for future explorations.

We start our analysis by considering only the first constraint. Without loss of generality, we assume that the page is bound at its top edge, and the fold is in the upper-right direction (as in Fig. 1). Then the objective is to have the folded page extend rightward as much as possible, while keeping the top edge fixed.

Using the coordinate system defined in Fig. 1, the fold maps the original rectangle Ω = [ 0 , 1 ] × [ 0 , l ] to a new region Ω . The optimal fold is the one that maximizes the excess length, which is denoted by e in Fig. 1. One may interpret Ω as one state of the system, and e as the negative of the state's potential energy. Then the optimal fold that maximizes e is the one that minimizes the system's energy.

Folding the page amounts to choose a straight crease in Ω and reflect the region below it—which we will call the “folded region”—with respect to this line. For a rectangular page, there are only two possible cases: either the crease crosses the left and bottom edges and the folded region is a triangle (Fig. 2, case 1), or it crosses the left and right edges in which case the folded region is a trapezoid (Fig. 2, cases 2a and 2b). Regardless of the crease location, e will be determined by the location of a corner point of Ω (otherwise, there is always a neighboring point that, once reflected, is further to the right). For this reason, we only need to consider the corners of the folded region.

Fig. 2.

The two cases arising when folding a page in the upper-right direction and with the top edge fixed. The dashed lines mark the crease line defined by the length parameters (a, b). The corners of the folded region are distinguished by adding primes to the symbol.

Fig. 2.

The two cases arising when folding a page in the upper-right direction and with the top edge fixed. The dashed lines mark the crease line defined by the length parameters (a, b). The corners of the folded region are distinguished by adding primes to the symbol.

Close modal

In both cases, we will first find the x coordinate of the rightmost point of the folded region, which is denoted by xe. xe may correspond to different page corners, depending on how we fold the page (refer to Fig. 2) and is related to e by e = x e 1.

Case 1:

The crease line crosses the left edge at a distance a from the lower-left corner and the bottom edge at a distance b from the lower-left corner (Fig. 2, case 1). In this case, the folded region is triangle Δ OBA, and the right most point is O . We have
(1)
x e ( a , b ) satisfies x e / a 0 always, and x e / b 0 conditionally for b [ 0 , a ]. The maxima of xe and e are achieved at the upper bounds of a and b: a = l and b = min ( l , 1 ) because b is also limited by 0 b 1.

When l 1, meaning the page is “short and wide” (including the square page shown in Fig. 1), substituting a = b = l in xe gives x e ( l , l ) = l ( denotes optimal values). As a result, e = l 1 0. This means that, for a short page, it is impossible to fold a bookmark using case 1 creases.

When l > 1, meaning the page is “long and narrow,” substituting a = l and b = 1 gives
Because e > 0, a bookmark can be folded for a long page. The optimal crease line then coincides with the page diagonal that connects the upper-left and lower-right corners. However, the limit lim l + e ( l , 1 ) = 1 is the maximum possible excess length, no matter how long the page is.

Case 2:

The crease line crosses the left edge at a distance a from the lower-left corner and the right edge at a distance b from the lower-right corner (Fig. 2, cases 2a and 2b). In this case, the folded region is trapezoid OCBA, with a b for the fold to be in the upper-right direction. We define θ as the angle between the crease line and the x-axis. Clearly, tan θ = a b. The rightmost point in Ω depends on the value of θ: It is C when θ 45 ° (Fig. 2, case 2a) and O when θ > 45 ° (Fig. 2, case 2b).

Case 2a. When θ 45 ° (i.e., a b 1), x e = x C and
(2)
With a b 1, Eq. (2) increases monotonically with a. Letting a = l and e ( l , b ) / b = 0 gives the optimal solution
(3)
Because e > 0, this fold leads to a valid bookmark regardless of the value of l. For example, if the page is square (l = 1), we have e ( 1 , 2 2 ) = 2 1 0.41421, whereas, in the same condition, case 1 leads to e = 0 (an invalid bookmark). One can easily show that, for a given height l, solution Eq. (3) always outperforms the solution of case 1.
When l is large, Eq. (3) is approximately
(4)
So for a thin paper strip (high aspect ratio), the optimal fold tends to overlap the left edge with the top edge using a crease line that passes through the upper-left corner and is at 45 ° with respect to the x-axis (the gray fold in Fig. 3). Certainly, in the limit of an infinitely long strip, the page resembles a line and this fold is optimal: It simply reorients the page horizontally.
Fig. 3.

The original (dashed), optimally folded (solid black), and nearly optimally folded (solid gray) pages with an aspect ratio of 4 (an example of Fig. 2, case 2a).

Fig. 3.

The original (dashed), optimally folded (solid black), and nearly optimally folded (solid gray) pages with an aspect ratio of 4 (an example of Fig. 2, case 2a).

Close modal
Case 2b. When θ 45 ° (i.e., a b 1), x e = x O and
(5)
With a b 1, Eq. (5) increases monotonically with b. Setting b = l 1 gives a = l, θ = 45 °, and e = l 1. Because θ = 45 °, this solution is also a special case of case 2a and is always sub-optimal compared with the solution (Eq. (3)) to case 2a.

All cases considered, the optimal fold is Eq. (3) of case 2a and resembles the illustration Fig. 2, case 2a.

We continue our analysis by including the second constraint: The folded page should stay within one side of the bind, that is to say, “below” the bind—so that it does not get accidentally creased when one closes the notebook. We start from case 2a, which gave us the optimal fold considering only the first constraint, and improve it to satisfy the second constraint. One will see that the additional constraint leads to bifurcation behaviors of the optimal solution.

In case 2a, the optimally folded page does exceed the top edge (i.e., the bind). This is shown by finding the maximum y coordinate of Ω , which is defined at point O (refer to Fig. 2, case 2a)
(6)
Substituting a and b from Eq. (3) gives
The constraint that limits the y extent reads y e l. Together with Eqs. (2) and (6), we finally have the complete problem statement as
(7)
It is important to keep in mind the physical picture behind Eq. (7), which is the maximization of the x coordinate of point C subject to the constraint imposed on the y coordinate of point O , being understood that C and O result from the isometric transformation (i.e., shape-preserving) of the page upon folding along the crease line.

The accessible region of the (a, b) phase space according to Eq. (7) is shown in Fig. 4 (shaded area) for l 1. Each point in the accessible region represents one way to fold the page without violating any of our constraints. At low a, the region is bounded by the condition b a. At higher values, it is determined by 2 a l ( 1 + ( a b ) 2 ) and is, therefore, more complicated. Still, one can show that e / a 0, so that e occurs on the curve boundary. One might also notice that e = 0 along the other two (straight) boundaries of the accessible region.

Fig. 4.

(a, b) phase diagram for l 1. The accessible region, determined in Eq. (7), corresponds to the shaded area.

Fig. 4.

(a, b) phase diagram for l 1. The accessible region, determined in Eq. (7), corresponds to the shaded area.

Close modal

The equation of the curved boundary is b = a 2 a l l 2 / l, which is defined over the interval a [ l / 2 , l ] (the other solution b = a + 2 a l l 2 / l violates b a and is discarded). Points on this boundary correspond to a special kind of folding: the ones for which the page's lower-left corner overlaps with its top edge. This physical picture explains the lower limit l / 2 for a. The remaining constraint, b a 1, is automatically satisfied on this boundary.

By substituting b = a 2 a l l 2 / l into Eq. (2), we remove the dependence on b and obtain
(8)
The aspect ratio l itself becomes a parameter that affects the optimal solution.
In Fig. 5, e(a, l) as given in Eq. (8) is plotted as a function of a for five different values of l. The local and global maxima of each curve are marked with open and solid circles. When l is small (1.0 and 1.1), we have a unique local maximum, which is also the global maximum. The maximum occurs at e / a = 0, which leads to an algebraic equation a 4 = ( 2 a l ) l. The solution a i n = a i n ( l ) (a complex algebraic expression in l that can be obtained using any computer algebra system), when substituted into Eq. (8), gives the local optimal solution
(9)
Fig. 5.

Plot of e(a, l) with different values of l (labeled by text). Each curve has its own horizontal extent a [ l / 2 , l ]. The small solid/open circles mark the local maxima and minima, and their trajectories are shown by the dotted/dashed curves, respectively. The crosses mark the boundary solutions.

Fig. 5.

Plot of e(a, l) with different values of l (labeled by text). Each curve has its own horizontal extent a [ l / 2 , l ]. The small solid/open circles mark the local maxima and minima, and their trajectories are shown by the dotted/dashed curves, respectively. The crosses mark the boundary solutions.

Close modal

Because the solutions are such that a i n ( l ) is within the corresponding range in a ( a [ l / 2 , l ]), we call them internal solutions. The crease lines intersect with the page's left edge but do not pass through either of its two end points.

As l increases, the right boundary point of the curve rises continuously (see Fig. 5). Substituting a = l in Eq. (8) gives the excess of those points as
(10)
In addition, b = l 1 and θ = 45 ° (refer to θ in Fig. 2). We call those points the boundary solutions. Their crease lines pass the top-left corner of the page as illustrated in Fig. 3 (gray fold).

The boundary and internal solutions become equal at the critical value lcr = 1.20711, which can only be found numerically by solving e ̂ i n ( l ) = e ̂ b d ( l ) or e ( a i n ( l ) , l ) = l 1. The solution steps for mathematica are provided in the  Appendix as an example.

When l is greater than lcr, the optimal solution e jumps from the internal solution to the boundary solution. Interestingly, the internal solution ceases to exist beyond a second critical value l c r 2. This is caused by the change of concavity of the e(a, l) curves in Fig. 5. The critical value l c r 2 is obtained by solving
(11)
with e(a, l) from Eq. (8). Numerically, l c r 2 = 3 3 / 4 1.29904.

To further understand the disappearance of the internal solution, Fig. 5 plots its trajectory as l is increased (dotted curve). In addition, Fig. 5 plots the trajectory of the local minima (dashed curve). It is clear that, as l increases, the two solutions approach each other. They coincide when l = l c r 2 and, as a pair, spontaneously disappear (Fig. 6 will also show this).

Fig. 6.

The optimal solution e ( l ) in the bookmark folding problem. Top: Three pages of increasing l are shown with a sequence of folds. The folds have a uniform distribution of a [ l / 2 , l ], with b chosen to have O overlap with the top edge. The thick black curves track the trajectory of C . The internal solution occurs internal to the trajectory, whereas the boundary solution occurs at its end; the rightmost of the two is where e occurs. Bottom: Variation of the internal (dotted) and boundary (dash-dotted) solutions with l; the maximum of two gives the global optimum e (solid). The two curves intersect at the critical value l c r 1.20711. Furthermore, the internal solution disappears at l c r 2 1.29904 together with the local minimum solution (dashed gray).

Fig. 6.

The optimal solution e ( l ) in the bookmark folding problem. Top: Three pages of increasing l are shown with a sequence of folds. The folds have a uniform distribution of a [ l / 2 , l ], with b chosen to have O overlap with the top edge. The thick black curves track the trajectory of C . The internal solution occurs internal to the trajectory, whereas the boundary solution occurs at its end; the rightmost of the two is where e occurs. Bottom: Variation of the internal (dotted) and boundary (dash-dotted) solutions with l; the maximum of two gives the global optimum e (solid). The two curves intersect at the critical value l c r 1.20711. Furthermore, the internal solution disappears at l c r 2 1.29904 together with the local minimum solution (dashed gray).

Close modal

In the top part of Fig. 6, three pages with increasing aspect ratios are each folded in many different ways, and the trajectories of the lower-right page corner ( C in Fig. 2) are recorded by thick black curves. The internal and boundary solutions are marked in the leftmost example. As l increases, the optimal solution clearly changes from the internal to the boundary solution.

Overall, the solution to e as a function of the page aspect ratio l has three distinct behaviors, which is surprising for such a simple problem. A graphical summary is found at the bottom of Fig. 6, which plots the internal solution e ̂ i n ( l ) and the boundary solution e ̂ b d ( l ); the maximum of the two gives e . The two curves intersect at the critical point l c r 1.20711, which signals an abrupt change of the optimal fold state. Figure 6 also plots the excess length given by the local minima (gray curve). This curve is the dashed curve in Fig. 5 but plotted as a function of l instead of a(l). It intersects with the internal solution at l = l c r 2, forming a cusp where the two solutions annihilate. The jump of the optimal solution from the internal solution branch to the boundary solution branch as well as the spontaneous disappearance of extrema in pair are all well-known manifestations of bifurcations. Interestingly, the behavior of e can be loosely likened to the (negative of) free energy in first-order phase transitions.

In this paper, we investigated an everyday physics question, “how to optimally fold a rectangular notebook page into a bookmark?” so that the folded page remains most visible but stays within one side of the notebook bind. The solution to this problem is unexpectedly rich given that the problem is governed by a single control parameter, the page's aspect ratio. In particular, this system exhibits a bifurcating behavior found in dynamical systems.

While related to origami, the problem has a fixed edge and makes a single fold. To distinguish the two, we may name ours as the “1-fold bookmark problem.” A few extensions can be readily proposed for future studies.

  1. Use, for example, circular or even concave pages. For a circular page, one may fix its topmost point, as there is no straight edge. Regardless of its shape, the page can be contained within a bounding rectangle. Or more intuitively, it can be glued on top of a larger rectangular “bounding page.” The maximum excess length of the bounding page is known from the current work and is always greater than that of the bounded page. Thus, the current work provides useful upper limits even for folding pages of arbitrary shapes.

  2. Extend to higher dimensions. This is challenging to represent. While the page in our problem is two-dimensional, during folding, it temporarily relies on the third dimension. Similarly, folding a three-dimensional shape relies on the fourth dimension.

  3. Allow more than one fold to tackle the n-fold bookmark problem. Because each fold yields a more complex shape, the geometric complexity increases rapidly. However, a quick look at the “2-fold” problem of a square page shows that a simple 2-fold solution easily beats our optimal 1-fold solution: first fold the square along its diagonal to have an isosceles right triangle and then fold the triangle to align its longer edge with its top edge. This gives e = 2 1 0.41421, while our optimal solution from Eq. (8) is only

Although the simple 2-fold solution outperforms the best 1-fold solution, it is unclear whether it is the optimal way to fold a square page twice into a bookmark, and whether its performance generalizes to pages of different aspect ratios.

The author would like to thank the reviewers for their comments and Dr. Rodolfo R. Rosales for the discussions.

The author has no conflicts to disclose.

As an example, we list the steps of finding lcr using mathematica:

  1. Find the roots by Solve[a4̂==(2a − l)l, a], which gives four roots.

  2. To pick the right root, one can substitute some test values of l. Suppose that root is named aSol.

  3. Substitute aSol to have the internal solution eIn =−2+l/a+Sqrt[(2a−l)l]/.{a->aSol}.

  4. Find lcr numerically by providing a guess of its value FindRoot[eIn == l − 1, {l, 1.2}].

1.
Fabian
Haas
and
Robin J.
Wootton
, “
Two basic mechanisms in insect wing folding
,”
Proc. R. Soc. London. Ser. B: Biol. Sci.
263
(
1377
),
1651
1658
(
1996
).
2.
Jakob A.
Faber
,
Andres F.
Arrieta
, and
André R.
Studart
, “
Bioinspired spring origami
,”
Science
359
,
1386
1391
(
2018
).
3.
Roger C.
Alperin
and
Robert J.
Lang
, “
One-, two-, and multi-fold origami axioms
,” in
Origami 4
, edited by
J.
Robert Lang
(
A K Peters/CRC Press
,
New York
,
2009
), Chap. 32, p.
371
.
4.
Robert J.
Lang
,
Origami Design Secrets: Mathematical Methods for an Ancient Art
(
AK Peters/CRC Press
,
Boca Raton, FL
,
2011
).
5.
Kazuo
Haga
,
Origamics: Mathematical Explorations Through Paper Folding
(
World Scientific
,
Singapore
,
2008
).
6.
Thomas
Hull
,
Project Origami: Activities for Exploring Mathematics
, 2nd ed. (
CRC Press
,
Boca Raton, FL
,
2012
).
7.
John
Jumper
,
Richard
Evans
,
Alexander
Pritzel
,
Tim
Green
,
Michael
Figurnov
,
Kathryn
Tunyasuvunakool
,
Olaf
Ronneberger
,
Russ
Bates
,
Augustin
Žídek
,
Alex
Bridgland
et al, “
High accuracy protein structure prediction using deep learning
,” in
Fourteenth Critical Assessment of Techniques for Protein Structure Prediction (Abstract Book)
(May–September
2020
), pp. 22 and 24.
8.
Timothy G.
Leong
,
Paul A.
Lester
,
Travis L.
Koh
,
Emma K.
Call
, and
David H.
Gracias
, “
Surface tension-driven self-folding polyhedra
,”
Langmuir
23
(
17
),
8747
8751
(
2007
).
9.
Edwin A.
Peraza Hernandez
,
Darren J.
Hartl
, and
Dimitris C.
Lagoudas
,
Active Origami: Modeling, Design, and Applications
(
Springer
,
Berlin
,
2018
).
10.
Daniel L.
Blair
and
Arshad
Kudrolli
, “
Geometry of crumpled paper
,”
Phys. Rev. Lett.
94
(
16
),
166107
(
2005
).
11.
Jovana
Andrejevic
,
Lisa M.
Lee
,
Shmuel M.
Rubinstein
, and
Chris H.
Rycroft
, “
A model for the fragmentation kinetics of crumpled thin sheets
,”
Nat. Commun.
12
(
1
),
1470
(
2021
).
Published open access through an agreement with Massachusetts Institute of Technology