Superatom Compounds under Oriented External Electric Fields

The first hyperpolarizabilities (β0) and natural bond orbital (NBO) charges of ... BF4 and negative fields (F < 0) with the direction of NLi4 ← BF4...
0 downloads 0 Views 1MB Size
Subscriber access provided by UNIVERSITY OF TOLEDO LIBRARIES

C: Plasmonics, Optical Materials, and Hard Matter

Superatom Compounds under Oriented External Electric Fields: Simultaneously Enhanced Bond Energies and Nonlinear Optical Responses Wei-Ming Sun, Chunyan Li, Jie Kang, Di Wu, Ying Li, Bi-Lian Ni, Xiang-Hui Li, and Zhi-Ru Li J. Phys. Chem. C, Just Accepted Manuscript • Publication Date (Web): 27 Mar 2018 Downloaded from http://pubs.acs.org on March 27, 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 25 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

Superatom Compounds under Oriented External Electric Fields: Simultaneously Enhanced Bond Energies and Nonlinear Optical Responses Wei-Ming Suna, Chun-Yan Lia, Jie Kang,a Di Wuc, Ying Lic, Bi-Lian Nia, Xiang-Hui Lib, Zhi-Ru Lic a

The Department of Basic Chemistry, The School of Pharmacy, Fujian Medical University, Fuzhou 350108, People’s Republic of China. b

Medical Technology and Engineering College, Fujian Medical University Fuzhou 350004, People’s Republic of China. c Laboratory of Theoretical and Computational Chemistry, Institute of Theoretical Chemistry, Jilin University, Changchun 130023, People’s Republic of China.

ABSTRACT: It has recently been demonstrated that oriented external electric fields (OEEFs) can serve as future smart reagents and effectors in chemistry. Consequently, the effect of OEEFs on the geometric structures, electronic properties, bonding properties, and NLO responses of three typical superatom compounds, i.e., (NLi4)(BF4) and (BLi6)X (X = BeF3 and BF4) has been systemically investigated by ab initio computations in this work. The computational results reveal that the (NLi4)(BF4) possesses almost the same performance as the traditional alkali metal halide under the same OEEFs. It is found that the vertical ionization potential, bond energy, and first hyperpolarizability of (NLi4)(BF4) can be gradually enhanced by increasing the imposed OEEF from zero to the critical external electric field (Fc) along the charge transfer direction (NLi4→BF4). In particular, the first hyperpolarizability of (NLi4)(BF4) is greatly enlarged from 2.84103 au to 1.36107 au (for 4790 fold) by increasing OEEF from 0 to 121×10−4 au. Moreover, a similar effect is also observed for (BLi6)X. Consequently, this work reveals an effective approach to simultaneously enhance the bond energies and NLO responses of superatom compounds through imposing OEEFs along the charge transfer direction.



Corresponding author.

E-mail address: [email protected] 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

Page 2 of 25

1. INTRODUCTION For more than a century, the Periodic Table of elements continues to serve as the foundation for the synthesis of materials with desired properties. Nevertheless, in a pioneering work, Khanna and Jena1 proposed that atomic clusters with unique electronic structures can mimic the chemical behavior of individual atoms. Hence, they named such unusual clusters as “superatoms” and considered them as a three-dimensional (3D) expansion of the conventional Periodic Table.2,3 Studies in the past decades have greatly expanded the scope of superatoms and yielded plentiful exciting results.4-12 More importantly, it is established that superatoms can maintain their identity, much in the same way as ordinary atoms, when assembled into an extended nanostructure,13 and thereby they offer the exciting prospect of serving as novel building blocks for cluster-assembled compounds and nanomaterials with tailored properties for a great variety of potential technologies.14-17 Up to now, a large number of novel compounds using superatoms as building blocks have been synthesized and characterized by experimental and theoretical workers.4,5,13,18-29 Among them, the “supersalts”,18-24 a special kind of superatom compounds formed by combining superalkalis30 and superhalogens31, have attracted great interests because of their potential applications in nonlinear optical (NLO) materials,19,20 high energy density materials,22 energy storage and energy conversion materials23. For example, Khanna and coworkers have theoretically constructed a new type of unique superatom compound by combining superhalogen Al13 and superalkalis (K3O and Na3O).18 Wu and coworkers have presented the theoretical evidence for a series of superatom compounds, including (BLi6)X (X = F, LiF2, BeF3, and BF4),19 (BF4)M (M = Li, FLi2, OLi3, NLi4),20 (Li3)(SH) (SH = LiF2, BeF3, and BF4)21. Among them, (BLi6)X and (BF4)(NLi4) have been predicted to be stable NLO molecules in view of their large bond energies and considerable first hyperpolarizabilities. In another respect, it is known that the oriented external electric fields (OEEFs) plays an important role in chemistry. For instance, OEEFs can affect the H-bonding and

π-stacking

interactions,32

cause

structural

distortions33

2

ACS Paragon Plus Environment

and

molecular

Page 3 of 25 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

isomerization,34-36 induce proton transfer37–40 and mutation of DNA,41,42 control metal-ligand binding,43,44 elicit spin-crossover transitions45 and spin-polarized conductivity,46,47 promote electron transfer and redox reactions,48–50 and activate click reactions.51 In particular, Shaik and coworkers have successfully utilized OEEFs to catalyze and control a variety of non-redox reactions, such as Diels–Alder reaction,52 enzymatic-like bond-activation reaction,53 hydrogen abstraction reaction,54 and the catalytic cycle of the enzyme cytochrome P450.55 They proposed that an OEEF along the reaction axis will catalyze nonpolar reactions by orders of magnitude, control regioselectivity and induce spin-state selectivity. More recently, the elegant experiment of an OEEF-catalyzed Diels–Alder reaction56 showed that these theoretical predictions are reachable with the help of already accessible techniques.37,56-58 Hence, Shaik et al.59 proposed that the OEEFs have the potential of becoming smart catalysts or inhibitors of non-redox reactions and as controllers of reaction mechanisms. Inspired by the above facts, we want to know whether OEEF can be utilized to modulate the electronic properties and NLO responses of superatom compounds. Consequently, the performance of three superatom compounds, including (NLi4)(BF4) and (BLi6)X (X = BeF3 and BF4), under varying OEEFs has been systemically investigated in the present work. Here, we will mainly solve the following problems: (1) Do the superatom compounds possess the same performance as traditional alkali metal halide under the same OEEFs? (2) How the geometric structures, electronic properties, bonding properties, and NLO responses of superatom compounds evolve with the varying strength of OEEFs? (3) Can OEEF be utilized to enhance the bond energies and NLO responses of superatom compounds? (4) If increasing OEEF can indeed enlarge the bond energies and NLO responses of superatom compounds, does a critical external electric field (Fc, the maximum electric field that the structure can withstand) exist for each superatom compound? (5) What is the detailed mechanism of the enhancement of bond energies and NLO responses of superatom compounds under the chosen OEEFs? 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 25

2. COMPUTATIONAL METHODS The MP2 method has been widely applied to study the structures and properties of superatom compounds.19-21,24 Thus, the equilibrium geometries with all real frequencies of the studied superatom compounds under a series of OEEFs were obtained by the second-order MP2 method in conjunction with the aug-cc-pVDZ basis set. As for the direction of OEEFs, the formal physics convention that positive field is the direction a positive test charge would move was used in this work. To be specific, the keyword “Field = z + n” (n > 0) was used for imposing positive electric field on (NLi4)(BF4), whereas “Field = z - n” (n > 0) was used for applying positive electric field on (BLi6)X (X = BeF3 and BF4). The first hyperpolarizabilities (β0) and natural bond orbital (NBO) charges of these compounds are calculated at the MP2/aug-cc-pVTZ level. Herein, the β0 values were calculated by the finite-field approach, in which the total energy of a molecular system in the presence of a homogeneous electric field can be expressed as follows: 1 1 E  E 0   F   F F   F F F  2 6

(1)

where E0 is the molecular total energy without the electric field, and Fα is the electric field component along α direction; μα, ααβ, and βαβγ are the dipole moment, the polarizability, and the first hyperpolarizability, respectively. Herein, the static mean first hyperpolarizability (β0) are defined as follows:

 (  x 2  y 2  z )2 1/2

(2)

0

3 where  i  (  iii   ijj   ikk ) , i, j, k = x, y, z. 5

In

addition,

the

TD-DFT

calculations

were

performed

at

the

M06-2X/aug-cc-pVTZ level to obtain the oscillator strengths (f0), transition energies (E), and the difference of the dipole moments (μ) between the ground state and crucial excited state for the studied systems. The vertical ionization potential (VIP) of each compound, defined as the total energy difference between the ionized and the neutral compound with the same geometry as the neutral, was obtained at the MP2/aug-cc-pVTZ level. The intramolecular interaction energies (Eint) between the 4

ACS Paragon Plus Environment

Page 5 of 25 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

superalkali and superhalogen subunits of these studied compounds were also calculated at the MP2/aug-cc-pVTZ level. We used the counterpoise procedure60 to eliminate the basis set superposition error (BSSE) effect given as follows,61 Eint  EAB ( X AB )  EA ( X AB )  EB ( X AB )

(3)

where the whole basis set, XAB, was used for the subunit energy (EA and EB) calculations. And then, the bond energies (Eb) of these studied compounds are defined as the negative of Eint. All the above calculations were performed by using the GAUSSIAN 09 program package.62 Dimensional plots of molecular configurations and orbitals were generated with the GaussView program.63

3. RESULTS AND DISCUSSION 3.1 Performance of (NLi4)(BF4) under OEEFs To detect how OEEFs affect the geometric structures, electronic properties, bonding properties, and NLO responses of superatom compounds, the reported global minimum of the supersalt (NLi4)(BF4)20 was firstly chosen as an example and has been optimized both in the absence of external perturbation and under the influence of twenty-one OEEFs with strengths ranging from -80×10−4 to 121×10−4 au (1 au = 5.142  109 V/cm). All the external electric fields have been applied parallel to the charge transfer direction (NLi4 → BF4) in both directions, that is, positive fields (F > 0) with the direction of NLi4 → BF4 and negative fields (F < 0) with the direction of NLi4 ← BF4 (see Figure 1). Our computational test demonstrates that the critical positive electric field (Fc) for (NLi4)(BF4) is 121×10−4 au because a larger positive electric field (≥122×10−4 au) will destroy its geometric structure. Under the above OEEFs, the variations of selective geometric parameters, NBO charges on the superalkali units, the highest occupied molecular orbital (HOMO)−lowest unoccupied molecular orbital (LUMO) gaps, Eb values, Wiberg bond index, VIP, and β0 values of (NLi4)(BF4) were gathered in Table 1, while the rest geometric parameters and total energies are listed in Table S1 of the Supporting Information. It is reported that the 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

superalkali and superhalogen subunits can behave as single alkali metal atoms and halogen atoms, respectively, when they are interacting with each other in the superatom compounds.19,20 Thus, for comparison, the behavior of the traditional alkali-halide LiF under the same OEEFs has also been studied. In the following subsections, the changing trends of geometric features, electronic properties, bonding properties, and NLO responses of (NLi4)(BF4) and LiF under the varying OEEFs are discussed in detail.

Figure 1 Geometric structures of (NLi4)(BF4) and LiF without an electric field. Cartesian axis, atomic numbering, the direction of inherent μz, and the direction of the positive and negative electric fields (F, in 10-4 au) are also shown.

6

ACS Paragon Plus Environment

Page 6 of 25

Page 7 of 25 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

Table 1. Changes of Symmetry Points, Bond Lengths of Li2-F3 (RLi2-F3, in Å), NBO Charges on Superalkali Units (QSA, in |e|), μz (in debye), HOMO-LUMO Gaps (in eV), Vertical Ionization Potentials (VIP, in eV), Bond Energies (Eb, in kcal/mol), Wiberg Bond Index (WBI) of the Li2-F3 Bond, and the First Hyperpolarizability (β0, in au) with the Oriented External Electric Field (OEEF, in 10−4 au) for (NLi4)(BF4). OEEF

Symmetry

RLi2-F3

QSA

μz

Gap

VIP

Eb

WBI

β0

-80 -70 -60 -50 -40 -30 -20 -10 0 10 20 30 40 50 60 70 80 90 100 110 121

C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v

1.985 1.974 1.962 1.952 1.942 1.933 1.924 1.916 1.908 1.899 1.893 1.886 1.879 1.872 1.865 1.859 1.853 1.847 1.841 1.835 1.828

0.866 0.864 0.859 0.861 0.858 0.855 0.852 0.851 0.847 0.847 0.846 0.845 0.843 0.838 0.838 0.835 0.833 0.832 0.834 0.831 0.832

10.200 9.859 9.521 9.189 8.861 8.536 8.213 7.888 7.568 7.245 6.945 6.641 6.315 5.981 5.659 5.333 5.003 4.669 4.328 3.979 3.557

4.21 4.56 4.91 5.27 5.62 5.97 6.31 6.58 6.45 6.14 5.73 5.50 4.79 4.31 3.83 3.34 2.85 2.36 1.88 1.39 0.85

5.55 5.65 5.75 5.85 5.96 6.06 6.16 6.27 6.38 6.48 6.59 6.71 6.82 6.93 7.05 7.16 7.28 7.41 7.53 7.66 7.80

106.6 108.7 110.7 112.6 114.6 116.6 118.5 120.5 122.4 124.4 126.3 128.2 130.1 132.1 134.1 136.1 138.1 140.1 142.2 144.2 146.6

0.0615 0.0626 0.0636 0.0647 0.0656 0.0666 0.0674 0.0682 0.0690 0.0699 0.0705 0.0711 0.0717 0.0724 0.0729 0.0734 0.0737 0.0741 0.0743 0.0744 0.0742

5.68102 7.91102 1.02103 1.26103 1.51103 1.79103 2.10103 2.44103 2.84103 3.31103 3.87103 4.55103 5.38103 6.42103 7.81103 9.72103 1.25104 1.70104 2.55104 4.93104 1.36107

Figure 2. Evolution of (a) the Li2-F3 bond lengths and (b) μz as a function of the oriented external electric field (OEEF, in 10-4 au) for the (NLi4)(BF4) and LiF (inset) compounds. 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

A. Geometric Characteristics Figure 1 displays the optimized C3v-symmetric (NLi4)(BF4) without any OEEF at the MP2/aug-cc-pVDZ level, where the bond length of Li2-F3 is 1.908 Å and the angles of θ and φ are 93.47°and 111.15°, respectively. Under different OEEFs, this structure preserves its C3v geometry, where the superalkali and superhalogen subunits retain their structural integrities. Even so, the applied OEEFs also have certain influence on the geometric structure of (NLi4)(BF4), which has been clearly displayed in Figures 2a and S2. From Figure S2, it can be observed that the bond lengths of N–Li2 and B-F3 of (NLi4)(BF4) are slightly elongated, whereas the Li2-F3 and B-F4 bonds are gradually reduced with the OEEF changing from -80×10−4 to 121×10−4 au. It is found that, under the external electric fields, the variation of bond length is more obvious for Li2-F3 (ca. 0.17 Å), while the change in bond length is less than 0.07 Å for the other bonds, indicating that the OEEFs have substantial effects on the Li-F bonds formed between the NLi4 and BF4 subunits of (NLi4)(BF4). As is revealed in Figure 2a, the Li-F bond of LiF exhibits almost the same dependence on OEEF. To be specific, the Li2-F3 and Li-F bonds are gradually shortened under the increasing positive OEEF along the charge transfer direction. This is reasonable considering the fact that the F− anion tends to move towards the positive pole while the Li+ cation tends to move towards the negative pole of the applied electric field, leading to the contraction of the bridging bonds between the (super)halogen anion and (super)alkali cation in the ionic compounds. When taking the other geometric parameters into account, as shown in Figure S2, it is found that the distance between the B and N atoms of (NLi4)(BF4) is generally increased with the applied OEEF changing from -80×10−4 to 121×10−4 au. This can be understood by the fact that the B center of BF4-unit carries positive charge and the N center of NLi4-unit has negative charge, which gradually makes them move to the opposite direction as the imposed OEEF varies from -80×10−4 to 121×10−4 au. As the Li2-F3 bonds are shortened, the elongation of B-N distances leads to the synchronous stretch of the N–Li2 and B-F3 bonds. Meanwhile, the angle of θ is gradually reduced, 8

ACS Paragon Plus Environment

Page 8 of 25

Page 9 of 25 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

whereas the angle of φ is slightly enlarged under the variation of OEEF. B. Electronic Properties To explore the effect of OEEFs on the charge transfer in (NLi4)(BF4) and LiF compounds, NBO analyses were performed and the related results are listed in Tables 1, S2 and S3 and plotted in Figure S3. From Table 1 and Figure S3a, it can be seen that the general changing trend of NBO charges on the NLi4 unit (QSA) of (NLi4)(BF4) is decreased from 0.866 to 0.831|e| with the OEEF changing from -80×10−4 to 121×10−4 au, which is similar to the varying tendency of QLi for LiF (see Figure S3b). This indicates that the charge transfer from (super)alkali to (super)halogen in these ionic compounds is gradually weakened by enlarging positive electric fields along the charge transfer direction, but is gradually enhanced by increasing the negative electric fields with the opposite direction. Moreover, from Table S2, it is observed that the NBO charges on the Li1 and F4 atoms of (NLi4)(BF4) change more obviously as compared with the other atoms in each superatom subunit. In addition, our computations reveal that only the μz component makes contribution to the dipole moment of (NLi4)(BF4) or LiF, as the μx and μy are zero. As shown in Figure 1, the inherent μz points from (super)halogen to (super)alkali with the same direction as the negative field. Consequently, as shown in Table 1 and Figure 2b, the μz of (NLi4)(BF4) is gradually enlarged from 7.568 to 10.200 D under the increasing negative OEEF, whereas it is gradually reduced to 3.557 D under the increasing positive field. As a result, an excellent linear relationship between μz and electric field (F) was obtained, and the resulting linear equation is μz = -0.032 F + 7.580 with R2 = 0.99978 (F in 10-4 au). The same case is also true for the LiF compound, whose μz decreases from 6.994 to 6.338 D with the linear equation of μz = -0.0032 F + 6.723 (F in 10-4 au, R2 = 0.99933) under the same OEEFs. This phenomenon can be attributed to the fact that the increasing positive electric field not only gradually shortens the Li-F bonds but also weakens the charge transfer from Li to F (see Figures 2a and S3b), jointly leading to the decrease of μz as the magnitude of dipole moment is equal to the separated charge multiplied by the distance between the 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

charges (μ = q· d). The HOMO–LUMO gap of cluster is considered to be an important electronic property in chemistry, which can represent its ability to participate in chemical reactions to some degree. Moreover, the HOMO–LUMO gap is also closely related to the band gap in solid-state physics, serving as an important indicator to predict the conductivity of cluster-assembled materials. As shown in Table 1, the HOMO-LUMO gap of (NLi4)(BF4) without external electric field is as large as 6.45 eV, exhibiting its intrinsic insulator properties. The dependence of HOMO–LUMO gaps of (NLi4)(BF4) and LiF compounds on OEEF is plotted in Figure S3. From Figure S3c, it is observed that the HOMO–LUMO gap of (NLi4)(BF4) increases from 4.21 to 6.58 eV with the negative OEEF decreasing from -80×10−4 to -10×10−4 au, and then gradually reduces from 6.45 to 0.58 eV with the increasing positive electric field from 0 to 121×10−4 au. This variation of HOMO–LUMO gap suggests that our work may provide an effective method to control the energy gaps of superatom compounds. Similar varying trend of HOMO–LUMO gap is also observed for LiF, but the difference is that the peak value of gap appears at OEEF = -40×10−4 au for LiF (see Figure S3d). To deeply understand the origin of the change of energy gaps, the variations of HOMO and LUMO energies of (NLi4)(BF4) as a function of OEEF are shown in Figure S4 of Supporting Information. From this figure, it is observed that, regardless of which direction of the OEEF is adopted, the increasing electric field strength leads to the quickly decrease of LUMO level. As the HOMO level decreases slowly and linearly under the same OEEFs, the change tendency of the HOMO-LUMO gap of (NLi4)(BF4) accords well with the variation of LUMO level. In other words, the change rule of HOMO-LUMO gap of (NLi4)(BF4) mainly depends on the evolution of LUMO level under the varying OEEFs. Among various electronic properties, the vertical ionization potential (VIP) is considered to be an important parameter for estimating the electronic stability of clusters. Thus, the VIP values of (NLi4)(BF4) and LiF compounds under the chosen OEEFs have been also calculated and are listed in Tables 1 and S3, respectively. To 10

ACS Paragon Plus Environment

Page 10 of 25

Page 11 of 25 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

better visualize the results, the dependences of VIP on OEEF for (NLi4)(BF4) and LiF are also plotted in Figure 3. It noted that the VIP values of (NLi4)(BF4) and LiF linearly increase with the imposed OEEF changing from -80×10−4 to 121×10−4 au, demonstrating that the electronic stability of such compounds can be gradually enhanced by applying increased positive electric fields along the charge transfer direction.

Figure 3. Dependences of (a) VIP and (b) Eb on the oriented external electric field (OEEF, in 10-4 au) for the (NLi4)(BF4) and LiF (inset) compounds.

C. Bonding Properties The bond energy (Eb) is considered to a crucial parameter to evaluate the strength of electrostatic interaction between the cation and anion in ionic superatom compounds.19,20 Herein, the dependences of the Eb values of (NLi4)(BF4) and LiF compounds on OEEF are plotted in Figure 3b. It is observed that the Eb value of (NLi4)(BF4) gradually increases with the imposed OEEF changing from -80×10−4 to 121×10−4 au, and an excellent linear relationship between Eb and electric field (F) is obtained with the linear equation of Eb = 0.197 F + 122.417 (R2 = 0.99989 and F in 10-4 au). This indicates that the strength of electrostatic interaction between the NLi4+ and BF4− subunits in (NLi4)(BF4) can be enhanced by employing a larger electric field along the charge transfer direction, resulting in a more stable supersalt where the superalkali and superhalogen subunits are more tightly bound. From the inset picture of Figure 3b, the same case is also observed for the traditional LiF compound. 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

How to understand the effect of OEEFs on the electrostatic interaction between the NLi4+ and BF4− units in the (NLi4)(BF4) compound? This may be related to the fact that the Li-F bonds linking the NLi4+ and BF4− subunits together in (NLi4)(BF4) are gradually shortened with changing OEEF from -80×10−4 to 121×10−4 au (see Figure 2a). According to the Coulomb’s law,64 the electrostatic interaction between two point charges is inversely proportional to the square of the distance between them. Hence, the electrostatic interaction between the Li+ and F− of the bridging Li-F bonds in (NLi4)(BF4) gradually increases when the Li-F bond lengths decrease. Also, this can be further confirmed by the increasing Wiberg bond index (WBI) of bridging Li-F bonds under the changing OEEF from -80×10−4 to 121×10−4 au (see Figure S3e). Similarly, the dependence of Eb of the LiF compound on OEEF can also be understood by the decreasing Li-F bonds and the increasing WBI values (see Figures 2 and S3f). D. NLO Responses It is reported that the (NLi4)(BF4) compound exhibits relatively large NLO response among the BF4-M (M = Li, FLi2, OLi3, and NLi4) series.20 However, the field-free first hyperpolarizability (β0) of 2.84103 au for (NLi4)(BF4) is not large enough; this suggests that this compound can not serve as a good NLO molecule in the absence of an external electric field. One interesting question therefore emerges: can the NLO response of this compound be enhanced by imposing appropriate OEEFs? Consequently, the β0 values of (NLi4)(BF4) under the chosen OEEFs have been calculated and listed in Table 1. To better visualize the results, the dependence of the β0 value of (NLi4)(BF4) on electric field is exhibited in Figure 4. As shown in Figure 4, the β0 value of (NLi4)(BF4) gradually increases along with the imposed OEEF changing from -80×10−4 to 121×10−4 au, indicating that the NLO response of (NLi4)(BF4) can be modulated by changing the magnitude and direction of an applied OEEF. To be specific, the field-free β0 value (2.84103 au) of (NLi4)(BF4) is greatly enlarged to 1.36107 au (4790 fold increase) via imposing an OEEF of 121×10−4 au along the charge transfer direction, whereas it is reduced to 12

ACS Paragon Plus Environment

Page 12 of 25

Page 13 of 25 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

5.68102 au under OEEF of -80×10−4 au against the charge transfer direction. It is also noted that, under positive electric fields, the β0 of (NLi4)(BF4) is slowly increased in the range of low fields, but is quickly enlarged in the range of high fields, indicating that the enhanced effect of positive fields on β0 is more obvious for the larger electric fields. This case is quite different to LiF, for which β0 almost linearly increases with the electric field changes from -80×10−4 to 121×10−4 au (see Figure S5).

Figure 4. Dependence of the first hyperpolarizability (β0) on the oriented external electric field (OEEF, in 10-4 au) for the (NLi4)(BF4) compound. The inset is the change of β0 under the fields ranging from 115×10−4 to 121×10−4 au with the step size of 1×10−4 au.

As shown in Table 1, the β0 of (NLi4)(BF4) is sharply increased from 4.93104 to 1.36107 au when the imposed OEEF is changed from 110×10−4 au to the critical electric field (121×10−4 au). Thus, we focused our attention on the variation of β0 values under this OEEF range. The β0 values of (NLi4)(BF4) for fields ranging from 115×10−4 to 121×10−4 au (with a step size of 1×10−4 au) have been calculated and are shown in Table S4 and Figure 4. Similar to the above case, the β0 of (NLi4)(BF4) is slowly increased in the range of low fields (115×10−4~120×10−4 au), but is greatly enlarged from 120×10−4 to 121×10−4 au, indicating that the enhanced effect of electric field on β0 is gradually enlarged with the increase of OEEF. Additionally, it can be 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

seen from Table S4 that there is a prominent tensor component (βz) that makes the dominant contribution to the β0 values of (NLi4)(BF4) under the applied fields ranging from 0 to 117×10−4 au. In contrast, when the applied fields is larger than 117×10−4 au, more than one component contributes to the β0 value. How to understand the dependence of β0 value on the chosen OEEFs? We may find some clues from the two-level model proposed by Oudar and Chemla65,66. For the static case, the expression of the two-level model is employed as follows:

0 

  f 0 E 3

where ΔE, f0, and Δµ are the transition energy, oscillator strength, and the difference in the dipole moment between the ground state and crucial excited state, respectively. According to the two-level expression, β0 is proportional to f0 and Δµ, whereas is inversely proportional to the third power of ΔE and thus ΔE is usually considered to be the decisive factor on β0. Herein, the ΔE, f0, and Δµ values of (NLi4)(BF4) under the applied electric fields are presented in Table S4 and the change trend of ΔE is illustrated in Figure S6. As shown in Figure S6, it can be clearly seen that the ΔE values of (NLi4)(BF4) exhibit a decreased trend as the OEEF changes from -80×10−4 to 121×10−4 au, which justifies the increasing tendency of corresponding β0 values. Moreover, the greatly large Δµ value of 50.273 D at OEEF = 121×10−4 au can also contribute to the remarkably large β0 of (NLi4)(BF4) under this electric field. Accordingly, the combination of the smallest ΔE and largest Δµ is responsible for the largest β0 value of (NLi4)(BF4) under its critical electric field. 3.2 Performance of (BLi6)X (X = BeF3 and BF4) under OEEFs Based on the above analyses, we gained a significative conclusion that imposing an increasing OEEF along the charge transfer direction of (NLi4)(BF4) is an effective approach to simultaneously enhance its bond energy and first hyperpolarizability. Can this conclusion be extended to the other superatom compounds? The performance of (BLi6)X (X = BeF3 and BF4) under the OEEFs ranging from 0 to Fc along the charge transfer direction has therefore also been investigated. The Fc values of (BLi6)(BeF3) 14

ACS Paragon Plus Environment

Page 14 of 25

Page 15 of 25 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

and (BLi6)(BF4) are found to be 95×10−4 and 93×10−4 au, respectively. The geometric structures of (BLi6)X (X = BeF3 and BF4) are shown in Figure 5, and the geometric parameters and relevant properties are summarized in Tables 2 and S5.

Figure 5 Geometric structures of (BLi6)X (X = BeF3 and BF4) without an external electric field. Cartesian axis, atomic numbering, the direction of inherent μz, and the direction of applied electric field (F, in 10-4 au) are also shown.

Computationl results reveal that both of (BLi6)(BeF3) and (BLi6)(BF4) preserve their C3v-symmetric structures under the chosen external electric fields, where the superalkali and superhalogen subunits always retain their structural integrities under the varying OEEF. The evolution of the geometric parameters of (BLi6)X (X = BeF3 and BF4) under the varying OEEFs has been plotted in Figures S7 and S8, respectively. Similar to (NLi4)(BF4), the Li2-F3 bonds linking the superalkali and superhalogen units together in (BLi6)X are gradually shortened with the increasing electric fields. Also, it is found that the B-Be and B-B distances of (BLi6)(BeF3) and (BLi6)(BF4) are generally increased with the increase of OEEF. Moreover, it can be observed from Figures S7 and S8 that the structural changes of BLi6 unit under the varying OEEFs are the same for (BLi6)X. That is, regardless of whether the superhalogen unit is BeF3 or BF4 for (BLi6)X, the B-Li2 bond length is gradually increased, whereas the bond length of B-Li1 and the angle of θ are reduced as the electric field increases.

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

Page 16 of 25

Table 2. Changes of Symmetry Points, Bond Lengths of Li2-F3 (RLi2-F3, in au), μz (in debye), NBO Charges on Superalkali Units (QSA, in |e|), HOMO-LUMO Gaps (in eV), Bond Energies (Eb, in kcal/mol), Wiberg Bond Index (WBI), Vertical Ionization Potentials (VIP, in eV), and the First Hyperpolarizability (β0, in au) with Oriented External Electric Field (OEEF, in 10−4 au) for (BLi6)X (X = BeF3 and BF4). Compounds (BLi6)(BeF3)

OEEF 0 10 20 30 40 50 60 70 80 95

Symmetry C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v

RLi2-F3 1.814 1.813 1.813 1.812 1.810 1.808 1.803 1.803 1.802 1.795

μz -5.083 -3.931 -2.744 -1.565 -0.378 0.8572 2.1459 3.5939 5.5795 21.277

QSA 0.883 0.871 0.861 0.854 0.846 0.841 0.838 0.830 0.827 0.817

Gap 5.10 4.85 4.43 3.98 3.53 3.09 2.64 2.21 1.82 1.56

Eb 132.4 134.6 137.2 140.1 143.1 146.3 149.4 152.9 156.3 161.9

WBI 0.0789 0.0786 0.0786 0.0786 0.0788 0.0792 0.0792 0.0797 0.0805 0.0823

VIP 5.54 5.51 5.63 5.72 5.82 5.92 6.04 6.13 6.25 6.43

β0

(BLi6)(BF4)

0 10 20 30 40 50 60 70 80 93

C3v C3v C3v C3v C3v C3v C3v C3v C3v C3v

1.881 1.874 1.865 1.858 1.851 1.845 1.839 1.833 1.829 1.820

-5.818 -4.925 -3.991 -3.025 -2.016 -0.965 0.163 1.4143 3.0009 11.616

0.819 0.819 0.817 0.815 0.815 0.814 0.813 0.812 0.811 0.807

5.03 4.87 4.52 4.09 3.66 3.22 2.79 2.37 1.97 1.60

121.1 124.5 127.8 131.1 134.5 137.9 141.3 144.8 148.4 153.1

0.0736 0.0737 0.0739 0.0741 0.0742 0.0744 0.0746 0.0749 0.0754 0.0764

5.48 5.59 5.69 5.84 5.98 6.12 6.26 6.41 6.56 6.76

5.93103 9.01103 1.35104 1.99104 2.87104 4.12104 6.21104 1.10105 2.56105 4.15105

16

ACS Paragon Plus Environment

1.01104 1.53104 2.13104 2.92104 3.96104 5.41104 8.41104 1.65105 3.63105 1.47106

Page 17 of 25 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 NBO computations show that, similar to (NLi4)(BF4), the NBO charges on the BLi6 subunits of (BLi6)X decrease with the increasing OEEF (see Figure S9a), indicating that the charge transfer from superalkali to superhalogen in (BLi6)X can be weakened by exposing them to electric fields along the charge transfer direction. It is noted that the decreasing tendency of QSA is more obvious for (BLi6)(BeF3) as compared with (BLi6)(BF4), implying that the effect of electric field on the charge transfer of the former is much larger than that of the latter. As for μz, a similar variation is observed for the (BLi6)X compounds (see Figure S9b). To be specific, the magnitudes of μz for (BLi6)X are firstly reduced with the increasing OEEF until the directions of inherent μz are reversed, and then they are gradually increased as the electric field increases. In particular, it can be seen from Table 2 that the μz values of (BLi6)X are enlarged by ca. 4 times with the OEEF changing from 80×10−4 au to their respective Fc, whereas they are almost linearly increased under the lower OEEFs from 0 to 80×10−4 au (see Figure S9b).

Figure 6. Dependence of (a) VIP and (b) Eb on the oriented external electric field (F, in 10-4 au) for the (BLi6)X (X = BeF3 and BF4) compounds.

From Figure S9c, it is observed that the HOMO–LUMO gaps of (BLi6)X gradually decrease with the increasing electric fields, suggesting that their energy gaps can also be manipulated by changing the imposed OEEF. When the variations of HOMO and LUMO energies of (BLi6)X are taken into consideration, it can be seen from Figure S10 that both of HOMO and LUMO energy levels are decreased with the increasing OEEF, but the decrease of the LUMO is much obvious than that of the 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

HOMO, leading to reduced tendency of the HOMO-LUMO gaps. It is worth mentioning that the gaps of (BLi6)(BeF3) and (BLi6)(BF4) are 1.56 and 1.60 eV, respectively, even when they are exposed to their Fc. These values are comparable to the experimental HOMO−LUMO gap of 1.57 eV for the kinetically stable C60,67 suggesting that the chemical stabilities of (BLi6)(BeF3) and (BLi6)(BF4) continue to be large under the Fc. In addition, as shown in Table 2 and Figure 6a, the calculated VIP values of (BLi6)X (X = BeF3 and BF4) increase with increasing electric field with only one exception at OEEF = 10×10−4 au for (BLi6)(BeF3); this further demonstrates that the electronic stability of such compounds can be enhanced by the application of electric field along the charge transfer direction. Turn to the bond energy (Eb), it is found that this parameter almost linearly increases with the increase of imposed OEEF for (BLi6)(BeF3) and (BLi6)(BF4), indicating that the strength of ionic bonds between the BLi6+ and BeF3−/BF4− subunits in these two compounds can be gradually enhanced by employing an increasing OEEF along the same direction as charge transfer. This is further confirmed by the generally increasing trend of Wiberg bond indexes (WBI) of the bridging Li2-F3 bonds under the increasing OEEF (see Figure S9d). Thus, this work provides an effective strategy to enhance the interaction between cationic superalkali and anionic superhalogen units in such supersalts.

Figure 7. Dependence of the first hyperpolarizability (β0) on the oriented external electric field (F, in 10-4 au) for the (BLi6)X (X = BeF3 and BF4) compounds.

18

ACS Paragon Plus Environment

Page 18 of 25

Page 19 of 25 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

We now turn to the key question of this study: can the NLO responses of (BLi6)X (X = BeF3 and BF4) also be enhanced by exposing them to OEEFs? Consequently, the β0 values of (BLi6)X under the varying OEEFs have been obtained and listed in Table 2. To clearly display the results, the relationship between β0 and OEEF for (BLi6)X is also illustrated in Figure 7. It is notable that the β0 values of (BLi6)X gradually increase with the increasing OEEF, indicating that their NLO responses can be substantially enhanced by introducing an OEEF along the charge transfer direction. Here, the increased trend of β0 values for (BLi6)X can also be well understood by the corresponding decrease in ∆E values when the applied electric field increases (see Table S6 and Figure S11); this demonstrates that the ∆E is the decisive factor determining the NLO responses of such superatom compounds. Similar to (NLi4)(BF4), the β0 values of (BLi6)X also slowly increase under the lower electric fields, while quickly increasing under relatively higher fields. This further confirms the conclusion that the enhanced effect of OEEFs on β0 is more obvious for larger electric fields. In particular, as shown in Table 2 and Figure 7, the β0 of (BLi6)(BeF3) is greatly increased from 1.01104 to 1.47106 au (146 fold), while the β0 of (BLi6)(BF4) is enlarged from 5.93103 to 4.15105 au (70 fold) when the imposed electric field increases from 0 to the critical electric field (Fc). Additionally, it can be seen from Table S6 that the tensor component βz always makes the dominant contribution to the β0 values of (BLi6)X under the applied OEEFs. 4. CONCLUSIONS In the present work, the evolution of geometric structures, electronic properties, bonding properties, and NLO responses of (NLi4)(BF4) and (BLi6)X (X = BeF3 and BF4) under the varying OEEFs has been studied in detail by ab initio computations. Our results demonstrate that the OEEFs have substantial effect on the geometric and electronic properties, bonding properties, and NLO responses of the (NLi4)(BF4), and that these effects are almost the same as for traditional alkali metal halides. It is found that the HOMO–LUMO gaps of the studied compounds can be modulated by changing the magnitude and direction of OEEFs. In particular, the bond energies and 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

first hyperpolarizabilities can be gradually enhanced by increasing OEEF along the charge transfer direction from zero to the critical external electric field (Fc) along the charge transfer direction. This indicates that imposing an increasing OEEF along the charge transfer direction of superatom compounds is an effective approach to simultaneously enhance their bond energies and NLO responses. Despite that imposing extremely high electric field on clusters is a great challenge, it is hoped that this work could not only offer theoretical insights into the effect of OEEFs on the structures, electronic properties, and NLO responses of superatom compounds, but also attract more research interests and efforts of experimentalist in the application of OEEFs to modulate the desired properties of superatom compounds.

Supporting Information Effect of OEEFs on geometric parameters, β0, f0, ΔE, and Δμ for (NLi4)(BF4) and (BLi6)X. Charge distribution and WBI within NLi4+ and BF4− of (NLi4)(BF4). Dependence of QSA, HOMO-LUMO gaps, and WBI on OEEF for LiF, (NLi4)(BF4), and (BLi6)X. Dependence of HOMO and LUMO energy levels as well as ∆E on OEEF for (NLi4)(BF4) and (BLi6)X. Dependence of β0 on OEEF as well as the changes of symmetry points, bond length, μz, QLi, HOMO-LUMO gaps, VIP, Eb, WBI, and β0, with OEEF for LiF. Optimized structures of BeF3−, BF4−, NLi4+, and BLi6+. Cartesian coordinates of field-free (NLi4)(BF4) and (BLi6)X.

Acknowledgements This work was supported by the National Natural Science Foundation of China (21303066, 21375017, 21573089, 21603032), Natural Science Foundation of Fujian Province (2016J05032, 2016J01771), State Key Development Program for Basic Research of China (2013CB834801), Academic Foundation for Professor of Fujian Medical University (JS14008, JS14009). The authors also acknowledge the National Supercomputing Center in Shenzhen for providing computational resources.

References (1) Khanna, S.; Jena, P. Atomic Clusters: Building Blocks for a Class of Solids. Phys. Rev. B 1995, 51, 13705-13716. (2) Castleman, A. W. From Elements to Clusters: The Periodic Table Revisited. J. Phys. Chem. Lett. 2011, 2, 1062-1069. (3) Jena, P. Beyond the Periodic Table of Elements: The Role of Superatoms. J. Phys. Chem. Lett. 2013, 4, 1432-1442. (4) Bergeron, D. E.; Castleman, A. W.; Morisato, T.; Khanna, S. N. Formation of Al 13I-: Evidence for the Superhalogen Character of Al13. Science 2004, 304, 84-87. (5) Bergeron, D. E.; Roach, P. J.; Castleman, A. W.; Jones, N. O.; Khanna, S. N. Al Cluster 20

ACS Paragon Plus Environment

Page 20 of 25

Page 21 of 25 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

Superatoms as Halogens in Polyhalides and as Alkaline Earths in Iodide Salts. Science 2005, 307, 231-235. (6) Reveles, J. U.; Khanna, S. N.; Roach, P. J.; Castleman, A. W. Multiple Valence Superatoms. Proc. Natl. Acad. Sci. U. S. A. 2006, 103, 18405-18410. (7) Reveles, J. U.; Clayborne, P. A.; Reber, A. C.; Khanna, S. N.; Pradhan, K.; Sen, P.; Pederson, M. R. Designer Magnetic Superatoms. Nat. Chem. 2009, 1, 310-315. (8) Sun, W.-M.; Li, Y.; Wu, D.; Li, Z.-R. Designing Aromatic Superatoms. J. Phys. Chem. C 2013, 117, 24618-24624. (9) Sun, W.-M.; Wu, D.; Li, X.-H.; Li, Y.; Chen, J.-H.; Li, C.; Liu, J.-Y.; Li, Z.-R. Quasi-chalcogen Characteristics of Al12Be: A New Member of Three-dimensional Periodic Table. J. Phys. Chem. C 2016, 120, 2464-2471. (10) Sun, W.-M.; Li, X.-H.; Li, Y.; Liu, J.-Y.; Wu, D.; Li, C.-Y.; Ni, B.-L.; Li, Z.-R. On the Feasibility of Designing Hyperalkali Cations Using Superalkali Clusters as Ligands. J. Chem. Phys. 2016, 145, 194303. (11) Hou, N.; Wu, D.; Li, Y.; Li, Z. R. Lower the Electron Affinity by Halogenation: An Unusual Strategy to Design Superalkali Cations. J. Am. Chem. Soc. 2014, 136, 2921-2927. (12) Sun, W. M.; Hou, D.; Wu, D.; Li, X. H.; Li, Y.; Chen, J. H.; Li, C. Y.; Li, Z. R. Theoretical Characterization of a Series of N5-based Aromatic Hyperhalogen Anions. Dalton Trans. 2015, 44, 19901-19908. (13) Castleman, A. W.; Khanna, S. N. Clusters, Superatoms, and Building Blocks of New Materials. J. Phys. Chem. C 2009, 113, 2664–2675. (14) CastlemanJr., A. W.; Khanna, S. N.; Sen, A.; Reber, A. C.; Qian, M.; Davis, K. M.; Peppernick, S. J.; Ugrinov, A.; Merritt, M. D. From Designer Clusters to Synthetic Crystalline Nanoassemblies. Nano Lett. 2007, 7, 2734-2741. (15) Claridge, S. A.; A. W. Castleman, J.; Khanna, S. N.; Murray, C. B.; Sen, A.; Weiss, P. S. Cluster-Assembled Materials. ACS Nano 2009, 3, 244–255. (16) Qian, M.; Reber, A. C.; Ugrinov, A.; Chaki, N. K.; Mandal, S.; H. M. Saavedra; Khanna, S. N.; Sen, A.; Weiss, P. S. Cluster-Assembled Materials: Toward Nanomaterials with Precise Control over Properties. ACS Nano 2010, 4, 235–240. (17) Mandal, S.; Reber, A. C.; Qian, M.; Weiss, P. S.; Khanna, S. N.; Sen, A. Controlling the Band Gap Energy of Cluster-Assembled Materials. Acc. Chem. Res. 2013, 46, 2385-95.. (18) Reber, A. C.; Khanna, S. N.; Castleman, A. W. Superatom Compounds, Clusters, and Assemblies: Ultra Alkali Motifs and Architectures. J. Am. Chem. Soc. 2007, 129, 10189-10194. (19) Li, Y.; Wu, D.; Li, Z. R. Compounds of Superatom Clusters: Preferred Structures and Significant Nonlinear Optical Properties of the BLi6-X (X = F, LiF2, BeF3, BF4) Motifs. Inorg. Chem. 2008, 47, 9773-9778. (20) Yang, H.; Li, Y.; Wu, D.; Li, Z.-r. Structural Properties and Nonlinear Optical Responses of Superatom Compounds BF4-M (M = Li, FLi2, OLi3, NLi4). Int. J. Quantum Chem. 2012, 112, 770-778. (21) Wang, F. F.; Li, Z. R.; Wu, D.; Sun, X. Y.; Chen, W.; Li, Y.; Sun, C. C. Novel Superalkali Superhalogen Compounds (Li3)+(SH)- (SH = LiF2, BeF3, and BF4) with Aromaticity: New Electrides and Alkalides. Chemphyschem 2006, 7, 1136-1141. (22) Behera, S.; Samanta, D.; Jena, P. Nitrate Superhalogens as Building Blocks of Hypersalts. J. Phys. Chem. A 2013, 117, 5428-5434. (23) Giri, S.; Behera, S.; Jena, P. Superalkalis and Superhalogens as Building Blocks of 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

Supersalts. J. Phys. Chem. A 2014, 118, 638-645. (24) Srivastava, A. K.; Misra, N. Novel Li3X3 Supersalts (X = F, Cl, Br & I) and Their Alkalide Characteristics New J. Chem. 2014, 38, 2890-2893. (25) Sun, W. M.; Fan, L. T.; Li, Y.; Liu, J. Y.; Wu, D.; Li, Z. R. On the Potential Application of Superalkali Clusters in Designing Novel Alkalides with Large Nonlinear Optical Properties. Inorg. Chem. 2014, 53, 6170-6178. (26) Sun, W. M.; Li, Y.; Li, X. H.; Wu, D.; He, H. M.; Li, C. Y.; Chen, J. H.; Li, Z. R. Stability and Nonlinear Optical Response of Alkalides that Contain a Completely Encapsulated Superalkali Cluster. Chemphyschem 2016, 17, 2672-2678. (27) Wang, J. J.; Zhou, Z. J.; Bai, Y.; Liu, Z. B.; Li, Y.; Wu, D.; Chen, W.; Li, Z. R.; Sun, C. C. The Interaction Between Superalkalis (M3O, M = Na, K) and a C20F20 Cage Forming Superalkali Electride Salt Molecules with Excess Electrons Inside the C20F20 Cage: Dramatic Superalkali Effect on the Nonlinear Optical Property. J. Mater. Chem. 2012, 22, 9652-9657. (28) Huang, S.; Liao, K.; Peng, B.; Luo, Q. On the Potential of Using the Al Superatom as an Excess Electron Acceptor To Construct Materials with Excellent Nonlinear Optical Properties. Inorg. Chem. 2016, 55, 4421-4427. (29) Reveles, J. U.; Baruah, T.; Zope, R. R. Al12Cu Superatom as Stable Building Block of Ionic Salts. J. Phys. Chem. C 2015, 119, 5129-5137. (30) Gutsev, G. L.; Boldyrev, A. I. DVM X-Alpha Calculations of the Electronic-structure of “Superalkali” Cations. Chem. Phys. Lett. 1982, 92, 262–266. (31) Gutsev, G. L.; Boldyrev, A. I. DVM X-α Calculations on the Ionization Potentials of Complex Anions and the Electron Affinities of MXk+1 “Superhalogens”. Chem. Phys. 1981, 56, 277-283. (32) Jissy, A. K.; Datta., A. Effect of External Electric Field on H-Bonding and π-stacking Interactions in Guanine Aggregates. Chemphyschem 2012, 13, 4163-4172. (33) Bhattacharyya, K.; Surendran, A.; Chowdhury, C.; Datta, A. Steric and Electric Field Driven Distortions in Aromatic Molecules: Spontaneous and Non-spontaneous Symmetry Breaking. Phys. Chem. Chem. Phys. 2016, 18, 31160-31167. (34) Darwish, N.; Aragones, A. C.; Darwish, T.; Ciampi, S.; Diez-Perez, I. Multi-responsive Photo- and Chemo-electrical Single-molecule Switches. Nano Lett. 2014, 14, 7064-7070. (35) Alemani, M.; Peters, M. V.; Hecht, S.; Rieder, K.-H.; Moresco, F.; Grill, L. Electric Field-Induced Isomerization of Azobenzene by STM. J. Am. Chem. Soc. 2006, 128, 14446-14447. (36) Oklejas, V.; Uibel, R. H.; Horton, R.; Harris, J. M. Electric-field Control of the Tautomerization and Metal Ion Binding Reactivity of 8-hydroxyquinoline Immobilized to an Electrode Surface. Anal. Chem. 2008, 80, 1891-1901. (37) Susarrey-Arce, A.; Tiggelaar, R. M.; Morassutto, M.; Geerlings, J.; Sanders, R. G. P.; Geerdink, B.; Schlautmann, S.; Lefferts, L.; van Houselt, A.; Gardeniers, J. G. E. A New ATR-IR Microreactor to Study Electric Field-driven Processes. Sensor Actuator B: Chem. 2015, 220, 13-21. (38) Zhou, Z.-J.; Li, X.-P.; Liu, Z.-B.; Li, Z.-R.; Huang, X.-R.; Sun, C.-C. Electric Field-Driven Acid−Base Chemistry: Proton Transfer from Acid (HCl) to Base (NH 3/H2O). J. Phys. Chem. A 2011, 115, 1418-1422. (39) Ceron-Carrasco, J. P.; Cerezo, J.; Jacquemin, D. How DNA Is Damaged by External Electric Fields: Selective Mutation vs. Random Degradation. Phys. Chem. Chem. Phys. 2014, 16, 8243-8246. (40) Arabi, A. A.; Matta, C. F. Effects of External Electric Fields on Double Proton Transfer 22

ACS Paragon Plus Environment

Page 22 of 25

Page 23 of 25 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

Kinetics in the Formic Acid Dimer. Phys. Chem. Chem. Phys. 2011, 13, 13738-13748. (41) Ceron-Carrasco, J. P.; Jacquemin, D. Electric Field Induced DNA Damage: An Open Door for Selective Mutations. Chem. Commun. 2013, 49, 7578-7580. (42) Ceron-Carrasco, J. P.; Jacquemin, D. Electric-Field Induced Mutation of DNA: A Theoretical Investigation of the GC Base Pair. Phys. Chem. Chem. Phys. 2013, 15, 4548-4553. (43) Biase, P. M. D.; Paggi, D. A.; Doctorovich, F.; Hildebrandt, P.; Estrin, D. A.; Murgida, D. H.; Mart, M. A. Molecular Basis for the Electric Field Modulation of Cytochrome c Structure and Function. J. Am. Chem. Soc. 2009, 131, 16248-16256. (44) Karafiloglou, P. Control of Delocalization and Structural Changes by Means of an Electric Field. J. Comput. Chem. 2006, 27, 1883-1891. (45) Harzmann, G. D.; Frisenda, R.; van der Zant, H. S. J.; Mayor, M. Single-Molecule Spin Switch Based on Voltage-Triggered Distortion of the Coordination Sphere. Angew. Chem. Int. Ed. 2015, 54, 13425-13430. (46) Son, Y. W.; Cohen, M. L.; Louie, S. G. Half-metallic Graphene Nanoribbons. Nature 2006, 444, 347-349. (47) Naaman, R.; Waldeck, D. H. Spintronics and Chirality: Spin Selectivity in Electron Transport Through Chiral Molecules. Annu. Rev. Phys. Chem. 2015, 66, 263-281. (48) Franzen, S.; Goldstein, R. F.; Boxer, S. G. Electric Field Modulation of Electron Transfer Reaction Rates in Isotropic Systems: Long Distance Charge Recombination in Photosynthetic Reaction Centers. J. Phys. Chem. 1990, 94, 5135-5149. (49) Lao, K.; Franzen, S.; Stanley, R. J.; Lambright, D. G.; Boxer, S. G. Effects of Applied Electric Fields on the Quantum Yields of the Initial Electron-transfer Steps in Bacterial Photosynthesis. 1. Quantum Yield Failure. J. Phys. Chem. 1993, 97, 13165-13171. (50) Wahadoszamen, M.; Nakabayashi, T.; Kang, S.; Imahori, H.; Ohta, N. External Electric Field Effects on Absorption and Fluorescence Spectra of a Fullerene Derivative and Its Mixture with Zinc-tetraphenylporphyrin Doped in a PMMA Film. J. Phys. Chem. B 2006, 110, 20354-20361. (51) Bhattacharyya, K.; Karmakar, S.; Datta, A. External Electric Field Control: Driving the Reactivity of Metal-free Azide-alkyne Click Reactions. Phys. Chem. Chem. Phys. 2017, 19, 22482-86. (52) Meir, R.; Chen, H.; Lai, W.; Shaik, S. Oriented Electric Fields Accelerate Diels-Alder Reactions and Control the Endo/exo Selectivity. Chemphyschem 2010, 11, 301-310. (53) Shaik, S.; de Visser, S. P.; Kumar, D. External Electric Field will Control the Selectivity of Enzymatic-Like Bond Activations. J. Am. Chem. Soc. 2004, 126, 11746-11749. (54) Hirao, H.; Chen, H.; Carvajal, M. A.; Wang, Y.; Shaik, S. Effect of External Electric Fields on the C−H Bond Activation Reactivity of Nonheme Iron−oxo Reagents. J. Am. Chem. Soc. 2008, 130, 3319-3327. (55) Lai, W.; Chen, H.; Cho, K.-B.; Shaik, S. External Electric Field Can Control the Catalytic Cycle of Cytochrome P450cam: A QM/MM Study. J. Phys. Chem. Lett. 2010, 1, 2082-2087. (56) Aragones, A. C.; Haworth, N. L.; Darwish, N.; Ciampi, S.; Bloomfield, N. J.; Wallace, G. G.; Diez-Perez, I.; Coote, M. L. Electrostatic Catalysis of a Diels-Alder Reaction. Nature 2016, 531, 88-91. (57) Gorin, C. F.; Beh, E. S.; Kanan, M. W. An Electric Field-induced Change in the Selectivity of a Metal Oxide-catalyzed Epoxide Rearrangement. J. Am. Chem. Soc. 2012, 134, 186-189. (58) Gorin, C. F.; Beh, E. S.; Bui, Q. M.; Dick, G. R.; Kanan, M. W. Interfacial Electric Field Effects on a Carbene Reaction Catalyzed by Rh Porphyrins. J. Am. Chem. Soc. 2013, 135, 11257-11265. 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

(59) Shaik, S.; Mandal, D.; Ramanan, R. Oriented Electric Fields as Future Smart Reagents in Chemistry. Nat. Chem. 2016, 8, 1091-1098. (60) Boys, S. F.; Bernardi, F. The Calculation of Small Molecular Interactions by the Differences of Separate Total Energies. Some Procedures with Reduced Errors. Mol. Phys. 1970, 19, 553-566. (61) Alkorta, I.; Elguero, J. Theoretical Study of Strong Hydrogen Bonds between Neutral Molecules:  The Case of Amine Oxides and Phosphine Oxides as Hydrogen Bond Acceptors. J. Phys. Chem. A 1999, 103, 272-279. (62) 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 A02; Gaussian, Inc.: Wallingford CT, 2009. (63) Dennington, R.; Keith, T.; Millam, J. GaussView, Version 5; Semichem Inc.: Shawnee Mission, KS, 2009. (64) Williams, E. R.; Faller, J. E.; Hill, H. A. New Experimental Test of Coulomb’s Law: A Laboratory Upper Limit on the Photon Rest Mass. Phys. Rev. Lett. 1971, 26, 721-724. (65) Oudar, J. L. Optical Nonlinearities of Conjugated Molecules. Stilbene Derivatives and Highly Polar Aromatic Compounds. J. Chem. Phys. 1977, 67, 446. (66) Oudar, J. L.; Chemla, D. S. Hyperpolarizabilities of the Nitroanilines and Their Relations to the Excited State Dipole Moment. J. Chem. Phys. 1977, 66, 2664. (67) Wang, X.-B.; Ding, C.-F.; Wang, L.-S. High Resolution Photoelectron Spectroscopy of C60−. J. Chem. Phys. 1999, 110, 8217-8220.

24

ACS Paragon Plus Environment

Page 24 of 25

Page 25 of 25 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

Table of Contents Graphic

25

ACS Paragon Plus Environment