First ‐ Principles Density Functional Theory Characterisation of the Adsorption Complexes of H 3 AsO 3 on Cobalt Ferrite (Fe 2 CoO 4 ) Surfaces

: The mobility of arsenic in aqueous systems can be controlled by its adsorption onto the surfaces of iron oxide minerals such as cobalt ferrite (Fe 2 CoO 4 ). In this work, the adsorption energies, geometries, and vibrational properties of the most common form of As(III), arsenous acid (H 3 AsO 3 ), onto the low ‐ index (001), (110), and (111) surfaces of Fe 2 CoO 4 have been investigated under dry and aqueous conditions using periodic density functional theory (DFT) calculations. The dry and hy ‐ droxylated surfaces of Fe 2 CoO 4 steadily followed an order of increasing surface energy, and thus decreasing stability, of (001) < (111) < (110). Consequently, the favourability of H 3 AsO 3 adsorption increased in the same order, favouring the least stable (110) surface. However, by analysis of the equilibrium crystal morphologies, this surface is unlikely to occur naturally. The surfaces were demonstrated to be further stabilised by the introduction of H 2 O/OH species, which coordinate the surface cations, providing a closer match to the bulk coordination of the surface species. The ad ‐ sorption complexes of H 3 AsO 3 on the hydroxylated Fe 2 CoO 4 surfaces with the inclusion of explicit solvation molecules are found to be generally more stable than on the dry surfaces, demonstrating the importance of hydrogen ‐ bonded interactions. Inner ‐ sphere complexes involving bonds between the surface cations and molecular O atoms were strongly favoured over outer ‐ sphere complexes. On the dry surfaces, deprotonated bidentate binuclear configurations ‐ cally favoured, whereas monodentate mononuclear configurations were typically more prevalent on the hydroxylated surfaces. Vibrational frequencies were analysed to ascertain the stabilities of the different adsorption complexes and to assign the As ‐ O and O ‐ H stretching modes of the ad ‐ sorbed arsenic species. Our results highlight the importance of cobalt as a potential adsorbent for arsenic contaminated water treatment. Compared with the dry surfaces, adsorption of H 3 AsO 3 the hydroxylated surfaces with the inclusion of solvation effects found to increase the stability of the H 3 AsO 3 adsorption complexes on the Fe 2 CoO 4 sur ‐ faces. Similar results were reported for H 3 AsO 3 and H 3 AsO 3 adsorption complexes on fer ‐ rihydrite (110) surface and for H 3 AsO 3 adsorption complexes on hydrated mackinawite (FeS) surfaces.


Introduction
Water is a fundamental resource for the sustainment of life, and so it is essential that access to clean water for drinking, sanitation, and other uses is feasible globally. Inorganic arsenic is the most significant contaminant in the world's drinking water. It is highly toxic, with a lethal range of inorganic arsenic estimated to be 1-3 mg As kg −1 , and is thus well known to be used for its toxic properties [1]. Arsenical compounds occur naturally in water due to volcanic activity, weathering of rocks, solubilisation, and transport of sediment. They are also becoming increasingly present in the environment as a result of human activities and the use of herbicides, pesticides, and waste products, accelerated by erosion of land and combustion of fossil fuels [2].
Adsorption is a widely used technique for the removal of arsenic and other species from water and other systems [3]. The adsorption of arsenic species onto iron oxide and iron oxyhydroxide mineral surfaces has been at the centre of many studies [4][5][6][7]. There remains, however, ambiguity in the possible reactions and adsorption mechanisms of such techniques. Molecular-level insight into the interactions of arsenic species with different iron oxide mineral surfaces is essential to understanding the adsorption mechanisms. Owing to the complex reactions occurring at the mineral surfaces, it is experimentally difficult to predict the nature of the adsorption complexes; whether monodentate or bidentate bonding is energetically favoured. However, using computer simulations based on density functional theory (DFT) calculations, it is now possible to unravel the adsorption complexes of pollutant species on mineral surfaces at the atomic level. For example, Dzade et al., employed DFT calculations to unravel the structure and properties of As(OH)3 adsorption complexes on mackinawite (FeS) [8] and ferrihydrite (110) [9] mineral surfaces. Goffinet and Mason employed spin-polarised DFT calculations to study innersphere As(III) complexes on hydrated α-Fe2O3(0001) surface models [10]. Corum et al., characterised the adsorption complexes of arsenate onto alumina surfaces using DFT calculations [11].
However, there is limited research dedicated to characterising the adsorption mechanism of arsenic on magnetic spinel B2CoO4 compounds (B = Cr, Mn, or Fe). Magnetic spinel materials have attracted a surge of interest recently for a wide range of applications including magnetic data storage, catalysis, and as contrast agents in magnetic resonance imaging, because of their high stability, resistance towards acid and alkali attacks, and magnetic properties [12]. Cobalt ferrite (Fe2CoO4) is a ferrimagnetic oxide like magnetite (Fe3O4). It exhibits remarkably unique and interesting properties such as high saturation magnetisation, high coercivity and anisotropy, good mechanical hardness, chemical stability, and has the highest magnetostriction of all 3d-element-based spinel oxides, with an additional benefit of being cheap and easy to make [13]. The preparation method and conditions of this spinel ferrite strongly influence the acquired degree of inversion and thus its properties [14][15][16][17][18]. In an investigation of the adsorption of both arsenite (H3AsO3) and arsenate (H3AsO4) on Fe2MnO4, Fe2CoO4 and Fe3O4, Zhang et al. synthesised the material via co-precipitation, an effective synthetic approach that permits good control of shape and particle size distribution [18,19]. It was demonstrated that the maximum adsorption capacities of arsenite and arsenate on Fe2CoO4 were 100 and 74 mg g −1 , respectively, which are higher than those for Fe3O4 (50 and 44 mg g −1 , respectively). This trend is in agreement with other experimental findings, indicating that cobalt ferrite has a greater potential for application in water purification [20].
Although there exist experimental studies about the magnetic properties and pollutants adsorption on Fe2CoO4 and its variants [13,[18][19][20], there are no atomistic first-principles studies dedicated to characterising the fundamental adsorption process and complexation of arsenite onto the Fe2CoO4 surface, which makes this study timely. In the present study, we have employed computational approaches based on DFT to characterise the surface structures, compositions, and stabilities of the inverse spinel Fe2CoO4 mineral. Subsequently, the structures and energetics of arsenous acid (H3AsO3) adsorption complexes have been predicted under anhydrous and hydrous conditions. Geometrical parameters and vibrational frequency assignment of the As-O and O-H stretching modes are reported and discussed. These results improve our understanding of the adsorption reactions and complexation of H3AsO3 on different cobalt ferrite surfaces at the molecular level.

Materials and Methods
All calculations were carried out using the VASP -Vienna Ab initio Simulation Package [21]. The projected augmented wave (PAW) method was used to describe the interactions between the valence electrons and the cores. The exchange-correlation potential was calculated using the Perdew-Burke-Ernzerhof (PBE) generalised gradient approximation (GGA) functional with a Hubbard correction (PBE+U), which accounts for the Coulomb interaction of localised d-electrons [22]. In this study, Ueff values of 4 eV for Fe and 5 eV for Co were found to provide sufficiently accurate lattice parameters, electronic, and magnetic properties of Fe2CoO4. Long-range van der Waals (vdW) interactions were accounted for using the method of Grimme (DFT-D3) [23]. A plane-wave basis set with a kinetic energy cut-off of 600 eV was tested to be sufficient to converge the total energy of Fe2CoO4 to within 10 −6 eV. Geometry optimisations were performed based on the conjugate-gradient algorithm until the residual Hellmann-Feynman forces on all relaxed atoms reached 10 −3 eV Å −1 . The Brillouin zone of the bulk Fe2CoO4 was sampled using a 5 × 5 × 5 Monkhorst-Pack K-point mesh [24].
The bulk Fe2CoO4 containing 56 atoms (Fe16Co8O32) was modelled in the inverse spinel cubic structure with space group Fd3 m. The crystallographic information file (CIF) of the normal spinel structure obtained from the Materials Project online database was used to create the inverse spinel structure, which was then subjected to geometry optimisation in order to determine its lattice constant and structural parameters [25]. From the fully relaxed bulk Fe2CoO4 material, the commonly observed (001), (110), and (111) surfaces were created using the METADISE code [26]. In each simulation, a vacuum of 15 Å was added in the c-axis to avoid interactions between periodic slabs. The surfaces were optimised using a k-point mesh of 3 × 3 × 1. From calculated surface energies, the equilibrium morphology of the dry and hydroxylated Fe2CoO4 nanoparticles was predicted based on Wulff's construction as implemented in the GDIS visualisation software [27]. Different coordination modes of H3AsO3 were studied including bidentate binuclear (BB) and monodentate mononuclear (MM) configurations in order to find the most stable coordination complexes. Large area surface slabs were used to characterise H3AsO3 adsorption: the (001), (110), and (111) surfaces have a surface area of 72.82, 102.98, and 126.12 Å 2 , respectively. The introduction of a single H3AsO3 molecule onto these large area surfaces ensured that the lateral interactions between the H3AsO3 species in neighbouring image cells, which may influence their adsorption geometries, are minimised. The strength of the adsorption for each adsorbate molecule (M) was determined by calculating the adsorption energy, (Eads) as follows: where Esurf+M, Esurf, and EM are the total energies of the whole system (dry or solvated adsorbate on the adsorbent slab), relaxed dry or solvated surface, and adsorbate molecule (M = H3AsO3), respectively. When optimising these systems, the adsorbate and the topmost four surface layers were allowed to relax unconstrainedly, while the remaining bottom layers were kept fixed to mimic the bulk structure. The stabilities of the different predicted arsenic adsorption complexes were ascertained by calculating their vibrational frequencies, for which stable structures are characterised by positive vibrational frequencies.

Bulk Properties
The inverse spinel Fe2CoO4 structure is characterised by half of the octahedral (Oh) sites being occupied by Co 2+ cations and the other half of the octahedral sites and all of the tetrahedral (Td) sites being occupied by Fe 3+ cations (Figure 1a). After full relaxation of the bulk Fe2CoO4 unit cell, the lattice constant is determined to be a = 8.435 Å, which is in good agreement with experimental (a = 8.533 Å) and previous GGA+U (8.41−8.46 Å) values [28,29]. The optimised Co-O, FeOh-O, and FeTd-O bond lengths are calculated at 2.10 Å, 2.04 Å, and 1.90 Å, respectively, in close agreement with previous GGA+U predictions [30]. In order to further investigate the material's electronic properties, the partial density of states (PDOS) of Fe2CoO4 was calculated as shown in Figure 1b. It is evident from the PDOS that the spin-up and spin-down density of states are asymmetric, confirming the characteristic ferrimagnetic nature of Fe2CoO4. The semiconducting behaviour of Fe2CoO4 is also well reproduced, predicting a bandgap of approximately 1.82 eV. The valence band is found to be dominated by the 2p states of O, whereas the conduction band is dominated by the 3d states of Fe. The features of the PDOS in the present study are quite comparable to those obtained by Das et al. [28] and Hou et al. [30] with GGA+U methods, although they observed smaller bandgaps estimated at 1.14 eV and 0.72 eV, respectively. The average magnetic moments of the tetrahedral and octahedral Fe 3+ ions in the bulk structure are predicted to be 4.016 and 4.154 μB, respectively (Table S1). The octahedral Co 2+ ions on other hand have a magnetic moment of 2.696 μB. Previous work by Das et al. [28] predicted the magnetic moments to be 3.98, 4.10, and 2.66 μB for Fe 3+ Td, Fe 3+ Oh, and Co 2+ , respectively, and Hou et al. [30] predicted similar values of 3.97, 4.10, and 2.61 μB, respectively.

Surface Characterisation
The low-index (001), (110), and (111) surfaces were created from the fully relaxed bulk Fe2CoO4 structure. Shown in Figure 2 are the optimised surface structures for the most stable terminations of the dry (naked) (001), (110), and (111) surfaces. The topmost layer of the dry (001) surface ( Figure 2a) is made up of two five-coordinate octahedral Co atoms (CoOh), two five-coordinate octahedral Fe (FeOh) atoms, and one three-coordinate tetrahedral Fe (FeTd) atom, which sits slightly higher than the FeOh atoms. The (110) surface ( Figure 2b), however, has two three-coordinate FeOh atoms on the topmost surface layer. The (111) surface ( Figure 2c) contains a three-coordinate Co and four three-coordinate Fe atoms, one of which comes from an FeOh in the bulk structure and is more raised than the others, which are FeTd. Another FeTd atom sits slightly lower than these, but still has all four of its bonds, and so is likely to be less reactive compared to the others nearby. All three surfaces contain a mixture of three-and four-coordinate oxygen atoms. Because of the decreased coordination number of the surface cations, we observed modification of their magnetic moments compared to the bulk to material (Table S1). In the case of the (001) surface, average values of 3.550 μB/FeTd, 4.092 μB/FeOh, and 2.701 μB/Co were obtained. The (110) surface had two equivalent terminating FeOh ions, both with a magnetic moment of 3.581 μB/FeOh, whereas FeTd ions in the second layer have magnetic moment of 3.851 μB. The (111) surface gave average values of 3.932 μB/FeTd, 3.526 μB/FeOh, and 2.599 μB/Co. The changes in these values compared with those in the bulk structure can be attributed to the decreased coordination numbers of the surface cations. This effect is less strong for the Co 2+ surface cations and remains similar to bulk values. As the coordination number of the FeOh decreased from 6 to 5 to 3, we see a coinciding change in magnetic moment from 4.154 μB in the bulk, to 4.092 μB for the five-coordinate FeOh on the (001) surface, to approximately 3.550 μB on the other two surfaces. The magnetic moment of the Co 2+ surface cations remains generally close those of the bulk magnetic moments. Considering that, in aqueous conditions, the surfaces may be hydrated (covered with water molecules) or hydroxylated (covered with OH species, resulting from water dissociation), we have first quantified the energetics of molecular and dissociative water adsorption on the three different surfaces. This is to ascertain whether the surfaces will preferentially be covered with molecular adsorbed water or OH species resulting from water dissociation. Generally, dissociative water adsorption is found to be energetically favoured over molecular water adsorption on all three surfaces (Supplementary Materials, Figure S1, Table S2), suggesting that the Fe2CoO4 surfaces will be hydroxylated under normal conditions. Having found that the water molecules will preferentially adsorb dissociatively on the Fe2CoO4 surfaces, we created a monolayer of OH-covered (001), (110), and (111) surfaces by introducing 7, 4, and 6 water molecules adsorbed dissociatively (Figure 3a-c). During energy minimisation, some of the OH/H species at the (001) surface recombined to form H2O on the (001), two of which were physisorbed. Seeing that the hydration/hydroxylation may have a stabilising effect on the Fe2CoO4 surfaces, the relative stabilities of the Fe2CoO4 surfaces were determined by calculating the surface energies of the dry surfaces (γr) and hydroxylated surfaces (γhydrox.) using Equations (2) and (3).

2
( 2) where Esurf is the energy of the slab, Ebulk is the energy of the bulk Fe2CoO4 material, n is the number of unit cells used in the slab, and A is the surface area of one side of the slab.
wherein Ehydrox. is the energy of the relaxed hydroxylated surface, x is the number of H2O molecules used to hydroxylate the surface, EH2O is the energy of the isolated H2O molecule, and γr is the relaxed surface energy of the corresponding surface, as determined by Equation (2). The relaxed surface energy of the dry surfaces is calculated at 1.454, 2.121, and 1.882 Jm −2 for the (001), (110), and (111) Fe2CoO4 surfaces, respectively. These results show that the order of increasing surface energy, and thus decreasing stability for the dry surfaces, follow the trend (001) < (111) < (110). When hydroxylated, the surface energy is predicted at −0.011, 0.993, and 0.596 Jm −2 for the (001), (110), and (111) Fe2CoO4 surfaces, respectively. In comparison with the dry surfaces, all surface energies were reduced considerably upon hydroxylation. Most notably, the (001) surface energy became negative, signifying that it is now significantly more stable than the other two hydroxylated surfaces.
The order of stability remained as (001) < (111) < (110), as was predicted for the dry surfaces. Based on the calculated surface energies, the equilibrium morphology of the dry and hydroxylated Fe2CoO4 nanoparticles was predicted in order to determine which facets are present and/or dominant. For the dry nanoparticle (Figure 4a), it is observed that the crystal adopts a truncated cubic shape, with the most stable (001) surface having the largest surface area, followed by the second most stable surface, (111), truncating the corners of the cube. The absence of the (110) surface in the dry nanoparticle may be attributed to its relatively high surface energy [31]. The hydroxylated nanoparticle (Figure 4b), on the other hand, contains only the (001) lowest energy surface, forming a cubic crystal structure, consistent with its significant stabilisation upon hydroxylation relative to the other two surfaces. Although the surface energies of the (110) and (111) surfaces were also lowered considerably, these reductions were not significant enough to compete with the (001) surface for expression in the nanoparticle.

Adsorption Complexes of H3AsO3 on Dry Fe2CoO4 Surfaces
Although contamination of water by arsenic species is a more globally prevalent concern, we can benefit from investigating the adsorption of H3AsO3 in the absence of water, as arsenic-containing species such as H3AsO3 can be found in soils and sediments and other non-aqueous environments. Hence, the adsorption complexes and properties of H3AsO3 were investigated on the dry Fe2CoO4 surfaces. In order to find the most stable adsorption geometries, a variety of monodentate mononuclear (MM) and bidentate binuclear (BB) modes were explored. The resulting unique low-energy adsorption configurations are shown in Figure 5, with the calculated adsorption energies summarised in Table 1   The most stable intact H3AsO3 geometry on the (111) surface was bound by an O-Fe bond, releasing an adsorption energy of −1.746 eV. The As-O bond length of the Omol atom bound to the surface increased to 1.999 Å. The equivalent structure with a monodentate O-Co bond instead had an adsorption energy of only −0.943 eV, indicating that this surface has more favourable interactions with the Fe surface atoms. This can be explained by the fact that the Fe ion involved here is tetrahedral, and so adsorption at this site would complete its coordination sphere, whereas, when adsorbing at Co sites, these are three-coordinate octahedral sites, so the coordination spheres remain incomplete upon adsorption with H3AsO3. This favourability can also be seen in comparing the adsorption energies of the two deprotonated bridging AsO geometries. The most thermodynamically stable deprotonated complex on the (111) surface is a BB-Fe-OO-Fedp geometry, in which one of the O atoms bound to Fe was deprotonated, releasing an adsorption energy of −3.262 eV. The other deprotonated BB-Co-AsO-Fedp and BB-Fe-AsO-Fedp released adsorption energies of −2.130 and −2.250 eV, respectively.
On the least stable (110) surface, the preferred H3AsO3 geometry turned out to be a bidentate mononuclear (BM) complex, BM-OO-Fe, in which two of the molecule's O atoms were bound to a single FeOh surface atom with O-Fe bond lengths of 2.158 and 2.434 Å. The adsorption energy of the BM-OO-Fe complex is −1.832 eV, the most stable among all the protonated complexes on all three surfaces. These results show that the strength of adsorption is related to the stability of the surface, where the H3AsO3 binds most strongly to the least stable (110) surface. As expected, this surface also produced the most stable deprotonated adsorption complex (BB-O-AsO-Fedp), which released an adsorption energy of −3.942 eV. The deprotonated As-O bond length is predicted at 1.653 Å, which is shorter than the protonated As-O bond (1.821 and 1.772 Å). The most stable coordination modes are found to be the deprotonated complexes in which there is a proton transfer to form surface hydroxyl species. These results suggest that, on the dry Fe2CoO4 surfaces, H3AsO3 species will preferentially exist in deprotonated states.

Adsorption Complexes of H3AsO3 on Hydroxylated Fe2CoO4 Surfaces
The unique optimised geometries of H3AsO3 on the hydroxylated (001), (110), and (111) Fe2CoO4 surfaces are shown in Figure 6 and their structural properties are summarised in Table 2. Solvation effects were accounted for through the inclusion of four explicit water molecules near the H3AsO3 adsorbate on the (001), (110), and (111) surfaces. The adsorption of H3AsO3 onto the hydroxylated surfaces is found to be significantly stronger than on the dry surfaces. This is expected as the surface hydroxyl species and the solvating water molecules form hydrogen-bonding interactions with the H3AsO3 species, which contribute to the stability of both the inner-and outer-sphere adsorption complexes on the hydroxylated surfaces. The most thermodynamically favourable of all adsorption geometries was the monodentate mononuclear MM-O-Co complex on the (001) surface, which was deprotonated and released an adsorption energy of −4.457 eV, which is almost twice the energy released by the MM-O-Co complex on the dry surface. The deprotonated As-O bond shortened to 1.737 Å and the remaining two As-O bonds lengthened compared with those in the non-aqueous system, owing to the additional hydrogen-bonded interactions introduced by the H2O molecules. The strength of the hydrogen-bonded interactions in H3AsO3 adsorption complexes on the hydroxylated Fe2CoO4 surfaces is quantified by measuring all possible hydrogen-bonded interactions with nearest O_H2O (solvated or adsorbed) and OHsurf, as shown in Table S3

Vibrational Properties
In order to gain further insight into the stabilities of each uniquely identified adsorption complex on the dry and hydroxylated Fe2CoO4 surfaces, the vibrational frequencies of the adsorbed H3AsO3 species were calculated and the As-O and O-H stretching modes were assigned. These were also compared to the vibrational frequencies of the gas-phase H3AsO3 molecule. All predicted adsorption structures are deemed stable considering that almost all normal modes were positive, except for a few cases in which low imaginary modes were observed, corresponding to translational modes, and not to the As-O or O-H stretching modes. A clear and predictable trend is observed on the dry and hydroxylated surfaces, in which the bonds of the shortest length, and thus greatest strength, give rise to the highest vibrational frequencies. This is amplified for the deprotonated As-O bonds. The bond lengths for deprotonated As-O bonds are shortest, thus the atoms are more tightly bound and stable, resulting in a more rigid vibration of higher frequency. Tables 3  and 4   The three most stable configurations on the (001) surface, as determined by their adsorption energies, are those that do not give rise to any imaginary vibrational frequencies in this analysis, confirming their thermodynamic stability. The most stable protonated adsorbing geometry is the MM-O-Co configuration, and it is the only protonated geometry on this surface with no imaginary frequencies. Its weakest As-O bond with a frequency of 522 cm −1 is given by the O molecule bound directly to the Co surface atom. It is seen that increased coordination onto the surface causes a weakening of the bonds directly involved in the coordination, and hence reduces the vibrational frequencies. This is seen for exam-  (111) surface contained the most imaginary frequencies of all configurations. Its adsorption energy of −0.943 eV is significantly lower than the other configurations at this surface, and so, although it is still sufficiently stable, it is least likely to occur in a real system. Furthermore, this configuration is the one that most resembles the frequencies of those in the isolated H3AsO3 molecule.
No imaginary frequencies were found for any of the solvated H3AsO3 adsorption complexes on the hydroxylated surfaces, suggesting that they are all stable, more so than those on the dry surfaces. As seen in Table 4, the As-O bond stretching frequencies are generally greater for the outer sphere, physisorbed, H3AsO3 complexes than the innersphere, chemisorbed complexes, except when those Omol atoms are deprotonated, again highlighting the weakening of these bonds upon adsorption onto the surfaces. As shown in Table S5, we observed elongation (>1 Å) of the O-H bonds (reduction in stretching frequencies), most especially for those for those forming hydrogen bonds with the surrounding O atoms of the surface hydroxyl groups and water molecules. The elongated O-H bonds and the deprotonated O-H ones are found to have stretching frequencies below 2600 cm −1 , whereas the shorter O-H bonds (<1 Å) are typically characterised by higher frequencies >3000 cm −1 . The change in the O-H bonds follows the same reasoning as above, in which longer bonds give rise to lower frequencies.

Conclusions
Using periodic DFT+U calculations, we have performed an atomistic-level investigation of the adsorption of H3AsO3 onto the low-index Fe2CoO4 (001), (110), and (111) sur-faces under dry and aqueous conditions, including analysis of the coordination geometries, electronic, and vibrational properties. The strength of H3AsO3 adsorption is demonstrated to be related to the stability of different surfaces; it increases in the order (001) < (111) < (110). Along with the Wulff plots of the crystal morphologies, the surface energies of both sets of surfaces suggest that the reactive (110) surface is unlikely to exist unless engineered to be present in the structure. Hydroxylated surfaces were highly favoured over hydrated ones, indicating that these systems would operate optimally under intermediate-to-high pH conditions. The adsorbing H2O/OH species also bring the coordination of the surface atoms closer to that in the bulk material, stabilising the surfaces. H3AsO3 also tended towards deprotonation, which can be attributed to the fact that less energy is required to break a weak O-H bond in H3AsO3 than what is released upon forming an Osurf-H bond. Monodentate mononuclear geometries were more prevalent than bidentate geometries and were largely the more thermodynamically favoured configurations in the solvated systems. In the absence of water, deprotonated bidentate binuclear geometries were generally more favourable. Solvation of the systems is found to enhance adsorption due to the introduction of hydrogen-bonding interactions between the H3AsO3 molecule and surrounding water species on the surface and in solution. The vibrational frequencies associated with each coordination mode and the crystal morphologies further support these findings and provide additional useful insight for future research into this system for arsenic adsorption. The adsorption characteristics of H3AsO3 on Fe2CoO4 surfaces in the present study are similar to, but are overall stronger than, those observed on the ironcontaining materials ferrihydrite and mackinawite, highlighting its capacity to be a potential adsorbent over other alternatives [9,33]. It is thus expected that these results should be comparable to those of other iron oxides and ferrites that may be investigated in the future. The work presented in this report provides a good foundation for further computational and experimental research into this H3AsO3-Fe2CoO4 system and provides a molecular-level understanding of the adsorption mechanisms. Future investigations may expand the work presented here to include the effects of coverage (larger surface areas and multiple H3AsO3 species), surface defects (oxygen vacancies), and the use of classical molecular dynamics (MD) simulations, which will provide a complete description of the dynamic processes occurring at the H3AsO3-water-Fe2CoO4 interfaces.
Supplementary Materials: The following are available online at www.mdpi.com/2075-163X/11/2/195/s1, Figure S1: Low-energy adsorption geometries for molecular and dissociative H2O on the Fe2CoO4 (001), (110), and (111) surfaces (Fe = brown, Co = blue, Osurf = red, Omol = purple, H = white), Table S1: Magnetic moments of FeTd, FeOh, and Co atoms in the bulk and surfaces of Fe2CoO4, Table S2: Adsorption energies and Omol-metal bond lengths of molecular and dissociated H2O geometries on the Fe2CoO4 surfaces,   Data Availability Statement: Information on the data that underpins the results presented here, including how to access them, can be found in the Cardiff University data catalogue at http://doi.org/10.17035/d.2021.0128750010. use of the facilities of ARCHER (http://www.archer.ac.uk), the U.K.'s national supercomputing service via our membership of the U.K.'s HEC Materials Chemistry Consortium, which is funded by EPSRC (EP/L000202).

Conflicts of Interest:
The authors declare no conflict of interest.