Конструкция Максвелла
В термодинамическом равновесии необходимым условием стабильности является давление, , не увеличивается с увеличением объема или молярного объема, ; математически это выражается как , где это температура. [1] Это основное требование стабильности, как и аналогичные требования для других сопряженных пар переменных , нарушается в аналитических моделях фазовых переходов первого рода . Самый известный случай — уравнение Ван дер Ваальса . [2] [3] где являются размерными константами. Это нарушение не является дефектом, а является причиной наблюдаемого разрыва свойств, отличающих жидкость от пара, и определяет фазовый переход первого рода.

На рис. 1 показана изотерма, построенная для , как непрерывно дифференцируемую сплошную черную, пунктирную черную и пунктирную серую кривые. Убывающая часть кривой справа от точки С на рис. 1 описывает газ, а убывающая часть слева от точки Е описывает жидкость. Эти две части разделены областью между локальным минимумом и локальным максимумом на кривой с положительным наклоном, что нарушает критерий устойчивости. Этот математический критерий выражает физическое состояние, которое Эпштейн [4] описано следующим образом:
«Очевидно, что эта средняя часть, отмеченная на наших кривых пунктиром [здесь пунктирная линия на рис.1], не может иметь никакой физической реальности. В самом деле, представим себе жидкость в состоянии, соответствующем этой части кривой, заключенной в проводящий вертикальный цилиндр, вершина которого образована поршнем. Поршень может скользить вверх и вниз в цилиндре, и мы приложим к нему нагрузку, точно уравновешивающую давление газа. Если мы снимем с поршня небольшой вес, этого не произойдет. больше не будет находиться в равновесии, и он начнет двигаться вверх. Однако по мере его движения объем газа увеличивается, а вместе с ним и его давление, действующее на поршень, увеличивается, поэтому поршень будет продолжать двигаться вверх. двигаться, и газ будет расширяться, пока не достигнет состояния, представленного максимумом изотермы. И наоборот, если мы добавим хоть немного к нагрузке сбалансированного поршня, газ схлопнется до состояния, соответствующего минимуму изотермы. ."
Эта ситуация аналогична телу, точно сбалансированному на вершине гладкой поверхности, которое при малейшем возмущении отойдет от своего положения равновесия и будет продолжать движение до тех пор, пока не достигнет локального минимума. В описанном виде такие состояния динамически неустойчивы и поэтому не наблюдаются. Разрыв является предшественником фактического фазового перехода от жидкости к пару. Точки Е и С , где , которые ограничивают максимально возможное состояние жидкости и наименьшее возможное состояние пара, называются спинодальными точками. Их локус образует спинодальную кривую, ограничивающую область, в которой не могут существовать однородные устойчивые состояния.
Эксперименты показывают, что если объем сосуда, содержащего фиксированное количество жидкости, нагревается и расширяется при постоянной температуре и определенном давлении, , пар, (обозначены точками в точках и на рис. 1) образуются пузырьки, поэтому жидкость больше не является однородной, а превращается в гетерогенную смесь кипящей жидкости и конденсирующегося пара. Гравитация отделяет кипящую (насыщенную) жидкость, , от менее плотного конденсирующегося (насыщенного) пара, которые сосуществуют при одинаковых температуре и давлении насыщения. По мере нагревания количество пара, увеличивается и жидкость, , уменьшается. Все это время давление, и температура, , остаются постоянными, а объем увеличивается. В этой ситуации молярный объем смеси представляет собой средневзвешенное значение ее компонентов. где , мольная доля пара, , непрерывно возрастает; однако молярный объем самого вещества имеет только максимально возможное стабильное значение для его жидкого состояния и наименьшее возможное стабильное значение для его парообразного состояния при данном . Повторим, хотя мольный объем смеси непрерывно переходит от к (обозначено пунктирной линией на рис. 1), подстилающая жидкость имеет разрыв по этому свойству, как и другие. Это уравнение состояния смеси называется правилом рычага . [5] [6] [7]
Пунктирные части кривой на рис. 1 — метастабильные состояния . В течение многих лет такие состояния были академической диковинкой; Каллен [8] привел в качестве примера,
вода, охлажденная ниже 0°С при давлении 1 атм. Удар по стакану с водой в этом состоянии вызывает внезапную резкую кристаллизацию системы.
Однако исследования теплопередачи при кипении ясно показали, что метастабильные состояния возникают обычно как неотъемлемая часть этого процесса. В нем температура поверхности нагрева выше температуры насыщения, часто значительно, поэтому прилегающая жидкость должна быть перегрета. [9] Кроме того, появление устройств, работающих с очень высокими тепловыми потоками, вызвало интерес к метастабильным состояниям и связанным с ними термодинамическим свойствам, в частности к состояниям перегретой жидкости. [10] Более того, тот факт, что они предсказываются уравнением Ван-дер-Ваальса и кубическими уравнениями в целом, является убедительным свидетельством его эффективности при описании фазовых переходов; Зоммерфельд описал это следующим образом: [11]
Весьма примечательно, что теория Ван-дер-Ваальса способна предсказать, по крайней мере качественно, существование нестабильных (называемых здесь метастабильными) состояний вдоль ветвей AA` или BB` (BC и FE на рис. 1). здесь].
Правило равной площади
[ редактировать ]Прерывистость в и другие свойства, например внутренняя энергия, и энтропия, , вещества, называется фазовым переходом первого рода. [12] [13] Чтобы указать уникальное экспериментально наблюдаемое давление, , при котором это происходит, необходимо другое термодинамическое условие, поскольку из рис.1 ясно, что это может иметь место при любом давлении в диапазоне . Такое условие было впервые сформулировано Максвеллом в остроумном термодинамическом рассуждении на лекции, которую он прочитал в Британском химическом обществе 18 февраля 1875 года. [14] (Рис.1, включая буквы BCDEF, представляет собой описанную им кривую):

«Участок кривой от С до Е представляет собой точки, которые по существу нестабильны и поэтому не могут быть реализованы.
Предположим теперь, что среда переходит из B в F по гипотетической кривой B C D E F в состоянии всегда однородном, а возвращается по прямому пути F B в виде смеси жидкости и пара. Поскольку температура повсюду была постоянной, никакое тепло не могло быть преобразовано в работу. Теперь тепло, преобразованное в работу, представлено превышением площади F D E над B C D . Следовательно, условием, определяющим максимальное давление пара при данной температуре, является то, что линия B F отсекает равные площади от кривой сверху и снизу».
Самый простой способ понять аргумент Максвелла — рассмотреть предложенный им цикл на плоскости температура — молярная энтропия. Каждый вводный текст по термодинамике [15] [16] представлен тот факт, что на такой плоскости площадь под любой кривой представляет собой теплоотдачу вещества на моль, положительную, идущую слева направо, и отрицательную, справа налево; более того, в циклическом процессе чистая передача тепла веществу представляет собой площадь, ограниченную замкнутой кривой цикла. Поскольку рассматриваемый им цикл состоит из двух серых пунктирных изотерм при одной и той же температуре: одна идет от B к F (через CD и E), а другая прямо обратно от F к B, эти две линии идентичны, только что пройдены в обеспечить регресс; заключена нулевая область, и, следовательно, .
Далее в тех же текстах площадь под этими кривыми, нанесенная на давление-молярный объем, см. рис. 1, описывается как работа, совершаемая веществом, положительная, идущая слева направо, и отрицательная, справа налево. Точно так же чистая работа, совершенная за цикл, представляет собой площадь, ограниченную замкнутой кривой. Поскольку первый закон термодинамики в частном случае цикла дает результат , для цикла, предложенного Максвеллом ; тогда, поскольку заключенная область равна I+II=0, см. рис.1, где I положительный, а II отрицательный, переходное давление должно быть таким, чтобы обе области были равны.
Записанное в виде математического уравнения с точки зрения работы, совершаемой в каждом процессе, это: Это уравнение вместе с уравнением состояния, записанным для каждого из состояний и представляют собой три уравнения для четырех переменных, , поэтому, учитывая любой из них, скажем , остальные три определены. Другими словами, существует уникальная ценность , а также и , при котором может произойти фазовый переход.
критерий Гиббса
[ редактировать ]В конце своей лекции, похвалив Ван дер Ваальса, назвав его работу «чрезвычайно гениальной диссертацией», Максвелл закончил ее словами:
Однако я не должен упустить из виду наиболее важный американский вклад в эту часть термодинамики, сделанный профессором Уиллардом Гиббсом из Йельского колледжа США, который дал нам удивительно простой и вполне удовлетворительный метод представления отношений между различными состояниями материи. с помощью модели. С помощью этой модели проблемы, которые долгое время сопротивлялись усилиям меня и других, могут быть решены сразу.
Это замечание оказалось пророческим, поскольку в 1876–1878 годах Гиббс опубликовал свою исчерпывающую работу по термодинамике. [17] в котором он показал, что термодинамическое равновесие гетерогенного вещества требует, чтобы, кроме механического равновесия (одинаковое давление для каждого компонента) и теплового равновесия (одинаковая температура для каждого компонента), должно существовать еще материальное равновесие (одинаковый химический потенциал для каждого компонента). В данном случае одно вещество и две фазы в дополнение к и , материальное равновесие требует (для частного случая одного вещества его химический потенциал представляет собой молярную функцию Гиббса, где ). [18] Это условие можно вывести с помощью простого физического аргумента следующим образом: энергия, необходимая для испарения моля, определяется вторым законом при постоянной температуре. , и из первого закона при постоянном давлении , затем приравнивание этих двух и перестановка дает результат, поскольку . Условия материального равновесия приводят к знаменитому правилу фаз Гиббса : , где количество веществ, количество фаз и количество независимых интенсивных переменных, необходимых для задания состояния. [19] [20] В случае одного вещества и двух фаз, обсуждаемых здесь, это дает: , экспериментально наблюдаемое число.
Сейчас – термодинамическая потенциальная функция , ее дифференциал равен [21]
Интегрирование этого при постоянной температуре дает здесь является константой интегрирования, но константа различна для каждой изотермы, поэтому она записана как функция . [22] Чтобы оценить нужно инвертировать чтобы получить . Однако природа явления фазового перехода заключается в том, что эта инверсия не уникальна; например, уравнение Ван дер Ваальса, написанное для является,

кубика с 1 или, в данном случае, 3 действительными корнями. Таким образом, как видно на рис. 2, существуют три кривые, состоящие из стабильного (показано сплошным черным), метастабильного (показано черным пунктиром) и нестабильного (показано пунктирным серым) состояний.
На самом деле эта цифра была получена не путем решения кубической задачи и интегрирования, а скорее было получено из его определения путем предварительного получения и , что легко сделать аналитически для уравнения Ван-дер-Ваальса, и построить его параметрически с помощью , с использованием в качестве параметра. Рассматривая только его устойчивые состояния является непрерывным с разрывными частными производными , и , в точке фазового перехода. В Эренфеста классификации фазовый переход первого рода означает разрыв первых частных производных в то время как фазовый переход второго рода будет включать разрывы вторых частных производных. [23]
Связь между критериями Гиббса и Максвелла
[ редактировать ]Вычисление интегрального выражения для заданное ранее между состояниями насыщенной жидкости и пара, и применение критерия Гиббса материального равновесия к этому процессу фазового перехода требует записи его как

Здесь интеграл разделен на три части с использованием трех действительных корней кубической степени, соответствующей жидкости: , нестабильный, и пар, , говорится соответственно. Эти интегралы лучше всего можно визуализировать, посмотрев на рис. 1 повернутым. против часовой стрелки в бумажной плоскости, затем о ось, так что появляется на левой ординате кривой, как показано на прилагаемом графике. С этой точки зрения функция явно многозначен; именно по этой причине ему требуются три реальные функции для описания его поведения между и . Теперь о разбиении среднего интеграла на два Первые два интеграла здесь представляют собой область I, а вторые два являются отрицательными для области II. Сумма двух площадей равна нулю, следовательно, их величины равны согласно этому критерию Гиббса. Это снова правило равных площадей Максвелла, конструкция Максвелла, и ее также можно показать аналитически. С , Интегрируя это для постоянной температуры из состояния к с условием Гиббса дает что является результатом Максвелла. Это правило равных площадей также можно вывести, используя свободную энергию Гельмгольца. [24] В любом случае конструкция Максвелла вытекает из условия Гиббса материального равновесия. Однако, хотя является более фундаментальным, оно более абстрактно, чем правило равной площади, которое понимается геометрически.
Общая касательная конструкция
[ редактировать ]
Другой метод определения точек сосуществования основан на принципе минимума потенциала Гельмгольца, который гласит, что в системе, находящейся в диатермическом контакте с тепловым резервуаром , и , а именно в состоянии равновесия потенциал Гельмгольца минимален. [25] Поскольку, как , молярная функция Гельмгольца также является потенциальной функцией, дифференциал которой равен: [26] этот принцип минимума приводит к условию устойчивости . [27] Это условие требует, чтобы в любом устойчивом состоянии системы функция строго выпукла , а именно, что в ее окрестности кривая лежит на ее касательной или выше нее. [28] Более того, для этих состояний обязательно выполняется и предыдущее условие устойчивости давления.
График этой функции для той же докритической изотермы уравнения ВДВ, что и на рис. 1 и 2 показаны на рис. 3. На этом рисунке изображена (пунктирная/сплошная) прямая линия, имеющая двойную (общую) касательную с кривой функции в B и F. Эта прямая линия: , с константа, которую можно записать как . Последнее равенство следует из соотношения , [29] вместе с . [18] Все это означает, что каждая точка прямой имеет одинаковые значения , в частности точки B и F, что дает условие Гиббса материального равновесия. а также равенство температуры и давления. [30] Поэтому эта конструкция эквивалентна как условиям Гиббса, так и конструкции Максвелла.
Эта конструкция, основанная на определенное ранее Гиббсом, [31] [32] первоначально использовался Ван дер Ваальсом (он называл его одновременно двойной и общей касательной), [33] потому что его можно легко расширить, включив в него бинарные смеси жидкостей , для которых изотерма , с переменная композиции образует поверхность, которая может иметь общую касательную плоскость. Впоследствии это стало популярным способом решения проблем с изменением фазы. [34] [35]
Приложение к уравнению Ван-дер-Ваальса
[ редактировать ]Из уравнения Ван-дер-Ваальса, примененного к насыщенной жидкости, и пар, , утверждает Эти два уравнения задают 4 переменные, поэтому их можно решить относительно с точки зрения . Это приводит к где , , и — характеристическое давление, мольный объём и температура, определяемые константами (обратите внимание, что ). Применение конструкции Максвелла к уравнению Ван-дер-Ваальса дает Эти три уравнения можно решить численно. Это было сделано с учетом значения либо или и представлены табличные результаты; [36] [37] однако уравнения также допускают аналитическое параметрическое решение, которое, по мнению Ленкнера, [38] был получен Гиббсом. Сам Ленкнер разработал простой и элегантный метод получения этого решения, исключив и из уравнений и записать их в терминах растянутой безразмерной плотности чисел, , это варьируется между и 0 как варьируется от к ; это производит Несмотря на трансцендентность, это уравнение имеет простое аналитическое параметрическое решение, полученное путем записи левой части уравнения, которая просто , как Затем и когда используется для устранения с правой стороны линейное уравнение для получено решение которого есть Соответственно, фундаментальная переменная, которая определяет все остальные в этом процессе фазового перехода, равна . Это решение проблемы насыщения легко расширить, чтобы охватить все ее переменные. где .


Значения всех других разрывов свойств на кривой насыщения также следуют из этого решения. [39]
Эти функции определяют кривую сосуществования, которая является местом расположения состояний насыщенной жидкости и насыщенного пара жидкости Ван-дер-Ваальса. На рис. 4 эта кривая изображена синим цветом вместе со спинодальной кривой черного цвета, рассчитанной по формуле где является параметром. Переменные, используемые при построении этих графиков, представляют собой приведенные (безразмерные) переменные, , , и где подстрочные величины представляют собой значения критической точки. Они определяются, , и в критической точке, [40] и являются измеримыми величинами. Отношения , , используются для перевода звездных величин в растворе в количества, использованные на рисунках. Кривая полностью согласуется с численными результатами, упомянутыми ранее. В области внутри спинодальной кривой в каждой точке имеется два состояния: одно стабильное и одно метастабильное: либо перегретая жидкость справа от синей кривой, либо недоохлажденный пар слева, тогда как вне спинодальной кривой имеется одно устойчивое состояние при каждая точка. На рис. 5 область под спинодальной кривой (черный штрихпунктир) не содержит однородных стабильных состояний, в то время как между кривыми сосуществования (красный штрихпунктир) и кривыми спинодали имеется по одному метастабильному состоянию в каждой точке, а вне кривой сосуществования имеется одно стабильное состояние. состояние в каждой точке. Два синих и два зеленых кружка обозначают состояния насыщенной жидкости и пара на соответствующих изотермах. Повсюду под кривой сосуществования также наблюдаются гетерогенные состояния, удовлетворяющие правилу рычага; однако они не являются однородными состояниями уравнения Ван-дер-Ваальса, поэтому их существование, обозначенное горизонтальными линиями, соединяющими точки насыщения на каждой докритической изотерме, не отображается. Кроме того, абсцисса на этом рисунке логарифмическая, а не линейная, чтобы показать большую часть области пара в целом. без чрезмерного сжатия жидкости и нестабильных областей при небольших ; однако это устройство искажает площади, поэтому две области I и II на рис.1 здесь не будут выглядеть равными.
В диапазоне параметров , монотонно убывает от и приближается к 0 как в пределе . Поэтому и в пределе , и . Поведение и следуют из уравнений. Оба эти свойства также монотонно убывают с и и приближаемся к 0 как и в пределе . Обратите внимание на то, что ; В этом пределе насыщенный пар Ван-дер-Ваальса является идеальным газом. Перефразируя Зоммерфельда, примечательно, что теория Ван дер Ваальса способна предсказать, что когда насыщенный пар ведет себя как идеальный газ; именно так ведут себя насыщенные пары реальных газов.
Кроме того, для жидкая спинодальная точка возникает при отрицательном давлении, а изотерма включен в рис. 4, чтобы проиллюстрировать этот момент. Это означает, что некоторая часть этих жидких метастабильных состояний находится в состоянии растяжения, и чем ниже температура, тем больше растягивающее напряжение. Хотя это кажется нелогичным, известно, что при некоторых обстоятельствах жидкости могут поддерживать напряжение. Тьен и Линхард [41] заметил это и написал:
Уравнение Ван-дер-Ваальса предсказывает, что при низких температурах жидкости выдерживают огромное напряжение — факт, который заставил некоторых авторов относиться к этому уравнению легкомысленно. В последние годы были проведены измерения, которые показали, что это совершенно верно. [42] Чистые жидкости, не содержащие растворенных газов, могут подвергаться напряжениям, превышающим по величине .
Это еще одна интересная особенность теории Ван дер Ваальса.
Примечания
[ редактировать ]- ^ Каллен, стр. 131-135, (1960}
- ^ Ван дер Ваальс, с. 174, (1988)
- ^ Эпштейн, стр. 9
- ^ Эпштейн, стр. 10.
- ^ Каллен, стр. 146-163.
- ^ Гудштейн, стр. 443-452.
- ^ Кондепуди, Дилип; Пригожин, Илья (31 декабря 2014 г.). Современная термодинамика: от тепловых двигателей к диссипативным структурам . Джон Уайли и сыновья. ISBN 9781118371817 .
- ^ Каллен, стр. 163.
- ^ Лиенхард IV и Линхард V, стр. 467-469.
- ^ Линхард, 169–187.
- ^ Зоммерфельд, с. 66
- ^ Каллен, стр. 146.
- ^ Гудштейн, стр. 443.
- ^ Максвелл, Дж. К., стр. 357-359.
- ^ Ван Вайлен и воскресенье, с. 211
- ^ Моран и Шапиро, HN с. 251
- ^ Гиббс, стр. 62-65
- ^ Jump up to: а б Каллен, стр. 99.
- ^ Гиббс, стр. 96-100.
- ^ Каллен, стр. 163-167.
- ^ Каллен, стр. 98-100.
- ^ Каллен, стр. 150.
- ^ Каллен, стр. 172-173.
- ^ Зоммерфельд, А., стр. 67-68, (1956)
- ^ Каллен, с. 105
- ^ Каллен, с. 120
- ^ Каллен, с. 135
- ^ Rectorys, с. 429
- ^ Каллен, стр.98
- ^ Ван дер Ваальс, с. 246
- ^ Де Бур, стр.8
- ^ Гиббс, стр. 87-90.
- ^ Ван дер Ваальс, стр. 246-247
- ^ Уэльс, Дэвид; Уэльс (2003). Энергетические ландшафты: приложения к кластерам, биомолекулам и стеклам . Издательство Кембриджского университета. п. 444. ИСБН 9780521814157 .
- ^ «Фазовые переходы первого рода и динамика спинодального распада» . www.mhkoepf.de . Проверено 12 ноября 2019 г.
- ^ Шамсундар и Линхард, 876-880 гг.
- ^ Барруфет и Юбенк, стр. 168-175.
- ^ Лекнер, с. 161
- ^ Джонстон, стр. 16-18.
- ^ Зоммерфельд, А., стр. 56-57.
- ^ Тьен и Линхард, стр.254.
- ^ Темперли, HNV, Поведение воды при гидростатическом растяжении: III , Proc. Физ. Соц. (Лондон), 59 , стр. 199, (1947)
Ссылки
[ редактировать ]- Барруфет, Массачусетс; Юбанк, ПТ (1989). «Обобщенные свойства насыщения чистых жидкостей посредством кубических уравнений состояния». Химико-технологическое образование . 23 (3): 168–175.
- Каллен, Х.Б. (1960). Термодинамика . Нью-Йорк: Джон Уайли и сыновья.
- ДеБоер, Дж. (1974). «Ван дер Ваальс в свое время и нынешнее вступительное слово возрождения». Физика . 73 : 1–27. дои : 10.1016/0031-8914(74)90223-7 .
- Эпштейн, PS (1937). Учебник термодинамики . Нью-Йорк: Джон Уайли и сыновья.
- Гиббс, Дж. В. (1928) [1876, 1878]. «О равновесии гетерогенных веществ». В Лонгли, WR; Ван Наме, Р.Г. (ред.). Собрание сочинений Дж. Уилларда Гиббса, том I по термодинамике (изд. 1948 г.). Нью-Хейвен: Издательство Йельского университета. стр. 55–353.
- Гудштейн, Д.Л. (1985) [1975]. Состояния материи . Нью-Йорк: Дувр.
- Джонстон, округ Колумбия (2014). Достижения в термодинамике жидкости Ван-дер-Ваальса . arXiv : 1402.1205 . дои : 10.1088/978-1-627-05532-1 . ISBN 978-1-627-05532-1 .
- Лекнер, Дж. (1982). «Параметрическое решение кривой сосуществования жидкости и пара Ван-дер-Ваальса». Являюсь. Дж. Физ . 50 (2): 161–163. дои : 10.1119/1.12877 .
- Линхард IV, Дж. Х .; Линхард В. (2019). Учебник по теплопередаче (Пятое изд.). Минеола, Нью-Йорк: Dover Publications Inc.
- Линхард, Дж. Х. (1986). «Свойства и поведение перегретых жидкостей». лат. Являюсь. J. Тепломассоперенос . 10 : 169–187.
- Максвелл, Дж. К. (1875). «О динамических доказательствах молекулярного строения тел». Природа . 11 (279): 357–359. дои : 10.1038/011357a0 .
- Моран, MJ; Шапиро, Х.Н. (2000). Основы инженерной термодинамики (4-е изд.). Нью-Йорк: МакГроу-Хилл.
- Ректорис, К. (1969). «Дифференциальное исчисление функций действительной переменной». В ректориях, К. (ред.). Обзор прикладной математики . Кембридж, Массачусетс: MIT Press. стр. 397–439.
- Шамсундар, Н.; Линхард, Дж. Х. (1983). «Насыщение и метастабильные свойства жидкости Ван дер Ваальса». Канада J Chem Eng . 61 (6): 876–880. doi : 10.1002/cjce.5450610617 .
- Зоммерфельд, А. (1956). Бопп, Ф.; Мейкснер, Дж. (ред.). Термодинамика и статистическая механика - Лекции по теоретической физике Том V. Перевод Кестина, Дж. Нью-Йорк: Academic Press.
- Тьен, CL ; Линхард, Дж. Х. (1979). Пересмотренная печать по статистической термодинамике . Нью-Йорк: Издательство Hemisphere.
- Ван Вайлен, Дж.Дж .; Зоннтаг, RE (1973). Основы классической термодинамики (второе изд.). Нью-Йорк: Джон Уайли и сыновья.
- Ван дер Ваальс, доктор юридических наук (1984). Роулинсон, Дж. С. (ред.). О непрерывности газообразного и жидкого состояний, под редакцией и с введением Дж. С. Роулинсона . Нью-Йорк: Dover Phoenix Editions.