after merge

This commit is contained in:
jcuny 2021-06-15 19:42:03 +02:00
commit debb6f898c
15 changed files with 7470 additions and 7614 deletions

File diff suppressed because it is too large Load Diff

View File

@ -127,7 +127,7 @@
\newlabel{Ehartree}{{2.21}{24}{Density Functional Theory}{equation.2.3.21}{}}
\newlabel{Exc}{{2.22}{24}{Density Functional Theory}{equation.2.3.22}{}}
\newlabel{Veff}{{2.24}{24}{Density Functional Theory}{equation.2.3.24}{}}
\newlabel{rho}{{2.25}{24}{Density Functional Theory}{equation.2.3.25}{}}
\newlabel{KS}{{2.25}{24}{Density Functional Theory}{equation.2.3.25}{}}
\citation{Kohn1965}
\citation{Kohn1965,Herman1969,Perdew1996}
\citation{Becke1993}
@ -136,26 +136,26 @@
\citation{Henderson2009,Krukau2006}
\citation{Zhao2008}
\citation{Li2018}
\newlabel{KS}{{2.26}{25}{Density Functional Theory}{equation.2.3.26}{}}
\newlabel{orthonormality}{{2.27}{25}{Density Functional Theory}{equation.2.3.27}{}}
\@writefile{brf}{\backcite{Kohn1965}{{25}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Kohn1965}{{25}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Herman1969}{{25}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Perdew1996}{{25}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Becke1993}{{25}{2.3.2}{equation.2.3.27}}}
\newlabel{Ehybrid}{{2.28}{25}{Density Functional Theory}{equation.2.3.28}{}}
\newlabel{orthonormality}{{2.26}{25}{Density Functional Theory}{equation.2.3.26}{}}
\@writefile{brf}{\backcite{Kohn1965}{{25}{2.3.2}{equation.2.3.26}}}
\@writefile{brf}{\backcite{Kohn1965}{{25}{2.3.2}{equation.2.3.26}}}
\@writefile{brf}{\backcite{Herman1969}{{25}{2.3.2}{equation.2.3.26}}}
\@writefile{brf}{\backcite{Perdew1996}{{25}{2.3.2}{equation.2.3.26}}}
\@writefile{brf}{\backcite{Becke1993}{{25}{2.3.2}{equation.2.3.26}}}
\newlabel{Ehybrid}{{2.27}{25}{Density Functional Theory}{equation.2.3.27}{}}
\citation{Foulkes1989}
\citation{Porezag1995}
\citation{Elstner1998,Porezag1995,Seifert1996,Elstner1998,Elstne2007,Oliveira2009}
\@writefile{brf}{\backcite{Becke1993}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Becke1988}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Lee1988}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Stephens1994}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Adamo1999}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Henderson2009}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Krukau2006}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Zhao2008}{{26}{2.3.2}{equation.2.3.28}}}
\@writefile{brf}{\backcite{Li2018}{{26}{2.3.2}{equation.2.3.28}}}
\citation{Gaus2011}
\@writefile{brf}{\backcite{Becke1993}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Becke1988}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Lee1988}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Stephens1994}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Adamo1999}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Henderson2009}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Krukau2006}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Zhao2008}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{brf}{\backcite{Li2018}{{26}{2.3.2}{equation.2.3.27}}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.3.3}Density Functional based Tight-Binding Theory}{26}{subsection.2.3.3}}
\newlabel{sec:DFTB}{{2.3.3}{26}{Density Functional based Tight-Binding Theory}{subsection.2.3.3}{}}
\@writefile{brf}{\backcite{Foulkes1989}{{26}{2.3.3}{subsection.2.3.3}}}
@ -166,67 +166,69 @@
\@writefile{brf}{\backcite{Seifert1996}{{26}{2.3.3}{subsection.2.3.3}}}
\@writefile{brf}{\backcite{Elstne2007}{{26}{2.3.3}{subsection.2.3.3}}}
\@writefile{brf}{\backcite{Oliveira2009}{{26}{2.3.3}{subsection.2.3.3}}}
\citation{Gaus2011}
\@writefile{brf}{\backcite{Gaus2011}{{27}{2.3.3}{subsection.2.3.3}}}
\newlabel{Edft}{{2.29}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.29}{}}
\newlabel{Edftb}{{2.30}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.30}{}}
\newlabel{Edftb-long}{{2.31}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.31}{}}
\newlabel{H0}{{2.32}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.32}{}}
\newlabel{Eterms}{{2.33}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.33}{}}
\newlabel{orbitals}{{2.34}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.34}{}}
\newlabel{Eband}{{2.35}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.35}{}}
\newlabel{matrix}{{2.36}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.36}{}}
\newlabel{Vdftb}{{2.37}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.37}{}}
\newlabel{Hamil}{{2.38}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.38}{}}
\newlabel{Hmatrix}{{2.39}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.39}{}}
\newlabel{Erep-infinity}{{2.40}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.40}{}}
\newlabel{Erep-right}{{2.41}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.41}{}}
\newlabel{Erep-final}{{2.42}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.42}{}}
\newlabel{Edft}{{2.28}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.28}{}}
\newlabel{Edftb}{{2.29}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.29}{}}
\newlabel{Edftb-long}{{2.30}{27}{Density Functional based Tight-Binding Theory}{equation.2.3.30}{}}
\newlabel{H0}{{2.31}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.31}{}}
\newlabel{Eterms}{{2.32}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.32}{}}
\newlabel{orbitals}{{2.33}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.33}{}}
\newlabel{Eband}{{2.34}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.34}{}}
\newlabel{matrix}{{2.35}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.35}{}}
\newlabel{Vdftb}{{2.36}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.36}{}}
\newlabel{Hamil}{{2.37}{28}{Density Functional based Tight-Binding Theory}{equation.2.3.37}{}}
\citation{Mulliken1955}
\newlabel{Hmatrix}{{2.38}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.38}{}}
\newlabel{relative_p}{{2.39}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.39}{}}
\newlabel{relative_pa}{{2.40}{29}{Density Functional based Tight-Binding Theory}{equation.2.3.40}{}}
\@writefile{brf}{\backcite{Mulliken1955}{{29}{2.3.3}{equation.2.3.40}}}
\citation{Elstner1998}
\citation{Mulliken1955}
\newlabel{relative_p}{{2.43}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.43}{}}
\newlabel{relative_pa}{{2.44}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.44}{}}
\@writefile{brf}{\backcite{Mulliken1955}{{30}{2.3.3}{equation.2.3.44}}}
\newlabel{E2nd}{{2.46}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.46}{}}
\newlabel{E2ndsimple}{{2.47}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.47}{}}
\@writefile{brf}{\backcite{Mulliken1955}{{31}{2.3.3}{equation.2.3.47}}}
\newlabel{Mulliken}{{2.48}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.48}{}}
\newlabel{Escc}{{2.49}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.49}{}}
\newlabel{secular-eq}{{2.50}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.50}{}}
\newlabel{Huv}{{2.51}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.51}{}}
\newlabel{E2nd}{{2.41}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.41}{}}
\newlabel{simpleE2nd}{{2.42}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.42}{}}
\@writefile{brf}{\backcite{Elstner1998}{{30}{2.3.3}{equation.2.3.42}}}
\@writefile{brf}{\backcite{Mulliken1955}{{30}{2.3.3}{equation.2.3.42}}}
\newlabel{Mulliken}{{2.43}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.43}{}}
\newlabel{Erep-right}{{2.44}{30}{Density Functional based Tight-Binding Theory}{equation.2.3.44}{}}
\citation{Bader1985,Bader1990,Glendening1998,Glendening2012}
\citation{Singh1984,Besler1990}
\citation{Li1998,Kalinowski2004,Kelly2005,Rapacioli2009corr}
\newlabel{Erep-final}{{2.45}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.45}{}}
\newlabel{Escc}{{2.46}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.46}{}}
\newlabel{secular-eq}{{2.47}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.47}{}}
\newlabel{Huv}{{2.48}{31}{Density Functional based Tight-Binding Theory}{equation.2.3.48}{}}
\citation{Mayer1983,Cances1997,Lu2012}
\citation{Lewis1995,Gianturco1999,Elstner2001,Wu2002,Grimme2004,Zimmerli2004,Neumann2005,Goursot2007}
\@writefile{brf}{\backcite{Bader1985}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Bader1990}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Glendening1998}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Glendening2012}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Singh1984}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Besler1990}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Li1998}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Kalinowski2004}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Kelly2005}{{32}{2.3.3}{equation.2.3.51}}}
\@writefile{brf}{\backcite{Rapacioli2009corr}{{32}{2.3.3}{equation.2.3.51}}}
\newlabel{CM3}{{2.52}{32}{Density Functional based Tight-Binding Theory}{equation.2.3.52}{}}
\@writefile{brf}{\backcite{Mayer1983}{{32}{2.3.3}{equation.2.3.52}}}
\@writefile{brf}{\backcite{Cances1997}{{32}{2.3.3}{equation.2.3.52}}}
\@writefile{brf}{\backcite{Lu2012}{{32}{2.3.3}{equation.2.3.52}}}
\newlabel{CM3-classical}{{2.53}{32}{Density Functional based Tight-Binding Theory}{equation.2.3.53}{}}
\@writefile{brf}{\backcite{Lewis1995}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Gianturco1999}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Elstner2001}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Wu2002}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Grimme2004}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Zimmerli2004}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Neumann2005}{{32}{2.3.3}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Goursot2007}{{32}{2.3.3}{equation.2.3.53}}}
\citation{Frenkel2001}
\@writefile{brf}{\backcite{Bader1985}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Bader1990}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Glendening1998}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Glendening2012}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Singh1984}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Besler1990}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Li1998}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Kalinowski2004}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Kelly2005}{{32}{2.3.3}{equation.2.3.48}}}
\@writefile{brf}{\backcite{Rapacioli2009corr}{{32}{2.3.3}{equation.2.3.48}}}
\newlabel{CM3}{{2.49}{32}{Density Functional based Tight-Binding Theory}{equation.2.3.49}{}}
\@writefile{brf}{\backcite{Mayer1983}{{32}{2.3.3}{equation.2.3.49}}}
\@writefile{brf}{\backcite{Cances1997}{{32}{2.3.3}{equation.2.3.49}}}
\@writefile{brf}{\backcite{Lu2012}{{32}{2.3.3}{equation.2.3.49}}}
\newlabel{CM3-classical}{{2.50}{32}{Density Functional based Tight-Binding Theory}{equation.2.3.50}{}}
\@writefile{brf}{\backcite{Lewis1995}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Gianturco1999}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Elstner2001}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Wu2002}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Grimme2004}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Zimmerli2004}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Neumann2005}{{32}{2.3.3}{equation.2.3.50}}}
\@writefile{brf}{\backcite{Goursot2007}{{32}{2.3.3}{equation.2.3.50}}}
\newlabel{dispersionE}{{2.51}{32}{Density Functional based Tight-Binding Theory}{equation.2.3.51}{}}
\citation{Jones1924,Lennard1924,Morse1929,Born1932,Stillinger1985,Kaplan2006}
\citation{Bahar1997}
\citation{Morse1929,Girifalco1959}
\newlabel{dispersionE}{{2.54}{33}{Density Functional based Tight-Binding Theory}{equation.2.3.54}{}}
\citation{Jones1924,Lennard1924}
\citation{Lennard1931}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.3.4}Force Field Methods}{33}{subsection.2.3.4}}
\@writefile{brf}{\backcite{Frenkel2001}{{33}{2.3.4}{subsection.2.3.4}}}
\@writefile{brf}{\backcite{Jones1924}{{33}{2.3.4}{subsection.2.3.4}}}
@ -236,43 +238,41 @@
\@writefile{brf}{\backcite{Stillinger1985}{{33}{2.3.4}{subsection.2.3.4}}}
\@writefile{brf}{\backcite{Kaplan2006}{{33}{2.3.4}{subsection.2.3.4}}}
\@writefile{brf}{\backcite{Bahar1997}{{33}{2.3.4}{subsection.2.3.4}}}
\newlabel{harmonic}{{2.55}{33}{Force Field Methods}{equation.2.3.55}{}}
\citation{Jones1924,Lennard1924}
\citation{Lennard1931}
\newlabel{harmonic}{{2.52}{33}{Force Field Methods}{equation.2.3.52}{}}
\@writefile{brf}{\backcite{Morse1929}{{33}{2.3.4}{equation.2.3.52}}}
\@writefile{brf}{\backcite{Girifalco1959}{{33}{2.3.4}{equation.2.3.52}}}
\newlabel{Morse}{{2.53}{33}{Force Field Methods}{equation.2.3.53}{}}
\@writefile{brf}{\backcite{Jones1924}{{33}{2.3.4}{equation.2.3.53}}}
\@writefile{brf}{\backcite{Lennard1924}{{33}{2.3.4}{equation.2.3.53}}}
\citation{Monticelli2013}
\@writefile{brf}{\backcite{Morse1929}{{34}{2.3.4}{equation.2.3.55}}}
\@writefile{brf}{\backcite{Girifalco1959}{{34}{2.3.4}{equation.2.3.55}}}
\newlabel{Morse}{{2.56}{34}{Force Field Methods}{equation.2.3.56}{}}
\@writefile{brf}{\backcite{Jones1924}{{34}{2.3.4}{equation.2.3.56}}}
\@writefile{brf}{\backcite{Lennard1924}{{34}{2.3.4}{equation.2.3.56}}}
\@writefile{brf}{\backcite{Lennard1931}{{34}{2.3.4}{equation.2.3.56}}}
\newlabel{LennarsJones}{{2.57}{34}{Force Field Methods}{equation.2.3.57}{}}
\@writefile{brf}{\backcite{Monticelli2013}{{34}{2.3.4}{equation.2.3.57}}}
\newlabel{potential}{{2.58}{34}{Force Field Methods}{equation.2.3.58}{}}
\@writefile{brf}{\backcite{Lennard1931}{{34}{2.3.4}{equation.2.3.53}}}
\newlabel{LennarsJones}{{2.54}{34}{Force Field Methods}{equation.2.3.54}{}}
\@writefile{brf}{\backcite{Monticelli2013}{{34}{2.3.4}{equation.2.3.54}}}
\newlabel{potential}{{2.55}{34}{Force Field Methods}{equation.2.3.55}{}}
\@writefile{toc}{\contentsline {section}{\numberline {2.4}Exploration of PES}{35}{section.2.4}}
\@writefile{lof}{\contentsline {figure}{\numberline {2.3}{\ignorespaces Schematic representation of some key points on a model potential energy surface.}}{35}{figure.caption.7}}
\newlabel{PES}{{2.3}{35}{Schematic representation of some key points on a model potential energy surface}{figure.caption.7}{}}
\citation{Metropolis1987}
\citation{Metropolis1949}
\citation{Kroese2014}
\@writefile{lof}{\contentsline {figure}{\numberline {2.3}{\ignorespaces Schematic representation of some key points on a model potential energy surface.}}{36}{figure.caption.7}}
\newlabel{PES}{{2.3}{36}{Schematic representation of some key points on a model potential energy surface}{figure.caption.7}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.4.1}Monte Carlo Simulations}{36}{subsection.2.4.1}}
\@writefile{brf}{\backcite{Metropolis1987}{{36}{2.4.1}{subsection.2.4.1}}}
\@writefile{brf}{\backcite{Metropolis1949}{{36}{2.4.1}{subsection.2.4.1}}}
\@writefile{brf}{\backcite{Kroese2014}{{36}{2.4.1}{subsection.2.4.1}}}
\newlabel{SingleIntegral}{{2.56}{36}{Monte Carlo Simulations}{equation.2.4.56}{}}
\newlabel{SingleAverage}{{2.57}{37}{Monte Carlo Simulations}{equation.2.4.57}{}}
\newlabel{limI}{{2.58}{37}{Monte Carlo Simulations}{equation.2.4.58}{}}
\newlabel{IntegralPx}{{2.59}{37}{Monte Carlo Simulations}{equation.2.4.59}{}}
\newlabel{INp}{{2.60}{37}{Monte Carlo Simulations}{equation.2.4.60}{}}
\newlabel{averaA}{{2.61}{37}{Monte Carlo Simulations}{equation.2.4.61}{}}
\citation{Rosenbluth1955,Binder1993,Baeurle2009}
\@writefile{brf}{\backcite{Rosenbluth1955}{{37}{2.4.1}{subsection.2.4.1}}}
\@writefile{brf}{\backcite{Binder1993}{{37}{2.4.1}{subsection.2.4.1}}}
\@writefile{brf}{\backcite{Baeurle2009}{{37}{2.4.1}{subsection.2.4.1}}}
\newlabel{thermalAverage}{{2.59}{38}{Monte Carlo Simulations}{equation.2.4.59}{}}
\newlabel{SingleIntegral}{{2.60}{38}{Monte Carlo Simulations}{equation.2.4.60}{}}
\newlabel{SingleAverage}{{2.61}{38}{Monte Carlo Simulations}{equation.2.4.61}{}}
\newlabel{limI}{{2.62}{38}{Monte Carlo Simulations}{equation.2.4.62}{}}
\newlabel{IntegralPx}{{2.63}{39}{Monte Carlo Simulations}{equation.2.4.63}{}}
\newlabel{INp}{{2.64}{39}{Monte Carlo Simulations}{equation.2.4.64}{}}
\newlabel{averaA}{{2.65}{39}{Monte Carlo Simulations}{equation.2.4.65}{}}
\newlabel{Markov1}{{2.66}{39}{Monte Carlo Simulations}{equation.2.4.66}{}}
\newlabel{Markov2}{{2.67}{39}{Monte Carlo Simulations}{equation.2.4.67}{}}
\newlabel{balanceEq}{{2.68}{39}{Monte Carlo Simulations}{equation.2.4.68}{}}
\newlabel{Markov1}{{2.62}{38}{Monte Carlo Simulations}{equation.2.4.62}{}}
\newlabel{Markov2}{{2.63}{38}{Monte Carlo Simulations}{equation.2.4.63}{}}
\newlabel{balanceEq}{{2.64}{38}{Monte Carlo Simulations}{equation.2.4.64}{}}
\newlabel{assess}{{2.65}{38}{Monte Carlo Simulations}{equation.2.4.65}{}}
\@writefile{brf}{\backcite{Rosenbluth1955}{{38}{2.4.1}{equation.2.4.65}}}
\@writefile{brf}{\backcite{Binder1993}{{38}{2.4.1}{equation.2.4.65}}}
\@writefile{brf}{\backcite{Baeurle2009}{{38}{2.4.1}{equation.2.4.65}}}
\citation{Alder1957}
\citation{Gibson1960,Rahman1964}
\citation{Allen2017,Cicccotti1987,Van1990,Berne1998,Tuckerman2000,Frenkel2001,Karplus2002,Rapaport2004,Tavan2005,Van2006}
@ -281,10 +281,10 @@
\citation{Charati1998,Yeh2004,Braga2014,Pranami2015}
\citation{Mondello1997,Hess2002,Yeh2004,Zhang2015}
\citation{Aragon1976,Vuilleumier1999,Pavone2006,Brancato2006}
\newlabel{thermalAverage}{{2.66}{39}{Monte Carlo Simulations}{equation.2.4.66}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.4.2}Classical Molecular Dynamics}{39}{subsection.2.4.2}}
\citation{Beck2000,Wang2003}
\citation{Wei1994,Light2000}
\newlabel{assess}{{2.69}{40}{Monte Carlo Simulations}{equation.2.4.69}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.4.2}Classical Molecular Dynamics}{40}{subsection.2.4.2}}
\@writefile{brf}{\backcite{Alder1957}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Gibson1960}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Rahman1964}{{40}{2.4.2}{subsection.2.4.2}}}
@ -319,48 +319,54 @@
\@writefile{brf}{\backcite{Vuilleumier1999}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Pavone2006}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Brancato2006}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Beck2000}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Wang2003}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Wei1994}{{40}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Light2000}{{40}{2.4.2}{subsection.2.4.2}}}
\citation{Batina1991,Butcher2008}
\citation{Gear1966,Diethelm2002}
\citation{Runge1895,Kutta1901,Butcher1996,Butcher2015}
\citation{Verlet1967}
\citation{Fraige2004}
\citation{Swope1982}
\@writefile{brf}{\backcite{Beck2000}{{41}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Wang2003}{{41}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Wei1994}{{41}{2.4.2}{subsection.2.4.2}}}
\@writefile{brf}{\backcite{Light2000}{{41}{2.4.2}{subsection.2.4.2}}}
\newlabel{Newton}{{2.70}{41}{Classical Molecular Dynamics}{equation.2.4.70}{}}
\newlabel{Force}{{2.71}{41}{Classical Molecular Dynamics}{equation.2.4.71}{}}
\@writefile{brf}{\backcite{Batina1991}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Butcher2008}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Gear1966}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Diethelm2002}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Runge1895}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Kutta1901}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Butcher1996}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Butcher2015}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Verlet1967}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Fraige2004}{{42}{2.4.2}{equation.2.4.71}}}
\@writefile{brf}{\backcite{Swope1982}{{42}{2.4.2}{equation.2.4.71}}}
\newlabel{DiscretisationError1}{{2.72}{42}{Classical Molecular Dynamics}{equation.2.4.72}{}}
\newlabel{DiscretisationError2}{{2.73}{42}{Classical Molecular Dynamics}{equation.2.4.73}{}}
\newlabel{Verlet}{{2.74}{42}{Classical Molecular Dynamics}{equation.2.4.74}{}}
\newlabel{Velo}{{2.75}{42}{Classical Molecular Dynamics}{equation.2.4.75}{}}
\newlabel{Newton}{{2.67}{41}{Classical Molecular Dynamics}{equation.2.4.67}{}}
\newlabel{Force}{{2.68}{41}{Classical Molecular Dynamics}{equation.2.4.68}{}}
\@writefile{brf}{\backcite{Batina1991}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Butcher2008}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Gear1966}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Diethelm2002}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Runge1895}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Kutta1901}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Butcher1996}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Butcher2015}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Verlet1967}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Fraige2004}{{41}{2.4.2}{equation.2.4.68}}}
\@writefile{brf}{\backcite{Swope1982}{{41}{2.4.2}{equation.2.4.68}}}
\newlabel{DiscretisationError1}{{2.69}{42}{Classical Molecular Dynamics}{equation.2.4.69}{}}
\newlabel{DiscretisationError2}{{2.70}{42}{Classical Molecular Dynamics}{equation.2.4.70}{}}
\newlabel{Verlet}{{2.71}{42}{Classical Molecular Dynamics}{equation.2.4.71}{}}
\newlabel{Velo}{{2.72}{42}{Classical Molecular Dynamics}{equation.2.4.72}{}}
\newlabel{VV}{{2.74}{42}{Classical Molecular Dynamics}{equation.2.4.74}{}}
\citation{Ray1981,Ray1999}
\citation{Nose1984,Nose1984M,Hoover1985,Hunenberger2005}
\newlabel{VV}{{2.77}{43}{Classical Molecular Dynamics}{equation.2.4.77}{}}
\newlabel{VV1}{{2.78}{43}{Classical Molecular Dynamics}{equation.2.4.78}{}}
\newlabel{VV2}{{2.79}{43}{Classical Molecular Dynamics}{equation.2.4.79}{}}
\newlabel{VV3}{{2.80}{43}{Classical Molecular Dynamics}{equation.2.4.80}{}}
\@writefile{brf}{\backcite{Ray1981}{{43}{2.4.2}{equation.2.4.80}}}
\@writefile{brf}{\backcite{Ray1999}{{43}{2.4.2}{equation.2.4.80}}}
\citation{Torrie1977,Hansmann1993,Marchi1999,Bartels2000,Darve2001}
\citation{Laio2002,Iannuzzi2003,Barducci2011}
\citation{Kirkpatrick1983,Van1987}
\@writefile{brf}{\backcite{Nose1984}{{44}{2.4.2}{equation.2.4.80}}}
\@writefile{brf}{\backcite{Nose1984M}{{44}{2.4.2}{equation.2.4.80}}}
\@writefile{brf}{\backcite{Hoover1985}{{44}{2.4.2}{equation.2.4.80}}}
\@writefile{brf}{\backcite{Hunenberger2005}{{44}{2.4.2}{equation.2.4.80}}}
\newlabel{VV1}{{2.75}{43}{Classical Molecular Dynamics}{equation.2.4.75}{}}
\newlabel{VV2}{{2.76}{43}{Classical Molecular Dynamics}{equation.2.4.76}{}}
\newlabel{VV3}{{2.77}{43}{Classical Molecular Dynamics}{equation.2.4.77}{}}
\@writefile{brf}{\backcite{Ray1981}{{43}{2.4.2}{equation.2.4.77}}}
\@writefile{brf}{\backcite{Ray1999}{{43}{2.4.2}{equation.2.4.77}}}
\@writefile{brf}{\backcite{Nose1984}{{43}{2.4.2}{equation.2.4.77}}}
\@writefile{brf}{\backcite{Nose1984M}{{43}{2.4.2}{equation.2.4.77}}}
\@writefile{brf}{\backcite{Hoover1985}{{43}{2.4.2}{equation.2.4.77}}}
\@writefile{brf}{\backcite{Hunenberger2005}{{43}{2.4.2}{equation.2.4.77}}}
\citation{Swendsen1986}
\citation{Geyer1991}
\citation{Hukushima1996}
\citation{Falcioni1999}
\citation{Earl2005}
\citation{Sugita1999}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.4.3}Parallel-Tempering Molecular Dynamics}{44}{subsection.2.4.3}}
\newlabel{sec:PTMD}{{2.4.3}{44}{Parallel-Tempering Molecular Dynamics}{subsection.2.4.3}{}}
\@writefile{brf}{\backcite{Torrie1977}{{44}{2.4.3}{subsection.2.4.3}}}
@ -373,35 +379,29 @@
\@writefile{brf}{\backcite{Barducci2011}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Kirkpatrick1983}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Van1987}{{44}{2.4.3}{subsection.2.4.3}}}
\citation{Swendsen1986}
\citation{Geyer1991}
\citation{Hukushima1996}
\citation{Falcioni1999}
\citation{Earl2005}
\citation{Sugita1999}
\@writefile{brf}{\backcite{Swendsen1986}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Geyer1991}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Hukushima1996}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Falcioni1999}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Earl2005}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Sugita1999}{{45}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Swendsen1986}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Geyer1991}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Hukushima1996}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Falcioni1999}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Earl2005}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{brf}{\backcite{Sugita1999}{{44}{2.4.3}{subsection.2.4.3}}}
\@writefile{lof}{\contentsline {figure}{\numberline {2.4}{\ignorespaces Schematics of the PTMD algorithm in its synchronous version. Replicas of the same system, numbered from $C_1$ to $C_N$, are simulated subject to different temperatures (from T$_1$ to T$_N$). Once $X$ MD steps (straight solid arrows) have been performed by each replica, configuration exchanges are attempted between neighbouring simulations according to the Metropolis criterion. Some of them undergo successful swapping (solid curved arrows) while other not (dashed curved arrows). MD simulations then proceed for $X$ additional MD steps before new attempts of exchange.\relax }}{45}{figure.caption.8}}
\newlabel{ptmd_s}{{2.4}{45}{Schematics of the PTMD algorithm in its synchronous version. Replicas of the same system, numbered from $C_1$ to $C_N$, are simulated subject to different temperatures (from T$_1$ to T$_N$). Once $X$ MD steps (straight solid arrows) have been performed by each replica, configuration exchanges are attempted between neighbouring simulations according to the Metropolis criterion. Some of them undergo successful swapping (solid curved arrows) while other not (dashed curved arrows). MD simulations then proceed for $X$ additional MD steps before new attempts of exchange.\relax }{figure.caption.8}{}}
\citation{Wales1997,Wales1999}
\citation{Hartke1993,Unger1993,Sivanandam2008,Toledo2014}
\newlabel{Metropolis}{{2.81}{46}{Parallel-Tempering Molecular Dynamics}{equation.2.4.81}{}}
\newlabel{newVelo}{{2.82}{46}{Parallel-Tempering Molecular Dynamics}{equation.2.4.82}{}}
\newlabel{Metropolis}{{2.78}{46}{Parallel-Tempering Molecular Dynamics}{equation.2.4.78}{}}
\newlabel{newVelo}{{2.79}{46}{Parallel-Tempering Molecular Dynamics}{equation.2.4.79}{}}
\@writefile{toc}{\contentsline {subsection}{\numberline {2.4.4}Global Optimization}{46}{subsection.2.4.4}}
\@writefile{brf}{\backcite{Wales1997}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Wales1999}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Hartke1993}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Unger1993}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Sivanandam2008}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Toledo2014}{{47}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Wales1997}{{46}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Wales1999}{{46}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Hartke1993}{{46}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Unger1993}{{46}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Sivanandam2008}{{46}{2.4.4}{subsection.2.4.4}}}
\@writefile{brf}{\backcite{Toledo2014}{{46}{2.4.4}{subsection.2.4.4}}}
\FN@pp@footnotehinttrue
\@setckpt{2_Introduction/introduction}{
\setcounter{page}{48}
\setcounter{equation}{82}
\setcounter{equation}{79}
\setcounter{enumi}{5}
\setcounter{enumii}{0}
\setcounter{enumiii}{0}

View File

@ -482,17 +482,16 @@ Combining eqs \ref{E_KS}, \ref{Tn}, and \ref{Ehartree}, the effective potential
\begin{align}
V_\mathrm{eff}[\rho(\mathbf{r})]&=V_\mathrm{ext}[\rho(\mathbf{r})]+V_\mathrm{H}[\rho(\mathbf{r})]+V_\mathrm{xc}[\rho(\mathbf{r})] \nonumber \\
&=V_\mathrm{ext}[\rho(\mathbf{r})]+\int\frac{\rho(\mathbf{r'})}{|\mathbf{r}-\mathbf {r'}|}d\mathbf{r'} +\frac{\delta E_\mathrm{xc}[\rho(\mathbf{r})]}{\delta \rho(\mathbf{r})} \\
&=\frac{\delta E_\mathrm{eff}[\rho(\mathbf{r})]}{\delta\rho(\mathbf{r})}
&=\frac{\partial E_\mathrm{eff}[\rho(\mathbf{r})]}{\partial\rho(\mathbf{r})}
\label{Veff}
\end{align}
where $V_{H}$[$\rho$(\textbf{r})] refers to the Hartree potential and $V_{xc}$[$\rho$(\textbf{r})] is the exchange-correlation potential.
Therefore, for a given $V_\mathrm{eff}$[$\rho$(\textbf{r})], one can obtain $\rho$(\textbf{r}) by making the right-hand side of equation~\ref{Lagrange}
independent of \textbf{r} which is done by introducing the molecular orbitals $\psi_{i}$(\textbf{r}) such that:
\begin{eqnarray}
\label{rho}
\rho(\mathbf{r})=\sum_{i}^{N} \left |\psi_i(\mathbf{r}) \right |^2
\end{eqnarray}
%Therefore, for a given $V_\mathrm{eff}$[$\rho$(\textbf{r})], one can obtain $\rho$(\textbf{r}) by making the right-hand side of equation~\ref{Lagrange} independent of \textbf{r} which is done by introducing the molecular orbitals $\psi_{i}$(\textbf{r}) such that:
%\begin{eqnarray}
%\label{rho}
%\rho(\mathbf{r})=\sum_{i}^{N} \left |\psi_i(\mathbf{r}) \right |^2
%\end{eqnarray}
%<<<<<<< HEAD
%where \textbf{x} denotes the four vector-containing space and spin variables, and the integration is implemented over the spin variable $\sigma$. Moreover, the molecular orbitals $\psi_{i}$[$\rho$(\textbf{r})] should satisfy the one-electron KS equation,
%\begin{eqnarray}
@ -500,19 +499,18 @@ independent of \textbf{r} which is done by introducing the molecular orbitals $\
%\underbrace{\left(-\frac{1}{2}\mathbf{\nabla}_i^2+V_\mathrm{eff}[\rho(\mathbf{r})]\right)}_\mathrm{KS ~ operator} \psi_i[\rho(\mathbf{r})]=E_i\psi_i[\rho(\mathbf{r})]
%=======
The molecular orbitals $\psi_{i}$(\textbf{r}) should satisfy the one-electron KS equations:
In addition, the molecular orbitals $\phi_{i}$(\textbf{r}) should satisfy the one-electron KS equations:
\begin{eqnarray}
\label{KS}
\underbrace{\left(-\frac{1}{2}\mathbf{\nabla}_i^2+V_\mathrm{eff}[\rho(\mathbf{r})]\right)}_\mathrm{KS ~ operator} \psi_{i}(\textbf{r})=E_i \psi_{i}(\textbf{r})
\underbrace{\left(-\frac{1}{2}\mathbf{\nabla}_i^2+V_\mathrm{eff}[\rho(\mathbf{r})]\right)}_\mathrm{KS ~ operator} \phi_{i}(\textbf{r})=E_i \phi_{i}(\textbf{r})
%>>>>>>> 92023a10c3aa8b7dc4ace43987c1d571fb99a738
\end{eqnarray}
This result can be reobtained within a variational context when looking for those orbitals minimizing the
energy functional of equation~\ref{HK}, subject to orthonormality conditions:
\begin{eqnarray}
\label{orthonormality}
\int \psi^*_i(\mathbf{r}) \psi_j(\mathbf{r}) d\mathbf{r}=\delta_{ij}
\int \phi^*_i(\mathbf{r})\phi_j(\mathbf{r}) d\mathbf{r}=\delta_{ij}
\end{eqnarray}
%<<<<<<< HEAD
@ -570,9 +568,9 @@ Non-self-consistent DFTB scheme is suitable to study systems in which polyatomic
This is the case ighly ionic and homonuclear covalent systems. However, uncertainty rises in the non-self-consistent DFTB scheme whem chemical bonds
are controlled by a subtler charge balance between atoms especially for polar, semi-conductor and heteronuclear molecules. The self-consistent-charge
extension of DFTB, SCC-DFTB, was born as an improvement of standard DFTB to provide a better description of electronic systems in which long-range
Coulomb interactions are significant.\cite{Elstner1998, Porezag1995, Seifert1996, Elstner1998, Elstne2007, Oliveira2009} However, this method is not
fully suitable for the calculation of molecular aggregates because of the dispersion and charge resonance present in such systems. Further corrections
were made to the model in order to address these problems. I know describe in more details how SCC-DFTB can be derived from DFT.
Coulomb interactions are significant.\cite{Elstner1998, Porezag1995, Seifert1996, Elstner1998, Elstne2007, Oliveira2009}
%However, this method is not fully suitable for the calculation of molecular aggregates because of the dispersion and charge resonance present in such systems. Further corrections were made to the model in order to address these problems.
I now describe in more details how SCC-DFTB can be derived from DFT.
\textbf{General principles of DFTB.} DFTB is derived from DFT based on the following approximations:
@ -593,7 +591,9 @@ were made to the model in order to address these problems. I know describe in mo
\textbf{First-order DFTB} (historically referred to as zeroth-order DFTB) takes into account the first term of Taylor's expansion and which is
equivalence with the other tight-binding model. \textbf{Second-order DFTB} (historically the SCC-DFTB) introduces a self-consistent procedure on
atomic charges. There is also a more recent \textbf{third-order extension of DFTB} (referred to as DFTB3).\cite{Gaus2011} DFTB3 was not used in
this thesis so the third-order expansion term will not be shown in the following equations. According to equation~\ref{E_KS}, $E_\mathrm{DFT}[\rho(\mathbf{r})]$
this thesis so the third-order expansion term will not be shown in the following equations.
According to equation~\ref{E_KS}, $E_\mathrm{DFT}[\rho(\mathbf{r})]$
can be written as:
%>>>>>>> 92023a10c3aa8b7dc4ace43987c1d571fb99a738
\begin{align}
@ -620,58 +620,61 @@ E_\mathrm{DFTB}[\rho_{_0}(\mathbf{r})+\delta \rho(\mathbf{r})]=& \overbrace{E_\m
\label{Edftb-long}
\end{align}
From the first line of equation~\ref{Edftb-long}, we can define that the reference Hamiltonian $\hat{H}^0$ only depends on the reference
electron density $\rho_{_0}$:
The right-hand side terms in the first line of equation~\ref{Edftb-long} only depends on with $\rho_0(\mathbf{r})$, and correspond to a repulsive contribution $E_\mathrm{rep}$. The sum of the terms in the second line is the so-called band energy $E_\mathrm{band}$.
We can define that a reference Hamiltonian $\hat{H}^0$, which only depends on the reference electron density $\rho_{_0}$:
\begin{eqnarray}
\label{H0}
\hat{H}^0=-\frac{1}{2}\nabla^2 + \underbrace{V_\mathrm{ext}[\rho(\mathbf{r})]+ \int \frac{\rho_{_0}'(\mathbf{r'})}{|\mathbf{r}-\mathbf{r'}|}d\mathbf{r'} +V_\mathrm{xc}[\rho_{_0}(\mathbf{r})] }_{V_\mathrm{eff}([\rho_{_0}(\mathbf{r})])}
\end{eqnarray}
where we combined the last three terms in the operator as $V_\mathrm{eff}([\rho_{_0}(\mathbf{r})])$.
The right-hand side terms in the first line of equation~\ref{Edftb-long} vary linearly with $\rho_0(\mathbf{r})$, and correspond to a
repulsive contribution $E_\mathrm{rep}$. The sum of the terms in the second line is the so-called band energy $E_\mathrm{band}$. The third
line is the second-order energy $E_\mathrm{2nd}$. Equation~\ref{Edftb-long} can be rewritten as follows:
The third line is the second-order energy $E_\mathrm{2nd}$. Equation~\ref{Edftb-long} can be rewritten as follows:
\begin{align}
E_\mathrm{DFTB}[\rho_{_0}(\mathbf{r})+\delta \rho(\mathbf{r})]=& E_\mathrm{rep}[\rho_{_0}(\mathbf{r})] + \overbrace{\sum_{i}^\mathrm{occ} n_i \left\langle \psi_i \left\vert \hat{H}^0 \right\vert \psi_i \right\rangle}^{E_\mathrm{band}} ~ + E_\mathrm{2nd}\left[\rho_{_0}(\mathbf{r}), (\delta\rho(\mathbf{r}))^2\right]
\label{Eterms}
\end{align}
\textbf{Band energy term.}
In DFTB, one relies on the use of LCAO for the description of the Kohn-Sham molecular orbitals $\psi_i$(\textbf{r}). Here the atomic
In DFTB, it relies on the use of LCAO for the description of the KS molecular orbitals $\phi_i$(\textbf{r}). Here the atomic
orbitals are limited to the valence orbitals of atoms:
\begin{align}
\psi_i(\mathbf{r}) = \sum_\nu C_{i\nu} \varphi_\nu(\mathbf{r}-\mathbf{R_\alpha})
\label{orbitals}
\end{align}
$E_\mathrm{band}$ is the sum over the energies of all occupied orbitals obtained by diagonalization of the parameterized Hamiltonian matrix.
For atoms $\alpha$ and $\beta$, the corresponding atomic orbitals are $\varphi_\mu$ and $\varphi_\nu$.
Then $E_\mathrm{band}$ can be rewritten with equations~\ref{Veff} and \ref{orbitals}:
where $\varphi_\mu$ is the orbital of atom $\alpha$.
%$E_\mathrm{band}$ is the sum over the energies of all occupied orbitals obtained by diagonalization of the parameterized Hamiltonian matrix.
%For atoms $\alpha$ and $\beta$, the corresponding atomic orbitals are $\varphi_\mu$ and $\varphi_\nu$.
$E_\mathrm{band}$ can be rewritten with equations~\ref{Veff} and \ref{orbitals}:
\begin{align}
\sum_{i}^\mathrm{occ} n_i \left\langle \psi_i \left\vert \hat{H}^0 \right\vert \psi_i \right\rangle=\sum_{i}^\mathrm{occ} n_i \sum_{\mu}^\mathrm{occ} \sum_{\nu}^\mathrm{occ} C_{i\mu}C_{i\nu} \underbrace{\left\langle \varphi_\mu \left\vert \hat{T}_e[\rho(\mathbf{r})] + V_\mathrm{eff}[\rho_{_0}(\mathbf{r})] \right\vert \varphi_\nu \right \rangle}_{H^0_{\mu\nu}}
\sum_{i}^\mathrm{occ} n_i \left\langle \psi_i \left\vert \hat{H}^0 \right\vert \psi_i \right\rangle=\sum_{i}^\mathrm{occ} n_i \sum_{\mu}^\mathrm{occ} \sum_{\nu}^\mathrm{occ} C_{i\mu}C_{i\nu} \underbrace{\left\langle \varphi_\mu \left\vert \hat{T}_e[\rho(\mathbf{r})] + V_\mathrm{eff}[\rho_{_0}(\mathbf{r})] \right\vert \varphi_\nu \right \rangle}_{H^0_{\mu\nu}}, \mathrm {with} ~ \mu \in {\alpha}, \nu \in {\beta}
\label{Eband}
\end{align}
$\forall\mu\in {\alpha}, \nu\in{\beta}$. The Hamiltonian matrix element $H^0_{\mu\nu}$ is defined as:
%where $\varphi_\nu$ is the orbital of atom $\beta$.
%$\forall\mu\in {\alpha}, \nu\in{\beta}$.
The Hamiltonian matrix element $H^0_{\mu\nu}$ is defined as:
\begin{align}
H^0_{\mu\nu}&= \left\langle \varphi_\mu \left\vert \hat{H}^0 \right\vert \varphi_\nu \right\rangle, \nonumber \\
\label{matrix}
\end{align}
The effective potential $V_\mathrm{eff}[\rho_{_0}(\mathbf{r})]$ is defined as the sum of potentials $V_\alpha(\textbf{r})$ centered on the atoms:
\begin{align}
V_\mathrm{eff}[\rho_{_0}(\mathbf{r})]=\sum_{\alpha}V
_\alpha(\mathbf{r}_\alpha)
_\alpha(\mathbf{r}-\mathbf{R}_\alpha)
\label{Vdftb}
\end{align}
where $\mathbf{r}_\alpha=\mathbf{r}-\mathbf{R}_\alpha$.
%where $\mathbf{r}-\mathbf{R}_\alpha$ equals to $\mathbf{r}_\alpha$.
The Hamiltonian matrix elements can be written as follows:
%\begin{align}
%H^0_{\mu\nu}=\left\langle\varphi_\mu \left\vert -\frac{1}{2}\nabla_\nu^2 + V_\alpha+(1-\delta_{\alpha\beta})V_\beta \right\vert \varphi_\nu \right\rangle, \mathrm {with} ~ \mu \in {\alpha}, \nu \in {\beta}
%\label{Hamil}
%\end{align}
%where $\delta_{\alpha\beta}$ is the Kronecker's delta.
\begin{align}
H^0_{\mu\nu}=\left\langle\varphi_\mu \left\vert -\frac{1}{2}\nabla_\nu^2 + V_\alpha+(1-\delta_{\alpha\beta})V_\beta \right\vert \varphi_\nu \right\rangle, \mathrm {with} ~ \mu \in {\alpha}, \nu \in {\beta}
H^0_{\mu\nu}=\left\langle\varphi_\mu \left\vert -\frac{1}{2}\nabla_i^2 + V_\alpha+ V_\beta\right\vert \varphi_\nu \right\rangle
\label{Hamil}
\end{align}
where $\delta_{\alpha\beta}$ is the Kronecker's delta.
%In practice, $H^0_{\mu\nu}$ are calculated as follows:
For diagonal elements, the energy level in the free atom is chosen, which ensures the correct dissociation limits.
@ -682,7 +685,7 @@ To summarize, within the electronic density superposition approach, the $H^0_{\m
H^0_{\mu\nu}=
\begin{cases}
\varepsilon_\mu^\mathrm{free~atom}, ~ & \mu=\nu \\
\left \langle \varphi_\mu \left\vert -\frac{1}{2}\nabla_\nu^2 + V_\alpha+(1-\delta_{\alpha\beta})V_\beta \right\vert \varphi_\nu \right\rangle, ~ & \mu \in {\alpha}, \nu \in {\beta}, \alpha \neq \beta \\
\left \langle \varphi_\mu \left\vert -\frac{1}{2}\nabla_i^2 + V_\alpha+V_\beta \right\vert \varphi_\nu \right\rangle, ~ & \mu \in {\alpha}, \nu \in {\beta}, \alpha \neq \beta \\
0, ~ & \mathrm{otherwise}
\end{cases}
\label{Hmatrix}
@ -691,9 +694,6 @@ H^0_{\mu\nu}=
It should be noted that the $H^0_{\mu\nu}$ elements only depend on atoms $\alpha$ and $\beta$. Therefore only the two-center matrix elements and
the two-center elements of the overlap matrix can be explicitly calculated, in other words, interactions at three or more centers are neglected as stated above.
\textbf{Repulsive energy term.} $E_\mathrm{rep}$ is a repulsive contribution obtained from the sum of atomic-pair terms, which only depend
on the reference electronic density $\rho_{_0}(\textbf{r})$.
%In order to obtain a good evaluation of $\rho_{_0}(\textbf{r})$, it is approximated as a overlap of atom-like densities centered on the nuclei $\alpha$:
%\begin{align}
%\rho_{_0}(\mathbf{r})=\sum_{\alpha}^{N}
@ -712,29 +712,6 @@ on the reference electronic density $\rho_{_0}(\textbf{r})$.
%\end{align}
%where $E_\mathrm{long}$ refers to the interactions energy of more than two centers.
Due to the screening of terms of more than two centers, one can assume that the two-center contributions are short ranged.
But $E_\mathrm{rep}$ doesn't decay to zero for long interatomic distances. Instead, it decays to a constant given by the atomic
contributions:
\begin{align}
\displaystyle \lim_{R_{\alpha\beta} \to \infty}E_\mathrm{rep}[\rho_{_0}(\mathbf{r})]=\sum_{\alpha}^{N}
E_\mathrm{rep}[\rho_{_0}^\alpha(\mathbf{r}_\alpha)]
\label{Erep-infinity}
\end{align}
The right side of equation~\ref{Erep-infinity} is assumed to make $E_\mathrm{rep}$ only rely on the two-center contributions:
\begin{align}
E_\mathrm{rep}[\rho_{_0}(\mathbf{r})]\approx\frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} V(\mathbf{R}_\alpha-\mathbf{R}_\beta)
\label{Erep-right}
\end{align}
In practice, it it possible to calculate $E_\mathrm{rep}$ with known values of $\rho_{_0}(\textbf{r})$, but it's more convenient to adjust
the expression of $E_\mathrm{rep}$ to $ab ~ initio$ calculations. Therefore, $E_\mathrm{rep}$ is determined by comparing the
difference between the DFT energy $E_\mathrm{DFT}$ and $E_\mathrm{band}$+$E_\mathrm{2nd}$ as a function of the interatomic distance $R_{\alpha\beta}$:
\begin{align}
E_\mathrm{rep}[\rho_{_0}(\mathbf{r})] \equiv E_\mathrm{rep}(R_{\alpha\beta})=E_\mathrm{DFT}(R_{\alpha\beta})-E_\mathrm{band}(R_{\alpha\beta})-E_\mathrm{2nd}(R_{\alpha\beta})
\label{Erep-final}
\end{align}
\textbf{Second-order term.}
In SCC-DFTB, the electronic density is corrected by including the second-order contribution $E_\mathrm{2nd}$ in equation~\ref{Eterms},
which is ignored in first-order DFTB.
@ -761,20 +738,21 @@ $E_\mathrm{2nd}$ can be rewritten with equations~\ref{Edftb-long} and \ref{Eterm
\begin{align}
E_\mathrm{2nd}\approx
&\frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} \Delta q_{\alpha} \Delta q_{\beta} \overbrace{\iint \left(\frac{1}{|\mathbf{r}-\mathbf {r'}|} + \frac{\delta^2 E_{xc}[\rho_{_0}(\mathbf{r})]}{\delta\rho(\mathbf{r})\delta\rho(\mathbf{r'})} \Big\vert _{\rho_{_0}, ~ \rho_{_0}'} \right) F(\alpha, \beta) d\mathbf{r}d\mathbf{r'}}^{\gamma_{\alpha\beta}} \nonumber \\
=& \frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} \Delta q_{\alpha} \Delta q_{\beta} \gamma_{\alpha\beta} \\
F(\alpha, \beta)= & F_{0}(\mathbf{r}-\mathbf{R}_\alpha)\times F_{0}(\mathbf{r'}-\mathbf{R}_\beta)
=& \frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} \Delta q_{\alpha} \Delta q_{\beta} \gamma_{\alpha\beta}
\label{E2nd}
\end{align}
where the two-electron integrals $\gamma_{\alpha\beta}$ is introduced for convenience.
where $F(\alpha, \beta)= F_{0}(\mathbf{r}-\mathbf{R}_\alpha)\times F_{0}(\mathbf{r'}-\mathbf{R}_\beta) \nonumber$ the two-electron integrals $\gamma_{\alpha\beta}$ is introduced for convenience.
To calculate equation~\ref{E2nd}, $\gamma_{\alpha\beta}$ must be analyzed. In the limit case, the interatomic distance is very large,
$|{\textbf{R}_\alpha - \textbf{R}_\beta}|$ = $|\textbf{r} - \textbf{r}'| \rightarrow \infty$ with GGA-DFT, the exchange-correlation term tends to zero.
$\gamma_{\alpha\beta}$ that describes the interaction of two normalized spherical electronic densities reduces to 1/$|\textbf{R}_\alpha-\textbf{R}_\beta|$, so $E_\mathrm{2nd}$ can be expressed as follows:
\begin{align}
E_\mathrm{2nd}=\frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} \frac{\Delta q_{\alpha} \Delta q_{\beta}} {|\mathbf{R}_\alpha-\mathbf{R}_\beta|}
\label{E2ndsimple}
E_\mathrm{2nd}=\frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} \frac{\Delta q_{\alpha} \Delta q_{\beta}} {|\mathbf{R}_\alpha-\mathbf{R}_\beta|} + f(U_\alpha, U_\beta, (\mathbf{R}_\alpha-\mathbf{R}_\beta))
\label{simpleE2nd}
\end{align}
where U$_\alpha$, U$_\beta$ are the Hubbard parameters: $\gamma_{\alpha\alpha} \approx I_\alpha-A_\alpha \approx 2\eta\alpha \approx U_\alpha$.
$I$ is the ionisation potential and $A$ is the electonic affinity of atom $\alpha$. $\eta_\alpha$ refers to the chemical hardness of atom $\alpha$.\cite{Elstner1998}
It is worth noting that the electronic density $\rho(\mathbf{r})$ influences explicitly the calculation of the electrostatic energy in DFT.
In the context of DFTB, point charges are used and the electronic density around the atom is condensed at a point. In practice, Mulliken's
@ -785,7 +763,30 @@ q_\alpha^\mathrm{Mull}=\frac{1}{2}\sum_{\alpha}^\mathrm{occ} n_i
\right)
\label{Mulliken}
\end{align}
where $\mu$ denotes the orbitals belonging to the atom $\alpha$. This definition does not allow the bond between two different atoms to polarize.
%where $\mu$ denotes the orbitals belonging to the atom $\alpha$. This definition does not allow the bond between two different atoms to polarize.
\textbf{Repulsive energy term.} $E_\mathrm{rep}$ is a repulsive contribution obtained from the sum of atomic-pair terms, which only depend
on the reference electronic density $\rho_{_0}(\textbf{r})$.
%Due to the screening of terms of more than two centers, one can assume that the two-center contributions are short ranged. But $E_\mathrm{rep}$ doesn't decay to zero for long interatomic distances. Instead, it decays to a constant given by the atomic contributions:
%\begin{align}
%\displaystyle \lim_{R_{\alpha\beta} \to \infty}E_\mathrm{rep}[\rho_{_0}(\mathbf{r})]=\sum_{\alpha}^{N}
%E_\mathrm{rep}[\rho_{_0}^\alpha(\mathbf{r}_\alpha)]
%\label{Erep-infinity}
%\end{align}
%The right side of equation~\ref{Erep-infinity}
It is assumed to make $E_\mathrm{rep}$ only rely on the two-center contributions:
\begin{align}
E_\mathrm{rep}[\rho_{_0}(\mathbf{r})]\approx\frac{1}{2}\sum_{\alpha}^{N}\sum_{\beta}^{N} V(\mathbf{R}_\alpha-\mathbf{R}_\beta)
\label{Erep-right}
\end{align}
In practice, it it possible to calculate $E_\mathrm{rep}$ with known values of $\rho_{_0}(\textbf{r})$, but it's more convenient to adjust
the expression of $E_\mathrm{rep}$ to $ab ~ initio$ calculations. Therefore, $E_\mathrm{rep}$ is determined by comparing the
difference between the DFT energy $E_\mathrm{DFT}$ and $E_\mathrm{band}$+$E_\mathrm{2nd}$ as a function of the interatomic distance $R_{\alpha\beta}$:
\begin{align}
E_\mathrm{rep}[\rho_{_0}(\mathbf{r})] \equiv E_\mathrm{rep}(R_{\alpha\beta})=E_\mathrm{DFT}(R_{\alpha\beta})-E_\mathrm{band}(R_{\alpha\beta})-E_\mathrm{2nd}(R_{\alpha\beta})
\label{Erep-final}
\end{align}
\textbf{Total energy.}
The total energy in SCC-DFTB can be written from the previous different contributions as follows:
@ -802,24 +803,26 @@ V(\mathbf{R}_\alpha-\mathbf{R}_\beta) \nonumber \\
%\textbf{Secular equation}
%=======
Here the contribution from $H^0$ is exactly the same with the one in standard DFTB scheme and the molecular orbitals $\psi_i$ are expanded in a LCAO model (equation~\ref{orbitals}).
%Here the contribution from $H^0$ is exactly the same with the one in standard DFTB scheme and the molecular orbitals $\psi_i$ are expanded in a LCAO model (equation~\ref{orbitals}).
%>>>>>>> 92023a10c3aa8b7dc4ace43987c1d571fb99a738
\textbf{Secular equations.}
From this LCAO model, we can get the secular equations:
From this energy expression model, we can get the secular equations:
\begin{align}
\sum_\nu C_{i\nu} \left({H}^0_{\mu\nu} - \varepsilon_i S_{\mu\nu} \right) =0, ~ \forall \mu, ~ \nu
\label{secular-eq}
\end{align}
where $H^0_{\mu\nu}$ are the Hamiltonian matrix elements and $S_{\mu\nu}$ are the overlap matrix elements.
where $H^0_{\mu\nu}$ is the DFTB operator matrix element and $S_{\mu\nu}$ are the overlap matrix elements.
The secular equations can be rewritten with modified Hamiltonian matrix elements, $H_{\mu\nu}$, defined by:
The secular equations can be rewritten using the matrix, $H_{\mu\nu}$, defined by:
\begin{align}
H_{\mu\nu}=H_{\mu\nu}^0 + \frac{1}{2} S_{\mu\nu} \sum_{\zeta}(\gamma_{\alpha\zeta} + \gamma_{\beta\zeta})\Delta q_\zeta
H_{\mu\nu}=H_{\mu\nu}^0 + \underbrace{\frac{1}{2} S_{\mu\nu} \sum_{\zeta}(\gamma_{\alpha\zeta} + \gamma_{\beta\zeta})\Delta q_\zeta}_{H_{\mu\nu}^\prime}
\label{Huv}
\end{align}
The $H_{\mu\nu}^0$ and $S_{\mu\nu}$ matrix elements are the same than the ones defined in the standard DFTB method (equation~\ref{matrix}).
$H_{\mu\nu}^0$ comes from the band energy and $H_{\mu\nu}^\prime$ is from the second order term.
%The $H_{\mu\nu}^0$ and $S_{\mu\nu}$ matrix elements are the same than the ones defined in the standard DFTB method (equation~\ref{matrix}).
The $H_{\mu\nu}$ elements rely on the atomic charges explicitly, and the atomic charges depend on the molecular orbitals (see equation~\ref{Mulliken}).
Then the resolution can be achieved in a self-consistent way.
First, from an initial set of charges
@ -843,20 +846,20 @@ charges were proposed by J. Li $et al.$ in 1998 and they have been considered in
and the electrostatic potential, partial atomic charges in molecules, and Coulombic intermolecular potential of polycyclic aromatic hydrocarbon clusters.\cite{Li1998, Kalinowski2004, Kelly2005, Rapacioli2009corr}
The CM3 charges are defined as:
\begin{align}
q_{\alpha}^\mathrm{CM3} = q_{\alpha}^\mathrm{Mull} + \sum_{\beta \neq \alpha} (C_{t_\alpha t_\beta} K_{\alpha\beta} + D_{t_\alpha t_\beta} K_{\alpha\beta}^2 )
q_{\alpha}^\mathrm{CM3} = q_{\alpha}^\mathrm{Mull} + \sum_{\beta \neq \alpha} (D_{t_\alpha t_\beta} K_{\alpha\beta} + D_{t_\alpha t_\beta} K_{\alpha\beta}^2 )
\label{CM3}
\end{align}
where $K_{\alpha\beta}$ is the Mayer bond order,\cite{Mayer1983, Cances1997, Lu2012} which depends on the density matrix of the orbitals belonging to
each of the two atoms $\alpha$ and $\beta$. $C_{t_\alpha t\beta}$ and $D_{t_\alpha t\beta}$ are empirical parameters which are related
to the nature of atoms $\alpha$ and $\beta$.
where $K_{\alpha\beta}$ is the Mayer bond order,\cite{Mayer1983, Cances1997, Lu2012} between
%which depends on the density matrix of the orbitals belonging to each of the two
atoms $\alpha$ and $\beta$. $C_{t_\alpha t\beta}$ and $D_{t_\alpha t\beta}$ are empirical parameters which are related to the nature of
atoms $\alpha$ and $\beta$.
In this thesis, in practice for the calculation of electronic energy, the definition of CM3 charges is simplified as:
\begin{align}
q_{\alpha}^\mathrm{CM3} = q_{\alpha}^\mathrm{Mull} + \sum_{\beta \neq \alpha} C_{t_\alpha t_\beta} K_{\alpha\beta}
q_{\alpha}^\mathrm{CM3} = q_{\alpha}^\mathrm{Mull} + \sum_{\beta \neq \alpha} D_{t_\alpha t_\beta} K_{\alpha\beta}
\label{CM3-classical}
\end{align}
%<<<<<<< HEAD
%\textbf{Dispersion energy}
%
@ -987,15 +990,12 @@ to understand the PES. Stable molecular structures correspond to the minima in t
of the valley around a minimum determines the vibrational spectrum.
The key points on a PES can be classified according to the first and all second derivatives of the energy with
respect nuclear coordinates, which correspond to the gradients and the curvatures, respectively.
When points have a zero gradient (stationary points) and their second derivatives are positive, this corresponds
to a \textbf{local minimum} (physically stable structure). Among these local minima, the lowest energy minima is called
the \textbf{global minimum}. When at least one of the second derivatives is negative, the point is a transition state (saddle point).
The points have a zero gradient (stationary points) and their second derivatives are positive, they correspond
to \textbf{local minima} (physically stable structure). Among these local minima, the lowest energy minimum is called the \textbf{global minimum}. When at least one of the second derivatives is negative, the point is a transition state (saddle point).
\figuremacro{PES}{Schematic representation of some key points on a model potential energy surface.}
Molecular structure, properties, chemical reactivity, dynamics, and vibrational spectra of molecules can be readily
understood in terms of PES. Only very simple PES can be obtained from experiment whereas computational chemistry
has developed different kinds of methods to efficiently explore PES. To survey the PES, the choice of an exploration
Molecular structure, properties, chemical reactivity, dynamics, and vibrational spectra of molecules can be readily understood in terms of PES. Only very simple PES can be obtained from experiment whereas computational chemistry has developed different kinds of methods to efficiently explore PES. To survey the PES, the choice of an exploration
method can be guided by the shape of the PES (the statistical set that one wishes to study) and the temporal aspect.
\textbf{Monte Carlo} and classical \textbf{MD} simulations are widely recognized approaches for the exploration of PES of systems containing a large number of degrees of freedom such as molecular aggregates.
@ -1019,60 +1019,13 @@ inherent randomness, and the computing power of the computer can directly simula
in nuclear physics, analysis of the transmission process of neutrons in a reactor.
The second category applies to problems that can be transformed into randomly distributed characteristic numbers, such as the
probability of a random event. Through the random sampling method, the probability of random events is estimated by the
frequency of occurrence, and this is used as the solution to the problem. This method is usually used to solve complicated
multi-dimensional integration problems.
frequency of occurrence, and this is used as the solution to the problem. This method is usually used to solve complicated multi-dimensional integration problems. Monte Carlo methods are mainly used in the field of optimization, numerical integration, and generating draws from a probability distribution.
Monte Carlo methods are mainly used in the field of optimization, numerical integration, and generating draws from a probability
distribution. In physics related problems, Monte Carlo methods are applied to simulat systems diplaying many coupled degrees
of freedom, such as fluids, strongly coupled solids, disordered materials, and cellular structures. In statistical physics not
related to thermodynamics, Monte Carlo molecular simulation is an alternative to MD simulations. Monte Carlo methods
can be used to compute statistical field theories of simple particles and polymer systems.\cite{Rosenbluth1955, Binder1993, Baeurle2009}
When solving practical problems using a Monte Carlo methods, it includes two main parts. First, random variables with various
probability distributions need to be generated to simulate a certain process. Second, statistical methods need to be used to
estimate the numerical characteristics of the model to obtain a numerical solutions to the actual problem. Monte Carlo methods
make it possible for the sampling of a PES by performing random shifts in order to correctly reproduce the probability distribution
of the configurations. One example for the calculation procedure of a molecular simulation is as follows:
\begin{enumerate}
\item A random molecular configuration is generated using a random number generator.
\item Random changes are made to the particle coordinates of this molecular configuration, resulting in a new molecular configuration.
\item Calculate the energy of the new molecular configuration.
\item Compare the energy change between the new molecular configuration and the former one to determine whether to accept the configuration change or not.
%>>>>>>> 92023a10c3aa8b7dc4ace43987c1d571fb99a738
\begin{itemize}
\item[$\bullet$] If the energy of the new molecular configuration is lower than that of the original one, the new configuration change is accepted, and the
new configuration is used for the next iteration.
\item[$\bullet$] If not, the Boltzmann factor is calculated and a random number is generated. If this random number is greater than the calculated
Boltzmann factor, then the new configuration is discarded. If not, the new configuration is conserved and is used for the next iteration.
\end{itemize}
\item If a new iteration is required, the process is repeated from step 2.
\end{enumerate}
\textbf{General principles.}
Monte Carlo calculations that lead to quantitative results may be regarded as attempts to estimate the value of a multiple integral.
This is particularly true for the applications in equilibrium statistical thermodynamics, where one hopes to calculate the thermal
average $\langle A \rangle _T$ of an observable $A(\textbf X)$ as an integral over phase space $\Omega$, where $\textbf X$ is a
point in $\Omega$:
\begin{align}
\langle A \rangle _T= \frac{1}{Z} \int_\Omega d\mathbf X A(\mathbf X)e^{-H(\mathbf X)/k_BT}
\label{thermalAverage}
\end{align}
in which $Z$ is the partition function, and $k_B$ is the Boltzmann constant. $T$ refers to the temperature, and $H(\textbf X)$ denotes the Hamiltonian of the system.
To illustrate the general application of Monte Carlo techniques, here the standard example of the one-dimensional integral $I$ over integration space $\Omega$ is taken:
\begin{align}
I=\int_\Omega f(\mathbf x)d\mathbf x
\label{SingleIntegral}
\end{align}
Such integral can be evaluated more efficiently using conventional numerical means than using Monte Carlo methods. However, the extension to higher
dimensions (number of degrees of freedom greater than 3) is always difficult in practice with conventional numerical integration. It is then possible to
use stochastic approaches that one explores the configuration space randomly and the computation of the integral is estimated in the form of an average
@ -1132,6 +1085,51 @@ In the case of symmetrical movements, $\eta (\mathbf x_i \rightarrow \mathbf x_j
\label{assess}
\end{align}
In physics related problems, Monte Carlo methods are applied to simulat systems diplaying many coupled degrees
of freedom, such as fluids, strongly coupled solids, disordered materials, and cellular structures. In statistical physics not
related to thermodynamics, Monte Carlo molecular simulation is an alternative to MD simulations. Monte Carlo methods
can be used to compute statistical field theories of simple particles and polymer systems.\cite{Rosenbluth1955, Binder1993, Baeurle2009}
When solving practical problems using a Monte Carlo methods, it includes two main parts. First, random variables with various probability distributions need to be generated to simulate a certain process. Second, statistical methods need to be used to
estimate the numerical characteristics of the model to obtain a numerical solutions to the actual problem. Monte Carlo methods
make it possible for the sampling of a PES by performing random shifts in order to correctly reproduce the probability distribution
of the configurations.
\textbf{General principles.}
Monte Carlo calculations that lead to quantitative results may be regarded as attempts to estimate the value of a multiple integral.
This is particularly true for the applications in equilibrium statistical thermodynamics, where one hopes to calculate the thermal
average $\langle A \rangle _T$ of an observable $A(\textbf X)$ as an integral over phase space $\Omega$, where $\textbf X$ is a
point in $\Omega$:
\begin{align}
\langle A \rangle _T= \frac{1}{Z} \int_\Omega d\mathbf X A(\mathbf X)e^{-H(\mathbf X)/k_BT}
\label{thermalAverage}
\end{align}
in which $Z$ is the partition function, and $k_B$ is the Boltzmann constant. $T$ refers to the temperature, and $H(\textbf X)$ denotes the Hamiltonian of the system.
5
One example for the calculation procedure of a molecular simulation is as follows:
\begin{enumerate}
\item A random molecular configuration is generated using a random number generator.
\item Random changes are made to the particle coordinates of this molecular configuration, resulting in a new molecular configuration.
\item Calculate the energy of the new molecular configuration.
\item Compare the energy change between the new molecular configuration and the former one to determine whether to accept the configuration change or not.
%>>>>>>> 92023a10c3aa8b7dc4ace43987c1d571fb99a738
\begin{itemize}
\item[$\bullet$] If the energy of the new molecular configuration is lower than that of the original one, the new configuration change is accepted, and the
new configuration is used for the next iteration.
\item[$\bullet$] If not, the Boltzmann factor is calculated and a random number is generated. If this random number is greater than the calculated
Boltzmann factor, then the new configuration is discarded. If not, the new configuration is conserved and is used for the next iteration.
\end{itemize}
\item If a new iteration is required, the process is repeated from step 2.
\end{enumerate}
\subsection{Classical Molecular Dynamics}
%<<<<<<< HEAD
@ -1167,13 +1165,12 @@ energies were studied using MD in 1960 and MD was also applied to simulate liqui
the description of phase equilibria, the relaxation of metastable states and the dynamics of processes
at interfaces. In addition, MD can also model chemical reactions in complex environments.
MD simulations allow to model real time evolution of particles, one can then access time-dependent properties. If the
As MD simulations allow to model real time evolution of particles, one can then access time-dependent properties. If the
time evolution is obtained by \textbf{integrating Newton's equations of motion} for a system of interacting particles, it is
referred to as \textbf{classical MD}.
A classical MD simulation needs the definition of a potential describing the interaction between particles in order to calculate
the PES. %This is often referred to as a force field in chemistry and biology, and as an interatomic potential in materials physics.
Potentials can be of different levels of accuracy as described in the previous section. The most commonly used potentials
in chemistry are \textbf{force fields}. In that case, one refers to \textbf{molecular mechanics}, which embodies a classical
Potentials can be of different levels of accuracy as described in the previous section. The most commonly used potentials in chemistry are \textbf{force fields}. In that case, one refers to \textbf{molecular mechanics}, which embodies a classical
mechanics treatment of the interactions between particles. In classical molecular dynamics, electrons and nuclei are not
distinguished and one refers only to particles. As already mentioned, the main drawback of force fields is that
they usually can not model chemical reactions. If the potential comes from a quantum chemical treatment of the
@ -1208,7 +1205,7 @@ position ($\textbf R_\alpha$) of $\alpha$:
\label{Force}
\end{align}
where $V$ depends on the positions of all atoms or particles. This leads to a system of $f \times N$ second-order differential equations
where $f$ is the dimension of space. In our case, the degree of freedom $f$ is equal to 3. According to the known initial positions of
where $f$ is the dimension of space. In our case, the number of degrees of freedom $f$ is equal to 3. According to the known initial positions of
particles, the potential energy can be obtained.
Then, a numerical resolution of the partial derivative equations provided in equation~\ref{Force} can be obtained using a suitable \textbf{integration algorithm}.
The integration algorithm gives access to the positions and velocities of atoms or particles and to the forces acting on these atoms
@ -1219,11 +1216,12 @@ Many high order integration algorithms have been proposed depending on the desir
Euler algorithm,\cite{Batina1991, Butcher2008} Verlet algorithms, predictor-corrector algorithm,\cite{Gear1966, Diethelm2002} and
Runge-Kutta algorithm.\cite{Runge1895, Kutta1901, Butcher1996, Butcher2015}
The Verlet algorithms include the Simple Verlet (SV),\cite{Verlet1967} the Leapfrog Verlet (LFV),\cite{Fraige2004} and the Velocity Verlet (VV).\cite{Swope1982}
The \textbf{Velocity Verlet algorithm} is the most widely used in many MD codes owing to its numerical stability and implementation simplicity.
Furthermore, movement constants are very well preserved over time. We briefly describe the basis of the Velocity Verlet algorithm below.
The \textbf{VV algorithm} is the most widely used in many MD codes owing to its numerical stability and simplicity implementation.
Furthermore, movement constants are very well preserved over time. We briefly describe the basis of the VV algorithm below.
%The time symmetry in the Verlet algorithm reduces the level of local errors introduced into the integration by the discretization through removing all the odd-degree terms which refers to the three-degree terms in $\delta t$.
The local error is quantified by inserting the exact values $\textbf R_\alpha(t_{n-1})$, $\textbf R_\alpha(t_n)$, and $\textbf R_\alpha(t_{n+1})$ into the iteration and calculating the Taylor expansions at time $t=t_n$ of the position vector $\textbf R_\alpha(t\pm\delta t)$ in different time directions:
The local error is quantified by inserting the exact values $\textbf R_\alpha(t_{n-1})$, $\textbf R_\alpha(t_n)$, and $\textbf R_\alpha(t_{n+1})$ into the iteration and calculating the Taylor expansions at time $t=t_n$ of the position vector $\textbf R_\alpha(t\pm\delta t)$
%in different time directions:
\begin{align}
\mathbf R_\alpha(t+\delta t) = \mathbf R_\alpha(t) + \mathbf v_\alpha(t)\delta t +\frac{\mathbf a_\alpha(t)\delta t^2}{2}+ \frac{\mathbf b_\alpha(t)\delta t^3}{6} + O(\delta t^4)
@ -1265,24 +1263,23 @@ The VV algorithm is often applied to solve this problem as it allows velocities
\end{align}
%The error in the VV is of the same order with the one in the basic Verlet. However, compared to the basic Verlet, the VV algorithm requires less memory because two vectors of position are tracked in basic Verlet while one position vector and one velocity vector are tracked in VV.
The standard implementation of the VV algorithm is a four steps scheme: firstly to calculate equation~\ref{VV1}, secondly to calculate equation~\ref{VV2}, thirdly to
derive $\mathbf a_\alpha(t+ \delta t)$ from the interaction potential with $\mathbf R_\alpha(t+ \delta t)$, finally to compute equation~\ref{VV3}.
The standard implementation of the VV algorithm is a four steps scheme: firstly to calculate the following equation,
\begin{align}
\mathbf v_\alpha(t+\frac{1}{2}\delta t) = \mathbf v_\alpha(t) + \frac{1}{2}\mathbf a_\alpha(t)\delta t
\label{VV1}
\end{align}
Secondly to calculate equation as follows:
\begin{align}
\quad \qquad \mathbf R_\alpha(t+ \delta t) =\mathbf R_\alpha(t) + \mathbf v_\alpha(t + \frac{1}{2} \delta t) \delta t
\label{VV2}
\end{align}
Thirdly to derive $\mathbf a_\alpha(t+ \delta t)$ from the interaction potential at $\mathbf R_\alpha(t+ \delta t)$, finally to compute the following equation.
\begin{align}
\qquad \qquad \qquad \mathbf v_\alpha(t+ \delta t) = \mathbf v_\alpha(t + \frac{1}{2} \delta t) + \frac{1}{2} \mathbf a_\alpha(t+ \delta t) \delta t
\label{VV3}
\end{align}
<%<<<<<< HEAD
%<%<<<<<< HEAD
%The VV algorithm needs to ensure two intrinsic properties of classical motion equations.
%One is temporal reversibility, the invariance of the trajectories at $t$ and $-t$. This symmetry leads to the independence of the dynamics of the microscopic system from the direction of time. The other one is the conservation of the Hamiltonian function over time. Because of the discretization of the trajectories, this conservation can not be assured. So a stable integration algorithm must impose this conservation for enough long time steps ($\delta t$) to allow the long simulation times. The VV algorithm is able to do this due to its sufficient numerical stability. At the end of each integration step, VV algorithm allows to access to $\mathbf R_\alpha (t+\delta t)$, $\mathbf v_\alpha (t+\delta t)$, $\mathbf F_\alpha (t+\delta t)$ directly.
@ -1303,9 +1300,11 @@ Temperature is not a state variable in the simulations, an average kinetic tempe
At the end of each integration step, the VV algorithm gives directly access to $\mathbf R_\alpha (t+\delta t)$, $\mathbf v_\alpha (t+\delta t)$,
and $\mathbf F_\alpha (t+\delta t)$.
The VV algorithm needs to ensure two intrinsic properties of the classical equations of motion. One is temporal reversibility,
The VV algorithm ensures two intrinsic properties of the classical equations of motion. One is temporal reversibility,
the invariance of the trajectories at $t$ and $-t$. This symmetry leads to the independence of
he dynamics from the direction of time. The other property is the conservation of the Hamiltonian function over time.
he dynamics from the direction of time. The other property is the conservation of the
%Hamiltonian function
total energy over time.
Because of the discretization of trajectories, this conservation can not be insured. A stable integration algorithm must
impose this conservation for long enough time steps ($\delta t$) to allow for sufficiently long simulation times. The VV algorithm
is able to do this due to its sufficient numerical stability.
@ -1322,10 +1321,7 @@ forces. \textbf{In this thesis, the SCC-DFTB method has been applied to compute
\subsection{Parallel-Tempering Molecular Dynamics} \label{sec:PTMD}
In a number of cases, it is necessary to explore the PES as thoroughly as possible for the study of dynamical,
thermodynamical, and structural properties of a given system. For atomic and molecular clusters, PES usually
presents a large number of stable configurations linked together by energy barriers. Unfortunately, MD simulations
cannot overcome these energy barriers in a reasonable simulation time even within the canonical ensemble ($N, V, T$).
In a number of cases, it is necessary to explore the PES as thoroughly as possible for the study of dynamical, thermodynamical, and structural properties of a given system. For atomic and molecular clusters, PES usually presents a large number of stable configurations linked together by energy barriers. Unfortunately, MD simulations cannot overcome these energy barriers in a reasonable simulation time even within the canonical ensemble ($N, V, T$).
This leads to non-ergodic simulations that cannot be used to extract meaningful statistical averages. Actually, when
a simulation explores the well of a PES, it may often be blocked in this well at low temperature because the energy
barriers are too high to be crossed. In this case, if $E_b$ refers to the energy barrier and $T$ is the temperature of

File diff suppressed because it is too large Load Diff

View File

@ -3,7 +3,7 @@
\FN@pp@footnotehinttrue
\citation{Brechignac1989,Brechignac1994}
\citation{Wong2004,Bush2008}
\citation{Holm2010,M. Gatchell2014,Gatchell2017}
\citation{Holm2010,Gatchell2014,Gatchell2017}
\citation{Boering1992,Wells2005,Zamith2019thermal}
\@writefile{toc}{\contentsline {chapter}{\numberline {4}Dynamical Simulation of Collision-Induced Dissociation}{97}{chapter.4}}
\@writefile{lof}{\addvspace {10\p@ }}
@ -33,6 +33,7 @@
\citation{Carl2013,Hofstetter2013,Coates2018}
\@writefile{brf}{\backcite{Holm2010}{{98}{4.1}{section.4.1}}}
\@writefile{brf}{\backcite{Gatchell2017}{{98}{4.1}{section.4.1}}}
\@writefile{brf}{\backcite{Gatchell2014}{{98}{4.1}{section.4.1}}}
\@writefile{brf}{\backcite{Zamith2019thermal}{{98}{4.1}{section.4.1}}}
\@writefile{brf}{\backcite{Boering1992}{{98}{4.1}{section.4.1}}}
\@writefile{brf}{\backcite{Wells2005}{{98}{4.1}{section.4.1}}}

View File

@ -30,7 +30,7 @@ The stability of cluster can be investigated from dissociation experiments. Clus
For instance, the sodium cluster ions and lithium cluster cation were dissociated with a pulsed UV laser source.\cite{Brechignac1989, Brechignac1994}
Gaseous hydrated trivalent metal ions were dissociated using blackbody infrared radiative dissociation (BIRD).\cite{Wong2004, Bush2008}
The collision between cluster and high or low energetic particles at different pressure also have been investigated.
Collisions between the high energetic projectile ions (such as 3 keV Ar$^+$, 22.5 keV He$^{2+}$) and neutral targets were investigated by Gatchell and A. Holm.\cite{Holm2010, M. Gatchell2014, Gatchell2017}
Collisions between the high energetic projectile ions (such as 3 keV Ar$^+$, 22.5 keV He$^{2+}$) and neutral targets were investigated by M. Gatchell and A. Holm.\cite{Holm2010, Gatchell2014, Gatchell2017}
Collisions between clusters and projectile have been also explored at low collision energy, which allows for the derivation of dissociation energies and the thermal evaporation and stability of clusters. \cite{Boering1992, Wells2005, Zamith2019thermal}
By colliding a molecule, or a molecular aggregate, with a non-reactive rare gas atom (neon, argon) or a small molecule such as H$_2$O or N$_2$, it is possible to monitor the parent ions and collision products by use, for instance, of tandem mass spectrometry (MS/MS).\cite{Ma1997, Chowdhury2009} The resulting mass spectra provide a wealth of information about the structure of the parent and product ions from which one can infer, for instance, dissociation mechanisms \cite{Nelson1994, Molina2015} or bond and hydration enthalpies \cite{Carl2007}.

View File

@ -21,7 +21,7 @@
\setcounter{ContinuedFloat}{0}
\setcounter{pp@next@reset}{1}
\setcounter{@fnserial}{0}
\setcounter{NAT@ctr}{482}
\setcounter{NAT@ctr}{483}
\setcounter{Item}{5}
\setcounter{Hfootnote}{0}
\setcounter{bookmark@seq@number}{65}

View File

@ -466,68 +466,69 @@
\bibcite{Brechignac1994}{{418}{}{{}}{{}}}
\bibcite{Wong2004}{{419}{}{{}}{{}}}
\bibcite{Bush2008}{{420}{}{{}}{{}}}
\bibcite{Boering1992}{{421}{}{{}}{{}}}
\bibcite{Wells2005}{{422}{}{{}}{{}}}
\bibcite{Ma1997}{{423}{}{{}}{{}}}
\bibcite{Chowdhury2009}{{424}{}{{}}{{}}}
\bibcite{Sleno2004ion}{{425}{}{{}}{{}}}
\bibcite{Cody1982}{{426}{}{{}}{{}}}
\bibcite{Olsen2007higher}{{427}{}{{}}{{}}}
\bibcite{Hart2011}{{428}{}{{}}{{}}}
\bibcite{Gauthier1991}{{429}{}{{}}{{}}}
\bibcite{Laskin2005}{{430}{}{{}}{{}}}
\bibcite{Mcquinn2009}{{431}{}{{}}{{}}}
\bibcite{Carl2013}{{432}{}{{}}{{}}}
\bibcite{Hofstetter2013}{{433}{}{{}}{{}}}
\bibcite{Coates2017}{{434}{}{{}}{{}}}
\bibcite{Haag2009}{{435}{}{{}}{{}}}
\bibcite{Li1992}{{436}{}{{}}{{}}}
\bibcite{Bobbert2002}{{437}{}{{}}{{}}}
\bibcite{Spasov2000}{{438}{}{{}}{{}}}
\bibcite{Armentrout2008}{{439}{}{{}}{{}}}
\bibcite{Rodgers1998}{{440}{}{{}}{{}}}
\bibcite{Armentrout2007}{{441}{}{{}}{{}}}
\bibcite{Braud2017}{{442}{}{{}}{{}}}
\bibcite{Chirot2006new}{{443}{}{{}}{{}}}
\bibcite{Frenzel2004}{{444}{}{{}}{{}}}
\bibcite{Spiegelman2020}{{445}{}{{}}{{}}}
\bibcite{Simon2017}{{446}{}{{}}{{}}}
\bibcite{Rapacioli2018}{{447}{}{{}}{{}}}
\bibcite{Simon2018}{{448}{}{{}}{{}}}
\bibcite{Warshel1976}{{449}{}{{}}{{}}}
\bibcite{Cui2001}{{450}{}{{}}{{}}}
\bibcite{Iftner2014}{{451}{}{{}}{{}}}
\bibcite{Kukk2015}{{452}{}{{}}{{}}}
\bibcite{Rapacioli2018atomic}{{453}{}{{}}{{}}}
\bibcite{Dontot2019}{{454}{}{{}}{{}}}
\bibcite{Chung2011}{{455}{}{{}}{{}}}
\bibcite{Saggese2015}{{456}{}{{}}{{}}}
\bibcite{Mao2017}{{457}{}{{}}{{}}}
\bibcite{Wang2018}{{458}{}{{}}{{}}}
\bibcite{Kyrtopoulos2001}{{459}{}{{}}{{}}}
\bibcite{Farmer2003}{{460}{}{{}}{{}}}
\bibcite{Aumaitre2019}{{461}{}{{}}{{}}}
\bibcite{Rapacioli2005}{{462}{}{{}}{{}}}
\bibcite{Berne2008}{{463}{}{{}}{{}}}
\bibcite{Eschenbach1998}{{464}{}{{}}{{}}}
\bibcite{Goulart2017}{{465}{}{{}}{{}}}
\bibcite{Lei2019}{{466}{}{{}}{{}}}
\bibcite{Roser2015}{{467}{}{{}}{{}}}
\bibcite{Lemmens2019}{{468}{}{{}}{{}}}
\bibcite{Delaunay2015}{{469}{}{{}}{{}}}
\bibcite{Piacenza2005}{{470}{}{{}}{{}}}
\bibcite{Birer2015}{{471}{}{{}}{{}}}
\bibcite{Zhao2008truhlar}{{472}{}{{}}{{}}}
\bibcite{Bowal2019}{{473}{}{{}}{{}}}
\bibcite{Ricca2013}{{474}{}{{}}{{}}}
\bibcite{Grafenstein2009}{{475}{}{{}}{{}}}
\bibcite{Dontot2016}{{476}{}{{}}{{}}}
\bibcite{Dontot2020}{{477}{}{{}}{{}}}
\bibcite{Rapacioli2011}{{478}{}{{}}{{}}}
\bibcite{Gatchell2016}{{479}{}{{}}{{}}}
\bibcite{Gatchell2016knockout}{{480}{}{{}}{{}}}
\bibcite{Levine1987}{{481}{}{{}}{{}}}
\bibcite{Chen2014}{{482}{}{{}}{{}}}
\bibcite{Gatchell2014}{{421}{}{{}}{{}}}
\bibcite{Boering1992}{{422}{}{{}}{{}}}
\bibcite{Wells2005}{{423}{}{{}}{{}}}
\bibcite{Ma1997}{{424}{}{{}}{{}}}
\bibcite{Chowdhury2009}{{425}{}{{}}{{}}}
\bibcite{Sleno2004ion}{{426}{}{{}}{{}}}
\bibcite{Cody1982}{{427}{}{{}}{{}}}
\bibcite{Olsen2007higher}{{428}{}{{}}{{}}}
\bibcite{Hart2011}{{429}{}{{}}{{}}}
\bibcite{Gauthier1991}{{430}{}{{}}{{}}}
\bibcite{Laskin2005}{{431}{}{{}}{{}}}
\bibcite{Mcquinn2009}{{432}{}{{}}{{}}}
\bibcite{Carl2013}{{433}{}{{}}{{}}}
\bibcite{Hofstetter2013}{{434}{}{{}}{{}}}
\bibcite{Coates2017}{{435}{}{{}}{{}}}
\bibcite{Haag2009}{{436}{}{{}}{{}}}
\bibcite{Li1992}{{437}{}{{}}{{}}}
\bibcite{Bobbert2002}{{438}{}{{}}{{}}}
\bibcite{Spasov2000}{{439}{}{{}}{{}}}
\bibcite{Armentrout2008}{{440}{}{{}}{{}}}
\bibcite{Rodgers1998}{{441}{}{{}}{{}}}
\bibcite{Armentrout2007}{{442}{}{{}}{{}}}
\bibcite{Braud2017}{{443}{}{{}}{{}}}
\bibcite{Chirot2006new}{{444}{}{{}}{{}}}
\bibcite{Frenzel2004}{{445}{}{{}}{{}}}
\bibcite{Spiegelman2020}{{446}{}{{}}{{}}}
\bibcite{Simon2017}{{447}{}{{}}{{}}}
\bibcite{Rapacioli2018}{{448}{}{{}}{{}}}
\bibcite{Simon2018}{{449}{}{{}}{{}}}
\bibcite{Warshel1976}{{450}{}{{}}{{}}}
\bibcite{Cui2001}{{451}{}{{}}{{}}}
\bibcite{Iftner2014}{{452}{}{{}}{{}}}
\bibcite{Kukk2015}{{453}{}{{}}{{}}}
\bibcite{Rapacioli2018atomic}{{454}{}{{}}{{}}}
\bibcite{Dontot2019}{{455}{}{{}}{{}}}
\bibcite{Chung2011}{{456}{}{{}}{{}}}
\bibcite{Saggese2015}{{457}{}{{}}{{}}}
\bibcite{Mao2017}{{458}{}{{}}{{}}}
\bibcite{Wang2018}{{459}{}{{}}{{}}}
\bibcite{Kyrtopoulos2001}{{460}{}{{}}{{}}}
\bibcite{Farmer2003}{{461}{}{{}}{{}}}
\bibcite{Aumaitre2019}{{462}{}{{}}{{}}}
\bibcite{Rapacioli2005}{{463}{}{{}}{{}}}
\bibcite{Berne2008}{{464}{}{{}}{{}}}
\bibcite{Eschenbach1998}{{465}{}{{}}{{}}}
\bibcite{Goulart2017}{{466}{}{{}}{{}}}
\bibcite{Lei2019}{{467}{}{{}}{{}}}
\bibcite{Roser2015}{{468}{}{{}}{{}}}
\bibcite{Lemmens2019}{{469}{}{{}}{{}}}
\bibcite{Delaunay2015}{{470}{}{{}}{{}}}
\bibcite{Piacenza2005}{{471}{}{{}}{{}}}
\bibcite{Birer2015}{{472}{}{{}}{{}}}
\bibcite{Zhao2008truhlar}{{473}{}{{}}{{}}}
\bibcite{Bowal2019}{{474}{}{{}}{{}}}
\bibcite{Ricca2013}{{475}{}{{}}{{}}}
\bibcite{Grafenstein2009}{{476}{}{{}}{{}}}
\bibcite{Dontot2016}{{477}{}{{}}{{}}}
\bibcite{Dontot2020}{{478}{}{{}}{{}}}
\bibcite{Rapacioli2011}{{479}{}{{}}{{}}}
\bibcite{Gatchell2016}{{480}{}{{}}{{}}}
\bibcite{Gatchell2016knockout}{{481}{}{{}}{{}}}
\bibcite{Levine1987}{{482}{}{{}}{{}}}
\bibcite{Chen2014}{{483}{}{{}}{{}}}
\FN@pp@footnotehinttrue
\@input{6_backmatter/declaration.aux}
\providecommand\NAT@force@numbers{}\NAT@force@numbers

View File

@ -2554,6 +2554,13 @@
trivalent metal ions}.
\newblock {\em J. Am. Chem. Soc.}, {\bf 130}(28):9122--9128, 2008.
\bibitem{Gatchell2014}
{\sc Michael Gatchell, Patrick Rousseau, A~Domaracka, Mark~H Stockett, Tao
Chen, Henning~T Schmidt, JY~Chesnel, A~M{\'e}ry, S~Maclot, L~Adoui, et~al.}
\newblock {\bf Ions colliding with mixed clusters of C$_{60}$ and coronene:
Fragmentation and bond formation}.
\newblock {\em Phys. Rev. A}, {\bf 90}(2):022713, 2014.
\bibitem{Boering1992}
{\sc Kristie~A Boering, Joseph Rolfe, and John~I Brauman}.
\newblock {\bf Low energy collision induced dissociation: phase-shifting

View File

@ -9,14 +9,10 @@ A level-1 auxiliary file: 1_GeneIntro/GeneIntro.aux
A level-1 auxiliary file: 2_Introduction/introduction.aux
A level-1 auxiliary file: 3/structure_stability.aux
A level-1 auxiliary file: 4/collision.aux
White space in argument---line 6 of file 4/collision.aux
: \citation{Holm2010,M.
: Gatchell2014,Gatchell2017}
I'm skipping whatever remains of this command
A level-1 auxiliary file: 5/general_conclusion.aux
The style file: Latex/Classes/PhDbiblio-url2.bst
A level-1 auxiliary file: 6_backmatter/declaration.aux
Reallocated field_info (elt_size=4) to 18014 items from 5000.
Reallocated field_info (elt_size=4) to 18041 items from 5000.
Database file #1: 6_backmatter/references.bib
I was expecting a `,' or a `}'---line 1252 of file 6_backmatter/references.bib
:
@ -121,45 +117,45 @@ Warning--empty year in Unger1993
Warning--empty title in GaussianCode
Warning--empty journal in GaussianCode
Warning--empty year in GaussianCode
You've used 482 entries,
You've used 483 entries,
2776 wiz_defined-function locations,
3548 strings with 95460 characters,
and the built_in function-call counts, 222153 in all, are:
= -- 23195
> -- 9218
3553 strings with 95767 characters,
and the built_in function-call counts, 222834 in all, are:
= -- 23255
> -- 9276
< -- 22
+ -- 3680
- -- 3142
* -- 17439
:= -- 33780
add.period$ -- 1451
call.type$ -- 482
change.case$ -- 2521
+ -- 3703
- -- 3164
* -- 17505
:= -- 33895
add.period$ -- 1454
call.type$ -- 483
change.case$ -- 2533
chr.to.int$ -- 0
cite$ -- 494
duplicate$ -- 6839
empty$ -- 18808
format.name$ -- 3142
if$ -- 50990
cite$ -- 495
duplicate$ -- 6853
empty$ -- 18846
format.name$ -- 3164
if$ -- 51131
int.to.chr$ -- 0
int.to.str$ -- 482
missing$ -- 504
newline$ -- 2409
num.names$ -- 964
pop$ -- 3749
int.to.str$ -- 483
missing$ -- 505
newline$ -- 2414
num.names$ -- 966
pop$ -- 3763
preamble$ -- 1
purify$ -- 2520
purify$ -- 2532
quote$ -- 0
skip$ -- 11290
skip$ -- 11314
stack$ -- 0
substring$ -- 13655
swap$ -- 1968
substring$ -- 13678
swap$ -- 1972
text.length$ -- 22
text.prefix$ -- 0
top$ -- 0
type$ -- 1874
type$ -- 1878
warning$ -- 12
while$ -- 1746
width$ -- 485
write$ -- 5269
(There were 22 error messages)
while$ -- 1749
width$ -- 486
write$ -- 5280
(There were 21 error messages)

View File

@ -242,20 +242,20 @@
\backcite {Thomas1927}{{22}{2.3.2}{equation.2.3.16}}
\backcite {Fermi1928}{{22}{2.3.2}{equation.2.3.16}}
\backcite {Kohn1965}{{23}{2.3.2}{equation.2.3.18}}
\backcite {Kohn1965}{{25}{2.3.2}{equation.2.3.27}}
\backcite {Kohn1965}{{25}{2.3.2}{equation.2.3.27}}
\backcite {Herman1969}{{25}{2.3.2}{equation.2.3.27}}
\backcite {Perdew1996}{{25}{2.3.2}{equation.2.3.27}}
\backcite {Becke1993}{{25}{2.3.2}{equation.2.3.27}}
\backcite {Becke1993}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Becke1988}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Lee1988}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Stephens1994}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Adamo1999}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Henderson2009}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Krukau2006}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Zhao2008}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Li2018}{{26}{2.3.2}{equation.2.3.28}}
\backcite {Kohn1965}{{25}{2.3.2}{equation.2.3.26}}
\backcite {Kohn1965}{{25}{2.3.2}{equation.2.3.26}}
\backcite {Herman1969}{{25}{2.3.2}{equation.2.3.26}}
\backcite {Perdew1996}{{25}{2.3.2}{equation.2.3.26}}
\backcite {Becke1993}{{25}{2.3.2}{equation.2.3.26}}
\backcite {Becke1993}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Becke1988}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Lee1988}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Stephens1994}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Adamo1999}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Henderson2009}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Krukau2006}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Zhao2008}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Li2018}{{26}{2.3.2}{equation.2.3.27}}
\backcite {Foulkes1989}{{26}{2.3.3}{subsection.2.3.3}}
\backcite {Porezag1995}{{26}{2.3.3}{subsection.2.3.3}}
\backcite {Elstner1998}{{26}{2.3.3}{subsection.2.3.3}}
@ -265,29 +265,30 @@
\backcite {Elstne2007}{{26}{2.3.3}{subsection.2.3.3}}
\backcite {Oliveira2009}{{26}{2.3.3}{subsection.2.3.3}}
\backcite {Gaus2011}{{27}{2.3.3}{subsection.2.3.3}}
\backcite {Mulliken1955}{{30}{2.3.3}{equation.2.3.44}}
\backcite {Mulliken1955}{{31}{2.3.3}{equation.2.3.47}}
\backcite {Bader1985}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Bader1990}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Glendening1998}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Glendening2012}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Singh1984}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Besler1990}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Li1998}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Kalinowski2004}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Kelly2005}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Rapacioli2009corr}{{32}{2.3.3}{equation.2.3.51}}
\backcite {Mayer1983}{{32}{2.3.3}{equation.2.3.52}}
\backcite {Cances1997}{{32}{2.3.3}{equation.2.3.52}}
\backcite {Lu2012}{{32}{2.3.3}{equation.2.3.52}}
\backcite {Lewis1995}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Gianturco1999}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Elstner2001}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Wu2002}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Grimme2004}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Zimmerli2004}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Neumann2005}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Goursot2007}{{32}{2.3.3}{equation.2.3.53}}
\backcite {Mulliken1955}{{29}{2.3.3}{equation.2.3.40}}
\backcite {Elstner1998}{{30}{2.3.3}{equation.2.3.42}}
\backcite {Mulliken1955}{{30}{2.3.3}{equation.2.3.42}}
\backcite {Bader1985}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Bader1990}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Glendening1998}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Glendening2012}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Singh1984}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Besler1990}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Li1998}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Kalinowski2004}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Kelly2005}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Rapacioli2009corr}{{32}{2.3.3}{equation.2.3.48}}
\backcite {Mayer1983}{{32}{2.3.3}{equation.2.3.49}}
\backcite {Cances1997}{{32}{2.3.3}{equation.2.3.49}}
\backcite {Lu2012}{{32}{2.3.3}{equation.2.3.49}}
\backcite {Lewis1995}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Gianturco1999}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Elstner2001}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Wu2002}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Grimme2004}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Zimmerli2004}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Neumann2005}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Goursot2007}{{32}{2.3.3}{equation.2.3.50}}
\backcite {Frenkel2001}{{33}{2.3.4}{subsection.2.3.4}}
\backcite {Jones1924}{{33}{2.3.4}{subsection.2.3.4}}
\backcite {Lennard1924}{{33}{2.3.4}{subsection.2.3.4}}
@ -296,18 +297,18 @@
\backcite {Stillinger1985}{{33}{2.3.4}{subsection.2.3.4}}
\backcite {Kaplan2006}{{33}{2.3.4}{subsection.2.3.4}}
\backcite {Bahar1997}{{33}{2.3.4}{subsection.2.3.4}}
\backcite {Morse1929}{{34}{2.3.4}{equation.2.3.55}}
\backcite {Girifalco1959}{{34}{2.3.4}{equation.2.3.55}}
\backcite {Jones1924}{{34}{2.3.4}{equation.2.3.56}}
\backcite {Lennard1924}{{34}{2.3.4}{equation.2.3.56}}
\backcite {Lennard1931}{{34}{2.3.4}{equation.2.3.56}}
\backcite {Monticelli2013}{{34}{2.3.4}{equation.2.3.57}}
\backcite {Morse1929}{{33}{2.3.4}{equation.2.3.52}}
\backcite {Girifalco1959}{{33}{2.3.4}{equation.2.3.52}}
\backcite {Jones1924}{{33}{2.3.4}{equation.2.3.53}}
\backcite {Lennard1924}{{33}{2.3.4}{equation.2.3.53}}
\backcite {Lennard1931}{{34}{2.3.4}{equation.2.3.53}}
\backcite {Monticelli2013}{{34}{2.3.4}{equation.2.3.54}}
\backcite {Metropolis1987}{{36}{2.4.1}{subsection.2.4.1}}
\backcite {Metropolis1949}{{36}{2.4.1}{subsection.2.4.1}}
\backcite {Kroese2014}{{36}{2.4.1}{subsection.2.4.1}}
\backcite {Rosenbluth1955}{{37}{2.4.1}{subsection.2.4.1}}
\backcite {Binder1993}{{37}{2.4.1}{subsection.2.4.1}}
\backcite {Baeurle2009}{{37}{2.4.1}{subsection.2.4.1}}
\backcite {Rosenbluth1955}{{38}{2.4.1}{equation.2.4.65}}
\backcite {Binder1993}{{38}{2.4.1}{equation.2.4.65}}
\backcite {Baeurle2009}{{38}{2.4.1}{equation.2.4.65}}
\backcite {Alder1957}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Gibson1960}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Rahman1964}{{40}{2.4.2}{subsection.2.4.2}}
@ -342,27 +343,27 @@
\backcite {Vuilleumier1999}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Pavone2006}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Brancato2006}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Beck2000}{{41}{2.4.2}{subsection.2.4.2}}
\backcite {Wang2003}{{41}{2.4.2}{subsection.2.4.2}}
\backcite {Wei1994}{{41}{2.4.2}{subsection.2.4.2}}
\backcite {Light2000}{{41}{2.4.2}{subsection.2.4.2}}
\backcite {Batina1991}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Butcher2008}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Gear1966}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Diethelm2002}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Runge1895}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Kutta1901}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Butcher1996}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Butcher2015}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Verlet1967}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Fraige2004}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Swope1982}{{42}{2.4.2}{equation.2.4.71}}
\backcite {Ray1981}{{43}{2.4.2}{equation.2.4.80}}
\backcite {Ray1999}{{43}{2.4.2}{equation.2.4.80}}
\backcite {Nose1984}{{44}{2.4.2}{equation.2.4.80}}
\backcite {Nose1984M}{{44}{2.4.2}{equation.2.4.80}}
\backcite {Hoover1985}{{44}{2.4.2}{equation.2.4.80}}
\backcite {Hunenberger2005}{{44}{2.4.2}{equation.2.4.80}}
\backcite {Beck2000}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Wang2003}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Wei1994}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Light2000}{{40}{2.4.2}{subsection.2.4.2}}
\backcite {Batina1991}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Butcher2008}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Gear1966}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Diethelm2002}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Runge1895}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Kutta1901}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Butcher1996}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Butcher2015}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Verlet1967}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Fraige2004}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Swope1982}{{41}{2.4.2}{equation.2.4.68}}
\backcite {Ray1981}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Ray1999}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Nose1984}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Nose1984M}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Hoover1985}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Hunenberger2005}{{43}{2.4.2}{equation.2.4.77}}
\backcite {Torrie1977}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Hansmann1993}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Marchi1999}{{44}{2.4.3}{subsection.2.4.3}}
@ -373,18 +374,18 @@
\backcite {Barducci2011}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Kirkpatrick1983}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Van1987}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Swendsen1986}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Geyer1991}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Hukushima1996}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Falcioni1999}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Earl2005}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Sugita1999}{{45}{2.4.3}{subsection.2.4.3}}
\backcite {Wales1997}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Wales1999}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Hartke1993}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Unger1993}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Sivanandam2008}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Toledo2014}{{47}{2.4.4}{subsection.2.4.4}}
\backcite {Swendsen1986}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Geyer1991}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Hukushima1996}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Falcioni1999}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Earl2005}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Sugita1999}{{44}{2.4.3}{subsection.2.4.3}}
\backcite {Wales1997}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {Wales1999}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {Hartke1993}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {Unger1993}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {Sivanandam2008}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {Toledo2014}{{46}{2.4.4}{subsection.2.4.4}}
\backcite {deMonNano2009}{{50}{3.1.1}{subsection.3.1.1}}
\backcite {Elstner1998}{{50}{3.1.1}{subsection.3.1.1}}
\backcite {Gaus2013para}{{50}{3.1.1}{subsection.3.1.1}}
@ -590,6 +591,7 @@
\backcite {Bush2008}{{97}{4.1}{section.4.1}}
\backcite {Holm2010}{{98}{4.1}{section.4.1}}
\backcite {Gatchell2017}{{98}{4.1}{section.4.1}}
\backcite {Gatchell2014}{{98}{4.1}{section.4.1}}
\backcite {Zamith2019thermal}{{98}{4.1}{section.4.1}}
\backcite {Boering1992}{{98}{4.1}{section.4.1}}
\backcite {Wells2005}{{98}{4.1}{section.4.1}}

View File

@ -1,4 +1,4 @@
This is pdfTeX, Version 3.14159265-2.6-1.40.19 (TeX Live 2018) (preloaded format=pdflatex 2018.9.24) 15 JUN 2021 19:30
This is pdfTeX, Version 3.14159265-2.6-1.40.19 (TeX Live 2018) (preloaded format=pdflatex 2018.9.24) 15 JUN 2021 19:41
entering extended mode
restricted \write18 enabled.
%&-line parsing enabled.
@ -946,36 +946,38 @@ File: 2_Introduction/figures/variables.png Graphic file (type png)
Package pdftex.def Info: 2_Introduction/figures/variables.png used on input li
ne 415.
(pdftex.def) Requested size: 330.05322pt x 79.99901pt.
[22] [23 <./2_Introduction/figures/variables.png>]
LaTeX Warning: Reference `Lagrange' on page 24 undefined on input line 489.
[24] [25] [26] [27] [28] [29] [30] [31] [32] [33] [34]
<2_Introduction/figures/PES.png, id=1222, 1016.79875pt x 551.05875pt>
[22] [23 <./2_Introduction/figures/variables.png>] [24] [25] [26]
[27] [28] [29] [30] [31] [32] [33] [34]
<2_Introduction/figures/PES.png, id=1216, 1016.79875pt x 551.05875pt>
File: 2_Introduction/figures/PES.png Graphic file (type png)
<use 2_Introduction/figures/PES.png>
Package pdftex.def Info: 2_Introduction/figures/PES.png used on input line 995
Package pdftex.def Info: 2_Introduction/figures/PES.png used on input line 997
.
(pdftex.def) Requested size: 330.05322pt x 178.87334pt.
[35] [36 <./2_Introduction/figures/PES.png>] [37] [38] [39] [40] [41] [42] [43
] [44]
<2_Introduction/figures/PTMD_schema.png, id=1452, 865.7946pt x 656.2116pt>
[35 <./2_Introduction/figures/PES.png>] [36] [37] [38] [39] [40] [41] [42] [43
]
<2_Introduction/figures/PTMD_schema.png, id=1428, 865.7946pt x 656.2116pt>
File: 2_Introduction/figures/PTMD_schema.png Graphic file (type png)
<use 2_Introduction/figures/PTMD_schema.png>
Package pdftex.def Info: 2_Introduction/figures/PTMD_schema.png used on input
line 1353.
line 1349.
(pdftex.def) Requested size: 341.43306pt x 258.7852pt.
[45 <./2_Introduction/figures/PTMD_schema.png>] [46]) [47]
LaTeX Warning: `h' float specifier changed to `ht'.
[44] [45 <./2_Introduction/figures/PTMD_schema.png>] [46]) [47]
\openout2 = `3/structure_stability.aux'.
(./3/structure_stability.tex [48
(./3/structure_stability.tex [48
]
3.
[49] [50]
<3/figures/uracil.pdf, id=1535, 853.1875pt x 422.57875pt>
<3/figures/uracil.pdf, id=1526, 853.1875pt x 422.57875pt>
File: 3/figures/uracil.pdf Graphic file (type pdf)
<use 3/figures/uracil.pdf>
Package pdftex.def Info: 3/figures/uracil.pdf used on input line 82.
@ -985,7 +987,7 @@ Package pdftex.def Info: 3/figures/uracil.pdf used on input line 82.
LaTeX Warning: `!h' float specifier changed to `!ht'.
[51] [52 <./3/figures/uracil.pdf>] [53] [54]
<3/figures/E-distance-nh4-w.png, id=1686, 735.74875pt x 766.865pt>
<3/figures/E-distance-nh4-w.png, id=1677, 735.74875pt x 766.865pt>
File: 3/figures/E-distance-nh4-w.png Graphic file (type png)
<use 3/figures/E-distance-nh4-w.png>
Package pdftex.def Info: 3/figures/E-distance-nh4-w.png used on input line 200
@ -995,7 +997,7 @@ Package pdftex.def Info: 3/figures/E-distance-nh4-w.png used on input line 200
LaTeX Warning: `!h' float specifier changed to `!ht'.
<3/figures/E-distance-nh3-w.png, id=1687, 758.835pt x 763.85374pt>
<3/figures/E-distance-nh3-w.png, id=1678, 758.835pt x 763.85374pt>
File: 3/figures/E-distance-nh3-w.png Graphic file (type png)
<use 3/figures/E-distance-nh3-w.png>
Package pdftex.def Info: 3/figures/E-distance-nh3-w.png used on input line 208
@ -1006,7 +1008,7 @@ LaTeX Warning: `!h' float specifier changed to `!ht'.
[55] [56 <./3/figures/E-distance-nh4-w.png>] [57 <./3/figures/E-distance-nh3-w.
png>]
<3/figures/dimers.png, id=1710, 489.83pt x 425.59pt>
<3/figures/dimers.png, id=1701, 489.83pt x 425.59pt>
File: 3/figures/dimers.png Graphic file (type png)
<use 3/figures/dimers.png>
Package pdftex.def Info: 3/figures/dimers.png used on input line 245.
@ -1097,17 +1099,17 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
(hyperref) removing `math shift' on input line 258.
[58 <./3/figures/dimers.png>]
<3/figures/nh3-nh4-1w.png, id=1727, 312.16624pt x 555.07375pt>
<3/figures/nh3-nh4-1w.png, id=1718, 312.16624pt x 555.07375pt>
File: 3/figures/nh3-nh4-1w.png Graphic file (type png)
<use 3/figures/nh3-nh4-1w.png>
Package pdftex.def Info: 3/figures/nh3-nh4-1w.png used on input line 282.
(pdftex.def) Requested size: 82.51172pt x 146.71269pt.
<3/figures/nh3-nh4-2w.png, id=1728, 695.59875pt x 513.92pt>
<3/figures/nh3-nh4-2w.png, id=1719, 695.59875pt x 513.92pt>
File: 3/figures/nh3-nh4-2w.png Graphic file (type png)
<use 3/figures/nh3-nh4-2w.png>
Package pdftex.def Info: 3/figures/nh3-nh4-2w.png used on input line 291.
(pdftex.def) Requested size: 206.28247pt x 152.40494pt.
<3/figures/nh3-nh4-3w.png, id=1729, 747.79375pt x 888.31876pt>
<3/figures/nh3-nh4-3w.png, id=1720, 747.79375pt x 888.31876pt>
File: 3/figures/nh3-nh4-3w.png Graphic file (type png)
<use 3/figures/nh3-nh4-3w.png>
Package pdftex.def Info: 3/figures/nh3-nh4-3w.png used on input line 300.
@ -1166,7 +1168,7 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
(hyperref) removing `math shift' on input line 341.
<3/figures/nh4-4-6w.png, id=1765, 1189.44376pt x 887.315pt>
<3/figures/nh4-4-6w.png, id=1756, 1189.44376pt x 887.315pt>
File: 3/figures/nh4-4-6w.png Graphic file (type png)
<use 3/figures/nh4-4-6w.png>
Package pdftex.def Info: 3/figures/nh4-4-6w.png used on input line 353.
@ -1178,7 +1180,7 @@ LaTeX Warning: `!h' float specifier changed to `!ht'.
LaTeX Warning: Float too large for page by 46.56264pt on input line 408.
[61] [62] [63 <./3/figures/nh4-4-6w.png>] [64]
<3/figures/nh4-7-10w.png, id=1805, 1093.08376pt x 902.37125pt>
<3/figures/nh4-7-10w.png, id=1796, 1093.08376pt x 902.37125pt>
File: 3/figures/nh4-7-10w.png Graphic file (type png)
<use 3/figures/nh4-7-10w.png>
Package pdftex.def Info: 3/figures/nh4-7-10w.png used on input line 423.
@ -1224,7 +1226,7 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
(hyperref) removing `math shift' on input line 444.
<3/figures/nh3-4-7w.png, id=1837, 911.405pt x 861.2175pt>
<3/figures/nh3-4-7w.png, id=1828, 911.405pt x 861.2175pt>
File: 3/figures/nh3-4-7w.png Graphic file (type png)
<use 3/figures/nh3-4-7w.png>
Package pdftex.def Info: 3/figures/nh3-4-7w.png used on input line 449.
@ -1233,7 +1235,7 @@ Package pdftex.def Info: 3/figures/nh3-4-7w.png used on input line 449.
LaTeX Warning: `!h' float specifier changed to `!ht'.
[68] [69 <./3/figures/nh3-4-7w.png>] [70]
<3/figures/nh3-8-10w.png, id=1862, 937.5025pt x 854.19125pt>
<3/figures/nh3-8-10w.png, id=1852, 937.5025pt x 854.19125pt>
File: 3/figures/nh3-8-10w.png Graphic file (type png)
<use 3/figures/nh3-8-10w.png>
Package pdftex.def Info: 3/figures/nh3-8-10w.png used on input line 472.
@ -1291,7 +1293,7 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
(hyperref) removing `math shift' on input line 488.
<3/figures/nh3-nh4-20w.jpeg, id=1884, 465.74pt x 253.94875pt>
<3/figures/nh3-nh4-20w.jpeg, id=1875, 465.74pt x 253.94875pt>
File: 3/figures/nh3-nh4-20w.jpeg Graphic file (type jpg)
<use 3/figures/nh3-nh4-20w.jpeg>
Package pdftex.def Info: 3/figures/nh3-nh4-20w.jpeg used on input line 495.
@ -1309,64 +1311,64 @@ File: omlptm.fd
LaTeX Font Info: Font shape `OML/ptm/m/n' in size <10.95> not available
(Font) Font shape `OML/cmm/m/it' tried instead on input line 549.
[76]
<3/figures/mass7w.pdf, id=1962, 804.00375pt x 546.04pt>
<3/figures/mass7w.pdf, id=1953, 804.00375pt x 546.04pt>
File: 3/figures/mass7w.pdf Graphic file (type pdf)
<use 3/figures/mass7w.pdf>
Package pdftex.def Info: 3/figures/mass7w.pdf used on input line 582.
(pdftex.def) Requested size: 330.05322pt x 224.1528pt.
[77]
[78 <./3/figures/mass7w.pdf>]
<3/figures/fragcrosssec.pdf, id=2006, 893.3375pt x 589.20125pt>
<3/figures/fragcrosssec.pdf, id=1997, 893.3375pt x 589.20125pt>
File: 3/figures/fragcrosssec.pdf Graphic file (type pdf)
<use 3/figures/fragcrosssec.pdf>
Package pdftex.def Info: 3/figures/fragcrosssec.pdf used on input line 599.
(pdftex.def) Requested size: 330.05322pt x 217.68637pt.
<3/figures/Uloss.pdf, id=2012, 695.59875pt x 512.91624pt>
<3/figures/Uloss.pdf, id=2003, 695.59875pt x 512.91624pt>
File: 3/figures/Uloss.pdf Graphic file (type pdf)
<use 3/figures/Uloss.pdf>
Package pdftex.def Info: 3/figures/Uloss.pdf used on input line 604.
(pdftex.def) Requested size: 330.05322pt x 243.37164pt.
[79] [80 <./3/figures/fragcrosssec.pdf>]
<3/figures/protonAffinity.pdf, id=2040, 681.54625pt x 458.71375pt>
<3/figures/protonAffinity.pdf, id=2031, 681.54625pt x 458.71375pt>
File: 3/figures/protonAffinity.pdf Graphic file (type pdf)
<use 3/figures/protonAffinity.pdf>
Package pdftex.def Info: 3/figures/protonAffinity.pdf used on input line 614.
(pdftex.def) Requested size: 330.05322pt x 222.13992pt.
[81 <./3/figures/Uloss.pdf>] [82 <./3/figures/protonAffinity.pdf>]
<3/figures/a-b.pdf, id=2084, 437.635pt x 272.01625pt>
<3/figures/a-b.pdf, id=2075, 437.635pt x 272.01625pt>
File: 3/figures/a-b.pdf Graphic file (type pdf)
<use 3/figures/a-b.pdf>
Package pdftex.def Info: 3/figures/a-b.pdf used on input line 629.
(pdftex.def) Requested size: 206.28247pt x 128.21706pt.
[83 <./3/figures/a-b.pdf>]
<3/figures/1a-f.pdf, id=2112, 907.39pt x 718.685pt>
<3/figures/1a-f.pdf, id=2103, 907.39pt x 718.685pt>
File: 3/figures/1a-f.pdf Graphic file (type pdf)
<use 3/figures/1a-f.pdf>
Package pdftex.def Info: 3/figures/1a-f.pdf used on input line 670.
(pdftex.def) Requested size: 412.56496pt x 326.76093pt.
<3/figures/1a-f-b3lyp.pdf, id=2113, 997.7275pt x 757.83125pt>
<3/figures/1a-f-b3lyp.pdf, id=2104, 997.7275pt x 757.83125pt>
File: 3/figures/1a-f-b3lyp.pdf Graphic file (type pdf)
<use 3/figures/1a-f-b3lyp.pdf>
Package pdftex.def Info: 3/figures/1a-f-b3lyp.pdf used on input line 672.
(pdftex.def) Requested size: 412.56496pt x 313.36089pt.
[84]
<3/figures/2a-f.pdf, id=2124, 1054.94125pt x 807.015pt>
<3/figures/2a-f.pdf, id=2115, 1054.94125pt x 807.015pt>
File: 3/figures/2a-f.pdf Graphic file (type pdf)
<use 3/figures/2a-f.pdf>
Package pdftex.def Info: 3/figures/2a-f.pdf used on input line 676.
(pdftex.def) Requested size: 412.56496pt x 315.59656pt.
<3/figures/3a-f.pdf, id=2125, 1054.94125pt x 807.015pt>
<3/figures/3a-f.pdf, id=2116, 1054.94125pt x 807.015pt>
File: 3/figures/3a-f.pdf Graphic file (type pdf)
<use 3/figures/3a-f.pdf>
Package pdftex.def Info: 3/figures/3a-f.pdf used on input line 678.
(pdftex.def) Requested size: 412.56496pt x 315.59656pt.
[85 <./3/figures/1a-f.pdf>] [86 <./3/figures/1a-f-b3lyp.pdf>]
<3/figures/4a-f.pdf, id=2160, 1021.8175pt x 793.96625pt>
<3/figures/4a-f.pdf, id=2151, 1021.8175pt x 793.96625pt>
File: 3/figures/4a-f.pdf Graphic file (type pdf)
<use 3/figures/4a-f.pdf>
Package pdftex.def Info: 3/figures/4a-f.pdf used on input line 682.
(pdftex.def) Requested size: 412.56496pt x 320.56114pt.
<3/figures/5a-f.pdf, id=2161, 964.60374pt x 824.07875pt>
<3/figures/5a-f.pdf, id=2152, 964.60374pt x 824.07875pt>
File: 3/figures/5a-f.pdf Graphic file (type pdf)
<use 3/figures/5a-f.pdf>
Package pdftex.def Info: 3/figures/5a-f.pdf used on input line 684.
@ -1374,22 +1376,22 @@ Package pdftex.def Info: 3/figures/5a-f.pdf used on input line 684.
[87 <./3/figures/2a-f.pdf>] [88 <./3/figures/3a-f.pdf>] [89 <./3/figures/4a-f.p
df>]
<3/figures/6a-f.pdf, id=2215, 1035.87pt x 826.08624pt>
<3/figures/6a-f.pdf, id=2206, 1035.87pt x 826.08624pt>
File: 3/figures/6a-f.pdf Graphic file (type pdf)
<use 3/figures/6a-f.pdf>
Package pdftex.def Info: 3/figures/6a-f.pdf used on input line 693.
(pdftex.def) Requested size: 412.56496pt x 329.01688pt.
<3/figures/7a-f.pdf, id=2216, 922.44624pt x 814.04124pt>
<3/figures/7a-f.pdf, id=2207, 922.44624pt x 814.04124pt>
File: 3/figures/7a-f.pdf Graphic file (type pdf)
<use 3/figures/7a-f.pdf>
Package pdftex.def Info: 3/figures/7a-f.pdf used on input line 695.
(pdftex.def) Requested size: 412.56496pt x 364.07935pt.
<3/figures/11a-f.pdf, id=2217, 930.47626pt x 761.84625pt>
<3/figures/11a-f.pdf, id=2208, 930.47626pt x 761.84625pt>
File: 3/figures/11a-f.pdf Graphic file (type pdf)
<use 3/figures/11a-f.pdf>
Package pdftex.def Info: 3/figures/11a-f.pdf used on input line 696.
(pdftex.def) Requested size: 412.56496pt x 337.7941pt.
<3/figures/12a-f.pdf, id=2218, 955.57pt x 780.9175pt>
<3/figures/12a-f.pdf, id=2209, 955.57pt x 780.9175pt>
File: 3/figures/12a-f.pdf Graphic file (type pdf)
<use 3/figures/12a-f.pdf>
Package pdftex.def Info: 3/figures/12a-f.pdf used on input line 699.
@ -1446,19 +1448,15 @@ s/12a-f.pdf>]
]
4.
Package natbib Warning: Citation `M. Gatchell2014' on page 97 undefined on inpu
t line 33.
[97] [98] [99]
<4/figures/experiment-setup.pdf, id=2400, 1187.43625pt x 274.02374pt>
<4/figures/experiment-setup.pdf, id=2393, 1187.43625pt x 274.02374pt>
File: 4/figures/experiment-setup.pdf Graphic file (type pdf)
<use 4/figures/experiment-setup.pdf>
Package pdftex.def Info: 4/figures/experiment-setup.pdf used on input line 94.
(pdftex.def) Requested size: 412.56496pt x 95.20729pt.
[100] [101 <./4/figures/experiment-setup.pdf>]
<4/figures/exp-setup.pdf, id=2429, 959.585pt x 327.2225pt>
<4/figures/exp-setup.pdf, id=2422, 959.585pt x 327.2225pt>
File: 4/figures/exp-setup.pdf Graphic file (type pdf)
<use 4/figures/exp-setup.pdf>
Package pdftex.def Info: 4/figures/exp-setup.pdf used on input line 106.
@ -1467,56 +1465,56 @@ Package pdftex.def Info: 4/figures/exp-setup.pdf used on input line 106.
LaTeX Font Info: Font shape `OT1/ptm/bx/it' in size <10.95> not available
(Font) Font shape `OT1/ptm/b/it' tried instead on input line 128.
[103]
<4/figures/howinputs.pdf, id=2487, 1150.2975pt x 584.1825pt>
<4/figures/howinputs.pdf, id=2480, 1150.2975pt x 584.1825pt>
File: 4/figures/howinputs.pdf Graphic file (type pdf)
<use 4/figures/howinputs.pdf>
Package pdftex.def Info: 4/figures/howinputs.pdf used on input line 135.
(pdftex.def) Requested size: 412.56496pt x 209.5211pt.
[104 <./4/figures/howinputs.pdf>] [105]
<4/figures/3b-sphere.png, id=2517, 787.94376pt x 759.83875pt>
<4/figures/3b-sphere.png, id=2510, 787.94376pt x 759.83875pt>
File: 4/figures/3b-sphere.png Graphic file (type png)
<use 4/figures/3b-sphere.png>
Package pdftex.def Info: 4/figures/3b-sphere.png used on input line 188.
(pdftex.def) Requested size: 330.05322pt x 318.27213pt.
[106]
<4/figures/12f-sphere.png, id=2524, 784.9325pt x 770.88pt>
<4/figures/12f-sphere.png, id=2517, 784.9325pt x 770.88pt>
File: 4/figures/12f-sphere.png Graphic file (type png)
<use 4/figures/12f-sphere.png>
Package pdftex.def Info: 4/figures/12f-sphere.png used on input line 190.
(pdftex.def) Requested size: 330.05322pt x 324.1438pt.
[107 <./4/figures/3b-sphere.png>] [108 <./4/figures/12f-sphere.png>] [109] [11
0] [111]
<4/figures/proporEachFrag-1a2a.pdf, id=2572, 1046.91125pt x 505.89pt>
<4/figures/proporEachFrag-1a2a.pdf, id=2565, 1046.91125pt x 505.89pt>
File: 4/figures/proporEachFrag-1a2a.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-1a2a.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-1a2a.pdf used on input line
307.
(pdftex.def) Requested size: 412.56496pt x 199.35733pt.
<4/figures/proporEachFrag-3a4a-zoom.pdf, id=2573, 932.48375pt x 889.3225pt>
<4/figures/proporEachFrag-3a4a-zoom.pdf, id=2566, 932.48375pt x 889.3225pt>
File: 4/figures/proporEachFrag-3a4a-zoom.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-3a4a-zoom.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-3a4a-zoom.pdf used on input
line 309.
(pdftex.def) Requested size: 330.05322pt x 314.76857pt.
<4/figures/proporEachFrag-5a6a-zoom.pdf, id=2574, 932.48375pt x 890.32625pt>
<4/figures/proporEachFrag-5a6a-zoom.pdf, id=2567, 932.48375pt x 890.32625pt>
File: 4/figures/proporEachFrag-5a6a-zoom.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-5a6a-zoom.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-5a6a-zoom.pdf used on input
line 311.
(pdftex.def) Requested size: 330.05322pt x 315.12386pt.
<4/figures/proporEachFrag-11a-zoom.pdf, id=2575, 966.61125pt x 469.755pt>
<4/figures/proporEachFrag-11a-zoom.pdf, id=2568, 966.61125pt x 469.755pt>
File: 4/figures/proporEachFrag-11a-zoom.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-11a-zoom.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-11a-zoom.pdf used on input l
ine 313.
(pdftex.def) Requested size: 330.05322pt x 160.39554pt.
<4/figures/proporEachFrag-7a12a-zoom.pdf, id=2576, 932.48375pt x 892.33376pt>
<4/figures/proporEachFrag-7a12a-zoom.pdf, id=2569, 932.48375pt x 892.33376pt>
File: 4/figures/proporEachFrag-7a12a-zoom.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-7a12a-zoom.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-7a12a-zoom.pdf used on input
line 315.
(pdftex.def) Requested size: 330.05322pt x 315.8344pt.
<4/figures/proporEachFrag-7d12c-zoom.pdf, id=2577, 972.63374pt x 881.2925pt>
<4/figures/proporEachFrag-7d12c-zoom.pdf, id=2570, 972.63374pt x 881.2925pt>
File: 4/figures/proporEachFrag-7d12c-zoom.pdf Graphic file (type pdf)
<use 4/figures/proporEachFrag-7d12c-zoom.pdf>
Package pdftex.def Info: 4/figures/proporEachFrag-7d12c-zoom.pdf used on input
@ -1531,36 +1529,36 @@ LaTeX Warning: Float too large for page by 87.84254pt on input line 375.
ag-5a6a-zoom.pdf> <./4/figures/proporEachFrag-11a-zoom.pdf>] [115] [116 <./4/fi
gures/proporEachFrag-7a12a-zoom.pdf>] [117 <./4/figures/proporEachFrag-7d12c-zo
om.pdf>]
<4/figures/fig-1a-3b.pdf, id=2677, 939.51pt x 735.74875pt>
<4/figures/fig-1a-3b.pdf, id=2670, 939.51pt x 735.74875pt>
File: 4/figures/fig-1a-3b.pdf Graphic file (type pdf)
<use 4/figures/fig-1a-3b.pdf>
Package pdftex.def Info: 4/figures/fig-1a-3b.pdf used on input line 385.
(pdftex.def) Requested size: 330.05322pt x 258.47017pt.
<4/figures/fig-4a-5d.pdf, id=2678, 955.57pt x 808.01875pt>
<4/figures/fig-4a-5d.pdf, id=2671, 955.57pt x 808.01875pt>
File: 4/figures/fig-4a-5d.pdf Graphic file (type pdf)
<use 4/figures/fig-4a-5d.pdf>
Package pdftex.def Info: 4/figures/fig-4a-5d.pdf used on input line 387.
(pdftex.def) Requested size: 330.05322pt x 279.08733pt.
[118 <./4/figures/fig-1a-3b.pdf>] [119 <./4/figures/fig-4a-5d.pdf>]
<4/figures/fig-6a-6f.pdf, id=2714, 1031.855pt x 770.88pt>
<4/figures/fig-6a-6f.pdf, id=2707, 1031.855pt x 770.88pt>
File: 4/figures/fig-6a-6f.pdf Graphic file (type pdf)
<use 4/figures/fig-6a-6f.pdf>
Package pdftex.def Info: 4/figures/fig-6a-6f.pdf used on input line 405.
(pdftex.def) Requested size: 330.05322pt x 246.56902pt.
<4/figures/fig-7a-7d.pdf, id=2715, 618.31pt x 635.37375pt>
<4/figures/fig-7a-7d.pdf, id=2708, 618.31pt x 635.37375pt>
File: 4/figures/fig-7a-7d.pdf Graphic file (type pdf)
<use 4/figures/fig-7a-7d.pdf>
Package pdftex.def Info: 4/figures/fig-7a-7d.pdf used on input line 407.
(pdftex.def) Requested size: 330.05322pt x 339.17023pt.
[120 <./4/figures/fig-6a-6f.pdf>] [121 <./4/figures/fig-7a-7d.pdf>]
<4/figures/neutralUloss-Ne-Ar.pdf, id=2749, 853.1875pt x 853.1875pt>
<4/figures/neutralUloss-Ne-Ar.pdf, id=2742, 853.1875pt x 853.1875pt>
File: 4/figures/neutralUloss-Ne-Ar.pdf Graphic file (type pdf)
<use 4/figures/neutralUloss-Ne-Ar.pdf>
Package pdftex.def Info: 4/figures/neutralUloss-Ne-Ar.pdf used on input line 4
16.
(pdftex.def) Requested size: 330.05322pt x 330.05997pt.
<4/figures/cross-section-Ne-Ar.pdf, id=2750, 853.1875pt x 853.1875pt>
<4/figures/cross-section-Ne-Ar.pdf, id=2743, 853.1875pt x 853.1875pt>
File: 4/figures/cross-section-Ne-Ar.pdf Graphic file (type pdf)
<use 4/figures/cross-section-Ne-Ar.pdf>
Package pdftex.def Info: 4/figures/cross-section-Ne-Ar.pdf used on input line
@ -1605,24 +1603,24 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
[122 <./4/figures/neutralUloss-Ne-Ar.pdf>] [123 <./4/figures/cross-section-Ne-A
r.pdf>]
<4/figures/fig-11a-f.pdf, id=2795, 849.1725pt x 572.1375pt>
<4/figures/fig-11a-f.pdf, id=2788, 849.1725pt x 572.1375pt>
File: 4/figures/fig-11a-f.pdf Graphic file (type pdf)
<use 4/figures/fig-11a-f.pdf>
Package pdftex.def Info: 4/figures/fig-11a-f.pdf used on input line 448.
(pdftex.def) Requested size: 330.05322pt x 222.3732pt.
<4/figures/fig-12a-f.pdf, id=2796, 994.71625pt x 767.86874pt>
<4/figures/fig-12a-f.pdf, id=2789, 994.71625pt x 767.86874pt>
File: 4/figures/fig-12a-f.pdf Graphic file (type pdf)
<use 4/figures/fig-12a-f.pdf>
Package pdftex.def Info: 4/figures/fig-12a-f.pdf used on input line 451.
(pdftex.def) Requested size: 330.05322pt x 254.78004pt.
[124] [125 <./4/figures/fig-11a-f.pdf> <./4/figures/fig-12a-f.pdf>]
<4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf, id=2832, 974.64125pt x 865.2325pt>
<4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf, id=2825, 974.64125pt x 865.2325pt>
File: 4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf Graphic file (type pdf)
<use 4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf>
Package pdftex.def Info: 4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf used on input
line 459.
(pdftex.def) Requested size: 330.05322pt x 293.00015pt.
<4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf, id=2833, 975.645pt x 872.25874pt>
<4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf, id=2826, 975.645pt x 872.25874pt>
File: 4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf Graphic file (type pdf)
<use 4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf>
Package pdftex.def Info: 4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf used on input
@ -1630,7 +1628,7 @@ Package pdftex.def Info: 4/figures/MS-BR-5w-11w-Ne-Ar-branch.pdf used on input
(pdftex.def) Requested size: 330.05322pt x 295.07336pt.
[126 <./4/figures/MS-BR-1w-4w-Ne-Ar-branch.pdf>]
<4/figures/MS-BR-12w-Ne-branch.pdf, id=2851, 647.41875pt x 435.6275pt>
<4/figures/MS-BR-12w-Ne-branch.pdf, id=2844, 647.41875pt x 435.6275pt>
File: 4/figures/MS-BR-12w-Ne-branch.pdf Graphic file (type pdf)
<use 4/figures/MS-BR-12w-Ne-branch.pdf>
Package pdftex.def Info: 4/figures/MS-BR-12w-Ne-branch.pdf used on input line
@ -1675,99 +1673,99 @@ Package hyperref Warning: Token not allowed in a PDF string (PDFDocEncoding):
(hyperref) removing `math shift' on input line 509.
[129] [130] [131] [132] [133] [134]
<4/figures/MDPST.pdf, id=3010, 896.34875pt x 760.8425pt>
<4/figures/MDPST.pdf, id=3003, 896.34875pt x 760.8425pt>
File: 4/figures/MDPST.pdf Graphic file (type pdf)
<use 4/figures/MDPST.pdf>
Package pdftex.def Info: 4/figures/MDPST.pdf used on input line 626.
(pdftex.def) Requested size: 412.56496pt x 350.18915pt.
[135] [136 <./4/figures/MDPST.pdf>]
<4/figures/expTOF.pdf, id=3037, 1172.38pt x 413.545pt>
<4/figures/expTOF.pdf, id=3030, 1172.38pt x 413.545pt>
File: 4/figures/expTOF.pdf Graphic file (type pdf)
<use 4/figures/expTOF.pdf>
Package pdftex.def Info: 4/figures/expTOF.pdf used on input line 645.
(pdftex.def) Requested size: 412.56496pt x 145.52538pt.
[137 <./4/figures/expTOF.pdf>]
<4/figures/collisions.pdf, id=3062, 727.71875pt x 815.045pt>
<4/figures/collisions.pdf, id=3055, 727.71875pt x 815.045pt>
File: 4/figures/collisions.pdf Graphic file (type pdf)
<use 4/figures/collisions.pdf>
Package pdftex.def Info: 4/figures/collisions.pdf used on input line 658.
(pdftex.def) Requested size: 412.56496pt x 462.08072pt.
[138] [139 <./4/figures/collisions.pdf>]
<4/figures/distriPerc-Etf-175eV-d-bin03.pdf, id=3087, 1137.24875pt x 881.2925pt
<4/figures/distriPerc-Etf-175eV-d-bin03.pdf, id=3080, 1137.24875pt x 881.2925pt
>
File: 4/figures/distriPerc-Etf-175eV-d-bin03.pdf Graphic file (type pdf)
<use 4/figures/distriPerc-Etf-175eV-d-bin03.pdf>
Package pdftex.def Info: 4/figures/distriPerc-Etf-175eV-d-bin03.pdf used on in
put line 671.
(pdftex.def) Requested size: 330.05322pt x 255.76999pt.
<4/figures/fragmentation.pdf, id=3089, 1040.88875pt x 624.3325pt>
<4/figures/fragmentation.pdf, id=3082, 1040.88875pt x 624.3325pt>
File: 4/figures/fragmentation.pdf Graphic file (type pdf)
<use 4/figures/fragmentation.pdf>
Package pdftex.def Info: 4/figures/fragmentation.pdf used on input line 678.
(pdftex.def) Requested size: 412.56496pt x 247.46126pt.
[140]
<4/figures/opacitycurves.pdf, id=3100, 1137.24875pt x 853.1875pt>
<4/figures/opacitycurves.pdf, id=3093, 1137.24875pt x 853.1875pt>
File: 4/figures/opacitycurves.pdf Graphic file (type pdf)
<use 4/figures/opacitycurves.pdf>
Package pdftex.def Info: 4/figures/opacitycurves.pdf used on input line 685.
(pdftex.def) Requested size: 330.05322pt x 247.61331pt.
[141 <./4/figures/distriPerc-Etf-175eV-d-bin03.pdf>] [142 <./4/figures/fragmen
tation.pdf> <./4/figures/opacitycurves.pdf>]
<4/figures/cross-section.pdf, id=3145, 1137.24875pt x 853.1875pt>
<4/figures/cross-section.pdf, id=3138, 1137.24875pt x 853.1875pt>
File: 4/figures/cross-section.pdf Graphic file (type pdf)
<use 4/figures/cross-section.pdf>
Package pdftex.def Info: 4/figures/cross-section.pdf used on input line 701.
(pdftex.def) Requested size: 330.05322pt x 247.61331pt.
[143] [144 <./4/figures/cross-section.pdf>]
<4/figures/transferredE-Ar-300.pdf, id=3172, 1137.24875pt x 853.1875pt>
<4/figures/transferredE-Ar-300.pdf, id=3165, 1137.24875pt x 853.1875pt>
File: 4/figures/transferredE-Ar-300.pdf Graphic file (type pdf)
<use 4/figures/transferredE-Ar-300.pdf>
Package pdftex.def Info: 4/figures/transferredE-Ar-300.pdf used on input line
713.
(pdftex.def) Requested size: 330.05322pt x 247.61331pt.
[145 <./4/figures/transferredE-Ar-300.pdf>]
<4/figures/Epartition-Ar-300-SP.pdf, id=3202, 1137.24875pt x 910.40125pt>
<4/figures/Epartition-Ar-300-SP.pdf, id=3194, 1137.24875pt x 910.40125pt>
File: 4/figures/Epartition-Ar-300-SP.pdf Graphic file (type pdf)
<use 4/figures/Epartition-Ar-300-SP.pdf>
Package pdftex.def Info: 4/figures/Epartition-Ar-300-SP.pdf used on input line
767.
(pdftex.def) Requested size: 330.05322pt x 264.21797pt.
<4/figures/Epartition-Ar-300-Tstep-01.pdf, id=3205, 1137.24875pt x 853.1875pt>
<4/figures/Epartition-Ar-300-Tstep-01.pdf, id=3197, 1137.24875pt x 853.1875pt>
File: 4/figures/Epartition-Ar-300-Tstep-01.pdf Graphic file (type pdf)
<use 4/figures/Epartition-Ar-300-Tstep-01.pdf>
Package pdftex.def Info: 4/figures/Epartition-Ar-300-Tstep-01.pdf used on inpu
t line 796.
(pdftex.def) Requested size: 330.05322pt x 247.61331pt.
[146]
<4/figures/prot-Ar-300.pdf, id=3215, 1137.24875pt x 995.72pt>
<4/figures/prot-Ar-300.pdf, id=3208, 1137.24875pt x 995.72pt>
File: 4/figures/prot-Ar-300.pdf Graphic file (type pdf)
<use 4/figures/prot-Ar-300.pdf>
Package pdftex.def Info: 4/figures/prot-Ar-300.pdf used on input line 801.
(pdftex.def) Requested size: 330.05322pt x 288.97931pt.
[147 <./4/figures/Epartition-Ar-300-SP.pdf>] [148 <./4/figures/Epartition-Ar-3
00-Tstep-01.pdf>]
<4/figures/Epartition-Ar-300-d-ud.pdf, id=3257, 1137.24875pt x 1023.825pt>
<4/figures/Epartition-Ar-300-d-ud.pdf, id=3250, 1137.24875pt x 1023.825pt>
File: 4/figures/Epartition-Ar-300-d-ud.pdf Graphic file (type pdf)
<use 4/figures/Epartition-Ar-300-d-ud.pdf>
Package pdftex.def Info: 4/figures/Epartition-Ar-300-d-ud.pdf used on input li
ne 810.
(pdftex.def) Requested size: 330.05322pt x 297.13597pt.
[149 <./4/figures/prot-Ar-300.pdf>]
<4/figures/figuretimescale.pdf, id=3295, 674.52pt x 497.86pt>
<4/figures/figuretimescale.pdf, id=3288, 674.52pt x 497.86pt>
File: 4/figures/figuretimescale.pdf Graphic file (type pdf)
<use 4/figures/figuretimescale.pdf>
Package pdftex.def Info: 4/figures/figuretimescale.pdf used on input line 818.
(pdftex.def) Requested size: 330.05322pt x 243.61168pt.
<4/figures/T-time-zoom_abcdef.pdf, id=3296, 642.4pt x 892.33376pt>
<4/figures/T-time-zoom_abcdef.pdf, id=3289, 642.4pt x 892.33376pt>
File: 4/figures/T-time-zoom_abcdef.pdf Graphic file (type pdf)
<use 4/figures/T-time-zoom_abcdef.pdf>
Package pdftex.def Info: 4/figures/T-time-zoom_abcdef.pdf used on input line 8
22.
(pdftex.def) Requested size: 330.05322pt x 458.47456pt.
<4/figures/E-time-abcdef.pdf, id=3297, 794.97pt x 843.15pt>
<4/figures/E-time-abcdef.pdf, id=3290, 794.97pt x 843.15pt>
File: 4/figures/E-time-abcdef.pdf Graphic file (type pdf)
<use 4/figures/E-time-abcdef.pdf>
Package pdftex.def Info: 4/figures/E-time-abcdef.pdf used on input line 824.
@ -1861,7 +1859,7 @@ Underfull \hbox (badness 10000) in paragraph at lines 2497--2497
[]
[175] [176]
Underfull \hbox (badness 6252) in paragraph at lines 2841--2841
Underfull \hbox (badness 6252) in paragraph at lines 2848--2848
[]\OT1/ptm/m/sc/6 Annette Es-chen-bach, Rein-hard Wien-berg, and Bernd
[]
@ -1871,12 +1869,10 @@ Underfull \hbox (badness 6252) in paragraph at lines 2841--2841
(./6_backmatter/declaration.tex [179
])
Package natbib Warning: There were undefined citations.
Package atveryend Info: Empty hook `BeforeClearDocument' on input line 307.
Package atveryend Info: Empty hook `AfterLastShipout' on input line 307.
(./thesis.aux (./0_frontmatter/abstract.aux) (./0_frontmatter/dedication.aux)
(./thesis.aux
(./0_frontmatter/abstract.aux) (./0_frontmatter/dedication.aux)
(./0_frontmatter/acknowledgement.aux) (./0_frontmatter/glossary.aux)
(./1_GeneIntro/GeneIntro.aux) (./2_Introduction/introduction.aux)
(./3/structure_stability.aux) (./4/collision.aux) (./5/general_conclusion.aux)
@ -1886,18 +1882,14 @@ Package atveryend Info: Executing hook `AtEndAfterFileList' on input line 307.
Package rerunfilecheck Info: File `thesis.out' has not changed.
(rerunfilecheck) Checksum: 54B844088F193135A8E4AE44D25629B8;5769.
Package rerunfilecheck Info: File `thesis.brf' has not changed.
(rerunfilecheck) Checksum: B4BA215824F3A1FEC194964220F76936;39256.
LaTeX Warning: There were undefined references.
(rerunfilecheck) Checksum: E8709CD4E67FBBA5E85D00C43F8AA4F5;39356.
Package atveryend Info: Empty hook `AtVeryVeryEnd' on input line 307.
)
Here is how much of TeX's memory you used:
14705 strings out of 492649
225960 string characters out of 6129622
339468 words of memory out of 5000000
17242 multiletter control sequences out of 15000+600000
14706 strings out of 492649
226010 string characters out of 6129622
339212 words of memory out of 5000000
17245 multiletter control sequences out of 15000+600000
84285 words of font info for 177 fonts, out of 8000000 for 9000
1141 hyphenation exceptions out of 8191
43i,20n,51p,2185b,2512s stack positions out of 5000i,500n,10000p,200000b,80000s
@ -1913,10 +1905,10 @@ fonts/type1/urw/times/utmb8a.pfb></usr/local/texlive/2018/texmf-dist/fonts/type
1/urw/times/utmbi8a.pfb></usr/local/texlive/2018/texmf-dist/fonts/type1/urw/tim
es/utmr8a.pfb></usr/local/texlive/2018/texmf-dist/fonts/type1/urw/times/utmri8a
.pfb>
Output written on thesis.pdf (193 pages, 31313699 bytes).
Output written on thesis.pdf (193 pages, 31312211 bytes).
PDF statistics:
4478 PDF objects out of 5155 (max. 8388607)
3941 compressed objects within 40 object streams
923 named destinations out of 1000 (max. 500000)
4470 PDF objects out of 5155 (max. 8388607)
3933 compressed objects within 40 object streams
921 named destinations out of 1000 (max. 500000)
876 words of extra memory for PDF output out of 10000 (max. 10000000)

Binary file not shown.

Binary file not shown.

View File

@ -11,7 +11,7 @@
\contentsline {subsection}{\numberline {2.3.4}Force Field Methods}{33}{subsection.2.3.4}
\contentsline {section}{\numberline {2.4}Exploration of PES}{35}{section.2.4}
\contentsline {subsection}{\numberline {2.4.1}Monte Carlo Simulations}{36}{subsection.2.4.1}
\contentsline {subsection}{\numberline {2.4.2}Classical Molecular Dynamics}{40}{subsection.2.4.2}
\contentsline {subsection}{\numberline {2.4.2}Classical Molecular Dynamics}{39}{subsection.2.4.2}
\contentsline {subsection}{\numberline {2.4.3}Parallel-Tempering Molecular Dynamics}{44}{subsection.2.4.3}
\contentsline {subsection}{\numberline {2.4.4}Global Optimization}{46}{subsection.2.4.4}
\contentsline {chapter}{\numberline {3}Exploration of Structural and Energetic Properties}{49}{chapter.3}