Next Article in Journal
Stress and Strain Prediction of Zirconium Nitride under Oxygen Doping and Vacancy Introduction
Previous Article in Journal
Horocycles of Light in a Ferrocell
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Electron Lone-Pairs Stereochemistry and Drastic van der Waals and Pressure Effects in AsF3 from First Principles

1
CNRS-LCTS, Université de Bordeaux, 33600 Pessac, France
2
Lebanese German University, LGU, Jounieh 206, Lebanon
*
Author to whom correspondence should be addressed.
Formerly, senior researcher at CNRS-ICMCB, Université de Bordeaux.
Condens. Matter 2021, 6(3), 31; https://doi.org/10.3390/condmat6030031
Submission received: 26 June 2021 / Revised: 5 August 2021 / Accepted: 6 August 2021 / Published: 12 August 2021
(This article belongs to the Section Condensed Matter Theory)

Abstract

:
AsF3E as a representative of a molecular crystal has been chosen to find the precise localization of the lone pair (LP) E centroid 4s2 of As3+ and to enlighten the behavior of lone pair triplets of fluorine atoms. Starting from stereochemistry rationale, Density Functional (DFT) electronic structure calculations yielding the electron localization (ELF) mapping led to precise large crystal structure evolutions from basic X-rays data (V = 267.2Å3 at 193K), to (V = 230.5Å3) and under Van der Waals forces (DEW) V = 206.4Å3, and then under pressure P, all illustrated with ELF maps and band structures. Calibrated pressures up to 100 GPa exhibit the remarkable shrinking of all inter-atomic distances including As-E from 0.94Å down to 0.46Å, while the major three bonds As-F1, As-F2 and As-F3 are continuously expanding. The resulting picture of the application of pressure on AsF3 molecular structure leads to the progressive immersion of the 4s2 doublet within the electronic cloud with an original new status of As with its 4s2 E expressed as [AsE]3+. The expanded size becomes close to big LP-free K+, Ba2+ cations or LP characterized bismuth cation: [BiE]3+, which appear inserted in cubes or square anti-prisms with coordination number 8 (CN8) or in tricapped trigonal prisms with CN9.

1. Introduction

Following our earlier research works on electron lone pair (LP) ns2 elements [1,2], some particular features on fluoride compounds of group VA trivalent elements prompted us to carefully settle their stereochemistry and to include the fluoride anions LP triplets. For this purpose, we focused on AsF3E, where E designates AsIII 4s2 lone pair, to precisely address the features pertaining firstly to the stereoactivity of AsIII lone pair E on one hand, and on the other hand to focus on the significant role of the LPs brought by the fluorine atoms of the three pending As-F bonds. Additionally, AsF3 is a molecular solid state in its crystalline form [3], therefore it was a good candidate for highlighting and precising the stereochemistry and metrics of the involved two natures of LPs and their implication in structural evolution such as with pressure. It needs to be mentioned here that a recent study dealing with SnII where the eminent role of the LP triplets in overlapping [SnF4]n layers was quantitatively shown [2].
Stressing further the original features above, the 4s2 LP of AsIII is not alone in influencing the atomic network. Indeed, in the present study, we show that the (2s22p4) configuration of fluorine LP’s triplets which are all non-bonding play a paramount role in the architecture of the molecule, and consequently on the shape of the three dimensional packing of the crystal.
For a complete picture of the stereochemistry of the title compounds, our investigations merge crystal chemistry rationale and electronic structure calculations. Specifically, the results allow extracting accurate information pertaining on one hand to the electron localization (cf. next Section 1), the stereochemistry and sterical effects of LPs and to provide quantitative LP metrics, and on the other hand to follow AsF3 evolution at the atomic scale by accounting for van der Waals effects and under pressure (GPa) influence [4]. Let us note that during condensation of cell volume due to the effects of simulated pressure, there is a formation of original [As2F6E2]n potential polymer chains along [100] and [001] around 10–15 GPa occurring through AsF4E triangular bipyramids units sharing axial summits, and resulting in a three-dimensional crystal compound at 70–100 GPa closely related to β BiF3E [5], characterized by a big bismuth cation.
Beside the LP metrics, further assessments of the electronic structure focusing on simulated van der Waals (vdW) and pressure (P) effects will be provided from the electronic band structures, mirrored by the site projected density of states and further illustrated through the qualitative assessment of the chemical bonding illustrating the interactions between pairs of atoms.
The paper is organized as follows: After this contextual Introduction, the theory and computational frameworks are presented in Section 2. The stereochemistry of AsF3 is described in details in Section 3, followed by the results of LP metrics together with the vdW and pressure (P) effects detailed in Section 4. Section 5 details the electronic band structures, the site projected DOS and properties of chemical bonding. The last section highlights the conclusions.

2. Theoretical and Computational Frameworks

The ab initio electronic structure calculations herein are based on the quantum density functional theory (DFT) [6,7], where the effects of exchange and correlation (XC) are equally taken into account, opposite to the Hartree–Fock method, which mainly treats exchange effects. In the founding papers, the first development of XC effects was operated thanks to local density approximation LDA [8] scheme based on the homogeneous electron gas. However, the need for introducing gradients to the electron density to correctly describe solid state chemical compounds rapidly led to the generalized gradient approximation GGA [9] used herein.
In the present work, we used two methods within DFT. The Vienna ab initio simulation package (VASP) code [10,11] with the projector augmented wave PAW method [11,12] for generating atomic potentials of As and F was used to obtain equilibrium crystal structures starting from experiment; those are usually found close to experimental findings. This method also allows accounting for van de Waals effects following Grimme scheme [4] within so-called DFT-D2 approach implemented in the MEDEA simulation package with VASP [13]. The main outcome of the calculations is to ‘build’ the electron localization around the chemical constituents, i.e., the atoms with their chemical trend to ionize either positively (decrease of localization) or negatively (increase of localization). Two different schemes are practiced for this purpose: the electron localization indicator (ELI-D) [14] and the electron localization function (ELF) [15]; the latter is used herein. ELF scheme is based on the kinetic energy in which the Pauli Exclusion Principle is included. ELF is a normalized function meaning, that its magnitude ranges from 0 to 1: 0 for no electron localization (blue zones) and 1 (red zones) for strong localization in the plane-slices the halfway value of ½ corresponds to a free electron gas–like behavior (green zones), cf. Figures. These color codes identified in slices of planes crossing atomic zones of interest, are also accompanied by 3D iso-surfaces, which are relevant when analyzing electron lone pairs E stereo activity for the quantitative discussion of the LP metrics.
Detailed calculations of the electronic structure and the chemical bonding properties were carried out using the full-potential augmented spherical wave (ASW) method in its scalar-relativistic implementation [16], wherein the effects of exchange-correlation are treated with GGA [9]. In the ASW method, the wave function is expanded in atom-centered augmented spherical waves, which are Hankel functions and numerical solutions of Schrödinger’s equation, respectively, outside and inside the so-called augmentation spheres. In order to optimize the basis set, additional augmented spherical waves were placed at carefully selected interstitial sites. In the minimal ASW basis set, the outermost shells were explicitly used to represent the valence states with partial waves up to lmax +1 for the different atomic species where lmax is the maximum secondary quantum number; here, lmax= 2 for As and lmax = 1 for F. Filled As 3d10 states were not considered in the valence basis set. Self-consistency was achieved until the variation of the atomic charges was smaller than 10−8 electrons and the variation of the total energy was smaller than 10−7 eV. Further, the ASW method implements the chemical bonding for pair interactions provided as based on the overlap integral Sij (i and j designating two chemical species) within the crystal orbital overlap population (COOP) by Hoffmann [17]. Here, we use the integrated COOP (iCOOP, unit less), whereby the larger the area underneath the curve, the larger the bonding is.

3. Stereochemical Considerations

3.1. Arsenic Trifluoride AsF3E

In view of the Introduction section where fluorine LP (E) role is accounted for besides that of trivalent arsenic, we now express AsF3 as As{EAs}F3{EF}9.
An electron diffraction study of AsF3 molecule was early reported [18], establishing an As-F bond length of 1.70Å. In 1970, an electron diffraction study refining AsF3 molecule parameters with As-F = 1.71 Å and F-As-F angle = 96.2° data was published [3]. Finally, in 1979, an X-ray structure determination on single crystal was reported [19] and discussed with Va element trihalides [20].

3.2. Crystal Formation and Structure

The crystal structure of AsF3 was reported by Enjalbert and Galy [19] followed by crystal structure determination of other MF3 compounds, where M is the Va element [20]. AsF3 shows a typical molecular network with an AsF3E triangular pyramidal model, arsenic sitting in a triangular pyramid with a base of three F and As 4s2 lone pair E at the apex as shown in Figure 1 see upper left.
Cell parameters determined by X-ray were used for basic features of the structure, the data involving E lone pair for As and the lone pair triplets {EF} for fluorine atoms were issued from DFT–ELF analyses (DE) + Van der Waals (VdW); they complete the panel of AsF3E data before applying increasing pressures from 3 up to 100 GPa (present work: pw). DE + VdW, designated by DEW, serve to create the CIF base (crystal information file) for structural investigations under chosen pressures (GPa). As{EAs}F3{EF}3 crystallizes in the orthorhombic system with space group Pna21 as reported in Table 1.

3.3. Description and Crystal Chemistry Rationale

The molecular network shows largely separated AsF3E molecules exhibiting with cationic AsIII and its lone pair E the classical AsF3E triangular pyramidal (TP) coordination shape (Figure 1a). The molecules coil around 21 screw axes (along [001]), and they are packed in planes parallel to (100) and (010). However, the packing reveals a remarkable organization of five molecules around the molecule model AsF1F2F3Ec and more precisely around its lone pair Ec giving six fluorine atoms to complete a tricapped triangular prism around [AsE]. The two bases, F1F2F3 and F1a2a3a, of the prism, completed by F2c, F3b and F1b capping the three rectangular faces, are shown in Figure 1b. The Ec-F1,F2,F3 average distance amounts to 2.38Å, while the six other fluorine atoms enclosing [AsE] show Ec-F distances ranging between 2.58Å and 3.15Å (average 2.74Å). In Figure 1b, these fluorine atoms surrounding [AsF3E] molecule make a closed polyhedron, a distorted tricapped triangular prism (TTP), which merges [AsE].

4. Results for Stereochemical Assessments

4.1. Electron Localization Function ELF

From the calculation results, selective ELF sections allow us to enter in a precise location and description of lone pairs. Coordination details of AsF3 are summarized in Figure 2.
In Figure 2a,b, the sections by AsF1F3 and AsF1F2 reveal an important electronic density (this term being used indifferently from localization herein) on the top of As with the slice of As lone pair 4s2 of volume E. By using special software (cf. ref. [2]), the coordinates of its maximum density were determined and refined. This clear Ec centroid location allowed us to calculate the section illustrated in Figure 2c determined by AsF3Ec plane.
Subsequently, the coordinates of the center of gravity G of F1, F2 and F3 was calculated allowing us to precisely find a better symmetric section of Ec in GEcF1 plane. A precise volume of Ec and of its associated electronic cloud completed by a section perpendicular to As-Ec in Ec (Figure 2e) allowed us to determine the parameters of an Ec ellipsoid corresponding to a sphere of influence (Ec) in the packing with a radius rEc = 0.94Å. This is detailed in forthcoming figures. Therefore, it is worth noting the particular shape of the figure when we focus on fluorine atoms: They show two maxima around F1, F2 and F3 immersed in a density rounded rectangle (Figure 2f).
To set up these tori searches of density, maxima were systematically performed in the vicinity of the fluorine atoms. Figure 3a illustrates this search by several points with marked density shown with blue-green circles. After the precise localization of their coordinates in the network, they show that they were at the same distance from F with F-{EF} = 0.56Å and characterized with a center of gravity CT (CT-{EF} = 0.53Å) (Table 2).
Figure 3a shows, for example, an ELF section based on these points {EF} around As-F axes. It has been possible to link these points by a circle (green) with a radius CT-{EF} representing the torus gyration circle of {EF}, LP’s triplet locus of F. Figure 3b is of paramount importance, looking at the highest ring designed by the two IDC = 0.835 density curves (the size of the ring is indicated by green arrow). A constant ring density is shown, establishing that it results from the {E3} triplet rotating around the As-CT axis. The triplets {EF1} and {EF2} show identical sizes.
The gyration diameter is drawn in the middle of a ring of stable IDC designated by a yellow circle, with a diameter estimated at r = 1.06Å. The three {EFa,b,c} LPs are represented by pale blue transparent circles, their center (little clear blue circles) being settled on the yellow gyration circle. Their radii determined in Figure 2 and Figure 3, r{EF} = 0.45Å, suggest they are separated by 120°. Figure 3a shows the coherent dynamic crowns formed by these tori capping F atoms with their centers aligned on yellow gyration circle perpendicular in CT to As-F direction (angles: As-F1-CT1 = 179.5°, As-F2-CT2 = 179.9° and As-F3-CT3 = 179.9°) and distant by F-CT = 0.20Å.
The particular “dynamic” surface generated by the fluorine triplets lets us suppose that they participate in the high reactivity of this fluoride, especially through its extreme hygroscopicity.

4.2. Calculations Accounting for van der Waals Forces and Simulation of Various Pressures

From the above developments, the AsF3E structure, perfectly organized, lets us strongly suggest that subtle molecular interactions, at short and long ranges, were likely to cause a phase transition accompanied by network reorganizations after being exposed to pressures with increasing magnitudes. It was decided to apply pressures up to 100 GPa.
Parameters relevant to the cell volume V and cell dimensions (a,b,c) are reported in Table 1 and depicted in Figure 4a–d. A first approach was accomplished with the simulation of structure network already treated by DFT-ELF (DE) as previously done by adding up van der Waals forces (vdW). Such a DFT + D2 interaction is treated via a pair-wise force field optimized for DFT functional according to Grimme’s development in the MEDEA package [4].
Coming to these high pressures a first test at isotropic pressure of only 3 GPa V (DEW + 3 GPa) shows again an extra cell condensation, its volume decreasing from 206.4Å3 to 193.4Å3. Two examples at 20 GPa and 70 GPa are given (Figure 5c,d), cell volume decreasing down to 133 2Å3. Some large interaction in [100] and [001] directions, sensitive to high pressure, are marked by dotted red sticks, and their lengths are indicated (Figure 5).
The cell shape indicates its exponential decreasing evolution up to 100 GPa. Concomitantly, at the level of AsF3E molecule, the pressure shows atom by atom that the AsF3E atomic network evolves to a progressive phase transition, via its condensation, largely engaged around 70 GPa (Figure 4a–d). Some molecular data evolution, As-E, and As-F1,2,3, are directly illustrated by diagrams versus pressure steps (Figure 4c,d). A remarkable impact is revealed through the fundamental bonds As-F1,2,3 which, astonishingly, are all increasing up to P = 100 GPa. Therefore, another fact of paramount importance concerns the lone pair E, which starts to coil the As atom following a large diminution of As-E distance from 0.97 to 0.46Å. This highlights that the lone electronic centroid pair Ec and its associate electronic cloud E are very sensitive to the pressure, which pushes As3+ to the creation of a new cation: As3+ with its lone pair, i.e., [AsE]3+.
In a perspective cell view, some large As-F interactions are designed in Figure 5 by dotted red sticks roughly oriented to the three directions [100], [010] and [001] their length being included within the range 2.89Å–3.18Å.
Evolution of the structure from AsF3E from (a) DFT-ELF to; (b) 3 GPa application indicating the coherent large As--F interactions contraction; (c) an example at 20 GPa; (d) 70 GPa situation which indicate that the structure evolution as raised a state where a polymer circuit is traced along [100] (as an example) marked by large blue sticks from …F1b-Asb-F1a-Asa-F1b-Asb…corresponding eventually to a potential phase {As2F6E2}n, therefore, also the beginning of phase [AsE]F3 formation isostructural to orthorhombic β [Bi.E]F3 phase with Bi-E = 0.28Å (SG Pnam, a = 6.56, b =4.84, c =7.15Å, V = 222.9Å3 [5]) or a cubic one the synthetic γ [Bi.E]F3 with Bi-E = 0.13Å (SG P 4 ¯ 3m, a= 5.853Å, V = 200.5Å3) [21].
Then, the cell parameters evolution up to 100 GPa associated to the lengthening of primary As-F and As-E distances show that the last pressure was decisive to settle a new phase [AsE]F3 with a tri-dimensional network and a volume of 123.8Å (Table 1). A perspective view of the structure is given in Figure 6a. The detail of the crystal unit [AsE]F3 in its CN = 9 number as a form of tri-capped triangular prism (TTP) is displayed in Figure 6b. Such virtual [AsE]F3 phase under 100 GPa is directly comparable to β [BiE]F3. Both formulae correspond astonishingly to an isostructural network. Following the As-E evolution (Figure 4d) and its formula of the curve; a value for its derivative is found at P = 102 GPa, and d(As-E) = 0.467Å is identified as very close to the data corresponding to P = 100 GPa.
Such evolution lets us propose that a supplement of pressure could transform this [AsE]F3 structure into isostructural cubic form of BiF3 P 4 ¯ 3m (215) [21] and perhaps α gananite form [22] that are related to typical YF3 crystal structure as shown by Hyde and Andersson [23].

5. Electronic Band Structures, Density of States and Chemical Bonding

Further assessment of the van der Waals (vdW) effects, on one hand, and of the influence of pressure, on the other hand, regarding the electronic structure in AsF3 were subsequently carried out using the full potential augmented spherical wave method (ASW) [17]. ASW method also implements the crystal orbital overlap population (COOP) criterion following Hoffmann [18] to address the chemical bonding as based on the overlap population. The calculations were done using the crystal parameters in Table 1.
Figure 7 and Figure 8 present the results at Normal Pressure (N.P.) and P = 3 GPa, respectively, for the band structure calculated along the main lines of the simple orthorhombic Brillouin zone (BZ), followed by the site projected density of states (DOS) and the bonding rationalized through the integrated COOP (iCOOP). In all panels, the zero of energy is with respect to the top of the (filled) valence band VB, separated from the (empty) conduction band (CB) by an energy gap at different magnitudes (vide infra) indication insulating character of the trifluoride.

5.1. Normal Pressure Results

The energy band structures in Figure 7 exhibit differences between “off” and “on” vdW effects. Without vdW effects, the band gap Egap is large, close to 6 eV. Upon including vdW effects, it decreases down to Egap(vdW) = 5 eV. Additionally, the localized character of flat bands within both VB and CB is characteristic of a highly ionic chemical compound, whereby there is a clear separation of s-like bands centered at −10 eV from p-like bands extending over 5 eV, up to EV; and finally a narrow CB of ~2 eV breadth is observed. Upon tuning on the vdW effects, drastic changes occur with a significant broadening of all of three band blocks, i.e., by 2 eV or s-like, 8 eV for p-like and 5 eV for the CB. The overall picture is that of a less-isolated chemical system with intermolecular electronic interactions.
The observations are mirrored by the electronic density of states where the low energy lying DOS correspond to a broadening of the sharp As s-like states mixing specifically with F1 within a broad ~2 eV band. The p-block is broadened with a differentiation between the F states, especially F1, and the CB is mainly composed of As empty states.
In so far that the bonding is ensured by the directional p-states, the iCOOP top panel shows significant intensity essentially within the p-block energy range. The As-F bonds are almost identical for the three kinds of As-Fi (i = 1,2,3) interactions, but large changes occur in the lower panel upon accounting for vdW effects where As-F1 iCOOP intensity prevails versus As-F2 and As-F3 within the broader p-like VB, in perfect agreement with Figure 5b) showing infinite chains: ….F1-As-F1….

5.2. Pressure Effects

The major effects discussed above are reproduced in Figure 8, illustrating the consequence of exerting pressure. Expectedly, increasing pressure reduces the interatomic distances with the consequence of a larger overlap between the valence states, i.e., with broader bands. This is actually observed in the band structure panels exhibiting smaller band gaps.
The DOS panels show the same trends as in Figure 7 regarding the relative energy extension magnitudes of the s- and p-blocks, as well as the gap magnitudes.
Lastly, the iCOOP shows the remarkable As-F1 differentiated larger magnitude characterizing As-F1 bonding.

6. Conclusions

In this work a thorough analysis highlighting the major changes brought into the molecular structure of arsenic trifluoride AsF3E submitted via DFT-ELF to van der Waals corrections and high pressure treatments was performed.
The first point concerns the valence ns2 E lone pair of As, i.e., 4s2, which assumes with the three F atoms and its centroid Ec and electronic cloud E, the basal plane and the apex of the AsF3E triangular pyramidal molecule. Around [AsE] we have noted the presence of six fluorine atoms which form, with the three fluorine atoms involved with the As-F bonds, a closed tricapped triangular prism [AsE]F9 polyhedron.
The lone pair triplets of each fluorine atoms linked to As form a torus crowning each of them, and it is believed that such a dynamic of fluorine atoms of pending bonds must have a major impact in their high reactivity.
Coming to Van der Waals contraction of the cell, large length As--F interactions are all clearly reduced; this phenomenon being enhanced with the pressure increase.
However, a most astonishing fact is the increasing lengths of the three bonds As-F1,2,3, which amounts to ~117% between room and 100 GPa pressures. Such a phenomenon is well-balanced by the contraction of As-Ec, which enclose them, after continuous bond diminution from 0.96Å to 0.46Å, into a tricapped triangular prism, noted [AsE]F9. Accordingly, the volume [AsE]F93+ can be seen as a large cation like K+, Ba2+ or more precisely to [BiE]3+ in the series of α, β, and γ [BiE]F3 fluorides. However, here we learn that with pressure increase, the distance As-E diminishes but [AsE] size increases, its combination being illustrated in another occurrence with bismuth, namely in orthorhombic γ [BiE]F3 with Bi-E = 0.28Å, γ [BiE]F3 with Bi-E = 0.13Å and β [BiE]F3 with Bi-E = 0 the Gananite mineral with {[BiE]F8 cube surrounded by six {F}F8 cubes.
This stereochemical study indicates that it is possible with programmed high pressure to master the lone pair and to press it to merge into the cation, here AsIII, thus inducing it to behave, structurally, like a heavy atom as Bismuth.

Author Contributions

Both authors J.G. and S.F.M. have contributed equally to the conceptualization, methodology, formal analysis and writing—original draft preparation. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data supporting reported results including CIF files can be made available on demand.

Acknowledgments

One of us (J.G.) thanks Gérard Vignoles, LCTS-CNRS, for providing facilities and enriching discussions. Computational facilities from the University of Bordeaux and from the Lebanese German University are gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Galy, J.; Matar, S.F. ns2np4 (n = 4, 5) lone pair triplets whirling in M*F2E3 (M* = Kr, Xe): Stereochemistry and ab initio analyses. Solid State Sci. 2017, 64, 114–124. [Google Scholar] [CrossRef] [Green Version]
  2. Galy, J.; Matar, S.F. Joint stereochemical and ab initio overview of SnII electron lone pairs (E) and F(E) triplets effects on the crystal networks, the bonding and the electronic structures in a family of tin fluorides. J. Prog. Solid State Chem. 2019, 56, 100252. [Google Scholar] [CrossRef]
  3. Clippard, F.B., Jr.; Bartell, L.S. Molecular structures of arsenic trifluoride and arsenic pentafluoride as determined by elec-tron diffraction. Inorg. Chem. 1970, 9, 805–811. [Google Scholar] [CrossRef]
  4. Grimme, S. Semi-empirical GGA-type density functional constructed with a long-range dispersion correction. J. Comput. Chem. 2006, 27, 1787–1799. [Google Scholar] [CrossRef] [PubMed]
  5. Greis, O.; Martinez Ripoll, M. Crystal structure of BiF3. Z. Anorg. Allg. Chem. 1977, 436, 105–112. [Google Scholar] [CrossRef]
  6. Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864–B871. [Google Scholar] [CrossRef] [Green Version]
  7. Kohn, W.; Sham, L.J. Self-Consistent Equations Including Exchange and Correlation Effects. Phys. Rev. 1965, 140, A1133–A1138. [Google Scholar] [CrossRef] [Green Version]
  8. Ceperley, D.M.; Alder, B.J. Ground State of the Electron Gas by a Stochastic Method. Phys. Rev. Lett. 1980, 45, 560–566. [Google Scholar] [CrossRef] [Green Version]
  9. Perdew, J.; Burke, K.; Ernzerhof, M. The Generalized Gradient Approximation Made Simple. Phys. Rev. Lett. 1996, 77, 3865–3868. [Google Scholar] [CrossRef] [Green Version]
  10. Kresse, G.; Furthmüller, J. Efficiency of ab-initio total-energy calculations for metals and semiconductors using a plane-wave basis set. Phys. Rev. B Condens. Matter Phys. 1996, 54, 11169–11186. [Google Scholar] [CrossRef]
  11. Kresse, G.; Joubert, J. From ultrasoft pseudopotentials to the projector augmented wave. Phys. Rev. B Condens. Matter Phys. 1999, 59, 1758–1775. [Google Scholar] [CrossRef]
  12. Blöchl, P.E. The projected augmented wave method. Phys. Rev. B Condens. Matter Mater. Phys. 1994, 50, 17953–17979. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. MEDEA simulation Package. Materials Design Co. Available online: https://www.materialsdesign.com/software-releases (accessed on 7 August 2021).
  14. Bezugly, B.V.; Kohout, P.; Wagner, F.R. Electron localizability indicators ELI-D and ELIA for highly correlated wavefunctions of homonuclear dimers. II. N2, O2, F2, and Ne2. J. Comput. Chem. 2010, 31, 1404–1519. [Google Scholar] [CrossRef] [PubMed]
  15. Becke, A.D.; Edgecombe, K.E. A simple measure of electron localization in atomic and molecular systems. J. Chem. Phys. 1990, 92, 5397. [Google Scholar] [CrossRef]
  16. Eyert, V. Basic notions and applications of the augmented spherical wave method. Int. J. Quant. Chem. 2000, 77, 1007–1012. [Google Scholar] [CrossRef]
  17. Hoffmann, R. How chemistry and physics meet in the solid state. Angew. Chem. Int. Ed. Engl. 1987, 26, 846–848. [Google Scholar] [CrossRef]
  18. Pauling, L.; Brockway, L.O. The radial distribution method of interpretation of electron diffraction photographs of gas molecules. J. Am. Chem. Soc. 1935, 57, 2684. [Google Scholar] [CrossRef]
  19. Enjalbert, R.; Galy, J. Structure cristalline du trifluorure d’arsenic AsF3 à 193K. C.R. Acad. Sc. 1979, C289, 441–443. [Google Scholar]
  20. Galy, J.; Enjalbert, R. Crystal chemistry of the Va element trihalides: Lone pair, stereochemistry, and structural relationships. J. Solid State Chem. 1982, 44, 1–23. [Google Scholar] [CrossRef]
  21. Hund Von, F.; Fricke, R. Der Kristallbau von α-BiF3. Z. Anorg. Allg. Chem. 1949, 258, 198–204. [Google Scholar] [CrossRef]
  22. Croatto, U.N. The gananite mineral. Gazs. Chim. Ital. 1944, 71, 20. [Google Scholar]
  23. Hyde, B.G.; Andersson, S. Inorganic Crystal Structures; Wiley-Interscience Publication: Hoboken, NJ, USA, 1989; pp. 122–123. [Google Scholar]
Figure 1. Perspective views of AsF3E: (a) AsF3E molecular network; (b) the distorted fluorine tricapped triangular prism (TTP) of nine fluorine atoms surrounding [AsE], i.e., [AsE]F9.
Figure 1. Perspective views of AsF3E: (a) AsF3E molecular network; (b) the distorted fluorine tricapped triangular prism (TTP) of nine fluorine atoms surrounding [AsE], i.e., [AsE]F9.
Condensedmatter 06 00031 g001
Figure 2. Arsenic trifluoride ELF. (a,b) Sections based on As, F1, F2, and F3 atoms showing the presence of As lone pair Ec. The (c,d) sections allow us to find some precise metrics of the lone pair Ec of the couple [AsE]; (e) this Ec section is directly obtained by calculating a plane perpendicular to As-Ec in Ec, allowing us to refine a and b parameters of Ec ellipsoid. (d) Ec centroid is marked by a little red-white circle and pale blue ellipse circling also the center (blue-white circle) of electronic cloud ellipsoid. (e) Plane of Ec perpendicular to As-Ec. (f) Section based on the three fluorine atoms F1, F2, F3 (G their gravity centre) slicing also densities due to fluorine lone pairs 2s22p4.
Figure 2. Arsenic trifluoride ELF. (a,b) Sections based on As, F1, F2, and F3 atoms showing the presence of As lone pair Ec. The (c,d) sections allow us to find some precise metrics of the lone pair Ec of the couple [AsE]; (e) this Ec section is directly obtained by calculating a plane perpendicular to As-Ec in Ec, allowing us to refine a and b parameters of Ec ellipsoid. (d) Ec centroid is marked by a little red-white circle and pale blue ellipse circling also the center (blue-white circle) of electronic cloud ellipsoid. (e) Plane of Ec perpendicular to As-Ec. (f) Section based on the three fluorine atoms F1, F2, F3 (G their gravity centre) slicing also densities due to fluorine lone pairs 2s22p4.
Condensedmatter 06 00031 g002
Figure 3. (a) Projection onto (010) of the complex AsEF3{EF3}3 molecules. The three LPs of each F atom, {EF3}, have been localized (blue little circles). These points are located in a circle, the center (CT red circle) of which is moved away by 0.2Å. Section plane in CT perpendicular to As-F direction shows the generation by the three LPs {EFa,b,c} of the electronic torus (green circle) around each fluorine. (b) The ring with maximum intensity (IDC = 0.835) and thickness (green arrow) is the locus of gyration circle indicated by yellow ring, the diameter of which diameter amounts 1.06Å.
Figure 3. (a) Projection onto (010) of the complex AsEF3{EF3}3 molecules. The three LPs of each F atom, {EF3}, have been localized (blue little circles). These points are located in a circle, the center (CT red circle) of which is moved away by 0.2Å. Section plane in CT perpendicular to As-F direction shows the generation by the three LPs {EFa,b,c} of the electronic torus (green circle) around each fluorine. (b) The ring with maximum intensity (IDC = 0.835) and thickness (green arrow) is the locus of gyration circle indicated by yellow ring, the diameter of which diameter amounts 1.06Å.
Condensedmatter 06 00031 g003
Figure 4. (a) Volume (V) of the AsF3E cell versus pressure p (GPa) decreasing up to 100 GPa; (b) cell parameters a,b,c vs. p; (c) As-F bonds of initial molecule AsF3E increasing with pressure, and (d) the drastic negative evolution of As-E spacing.
Figure 4. (a) Volume (V) of the AsF3E cell versus pressure p (GPa) decreasing up to 100 GPa; (b) cell parameters a,b,c vs. p; (c) As-F bonds of initial molecule AsF3E increasing with pressure, and (d) the drastic negative evolution of As-E spacing.
Condensedmatter 06 00031 g004
Figure 5. AsF3E. (a) DFT-ELF (DE) situation where large As--F interactions are marked by dotted red sticks, (b) vdW and 3 GPa applied drastic contraction of As--F interactions was noted; (c) situation onto 20 GPa and (d) under 70 GPa pressure where contraction of these As--F interactions perpetuate a polymeric chain can be outlined (highlighted by a thick pale blue stick (4d)). All intermolecular interactions are drastically reduced (some indicated by dotted red sticks and their values in red numbers) but, concomitantly, As-F1,2,3 original bonds increase.
Figure 5. AsF3E. (a) DFT-ELF (DE) situation where large As--F interactions are marked by dotted red sticks, (b) vdW and 3 GPa applied drastic contraction of As--F interactions was noted; (c) situation onto 20 GPa and (d) under 70 GPa pressure where contraction of these As--F interactions perpetuate a polymeric chain can be outlined (highlighted by a thick pale blue stick (4d)). All intermolecular interactions are drastically reduced (some indicated by dotted red sticks and their values in red numbers) but, concomitantly, As-F1,2,3 original bonds increase.
Condensedmatter 06 00031 g005
Figure 6. (a) Detail of the AsF3E structure under 100 GPa pressure; all 9 fluorine atoms setting As and E are the green-red circles. (b) The formed fluorine tri-capped triangular prism (TTP) surrounding [AsE]3+, which behaves like large K+, Ba2+ or more precisely like [BiE]3+ cations in γ and α [BiE]F3 forms.
Figure 6. (a) Detail of the AsF3E structure under 100 GPa pressure; all 9 fluorine atoms setting As and E are the green-red circles. (b) The formed fluorine tri-capped triangular prism (TTP) surrounding [AsE]3+, which behaves like large K+, Ba2+ or more precisely like [BiE]3+ cations in γ and α [BiE]F3 forms.
Condensedmatter 06 00031 g006
Figure 7. AsF3 electronic structure at normal pressure: band dispersion (left), DOS (middle) and bonding (right), without (top) and with (bottom) van der Waals effects.
Figure 7. AsF3 electronic structure at normal pressure: band dispersion (left), DOS (middle) and bonding (right), without (top) and with (bottom) van der Waals effects.
Condensedmatter 06 00031 g007aCondensedmatter 06 00031 g007b
Figure 8. AsF3 electronic structure at P = 3 GPa pressure: band dispersion (left), DOS (middle) and bonding (right), without (top) and with (bottom) van der Waals effects.
Figure 8. AsF3 electronic structure at P = 3 GPa pressure: band dispersion (left), DOS (middle) and bonding (right), without (top) and with (bottom) van der Waals effects.
Condensedmatter 06 00031 g008
Table 1. X-rays [19,20], DFT-ELF (DE), DE + VdW (DEW) and some DEW + GPa data of AsF3E{EF3}3. Space group Pna21. Data: Bonds, interatomic distances (Å) and angles (°) from RX [19] and DEW + GPa [data pw]. (Ec centroid of As 4s2 lone pair).
Table 1. X-rays [19,20], DFT-ELF (DE), DE + VdW (DEW) and some DEW + GPa data of AsF3E{EF3}3. Space group Pna21. Data: Bonds, interatomic distances (Å) and angles (°) from RX [19] and DEW + GPa [data pw]. (Ec centroid of As 4s2 lone pair).
AsF3E/Cell (RX [21])a (Å)b (Å)c (Å)V (Å3)ZrV (F,E) (Å3)
RX [21] AsF3E7.0185.2057.315267.2416.7
Atom coordinates of the molecule and E localization (Figure 1)AtomsxyzWyckoff
As0.78620.53380.254a
F10.57360.39950.30654a
F20.80770.37400.04904a
F30.67920.79560.14934a
E0.88850.55290.33714a
AsF3E/Cellsa (Å)b (Å)c (Å)V (Å3)ZrV (F,E) (Å3)
AsF3E (DE)6.6324.9577.013230.5414.4
AsF3E (DEW)6.3464.7806.806206.4412.9
DEW+ 3 GPa6.1744.6816.691193.4412.1
DEW+ 30 GPa5.7014.3536.339157.349.8
DEW+ 100 GPa5.3193.9105.954123.847.7
Data
(Å, °, GPa)
RX [21]
(Å, °)
DEW + GPa
(Å, °, 3 GPa)
DEW + GPa
(Å, °, 30 GPa)
DEW + GPa
(Å, °, 100 GPa)
As-Ec0.970.940.690.46
As-F1 1.6991.7981.8341.966
As-F2 1.6961.8051.8341.966
As-F31.7221.8151.8612.036
F1-F2 2.5032.5932.5352.380
F2-F3 2.4832.5012.3952.174
F3-F1 2.4742.5432.4912.381
Asb-F1a3.0542.5462.2531.918
As-F2c2.8902.4222.2182.032
F-As-F angles (°)95.0, 93.2, 92.789.5, 87.4, 92.187.4, 80.8, 84.874.5, 65.8, 73.0
Table 2. AsEcF3{EF3}. Data involving E or [EF] lone pairs as issued from DFT–ELF analyses.
Table 2. AsEcF3{EF3}. Data involving E or [EF] lone pairs as issued from DFT–ELF analyses.
As-Ec0.97Ec-F1,2,32.38∠FAsF *93.6 *
F-{EF} (1,2,3)0.57F -CT (1,2,3)0.20Ec a/b/c ellipsoid 0.94/1.03/0.85
CT-{EF}(1,2,3)0.53As-CT (1,2,3)1.91{EF}a/b/c ellipsoid0.41/0.41/0.55
T. Rd. rect. {EF}
(Figure 2a,b) L/H/R
1.98/1.10/0.75ϕT.ext.{EF3}
blue dotted circle
1.98T. gyr. ϕ yellow stick or ϕ circle 1.06
rEc0.94 r{EF}0.45
Ec coordinates: Wyckoff 4a site x,y,z: 0.112, 0.447, 0.837. [pw]: present work. Ec: Oblate ellipsoid a/b/c: 0.94/1.03/0.85Å. * average. {EF} density maxima found around each F1,2,3 atoms. T. rd.rect.[EF] L/H/R: Rounded rectangle embedding {EF} (red line) lone pair triplet which generates a fluorine electronic torus by LP {EF}3 rotation around As-F bonds capping F atoms. CT1,2,3 tori centers of {EF1,2,3} triplets. T gyr.ϕ {EF}3: torus gyration diameter = 1.06Å.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Galy, J.; Matar, S.F. Electron Lone-Pairs Stereochemistry and Drastic van der Waals and Pressure Effects in AsF3 from First Principles. Condens. Matter 2021, 6, 31. https://doi.org/10.3390/condmat6030031

AMA Style

Galy J, Matar SF. Electron Lone-Pairs Stereochemistry and Drastic van der Waals and Pressure Effects in AsF3 from First Principles. Condensed Matter. 2021; 6(3):31. https://doi.org/10.3390/condmat6030031

Chicago/Turabian Style

Galy, Jean, and Samir F. Matar. 2021. "Electron Lone-Pairs Stereochemistry and Drastic van der Waals and Pressure Effects in AsF3 from First Principles" Condensed Matter 6, no. 3: 31. https://doi.org/10.3390/condmat6030031

APA Style

Galy, J., & Matar, S. F. (2021). Electron Lone-Pairs Stereochemistry and Drastic van der Waals and Pressure Effects in AsF3 from First Principles. Condensed Matter, 6(3), 31. https://doi.org/10.3390/condmat6030031

Article Metrics

Back to TopTop