Absorption and Mid-IR SHG in Two-Dimensional ... - ACS Publications

Nov 14, 2017 - ABSTRACT: Developing nonlinear optical nanodevices compatible for the current silicon chip industry is of great interest. Modified sili...
0 downloads 0 Views 1MB Size
Subscriber access provided by READING UNIV

Article

Absorption and Mid-IR SHG in Two-Dimensional Halogen and Hydrogen Saturated Silicene Series Guoyu Yang, and Kechen Wu J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.7b08810 • Publication Date (Web): 14 Nov 2017 Downloaded from http://pubs.acs.org on November 17, 2017

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 free 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 accessible to all readers and 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.

The Journal of Physical Chemistry C 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 32

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

Absorption and Mid-IR SHG in Two-Dimensional Halogen and Hydrogen Saturated Silicene Series Guoyu Yang†,‡ and Kechen Wu∗,† †

State Key Laboratory of Structural Chemistry, Fujian Institute of Research on the Structure of

Matter, Chinese Academy of Sciences, Fuzhou 350002, People’s Republic of China , and ‡

University of Chinese Academy of Sciences, Beijing 100049,People’s Republic of China

E-mail: [email protected]

*To whom correspondence should be addressed †Fujian Institute of Research on the Structure of Matter ‡University of Chinese Academy of Sciences

ACS Paragon Plus Environment

1

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 32

ABSTRACT

Developing nonlinear optical nano-devices compatible to current silicon chip industry is of great interests. Modified silicene, as its single atomic silicon thickness, has become a rising material. Here, based on first-principle calculations, two-dimentional Halogen and Hydrogen saturated silicene is proposed to be competing candidates for strong second harmonic gener ation (SHG) material working in mid-IR range. The SHG coefficient is more than ten times larger than KDP with high electron cloud deformability around halogen atoms. And the Si-H absorption peak is within the range of 2102 - 2140 cm-1 (4673 - 4757 nm), red shifted com pared to 3300 cm-1 (3030 nm) of C-H in graphane. In addition, the two-dimensional material construction method can be extended to other novel two-dimensional single-element materials including borophene, phosphorene, germanene, arsenene, antimonene and bismuthene.

Introduction The nonlinear opitcal (NLO) effect in IR region is crucial with wide application including bio and medical system1 , encrypted communication2 , high energy laser device3. However, corresponding high performance material has yet been well developed as failing to meet the five requirements at the same time, i.e. phase matching behavior, low absorption in working optical wavelength area, high laser induced damage thresholds(LIDTs), stable chemical and mechanic properties, large Second Harmonic Generation (SHG) coefficient4,5. Therefore, exploring new crystalline IR NLO material is of broad scientific and technological importance. If in crystalline NLO material6, the cations can be considered as frameworks mainly to stabilize the material, while the anions as building blocks to exhibit large electron cloud deformability and

ACS Paragon Plus Environment

2

Page 3 of 32

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

strong SHG effect5,7. Frequently discussed IR NLO material can be categorized into three groups based on their anions: chalcogenide (S group), including AgGaS2,AgGaSe28–12, VA group (P group), including ZnGeP28,13,14, halogen (X group), including Cs2Hg3I8 15–17,and HgBr218. Some recent work implies central symmetric material can exhibit NLO effect after special preparation into nano size19, even though, as a natural result of phase mismatching, the SHG coefficient di jk decreases as the crystal size going down into micrometer for most three dimensional (3D) NLO materials18,20,21. This defect, together with the requirement of equipment miniaturization, inspires scientists to search low dimension (LD) NLO materials. Among those LD materials investigated, two dimension (2D) single/few layer(s), such as MoS222–25, WSe226,27, h-BN22, and GaSe 28 , with their relatively low-cost top-down preparation technique, (i.e. exfoliation from the bulk) is more reported. However probably due to the limitation of laser generation and detection, most work is done in the visible and near-IR region. Silicene, with its atomic thickness 29–31 and compatibility with other electronics material32,33, was expected to be a promising nano material in silicon based chip industry34,35. It is one of the analogues of graphene36 has been synthesized and/or proposed, including borophene 37,38 , phosphorene 39 , germanene 40–42 , arsenene43,44, stanene 45–47 , antimonene 43 and bismuthene48.These single-element systems correspond to elements that are akin to carbon on periodic table of elements, and usually, as perfect 2D crystal, is central symmetric thus no strong SHG effect. To break central symmetry, methods like stacking 22,32 , charge accumulation 26 and electric field 27

, adding substrate 49,50 , strain and stress 33 has been adopted. Here, based on first principle

calculation, Si2XH (X refers to Halogen), a new 2D crystal, halogen and hydrogen saturated

ACS Paragon Plus Environment

3

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 32

silicene, or halogen substituted silicane51 , is proposed to be a promising 2D SHG material in IR region fulfilling the five requirements mentioned in paragraph one. The calculation results reveals that the silicon backbone plays the role of ’cation’ in the functional motif theory5,6, including stabilizing the material, red shifting the absorption away from mid-IR area, while the halogen atomsas ’anion’ to generate large SHG signal by offering highly deformable electron cloud.

Computational Methods Here, using Vienna ab initio simulation package (VASP)52 , density functional theory (DFT) based first principle calculation is performed. Projector augmented-wave (PAW) pseudopotential 53

is used at the general gradient approximation (GGA) level in the scheme of Perdew-Burke-

Ernzerhof (PBE) functional54. A vacuum layer of 50 Å along the z-direction is employed in order to avoid interactions between two Si2XH images in nearest-neighbor unit cell. The structure optimization is performed with 3x3x1 supercell, cutoff energy 550 eV, unitcell energy difference smaller than 10-4 eV. The atom position is relaxed with fixed lattice parameters until the force is less than 10-3 eV/Å and the lattice constant scanned to the accuracy of 0.01 Å. A Gammacentered Monkhorst-Pack scheme55 with k-point mesh of 7x7x1 is used. For static calculation, the convergence criterion of total energy is 10-6 eV, cutoff energy 750 eV. The k-point mesh is 7x7x1 for density of state (DOS) and optical properties. And for band calculation, 20 k points are intercepted between two symmetric points. The SHG coefficient was calculated using length gauge method56–58 without scissor operator correction.59 And the absorption spectrum is also confirmed by using many-body perturbation theory including quasiparticle GW0 and the Bethe-Salpeter equation (BSE).60,61 The total energy cutoff

ACS Paragon Plus Environment

4

Page 5 of 32

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

was set to 550 eV, and for the response function was set to 200 eV. A k-mesh of 5x5x1 with unit cell was used. A total of 160 bands and 96 frequency points were included in the GW0 calculation which was followed by the BSE spectrum calculation using the GW0 quasiparticle results. The 32 highest valence and 32 lowest conduction bands were included in the calculation of the excitonic states. The IR absorption wavenumber is calculated on crystal unit cell using the VASP Density Functional Perturbation Theory (DFPT), and intensity calculated by summing up the dipole derivative of rectified coordinate.62–65 The phonon dispersion spectrum is calculated using PHONOPY package66 with finite displacement method on supercell as or larger than 3x3x1.

Results and discussion Structure Monolayer Silicene, though well known as single layer (111) cut of fcc Si31, has rarely been prepared by liquid or mechanical exfoliation from the bulk like graphene36, antimonene67, MoS2 and h-BN22. Instead, reported synthesyzing methods is mainly bottom-up vacuum deposition31,68– 70

, like borophene37. To further saturize silicene, soft chemical synthesis hydrolization71,72,

followed by functional group substitution73–75 is a promising strategy. And substrate could be used to control which side to be covered or exposed for reaction109, 110. Though, possibly because that normal halogenated reagents (such as CCl4, Cl2 gas, PCl5) are highly toxic, to the best of our knowledge, Si2XH (Figure 1)has not been really synthesized by now. Therefore, first-principle based principle calculation should be serving as a very useful tool here.76

ACS Paragon Plus Environment

5

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 32

Figure 1: Structure of Si2XH (X = F, Cl, Br, I) Here the X;Y;Z directions were defined in this way: the Z direction was perpendicular to the Si2XH plane; X is the zigzag direction in the plane; Y is the armchair direction in the plane.

As seen from our phonon calculation result (See Supporting Information, Phonon Spectra), among the four Si2XH (X=F, Cl, Br, I) models, Si2FH and Si2IH are both stable with small imaginary frequency within the acceptable range compared with borophene77–79 and bismuthene80, while Si2ClH and Si2BrH has no imaginary frequency. The binding energy of each model (defined in (Eq. 1)) is respectively -0.29, -0.49, -0.42, -0.42, -0.58 eV/atom(Si2H2, Si2FH, Si2ClH, Si2BrH, Si2IH), which is consistent with previous calculation result in literature.81,82 (1)

ACS Paragon Plus Environment

6

Page 7 of 32

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

Where

and

crystal.

are respectively the total energy per atom of silicene and that of Si 2XH

is the total energy per atom in XH molecular.

and

denote the molar

fractions of the corresponding atoms.

For the structure, the Si2XH crystals belong to P3M1 symmetric group (group 156), and Si2H2 to P3-M1 symmetric group (group 164), with the Si sp3 hybrid orbital remains three fold symmetry (Figure 1).

Linear Optical Properties (Absorption) Since the three fold symmetry were kept in the Si2XH series, the relationship (Eq. 2,Eq. 3) in linear dielectric function hold. (2)

(3)

Also, the optical absorption in monolayer structure in different direction varies:83,84 in the direction perpendicular to the 2D surface is principally excitonic absorption, but the parallel com ponent not affected by that. The absorption coefficient is obtained by using Eq. 4:83,84

ACS Paragon Plus Environment

7

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 8 of 32

(4)

First, the electronic absorption calculated with DFT (Figure 2), (A.1) (A.2)) and GW0 (Figure 2), (B.1) (B.2)) was discussed. DFT calculation, though neglect many-body effect and exciton effect (which is usually illegible in bulk semiconductor material), was sometimes used to estimate the electronic optical absorption gap for being affordable.85 Here the DFT calculated absorption low energy limit is about 1 eV (1240 nm), which is without the mid-IR range (1500-6000nm).

Figure 2: Electronic and Vibrational Absorbtion of Si2H2 and Si2XH (A.1) (A.2): Electronic absorption coefficient calculated with DFT in the unit of eV and nm; (B.1) (B.2): Electronic absorption coefficient calculated with GW0 in the unit of eV and nm; (C.1) (C.2): Vibration Absorbtion in IR in the unit of cm-1 and nm.

ACS Paragon Plus Environment

8

Page 9 of 32

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

Here, the IR absorption was plotted both in wavenumber (cm-1) and wavelength (nm) with the Lorentz broadening of FWHM=100 cm-1 using Multiwfn.86

And for GW0, previous work usually stressed on two gaps in semiconductor: fundamental (or quasi-particle) and optical gaps, which were each calculated by ab initio GW0 to include electron-electron interaction, and by GW0 approach plus the Bethe-Salpeter equation (BSE) to include electron-hole interaction (excitonic effect).85 The GW calculated absorption result (Figure 2), (B.1) (B.2)) shares the same low-energy limit with DFT result (Figure 2), (A.1) (A.2)) (1 eV, 1240 nm), especially for Si2H2, the absorption low limit is consistent with previous work.87,88 However, the relative absorption strength of in-plane and perpendicular differs in these two methods. These result is in accordance with the conclusion proved in several groups that the perpendicular direction was principally affected by exciton, but the parallel was not83,84, so our calculation result should be considered as reliable in this aspect. Also, this validated our accuracy parameter choices in GW0 method. As the result, the electronic absorption in mid-IR range is eligible, so we turn attention to the vibration absorption in IR range. The IR region absorption is dominated by vibration. Comparison between the IR spectrum of Si2XH and Si2H2 (Figure 2), (C.1) (C.2)) gives two points: every model shares a small peak in the place about 2100 cm-1 (4761 nm), the peaks were assigned to Si-H bond vibration (see Supporting Information, IR Absorption); while in the broad and strong peak groups around 500 nm, every Si2XH has multiple/shoulder peaks, while Si2H2 does not. In Si2XH, part of these peaks corresponds to Si-X bond vibration, some were significantly coupled with Si-Si bond vibration, so the peak groups would be complex; while for Si2H2, for the structure is simpler with central symmetry, the peaks’ shape is plain. Concentrated on the absorption range,

ACS Paragon Plus Environment

9

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 10 of 32

elaboration on vibration (Figure 2), (C.1) (C.2)) reveals that on 1000-4000 nm, the Si2XH material could be considered as transparent, which takes a large part on mid-IR region (15006000 nm). So, the low absorption makes Si2XH series of strong candidates as SHG material in mid-IR.

Nonlinear Optical Properties Just like dielectric function, symmetry also simplifies the SHG coefficient analysis. In Si 2XH (P3M1, group 156) and Si2H2 (P3-M1, group 164), Si sp3 hybrid orbital remains three fold symmetry (Figure 1). So among the 10 reduced SHG coefficients di jk (i, j,k = X,Y,Z), with X,Y,Z in Figure 1 also denoting the direction of dielectric constant coordinate system,89 only four (showing in (Eq. 7) and (Eq. 8)) need to be paid attention to when calculating effective SHG coefficient de f f (Eq. 5) (Eq. 6).90 (5)

(6)

Where de f f (ooe) and de f f (eeo) denotes de f f in ooe and eeo mode (here o means ordinary and e means extraordinary ). And

and

are defined by optic polarization vector in polar

coordinates. group1 :

-dXXY = dYYY (7)

group2 :

dXXZ = dYYZ (8)

Nonlinear Optical Result. How to quantitatively evaluate the SHG effect of 2D monolayer (ML) material is another important issue to be considered. In previous literature, two methods were proposed: (i) to

ACS Paragon Plus Environment

10

Page 11 of 32

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

compare with a known freshly cleaved bulk surface (say, alpha-quartz22). This method is experimentally convenient, and the theory of bulk surface is well established91, yet it ignored the difference of thickness between the two models, ((sub)nano vs. (semi)infinite), thus hindered further analyzation. (ii) to compare with 3D bulk result92. This method adds a ’thickness’ (l) to define the ’volume’ of 2D ML material. Here l is decided by the layer distance of corresponding 3D crystals92, or microscope image of several-layer structure24. In this paper, SHG coefficient was defined within the first principle calculation scheme. To compare with the first method, the SHG coefficient direct calculation result in 3D vacuum layer d(vac) was multiplied with the vacuum layer height, (i.e. h = 50 Å here) (Eq. 9). The result d(2D) means SHG effect per unit area. (Figure 3, left Y axis). While to compare with the second method, as corresponding Si2XH 3D crystals do not actually exist, thickness l has to be newly defined. Here four two-layer (2L) models were built to include different stacking modes93. The layer distance of the four varied within 5%, so the average was used in (Eq. 10)(See Supporting Information, Geometry). The d(3D) here means the SHG coefficient per ’volume’. d(2D) = d(vac) ·h d(3D) = d(vac) ·h / l

(9) (10)

Here, d(2D) and d(3D) are ML SHG coefficients presented in 2D (per area) and 3D (per volume), and d(vac) is the result calculated with vacuum layer in models.

ACS Paragon Plus Environment

11

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 12 of 32

Figure 3: SHG Coefficient of Si2H2 and Si2XH Here, the black square, red dot, green up triangle, and blue down triangle is respectively for -dXXY , dYYY , dXXZ, dYYZ. And the solid characters are for d(2D) (left Y-axis), hollow for d(3D) (right Y-axis). The five horizontal lines, as the notation, denotes several large SHG 3D crystal in IR range, respectively: HgBrI, HgI2 21,94 in X group; AgGaS2, AgGaSe2 8–12 in S group; ZnGeP2 8,13,14 in P group.

The result of di jk is plotted in Figure 3. First, the SHG of Si2H2, a central symmetric structure, is very close to zero, which is consistent with the second order nonlinear optics principle.90 Second, it is shown that the two doubles -dXXY and dYYY , dXXZ and dYYZ data matched well, which is in accordance with previous symmetry analysis (Eq. 7) (Eq. 8). Thirdly, the corrected di jk (3D), is comparable with, if not stronger than, previous IR crystal material HgBrI,HgI2

21,94

(X group) ,

AgGaS2, AgGaSe2 8–12 (S group) ; and ZnGeP2 8,13,14 (P group) (Figure 3, right Y axis). Fourthly, the di jk increases as the X atom goes heavier, following the sequence of F, Cl, Br, I. This can be explained by two sides, heavier X makes the two sides of silicene more asymmetric, thus larger hyper-polarizabiliy90 (here hyper-polarizabiliy equals half of di jk 95); and halogen electron cloud d formability increases as the atoms go heavier96.

ACS Paragon Plus Environment

12

Page 13 of 32

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

Nonlinear Optical Band and DOS Analysis. To further specify the structure - property - function relationships, detailed band, density of state and partial charge analysis is performed. Since SHG means electron hopping from valence bands (VB) to conductive bands (CB), band analysis should be a good starting point. By adjusting the upper and lower boundary of the bands (Supporting Information, Define SHG of a Single Band) in NLO calculation, the contribution to SHG coefficient in certain direction from each band can be specified. For any di jk among the four discussed, it is found that the sum of the five most contributing VBs’ makes more than 80% of the total di jk, and that also applies to CB. Here Si2ClH is taken as an example; others are shown in (Supporting Information, Band and DOS). As shown in Figure 4, the five VBs and five CBs with most contribution of the each four effective di jk of Si2ClH were highlighted in the total bands. The ten bands are near Fermi level but not the closest ten, which means that the two-electron process selectivity is not just energy dependent. Careful investigation into the VB shows that the group -dXXY =dYYY share the same important five VBs and five CBs, while the dXXZ=dYYZ makes another group, which is in accordance with the P3M1 symmetry. And for the convenience of discussion, the bands were sorted by the energy at Gamma point and numbered. The criteria is that the top VB is numbered half the total number of the valence shell electrons (for example, 3x3 model Si2XH top VB corresponds to

× 3×

3×(2× 4 + 7 + 1) = 72). So here, as shown in Figure 4 Si2ClH, band No.71, 72 are important only for -dXXY and dYYY ; band No.49, 50 are for dXXZ and dYYZ; band No.52, 53 and 54 are important for both groups of di jk.

ACS Paragon Plus Environment

13

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 14 of 32

Figure 4: Electronic Properties of Si2ClH In the left figure for band, red lines are band No. 71,72,89,90, only important to -dXXY and dYYY; blue lines are band No. 49,50,76,76, only important to dXXZ and dYYZ; green lines are band No. 52,53,54,73,74,75, important to -dXXY and dYYY, dXXZ and dYYZ. This means for -dXXY and dYYY, the bands concerned are No. 52,53,54,71,72,73,74,75,89,90; for dXXZ and dYYZ, the bands concerned are 49, 50, 52, 53, 54, 73, 74, 75, 76, 77. All the bands are sorted and numbered by the energy of Gamma(G) points. The Valence Bands are Band No. 1 to No. 72.

Nonlinear Optical Charge Analysis. Partial charge density of the each seven bands was plotted in 3D real space (Supporting Information, Partial Charge Analysis). Band No.71 and 72 are basically the same situation, with charge distributed on X and Si atoms, but not bonded; No.49 and 50, 52, 53 and 54 are the similar, electron clouds is around X atoms extending towards Si (Figure 5). This comes to the

ACS Paragon Plus Environment

14

Page 15 of 32

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

conclusion that for -dXXY and dYYY , Cl atoms takes the main contribution, Si atoms are the auxiliary at best; and for dXXZ and dYYZ, the SHG signal is almost attributed to Cl atoms. This is in consistent with the previous analysis on NLO crystals binary and ternary metal halides.96

Figure 5: Partial Charge Analysis of Si2ClH Other Issues in Nonlinear Effect. For NLO crystals, proper transparency usually means low electron absorption and vibration absorption in working region. The Si2XH are semiconductor with band gap >1.8eV, calculated electron absorption is mainly in UV region (Supporting Information, Electronic Absorption) 97,98. And for vibration part in IR region(Supporting Information, IR Absorption), as the mid-IR region is 1500-6000 nm, the Si-H bond absorption peak ranging in 2102 - 2140 cm-1 (4673 4757 nm) indicates that Si2XH is better IR material compared to graphane, where C-H bond vibration is around 3300 cm-1 (3030 nm). The Si-X wavenumber decreases from 700 to 300 cm-1 (14286 nm to 33333 nm, or 14 mm to 33 mm) (Supporting Information, IR absorption). These Si-H and Si-X IR peaks are consistent with previous spectrum results99,100.

ACS Paragon Plus Environment

15

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 16 of 32

Si2XH, belonging to P3-M1 symmetry, is uniaxial crystal. But for monolayer, phase matching issue is trivial with atomical thickness perpendicular to the plane,24 so birefringence was illegible here. As to the laser damage threshold, many issues should be considered,

4

such as melting

temperature by dynamic simulation101, thermo-conductivity66,102, heat capacity 66 light absorption and heat conversion. These would be further complimented in our future work. To the best of our knowledge, at present, GW approach in low dimension material optical properties were mostly applied in linear optical aspect.85,103,104 In this paper, SHG coefficients were calculated by DFT method, which just include single-electron approximation in the twophoton process, and the actual multi-exciton process was indicated to be more prominent if correlation was taken into consideration. 105 So it is very likely our result could survive future reliable full GW calculation in many-electron process.106–108

Conclusion Si2XH, halogen and hydrogen saturated silicene, was predicted to be 2D monolayer mid-IR nonlinear optical material by first principle calculation. Electron properties analysis revealed that SHG signals mainly attributed to electron cloud of halogen atoms, and the transparency in midIR region attributed to the vibration red shifted by heavier silicon atoms. We expect our work not only inspire optoelectronic device design in silicon-based nano chip technology, but also be generalized to other graphene analogues.

ASSOCIATED CONTENT Supporting Information. Geometry parameters;

ACS Paragon Plus Environment

16

Page 17 of 32

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

Phonon Spectra; Define SHG of a Single Band; BAND and DOS; Partial Charge Analysis; IR Absorption.

AUTHOR INFORMATION Corresponding Author Kechen Wu, E-mail: [email protected] Author Contributions The manuscript was written with contributions from all authors. All authors have given approval to the final version of the manuscript. Funding Sources This study was supported by the National Natural Science Foundation of China (No. 21673240), Foreign Cooperation Project of Fujian Province (No. 2017I0019), and the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase). Notes The authors declare no competing financial interest.

ACKNOWLEDGMENT

ACS Paragon Plus Environment

17

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 18 of 32

This study was supported by the National Natural Science Foundation of China (No. 21673240), Foreign Cooperation Project of Fujian Province (No. 2017I0019), and the Special Program for Applied Research on Super Computation of the NSFC-Guangdong Joint Fund (the second phase).

References (1) Jung, Y.; Tong, L.; Tanaudommongkon, A.; Cheng, J.X.; Yang, C. In Vitro and in Vivo Nonlinear Optical Imaging of Silicon Nanowires. Nano Letters 2009, 9, 2440–2444. (2) Walter, F.; Li, G.; Meier, C.; Zhang, S.; Zentgraf, T. Ultrathin Nonlinear Metasurface for Optical Image Encoding. Nano Lett 2017, 17, 3171–3175. (3) Dergachev, A. In Solid State Lasers Xxii: Technology and Devices; Clarkson, W. A., Shori, R. K., Eds.; Proceedings of SPIE; 2013; Vol. 8599; Article Number 85990B. (4) Lin, Z.; Jiang, X.; Kang, L.; Gong, P.; Luo, S.; Lee, M.-H. First-Principles Materials Applications and Design of Nonlinear Optical Crystals. Journal of Physics D: Applied Physics 2014, 47, 253001. (5) Liu, B.W.; Zeng, H.Y.; Jiang, X.M.; Wang, G.-E.; Li, S.F.; Xu, L.; Guo, G.-C. [A3X][Ga3PS8] (A = K, Rb; X = Cl, Br): Promising IR Non-Linear Optical Materials Exhibiting Concurrently Strong Second-Harmonic Generation and High Laser Induced Damage Thresholds. Chem. Sci. 2016, 7, 6273–6277. (6) Guo, G.; Yao, Y.; Wu, K.; Wu, L.; Huang, J. Studies on the Structure-Sensitive Functional Materials. Progress In Chemistry 2001, 13, 151.

ACS Paragon Plus Environment

18

Page 19 of 32

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

(7) Jiang, X. M.; Zhang, M. J.; Zeng, H. Y.; Guo, G. C.; Huang, J. S. Inorganic Supermolecular Compounds with 3-D Chiral Frameworks Show Potential as Both Mid-IR Second-Order Nonlinear Optical and Piezoelectric Materials. J Am Chem Soc 2011, 133, 3410–3418. (8) Roberts, D. A. Simplified Characterization of Uniaxial and Biaxial Nonlinear Optical Crystals - A Plea for Standardization of Nomenclature and Conventions. IEEE Journal of Quantum Electronics 1992, 28, 2057–2074. (9) Kupecek, P. J.; Schwartz, C. A.; Chemla, D. S. Silver Thiogallate (AgGaS2) .1. Nonlinear Optical-Properties. IEEE Journal of Quantum Electronics 1974, QE10, 540–545. (10) Bhar, G. C.; Smith, R. C. Silver Thiogallate (AgGaS2) .2. Linear Optical-Properties. IEEE Journal of Quantum Electronics 1974, QE10, 546–550. (11) Boyd, G. D.; Kasper, H.; McFee, J. H. Linear and Nonlinear Optical Properties of AgGaS2, CuGaS2, and CuInS2, and Theory of Wedge Technique for Mearsurement of Nonlinear Coefficients. IEEE Journal of Quantum Electronics 1971, QE 7, 563–573. (12) Russell, D. A.; Ebert, R. Efficient Generation and Heterodyne-Detection of 4.75-mu-m Light with 2nd-Harmonic Generation. Applied Optics 1993, 32, 6638–6644. (13) Boyd, G. D.; Buehler, E.; Storz, F. G. Linear and Nonlinear Optical Properties of ZnGeP2 and CdSe. Applied Physics Letters 1971, 18, 301–304. (14) Mason, P. D.; Jackson, D. J.; Gorton, E. K. CO2-Laser Frequency-Doubling in ZnGeP2. Optics Communications 1994, 110, 163–166.

ACS Paragon Plus Environment

19

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 20 of 32

(15) Zhang, G.; Qin, J.; Liu, T.; Zhu, T.; Fu, P.; Wu, Y.; Chen, C. Synthesis, Characterization, and Crystal Growth of Cs2Hg3I8: A New Second-Order Nonlinear Optical Material. Crystal Growth & Design 2008, 8, 2946–2949. (16) Wu, Q.; Huang, Y.; Meng, X.; Zhong, C.; Chen, X.; Qin, J. Exploration of New Second Order Nonlinear Optical Materials of The Cs-Hg-Br-I System. Dalton Trans 2014, 43, 8899–904. (17) Sathiskumar, S.; Kathiravan, P.; Balakrishnan, T. In Proceedings of the 59th Dae Solid State Physics Symposium 2014, AIP Conference Proceedings; Bhattacharyya, D., Chitra, R., Sahoo, N. K., Eds.;2015; Vol. 1665; Article Number: 100014. (18) Liu, T.; Qin, J.; Zhang, G.; Zhu, T.; Niu, F.; Wu, Y.; Chen, C. Mercury Bromide (HgBr2): A Promising Nonlinear Optical Material In IR Region with a High Laser Damage Threshold. Applied Physics Letters 2008, 93, 091102. (19) Makarov, S. V.; Petrov, M. I.; Zywietz, U.; Milichko, V.; Zuev, D.; Lopanitsyna, N.; Kuksin, A.; Mukhin, I.; Zograf, G.; Ubyivovk, E. et al. Efficient Second-Harmonic Generation in Nanocrystalline Silicon Nanoparticles. Nano Lett 2017, 17, 3047–3053. (20) Dang, Y.; Meng, X.; Jiang, K.; Zhong, C.; Chen, X.; Qin, J. A Promising Nonlinear Optical Material in the Mid-IR Region: New Results on Synthesis, Crystal Structure and Properties of Noncentrosymmetric beta-HgBrCl. Dalton Trans 2013, 42, 9893–9897. (21) Wu, Q.; Li, Y.; Chen, H.; Jiang, K.; Li, H.; Zhong, C.; Chen, X.; Qin, J. HgBrI: A Promising Nonlinear Optical Material in IR Region. Inorganic Chemistry Communications 2013, 34, 1–3.

ACS Paragon Plus Environment

20

Page 21 of 32

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

(22) Li, Y.; Rao, Y.; Mak, K. F.; You, Y.; Wang, S.; Dean, C. R.; Heinz, T. F. Probing Symmetry Properties of Few-Layer MoS2 and h-BN by Optical Second-Harmonic Generation. Nano Lett 2013, 13, 3329–33. (23) Hsu, W.T.; Zhao, Z.-A.; Li, L.J.; Chen, C.-H.; Chiu, M.H.; Chang, P.S.; Chou, Y.C.; Chang, W.-H. Second Harmonic Generation from Artificially Stacked Transition Metal Dichalcogenide Twisted Bilayers. ACS Nano 2014, 8, 2951–2958. (24) Kumar, N.; Najmaei, S.; Cui, Q.; Ceballos, F.; Ajayan, P. M.; Lou, J.; Zhao, H. Second Harmonic Microscopy of Monolayer MoS2. Physical Review B 2013, 87. (25) Malard, L. M.; Alencar, T. V.; Barboza, A. P. M.; Mak, K. F.; de Paula, A. M. Observation of Intense Second Harmonic Generation from MoS2 Atomic Crystals. Physical Review B 2013, 87. (26) Seyler, K. L.; Schaibley, J. R.; Gong, P.; Rivera, P.; Jones, A. M.; Wu, S.; Yan, J.; Mandrus, D. G.; Yao, W.; Xu, X. Electrical Control of Second-Harmonic Generation in a WSe2 Monolayer Transistor. Nature Nanotechnology 2015, 10, 407–411. (27) Yu, H.; Talukdar, D.; Xu, W.; Khurgin, J. B.; Xiong, Q. Charge-Induced Second-Harmonic Generation in Bilayer WSe2. Nano Lett. 2015, 15, 5653–5657. (28) Zhou, X.; Cheng, J.; Zhou, Y.; Cao, T.; Hong, H.; Liao, Z.; Wu, S.; Peng, H.; Liu, K.; Yu, D. Strong Second-Harmonic Generation in Atomic Layered GaSe. J. Am. Chem. Soc. 2015, 137, 7994–7997.

ACS Paragon Plus Environment

21

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 22 of 32

(29) Cahangirov, S.; Topsakal, M.; Akturk, E.; Sahin, H.; Ciraci, S. Two- and One-Dimensional Honeycomb Structures of Silicon and Germanium. Phys Rev Lett 2009, 102, 236804. (30) Sahin, H.; Cahangirov, S.; Topsakal, M.; Bekaroglu, E.; Akturk, E.; Senger, R. T.; Ciraci, S. Monolayer Honeycomb Structures of Group-IV Elements and III-V Binary Compounds: FirstPrinciples Calculations. Physical Review B 2009, 80. (31) Vogt, P.; De Padova, P.; Quaresima, C.; Avila, J.; Frantzeskakis, E.; Asensio, M. C.; Resta, A.; Ealet, B.; Le Lay, G. Silicene: Compelling Experimental Evidence for Graphene like TwoDimensional Silicon. Phys Rev Lett 2012, 108, 155501. (32) Hu, W.; Li, Z.; Yang, J. Structural, Electronic, and Optical Properties of Hybrid Silicene and Graphene Nanocomposite. J Chem Phys 2013, 139, 154704. (33) Schriever, C.; Bianco, F.; Cazzanelli, M.; Ghulinyan, M.; Eisenschmidt, C.; de Boor, J.; Schmid, A.; Heitmann, J.; Pavesi, L.; Schilling, J. Second-Order Optical Nonlinearity in Silicon Waveguides: Inhomogeneous Stress and Interfaces. Advanced Optical Materials 2015, 3, 129– 136. (34) Schriever, C.; Bianco, F.; Cazzanelli, M.; Ghulinyan, M.; Eisenschmidt, C.; de Boor, J.; Schmid, A.; Heitmann, J.; Pavesi, L.; Schilling, J. Second-Order Optical Nonlinearity in Silicon Waveguides: Inhomogeneous Stress and Interfaces. Advanced Optical Materials 2015, 3, 129– 136. (35) Tao, L.; Cinquanta, E.; Chiappe, D.; Grazianetti, C.; Fanciulli, M.; Dubey, M.; Molle, A.; Akinwande, D. Silicene Field-Effect Transistors Operating at Room Temperature. Nature Nanotechnology 2015, 10, 227–31.

ACS Paragon Plus Environment

22

Page 23 of 32

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

(36) Novoselov, K. S.; Geim, A. K.; Morozov, S. V.; Jiang, D.; Zhang, Y.; Dubonos, S. V.; Grigorieva, I. V.; Firsov, A. A. Electric Field Effect in Atomically Thin Carbon Films. Science 2004, 306, 666–669. (37) Mannix, A. J.; Zhou, X.-F.; Kiraly, B.; Wood, J. D.; Alducin, D.; Myers, B. D.; Liu, X.; Fisher, B. L.; Santiago, U.; Guest, J. R.; Yacaman, M. J.; Ponce, A.; Oganov, A. R.; Hersam, M. C.; Guisinger, N. P. Synthesis of Borophenes: Anisotropic, Two-Dimensional Boron Polymorphs. Science 2015, 350, 1513–1516. (38) Wu, X.; Dai, J.; Zhao, Y.; Zhuo, Z.; Yang, J.; Zeng, X. C. Two-Dimensional Boron Monolayer Sheets. ACS Nano 2012, 6, 7443–7453. (39) Liu, H.; Neal, A. T.; Zhu, Z.; Luo, Z.; Xu, X.; Tomanek, D.; Ye, P. D. Phosphorene: An Unexplored 2D Semiconductor with a High Hole Mobility. ACS Nano 2014, 8, 4033–4041. (40) Ni, Z.; Liu, Q.; Tang, K.; Zheng, J.; Zhou, J.; Qin, R.; Gao, Z.; Yu, D.; Lu, J. Tunable Bandgap in Silicene and Germanene. Nano Letters 2012, 12, 113–118. (41) Bianco, E.; Butler, S.; Jiang, S.; Restrepo, O. D.; Windl, W.; Goldberger, J. E. Stability and Exfoliation of Germanane: A Germanium Graphane Analogue. ACS Nano 2013, 7, 4414–4421. (42) Davila, M. E.; Xian, L.; Cahangirov, S.; Rubio, A.; Le Lay, G. Germanene: a Novel Two Dimensional Germanium Allotrope Akin to Graphene and Silicene. New Journal of Physics 2014, 16.

ACS Paragon Plus Environment

23

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 24 of 32

(43) Zhang, S.; Yan, Z.; Li, Y.; Chen, Z.; Zeng, H. Atomically Thin Arsenene and Antimonene: Semimetal-Semiconductor and Indirect-Direct Band-Gap Transitions. Angewandte Chemie International Edition 2015, 54, 3112–3115. (44) Kamal, C.; Ezawa, M. Arsenene: Two-Dimensional Buckled and Puckered Honeycomb Arsenic Systems. Physical Review B 2015, 91. (45) Balendhran, S.; Walia, S.; Nili, H.; Sriram, S.; Bhaskaran, M. Elemental Analogues of Graphene: Silicene, Germanene, Stanene, and Phosphorene. Small 2015, 11, 640–652. (46) Tang, P.; Chen, P.; Cao, W.; Huang, H.; Cahangirov, S.; Xian, L.; Xu, Y.; Zhang, S.-C.; Duan, W.; Rubio, A. Stable Two-Dimensional Dumbbell Stanene: A Quantum Spin Hall Insulator. Physical Review B 2014, 90. (47) Zhu, F.F.; Chen, W.J.; Xu, Y.; Gao, C.L.; Guan, D.D.; Liu, C.H.; Qian, D.; Zhang, S.C.; Jia, J.F. Epitaxial Growth of Two-Dimensional Stanene. Nature Materials 2015, 14, 1020–1025. (48) Zhang, S.; Xie, M.; Li, F.; Yan, Z.; Li, Y.; Kan, E.; Liu, W.; Chen, Z.; Zeng, H. Semiconducting Group 15 Monolayers: A Broad Range of Band Gaps and High Carrier Mobilities. Angewandte Chemie-International Edition 2016, 55, 1666–1669. (49) Feng, Y. et al. Direct Evidence of Interaction-Induced Dirac Cones in a Monolayer Silicene/Ag(111) System. Proceedings of the National Academy of Sciences of the United States of America 2016, 113, 14656–14661. (50) Dean, J. J.; van Driel, H. M. Graphene and Few-Layer Graphite Probed by Second Harmonic Generation: Theory and Experiment. Physical Review B 2010, 82.

ACS Paragon Plus Environment

24

Page 25 of 32

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

(51) Pulci, O.; Gori, P.; Marsili, M.; Garbuio, V.; Del Sole, R.; Bechstedt, F. Strong Excitons in Novel Two-Dimensional Crystals: Silicane and Germanane. EPL 2012, 98. (52) Kresse, G.; Hafner, J. Ab Initio Molecular Dynamics for Open-Shell Transition Metals. Physical Review B 1993, 48, 13115–13118. (53) Blochl, P. E. Projector Augmented-Wave Method. Physical Review B 1994 50, 17953– 17979. (54) Perdew, J. P.; Burke, K.; Ernzerhof, M. Generalized Gradient Approximation Made Simple. Physical Review Letters 1996, 77, 3865–3868. (55) Monkhorst, H. J.; Pack, J. D. Special Points for Brillouin-Zone Integrations. Physical Review B 1976, 13, 5188–5192. (56) Aversa, C.; Sipe, J. E. Nonlinear Optical Susceptibilities of Semiconductors: Results with a Length-Gauge Analysis. Physical Review B 1995, 52, 14636–14645. (57) Rashkeev, S. N.; Lambrecht, W. R. L.; Segall, B. Efficient Ab Initio Method for the Calculation of Frequency-Dependent Second-Order Optical Response in Semiconductors. Phys. ev. B 1998, 57, 3905–3919. (58) Fang, Z.; Lin, J.; Liu, R.; Liu, P.; Li, Y.; Huang, X.; Ding, K.; Ning, L.; Zhang, Y. Computational Design of Inorganic Nonlinear Optical Crystals Based on a Genetic Algorithm. CrystEngComm 2014, 16, 10569–10580. (59) Levine, Z. H.; Allan, D. C. Calculation of the Nonlinear Susceptibility for Optical Second Harmonic Generation in III-V Semiconductors. Phys Rev Lett 1991, 66, 41–44.

ACS Paragon Plus Environment

25

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 26 of 32

(60) Rödl, C.; Bechstedt, F. Optical and Energy-Loss Spectra of the Antiferromagnetic Transition Metal Oxides MnO, FeO, CoO, and NiO Including Quasiparticle and Excitonic Effects. Phys. Rev. B 2012, 86, 235122. (61) Rinke, P.; Schleife, A.; Kioupakis, E.; Janotti, A.; Rödl, C.; Bechstedt, F.; Scheffler, M.; Van de Walle, C. G. First-Principles Optical Spectra for F Centers in MgO. Phys. Rev. Lett. 2012, 108, 126404. (62) Karhanek, D.; Bucko, T.; Hafner, J. A Density-Functional Study of the Adsorption of Methane-Thiol on the (111) Surfaces of the Ni-Group Metals: II. Vibrational Spectroscopy. J Phys Condens Matter 2010, 22, 265006. (63) Giannozzi, P.; Baroni, S. Vibrational and Dielectric Properties of C60 from Densityl’ Perturbation Theory. The Journal of Chemical Physics 1994, 100, 8537–8539. (64) Esfarjani, K.; Hashi, Y.; Onoe, J.; Takeuchi, K.; Kawazoe, Y. Vibrational Modes and IR Analysis of Neutral Photopolymerized C-60 Dimers. Physical Review B 1998, 57, 223–229. (65) Baroni, S.; de Gironcoli, S.; Dal Corso, A.; Giannozzi, P. Phonons and Related Crystal Properties from Density-Functional Perturbation Theory. Reviews of Modern Physics 2001, 73, 515–562. (66) Togo, A.; Tanaka, I. First Principles Phonon Calculations in Materials Science. Scripta Materialia 2015, 108, 1–5.

ACS Paragon Plus Environment

26

Page 27 of 32

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

(67) Huo, C.; Sun, X.; Yan, Z.; Song, X.; Zhang, S.; Xie, Z.; Liu, J.; Ji, J.; Jiang, L.; Zhou, S.; Zeng, H. Few-Layer Antimonene: Large Yield Synthesis, Exact Atomical Structure, and Outstanding Optical Limiting. Journal of the American Chemical Society 2017, 139, 3568– 3568. (68) Aufray, B.; Kara, A.; Vizzini, S.; Oughaddou, H.; L´leandri, C.; Ealet, B.; Lay, G. L. Graphene-Like Silicon Nanoribbons on Ag(110): A Possible Formation of Silicene. Applied Physics Letters 2010, 96, 183102. (69) Fleurence, A.; Friedlein, R.; Ozaki, T.; Kawai, H.; Wang, Y.; Yamada-Takamura, Y. Experimental Evidence for Epitaxial Silicene on Diboride Thin Films. Phys. Rev. Lett. 2012, 108, 245501. (70) Meng, L.; Wang, Y.; Zhang, L.; Du, S.; Wu, R.; Li, L.; Zhang, Y.; Li, G.; Zhou, H.; Hofer, W. A.; Gao, H.J. Buckled Silicene Formation on Ir(111). Nano Letters 2013, 13, 685–690. (71) Dahn, J. R.; Way, B. M.; Fuller, E.; Tse, J. S. Structure of Siloxene and Layered Polysilane (Si6H6). Phys. Rev. B 1993, 48, 17872–17877. (72) Yamanaka, S.; Matsuura, H.; Ishikawa, M. New Deintercalation Reaction of Calcium from Calcium Disilicide. Synthesis of Layered Polysilane. Materials Research Bulletin 1996, 31, 307– 316. (73) Sugiyama, Y.; Okamoto, H.; Mitsuoka, T.; Morikawa, T.; Nakanishi, K.; Ohta, T.; Nakano, H. Synthesis and Optical Properties of Monolayer Organosilicon Nanosheets. Journal of the American Chemical Society 2010, 132, 5946–5947.

ACS Paragon Plus Environment

27

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 28 of 32

(74) Buriak, J. M. Organometallic Chemistry on Silicon and Germanium Surfaces. Chemical Reviews 2002, 102, 1271–1308. (75) Bansal, A.; Li, X.; Yi, S. I.; Weinberg, W. H.; Lewis, N. S. Spectroscopic Studies of the Modification of Crystalline Si(111) Surfaces with Covalently-Attached Alkyl Chains Using a Chlorination/Alkylation Method. The Journal of Physical Chemistry B 2001, 105, 10266–10277. (76) Wu, C., Inorganic Two-dimensional Nanomaterials; Smart Materials Series; The Royal Society of Chemistry, 2017; (77) Chen, X.; Liu, J.; Zhang, W.; Xiao, B.; Zhang, P.; Li, L. First-Principles Study on the Mechanism of Hydrogen Decomposition and Spillover on Borophene. The Journal of Physical Chemistry C 2017, 121, 17314–17320. (78) Xu, L.C.; Du, A.; Kou, L. Hydrogenated Borophene as a Stable Two-Dimensional Dirac Material with an Ultrahigh Fermi Velocity. Phys. Chem. Chem. Phys. 2016, 18, 27284–27289. (79) Kou, L.; Ma, Y.; Tang, C.; Sun, Z.; Du, A.; Chen, C. Auxetic and Ferroelastic Borophane: A Novel 2D Material with Negative Possion Ratio and Switchable Dirac Transport Channels. Nano Letters 2016, 16, 7910–7914, (80) Akturk, E.; Akturk, O. U.; Ciraci, S. Single and Bilayer Bismuthene: Stability at High Temperature and Mechanical and Electronic Properties. Physical Review B 2016, 94. (81) Lu, N.; Li, Z.; Yang, J. Electronic Structure Engineering via On-Plane Chemical Functionalization: A Comparison Study on Two-Dimensional Polysilane and Graphane. The Journal of Physical Chemistry C 2009, 113, 16741–16746.

ACS Paragon Plus Environment

28

Page 29 of 32

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

(82) Gao, N.; Zheng, W. T.; Jiang, Q. Density Functional Theory Calculations for Two Dimensional Silicene with Halogen Functionalization. Phys Chem Chem Phys 2012, 14, 257–61. (83) Wang, R.; Pi, X.; Ni, Z.; Liu, Y.; Yang, D. Density Functional Theory Study on Organically Surface-Modified Silicene. RSC Adv. 2015, 5, 33831–33837. (84) Matthes, L.; Gori, P.; Pulci, O.; Bechstedt, F. Universal Infrared Absorbance of Two Dimensional Honeycomb Group-IV Crystals. Physical Review B 2013, 87. (85) Wang, Y.; Huang, P.; Ye, M.; Quhe, R.; Pan, Y.; Zhang, H.; Zhong, H.; Shi, J.; Lu, J. Many-Body Effect, Carrier Mobility, and Device Performance of Hexagonal Arsenene and Antimonene. Chemistry of Materials 2017, 29, 2191–2201. (86) Lu, T.; Chen, F. Multiwfn: A Multifunctional Wavefunction Analyzer. Journal of Computational Chemistry 2012, 33, 580. (87) Wei, W.; Dai, Y.; Huang, B.; Jacob, T. Many-Body Effects in Silicene, Silicane, Germanene and Germanane. Phys. Chem. Chem. Phys. 2013, 15, 8789–8794. (88) Pulci, O.; Gori, P.; Marsili, M.; Garbuio, V.; Sole, R. D.; Bechstedt, F. Strong Excitons in Novel Two-Dimensional Crystals: Silicane and Germanane. EPL (Europhysics Letters) 2012, 98, 37004. (89) Kleinman, D. A. Nonlinear Dielectric Polarization in Optical Media. Physical Review 1962, 126, 1977–1979. (90) Valentin, G. D. P. D.; Sci., G. G. G. P. D.; Lotsch, D. I. H. K. V. Handbook of Nonlinear Optical Crystals; Springer„ 1997; pp 3225–3231.

ACS Paragon Plus Environment

29

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 30 of 32

(91) Bloembergen, N.; Pershan, P. S. Light Waves at the Boundary of Nonlinear Media. Physical Review 1962, 128, 606–622. (92) Wang, C.Y.; Guo, G.Y. Nonlinear Optical Properties of Transition-Metal Dichalcogenide MX2(M = Mo, W; X = S, Se) Monolayers and Trilayers from First-Principles Calculations. The Journal of Physical Chemistry C 2015, 119, 13268–13276. (93) Li, Y.; Li, F.; Chen, Z. Graphane/Fluorographene Bilayer: Considerable C-H...F-C Hydrogen Bonding and Effective Band Structure Engineering. J. Am. Chem. Soc. 2012, 134, 11269–75. (94) Shi, X.; Ma, Z.; He, C.; Wu, K. Strong SHG Responses Predicted in Binary Metal Halide Crystal HgI2. Chemical Physics Letters 2014, 608, 219–223. (95) Shen, Y. R. The Principles of Nonlinear Optics; Wiley„ 1984. (96) Wu, X. T.; Chen, L. Structure-Property Relationships in Non-Linear Optical Crystals II; Springer Berlin Heidelberg, 2012; pp 10171–10178. (97) Li, J.; Duan, C. G.; Gu, Z. Q.; Wang, D.S. Linear Optical Properties and Multiphoton Absorption of Alkali Halides Calculated from First Principles. Phys. Rev. B 1998, 57, 2222–2228. (98) Guo, G. Y.; Chu, K. C.; Wang, D.-S.; Duan, C.-G. Linear and Nonlinear Optical Properties of Carbon Nanotubes from First-Principles Calculations. Phys. Rev. B 2004, 69, 205416. (99) Bürger, H.; Kilian, W. Spektroskopische Untersuchungen an Tris(Tri-Methylsilyl)Silan und Silan-D1. Journal of Organometallic Chemistry 1969, 18, 299 – 306.

ACS Paragon Plus Environment

30

Page 31 of 32

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

(100) Bürger, H.; Kilian, W.; Burczyk, K. Darstellung, Eigenschaften, Kmr-Spektren, Schwingungs-Spektren und Normal Koordinaten Analyse von Tris(Tri-Methylsilyl) Halogen Silanen. Journal of Organometallic Chemistry 1970, 21, 291 – 301. (101) Bucher, D.; Pierce, L. C. T.; McCammon, J. A.; Markwick, P. R. L. On the Use of Accelerated Molecular Dynamics to Enhance Configurational Sampling in Ab Initio Simulations. Journal of Chemical Theory and Computation 2011, 7, 890–897. (102) Feng, J.; Wan, C.; Xiao, B.; Zhou, R.; Pan, W.; Clarke, D. R. Calculation of the Thermal Conductivity of L2SrAl2O7 (L = La, Nd, Sm, Eu, Gd, Dy). Phys. Rev. B 2011, 84, 024302. (103) Shahrokhi, M. Quasi-Particle Energies and Optical Excitations of Novel Porous Graphene Phases from First-Principles Many-Body Calculations. Diamond and Related Materials 2017, 77, 35–40. (104) Lu, P.; Wu, L.; Yang, C.; Liang, D.; Quhe, R.; Guan, P.; Wang, S. Quasiparticle and Optical Properties of Strained Stanene and Stanane. Sci. Rep. 2017, 7, 3912. (105) Baer, R.; Rabani, E. Expeditious Stochastic Calculation of Multiexciton Generation Rates in Semiconductor Nanocrystals. Nano Lett 2012, 12, 2123–2128. (106) Hybertsen, M. S.; Louie, S. G. Electron Correlation in Semiconductors and Insulators: Band Gaps and Quasiparticle Energies. Phys. Rev. B 1986, 34, 5390–5413. (107) Rohlfing, M.; Louie, S. G. Electron-Hole Excitations and Optical Spectra from First Principles. Phys. Rev. B 2000, 62, 4927–4944.

ACS Paragon Plus Environment

31

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 32 of 32

(108) Deslippe, J.; Samsonidze, G.; Strubbe, D. A.; Jain, M.; Cohen, M. L.; Louie, S. G. BerkeleyGW: A Massively Parallel Computer Package for the Calculation of the Quasiparticle and Optical Properties of Materials and Nanostructures. Computer Physics Communications 2012, 183, 1269–1289. (109) Qiu, J.; Fu, H.; Xu, Y.; Oreshkin, A. I.; Shao, T.; Li, H.; Meng, S.; Chen, L.; Wu, K.; Ordered and Reversible Hydrogenation of Silicene Phys. Rev. Lett. 2015, 114, 126101 (110) W. Wang, W. Olovsson, and R. I. G. Uhrberg; Band Structure of Hhydrogenated Silicene on Ag(111): Evidence for Half-Silicane , Phys. Rev. B, 2016, 93, 081406

TOC Graphic

ACS Paragon Plus Environment

32