Rational design of mono- and di-anions as super acids through π-hole

Rational design of mono- and di-anions as super acids through π-hole interaction: Implications for Lithium and ... Publication Date (Web): July 2, 20...
0 downloads 0 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

C: Energy Conversion and Storage; Energy and Charge Transport

Rational design of mono- and di-anions as super acids through #hole interaction: Implications for Lithium and Magnesium Ion Batteries Mrinal Kanti Si, and Bishwajit Ganguly J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.8b01286 • Publication Date (Web): 02 Jul 2018 Downloaded from http://pubs.acs.org on July 3, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Rational Design of Mono- and Di-anions as Super Acids Through π-hole Interaction: Implications for Lithium and Magnesium Ion Batteries Mrinal Kanti Si1,2, Bishwajit Ganguly*1,2 1

Computation and Simulation Unit (Analytical Discipline and Centralized Instrument Facility), CSIR-Central Salt & Marine Chemicals Research Institute 2

Academy of Scientific and Innovative Research, CSIR-CSMCRI, Bhavnagar, Gujarat, India364 002. *Corresponding Author. Fax: (+91)-278-2567562, E-mail:[email protected]; [email protected] Abstract: Nature helps to generate multiply charged negative ions in abundance. Such negatively charged ions are stable as crystals with charge compensating cations, but in solution, they are protected with solvent molecules. In the gas phase, such multi-charged anions are not stable due to strong electrostatic repulsion between the extra electrons. Therefore, the generation of multi-charged anions without the influence of environment has been of great interest. The DFT calculations reveal that the π-hole interactions can stabilize the anion of acid using selenoaldehydic group attached to the cyclopentadiene ring to make them super acids. This strategy helped to generate the stability of dianions without any influence of environmental conditions. It is known that most electrolytes used in Li-ion batteries contain halogens, which are however toxic in nature. The designed super-acids can also be exploited as lithium ion batteries. These stable dianions showed the binding energies of lithium ions (5.8 and 63.6 kcal/mol) at M06-2X/6-311+G(2d,p) level of theory and can be considered as an efficient electrolyte in lithium ion battery. The study extended with magnesium ion also showed promise to function as efficient electrolyte in magnesium ion battery using such anions.

1 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction: The proton is a tiny particle that plays a pivotal role in many chemical and biological phenomena. The principal importance of proton is in acid-base chemistry. According to Arrhenius theory, an acid produces H+ ion in solution. In the dissociation of an acid HA, one important consequence in equilibria (HA

H+ + A-) is that every acid (HA) has a conjugate

base (A-) and vice-versa. The strong acids are molecular compounds that disassociates into H+ ions and anions in aqueous solution. The equilibrium shifts to the right hand side in the (HA

H+ + A-) means the HA can be the strong acid and the A- is the stable conjugate

base. Recently, efforts have been laid in the design of potential neutral organic superacids.1-5 Polycyanated hydrocarbon derivatives have been designed and examined, which behaves as hyperacids.6,7 The parent hydrocarbons of such polycyanated hydrocarbons exhibit tautomerism, which results a large number of tautomers. Some of these tautomers are less stable and possess lower acidity. The substitution of cyano groups in these hydrocarbons generate more number of tautomers where the hydrogen moves around the molecule leads to stable ketene imine species. Such phenomena yields hyper strong superacids in some cases.1 The tautomerism was also seen to play an important role in organic tri- and hexacyclic molecules to generate super- to hyper acids.8 A few strategies in preparing superacids have been reported. In the first strategy, electronic super acceptor substituents were used, however, in second strategy stabilization of the conjugate base through strong anionic resonance triggered by strong electron withdrawing groups.1,6,7,9 Yagupolskii has given the concept of electronic superacceptor substituents that can be showed by =NSO2CF3, where replacing doubly bonded oxygen to S, P and I atoms gives very strong superacids in the gas phase and in acetonitrile. This strategy is based on the use of one or more highly dipolar superacceptor groups and some strongly polarizable

2 ACS Paragon Plus Environment

Page 2 of 28

Page 3 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

substituents yields an extended conjugated system in their corresponding conjugate bases.10-12 In the second strategy, the anionic resonance can stabilize the conjugate base to increase the gas phase acidity of molecules. The increase in the acidity can relate to the increase in ionization potential of an anionic conjugate base, which can be calculated in Koopmans’ frozen electron density + clamped nuclei approximation.13 The large IPKoop ionization potentials indicate strong anionic resonance or aromatic stabilization upon deprotonation.10,13 These strategies have generated a large number of extremely powerful superacids in the gas phase and in solution.1,6,7,9 Reports on the stability of dianions are known in the literature, though such anions are relatively unstable due to electronic repulsions.14-17 Recently, the dianion of ortho-diethynylbenzene has been synthesized and the proton affinity was calculated using DFT methods.14 A stable and flexible dianion 2-dicyanomethylene1,1,3,4,5,5-hexacyanopentenediide

(DHCP2-)

has

also

been

prepared

as

a

tetraalkylammonium salt from hexacyanobutadiene.15,18-21 Such strong acids can be used for organic synthesis and catalyst for coal liquefaction, hydrocarbon isomerization, cracking, alkylations used in olefin polymerization and also isolation of many short-lived cations.6,22-29 Carboranes and boron clusters have been exploited to make an interesting class of strong acids.6,25,30 Researchers have also shown that tautomeric fulvene derivatives can also act as hyperacids.31 The stable anions have been shown to be important as materials for batteries.16,17,32 The reports suggest that the stability of multicharged anions can be exploited as electrolyte for metal ion batteries.16,17,32 It is reported in the literature that the halogen electrolytes are used in Li-ion battery, however, halogens are toxic in nature.32 Therefore, there are efforts to design the halogen free electrolyte for Li-ion battery.32 The DFT studies reveal that LiCB11H12 has the potential among all halogen-free electrolytes. The metal borohydrides such as LiBH4 can be used as potential candidates for halogen free electrolytes in Li-ion battery.32

3 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 28

Scheme 1 Structures of substituted selenoaldehydic cyclopentadiene derivatives (1-6) .

H

R1 H

Se

H

H

H H

R2

Se R4

R3 H

R5

H

Se

H R9 H

H

R6

R8

H R7 H

1

2-6

2. R1= H, R2= H, R3= H, R4= H, R5= H, R6= H, R7= H, R8= H, R9= H 3. R1= CN, R2= H, R3= H, R4= H, R5= H, R6= H, R7= H, R8= H, R9= CN 4. R1= CN, R2= CN, R3= H, R4= H, R5= H, R6= H, R7= H, R8= CN, R9= CN 5. R1= CN, R2= CN, R3= CN, R4= H, R5= H, R6= H, R7= CN, R8= CN, R9= CN 6. R1= CN, R2= CN, R3= CN, R4= OCH3, R5= H, R6= H, R7= CN, R8= CN, R9= CN

In this article, we have given a different concept to design superacids using non-covalent interactions using DFT calculations. Further, we have shown that the conjugated system with electron donating group instead of electron acceptor with π-hole bonds can arrive in the range of superacids. Reports reveal that the anion can be stabilized using non-covalent interactions such as anion…π, π…π, σ-hole and π-hole interaction.28 The π-hole bond is another depiction of the anion/lone pair…π interaction.33 The description of π-hole bond is extensive and includes organic and inorganic π-hole systems involving π-bond or conjugated π-bond.33-36 Generally, anion…π interaction conceives of a negative…negative interaction, nevertheless, the intrinsic characteristic of the “π” is described well in the π-hole. Such π-holes describe a region of positive surface electrostatic potential in a molecular entity, and the π-hole bond is also distinct. Generally, the π-molecular orbital or p orbital of one central atom is involved in π-hole bond. It is known that such non covalent interactions i.e., σ- and π-hole interactions are fully directional.33 We have incorporated for the first time the role of π-hole interactions to augment the acidity of cyclopentadiene (Cp) derivatives. Such strategy may lessen the synthetic work to introduce 4 ACS Paragon Plus Environment

Page 5 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

a number of substituents in the organic molecules.1,7 We have also shown that the di-anions can also be stabilized through non-covalent interactions. The stability of these dianions can be used as electrolyte in lithium and magnesium ion batteries. Such small molecules with such properties would be promising for various applications in chemistry and materials including batteries as storage devices. There are reports where non-covalent interactions have been exploited to influence the proton affinity of YHXOH (X = 0-3; Y = Li, Be, B, Na, Mg, Al etc.) in presence of Lewis bases. The Lewis bases form the complexes with YHXOH (X = 0-3; Y = Li, Be, B, Na, Mg, Al etc.) through non-covalent interaction which in turn enhances the proton affinity.37 Computational Methods: We have calculated the acidity of the compounds using the enthalpy change (∆Hacid) for the proton dissociation reaction. The proton dissociation reaction is given below AH(g) → A- (g) + H+(g)

(1)

The ∆Hacid are calculated using the following equation 2 ∆Hacid = ∆Eacid + ∆(pV)

(2)

Where ∆Eacid denote the change of energies of the species involved in the equation 1 where the energy includes total electronic energy and thermal enthalpy correction at 298.15 K. The ∆(pV) is pressure-volume work. Gibbs free energy change (∆Gacid) has been calculated using the species involved in eq. (1) at the same level of theory. We have considered the enthalpy of proton in gas phase as 1(5/2)RT or 6.197 kJ/mol at 298.15 K and Gibbs free energy of proton of “(5/2) RT - TS,” which is -26.25 kJ/mol at 298.15 K.31 We have optimized the structures at M06-2X/6-311+G(2d,p) level of theory in gas phase to calculate the electronic energies and Hcorrection.38-40 The optimized structures are minima which are confirmed by the positive harmonic vibrational frequency. The DFT M06-2X is 5 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 6 of 28

recognized to be well described with dispersion forces.41 Further, we have performed the single point of these molecules in DMSO solvent using the IPCM (Isodensity Polarized continuum model) solvent model.42,43 The acidity of the molecules can be calculated using the equation 3 AH + DMSO → A- + DMSOH+ + ∆rHDMSO

(3)

∆rH is the enthalpy change in the above equation. We have performed the calculations at M06-2X/6-311+G(2d,p)//M06-2X/6-311+G(2d,p) level of theory in solvent using isodensity polarized continuum model (IPCM) by using the thermal correction at M06-2X/6311+G(2d,p) level of theory. The pKa values for compound cyclopentadiene and 6 are calculated by using following equation which is previous reported to calculate the pKa for neutral organic acid.44 pKa(exp) = 0.661. ∆rHDMSO – 7.7

(4)

We have performed the AIM (The Bader’s theory of atom in molecule) analysis taking the wave function generated at M06-2X/6-311+G(2d,p)//M06-2X/6-311+G(2d,p) level of theory in GAUSSIAN09.45-47 The AIM has been performed using multiwfn software.48 The MESP (Molecular electrostatic potential) has been performed using Molekel-4.3 free software. The following equation has been used to calculate the electrostatic potential in space around molecules.49 V(r) = Σ ZA /|RA-r| - ʃρ(ŕ)d ŕ/|ŕ -r|

(5)

The V(r) is the potential at the distance r by the nuclei and electrons of the molecule; and ZA is the charge of nucleus A at the point RA. The ρ(r) is electron density. We have optimized the structure of Li-82-, Li2-82- and Mg-82- (Figure 5) at M06-2X/6311+G(2d,p) level of theory to get a quantitative estimate of the cation binding energy. The

6 ACS Paragon Plus Environment

Page 7 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

energies necessary to remove the first and second lithium ion as well as the magnesium ion have been calculated from the following equations.17 ∆ELi1 = E(Li) + E(Li18-) – E(Li28)

(6)

∆ELi2 = E(Li+) + E(82-) – E(Li18-)

(7)

∆EMg = E(Mg2+) + E(82-) – E(Mg8)

(8)

Where ∆ELi1 and ∆ELi2 are the dissociation energies for removal of first and second Li+ ions from the Li2-82- respectively (Figure 5). ∆EMg is the dissociation energy of magnesium of Mg-82-. The binding energies of 1st and 2nd electron of halogen-free electrolytes (8) are calculated in the following equation 9 and 10.12 ∆E1 = E(8) – E(8-)

(9)

∆E2 = E(8-) – E(82-)

(10)

The structure of 8 is given in Figure 4. Results and Discussion The stability of cyclopentadienyl anion with substituted selenoaldehyde was examined using density functional theory (DFT) calculations (Scheme 1). The (-CH2-CH2-) chain length is appropriate to achieve the interaction between the selenoaldehyde and the cyclopentadienyl anion. The smaller or longer chain lengths, however, do not align properly for any interaction (Figure 1). The selenoaldehyde attached with cyclopentadiene unit showed the remarkable stability of the corresponding conjugate base. The acidity of cyclopentadiene measured experimentally is 352.6 kcal/mol.31,44 The calculated ∆Hacid for cyclopentadiene with M062X/6-311+G(2d,p) level of theory is 351.0 kcal/mol which is in very good agreement with the experimental results. Compound 1 with selenoaldehyde showed the deprotonation energy computed at the same level of theory is 344.0 kcal/mol. The selenoaldehyde group can

7 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

possess σ- and π-holes for effective interaction with the electron rich systems.33,50 The σ-hole is present along the extension of the C=Se bond and π-holes are present above and below of the F2C=Se plane formed by the carbon and fluorine atoms.33,50 The σ- and the π-holes enable a molecule to interact attractively with a negative site such as the lone pairs, π electrons or anions.33 The strength of σ- and π-hole bonds can be modulated with the substituents attached to it. The electron withdrawing nature of the substituents enhances the ability of σ- and π-hole to interact strongly with the electron rich sites of the system.33,50 The presence of selenoaldehyde group in 1 remarkably eased the deprotonation by ~8.0 kcal/mol compared to parent cyclopentadiene. The substituted selenoaldehyde group aligns to interact with the cyclopentadienyl anion and resides at a distance of 3.64 Å (Figure 1). The Molecular electrostatic potential (MESP) calculations reveal that the π-hole of selenoaldehyde interacts with the negative site of cyclopentadienyl anion (Figure 2). The selenoaldehyde aligns properly to interact with the conjugate base, while in the case of parent cyclopentadiene, it is away from the ring (Figure 2). We have performed the MESP analysis to examine the stability arises by the π-hole interactions in conjugate bases of cyclopentadiene derivatives (Figure 2).49 The negatively charged cyclopentadienyl anion in 1 (indicated by red colour) is partially neutralized by the π-hole generated on selenoaldehyde. The negative region of cyclopentadienyl unit is dissipated in the MESP of 1- (Figure 2) causes the depletion of positive potential around -C=Se unit. This analysis clearly depicts the stabilization of charge within the system. To augment the stability of carbanion, additional selenoaldehyde was substituted in 2. The calculated ∆Hacid for 2 is 340.4 kcal/mol (Table 1). This result suggests that the cyclopentadienyl anion is stabilized by the π-holes of both the selenoaldehydic groups, which is appropriately oriented for suitable interactions. The MESP suggests that two selenoaldehyde groups interact with cyclopentadiene anion (Figure 2). To improve the effect

8 ACS Paragon Plus Environment

Page 8 of 28

Page 9 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

of π-hole interactions, compound 3 has been computed with the same level of theory. The cyano group was substituted at the selenoaldehyde carbon centers in 3 (Figure 1). The ∆Hacid calculated for 3 is 317.6 kcal/mol at M06-2X/6-311+G(2d,p) level of theory. This result shows significant enhancement in the acidity compared to parent cyclopentadiene. The influence of the electron withdrawing cyano group in the selenoaldehydic chain has been found to be significant to enhance the acidity of cyclopentadiene unit. The designed compounds 4 and 5 with more cyano groups substituted to the alkyl chains falls in the range of superacid (Table 1). It is known that electron withdrawing group like cyano and nitro groups can stabilize the conjugate base of acid through resonance stabilization and the electron donating group can have an opposing effect on the stabilization of a conjugate base.1, 6,7,9

On the contrary, we have employed electron donating group in the cyclopentadiene unit

to augment further the basicity of selenoaldehyde substituted cyclopentadiene derivative 6.

Figure 1 The structure of 3-(cyclopenta-1,3-dien-1-yl)propaneselenal (1) and it’s hexacyano (5) and methoxy substituted hexacyano (6) compounds and their conjugate bases (1-, 5-, 6-) 9 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

optimized at M06-2X/6-311+G(2d,p) level of theory in gas phase. (distances are given in Å and distance taken between selenium and nearest carbon of Cp ring, grey = carbon, white = hydrogen, blue = nitrogen, red = oxygen, orange = selenium).

The methoxy group was introduced at the R4 position in the cyclopentadiene ring (Scheme 1). The ∆Hacid calculated for 6 with M06-2X/6-311+G(2d,p) level of theory is 279.0 kcal/mol. The compound 6 could act as a stronger super acid compare to strong mineral Brөnsted acids like HNO3, H2SO4 and HClO4. The experimental gas phase ∆Hacid values of these inorganic acids are 324.5, 306.0 and 288.0 kcal/mol, respectively.1 We have introduced the electron donating group rather than electron withdrawing group in the cyclopentadienyl anion unit to increase the electron density in the cyclopentadienyl unit, which in turn would experience stronger attractive interaction with the π-hole of seleno-aldehyde group. On the other hand, if the electron-withdrawing group is attached to the cyclopentadienyl anion then the electron density would deplete in the cyclopentadienyl moiety, and hence the interaction would be weaker with the π-hole of seleno-aldehyde group. To examine this, we have calculated the deprotonation energy (∆Hacid) with electron withdrawing group (-CN) in cyclopentadiene ring of 6 instead of –OMe (6-CN) (Figure S1, Table TS1). The ∆Hacid calculated 6-CN is (284.3 kcal/mol) suggests that -CN substituted cyclopentadienyl anion would be less stable compared to –OMe substituted cyclopentadienyl anion. The calculated MESP for 6 reveals that there is more positive potential around -C=Se region compared to 1 which is neutralized by the negatively charged cyclopentadienyl anion in 6- (Table TS2, Figure 2). The MESP calculation reveals that the positive potential of π-hole shifted toward the Se atom of C=Se on cyano group substitution which facilities the stabilization of conjugate bases (Figure 2). It is important to note that other chalcogens, i.e., -C=O and -C=S can also produce similar stabilizing effect with cyclopentadienyl anion systems. It appears that the longer bond length of –C=Se can approach much closer to the cyclopentadienyl anion compared to its corresponding -C=O and -C=S systems (Figure S1 and Table TS1). We have examined the 10 ACS Paragon Plus Environment

Page 10 of 28

Page 11 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

acidity of 6 using aldehyde and thioaldehyde group (6-C=O and 6-C=S in Figure S1) at the same level of theory (Figure S1 and Table TS1). The calculated results suggest that –C=Se is more suitable to influence the acidity of these cyclopentadiene derivatives (Table 1 and Table TS1). We have calculated free energy change (∆Gacid) for all molecules studied here. The free energy change (∆Gacid) for 6 is lowest compared to other systems examined, which is in good agreement with the (∆Hacid) (Table 1). Table 1 The calculated ∆Hacid and ∆Gacid for compounds (1-6) at M06-2X/6-311+G(2d,p) level of theory in gas phase. (∆Hacid and ∆Gacid are given in kcal/mol) Molecules

∆Hacid

∆Gacid

Molecules

∆Hacid

∆Gacid

1

344.0

337.5

4

299.4

297.0

1-

-

-

4-

-

-

2

340.3

336.0

5

287.0

284.5

2-

-

-

5-

-

-

3

317.6

315.7

6

279.0

275.8

3-

-

-

6-

-

-

11 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 2 MESP calculation of 1, 5, and 6 and their conjugate bases. (Blue and light blue colors indicate the positive potential, and green, yellow, and red colors indicate the negative potential) The Gibbs free energies (∆Gacid) indicate that 6 is more acidic in this series (Table 1 and Table TS1). The Gibbs free energy for selenoaldehyde (6) is also less than that of aldehyde and thioaldehyde substituted in 6-C=O and 6-C=S, respectively (Table TS1). The pKa calculated for the 6 using M06-2X/6-311+G(2d,p) level of theory shows the value of -9.6, whereas the pKa calculated for the reference cyclopentadiene at the same level of theory is 19.0.48 The higher homologue of Se i.e., Te can have an even more significant impact on the acidity of similar derivatives. The calculations performed using telluro aldehydic group with the 5 and 6 show even lower ∆Hacid values than the corresponding C=Se group (Table TS3). These results suggest that the acidity can be modulated with even chalcogens of choice. The stabilization of a conjugate base of an acid with π-hole bonds can have an additional advantage, where the acidity of the system can be modulated with solvents. The non-covalent interactions can be influenced with the change in the environmental conditions.51,52

12 ACS Paragon Plus Environment

Page 12 of 28

Page 13 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 3: The topological analysis of compounds 1, 5, and 6 and their conjugate bases. (orange coloured point indicates the bond critical points (3, -1) appears between two interacting atoms) We have further examined the AIM (Atoms in Molecules) analysis for these systems (Figure 3).47,53 The AIM analysis for 1 reveals that one critical point appears between selenium atom and the carbon atom of Cp- ring after deprotonation, which suggests an interaction between these atoms. The electron density (ρ), Laplacian electron density (∇2ρ) values corroborate the MESP analysis to achieve the stability of the conjugate base of acids (Table TS4).47-53 The V(r)/G(r) ratio indicates that the non-covalent interaction between the C=Se group and the CP- of 1 is purely electrostatic in nature. The electron density values calculated for 6 reveals that the interaction between the C=Se group and the CP- of 6 is predominantly covalent in nature (Figure 3 and Table TS4).53 The calculated bond strength reveals the Se…CP- is stronger in compound 6- (Table TS4). The attraction between the Se group and cyclopentadiene is so strong that these groups fall in the region of the bond formation. The AIM calculated results qualitatively in agreement with the acidity of the systems studied here. The stronger interaction calculated with AIM shows that the conjugate base of that acid is more stable than the conjugate bases of other acids examined (Table 1 and Table TS4). 13 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

We have examined this strategy to examine the stability of dianions of similar systems. We have computed the di anions of cyclopentadiene derivatives.

Figure 4: Structure of di anions of Cyclopentadiene derivatives with selenoaldehyde, optimized at M06-2X/6-311+G(2d,p) level of theory. (∆Hacid is given in kcal/mol) The calculated acidity reveals that ∆Hacid of the parent compound (2Cp in Figure 4) are 349.8 kcal/mol and 467.0 kcal/mol, respectively, which suggests that the 2nd deprotonation would be unlikely in this case. However, the cyclopentadiene derivatives exhibit as super acids with seleno aldehydic group. The ∆Hacid of 1st and 2nd deprotonation of 7 are 291.5 kcal/mol and 353.9 kcal/mol, whereas, with electron donating substituents 8 are 289.6 kcal/mol and 348.4 kcal/mol, respectively (Figure 4). The generation of cyclopentadienyl anion from the cyclopentadiene has been reported in the literature.54 Therefore, the experimental condition maintained to prepare the cyclopentadienyl anion can as well lead to the formation of dianion.

14 ACS Paragon Plus Environment

Page 14 of 28

Page 15 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The binding energies of the first and second electron with the ligand system can also provide the stability of the anions formed.17 The calculations performed for the incorporation of first and second electrons to the ligand molecule (8) using Gaussian 09 program is discussed in the computation section. We have calculated the binding energies (in kcal/mol) of the first and second electron of 8 at M062X/6311+G(2d,p) level of theory. The calculated results suggest that the binding energy of the first electron is ∆E1 (64.3 kcal/mol) and for the second electron is ∆E2 (27.9 kcal/mol), respectively as obtained in other multi-charged anions.17 The stability predicted for the di-anions suggests that such ligands can be potentially used for applications in batteries.16,17 It has been reported that the potentially stable di-anions can be used as lithium and magnesium ion batteries. Therefore, we have extended our study with the dianion 82- to examine its potential as an electrolyte in lithium and magnesium ion batteries. Li-ion batteries play a pivotal role in today’s portable electronics, because they weigh light and possess high energy density.55-57 Li-ion battery comprises of three components such as anode, cathode, and electrolyte.58-62 The Li+ ions transfer from anode to cathode while discharging and reverse when charging and this is supplied by the electrolytes. Research efforts are underway to augment the efficiency, stability and safety of Li-ion batteries.63-70 The electrolytes function as ion carrier between the anode and cathode when current flows through an external circuit. There are electrolytes available consist of lithium salts such as LiAsF6, LiBF4, LiPF6, LiFePO4, LiClO4, LiN-(SO2F)2, and LiN(SO2CF3)2, combined with organic solvents like ethylene carbonate and dimethyl carbonate.71-74 However, there are disadvantages with such electrolytes.32 To resolve the problems associated with the electrolytes known attempts have been made by introducing new solvents or additives.32 One of the important characteristics is it to make the electrolytes halogen free for safety purposes.17,32 The designed dianions may be useful to fulfill the characteristics for Li+ ion batteries.

15 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

These di-anions were also used to examine for their suitability to function as an electrolyte in magnesium ion batteries. The magnesium ion batteries are used as replacement for lithium ion based battery due to their high-oxidative-stability and higher charge density.75 We have calculated the binding energies of metal ions like lithium (Li+) and magnesium (Mg2+) with the stable di-anions 82- using M06-2X/6-31G(2d,p) levels of theory (Figure 5, Table 2).

Figure 5: Structure of metal (Li+ and Mg2+) complexes with compound 82-. The optimized complex structure of Li-82- reveals that the lithium ion interacts with -CN group and the oxygen atom of the methoxy group (Figure 5). The weaker cation…π interaction of the Li+ ion was also observed with the negatively charged 2Cp2- ring (Figure 5). The binding affinity of Li+ with 82- become less as there is weaker cation…π interaction to partially quench the negative charge of anion. However, such cation…π interactions were not observed in Li2-82- complex structure. In Li2-82- complex each lithium ion interacts with CN group and the methoxy oxygen atoms (Figure 5). The interaction of Mg2+ with 82- revealed an interesting trend while complexing (Mg-82-), the magnesium ion interacts with Cp rings rather than complexing with CN and the methoxy group (Figure 5). The cation…π interaction between Mg2+ and Cp rings is weaker as the seleno aldehyde groups partially quench the negative charge of the anion (Figure 5). The calculated binding energies of Li+ and Mg2+ ions with 82- i.e., Li-82-, Li2-82- and Mg-82- at M06-2X/6-31G(2d,p) level of theory in gas phase are given in (Table 2).

16 ACS Paragon Plus Environment

Page 16 of 28

Page 17 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The high performance electrolytes in metal ion batteries should have lower binding affinity towards the metal ions.75 The reports suggest that the dianions such as B12H122- and its cyano derivatives can function as an excellent ligand for lithium and magnesium ion batteries. These studies also suggest that the binding of the metal ions with the anions depend on the size of the ligand anions. The binding energies of metal ions decrease with the increasing size of anions. Here, we have employed a different strategy for weaker binding of metal ions (Li+ and Mg2+) with anions where the anion can be stabilized through the π-hole interaction of seleno-aldehyde group. Therefore, the developed negative charge of the anion is partially quenched with the positive hole of seleno-aldehyde groups. The calculated binding energies for the 1st and 2nd Li+ ion with 82- and 8- are given (Table 2). The binding energies of lithium complexed with 8 (∆ELi1 and ∆ELi2) are 5.8 and 63.6 kcal/mol, respectively using M06-2X/6311+G(2d,p) level of theory. To compare the calculated dissociation energy values for 1st and 2nd Li ion with 82-/8- and the reported B12H122-, additional calculations were performed using B3LYP/6-31+G(d) level of theory in gas phase (Table TS5). The results show that the dissociation energy values are lower compared to the corresponding the lithiated B12H122-. The study extended with the Mg2+ ion with 82- also showed a similar trend as obtained for Li+ ions (Table 2 and Table TS5). These calculated results suggest that the dianion 82- can be suitable to function as an electrolyte in lithium and magnesium ion batteries.16,17 Table 2: Dissociation energies in lithium salts (∆ELi1, ∆ELi2) and magnesium salts (∆EMg) at M06-2X/6-311+G(2d,p) level of theory

Li2-8 Mg-8

M06-2X/6-311+G(2d,p) ∆ELi1 ∆ELi2/∆EMg (kcal/mol) 5.8 63.6 137.9

Conclusions

17 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

In this work, we have reported that π-hole interactions can be exploited to stabilize the anions that can turn an organic acid to super acid. Cyclopentadiene moiety appropriately substituted with selenoaldehydic group 6 can be turned into super acid (∆Hacid = 279.01 kcal/mol). The positive potential on the selenoaldehydic group can effectively interact with the negative region of anions which is formed after deprotonation of selenoaldehyde derivatives. The super acidic behaviour can be augmented with electron donating methoxy group placed to the cyclopentadiene contrary to the conceived opinion of electron withdrawing groups. This strategy has been successfully employed to stabilize the dianions as reports on the stabilization of such species are rare in the literature. The AIM analyses shed light on the nature of non-covalent interactions occurred to stabilize the conjugate base of the substituted cyclopentadiene derivatives. The stability of such anions was further explored for their use in lithium and magnesium ion batteries as this strategy helped to generate the stability of dianions without any influence of environmental conditions. The calculated results suggest that the di-anions can be potentially used for the preparation of lithium and magnesium ion battery electrolytes. The binding energies of lithium complexed with 8 (∆ELi1 and ∆ELi2) are 5.8 and 63.6 kcal/mol respectively using M06-2X/6-311+G(2d,p) level of theory. The stabilization of the negative charge of anions with π-hole interaction of selenoaldehyde groups facilitate the weaker interaction of lithium and magnesium ions with the ligand molecule, which is prerequisite for preparation of electrolyte in effective lithium and magnesium ion batteries. ASSOCIATED CONTENT Supporting Information The AIM calculation and table for ∆Hacid for tellurium substituted cyclopentadiene derivatives are given in SI.

18 ACS Paragon Plus Environment

Page 18 of 28

Page 19 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

AUTHOR INFORMATION Corresponding Author Corresponding Author. Fax: (+91)-278-2567562, E-mail: [email protected]; [email protected] Author Contributions B.G. and M.K.S. designed and performed the study and analyzed data. B.G. and M.K.S. have written the manuscript. ACKNOWLEDGMENT CSIR-CSMCRI registration number– 098/2018. The B.G. thanks DST (New Delhi), DBT, New Delhi, DAE-BRNS, Mumbai) for financial support. M.K.S. acknowledges UGC (New Delhi, India) for an SRF and Academy of Scientific & Innovative Research (AcSIR) for enrolment

in

the

PhD

program.

We

thank

the

reviewers

for

their

valuable

comments/suggestions that have helped us to improve the paper. References: 1. Vianello, R.; Maksic´, Z. B. Strong Acidity of Some Polycyclic Aromatic Compounds Annulated to a Cyclopentadiene Moiety and Their Cyano Derivatives – A Density Functional B3LYP Study. Eur. J. Org. Chem. 2005, 3571-3580. 2. Olah, G. A.; Surya Prakash, G. K.; Molnar, A.; Sommer, J. Superacid Chemistry, 1985, ISBN 978-0-471-59668-4(cloth) 3. Kuck, D. Thermochemical Data of Organic Ions Obtained from Investigations in the More or Less "Diluted" Gas Phase. Angew. Chem. Int. Ed. 2000, 39, 125-130. 4. Alcamı´, M.; Mo´, O.; Ya´nˇez, M. Computational chemistry: a useful (sometimes mandatory) tool in mass spectrometry studies. Mass Spectrom. Rev. 2001, 20, 195-245.

19 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

5. Alcamı´, M.; Mo´ O.; Ya´nˇez, M. Modeling intrinsic basicities and acidities. J. Phys. Org. Chem. 2002, 15, 174-186. 6. Vianello, R.; Maksic, Z. B. Rees polycyanated hydrocarbons and related compounds are extremely powerful Brønsted superacids in the gas-phase and DMSO a density functional B3LYP study. New J. Chem. 2008, 32, 413-427. 7. Si, M. K.; Ganguly, B. A DFT study to design super- and hyperacids with 1-(cyclopenta2,4-dien-1-yl)-4-nitrobenzene

and

3-(cyclopenta-2,4-dien-1-ylmethylene)-6-

methylenecyclohexa-1,4-diene molecules. New J. Chem. 2017, 41, 1425-1429. 8. Vianello, R.; Maksic, Z. B. Polycyano Derivatives of Some Organic Tri- and Hexacyclic Molecules Are Powerful Super- and Hyperacids in the Gas Phase and DMSO: Computational Study by DFT Approach. J. Org. Chem., 2010, 75, 7670-7681. 9. Vianello, R.; Liebman, J. F.; Maksic, Z. B. In search of ultrastrong brønsted neutral organic superacids: A DFT study on some cyclopentadiene derivatives. Chem.–Eur. J., 2004, 10, 5751-5760. 10. Vianello, R.; Maksic, Z. B. Tailoring of strong neutral organic superacids: DFT-B3LYP calculations on some fulvene derivatives. New. J. Chem., 2004, 28, 843–846. 11. Koppel, I. A.; Taft, R. W.; Anvia, F.; Zhu, S.-Z.; Hu, L.-Q.; Sung, K.-S.; DesMarteau, D. D.; Yagupolskii, L. M.; Yagupolskii, Y. L.; Ignat’ev, N. V.; Kondratenko, N. V.; Volkovskii, A. Y.; Vlasov, V. M.; Notario, R.; Maria, P.-C. The Gas-Phase Acidities of Very Strong Neutral Bronsted Acids. J. Am. Chem. Soc., 1994, 116, 3047-3057. 12. Koppel, I. A.; Burk, P.; Koppel, I.; Leito, I. Generalized Principle of Designing Neutral Superstrong Brønsted Acids. J. Am. Chem. Soc., 2002, 124, 5594-5600. 13. Maksic´, Z. B.; Vianello, R. Anionic versus Neutral Protonation—An Ab Initio Analysis with a Triadic Formula. Chem. Phys. Chem., 2002, 3, 696-700.

20 ACS Paragon Plus Environment

Page 20 of 28

Page 21 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

14. Poad, B. L. J.; Reed, N. D.; Hansen, C. S.; Trevitt, A. J.; Blanksby, S. J.; Mackay, E. G.; Sherburn, M. S.; Chan, B.; Radomd, L. Preparation of an ion with the highest calculated proton affinity: ortho-diethynylbenzene dianion. Chem. Sci. 2016, 7, 6245-6250. 15. Saito, G.; Sekizaki, S.; Konsha, A.; Yamochi, H.; Matsumoto, K.; Kusunokid, M.; Sakaguchi,

K.

A

stable

and

flexible

dianion:

2-dicyanomethylene-1,1,3,4,5,5-

hexacyanopentenediide (DHCP2−), and its complex formation. J. Mater. Chem. 2001, 11, 364. 16. Zhao, T.; Zhou, J.; Wang, Q.; Jena, P. Colossal Stability of Gas-Phase Trianions: SuperPnictogens, Angew. Chem. Int. Ed. 2017, 56, 13421-13425. 17. Zhao, H.; Zhou, J.; Jena, P. Stability of B12(CN)122−: Implications for Lithium and Magnesium Ion Batteries. Angew. Chem. Int. Ed. 2016, 55, 3704-3708. 18. Satumov, E. T.; Medvedev, J. J.; Nilov, D. I.; Sandzhieva, M. A.; Boyarskaya, I. A.; Nikolaev, V. A.; Vasilyev, A. V. Protonation and transformations of α-diazo-β-dicarbonyl compounds in superacids: generation of the strongest carbon-centered cationic electrophiles at the protonation of diazomalonates in Friedel–Crafts reactions. Tetrahedron 2016, 72, 48354844. 19. Chakarawet, K.; Knopf, I.; Nava, M.; Jiang, Y.; Stauber, J. M.; Cummins, C. C. Crystalline Metaphosphate Acid Salts: Synthesis in Organic Media, Structures, HydrogenBonding Capability, and Implication of Superacidity. Inorg. Chem. 2016, 55, 6178-6185. 20. Yanai, H.; Taguchi, T. Synthesis of superacidic carbon acid and its derivatives. J. Fluorine Chem. 2015, 174, 108-119. 21. Patel, H. A.; Mansor, N.; Gadipelli, S.; Brett, D. J. L.; Guo, Z. Superacidity in Nafion/MOF Hybrid Membranes Retains Water at Low Humidity to Enhance Proton Conduction for Fuel Cells. ACS Appl. Mater. Interfaces 2016, 8, 30687-30691.

21 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

22. Amani, M.; Burke, R. A.; Ji, X.; Zhao, P.; Lien, D.; Taheri, P.; Ahn, G. H.; Kirya, D.; Ager, J. W.; Yablonovitch, E.; Kong, J.; Dubey, M.; Javey, A. High Luminescence Efficiency in MoS2 Grown by Chemical Vapor Deposition. ACS Nano 2016, 10, 6535-6541. 23. Prajapati, R.; Kohli, K.; Maity, S. K.; Garg, M. O. Slurry Phase Hydrocracking of Residue by Phosphomolybdic and Phosphotungstic Acids. J. Pet. Environ. Biotechnol. 2016, 7, 1000281. 24. Dalpozzo, R.; Bartoli, G., Sambri, L.; Melchiorre, P. Perchloric Acid and Its Salts: Very Powerful Catalysts in Organic Chemistry. Chem. Rev. 2010, 110, 3501-3551. 25. Olah, G. A.; Rasul, G.; York, C.; Prakash, G. K. S. Superacid-Catalyzed Condensation of Benzaldehyde with Benzene. Study of Protonated Benzaldehydes and the Role of Superelectrophilic Activation. J. Am. Chem. Soc. 1995, 117, 11211-11214. 26. Lipping, L.; Leito, I.; Koppel, I.; Koppel, I. A. Gas-Phase Brønsted Superacidity of Some Derivatives of Monocarba-closo-Borates: a Computational Study. J. Phys. Chem. A 2009, 113, 12972-12978. 27. Thayer, A. M. Well-defined catalysts now allow producers to design polymers with exact properties and to create as yet unknown materials. Chem. Eng. News 1995, 73, 15-20. 28. Reed, C. A.; Kim, K. C.; Bolskar, R. D.; Taming, L.; Mueller, J. Superacids: Stabilization of the Fullerene Cations HC60+ and C60·+. Science 2000, 289, 101–104. 29. Reed, C. A.; Kim, K. C.; Stoyanov, E. S.; Stasko, D.; Tham, F. S.; Mueller, L. J.; Boyd, P. D. W. Isolating Benzenium Ion Salts. J. Am. Chem. Soc. 2003, 125, 1796-1804. 30. Reed, C. A. Carboranes:  A New Class of Weakly Coordinating Anions for Strong Electrophiles, Oxidants, and Superacids. Acc. Chem. Res. 1998, 31, 133-139. 31. Valadbeigi, Y. Design of neutral organic superacids using fulvene derivatives with dienol substituent. Int. J. Quantum Chem. 2017, 117, 190-196.

22 ACS Paragon Plus Environment

Page 22 of 28

Page 23 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

32. Giri, S.; Behera, S.; Jena, P. Superhalogens as Building Blocks of Halogen-Free Electrolytes in Lithium-Ion Batteries. Angew. Chem. Int. Ed. 2014, 53, 13916 –13919. 33. Politzer, P.; Murray, J. S. σ-Holes and π-Holes: Similarities and Differences, J. Comp. Chem. 2017, DOI: 10.1002/jcc.24891 34. Politzer, P.; Murray, J. S.; Clark, T. Halogen bonding: an electrostatically-driven highly directional noncovalent interaction. Phys. Chem. Chem. Phys. 2010, 12, 7748-7757. 35. Politzer, P.; Murray, J. S.; Clark, T. Halogen bonding and other σ-hole interactions: a perspective. Phys. Chem. Chem. Phys. 2013, 15, 11178-11189. 36. Murray, J. S.; Lane, P.; Clark, T.; Riley, K. E.; Politzer, P. Σ-holes, π-holes and electrostatically-driven interactions. J. Mol. Model. 2012, 18, 541-548. 37. Martín-Fernández, C.; Merced Montero-Campillo, M.; Alkorta, I.; Yáñez, M.; Mó,c, O.; Elguerob, J. Large Proton Affinity Enhancements Triggered by Non-Covalent Interactions. Chem. Eur. J. 2018, 24, 1971 –1977. 38. Zhao, Y.; Truhlar, D. G. The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals. Theor. Chem. Acc. 2008, 120, 215-241. 39. Woon, D. E.; Dunning Jr., T. H. Gaussian basis sets for use in correlated molecular calculations. III. The atoms aluminium through argon. J. Chem. Phys. 1993, 98, 1358. 40. Frisch, M. J.; pople, J.A.; Binkley, J. S. Self-consistent molecular orbital methods Supplementary functions for Gaussian basis sets. Chem. Phys. 1984, 80, 3265. 41. Walker, M.; Harvey, A. J. A.; Sen, A.; Dessent, C. E. H. Performance of M06, M06-2X, and M06-HF Density Functionals for Conformationally Flexible Anionic Clusters: M06 Functionals Perform Better than B3LYP for a Model System with Dispersion and Ionic Hydrogen-Bonding Interactions. J. Phys. Chem. A 2013, 117, 12590-12600.

23 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

42. Miertuš, S.; Scrocco, E.; Tomasi, J. Electrostatic interaction of a solute with a continuum. A direct utilizaion of AB initio molecular potentials for the prevision of solvent effects. J. Chem phys. 1981, 55, 117. 43. Miertuš, S.; Tomasi, J. Approximate evaluations of the electrostatic free energy and internal energy changes in solution processes. Chem phys. 1982, 65, 239-245. 44. Vianello, R.; Maksić, Z. B. The Acidity of Brønsted CH Acids in DMSO − The Extreme Acidity of Nonacyanocyclononatetraene. Eur. J. Org. Chem. 2004, 2004, 5003-5010. 45. Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al. Gaussian 09, Revision B01, Gaussian, Inc, Wallingford CT 2010. 46. Bader, R. F. W. Atoms in Molecules: A Quantum Theory; Oxford University Press: Oxford, 1990. 47. Bader, R. F. W. Are There Atoms in Molecules? Introduction to Advanced Topics of Computational Chemistry; Universidadde La Habana: La Habana, (2003) 1-28 48. Lu, T.; Chen, F. Multiwfn: A multifunctional wavefunction analyser. J. Comput. Chem. 2012, 33, 580. 49. Clark, T.; Politzer, P.; Murray, J. S. Correct electrostatic treatment of noncovalent interactions: the importance of polarization. WIREs Comput. Mol. Sci. 2015, 5, 169-177. 50. Wang, H.; Wang, W.; Jin, W. J. σ-Hole Bond vs π-Hole Bond: A Comparison Based on Halogen Bond. Chem. Rev. 2016, 116, 5072-5104. 51. Rueping, M.; Nachtsheim, B. J.; Ieawsuwan, W.; Atodiresei, I. Modulating the Acidity: Highly Acidic Brønsted Acids in Asymmetric Catalysis. Angew. Chem. Int. Ed. 2011, 50, 6706-6720.

24 ACS Paragon Plus Environment

Page 24 of 28

Page 25 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

52. Cavalu, S.; Damian, G. Rotational Correlation Times of 3-Carbamoyl-2,2,5,5tetramethyl-3-pyrrolin-1-yloxy Spin Label with Respect to Heme and Nonheme Proteins. Biomacromolecules 2003, 4, 1630-1635. 53. Si, M. K.; Lo, R. N.; Ganguly, B. The origin and magnitude of intramolecular quasicyclic S⋯ O and S⋯ S interactions revisited: A computational study. J. Chem. Phys. Lett. 2015, 631, 6-11. 54. Breton, G. W. Generation and Observation of the Cyclopentadienyl Anion: A Negatively Charged Aromatic Molecule. The Chemical Educator, 1997, 2 (http://journals.springerny.com/chedr). 55. Tarascon, J. M.; Armand, M. Issues and challenges facing rechargeable lithium batteries. Nature 2001, 414, 359-367. 56. Special Issue on Li Insertion into Host Materials (Ed.: D. Aurbach), Electrochim. Acta 1999. 57. Ehrlich, G. M. Handbook of Batteries, 3edrd edMcGraw-Hill, NY, 2002. 58. Gabano, J.-P. Lithium batteries, Academic Press, London, 1983. 59. Nazri, G.-A.; Pistoia, G. Lithium batteries, Science and Technology, Springer, Heidelberg, 2009. 60. Besenhard, J. O. Handbook of Battery Materials, Wiley-VCH, Weinheim, 1999. 61. Broussley, M.; Pistoia, G. Industrial Applications of Batteries, Elsevier, Amsterdam, 2007. 62. Dan, P.; Mengeritsky, E.; Geronov, Y.; Aurbach, D.; Weissman, I. Performances and safety behaviour of rechargeable AA-size Li/LixMnO2 cell. J. Power Sources 1995, 54, 143145.

25 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

63. Wrodnigg, G. H.; Wrodnigg, T. M.; Besenhard, J. O.; Winter, M. Propylene sulfite as film-forming electrolyte additive in lithium ion batteries. Electrochem. Commun. 1999, 1, 148-150. 64. Nakajima, T.; Dan, K.; Koh, M.; Ino, T.; Shimizu, T.; Effect of addition of fluoroethers to organic solvents for lithium ion secondary batteries. J. Fluorine Chem. 2001, 111, 167-174. 65. Schmidt, M.; Heider, U.; Kuehner, A.; Oesten, R.; Jungnitz, M.; Ignatev, N.; Sartori, P. Lithium fluoroalkylphosphates: a new class of conducting salts for electrolytes for high energy lithium-ion batteries. J. Power Sources 2001, 97 – 98, 557-560. 66. Gnanaraj, J. S.; Levi, M. D.; Gofer, Y.; Aurbach, D.; Schmidt, M. LiPF3(CF2 CF3)2 A Salt for Rechargeable Lithium Ion Batteries. J. Electrochem. Soc. 2003, 150, A445-A454. 67. Ein-Eli, Y.; Thomas, S. R.; Koch, V. J. New Electrolyte System for Li-Ion Battery, Electrochem. Soc. 1996, 143, L195-L197. 68. Aurbach, D.; Gofer, Y.; Ben-Zion, M.; Aped, P. The behaviour of lithium electrodes in propylene and ethylene carbonate: Te major factors that influence Li cycling efficiency. J. Electroanal. Chem. 1992, 339, 451-471. 69. Sun, X.; Lee, H. S.; Yang, X. Q.; McBreen, J. A New Additive for Lithium Battery Electrolytes Based on an Alkyl Borate Compound. J. Electrochem. Soc. 2002, 149, A355A359. 70. Aurbach, D.; Gamolsky, K.; Markovsky, B.; Gofer, Y. On the use of vinylene carbonate (VC) as an additive to electrolyte solutions for Li-ion batteries. Electrochim. Acta 2002, 47, 1423-1439. 71. Deepa, M.; Agnihotry, S. A.; Gupta, D.; Chandra, R. Ion-pairing effects and ion– solvent–polymer interactions in LiN(CF3SO2)2–PC–PMMA electrolytes: a FTIR study. Electrochim. Acta 2004, 49, 373-383.

26 ACS Paragon Plus Environment

Page 26 of 28

Page 27 of 28 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

72. Jung, C. Electrochemical absorption effect of BF4 anion salt on SEI layer formation. Solid State Ionics 2008, 179, 1717-1720. 73. Chen, Z.; Lu, W. Q.; Liu, J.; Amine, K. LiPF6/LiBOB blend salt electrolyte for highpower lithium-ion batteries. Electrochim. Acta 2006, 51, 3322-3326. 74. Zhang, Z.; Chen, X.; Li, F.; Lai, Y.; Li, J.; Liu, P.; Wang, X. LiPF6 and lithium oxalyldifluoroborate blend salts electrolyte for LiFePO4/artificial graphite lithium-ion cells. J. Power Sources 2010, 195, 7397-7402. 75. Carter, T. J.; Mohtadi, R.; Arthur, T. S.; Mizuno, F.; Zhang, R.; Shirai, S.; Kampf, J. W. Boron Clusters as Highly Stable Magnesium-Battery Electrolytes. Angew. Chem. Int. Ed. 2014, 53, 3173-3177.

27 ACS Paragon Plus Environment

The Journal of Physical Chemistry 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Graphical Abstract :

28 ACS Paragon Plus Environment

Page 28 of 28