transfer integral is an essential parameter in several theoretical models describing charge carrier transport in organic materials [ 17 ] . it is very often assumed that charge carrier mobility is proportional to the square of the charge transfer integral ( j ) which describes the transport of a charge between adjacent molecular sites . currently , the dft framework is commonly used to model charge transport in organic materials [ 3 , 814 ] . certain exchange correlation potentials are recognized to predict accurately geometries of molecules and shapes of molecular orbitals . however , it is well known , that wrong asymptotic behavior of conventional functionals create a real problem in calculations of some properties , especially for molecular complexes [ 15 , 16 ] . a recent systematic study of peach and co  workers may serve as an illustrative example . correlation functionals have difficulties with reliable description of excitation energies to charge  transfer states in molecules and molecular complexes . the charge for this reason , its evaluation might also present a challenge for the conventional exchange the primary goal of this study is to shed some light on this issue by employing recently proposed long range corrected density functional theory ( hereafter denoted as lrc  dft ) to compute charge dft is still being extensively tested primarily with an eye toward electric dipole ( hyper)polarizabilities , linear and nonlinear optical spectra [ 1729 ] . here , the lrc functionals employ the ewald split of \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ r_{{12}}^{{^ { { - { 1 } } } } } $ $ \end{document } operator which , in the case of the cam - b3lyp functional , takes the following form : 
 1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{1}{{{r_{{12 } } } } } = \frac{{1 - \left [ { \alpha + \beta \cdot { \text{erf}}\left ( { \mu { r_{{12 } } } } \right ) } \right]}}{{{r_{{12 } } } } } + \frac{{\alpha + \beta \cdot { \text{erf}}\left ( { \mu { r_{{12 } } } } \right)}}{{{r_{{12 } } } } } , $ $ \end{document}where the first and the second term on the rhs of the equation account for the short and the long  range interactions , respectively , and  and  are constants . here , we use the models based on the above equation with   =  0.33 for cam  b3lyp and   =  0.47 for lc  blyp . table  1the absolute values of effective charge transfer integrals ( |jeff| , given in ev ) computed with the aid of eqs . ( 3)(6 ) and the charge transfer integrals calculated using energy splitting in dimer approach ( |/2| ) . r is the intermolecular distance the absolute values of effective charge transfer integrals ( |jeff| , given in ev ) computed with the aid of eqs . ( 3)(6 ) and the charge transfer integrals calculated using energy splitting in dimer approach ( |/2| ) . r is the intermolecular distance as a model system to evaluate the performance of conventional exchange phthalocyanines are often considered as conductive materials with potential applications in organic electronics [ 3336 ] . in crystalline phase phthalocyanine molecules usually form regular columns and liquid crystals composed of phthalocyanines are promising materials for organic electronics . the liquid crystals in question are usually built from flat aromatic phthalocyanine center and aliphatic side groups . likewise , aromatic core of molecules in liquid crystal state form regular columns with molecules in stacked conformations and the fastest charge transport is observed inside a column with much smaller probability of charge transport between columns . the charge transfer integral between monomers in dimer can be used to describe charge transport inside of column composed of phthalocyanine molecules and as a first approximation of charge transfer in phthalocyanine based liquid crystals . in this work only charge  transfer integrals between highest occupied molecular orbitals ( homos ) of adjacent monomers are considered . the initial structure of the dimer was composed of two phthalocyanine monomers in stacked conformation with rotation axis passing through the centers of mass of monomers and perpendicular to the plane of monomer . the structure of the phthalocyanine monomer was optimized at the b3lyp/6 - 31g(d ) level of theory and was not reoptimized in dimer . the configurational space was spanned by a set of three parameters , namely intermolecular distance , twist angle around the symmetry axis and lateral slide of one of the monomers in two directions ( see fig .  1 ) . transfer integrals which are defined as : 
 2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { j_{{s , s\prime } } } = \left\langle { { \psi_s}\left| { { h_{{ks } } } } \right|{\psi_{{s\prime } } } } \right\rangle , $ $ \end{document}where hks is the kohn sham hamiltonian , s , s denote wave functions of the charge carrier localized on the sites s and s , respectively . in many organic systems zero spatial overlap between orbitals of the molecular sites . to account for this effect in calculations of charge transfer integral ( jeff ) may be introduced : 
 3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { j_{{eff } } } = { j_{{s , s\prime } } } - \frac{1}{2}{s_{{s , s\prime } } } \left ( { { \varepsilon_s } + { \varepsilon_{\prime } } } \right ) $ $ \end{document}ss , s denotes overlap integral of orbitals s and s ; s and s stand for energies of the sites s and s , respectively , and hereafter will be referred to as site energies . 1the schematic representation of three geometrical parameters used to describe the structure of phthalocyanine dimer :  denotes twist angle , r is intermolecular distance and l stands for lateral slide the schematic representation of three geometrical parameters used to describe the structure of phthalocyanine dimer :  denotes twist angle , r is intermolecular distance and l stands for lateral slide equation  2 can be directly used to compute the charge  transfer integral . in doing so , we express the hamiltonian of the dimer ( hks ) in the basis of molecular orbitals of the monomers [ 3 , 38 ] : 
 4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { h_{{ks } } } = sce{c^ { { - 1}}}. $ $ \end{document}in order to compute the charge transfer integral one needs to determine the eigenvalues for a dimer ( e ) , the eigenvectors in the basis of atomic orbitals ( ao ) for a dimer ( cao ) , the spatial overlap integrals in ao representation for a dimer ( sao ) and eigenvectors for monomers in ao representation \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( c_{{ao}}^{\prime } { \hbox{and } } c_{{ao}}^{{\prime\!\,\prime } } ) $ $ \end{document}.the eigenvector matrix for a dimer and spatial overlap matrix was transformed from ao representation to molecular orbital representation of monomers as follows : 
 5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ s = a{s_{{ao}}}{a^t } , $ $ \end{document}6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ c = a{c_{{ao}}}{a^t } , $ $ \end{document}where a denotes transformation matrix which is diagonal block matrix with monomer eigenvector matrices \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ ( c_{{ao}}^{\prime } { \hbox{and } } c_{{ao}}^{{\prime\!\,\prime } } ) $ $ \end{document } on the diagonal and a is transposed transformation matrix . transfer integrals , even more frequently used than the one described above , is the method which introduces energy splitting in dimer [ 3941 ] . in this approach , the eigenvalue of sth molecular orbital of the hamiltonian for dimer is given by : 
 7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { e_s } = \frac{{{\alpha_s}\pm { j_s}}}{{1\pm { s_s } } } , $ $ \end{document}where s is the s  th site energy integral , js and ss are the charge  transfer and the overlap integral between orbitals denoted by index s of two molecules forming dimer , respectively . it is assumed that the mixing between s  th orbital and the other molecular orbitals is not significant . the splitting in energy of s  th molecular orbital of monomer in a dimer can be written as : 
 8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left| { { e_{{s , a } } } - { e_{{s , b } } } } \right| = \left| { \frac{{2{j_s } - 2{s_s}{\alpha_s}}}{{1 - s_s^2 } } } \right| . ( 8) es , a and es , b denote eigenvalues for s  th molecular orbital of molecules a and b , respectively . transfer integral , site energy and overlap integral of s  th molecular orbital of molecules a and b , respectively . since we consider a dimer composed of two identical monomers , it is further assumed that s is the same for monomers a and b. for a system composed of two nonequivalent molecules , expressions for the eigenvalue es and the charge transfer integral js take more complicated form . usually , it is also assumed that spatial overlap integral is equal zero ( ss  =  0 ) , which is reasonable assumption considering charge transport between two organic molecules . in organic materials ( 8) can be rewritten as : 
 9\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \frac{1}{2}\left| { { e_{{s , a } } } - { e_{{s , b } } } } \right| \equiv \delta /2 = \left| { { j_s } } \right| . $ $ \end{document}in the present study we compute charge  transfer integrals using both above described methods . therefore , once the splitting in energy is known , the charge - transfer integral j can be calculated according to this relation . correlation potentials using different basis sets , including dunning s correlation consistent cc - pvdz basis set   as well as recently proposed jensen s basis set . the results of calculations presented in this work were carried out using the gaussian 09 program . most quantum  chemistry programs use the gaussian functions , introduced by boys in the mid 1950s , to approximate a wave function . transfer integrals using various basis sets and proved that basis set limit might be reached quite rapidly for the property in question . it should be underlined that the smooth convergence of charge  transfer integrals with the basis set size is observed for effective charge since the phthalocyanine complexes studied in the present work do not fulfill the latter criteria , we find it interesting to explore this topic more deeply by including the spatial overlap in computations of charge transfer integrals for different intermonomer separations are presented . indeed , as seen in the figure , the values of charge transfer integrals are quite insensitive to the basis set , provided the spatial overlap is taken into account . 2(a ) spatial overlap integrals ( s ) , charge  transfer integrals ( j ) and effective charge ( b ) dependence of the spatial overlap on the geometrical parameters calculated using various exchange pvdz basis set was used ( a ) spatial overlap integrals ( s ) , charge  transfer integrals ( j ) and effective charge ( b ) dependence of the spatial overlap on the geometrical parameters calculated using various exchange pvdz basis set was used the values of the spatial overlap for different geometrical parameters calculated with the aid of different exchange correlation potentials as well as using the hartree  fock method are presented in fig . as seen in the figure , the values of the spatial overlap calculated with use of different dft functionals are comparable for all considered geometrical parameters . the values of spatial overlap integral calculated using the hf wavefunction seem to be overestimated and differ substantially from the values determined within the dft framework . as seen in fig . 3 the choice of exchange  correlation potential is quite important as far as the magnitude of charge  transfer integral ( j ) and effective charge transfer integrals in conjugated systems in stacked configurations [ 3 , 41 ] , gives comparable results to those determined with the aid of the blyp and the b3lyp potentials . the values of charge  transfer integrals predicted by long  range corrected functionals , namely cam b3lyp and lc - blyp , lie between hf and conventional dft results ( see fig .  3 ) . it was shown by peach and co - workers that long  range corrected functionals improve substantially upon their traditional counterparts as far as excitation energies to charge  transfer states are concerned . it is a particularly notable observation for the coulomb  attenuated model ( cam - b3lyp ) . since both quantities in principle might be similar in nature , the lrc potentials should give better results also in the case of charge transfer integrals . for this reason , with a bit of scepticism due to the lack of more solid quantitative basis transfer integrals in comparison with their lrc counterparts . 
 fig .  3the dependence of the charge transfer integral ( j , left side ) and the effective charge  transfer integral ( jeff , right side ) on the geometrical parameters for different dft functionals using the cc - pvdz basis set the dependence of the charge transfer integral ( j , left side ) and the effective charge  transfer integral ( jeff , right side ) on the geometrical parameters for different dft functionals using the cc - pvdz basis set one of the theories most commonly used for describing charge transport in organic system is the one proposed by marcus [ 44 , 45 ] . in this approach , the hopping rate constant between sites i and j is defined as : 
 10\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { k_{{ij } } } = \frac{{2\pi { { \left| { { j_{{ij } } } } \right|}^2}}}{h}\sqrt { { \frac{\pi } { { 2{k_b}\lambda t } } } } \exp \left ( { - \frac{{{{\left ( { \lambda + { \varepsilon_j } - { \varepsilon_i } } \right)}^2}}}{{4\lambda { k_b}t } } } \right ) , $ $ \end{document}where h , kb and t are planck constant , boltzmann constant and temperature , respectively . i and j denote energies of the charge carrier localized on the sites i and j and  stands for reorganization energy [ 4 , 13 , 46 ] . in order to estimate the charge carrier mobility in the system without structural disorder ( assuming that each hopping rate is the same ) one can use the relation [ 47 , 48 ] : 
 11\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \mu = \frac{e}{{{k_b}t}}{k_{{ij}}}{a^2 } , $ $ \end{document}where e is elementary charge , and a denotes distance between molecular sites . for the studied phthalo - cyanine dimers , the internal reorganization energy calculated at the b3lyp / cc  pvdz level of theory is 0.043  ev , which is similar to the results presented in literature [ 13 , 49 , 50 ] . for the intermolecular distance 3.5   , the rotation angle 0 and lateral slide 1.5   ( this is the structure similar to the crystal structure of the phthalocyanine ) the charge thus , one can easily see , that dft functional has a significant influence on the mobility value 2 leads to the conclusion that for certain areas of conformational space the differences might be even higher . a comparison of the effective charge - transfer integrals calculated using eq . ( 3 ) and the charge - transfer integrals determined from the energy splitting ( eq . ( 9 ) ) is shown in table  1 . transfer integrals calculated based on the two approaches are insignificant and do not exceed a few thousandths of ev . at first glance , it appears that it is sufficient to employ less accurate method , based on the energy splitting in dimer with assumption of zero spatial overlap , to compute j between molecules in  interacting system . however , as it has already been mentioned , it is important to include spatial overlap in calculations of charge otherwise , the values of j might strongly depend on the size of the basis set used in calculations . the other drawback of the method based on energy splitting in dimer is the lack of information about the sign of charge however , if the knowledge of the sign is important , it can be subsequently determined from the bonding the primary aim of the present study was to evaluate the performance of commonly employed conventional exchange , we apply the recently proposed coulomb  attenuated model as a reference as this approach is proven to be very successful in predicting excitation energies to charge  transfer states . it is shown that for certain areas of conformational space in phthalocyanine dimer the differences in values of charge - transfer integrals between the conventional schemes and the cam - b3lyp functional in values of charge  transfer integrals might be quite significant . the same is revealed for triphenylene dimer . as a result , the values of charge carrier mobilities estimated using marcus formula might differ by 20% and more . likewise , theoretical predictions of peaks intensity in electro - absorption spectrum of molecular crystals and molecular aggregates [ 53 , 54 ] might be determined to a large extent by the accuracy of charge - transfer integrals ( kulig w , petelenz p , ( 2010 ) . we have also confirmed the findings reported by other authors   that the size of the basis set used in calculations of charge transfer integrals plays only a minor role provided the spatial overlap is included in the theoretical model .
an assessment of several widely used exchange -- correlation potentials in computing charge - transfer integrals is performed . in particular , we employ the recently proposed coulomb - attenuated model which was proven by other authors to improve upon conventional functionals in the case of charge - transfer excitations . for further validation , two distinct approaches to compute the property in question are compared for a phthalocyanine dimer .