\documentclass[11pt,twoside]{article}\makeatletter

\IfFileExists{xcolor.sty}%
  {\RequirePackage{xcolor}}%
  {\RequirePackage{color}}
\usepackage{colortbl}
\usepackage{wrapfig}
\usepackage{ifxetex}
\ifxetex
  \usepackage{fontspec}
  \usepackage{xunicode}
  \catcode`⃥=\active \def⃥{\textbackslash}
  \catcode`❴=\active \def❴{\{}
  \catcode`❵=\active \def❵{\}}
  \def\textJapanese{\fontspec{Noto Sans CJK JP}}
  \def\textChinese{\fontspec{Noto Sans CJK SC}}
  \def\textKorean{\fontspec{Noto Sans CJK KR}}
  \setmonofont{DejaVu Sans Mono}
  
\else
  \IfFileExists{utf8x.def}%
   {\usepackage[utf8x]{inputenc}
      \PrerenderUnicode{–}
    }%
   {\usepackage[utf8]{inputenc}}
  \usepackage[english]{babel}
  \usepackage[T1]{fontenc}
  \usepackage{float}
  \usepackage[]{ucs}
  \uc@dclc{8421}{default}{\textbackslash }
  \uc@dclc{10100}{default}{\{}
  \uc@dclc{10101}{default}{\}}
  \uc@dclc{8491}{default}{\AA{}}
  \uc@dclc{8239}{default}{\,}
  \uc@dclc{20154}{default}{ }
  \uc@dclc{10148}{default}{>}
  \def\textschwa{\rotatebox{-90}{e}}
  \def\textJapanese{}
  \def\textChinese{}
  \IfFileExists{tipa.sty}{\usepackage{tipa}}{}
\fi
\def\exampleFont{\ttfamily\small}
\DeclareTextSymbol{\textpi}{OML}{25}
\usepackage{relsize}
\RequirePackage{array}
\def\@testpach{\@chclass
 \ifnum \@lastchclass=6 \@ne \@chnum \@ne \else
  \ifnum \@lastchclass=7 5 \else
   \ifnum \@lastchclass=8 \tw@ \else
    \ifnum \@lastchclass=9 \thr@@
   \else \z@
   \ifnum \@lastchclass = 10 \else
   \edef\@nextchar{\expandafter\string\@nextchar}%
   \@chnum
   \if \@nextchar c\z@ \else
    \if \@nextchar l\@ne \else
     \if \@nextchar r\tw@ \else
   \z@ \@chclass
   \if\@nextchar |\@ne \else
    \if \@nextchar !6 \else
     \if \@nextchar @7 \else
      \if \@nextchar (8 \else
       \if \@nextchar )9 \else
  10
  \@chnum
  \if \@nextchar m\thr@@\else
   \if \@nextchar p4 \else
    \if \@nextchar b5 \else
   \z@ \@chclass \z@ \@preamerr \z@ \fi \fi \fi \fi
   \fi \fi  \fi  \fi  \fi  \fi  \fi \fi \fi \fi \fi \fi}
\gdef\arraybackslash{\let\\=\@arraycr}
\def\@textsubscript#1{{\m@th\ensuremath{_{\mbox{\fontsize\sf@size\z@#1}}}}}
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\def\abbr{}
\def\corr{}
\def\expan{}
\def\gap{}
\def\orig{}
\def\reg{}
\def\ref{}
\def\sic{}
\def\persName{}\def\name{}
\def\placeName{}
\def\orgName{}
\def\textcal#1{{\fontspec{Lucida Calligraphy}#1}}
\def\textgothic#1{{\fontspec{Lucida Blackletter}#1}}
\def\textlarge#1{{\large #1}}
\def\textoverbar#1{\ensuremath{\overline{#1}}}
\def\textquoted#1{‘#1’}
\def\textsmall#1{{\small #1}}
\def\textsubscript#1{\@textsubscript{\selectfont#1}}
\def\textxi{\ensuremath{\xi}}
\def\titlem{\itshape}
\newenvironment{biblfree}{}{\ifvmode\par\fi }
\newenvironment{bibl}{}{}
\newenvironment{byline}{\vskip6pt\itshape\fontsize{16pt}{18pt}\selectfont}{\par }
\newenvironment{citbibl}{}{\ifvmode\par\fi }
\newenvironment{docAuthor}{\ifvmode\vskip4pt\fontsize{16pt}{18pt}\selectfont\fi\itshape}{\ifvmode\par\fi }
\newenvironment{docDate}{}{\ifvmode\par\fi }
\newenvironment{docImprint}{\vskip 6pt}{\ifvmode\par\fi }
\newenvironment{docTitle}{\vskip6pt\bfseries\fontsize{22pt}{25pt}\selectfont}{\par }
\newenvironment{msHead}{\vskip 6pt}{\par}
\newenvironment{msItem}{\vskip 6pt}{\par}
\newenvironment{rubric}{}{}
\newenvironment{titlePart}{}{\par }

\newcolumntype{L}[1]{){\raggedright\arraybackslash}p{#1}}
\newcolumntype{C}[1]{){\centering\arraybackslash}p{#1}}
\newcolumntype{R}[1]{){\raggedleft\arraybackslash}p{#1}}
\newcolumntype{P}[1]{){\arraybackslash}p{#1}}
\newcolumntype{B}[1]{){\arraybackslash}b{#1}}
\newcolumntype{M}[1]{){\arraybackslash}m{#1}}
\definecolor{label}{gray}{0.75}
\def\unusedattribute#1{\sout{\textcolor{label}{#1}}}
\DeclareRobustCommand*{\xref}{\hyper@normalise\xref@}
\def\xref@#1#2{\hyper@linkurl{#2}{#1}}
\begingroup
\catcode`\_=\active
\gdef_#1{\ensuremath{\sb{\mathrm{#1}}}}
\endgroup
\mathcode`\_=\string"8000
\catcode`\_=12\relax

\usepackage[a4paper,twoside,lmargin=1in,rmargin=1in,tmargin=1in,bmargin=1in,marginparwidth=0.75in]{geometry}
\usepackage{framed}

\definecolor{shadecolor}{gray}{0.95}
\usepackage{longtable}
\usepackage[normalem]{ulem}
\usepackage{fancyvrb}
\usepackage{fancyhdr}
\usepackage{graphicx}
\usepackage{marginnote}

\renewcommand{\@cite}[1]{#1}


\renewcommand*{\marginfont}{\itshape\footnotesize}

\def\Gin@extensions{.pdf,.png,.jpg,.mps,.tif}

  \pagestyle{fancy}

\usepackage[pdftitle={Vibrational Properties of Disorder Silicene},
 pdfauthor={}]{hyperref}
\hyperbaseurl{}

	 \paperwidth210mm
	 \paperheight297mm
              
\def\@pnumwidth{1.55em}
\def\@tocrmarg {2.55em}
\def\@dotsep{4.5}
\setcounter{tocdepth}{3}
\clubpenalty=8000
\emergencystretch 3em
\hbadness=4000
\hyphenpenalty=400
\pretolerance=750
\tolerance=2000
\vbadness=4000
\widowpenalty=10000

\renewcommand\section{\@startsection {section}{1}{\z@}%
     {-1.75ex \@plus -0.5ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large\bfseries}}
\renewcommand\subsection{\@startsection{subsection}{2}{\z@}%
     {-1.75ex\@plus -0.5ex \@minus- .2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\Large}}
\renewcommand\subsubsection{\@startsection{subsubsection}{3}{\z@}%
     {-1.5ex\@plus -0.35ex \@minus -.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\large}}
\renewcommand\paragraph{\@startsection{paragraph}{4}{\z@}%
     {-1ex \@plus-0.35ex \@minus -0.2ex}%
     {0.5ex \@plus .2ex}%
     {\reset@font\normalsize}}
\renewcommand\subparagraph{\@startsection{subparagraph}{5}{\parindent}%
     {1.5ex \@plus1ex \@minus .2ex}%
     {-1em}%
     {\reset@font\normalsize\bfseries}}


\def\l@section#1#2{\addpenalty{\@secpenalty} \addvspace{1.0em plus 1pt}
 \@tempdima 1.5em \begingroup
 \parindent \z@ \rightskip \@pnumwidth 
 \parfillskip -\@pnumwidth 
 \bfseries \leavevmode #1\hfil \hbox to\@pnumwidth{\hss #2}\par
 \endgroup}
\def\l@subsection{\@dottedtocline{2}{1.5em}{2.3em}}
\def\l@subsubsection{\@dottedtocline{3}{3.8em}{3.2em}}
\def\l@paragraph{\@dottedtocline{4}{7.0em}{4.1em}}
\def\l@subparagraph{\@dottedtocline{5}{10em}{5em}}
\@ifundefined{c@section}{\newcounter{section}}{}
\@ifundefined{c@chapter}{\newcounter{chapter}}{}
\newif\if@mainmatter 
\@mainmattertrue
\def\chaptername{Chapter}
\def\frontmatter{%
  \pagenumbering{roman}
  \def\thechapter{\@roman\c@chapter}
  \def\theHchapter{\roman{chapter}}
  \def\thesection{\@roman\c@section}
  \def\theHsection{\roman{section}}
  \def\@chapapp{}%
}
\def\mainmatter{%
  \cleardoublepage
  \def\thechapter{\@arabic\c@chapter}
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \pagenumbering{arabic}
  \setcounter{secnumdepth}{6}
  \def\@chapapp{\chaptername}%
  \def\theHchapter{\arabic{chapter}}
  \def\thesection{\@arabic\c@section}
  \def\theHsection{\arabic{section}}
}
\def\backmatter{%
  \cleardoublepage
  \setcounter{chapter}{0}
  \setcounter{section}{0}
  \setcounter{secnumdepth}{2}
  \def\@chapapp{\appendixname}%
  \def\thechapter{\@Alph\c@chapter}
  \def\theHchapter{\Alph{chapter}}
  \appendix
}
\newenvironment{bibitemlist}[1]{%
   \list{\@biblabel{\@arabic\c@enumiv}}%
       {\settowidth\labelwidth{\@biblabel{#1}}%
        \leftmargin\labelwidth
        \advance\leftmargin\labelsep
        \@openbib@code
        \usecounter{enumiv}%
        \let\p@enumiv\@empty
        \renewcommand\theenumiv{\@arabic\c@enumiv}%
	}%
  \sloppy
  \clubpenalty4000
  \@clubpenalty \clubpenalty
  \widowpenalty4000%
  \sfcode`\.\@m}%
  {\def\@noitemerr
    {\@latex@warning{Empty `bibitemlist' environment}}%
    \endlist}

\def\tableofcontents{\section*{\contentsname}\@starttoc{toc}}
\parskip0pt
\parindent1em
\def\Panel#1#2#3#4{\multicolumn{#3}{){\columncolor{#2}}#4}{#1}}
\newenvironment{reflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemsep}{0pt}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\itshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{sansreflist}{%
  \begin{raggedright}\begin{list}{}
  {%
   \setlength{\topsep}{0pt}%
   \setlength{\rightmargin}{0.25in}%
   \setlength{\itemindent}{0pt}%
   \setlength{\parskip}{0pt}%
   \setlength{\itemsep}{0pt}%
   \setlength{\parsep}{2pt}%
   \def\makelabel##1{\upshape ##1}}%
  }
  {\end{list}\end{raggedright}}
\newenvironment{specHead}[2]%
 {\vspace{20pt}\hrule\vspace{10pt}%
  \phantomsection\label{#1}\markright{#2}%

  \pdfbookmark[2]{#2}{#1}%
  \hspace{-0.75in}{\bfseries\fontsize{16pt}{18pt}\selectfont#2}%
  }{}
      \def\TheFullDate{2018-01-15 (revised: 15 January 2018)}
\def\TheID{\makeatother }
\def\TheDate{2018-01-15}
\title{Vibrational Properties of Disorder Silicene}
\author{}\makeatletter 
\makeatletter
\newcommand*{\cleartoleftpage}{%
  \clearpage
    \if@twoside
    \ifodd\c@page
      \hbox{}\newpage
      \if@twocolumn
        \hbox{}\newpage
      \fi
    \fi
  \fi
}
\makeatother
\makeatletter
\thispagestyle{empty}
\markright{\@title}\markboth{\@title}{\@author}
\renewcommand\small{\@setfontsize\small{9pt}{11pt}\abovedisplayskip 8.5\p@ plus3\p@ minus4\p@
\belowdisplayskip \abovedisplayskip
\abovedisplayshortskip \z@ plus2\p@
\belowdisplayshortskip 4\p@ plus2\p@ minus2\p@
\def\@listi{\leftmargin\leftmargini
               \topsep 2\p@ plus1\p@ minus1\p@
               \parsep 2\p@ plus\p@ minus\p@
               \itemsep 1pt}
}
\makeatother
\fvset{frame=single,numberblanklines=false,xleftmargin=5mm,xrightmargin=5mm}
\fancyhf{} 
\setlength{\headheight}{14pt}
\fancyhead[LE]{\bfseries\leftmark} 
\fancyhead[RO]{\bfseries\rightmark} 
\fancyfoot[RO]{}
\fancyfoot[CO]{\thepage}
\fancyfoot[LO]{\TheID}
\fancyfoot[LE]{}
\fancyfoot[CE]{\thepage}
\fancyfoot[RE]{\TheID}
\hypersetup{citebordercolor=0.75 0.75 0.75,linkbordercolor=0.75 0.75 0.75,urlbordercolor=0.75 0.75 0.75,bookmarksnumbered=true}
\fancypagestyle{plain}{\fancyhead{}\renewcommand{\headrulewidth}{0pt}}

\date{}
\usepackage{authblk}

\providecommand{\keywords}[1]
{
\footnotesize
  \textbf{\textit{Index terms---}} #1
}

\usepackage{graphicx,xcolor}
\definecolor{GJBlue}{HTML}{273B81}
\definecolor{GJLightBlue}{HTML}{0A9DD9}
\definecolor{GJMediumGrey}{HTML}{6D6E70}
\definecolor{GJLightGrey}{HTML}{929497} 

\renewenvironment{abstract}{%
   \setlength{\parindent}{0pt}\raggedright
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
   \textcolor{GJBlue}{\large\bfseries\abstractname\space}
}{%   
   \vskip8pt
   \textcolor{GJMediumGrey}{\rule{\textwidth}{2pt}}
   \vskip16pt
}

\usepackage[absolute,overlay]{textpos}

\makeatother 
      \usepackage{lineno}
      \linenumbers
      
\begin{document}

             \author[1]{Md Shamim  Mia}

\renewcommand\Authands{ and }

\date{\small \em Received: 12 December 2017 Accepted: 5 January 2018 Published: 15 January 2018}

\maketitle


\begin{abstract}
        


A large amount of work has been devoted to the calculation of the phonon density of states (PDOSs) in disordered systems. However, very little is known about the actual nature of phonon states. By analogy with the corresponding problem of electrons in disordered systems, one expect that for strong disorder, phonon states will be localized excitations in the Anderson sense .Phonon localization in one-dimensional solids has been considered by Ishii and Jackle .Numerical studies has been considered by Nagel, Rahman, and Grest, for Lenrand-Jones glasses The glasses were prepared by means of molecular dynamics simulations of rapid cooling from the liquid state. One a glass was formed, the normal modes of vibration were studied numerically. It was found that above a certain frequency the modes were localized, as expected from the general ideas of Anderson localization. The structural and vibrational properties of two-dimensional hexagonal silicon (silicene) are investigated by means of first-principles calculations. It is predicted that the silicene structure with a small buckling of 0.44 Å (0.7 Å) and bond lengths of 2.28 Å (2.44 Å) is energetically the most favourable, and it does not exhibit imaginary phonon mode.

\end{abstract}


\keywords{phonon density of state; defects; raman spectrum; lattice structure; brillouin zone.}

\begin{textblock*}{18cm}(1cm,1cm) % {block width} (coords) 
\textcolor{GJBlue}{\LARGE Global Journals \LaTeX\ JournalKaleidoscope\texttrademark}
\end{textblock*}

\begin{textblock*}{18cm}(1.4cm,1.5cm) % {block width} (coords) 
\textcolor{GJBlue}{\footnotesize \\ Artificial Intelligence formulated this projection for compatibility purposes from the original article published at Global Journals. However, this technology is currently in beta. \emph{Therefore, kindly ignore odd layouts, missed formulae, text, tables, or figures.}}
\end{textblock*}


\let\tabcellsep& 	 	 		 
\section[{I. Introduction}]{I. Introduction}\par
he field of nanomaterials has become one of the most quickly growing areas in science due to the unique properties and potential applications of these materials in electronics, medicine, consumer goods, defence, amongst others. Nanomaterials come in different shapes including zero-dimensional (0D), one-dimensional (1D) and two-dimensional (2D) forms. They are classified by having at least one of their dimensions less than 100 nm in size. The class of 2D nanomaterials are characterised by large lateral dimensions and small thicknesses of the order of less than, typically, several nanometres. They are often referred to as Nano sheets and can be considered akin to extremely thin sheets of paper. A wide variety of materials can be grown as 2D nanomaterials and can be composed of one or multiple elements. The elemental 2D nanomaterials usually have names ending in 'ene', namely silicene. Silicene is a crystalline twodimensional (2D) nanomaterial composed entirely of silicon (Si) atoms. The atoms in this single layer are arranged in a hexagonal honeycomb structure which, when viewed from the side, are buckled \hyperref[b0]{[1]}. In fact, in view of the unique electronic properties found for the widest studied 2D material, graphene, researchers are looking for similar or even better properties in other layered materials, which could be eventually integrated into the current Nano electronics technology. Silicon is nowadays the most important semiconducting element used in this technology. Thus, the idea of having the silicon counterpart of graphene (silicene), is at the same time fascinating from the physical point of view, but also very attractive by virtue its possible integration into Nano electronics devices. 
\section[{a) Objectives}]{a) Objectives}\par
We use force vibrational method for numerical analysis induced by Williams and Marris [121.s] and it is applicable for any disordered 2D material systems. We can apply this method as the field of ? A large complex disordered system. ? In low frequency regime we can calculate the DOS.\par
The main objectives of our thesis are:    ?? ? 1 = ????,\textbf{(1)}?? ? 2 = ? ?? 2 ???+ ?3?? 2 ??? (\textbf{2})?? ? 3 = ?? ?\textbf{(3)}\par
The distance of this buckling parameter (Î?"??) is often given as 0.44 °A. [17] While silicene may seem three-dimensional due to the buckling parameter, it is more like a pseudo two-dimensional structure. The structure is actually two-dimensional except for that buckling parameter. Vector ?? 3 is only present for calculation reasons in the case of the two-dimensional structures. The buckling amount will be verified in the calculations section. Notice the buckling amount will vary when introduced into the bilayer and super lattice structure. But the buckling amount of stand-alone silicene will be verified so that Quantum Espresso can be proven effective. 
\section[{d) Band Structure}]{d) Band Structure}\par
Once the lattice structure is determined, the next step is to determine the band diagrams of these silicene-graphene composites. To accomplish that one needs to construct the Brillouin zone of this hexagonal The high symmetry points in the Brillouin zone will be written in this frame of reference. The units of 2?? ?? are pulled out to the front.Brillouin zone of the two materials will be the same. ?? ?? 1 = 2?? ???3 ( ? 3?? ? ?? + ?? ? ??)\textbf{(4)}?? ?? 2 = 4?? ???3 ?? ? ??\textbf{(5)}?? ?? 3 = 2???? ? ??\textbf{(} 
\section[{e) Phonons}]{e) Phonons}\par
Silicene has a flexural mode, as does graphene. However, because of the symmetry reduction due to the buckling, this mode for silicene has both a z and a xy component. When silicene-based materials have been fabricated experimentally, they will have some vacancy and atoms-vacancy defects  {\ref [27]}. Vacancies have a striking effect on the electronic structure. In particular, it is well known that vacancies induce localized states with associated resonant peaks at the Dirac point, and that state in the vicinity of the Dirac point has an enhanced tendency for localization, as revealed by an enhancement of the electronic inverse participation ratio  {\ref [28]}.Variety of defects can exist in silicene and the impact on certain physical properties.\par
Single vacancies [33] and di-vacancies induce small gaps in silicene. Si atoms induce long-range spin polarization and a band gap, thus achieving an allsilicon magnetic semiconductor. Small defects were found to have a tendency to coalescence forming highly stable vacancy clusters  {\ref [34]}.Stone-Wales defect is a topological defect formed by the 90° rotation of a dimer bond which results in four hexagons turning into two pairs of pentagon-heptagon rings. In particular, it was found that the formation energy and kinetic barrier are lower in silicene than in graphene. A band gap of 0.01 eV is created. The effect of vacancies and Stone-Wales defects on the mechanical properties has been studied using molecular dynamic finite element method with Turnoff potential  {\ref [33]}. They found that pristine and lowly defective silicene sheets exhibit almost the same elastic nature up to the fracture points. However, a single defect significantly weakened the silicene sheet, leading to a considerable reduction in the fracture strength. Thus, one 2-atom vacancy in the reduced the fracture stress by 1820 \% and the fracture strain by 3335 \%. The weakening effects of Stone-Wales defects varied with the tensile direction and the orientation of these defects.\par
A defect that is particular to 2D materials but absent for graphene is the buckling interface formed between two pieces of silicene with oppositely oriented buckling  {\ref [34]}.This leads to a line defect that has a low formation energy but has a higher reactivity than the pristine silicene itself. The latter was deduced by studying the adsorption of a single gold atom and they turn out to have a binding energy of ?3.50 eV. 
\section[{II. Methodology a) Numerical method}]{II. Methodology a) Numerical method}\par
The equation of motion of the masses is  ?? ?? ?????(??) = ? ? ? ???? ????? ?(??) ?? ? (7)?? ?? (t)= ? ?? ?? ?? (??) ?? ?? (??) ??? ?? , (\textbf{8})\par
Where Q\textunderscore ? is the amplitude of the normal mode ? and e\textunderscore l (?) is the displacement pattern or "polarization" vector of the mode ?. The \{e\textunderscore l (?)\} and the frequencies \{?\textunderscore ? \} satisfy the equations? ?? ???? ?(?? ?? ?? ?? ?)? 1 2 ?? ????? ?(??) = ð??"ð??" ?? 2 ?? ?? (??)\textbf{(9)}\par
Thus, to find the frequencies and displacement patterns of the normal modes one has to find the eigenvalues and eigenvectors of an N×N matrix. Conventional methods require a large amount of computer time as N becomes large. Thus, one has to find another approach. 
\section[{b) Method for the Density of States}]{b) Method for the Density of States}\par
Start with each atom at rest and with zero displacement then apply a force on each atom given by ?? ?? cos(?t)Where ?? ?? is independent of time. After a time, large compared to the typical period of oscillation of the atoms, the total energy of the atom is,?? = 1 2 ? [? ?? ?? ?? ?? (??) ??? ?? ?? ] 2 ?? sin \{[(ð??"ð??" ?? ??)/2]??\} 2 (ð??"ð??" ?? ??) 2\textbf{(10)}\par
Thus, for large times the periodic force excites only those modes whose frequency is close to ?.?? ?? = ?? 0 ??? ?? cos ? ??\textbf{(11)}\par
Where ?? 0 is a constant and ?? ?? is a random quantity. If we average over all possible values of ? ?? we find that the average value of ?? is????= ?? 0 2 4 ? sin \{[(ð??"ð??" ?? ??)/2]??\} 2 (ð??"ð??" ?? ??) 2 ?? (12)\par
We have used the orthonormality of the mode patterns:? ?? ?? (??)?? ?? ??? ?? = ?? ???? ??\textbf{(13)}\par
The modes which contribute to the sum in equation 3.6 are those whose frequencies lie within about ± ?2?? ?? ? of ? . Let us choose t such that ??? ? 1 , (\textbf{14}) 4???? ð??"ð??" ?? ?? ? 1 (15)\par
Where ? ?? is the maximum frequency of the system. Equation (3.8) means that only modes in a narrow band of frequency on the scale of ? contribution to the sum, equation (3.9) ensures that the number of such modes is much larger than unity. Thus, if there conditions are met we have???? ? ???? ?? 0 2 8 ? ??(ð??"ð??" ?? ? ?) ?? = ???? ?? ?? 2 ????(?) 8 (\textbf{16})\par
Where g(?) is the phonon density of states. Thus, provided that a way can be found to carry out the time development in the presence of the periodic force, we can findg(?) = 8???? ???? ?? ?? 2 (17)\par
For a sufficiently large system, one expects that it will not be necessary to carry out explicitly the average over all values of the make one random choice of the \{?? ?? \}. 
\section[{c) Numerical Algorithm}]{c) Numerical Algorithm}\par
The problem thus reduces to the solution of the equations of motion of the system in the presence of a periodic force. These equations can be written as,?? ?? = ????? ?1 [?? ???? . ?? ?? ?(??) + ?? ?? cos(???)]\textbf{(18)}?????(??) = ?? ?? (??)\textbf{(19)}\par
The standard approach to the time development is to replace t by n?? where ?? is a small time step, and n an integer. Then a time development algorithm is?? 1 (?? + 1) = ?? ?? (??) + ?? ?? ?1 [?? ???? . ?? ?? ?(??) + ?? ?? ??????(????)]?? (20) ?? 1 (?? + 1) = ?? ?? + ?? ?? (??) ??\textbf{(21)}\par
Normally one would choose ?? to be small scale of the system, i.e. the period ?? ?? of the highest frequency mode.\par
The initial set at t= 0 is prepared in which all atoms are at rest and have zero displacements. A random force is applied to each atom at rest, which is given by:?? ?? cos(???) = ?? 0 ??? ?? cos(????) cos(???)\textbf{(22)}\par
Where ?? 0 is a constant and ???? varies from 0 to 2??. The corresponding equation of motion of the system becomes:?? ?? ????? + ? ????? ?? ????? ?(??) = ?? 0 ??? ?? cos(???) cos(???) (23)\par
The total energy of the system can be written as the summation of kinetic and potential energy as follows:??(??) = ?? + ?? = 1 2 ? ?? ?? ?? ?? ?? ?2(??) + 1 2 ? ? ?? ?? (??) ?? ?? ????? ????? ?(??) = 1 2 (? ?? ?? ?? ?? ?? ?2(??) + ? ? ?? ?? (??) ?? ?? ????? ????? ?(??))\textbf{(24)}\par
Thus, for large times the periodic external force excites only those modes whose frequency is close to. When average all possible values of ?? ?? and used the ortho-normality of the eigen-vector??? ?? (??), the average value of energy ????becomes, Therefore, from Eq. 3.10, we can write:???(?)?= ???? ?? 0 2 8 ? ??(ð??"ð??" ?? ? ?) ?? (27)\par
The density of states are related with the delta function by the following equation:??(?)= 1 ?? ? ??(ð??"ð??" ?? ? ?) ?? (28)\par
Therefore, we can write the density of the states of the system is???(?)?= ???? ?? 0 2 ????(?) 8\textbf{(29)}??(?) = 8???(?)? ???? ?? 0 2 ??\textbf{(30)}\par
Thus, it could be obtained phonon density of states to carry out the time development in the presence of external force. 
\section[{d) Raman Spectroscopy}]{d) Raman Spectroscopy}\par
Raman spectroscopy is a widely used experimental technique to identify the characteristic vibrational modes of materials. It is one of the most used optical techniques for the characterization of 2D materials like graphene. Several properties can be derived from the analysis of the Raman spectra, e.g., the shift in frequency of the Raman peaks has been successfully used to analysis flakes of layered materials such as ?MoS?\textunderscore 2and quantify the number of layers  {\ref [40,} {\ref 41]}; similarly the Raman spectrum of graphene flakes has been used not only for the identification of the number of layers, but also to estimate the quality (in term of defects) of the graphene sheets  {\ref [42,} {\ref 43]}.\par
Note however that the splitting of the iTO and iLO at Î?" could be due to an artifact of the DFT calculation and particularly to the slightly different bond length a ? b predicted by using the nonlinear core correction. Nevertheless, the iTO mode is very close to the LiO mode and its Raman spectra show only one peak at around 575 cm-1, as shown in Fig. \hyperref[fig_2]{3}.1 below.  
\section[{III. Result And Discussion}]{III. Result And Discussion}\par
We have calculated phonon density of states of silicene, g(?) for lattice of 1000 atoms in the array respectively 40×40 and we analysis for 5\%,10\% and 15\% defects.We also have calculated phono density of states of silicene, g(?) for lattice of 2000 atoms in the array respectively 50×50. We consider same 5\%,10\% and 15\% defects for this lattice structure. 
\section[{a) Force constants calculation}]{a) Force constants calculation}\par
As the force constants are effective in three planes such that in the radius plane, in the inner plane, and in the outer plane. The force constants are as given following table  {\ref 4}.  {\ref 1[44]}.     
\section[{Conclusions And Future Work}]{Conclusions And Future Work}\par
Silicene is a single layer of silicon atoms with sp3 bonds in the honeycomb lattice structure. More compatibility of the Silicene with the current semiconducting technology and some other advantages such as a new and promising alternative for the spintronic and the Nano-magnetic materials have attracted considerable scientific attention in very recent years. In addition, despite the single layer Graphene which owns zero band gap in the Dirac point, owns a tunable band gap from zero to semiconducting region, which make it better choice for the FET technology. However, stability of this material was an We studied the Phonon Density of states (PDOSs) of disorder silicene, we belief that this thesis work will help to investigate the electronic properties of silicene. Of course, the work on silicene is still a growing area of research and there are many questions yet to be answered about the potential of this material. We are very excited to be part of the discoveries in this area and look forward to seeing where this material will go in the future. We hope that you as the reader will share in our passion for science and 2D nanomaterials!\par
The simulating analysis of vibrational properties like phonon density of states have been studied and found phonon density of states for lattice atoms disorder. There have been considered the normal mode of vibration. At the same time, there have been several important issues in the field that remain to be addressed.\par
We have performed the mathematical calculation on basis of a 2d dimensional model. The method developed in this dissertation can be extended to the three dimension for calculating the vibrational properties of multi-layer silicene. We can further calculate the layered structure of silicene and, further investigation of many-body interactions and electronphonon scattering effect, calculation of Raman intensities, quantitative analysis of specific structural defects such as Stone-walse defects in silicene samples, and an harmonic effects responsible for the thermal conductivity in the silicene sheet. The electronphonon interactions can be further calculated. The harmonic effects and thermal conductivity in the silicene sheet can be studied in a more systematic way. issue that the recent experimental studies showed that it could be stable with a small buckling. However, its electronic properties such as the density of states, the total carriers and the conductance as well as their temperature dependence and comparison with the electronic properties of the Graphene have been being investigating analytically.\begin{figure}[htbp]
\noindent\textbf{1}\includegraphics[]{image-2.png}
\caption{\label{fig_0}Figure 1 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{2}\includegraphics[]{image-3.png}
\caption{\label{fig_1}Figure 2 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{3}\includegraphics[]{image-4.png}
\caption{\label{fig_2}Figure 3 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-5.png}
\caption{\label{fig_3}6 )F}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{5}\includegraphics[]{image-6.png}
\caption{\label{fig_4}Figure 5 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{6}\includegraphics[]{image-7.png}
\caption{\label{fig_5}Figure 6 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{7}\includegraphics[]{image-8.png}
\caption{\label{fig_6}Figure 7 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{}\includegraphics[]{image-9.png}
\caption{\label{fig_7}}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{4}\includegraphics[]{image-10.png}
\caption{\label{fig_8}Figure 4 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{42}\includegraphics[]{image-11.png}
\caption{\label{fig_9}4 ? 2 (}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{8}\includegraphics[]{image-12.png}
\caption{\label{fig_10}Figure 8 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{9}\includegraphics[]{image-13.png}
\caption{\label{fig_11}Figure 9 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{10}\includegraphics[]{image-14.png}
\caption{\label{fig_12}Figure 10 :}\end{figure}
 \begin{figure}[htbp]
\noindent\textbf{1314}\includegraphics[]{image-15.png}
\caption{\label{fig_13}Figure 13 :Figure 14 :}\end{figure}
  \begin{figure}[htbp]
\noindent\textbf{1} \par 
\begin{longtable}{P{0.12395833333333334\textwidth}P{0.33055555555555555\textwidth}P{0.2125\textwidth}P{0.18298611111111113\textwidth}}
Neighbour\tabcellsep ?? ð??"ð??"ð??"ð??"\tabcellsep ?? ??????\tabcellsep ?? ????\\
1st\tabcellsep 2.032544\tabcellsep 15.9965\tabcellsep 0.3814\\
2nd\tabcellsep -0.882486\tabcellsep 0.9010\tabcellsep 0.0683\\
3rd\tabcellsep 0.249645\tabcellsep -0.9737\tabcellsep 0.1396\\
4th\tabcellsep 0.295442\tabcellsep -0.1067\tabcellsep 0.1006\end{longtable} \par
  {\small\itshape [Note: b) ResultsPhonon Density of states of Silicene for different atom numbers:]} 
\caption{\label{tab_2}Table 1 :}\end{figure}
 			\footnote{© 2018 Global Journals} 		 		\backmatter  			  				\begin{bibitemlist}{1}
\bibitem[Lee et al. ()]{b7}\label{b7} 	 		‘A study on the nonlinear dynamics of PR interval variability using surrogate data’.  		 			Jongmin Lee 		,  		 			Kwangsuk Park 		,  		 			Insun Shin 		.  	 	 		\textit{18th Annual International Conference of the IEEE Engineering in Medicine and Biology Society},  				 (Amsterdam)  		1996.  	 
\bibitem[Thakor and Zhu (1991)]{b4}\label{b4} 	 		‘Applications of adaptive filtering to ECG analysis: Noise cancellation and arrhythmia detection’.  		 			N V Thakor 		,  		 			Y Zhu 		.  	 	 		\textit{IEEE Trans. Biomed. Eng}  		Aug. 1991. 38 p. .  	 
\bibitem[Uzun et al. (2001)]{b8}\label{b8} 	 		‘Nonlinear analysis of heart rate variability’.  		 			I S Uzun 		,  		 			M H Asyali 		,  		 			G Celebi 		,  		 			M Pehlivan 		.  	 	 		\textit{Proceedings of the 23 rd Annual International Conference of the IEEE},  				 (the 23 rd Annual International Conference of the IEEE)  		2001. Oct. 2001. 2 p. .  	 
\bibitem[Signorini et al. (2006)]{b0}\label{b0} 	 		‘Nonlinear analysis of Heart Rate Variability signal for the characterization of Cardiac Heart Failure patients’.  		 			M G Signorini 		,  		 			M Ferrario 		,  		 			M Marchetti 		,  		 			A Marseglia 		.  	 	 		\textit{Proceedings of the 28th IEEE, EMBS Annual International Conference},  				 (the 28th IEEE, EMBS Annual International ConferenceNew York City, USA)  		Aug 30-Sept 3, 2006.  	 
\bibitem[Minami et al. (1999)]{b2}\label{b2} 	 		‘Realtime discrimination of ventricular tachyarrhythmia with Fourier-transform neural network’.  		 			K Minami 		,  		 			H Nakajima 		,  		 			T Toyoshima 		.  	 	 		\textit{IEEE Trans. Biomed. Eng}  		Feb. 1999. 46 p. .  	 
\bibitem[Signorhi et al. ()]{b1}\label{b1} 	 		‘Regularity Patterns in heart rate variability signal: the Approximate Entropy Approach’.  		 			Maria G Signorhi 		,  		 			Roberto Sassi 		,  		 			Federico Lombardi 		,  		 			Sergio Cerutti 		.  	 	 		\textit{Proceedings of the 20 th Annual International Conference of the IEEE Engineering in Medicine and Biology Society},  				 (the 20 th Annual International Conference of the IEEE Engineering in Medicine and Biology Society)  		1998. 20.  	 
\bibitem[Owis et al. ()]{b6}\label{b6} 	 		‘Study of Features Based on Nonlinear Dynamical Modelingin ECG Arrhythmia Detection and Classification’.  		 			Mohamed I Owis 		,  		 			Ahmed H Abou-Zied 		,  		 			M Abou-Bakr 		,  		 			Yasser M Youssef 		,  		 			Kadah 		.  	 	 		\textit{IEEE Trans. on Biomedical Engineering}  		JULY 2002. 49  (7) .  	 
\bibitem[Thakor and Pan (1990)]{b3}\label{b3} 	 		‘Tachycardia and fibrillation detection by automatic implantable cardioverterdefibrillators: sequential testing in time domain’.  		 			N V Thakor 		,  		 			K Pan 		.  	 	 		\textit{IEEE Trans. Biomed. Engg}  		March.1990. 09 p. .  	 
\bibitem[Jun et al.]{b9}\label{b9} 	 		‘Time series prediction using Lyapunov exponents in embedding phase space’.  		 			Jun 		,  		 			K F Zhang 		,  		 			Man 		.  	 	 		\textit{Proceedings of ICSP "98},  				 (ICSP "98)  		 	 
\bibitem[Thakor et al. (1990)]{b5}\label{b5} 	 		‘Ventricular tachycardia and fibrillation detection by a sequential hypothesis testing algorithm’.  		 			N V Thakor 		,  		 			Y Zhu 		,  		 			K Pan 		.  	 	 		\textit{IEEE Trans. Biomed. Eng}  		Sept. 1990. 37 p. .  	 
\end{bibitemlist}
 			 		 	 
\end{document}
