Население мира N (Т) будет описываться функцией от времени Т, определяющей состояние демографической системы Земли. Тогда параметром порядка — ведущей переменной, подчиняющей все остальные переменные, — станет полное число людей N. Таким образом, в этом приближенном функциональном соотношении не учитывается ни распределение населения по нашей планете, ни его экономическое и возрастное состояние или расовый и национальный состав. Сам процесс роста также будет рассматриваться на значительном интервале времени T — большом числе поколений. Иными словами, мы будем рассматривать усредненные значения переменных и усредненные функции. Этим, в частности, вносится в уравнения память о прошлом, определяемая временем усреднения переменных.

Такое выделение главных переменных N и T и их усреднение характерно для системного подхода. Оно получило развитие в синергетике и лежит в основе асимптотических методов, разработанных для решения задач большой сложности, появляющихся при рассмотрении систем со многими степенями свободы. Существенно то, что эти переменные, которые представляют все социально значимые факторы о возрасте и поле, образованию и развитию, доходам и т. д., описываются статистическими распределениями. Поэтому, когда рассматриваются такие многофакторные проблемы, то можно полагать, что в известных пределах развитие системы статистически стационарно и потому происходит динамически самоподобно. Это сильное предположение означает, что остаются неизменными пропорции между относительными изменениями времени и населения.

Смысл этой основной гипотезы автомодельности состоит в том, что утверждается постоянство относительной скорости изменения системы аналогично принципу инерции. В таком случае можно показать, что такой самоподобный рост должен описываться степенной функцией без характерного параметра, такого как масштаб времени. Такие процессы обладают масштабной инвариантностью — скейлингом — аналогично развитой турбулентности в потоке жидкости. Эти понятия мало знакомы историкам и обществоведам, однако они должны помочь в расширении тех образов, которыми мы описываем исторический процесс.

Выше (см. с. 29) изложено, как данные демографии приводят к формуле:

Это выражение как степенная функция обладает масштабной инвариантностью — отсутствием собственного масштаба времени, свойством, открытым еще Эйлером и указывающим на автомодельность роста. В нашей задаче о росте населения эта формула является лишь первым приближением. Как асимптотическое выражение оно ограничено областью применения, и задача теории в первую очередь состоит в установлении этих пределов как вблизи особенности, когда эта функция устремляется в бесконечность, так и в далеком прошлом, когда ее уменьшение происходит слишком медленно. Иными словами, асимптотика ограничена в прошлом нулем и полюсом в настоящее время.

Чтобы описать переход, следует учесть время, характеризующее внутренние процессы, определяемое продолжительностью жизни человека и его репродуктивной деятельности — тех факторов, которые при прохождении через демографический переход ограничивают скорость роста по мере приближения к моменту, когда скорость роста приближается к своему пределу. Для этого следует обратиться к выражению для скорости роста в зависимости от времени, продифференцировав (1):

и затем ввести в это расходящееся выражение характерное время τ, ограничивающее скорость роста:

Этот прием может показаться произвольным шагом, однако полученное выражение очень хорошо описывает глобальный демографический переход. Так, мы обратились к методам в теоретической физике, которые развиты для регуляризации расходимости, появившейся при анализе демографического перехода.

Интегрируя (3), получим выражение для описания перехода:

При обращении к последним данным демографии (см. рис. 15) были получены уточненные значения постоянных, что учтено в последующих вычислениях:

Из-за введения конечного τ полюс в Т₁ сдвигается к новому значению Т₁ = 1995 г., которое и принято при расчетах, описывающих как демографический переход, так и рост населения мира за пределы Т₁ в выражении (4) (см. табл. 1).

В недалеком прошлом выражение (4) асимптотически непосредственно переходит в автомодельный гиперболический рост (1). Однако применительно к очень далекому прошлому скорость роста должна быть ограничена снизу. Этого предположения достаточно для того, чтобы приписать далекому прошлому линейный рост, при котором в первом приближении скорость роста не может быть меньше появления одного гоминида за характерное время τ, пока численность населения не достигает порядка 100 тыс. В популяционной генетике это число К характерно для численности стабильного вида, биологически подобного человеку, и именно с него 1,6 млн лет начинается квадратичный рост, который с тех пор становится доминирующим до эпохи перехода.

Величина К определяет не только масштаб численности человечества в начальную эпоху роста, но и дает оценку численности когерентной группы людей или племени — самодостаточной единицы населения. Как большой параметр задачи постоянная К определяет все соотношения между населением и длительностью процессов роста, а значительная величина константы К приводит к высокой эффективности асимптотических решений. В результате скорость роста населения Земли определяется нелинейным автономным дифференциальным уравнением:

где время t = Т/τ выражено в единицах времени τ, и в решениях уравнения (6) отсчитывается от момента прохождения через демографический переход Т₁ Это характерное время одинаково для фазовых переходов в прошлом и настоящем. Сделанные предположения упрощают задачу, сводя все к одной переменной N (Т), рост которой зависит от состояния системы N в момент времени Т. Формула роста (6) выражает природу того коллективного нелинейного взаимодействия, которое ответственно за рост человечества в эпоху его взрывного развития между двумя сингулярностями. В этом уравнении для усредненных переменных T и N скорость роста приравнена к развитию, которое равно квадрату численности населения мира, как выражение меры системной сложности населения планеты.

Население планеты также можно рассматривать как результат парного взаимодействия N человек или как некое эффективное поле, феноменологически определяющее рост. Полное решение должно описывать рост человечества в течение трех эпох. Первая эпоха А — антропогенеза — начинается с линейного роста с указанной выше минимальной скоростью. Когда население достигает величины порядка ста тысяч, наступает эпоха В — взрывного роста — со скоростью роста, пропорциональной квадрату населения Земли. Начиная с этого момента, человек заселяет всю планету.

Когда скорость квадратичного роста достигла своего предела при удвоении за характерное время τ, наступил кризис мирового демографического развития, переход в эпоху С — стабилизации населения мира в рамках приближений теории.

1 — абсолютный прирост населения, усредненный за декаду, млн; 2 — относительный прирост,% в год (данные ООН)

В результате на основании (3) максимальная абсолютная скорость глобального роста во время демографического перехода равна:

при относительном росте:

достигнутом в 1995 г., что согласуется с данными ООН, но дает несколько меньшее значение для абсолютной скорости роста при сравнении с табл. 1 (см. рис. 18).

Население нашей планеты в этот критический момент перехода Т₁ равно:

На основе этих выражений легко определить предел N_∞, в два раза больший, чем N₁ к которому в эпоху С асимптотически стремится население Земли:

В рамках сделанных предположений это число представляет верхнюю оценку населения Земли в предвидимом будущем. Таким образом, глобальное взаимодействие приводит к ускорению и синхронизации процессов и на заключительной стадии демографического перехода — к сужению перехода и тем самым к снижению предела для населения Земли. Этот вывод находится в согласии с эмпирическими наблюдениями демографов. Рассмотрение N (Т) как аналитической функции указывает на асимптотическое поведение при T → ∞, когда N → N_∞, в предположении об отсутствии особенностей — полюсов или нулей — в обозримом будущем.

Начальный линейный рост дает оценку времени для эпохи антропогенеза — критической сингулярности в начале предыстории человечества, которая случилась:

если использовать известное значение N₁ и одно и то же значение τ = 45 лет для сингулярности в далеком прошлом и в настоящем. Несмотря на сделанные упрощения, данная оценка вполне согласуется с оценками времени Т₀ в антропологии.

Интересно определить полное число людей, живших на Земле. Если переставить переменные в (6) и проинтегрировать:

то получим число людей, живших от Т₀ до нашего времени Т₁ В оценках других авторов длительность поколения принята равной 20 годам, что ведет к оценке Р_0,1 = 106 млрд [10]. Поэтому необходимо введение в (12) множителя 45/20 = 2,25:

Таким образом, в течение каждого из ln К = 11,0 выделенных периодов жило по 2,25 K² = 8 млрд людей. Это число является инвариантным для числа людей, живших в экспоненциально сокращающихся циклах, а ln K указывает на число циклов.

Циклы можно получить, обобщая решение (6) в область комплексных переменных или же просуммировав экспоненциально сокращающиеся циклы:

где θ = |ln t| — номер цикла, определить длительность развития при К >> 1:

и сравнить ее с (11), где длительность равна Т₁ - Т₀ = π/2·Kτ = 1,571. В первом случае рост суммируется по гиперболической траектории, во втором — по (4) — N = K tan t/K.

Демографические циклы определяют периодичность развития всего человечества за 4–5 млн лет, включая проходящий по гиперболическому закону рост от конца антропогенеза до наших дней.

Для дальнейшего обзора результатов перейдем к переменной n = N/K:

когда мерой численности становится К. Тогда уравнения для роста приобретают симметричный вид и видно сопряжение переменных n и t. Смена зависимой переменной в (16а) и (16d) видна при прохождении перехода, когда n становится независимой переменной вместо времени t, что выражено в уравнении роста (3).

Рост населения можно иллюстрировать геометрическим построением функции тангенса:

где угол Δφ = τ отображает течение времени, а приращение населения ΔN = 1 (рис. 16).

Линейный рост будет продолжаться до φ_A,B = Кτ = 1 и N_B = tan 1 в точке В на касательной АС. Дальнейший рост N = К(π/2 — φ)^-1 будет проходить по гиперболе, при которой время асимптотически стремится к π/2, а население достигнет значения N_c = К². Когда система приближается к моменту особенности, то от уравнения (16а) следует переходить к уравнению (16d), чтобы описать рост при прохождении особенности в течение эпохи С. Построение показывает, что после перехода от линейного к гиперболическому росту на эпоху В остается в два раза меньше времени, чем на начальную эпоху А. Вывод этого соотношения для всей эпохи В (см. рис. 19) построен при К = 7, когда время от Т₀ до Т₁ разделено на 11 интервалов, и поскольку к/2 = 1/7, то N_c = К² =49. Однако даже при таком малом значения К, когда In K = In 1,95 дает хорошую оценку для числа демографических циклов, 1 + In К ≈ 3. Таким образом нулевой цикл антропогенеза продолжался 7 единиц времени, первый цикл длился 3 и последний — одну единицу времени. Это построение показывает, как дискретность времени и населения приводит к появлению периодичности роста, выраженной в демографических циклах как главных эпохах развития человечества.

Линейный рост описывает поведение системы вблизи начальной сингулярности роста, начинающейся с N₀ = 1 и положительных значений N. Далее следует рост по гиперболе и в конце — сингулярность демографического взрыва. Построение, когда переменные n и t при прохождении перехода меняются местами, мы оставляем читателю.

На рис. 18 показаны функции, описывающие рост системы при К = 1, которые появляются при построении решения, начинающегося с сингулярности в эпоху А, переходящего затем в эпоху В гиперболического роста и завершающегося эпохой С. Асимптотический переход решений, описывающий рост в начале развития и на его конечном участке, получим, обратившись к рядам для функции cot (t/K) и cot^-1(t/K):

Эти функции пересекаются в точке А, посередине роста при логарифмическом представлении между временем T₀ и T₁ соответствующей наступлению неолита:

под углом 2/(3K) практически гладко при больших значениях К.

Очевидно, что решение можно строить, отсчитывая время от T₀ — от эпохи антропогенеза А при t₀ = 0. Тогда, исключив t из (15с), получим одно автономное дифференциальное уравнение, описывающее рост в зависимости от состояния системы, которое определяется населением Земли и где последний член добавлен с тем, чтобы рост в эпоху А никогда не был меньше одного гоминида при Δt = τ.

Интегрируя (20) при значениях K > 1 и начальных условиях t₀ = n₀ = 0, получим решение:

Это решение показывает симметрию переменных N и T — населения и времени. Для развития в течение эпохи В вдали от особенностей роста это выражено в (16в) и следует из сложности причинных связей в рамках развитых представлений о нелинейной динамике глобальной системы населения нашей планеты.

Для того чтобы выяснить устойчивость развития, следует обратиться к уравнению роста человечества (20). На основании (15) в линейном приближении устойчивость роста к возмущениям

определит показатель Ляпунова λ развития неустойчивости в системе населения:

По этому критерию при λ > О движение неустойчиво до перехода. Только после него развитие системы становится асимптотически устойчивым и впредь таким и остается. Отметим, что в этих решениях значение констант роста К и τ не эволюционируют. Более полное определение устойчивости потребует введения распределений для n и обращения к методам статистической физики при обобщении развитой выше модели.

Рис. 21. Переходные процессы и устойчивость роста в линейном приближении

1 — логистический переход ν = 1/1+е^-r; 2 — демографический переход η = 1/π соt^-1 T и λ (ν).

При гиперболическом росте мгновенное значение экспоненциального роста равно древности,

что и определяет скорость процессов развития в момент времени Т.

В гиперболической хронологии мгновенный экспоненциальный масштаб времени роста линейной неустойчивости по Ляпунову зависит от древности и до демографического перехода равен удвоенному времени роста неустойчивости:

Однако наличие выделенных антропологами и историками демографических циклов указывает на глобальную устойчивость с малыми отклонениями системы от предельной траектории роста.

Наконец, из (15) следует, что после каждого цикла до демографического перехода остается приблизительно половина времени длительности цикла:

что вполне подтверждается данными истории и антропологии (см. табл. 2, с. 74).

Асимптотические решения для нелинейных задач приобретают, как заметил Я. Б. Зельдович, то же значение, что частные решения для линейных задач, где действует принцип суперпозиции. В нашей задаче автомодельные решения имеют асимптотический характер, при котором значение некоторых параметров оказываются несущественным, подобно тому как в частных решениях линейных задач происходит вырождение по некоторым параметрам. Так, частота колебаний в линейной задаче не зависит от их амплитуды.

В некоторых случаях независимость от параметра позволяет от уравнений в частных производных перейти к обыкновенным дифференциальным уравнениям. Именно это происходит в задаче о росте человечества, когда в нелинейном уравнении (6) можно пренебречь пространственным распределением населения, поскольку в первом приближении перемещение на конечной по размерам Земле — миграция населения — не влияет на само число людей.

В заключение заметим, что изложенная теория рассматривает задачу о росте и развитии демографической системы в асимптотическом приближении. Из этого следует, что в этом приближении ресурсы не влияют на глобальное взрывное развитие. Поэтому ограничение роста и переход к стабилизации населения мира обязаны внутренним процессам, выраженным в принципе демографического императива, и не подчиняется внешним, в первую очередь ресурсным факторам.

Здесь, несомненно, открывается обширное поле для исследования этой проблемы во всеоружии современных методов статистической физики, теории фазовых переходов и неравновесных систем. Более того, эта задача во многом подобна задачам космологии. Но она гораздо ближе к нам, поскольку имеет огромное значение для фундаментального понимания истории и экономики. Выразим поэтому уверенность, что эта проблематика приведет к развитию своих методов и в теоретической физике.