therefore , they continuously remain of interest for non - linear optics , nanoelectronics , biotechnology and many other fields of science . from the point of view of material science , interaction - induced electric properties for hydrogen - bonded systems
however , systems that become relevant as novel materials are usually composed of more than several light atoms and therefore high - level calculations of their properties are prohibitively expensive .
the great abundance of the various available dft functionals allows to extend the investigations to the systems of hundreds of atoms and the wide range of analyzed properties .
it is known that the molecular electric properties can be reproduced by dft methods with good accuracy [ 112 ] .
however , one should bear in mind that for several classes of systems of importance to the non - linear optics , for instance long -conjugated polymer chains as well as push - pull systems with significant charge transfer , such calculations lead to significant deterioration in the accuracy of the results [ 1315 ] .
still , not much care has been devoted to the interaction - induced increments to dipole moments and ( hyper)polarizabilities calculated beyond the wave function theory .
, various dft functionals have been tested for the series of hydrogen - bonded dimers : water dimer , hf dimer and h2 co  hf .
although there is no unambiguous answer to the question which functional is best for interaction - induced electric properties , it is quite easy to reproduce the ccsd(t ) results for investigated systems with dft methods with satisfactory accuracy . regarding the gathered information , further studies seem necessary in order to generalize such an observation .
the aim of the present study is to verify the applicability of the dft formalism also for the interaction - induced electric properties of the longer hydrogen - bonded chains .
the hf linear chain , h2co crystal structure and h3n  (hf)n complex are investigated up to five units in the system .
the chosen structures were previously investigated at an angle of the many - body contributions to the ( non)linear electric response with second order mller - plesset perturbation theory and the many - body interactions in these systems proved to be significant [ 21 , 22 ] .
the performance of the dft functionals is tested against the ccsd(t ) data that provide the reference level . in order to keep the consistency with the previous study
, the same exchange - correlation functionals have been applied here : b3lyp [ 2326 ] and pbe0   as conventional hybrids , lc - blyp   as a long - range corrected functional and m06 - 2x hybrid meta functional .
additionally , this set was appended by the b3lyp long - range corrected modification ( cam - b3lyp )   for which we expect to improve the accuracy of the results toward reference data .
the subject of the present study are hydrogen - bonded complexes , namely : linear hf(hf)n , crystalline h2co(h2co)n and linear h3n(hf)n , for which n  =  14 . the structures of linear hf(hf)n   and h3n(hf)n   chains have been obtained based on gas - phase monomer geometries optimized at the mp2/6 - 31g(d , p ) level of theory .
moreover , the intermonomer distance which corresponds to the separation in the linear dimer has been optimized at the very same level of theory . on the other hand
, the geometry of the h2co(h2co)n chain refers to that found in the formaldehyde crystal .
all calculations of the static electric molecular properties were carried out at the scf - hf , mp2 , ccsd(t ) , b3lyp , pbe0 , lc - blyp , m06 - 2x and cam - b3lyp level of theory .
values of the electric dipole ( hyper)polarizability tensors have been obtained using the finite field ( ff ) method [ 33 , 34 ] by numerical differentiation of energy with respect to external electric field ( f ) . relying on our previous research [ 21 , 22 ] the electric field of the strength 0.001 a.u . has been chosen for further calculations to ensure the numerical stability of the ff procedure .
interaction - induced properties p for the investigated systems were obtained in a supermolecular manner:1\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \delta p = { { p}_{{{12 } \ldots n } } } - \sum\limits_{{i = 1}}^n { p_i , } $ $ \end{document}where p12
 n and pi stands for the magnitude of the property for the n - body system and its constituents , respectively . in order to make the consistent comparison of electric properties for the whole set of studied systems , the following values of the electric properties are discussed:2\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \left| \mu \right| = \sqrt { { \mu_x^2 { + } \mu_y^2 { + } \mu_z^2 , } } $ $ \end{document}3\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { { \alpha } _ 0 } = \frac{1}{3}\left ( { { { \alpha } _ { { xx } } } + { { \alpha } _ { { yy } } } + { { \alpha } _ { { zz } } } } \right ) , $ $ \end{document}4\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { { \beta } _ { \mu } } = \frac{{{{\beta } _ x}{{\mu } _ x } + { { \beta } _ y}{{\mu } _ y } + { { \beta } _ z}{{\mu } _
z}}}{{\left| \mu \right| } } , $ $ \end{document}where5\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \beta _ { i } = \frac{3 } { 5 } \cdot { \sum\limits_{j = x , y , z } { \beta _ { { ijj } } } } $ $ \end{document}and the kleinman symmetry is satisfied .
the values of static electric properties obtained within the supermolecular approach together with basis set of finite size might suffer from the basis set superposition error [ 35 , 36 ] . due to the fact that the first - order hyperpolarizability is the most sensitive to the choice of the basis set of all the properties analyzed here [ 35 , 36 ]
, the correction for bsse has been included according to the site - site function counterpoise scheme .
the total interaction - induced counterpoise - corrected electric property totalp has been decomposed into:6\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ { { \delta } _ { { total}}}{{p}^{{cp } } } = \sum\limits_{{i < j } } { { { \delta } _ { { ij } } } } p\left ( { 1 \ldots n } \right ) + \sum\limits_{{i < j < k } } { { { \delta } _ { { ijk } } } } p\left ( { 1 \ldots
n } \right ) + \ldots + { { \delta } _ { { 1 \ldots n}}}p\left ( { 1 \ldots n } \right ) , $ $ \end{document}where the first summation includes all the two - body terms ijp(1 
the relative error of the investigated properties in different methods has been estimated with respect to the ccsd(t ) calculations via:7\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \delta \left [ \% \right ] = \frac{{\delta { { p}^{{method } } } - \delta { { p}^{{ccsd(t)}}}}}{{\delta { { p}^{{ccsd(t ) } } } } } \cdot 100\% .
$ $ \end{document } similarly , the relative error for the bsse - uncorrected data has been predicted as:8\documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ bsse\left [ \% \right ] = \frac{{\delta { { p}^{{uncorr } } } - \delta { { p}^{{corr}}}}}{{\delta { { p}^{{corr } } } } } \cdot 100\% .
$ $ \end{document } all calculations have been performed with gamess - us   and gaussian09   program packages with the aug - cc - pvdz basis set [ 40 , 41 ] .
tables  1 , 2 and 3 present the values of the electric properties and their interaction - induced contributions for all the investigated systems in different approaches . the corresponding relative errors ( [% ] ) for \documentclass[12pt]{minimal } 
				 \usepackage{amsmath } 
				 \usepackage{wasysym } 
				 \usepackage{amsfonts } 
				 \usepackage{amssymb } 
				 \usepackage{amsbsy } 
				 \usepackage{mathrsfs } 
				 \usepackage{upgreek } 
				 \setlength{\oddsidemargin}{-69pt } 
				 \begin{document}$$ \delta \left| \mu \right|(0 ) $ $ \end{document } , 0(0 ) , and (0 ) are plotted in figs .  2 , 3 and 4 , respectively . for comparison the electric properties of the isolated units ( hf , h2co and h3n )
2the percentage relative errors of the approximation applied for estimation of the interaction - induced dipole moment with respect to the ccsd(t ) resultsfig .
3the percentage relative errors of the approximation applied for estimation of the interaction - induced polarizability with respect to the ccsd(t ) resultsfig .
4structure of investigated complexes of maximal lengthtable 1the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u . ] for the hydrogen - bonded linear hf(hf)n complexes .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body components||(0)||(0)0(0)0(0)(0)(0)ntotal2-bodym - bodytotal2-bodym - bodytotal2-bodym - bodyccsd(t)11.58580.17300.173010.0000.0480.04811.674.324.3222.49850.37890.36190.017015.1390.1770.1190.05811.6511.409.451.9533.42180.59520.55540.039820.3050.3340.1940.14010.6319.6414.794.8544.34940.81590.75100.064925.4820.5030.2700.2339.1628.3720.218.16mp211.59820.17390.173910.0480.0340.03410.974.184.1822.51760.38100.36400.017015.2010.1460.0900.05610.6111.289.451.8333.44770.59850.55890.039620.3790.2840.1490.1359.4319.3114.834.4844.38200.82030.75570.064625.5680.4340.2100.2247.8327.8020.297.51scf - hf11.69310.16770.16778.8930.0190.0199.613.173.1722.65270.36480.35090.013913.4320.0940.0550.03910.498.126.971.1533.62160.57110.53880.032317.9890.1880.0940.09410.8913.6210.852.7744.59420.78120.72860.052622.5550.2900.1340.15611.0619.3614.764.60b3lyp11.59760.18270.182710.4580.0960.09612.033.443.4422.52510.40180.38180.020015.8540.2770.2110.06610.4911.868.083.7833.46420.63230.58580.046521.2800.4900.3270.1637.1422.1712.719.4644.40790.86740.79180.075626.7200.7180.4450.2733.1733.2918.0315.26lc - blyp11.61520.18570.185710.8330.0810.0819.335.005.0022.55460.40870.38840.020316.4230.2590.1820.0775.9914.4011.552.8533.50600.64360.59610.047522.0460.4720.2870.1851.2425.4718.337.1444.46240.88340.80590.077527.6830.7020.3920.3104.1937.3025.2812.02pbe011.59890.17910.179110.2370.1010.10111.143.693.6922.52360.39340.37430.019115.5210.2830.2210.0629.7211.458.113.3433.45960.61890.57430.044620.8350.4970.3430.1546.7621.4113.398.0244.40010.84880.77620.072626.1630.7240.4660.2583.2031.8018.4213.38m06 - 2x11.60970.17520.17529.8700.0960.0968.973.093.0922.53630.38410.36610.018014.9580.2680.2140.0547.659.907.092.8133.47350.60350.56170.041820.0740.4690.3350.1345.0218.1011.696.4144.41510.82730.75930.068025.2020.6820.4580.2241.7626.9615.9910.97cam - b3lyp11.61610.18140.181410.3550.0790.07910.154.224.2222.55160.39840.37930.019115.6910.2430.1770.0668.1512.6310.062.5733.49850.62670.58210.044621.0550.4400.2780.1624.7622.6616.286.3844.45000.85950.78690.072626.4340.6500.3800.2700.7832.9421.2411.70table 2the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u . ] for the hydrogen - bonded linear h3n(hf)n complexes .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body components||(0)||(0)0(0)0(0)(0)(0)ntotal2-bodym - bodytotal2-bodym - bodytotal2-bodym - bodyccsd(t)11.72560.42940.429418.3660.6930.69339.358.388.3822.69840.69480.64430.050523.6400.4360.5380.10225.0314.2913.370.9233.63790.92730.84340.083928.8750.2140.4510.23717.0229.6622.137.5344.57141.15381.04100.112834.0760.0210.3720.35112.9141.1928.8112.38mp211.73680.42710.427118.4150.7700.77039.619.859.8522.71590.69300.64360.049423.6870.5450.6340.08926.2111.2911.390.1033.66200.92630.84410.082228.9280.3480.5620.21418.7925.6020.225.3844.60221.15381.04300.110834.1390.1770.4980.32114.8936.5326.969.57scf - hf11.79390.39870.398716.6550.6210.62123.7216.7716.7722.80310.64460.60450.040121.2820.4590.5320.07316.457.437.210.2233.78500.86360.79720.066425.8810.3250.4840.15913.2316.4713.173.3044.76241.07830.98890.089430.4610.2090.4410.23212.1023.5717.875.70b3lyp11.76040.45750.457518.9350.7210.72143.4811.4511.4522.75260.74070.68330.057424.4970.3820.4830.10118.9321.8423.241.4033.70940.98890.89300.095930.0230.0750.3500.2753.3444.5432.3912.1544.65981.23051.10150.129035.5000.1900.2560.4465.7461.2738.4322.84lc - blyp11.79990.45790.457919.1000.7500.75033.6113.7413.7422.80260.74360.68790.055724.8270.4400.5370.09715.0412.6615.182.5233.77140.99580.90130.094530.5330.1490.4160.2672.1732.2526.106.1544.73421.24201.11330.128736.2010.1100.3060.4166.4047.4734.4513.02pbe011.76390.45720.457218.5970.6410.64139.459.609.6022.75250.73460.67880.055824.0490.2970.3970.10017.4120.4621.901.4433.70590.97720.88420.093029.4550.0040.2600.2643.9141.0631.239.8344.65281.21351.08820.125334.8160.2650.1320.3973.9255.9037.6118.29m06 - 2x11.76790.44660.446618.0540.5940.59432.8511.9011.9022.75430.71460.66330.051323.2790.2860.3650.07916.5210.9012.871.9733.70770.95000.86410.085928.4670.0160.2270.2116.6526.4120.615.8044.65521.17951.06370.115833.6200.2230.1020.3250.7438.0225.3212.70cam - b3lyp11.78080.44930.449318.6020.6960.69634.7410.6710.6722.77780.72720.67310.054124.0780.3960.4970.10116.4615.2715.760.4933.74140.97210.88140.090729.5250.1240.3820.2584.7133.8625.558.3144.69911.21111.08840.122734.9320.1160.2750.3912.5347.7031.8415.86table 3the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u . ] for the hydrogen - bonded crystalline h2co(h2co)n complexes .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body componentsn||(0)||(0)0(0)0(0)(0)(0)total2-bodym - bodytotal2-bodym - bodytotal2-bodym - bodyccsd(t)12.03580.14320.143234.8680.3970.39793.099.409.4023.15400.31580.30270.013152.7320.9590.8460.113129.2323.5921.212.3834mp212.03960.14180.141834.9010.3910.39171.8412.3712.3723.15920.31310.29990.013252.7870.9560.8370.11995.6129.7727.722.0534.28810.49370.46280.030970.7281.5771.2970.280117.4949.0543.625.4345.42080.67810.62770.050488.6932.2231.7640.459138.5769.1259.949.18scf - hf12.46310.16510.165133.0780.3220.32293.101.821.8223.80930.36310.34920.013949.9540.7740.6890.085136.255.064.620.4435.16580.57130.53870.032666.8711.2681.0690.199178.808.947.611.3346.52650.78370.73070.053083.8051.7801.4550.325221.0613.1110.682.43b3lyp12.09930.15270.152735.7940.4390.439127.426.306.3023.25620.33710.32250.014654.1481.0440.9210.123176.6023.6616.387.2834.42310.53140.49740.034072.5571.7071.4160.291222.0944.7027.1417.5645.59410.72980.67430.055590.9902.3941.9170.477266.0267.3138.2829.03lc - blyp12.16470.15570.155735.8290.4440.44494.8110.2910.2923.35660.34370.32900.014754.2141.0630.9370.126129.6327.1823.883.3034.55870.54200.50750.034572.6561.7431.4450.298161.9046.6338.198.4445.76490.74450.68820.056391.1212.4471.9590.488193.0367.2252.7714.45pbe012.09490.14990.149935.2360.4530.453118.826.066.0623.24760.33070.31660.014153.3081.0680.9500.118165.9920.8615.135.7334.41000.52130.48820.033171.4331.7391.4600.279210.1338.6724.8313.8445.57640.71580.66180.054089.5812.4331.9750.458252.9857.7834.7922.99m06 - 2x12.13650.14730.147334.5930.4400.440100.177.747.7423.30830.32510.31140.013752.3251.0400.9250.115139.8321.3217.783.5434.48970.51250.48040.032170.1081.6941.4240.270177.4037.0128.408.6145.67500.70390.65160.052387.9132.3691.9280.441214.0953.5639.1114.45cam - b3lyp12.17130.15470.154735.2320.4170.417100.599.149.1423.36550.34110.32680.014353.2880.9990.8800.119139.1224.7521.293.4634.56970.53760.50400.033671.3971.6381.3570.281175.2942.7333.928.8145.77800.73820.68340.054889.5292.3001.8400.456210.3661.8346.8315.00table 4the values of the electric properties [ a.u . ] of analyzed monomershfh3nh2co||(0)ccsd(t)0.70590.59190.9470mp20.71180.59990.9496scf - hf0.76310.63501.1498b3lyp0.70630.59490.9740lc - blyp0.71290.62491.0053pbe00.70940.59930.9733m06 - 2x0.71680.60510.9954cam - b3lyp0.71610.61291.00900 ( 0)ccsd(t)4.94113.81517.152mp24.97113.89817.169scf - hf4.40912.65216.321b3lyp5.14514.19417.584lc - blyp5.33714.23117.595pbe05.03113.89317.307m06 - 2x4.85613.54217.005cam - b3lyp5.10313.91817.321 ( 0)ccsd(t)8.8922.2452.11mp28.0520.5642.98scf - hf6.9010.7148.60b3lyp8.7224.5566.95lc - blyp8.4913.7653.23pbe08.1821.4562.74m06 - 2x6.7215.9154.62cam - b3lyp8.3417.0855.43 structure of investigated complexes of maximal length the percentage relative errors of the approximation applied for estimation of the interaction - induced dipole moment with respect to the ccsd(t ) results the percentage relative errors of the approximation applied for estimation of the interaction - induced polarizability with respect to the ccsd(t ) results structure of investigated complexes of maximal length the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body components the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u . ] for the hydrogen - bonded linear h3n(hf)n complexes .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body components the values of the electric properties and their counterpoise corrected interaction - induced counterparts [ a.u . ] for the hydrogen - bonded crystalline h2co(h2co)n complexes .
n , the absolute value of the static dipole moment for the whole chain ||(0 ) , the total interaction - induced dipole moment ||(0 ) and its contributions : two- and many - body , static polarizability 0(0 ) , the total interaction - induced polarizability 0(0 ) and its two- and many - body components , first order hyperpolarizability (0 ) , the total interaction - induced hyperpolarizability (0 ) and its two- and many - body components the values of the electric properties [ a.u . ] of analyzed monomers the careful analysis of the gathered data for the interaction - induced dipole moments reveals quite a good performance of the tested dft functionals .
for most cases , the dft results are of better quality than scf - hf data .
however , the exceptions are observed for linear hf and h3n  (hf)nchains . in both cases lc - blyp performs noticeably worse than scf - hf approach with the relative error almost twice as large for the hf chain .
additionally , for the hf chain also the b3lyp functional and its long - range corrected version cam - b3lyp demonstrate slight aggravation with respect to the scf - hf approach .
the relative error for the dft calculations in all the investigated cases does not exceed 9% .
the worst data are obtained with lc - blyp , where the relative errors are between 7.5% and 8.5% and the best agreement with the ccsd(t ) interaction - induced dipole moments is observed for the m06 - 2x functional ( errors 1.5%-4% ) .
it is worth mentioning that the relative error for the mp2 calculations does not surpass 1% .
2 shows that the relative error of the interaction - induced dipole moment remains approximately constant with the chain elongation for the hf and h2co complexes . on the other hand , the lengthening of the h3n 
(hf)n system leads to the improvement of the m06 - 2x and pbe0 results - the errors decrease from 4% to 2% and from 6.5% to 5% , respectively .
the opposite behavior is observed for the lc - blyp data - the results worsen by more than 1% with increasing chain length .
the decomposition of the total interaction - induced dipole moment for the inspected complexes does not disclose any abrupt demeanor of the dft results : the quality of the obtained two- and more - body components does not raise many reservations and therefore the dft functionals tested in this study can be recommended for the many - body analysis of the interaction - induced dipole moment of the investigated complexes . due to the almost negligible non - additivity in the case of the polarizability of the analyzed molecular chains ( interaction - induced polarizability less than 1 a.u . for linear chains and less than 2.5 a.u . for formaldehyde
the dramatic increase of the error for the longest h3n  (hf)n complex should not be taken into account , since the interaction - induced increment for the ccsd(t ) calculations of the five - membered chain accounts for only 0.06% of the total polarizability value .
therefore , the interaction - induced polarizability within for instance pbe0 approach , that makes up 0.76% of the total pbe0 polarizability for the longest chain , appears to be in large error with respect to the ccsd(t ) data ; however this still gives a very good predictive power for the total polarizability values .
the contribution from many - body terms increases with the chain elongation and for long chains it may exceed the two - body terms . in the case of the h3n 
(hf)n complexes , where the two- and many - body terms have opposite signs , the expansion of the chain leads from the negative interaction - induced polarizability that weakens the total response of the system , to the positive 0 for larger number of hf units , where the interaction - induced contribution additionally strengthens the total polarizability of the system .
this tendency would be observed for all the investigated methods , however the applied dft functionals predict the length of the chain for which 0 changes sign as smaller than respective ccsd(t ) , mp2 or scf calculations .
in other words , dft overestimates the many - body terms with respect to the two - body contributions .
this is however not a general observation for dft functionals , since for the hf chains the propensity is opposite : dft seems to underestimate slightly the many - body terms in comparison to the two - body contributions , and in the case of the h2co chains the agreement between all the applied methods for the ratio of the many - body to two - body contributions to interaction - induced polarizability is very good .
all the calculations lead to the conclusion that again the mp2 method gives good quality results . additionally , the best - performing dft functional is cam - b3lyp . according to the author s expectations and experiences [ 20 , 42 ]
the worst quality results have been obtained definitely with the scf - hf approach - the errors ranging from 25% to 100% . for the formaldehyde complex
the dft functionals in general perform better than mp2 and the only exception is the pbe0 interaction - induced hyperpolarizability of the dimer .
taking mp2 approach as the suitable method for the incremental values estimation , one can treat its largest error ( 32% for the formaldehyde dimer ) as an indicator for the satisfactory performance of the dft functionals . for the linear hf chain all the investigated functionals achieve smaller errors than this limit . in the case of the formaldehyde complexes
the addition of one unit to the investigated dimers has a significant influence on the quality of the obtained results . with the chain elongation for the hf complexes the relative errors of the dft interaction - induced hyperpolarizability increase and in the case of pbe0 and b3lyp changes sign for trimer .
similar behavior might be expected for the m06 - 2x functional for chains longer than investigated . in other words , these functionals underestimate  for hf dimer and overestimate it for longer chains .
the dramatic decrease of the relative error value is observed for the lc - blyp , m06 - 2x and cam - b3lyp functionals when going from h3n  hf to h3n  (hf)2 .
thus , b3lyp , cam - b3lyp and pbe0 give the overestimation of  for all the analyzed lengths .
lc - blyp exhibits the underestimation of the  value for the h3n  (hf)2 and its overestimation besides .
moreover , m06 - 2x for three of four investigated h3n  (hf)n complexes underestimates  , however its relative errors could be expected to change sign again for the chains longer than investigated .
additionally the correct qualitative description of the total nonlinear response of the hf chain was achived by the dft calculations , which is in contrary to the scf - hf results ( observe the increase of the  to the absolute value ) .
the decomposition of the interaction - induced hyperpolarizability into the two- and more - body components reveals imperfections of the commonly applied hybrid functionals : b3lyp and pbe0 .
for all the investigated systems the many - body component calculated with these functionals are of poor quality with respect to the ccsd(t ) reference data . additionally , for the h3n 
(hf)n systems also the two - body terms are not reproduced well . on the other hand
, dft reproduces well the sign change of the two - body terms in the h3n 
(hf)n chains with the system elongation . for the many - body terms dft ( and mp2 ) also predicts the change of sign going from n  =  2 to n  =  3 , however this effect is not present in the ccsd(t ) reference results , when all the many - body terms are positive .
comparison of the ratio of the many - body to the two - body terms in the case of the interaction - induced hyperpolarizability indicates its overestimation for the hf and h2co chains and good correspondence in the case of h3n  (hf)n complexes .
additionally , the mutual relations of the monomer properties and interaction - induced contributions causes the total hyperpolarizability  to change signs for shorter hf and h3n
the careful analysis of the gathered data for the interaction - induced dipole moments reveals quite a good performance of the tested dft functionals . for most cases ,
however , the exceptions are observed for linear hf and h3n  (hf)nchains . in both cases lc - blyp performs noticeably worse than scf - hf approach with the relative error almost twice as large for the hf chain .
additionally , for the hf chain also the b3lyp functional and its long - range corrected version cam - b3lyp demonstrate slight aggravation with respect to the scf - hf approach .
the relative error for the dft calculations in all the investigated cases does not exceed 9% .
the worst data are obtained with lc - blyp , where the relative errors are between 7.5% and 8.5% and the best agreement with the ccsd(t ) interaction - induced dipole moments is observed for the m06 - 2x functional ( errors 1.5%-4% ) .
it is worth mentioning that the relative error for the mp2 calculations does not surpass 1% . the comparison of the data presented in fig . 
2 shows that the relative error of the interaction - induced dipole moment remains approximately constant with the chain elongation for the hf and h2co complexes . on the other hand , the lengthening of the h3n 
(hf)n system leads to the improvement of the m06 - 2x and pbe0 results - the errors decrease from 4% to 2% and from 6.5% to 5% , respectively .
the opposite behavior is observed for the lc - blyp data - the results worsen by more than 1% with increasing chain length .
the decomposition of the total interaction - induced dipole moment for the inspected complexes does not disclose any abrupt demeanor of the dft results : the quality of the obtained two- and more - body components does not raise many reservations and therefore the dft functionals tested in this study can be recommended for the many - body analysis of the interaction - induced dipole moment of the investigated complexes .
due to the almost negligible non - additivity in the case of the polarizability of the analyzed molecular chains ( interaction - induced polarizability less than 1 a.u . for linear chains and less than 2.5 a.u . for formaldehyde
the dramatic increase of the error for the longest h3n  (hf)n complex should not be taken into account , since the interaction - induced increment for the ccsd(t ) calculations of the five - membered chain accounts for only 0.06% of the total polarizability value .
therefore , the interaction - induced polarizability within for instance pbe0 approach , that makes up 0.76% of the total pbe0 polarizability for the longest chain , appears to be in large error with respect to the ccsd(t ) data ; however this still gives a very good predictive power for the total polarizability values .
the contribution from many - body terms increases with the chain elongation and for long chains it may exceed the two - body terms . in the case of the h3n 
(hf)n complexes , where the two- and many - body terms have opposite signs , the expansion of the chain leads from the negative interaction - induced polarizability that weakens the total response of the system , to the positive 0 for larger number of hf units , where the interaction - induced contribution additionally strengthens the total polarizability of the system .
this tendency would be observed for all the investigated methods , however the applied dft functionals predict the length of the chain for which 0 changes sign as smaller than respective ccsd(t ) , mp2 or scf calculations . in other words
, dft overestimates the many - body terms with respect to the two - body contributions .
this is however not a general observation for dft functionals , since for the hf chains the propensity is opposite : dft seems to underestimate slightly the many - body terms in comparison to the two - body contributions , and in the case of the h2co chains the agreement between all the applied methods for the ratio of the many - body to two - body contributions to interaction - induced polarizability is very good .
all the calculations lead to the conclusion that again the mp2 method gives good quality results .
according to the author s expectations and experiences [ 20 , 42 ] the analysis of the dft functionals assesment appeared to be most complex . figure  4 depicts the relative error for the interaction - induced hyperpolarizability .
the worst quality results have been obtained definitely with the scf - hf approach - the errors ranging from 25% to 100% . for the formaldehyde complex
the dft functionals in general perform better than mp2 and the only exception is the pbe0 interaction - induced hyperpolarizability of the dimer .
taking mp2 approach as the suitable method for the incremental values estimation , one can treat its largest error ( 32% for the formaldehyde dimer ) as an indicator for the satisfactory performance of the dft functionals .
for the linear hf chain all the investigated functionals achieve smaller errors than this limit . in the case of the formaldehyde complexes
the addition of one unit to the investigated dimers has a significant influence on the quality of the obtained results . with the chain elongation for the hf complexes the relative errors of the dft interaction - induced hyperpolarizability increase and in the case of pbe0 and b3lyp changes sign for trimer .
similar behavior might be expected for the m06 - 2x functional for chains longer than investigated . in other words ,
the dramatic decrease of the relative error value is observed for the lc - blyp , m06 - 2x and cam - b3lyp functionals when going from h3n  hf to h3n  (hf)2 .
thus , b3lyp , cam - b3lyp and pbe0 give the overestimation of  for all the analyzed lengths .
lc - blyp exhibits the underestimation of the  value for the h3n  (hf)2 and its overestimation besides . moreover , m06 - 2x for three of four investigated h3n 
(hf)n complexes underestimates  , however its relative errors could be expected to change sign again for the chains longer than investigated .
additionally the correct qualitative description of the total nonlinear response of the hf chain was achived by the dft calculations , which is in contrary to the scf - hf results ( observe the increase of the  to the absolute value ) .
the decomposition of the interaction - induced hyperpolarizability into the two- and more - body components reveals imperfections of the commonly applied hybrid functionals : b3lyp and pbe0 .
for all the investigated systems the many - body component calculated with these functionals are of poor quality with respect to the ccsd(t ) reference data . additionally , for the h3n 
on the other hand , dft reproduces well the sign change of the two - body terms in the h3n 
(hf)n chains with the system elongation . for the many - body terms dft ( and mp2 ) also predicts the change of sign going from n  =  2 to n  =  3 , however this effect is not present in the ccsd(t ) reference results , when all the many - body terms are positive .
comparison of the ratio of the many - body to the two - body terms in the case of the interaction - induced hyperpolarizability indicates its overestimation for the hf and h2co chains and good correspondence in the case of h3n  (hf)n complexes .
additionally , the mutual relations of the monomer properties and interaction - induced contributions causes the total hyperpolarizability  to change signs for shorter hf and h3n  (hf)n chains within the dft approaches than within the ccsd(t ) calculations .
the estimated bsse values for the investigated complexes are summarized in tables  5 , 6 and 7 .
it should be repeated that all the calculations have been carried out with the aug - cc - pvdz basis set of the medium size . the relative error ( bsse[% ] )
arising from the lack of the bsse correction for the dipole moment is usually smaller than 1% and only in the case of the linear hf chain for the lc - blyp functional it reaches 2% . on these basis
one can conclude that the counterpoise procedure has a negligible influence on the quality and quantity of the obtained interaction - induced dipole moment values for the studied complexes.table 5interaction - induced electric properties for the hydrogen - bonded linear hf 
(hf)n complexes with and without counterpoise correction and the basis set superposition error [ % ] n||(0)0(0)(0)cp - corcp - uncorbssecp - corcp - uncorbssecp - corcp - uncorbsseccsd(t)10.17300.17390.52%0.0480.118145.83%4.326.1041.20%20.37890.38070.48%0.1770.31678.53%11.4015.0131.67%30.59520.59800.47%0.3340.54061.68%19.6424.9226.88%40.81590.81960.45%0.5030.77754.47%28.3735.2824.36%mp210.17390.17460.40%0.0340.106211.76%4.185.1222.49%20.38100.38230.34%0.1460.28897.26%11.2813.5219.86%30.59850.60050.33%0.2840.49574.30%19.3122.7517.81%40.82030.82310.34%0.4340.71364.29%27.8032.3916.51%scf - hf10.16770.16700.42%0.0190.075294.74%3.174.1831.86%20.36480.36360.33%0.0940.206119.15%8.1210.2025.62%30.57110.56940.30%0.1880.35488.30%13.6216.6922.54%40.78120.77890.29%0.2900.51176.21%19.3623.4220.97%b3lyp10.18270.18501.26%0.0960.16773.96%3.445.4157.27%20.40180.40621.10%0.2770.41850.90%11.8615.6632.04%30.63230.63891.04%0.4900.69842.45%22.1727.7425.12%40.86740.87641.04%0.7180.99338.30%33.2940.4321.45%lc - blyp10.18570.18941.99%0.0810.15895.06%5.007.6553.00%20.40870.41591.76%0.2590.41158.69%14.4019.4735.21%30.64360.65441.68%0.4720.69647.46%25.4732.7228.46%40.88340.89781.63%0.7020.99641.88%37.3046.6325.01%pbe010.17910.18010.56%0.1010.17472.28%3.695.2341.73%20.39340.39540.41%0.2830.42750.88%11.4514.8329.52%30.61890.62190.48%0.4970.70942.66%21.4125.9821.35%40.84880.85300.49%0.7241.00538.81%31.8037.7218.62%m06 - 2x10.17520.17620.57%0.0960.15763.54%3.094.4744.66%20.38410.38600.49%0.2680.38844.78%9.9012.5026.26%30.60350.60650.50%0.4690.64838.17%18.1021.8620.77%40.82730.83130.48%0.6820.91934.75%26.9631.8318.06%cam - b3lyp10.18140.18401.43%0.0790.14887.34%4.226.5354.59%20.39840.40341.26%0.2430.38056.38%12.6316.8733.57%30.62670.63421.20%0.4400.64245.91%22.6628.5926.17%40.85950.86971.19%0.6500.91741.08%32.9440.9124.20%table 6interaction - induced electric properties for the hydrogen - bonded h3n  (hf)n complexes with and without counterpoise correction and the basis set superposition error [ % ] n||(0)0(0)(0)cp - corcp - uncorbssecp - corcp - uncorbssecp - corcp - uncorbsseccsd(t)10.42940.42770.40%0.6930.39043.72%8.388.221.91%20.69480.69450.04%0.4360.05786.93%14.2914.984.83%30.92730.92820.10%0.2140.236210.28%29.6631.887.48%41.15381.15570.16%0.0210.4972466.67%41.1944.888.96%mp210.42710.42510.47%0.7700.45341.17%9.8511.0011.68%20.69300.69240.09%0.5450.15371.93%11.2910.447.53%30.92630.92670.04%0.3480.117133.62%25.6025.901.17%41.15391.15520.11%0.1770.358302.26%36.5337.853.61%scf - hf10.39870.39580.73%0.6210.40634.62%16.776.1163.57%20.64460.64200.40%0.4590.18859.04%7.438.058.34%30.86360.86080.32%0.3250.003100.92%16.4718.1710.32%41.07831.07520.29%0.2090.174183.25%23.5726.2011.16%b3lyp10.45750.45920.37%0.7210.40543.83%11.4510.2110.83%20.74070.74500.58%0.3820.012103.14%21.8423.065.59%30.98890.99560.68%0.0750.393624.00%44.5447.376.35%41.23051.23960.74%0.1900.725281.58%61.2765.176.37%lc - blyp10.45790.46210.92%0.7500.46837.60%13.7411.3617.32%20.74360.75191.12%0.4400.07982.05%12.6615.7024.01%30.99581.00781.21%0.1490.290294.63%32.2537.0614.91%41.24201.25771.26%0.1100.620463.64%47.4754.1113.99%pbe010.45720.45520.44%0.6410.32748.99%9.609.822.29%20.73460.73440.03%0.2970.094131.65%20.4620.400.29%30.97720.97830.11%0.0040.46811600.00%41.0642.102.53%41.21351.21580.19%0.2650.797200.75%55.9058.103.94%m06 - 2x10.44660.44610.11%0.5940.34541.92%11.9010.2214.12%20.71460.71580.17%0.2860.024108.39%10.9012.8317.71%30.95000.95230.24%0.0160.3552318.75%26.4129.4211.40%41.17951.18310.31%0.2230.652192.38%38.0242.0510.60%cam - b3lyp10.44930.45180.55%0.6960.41939.80%10.679.3212.65%20.72720.73270.75%0.3960.04688.38%15.2717.3013.29%30.97210.98030.84%0.1240.297339.52%33.8637.3910.43%41.21111.22190.89%0.1160.601418.10%47.7052.9711.05%table 7interaction - induced electric properties for the hydrogen - bonded crystalline h2co  (h2co)n complexes with and without counterpoise correction and the basis set superposition error [ % ] n||(0)0(0)(0)cp - corcp - uncorbssecp - corcp - uncorbssecp - corcp - uncorbsseccsd(t)10.14320.14180.98%0.3970.56341.81%9.4011.1218.30%20.31580.31300.89%0.9591.27532.95%23.5927.0914.84%34mp210.14180.14040.99%0.3910.56343.99%12.3714.1314.23%20.31310.31030.89%0.9561.28134.00%29.7733.3411.99%30.49370.48950.85%1.5782.05330.10%49.0554.4410.99%40.67810.67260.81%2.2232.84928.16%69.1276.3510.46%scf - hf10.16510.16360.91%0.3220.43635.40%1.824.10125.27%20.36310.36000.85%0.7740.99228.17%5.069.5488.54%30.57130.56670.81%1.2681.58825.24%8.9415.5974.38%40.78370.77770.77%1.7802.20223.71%13.1121.9267.20%b3lyp10.15270.15120.98%0.4390.62742.82%6.306.482.86%20.33710.33410.89%1.0441.39733.81%23.6624.262.54%30.53140.52700.83%1.7072.22230.17%44.7045.722.28%40.72980.72390.81%2.3943.07228.32%67.3168.742.12%lc - blyp10.15570.15411.03%0.4440.63943.92%10.2911.6513.22%20.34370.34070.87%1.0631.42934.43%27.1830.0510.56%30.54200.53750.83%1.7432.27630.58%46.6351.019.39%40.74450.73840.82%2.4473.14728.61%67.2273.118.76%pbe010.14990.14841.00%0.4530.62237.31%6.066.669.90%20.33070.32780.88%1.0681.38729.87%20.8622.256.66%30.52130.51700.82%1.7392.20426.74%38.6740.855.64%40.71580.71010.80%2.4333.04525.15%57.7860.745.12%m06 - 2x10.14730.14571.09%0.4400.58432.73%7.749.0617.05%20.32510.32210.92%1.0401.31126.06%21.3224.0212.66%30.51250.50800.88%1.6942.08823.26%37.0141.0710.97%40.70390.69790.85%2.3692.88821.91%53.5658.9910.14%cam - b3lyp10.15470.15330.90%0.4170.59141.73%9.1410.2612.25%20.34110.33840.79%0.9991.32532.63%24.7527.169.74%30.53760.53350.76%1.6382.11429.06%42.7346.428.64%40.73820.73280.73%2.3002.92527.17%61.8366.777.99% interaction - induced electric properties for the hydrogen - bonded linear hf  (hf)n complexes with and without counterpoise correction and the basis set superposition error [ % ] interaction - induced electric properties for the hydrogen - bonded h3n 
(hf)n complexes with and without counterpoise correction and the basis set superposition error [ % ] interaction - induced electric properties for the hydrogen - bonded crystalline h2co  (h2co)n complexes with and without counterpoise correction and the basis set superposition error [ % ] the analysis of the relative bsse for the linear hf and h3n  (hf)n complex is unreasonable due to the very small 0 values .
however , for the formaldehyde chains the obtained results become more stable owing to the larger interaction - induced polarizability values than in the case of the remaining systems . here , the relative bsse changes from above 40% for short chains to almost 20% for longer chains .
nevertheless , again the influence of the interaction - induced bsse for the total value of the polarizability of the whole complex is negligible and the qualitative picture of the 0 changes with the chain elongation remains the same with or without the counterpoise correction .
the changes of the relative bsse for the interaction - induced hyperpolarizability are parallel to these observed for the dipole moment , however here the absolute values of the errors decrease from 60% to 20% for hf chains and from 20% to 2% for formaldehyde chains with the chain elongation .
again , similarly as in the case of dipole moment and polarizability , the counterpoise procedure does not affect much the total qualitative picture .
however , it can not be neglected for the correct quantitative description of the interaction - induced hyperpolarizability of the investigated systems [ 35 , 36 ] .
the investigation of the electric properties of the hydrogen - bonded dimers   has shown that the m06 - 2x functional performed worst mainly for hyperpolarizability ( compare hf dimer ) .
however in the case of longer chains , it presents a serious advantage over any other tested functionals .
therefore , m06 - 2x appears to be the best for the description of the influence of the mutual interaction on the electric properties in the hydrogen - bonded chains .
although these conclusions were verified in the present contribution only for one basis set of the medium size ( aug - cc - pvdz ) , they should not differ dramatically for other basis set choices .
the ordering of the functionals may vary , however the main conclusion of this work remain unchanged : dft functionals can be applied succesfully for the investigation of the interaction - induced electric properties .
one should however keep in mind that the dft functional performance can also significantly depend on the nature of the intermolecular interactions and the conventional functionals are known to fail for systems governed by dispersion interaction [ 4346 ] .
it is also of importance for the complex systems , that the basis set superposition error in the case of interaction - induced dipole moment and polarizability is negligible .
what is more for the first - order hyperpolarizability bsse only has a noticable qualitative meaning , but the quantitative picture remains the same .
this observation allows for instance to save a lot of effort in a rational design of novel materials .
the earlier studies suggest the possibility of estimating the electric properties for the large molecular aggregates based on the accurate description of molecular properties and approximate descripton of the interaction - induced contributions [ 21 , 47 ] .
current findings allow to expect that calculations of molecular properties on high theory level ( mp2 or ccsd(t ) ) supplemented by the interaction - induced increments accounted for in dft could provide a relatively accurate and cheap alternative for regular post - hartree - fock predictions .
 a detailed analysis of the selected dft functionals for the calculations of interaction - induced dipole moment , polarizability and first - order hyperpolarizability has been carried out . 
 the hydrogen - bonded model chains consisting of hf , h2co and h3n molecules have been chosen as a case study . 
 the calculations of the components of the static electric properties using the diffuse dunning s basis set ( aug - cc - pvdz ) have been performed employing different types of density functionals ( b3lyp , lc - blyp , pbe0 , m06 - 2x and cam - b3lyp ) . 
 obtained results have been compared with those gained at the ccsd(t ) level of theory . 
 the counterpoise correction scheme , namely site - site function counterpoise , has been applied in order to eliminate basis set superposition error . 
 the performed tests allow to conclude that the dft functionals can provide a useful tool for prediction of the interaction - induced electric properties , however a caution has to be urged to their decomposition to the two- and many - body terms.figure  hydrogen - bonded model chains consisting of hf , h2co and h3n molecules 