1.2模擬參數設置

1）模型參數與構建

采用Materials Studio中的Visualizer模塊構建模擬所需的分子及離子，並利用Forcite模塊分配立場與電荷，圖1a為各種分子及離子優化後的模型。其次，利用Crystals模塊創建含有優化好的1個MSDS陰離子與1個AEO-5分子的9×18×30Å3（a×b×c）晶胞結構，通過Symmetry中的Supercell功能，將構建好的晶胞結構沿a、b和c方向分別擴展6、3和1個，得到1個54×54×30Å3的表麵活性劑分子層，如圖1b所示。然後，利用Amorphous Cell模塊構建1個由1 500個水分子和無機鹽離子組成，尺寸為54×54×30Å3的水盒子，如圖1c所示。最後，依據甘建建立的泡沫液膜模型，利用Build Layer功能，得到54×54×200Å3的“三明治”泡沫模型，如圖1d所示。基於截斷半徑要求，模擬體係Z軸方向水層上下各添加5 nm厚度的真空層。

圖1模型參數與構建

2）動力學參數設置

分子動力學模擬過程由Forcite模塊完成。初始泡沫體係的結構並不穩定，需對泡沫初始模型進行結構優化以獲得穩定的低能量泡沫結構。在分子動力學模擬開始前，利用Forcite模塊中的Anneal對幾何優化後的泡沫模型進行退火操作。利用Dynamic進行分子動力學模擬，係綜、力場、溫度等相關參數設置情況見表1。

表1模擬參數設置

2.結果與討論

2.1界麵張力

通常在分子動力學模擬中，研究人員通過計算體係壓力張量的分量來得到氣液體係的界麵張力，具體可參見式（1）：

式中：γ為泡沫液膜的表麵張力；Px，Py和Pz分別為壓力張量在x，y和z方向的分量；Lz為模型中z軸方向的長度。

通過腳本文件統計各泡沫體係平衡後200幀的軌跡文件，得到不同濃度含鹽泡沫體係界麵張力的變化情況，如圖2所示。由圖2可知，不同濃度含鹽泡沫體係的界麵張力隨濃度的增加而增大，低、中濃度時泡沫體係界麵張力從大到小依次為Ca2+>Mg2+>Na+，而高濃度時泡沫體係界麵張力從大到小依次為Mg2+>Ca2+>Na+。界麵張力是體現泡沫穩定性的相關參數之一。界麵張力越小，吸附在氣液界麵處的表麵活性劑活性越高，泡沫越穩定。各泡沫體係界麵張力的變化表明，硬水離子的存在會對泡沫穩定性產生不利影響，且Ca2+在低、中濃度時產生的負麵效果較Mg2+更明顯，而高濃度時，Mg2+較Ca2+產生的負麵效果更明顯。

圖2不同濃度含鹽泡沫體係的界麵張力

2.2界麵結構

界麵結構可通過相對數密度分布表征，計算體係中水分子沿界麵垂直方向的分布可得出泡沫界麵結構。由於泡沫模型關於水層中心對稱，故在計算水分子沿Z軸相對數密度分布時，以水層中心位置為原點，取水層上、下層數據的均值進行作圖，得到不同濃度含鹽泡沫體係中水分子相對數密度沿Z軸的分布情況，如圖3所示。

圖3不同濃度含鹽泡沫體係中水分子的相對密度曲線

由圖3，水分子相對數密度從液相到氣相逐漸減少到0，在液相區域內，離子濃度越高水分子的相對數密度越低。其中Ca2+與Mg2+泡沫體係的變化規律大體一致，而Na+泡沫體係則有所差異。對於Ca2+、Mg2+泡沫體係而言，低濃度時水分子的相對數密度在除氣相區域外的任何區域內均大於中、高濃度時的相對數密度，中、高濃度曲線在氣液界麵區域內交於1點後，中濃度時水分子的相對數密度開始小於高濃度時水分子的相對數密度。但中、高濃度曲線的交匯處明顯不同，Ca2+泡沫體係在9.2Å附近產生交匯，Mg2+泡沫體係在7.2Å附近產生交匯，兩者交匯時水分子的相對數密度均為8.1。Na+泡沫體係中，高濃度時水分子的相對數密度在除氣相區域外的任何區域內均小於中、低濃度時的相對數密度，中、高濃度曲線在氣液界麵區域內交於1點後，中濃度時水分子的相對數密度開始大於高濃度中水分子的相對數密度。

溶劑密度從10%升到90%的範圍為界麵層，通過相對數密度分布曲線可估算出界麵厚度。不同濃度含鹽泡沫體係界麵厚度變化的火柴圖，如圖4所示。由圖4可知，不同濃度含鹽泡沫體係界麵厚度隨濃度的變化規律一致。即當離子種類一定時，高濃度含鹽泡沫體係的界麵厚度最小，低濃度次之，中濃度最大。對於低濃度含鹽泡沫體係而言，Mg2+的界麵厚度最大，Ca2+次之，Na+最小；中濃度含鹽泡沫體係中，Ca2+的界麵厚度最大，Mg2+次之，Na+最小；高濃度含鹽泡沫體係中，Mg2+的界麵厚度最大，Ca2+次之，Na+最小。綜上可知，硬水離子的存在使得泡沫界麵厚度有所增加，這意味著界麵中水分子內部的相互作用減弱，使得界麵水分子受到的指向泡沫內部的不平衡拉力減弱。

圖4不同濃度含鹽泡沫體係界麵厚度的火柴圖

3.結論

礦井水質硬度越高，泡沫體係的界麵張力越大，高硬水條件下鈣、鎂泡沫體係的界麵張力平均值大於100 mN/N，該條件下的界麵厚度最薄，抗幹擾能力最低，泡沫穩定性最差。