Xem mã nguồn của Writhe

Trong lý thuyết nút thắt, có nhiều cách khác nhau để hiểu đại lượng '''writhe''', hay <math>\operatorname{Wr}</math>. Một mặt, writhe là một đại lượng của link, có hướng và nhận giá trị nguyên. Một cách hiểu khác của writhe là số thực mô tả độ "xoắn" của nút thắt theo nghĩa toán học (hay bất kỳ đường cong đơn khép kín nào) trong không gian ba chiều. Trong cả hai trường hợp, nó là một đại lượng hình học, tức là không được bảo toàn dưới các phép biến hình bảo toàn cấu trúc tô pô.

== Writhe của biểu đồ link ==
Trong lý thuyết nút thắt, writhe là tổng số giao điểm dương trừ đi tổng số giao điểm âm.

Mỗi một nút thành phần trong link được gán một hướng (trái - phải) xuyên suốt. Với mỗi giao điểm ta gặp khi đi theo hướng được gán cho, nếu sợi phía dưới đi từ phải qua trái, giao điểm ấy là giao điểm dương, và ngược lại. Ta có thể nhớ quy tắc này bằng một biến thể của quy tắc bàn tay phải.
{| style="margin:1em auto;"
|[[File:Knot-crossing-plus.svg|liên_kết=https://en.wikipedia.org/wiki/File:Knot-crossing-plus.svg|64x64px]]
| width="64px;" |
|[[File:Knot-crossing-minus.svg|liên_kết=https://en.wikipedia.org/wiki/File:Knot-crossing-minus.svg|64x64px]]
|-
|Giao điểm
dương
|
|Giao điểm
âm
|}
Với một biểu đồ nút bất kì, quy tắc bàn tay phải cho ra kết quả không đổi khi ta gán cho biểu đồ nút ấy hai hướng khác nhau, từ đây ta suy ra writhe được định nghĩa tốt trên các biểu đồ nút vô hướng.
[[File:Reidemeister_move_1.svg|liên_kết=https://en.wikipedia.org/wiki/File:Reidemeister_move_1.svg|nhỏ|150x150px|Một [[thao tác Reidemeister]] loại I tăng hoặc giảm writhe bởi 1.]]
[[Thao tác Reidemeister]] loại II và III không ảnh hưởng đến writhe, nhưng thao tác loại I tăng hoặc giảm writhe 1 đơn vị. Điều này có nghĩa là writhe của một nút không phải là một bất biến đẳng cấu của một nút, mà chỉ đóng vai trò đó với biểu đồ nút mà thôi. Bằng một chuỗi các thao tác loại I, ta có thể thay đổi writhe của một biểu đồ nút thành một số nguyên dương bất kì.

== Writhe của một đường cong khép kín ==
Trong một bài báo năm 1961, Gheorghe Călugăreanu chứng minh định lý sau: cho một ruy băng trong <math>\R^3</math>, gọi <math>\operatorname{Lk}</math> là số linking của các nút nằm ở rìa ruy băng đó, và <math>\operatorname{Tw}</math> là tổng độ xoắn. Khi đó, hiệu <math>\operatorname{Lk}-\operatorname{Tw}</math> chỉ phụ thuộc vào đường cong định nghĩa của ruy băng, và

: <math>\operatorname{Wr}=\operatorname{Lk}-\operatorname{Tw}</math>.

Trong một bài báo từ năm 1959, Călugăreanu cũng tìm ra cách để tính writhe <math>\operatorname{Wr}</math> thông qua một tích phân. Cho <math>C</math> là một đường cong khép kín, đơn, trơn, và cho <math>\mathbf{r}_{1}</math> và <math>\mathbf{r}_{2}</math> là hai điểm trên <math>C</math>. Khi đó writhe của đường cong bằng với tích phân Gauss

: <math>
\operatorname{Wr}=\frac{1}{4\pi}\int_{C}\int_{C}d\mathbf{r}_{1}\times d\mathbf{r}_{2}\cdot\frac{\mathbf{r}_{1}-\mathbf{r}_{2}}{\left|\mathbf{r}_{1}-\mathbf{r}_{2}\right|^{3}}
</math>.

== Xấp xỉ tích phân Gauss cho writhe của một đường cong trong không gian ==
Vì writhe của một đường cong trong không gian được định nghĩa qua một tích phân đôi, ta có thể tính ước tính được giá trị của nó bằng cách xấp xỉ đường cong bằng một chuỗi gồm <math>N</math> đoạn thẳng. Vốn được Michael Levitt phát minh để mô tả quá trình gập protein,  quy trình dưới đây về sau được Konstantin Klenin và Jörg Langowski dùng cho các ADN siêu xoắn.

: <math>
\operatorname{Wr}=\sum_{i=1}^{N}\sum_{j=1}^{N}\frac{\Omega_{ij}}{4\pi}=2\sum_{i=2}^{N}\sum_{j<i}\frac{\Omega_{ij}}{4\pi}</math>,

trong đó <math>\Omega_{ij}/{4\pi}</math> là tính toán chính xác của tính phân đôi trên đoạn i và j; để ý rằng <math>\Omega_{ij}=\Omega_{ji}</math> và <math>\Omega_{i,i+1}=\Omega_{ii}=0</math>

Để tính toán được <math>\Omega_{ij}/{4\pi}</math> cho các đoạn được đánh số <math>i</math> và <math>j</math>, đánh số các điểm cuối của đoạn là 1, 2, 3 và 4. Gọi <math>r_{pq}</math> là vectơ bắt đầu từ điểm cuối <math>p</math> và kết thúc ở điểm cuối <math>q</math>. Định nghĩa các đại lượng sau:

: <math>
n_{1}=\frac{r_{13}\times r_{14}}{\left|r_{13}\times r_{14}\right|},\; n_{2}=\frac{r_{14}\times r_{24}}{\left|r_{14}\times r_{24}\right|},\; n_{3}=\frac{r_{24}\times r_{23}}{\left|r_{24}\times r_{23}\right|},\; n_{4}=\frac{r_{23}\times r_{13}}{\left|r_{23}\times r_{13}\right|}
</math>

Sau đó ta tính

: <math>
\Omega^{*}=\arcsin\left(n_{1}\cdot n_{2}\right)+\arcsin\left(n_{2}\cdot n_{3}\right)+\arcsin\left(n_{3}\cdot n_{4}\right)+\arcsin\left(n_{4}\cdot n_{1}\right).
</math>

Cuối cùng, ta điều chỉnh dấu và chia cho <math>4\pi</math> để thu được

: <math>
\frac{\Omega}{4\pi}=\frac{\Omega^{*}}{4\pi}\text{sign}\left(\left(r_{34}\times r_{12}\right)\cdot r_{13}\right).
</math>

Thêm vào đó, các phương pháp khác để tính writhe có thể được mô tả hoàn toàn bằng thuật toán và toán học, một vài trong số đó vượt trội hơn phương pháp trên (phương pháp trên có độ phức tạp đa thức bậc hai, các phương pháp khác có thể có độ phức tạp bậc nhất)
[[File:Simulation_of_an_elastic_rod_relieving_torsional_stress_by_forming_coils.ogv|nhỏ|Một mô phỏng gậy dẻo giảm tải sức căng xoắn bằng cách cuộn lại.]]

== Ứng dụng trong ngành tô pô ADN ==
Vì ADN cuộn lại khi bị xoắn, các nhà toán sinh học sử dụng đại lượng writhe để mô tả độ biến dạng của một ADN dưới sức căng xoắn. Tổng quát hơn, hiện tượng tạo cuộn thông qua writhe được gọi là ADN siêu xoắn. HIện tượng này khá phổ biến, thực tế thì trong hầu như mọi cơ quan sống, DNA bị siêu xoắn âm.

Mọi gậy dẻo, không chỉ mỗi DNA, giảm tải sức căng xoắn bằng cách cuộn lại. Việc cuộn vừa tháo xoắn vừa bẻ cong cây gậy. F. Brock Fuller đã chứng minh toán học cách "giảm thế năng đàn hồi sinh ra do sự xoắn địa phương của một gậy bằng cách cuộn đường cong trung tâm của gậy sao cho số writhe tăng".