The Scientific World Journal

Volume 2013 (2013), Article ID 832682, 9 pages

http://dx.doi.org/10.1155/2013/832682

Research Article

Static and Dynamic Electronic (Hyper)polarizabilities of Dimethylnaphthalene Isomers: Characterization of Spatial Contributions by Density Analysis

Andrea Alparone

Department of Chemistry, University of Catania, Viale A. Doria 6, 95125 Catania, Italy

Received 12 August 2013; Accepted 4 September 2013

Academic Editors: N. Lisitza and P. Pietrzyk

Copyright © 2013 Andrea Alparone. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Static and frequency-dependent electronic (hyper)polarizabilities of the dimethylnaphthalene (DMN) isomers were computed in vacuum using the Coulomb-attenuating Density Functional Theory method. The nonlinear optical Second Harmonic Generation (SHG) and Electro-Optical Pockels Effect (EOPE) were investigated at the characteristic Nd:YAG laser wavelength of 1064 nm. The response electric properties especially the longitudinal polarizability, polarizability anisotropy, and first-order hyperpolarizability are significantly affected by the position of the methyl groups. The SHG and EOPE techniques can be potentially useful to discriminate the ,-DMN isomers (2,6-DMN < 2,7-DMN < 2,3-DMN) as well as the ,-DMN isomers (1,5-DMN < 1,4-DMN < 1,8-DMN). The (hyper)polarizability differences among the investigated DMNs were elucidated through density analysis calculations. The predicted polarizabilities exhibit good linear relationships with the experimental first-order biomass-normalized rate coefficient, a physicochemical property connected to the rates of biodegradation processes of polycyclic aromatic hydrocarbons.

1. Introduction

Polycyclic aromatic hydrocarbons (PAHs) are very stable organic lipophilic pollutants principally produced during incomplete combustion processes, often exhibiting high levels of toxicity, mutagenicity, and carcinogenicity to humans and/or other living creatures [1–5]. Since PAHs are ubiquitous contaminants, there is much interest in searching practical strategies for their identification and removal from environments. Generally, alkylated derivatives are found in prevalence in the environmental samples of PAHs [6–8]. Dimethylnaphthalenes (DMNs) are substituted PAHs of great environmental interest [9–15]. DMNs can exist as ten isomeric forms (Figure 1): 1,2-DMN, 1,3-DMN, 1,4-DMN, 1,5-DMN, 1,6-DMN, 1,7-DMN, 1,8-DMN, 2,3-DMN, 2,6-DMN, and 2,7-DMN. It has been previously recognized that the enzymatic biodegradation of DMNs in aqueous media is strongly affected by the specific position of the methyl substituents [12]. Indeed the first-order biomass-normalized rate coefficient () of DMNs varies within one order of magnitude, the ,-disubstituted isomers showing the maximum values along the series [12]. Therefore, the detection of the DMN isomers in the environmental mixtures is of fundamental importance from environmental and biochemical viewpoints. On the basis of experimental and theoretical investigations, the active site of the enzyme which controls the biodegradative mechanism is mainly characterized by hydrophobic residues [12, 13], involving contributions from dispersive and/or inductive forces in enzyme-substrate complex formation. This result has been corroborated by recent computational studies on the electronic polarizabilities () [14, 16] and Raman spectra [15] of alkylated-naphthalenes, the average polarizabilities [14], and summation of the Raman activities [15] of DMNs being found to be linearly related to the experimental values. However, although the average polarizabilities () of DMNs regularly increase in the order (,-DMNs) < (,-DMNs) < (,-DMNs), they slightly vary among the ten isomers (up to 3.5%). Thus, the average polarizability is little useful to discriminate unambiguously the DMN isomers. A physicochemical property much more affected by the structural features is the electronic first-order hyperpolarizability () and the related nonlinear optical (NLO) properties Second Harmonic Generation (SHG) and Electro-Optical Pockels Effect (EOPE) [17–32]. SHG measurements are currently employed for the structural identification of biomolecules such as peptides and proteins [28–32].

Figure 1: Coordinate system and atom numbering of dimethylnaphthalene (DMN) isomers.

The main focus of this study is to determine the static and frequency-dependent electronic and values of the series of the DMN isomers, aiming to explore the effects of the position of the CH₃ groups on these electric properties, potentially helpful for the isomeric discrimination. The electronic (hyper)polarizabilities, are commonly predicted by means of ab initio and/or Density Functional Theory (DFT) computations. However, as well-known in the literature for an accurate determination of the electronic (hyper)polarizabilities, the choice of the functional is critical, especially in the case of -conjugated compounds [33]. In fact, on the whole, the conventional DFT methods tend to systematically overestimate the electronic (hyper)polarizabilities obtained by high-level correlated ab initio levels. On the other hand, the long-range corrected DFT methods incorporating nonlocal effects [34, 35], describe adequately the diffuse regions of the charge distributions, giving much more satisfactory performances for the prediction of the response electric properties. In the present study we used the Coulomb-attenuating hybrid exchange-correlation functional (CAM-B3LYP) [36], which has been recently employed with success for computing electronic (hyper)polarizabilities of organic compounds [37–48].

2. Computational Details

All calculations were performed with the Gaussian 09 program [49]. We used the molecular geometries previously optimized at the DFT-B3LYP level with the 6-31G* basis set [14]. Static and frequency-dependent electronic (hyper)polarizability tensor components and were obtained through the Coupled-Perturbed Hartree-Fock procedure [50, 51] using the DFT-CAM-B3LYP method with the polarized and diffuse 6-31+G*  basis set. The CAM-B3LYP functional and the 6-31+G*  basis set can be considered suitable choices especially for the prediction of the electronic (hyper)polarizabilities of organic molecules [37–48, 52–54]. However, we checked the performances of the basis set (6-31+G*  versus POL Sadlej’s basis set [55]) and level of calculation (CAM-B3LYP versus second-order Møller-Plesset perturbation theory (MP2)) on the related compound toluene. Dynamic (hyper)polarizabilities were evaluated at the experimental Nd:YAG laser wavelength of 1064 nm ( = 0.04282 a.u.) for the SHG and EOPE NLO processes.

We report the dipole moments (), the isotropically averaged polarizabilities (), the polarizability anisotropies (), and the isotropically first-order hyperpolarizabilities (), which are defined, respectively, as [56]

where is given by .

3. Results and Discussion

3.1. Basis Set and Level of Calculation Effects: Response Electric Properties of Toluene

For an accurate prediction of the electronic (hyper)polarizabilities, the choice of the basis set and level of computation is of crucial importance [57–67]. In the present study we explored the effects of the basis set and theoretical level on toluene as a test case. Table 1 reports the , , and values of toluene calculated using the CAM-B3LYP and MP2 levels with the 6-31+G* and Sadlej’s POL basis set. The latter basis set was specifically constructed for polarizability computations and has been recently employed with success to predict the electronic polarizabilities of naphthalene (N) [68]. However, it is well-demonstrated that for -conjugated compounds the smaller 6-31+G* basis set furnishes an adequate alternative to the POL as well as more extended basis sets for predicting response electric properties, but at significantly minor computational costs [39, 45–48, 52–54].

Table 1: Static electronic (), (), and (10⁻⁵³ C³m³J⁻²) of toluenea.

Δα

CAM-B3LYP/POL
12.25
6.68
279.3

CAM-B3LYP/6-31+G*
11.50
6.43
316.6

MP2/6-31+G*
11.43
6.24
312.4

All calculations are carried out on the B3LYP/6-31G* geometry.

The present results show that when passing from the 6-31+G* to the POL basis set, only marginal effects are observed. In fact, the and values increase by 0.75  (+6.5%) and 0.25  (+3.9%), respectively, whereas the decreases by 37.3 × 10⁻⁵³ C³m³J⁻² (−13.3%). Note that the (hyper)polarizability calculations carried out using the 6-31+G* basis set require noticeably minor CPU resources than those with the POL basis set (by about a factor of twenty!). In addition, we investigated the effects of the computational method, by comparing the CAM-B3LYP and MP2 (hyper)polarizability data. In line with the recent literature [37–39, 45–48, 69], the differences between the two levels are further smaller than those found for the basis sets, the , , and values being calculated within 0.07  (0.6%), 0.19  (3.0%), and 4.2 × 10⁻⁵³ C³m³J⁻² (1.3%), respectively. Thus, considering the above results, the CAM-B3LYP/6-31+G* level can be judged as an acceptable compromise between accuracy and computational cost and has been entirely employed for the subsequent calculations on the static and frequency-dependent (hyper)polarizabilities of the DMN isomers.

3.2. Static and Dynamic Polarizabilities of the DMN Isomers

Table 2 lists the static and frequency-dependent polarizabilities of the DMNs calculated in the gas phase at the CAM-B3LYP/6-31+G* level. For all the isomers, is the largest component, giving 43–49% of the total polarizabilities (). The dispersion effects here evaluated at the = 0.04282 a.u. are rather modest, increasing the static , and values by 0.54–0.70  + (2%), 0.34–0.36  (+2%) and 0.40–0.55  (+3%), respectively. Table 2 also reports the data of the unsubstituted compound N for which some experimental and high-level correlated ab initio values are available in the literature [68, 70]. The static CAM-B3LYP/6-31+G*  , and values of N agree satisfactorily with both the observed (within −0.8, −2.8, and +2.3%, resp.) [70] and CCSD/POL data (within −2.0, −4.0, and +2.4%, resp.) [68]. Not surprisingly, the static CAM-B3LYP/6-31+G* values of the DMN isomers underestimate the previously calculated B3LYP/6-31+G* figures [14] by 0.42–0.51  (2.0–2.4%), principally owing to the introduction of nonlocal effects.

Table 2: Static and dynamic  a.u.) electronic (), () and () of the dimethylnaphthalene isomers and naphthalenea.

0
0.04282
0
0.04282
0
0.04282

1,2-DMN
28.53
29.14
20.65
21.00
15.00
15.46

1,3-DMN
28.74
29.36
20.77
21.12
15.03
15.50

1,4-DMN
26.55
27.09
20.54
20.89
13.80
14.20

1,5-DMN
26.58
27.12
20.55
20.90
13.89
14.30

1,6-DMN
28.76
29.38
20.78
21.13
15.19
15.65

1,7-DMN
28.84
29.46
20.76
21.11
15.08
15.55

1,8-DMN
26.77
27.32
20.42
20.76
13.97
14.38

2,3-DMN
30.58
31.27
20.90
21.25
16.55
17.08

2,6-DMN
31.06
31.76
21.07
21.43
17.12
17.67

2,7-DMN
31.09
31.79
21.06
21.42
16.90
17.45

N
24.20(24.39)c
16.91(17.40)c
13.21(12.91)c

All calculations were carried out at the CAM-B3LYP/6-31+G* level on the B3LYP/6-31G* geometry taken from [14].

Experimental first-order biomass-normalized rate coefficient ((mg of protein/L)⁻¹(h)⁻¹) in aqueous systems, [12].

[70].

The order of the static and dynamic CAM-B3LYP/6-31+G* values is the following:

 1,4-DMN ~ 1,5-DMN < 1,8-DMN < 1,2-DMN < 1,3-DMN ~ 1,6-DMN < 1,7-DMN < 2,3-DMN < 2,6-DMN ~ 2,7-DMN.

For the values, the above order is slightly modified, with the 1,8-DMN and 2,6-DMN isomers being, respectively, the less and more polarizable along the series:

 1,8-DMN < 1,4-DMN ~ 1,5-DMN < 1,2-DMN < 1,7-DMN ~ 1,3-DMN ~ 1,6-DMN < 2,3-DMN < 2,7-DMN ~ 2,6-DMN.

The order of the values is rather similar to that found for the data, except for the inversions between 1,6-DMN and 1,7-DMN and between 2,6-DMN and 2,7-DMN:

 1,4-DMN < 1,5-DMN < 1,8-DMN < 1,2-DMN ~ 1,3-DMN ~ 1,7-DMN < 1,6-DMN 2,3-DMN < 2,7-DMN < 2,6-DMN.

All the predicted polarizabilities concordantly increase on passing from the ,-DMNs to the ,-DMNs and then to the ,-DMNs, in agreement with the previous data obtained at the B3LYP/6-31+G* level in gaseous and aqueous phases [14]. Specifically, when we compare the 1,4-DMN and 2,6-DMN isomers, the static CAM-B3LYP/6-31+G* , and values enhance by 4.51  (+17.0%), 0.53  (+2.6%), and 3.32  (+24.1%), respectively. However, whereas the static data slightly change along the series of isomers (within 0.65 , 3.2%), the and values are distributed over larger ranges being within 3.32  (24.1%) and 4.54  (17.1%), respectively. On the whole, these results suggest that, in comparison to the average polarizabilities, the and properties are much more affected by the position of the methyl substituent, being potentially useful to identify the DMN isomers. Additionally, in agreement with a previous study on the average polarizabilities [14], the present and data of DMNs are found to be linearly related to the biodegradation experimental biomass-normalized first-order rate coefficients [12], confirming the crucial role of the polarizabilities in the biodegradation process of this group of organic pollutants. The , /, and linear relationships are displayed in Figure 2 showing good statistics since the value is predicted between 0.97 and 1.00 (following the discussion reported in [14], the 2,7-DMN isomer was excluded from the relationships).

Figure 2: Relationships between the experimental biomass-normalized first-order rate coefficient [12] and the gas phase CAM-B3LYP/6-31+G* polarizabilities of the DMN isomers. (a)   ; (b)   ; (c)   .

In order to explain the polarizability differences among the DMN isomers, we analyzed the spatial contributions of electrons to the polarizabilities by using the concept of density of polarizability [71, 72]. The is defined as derivative of the charge density function with respect to the applied fields ( is the position vector). The is commonly expanded in powers of as

For a certain positive-negative pair, the sign is positive when the direction of the positive to negative densities coincides with the positive direction of the chosen coordinate system, whereas the magnitude is proportional to the distance between the two densities. In the present study, we determined the densities at the CAM-B3LYP/6-31+G* level for the longitudinal component . As representative cases of the ,-DMN and ,-DMN isomers, we investigated the densities for 1,4-DMN and 2,3-DMN. The distributions are displayed in Figure 3. From the CAM-B3LYP/6-31+G* computations, the (2,3-DMN) value is higher than the corresponding datum of the 1,4-DMN isomer by about 4  (Table 2). For both the isomers, the largest positive to negative contributions are mainly placed over the N moiety. However, as can be appreciated from the plots in Figure 3, the 2,3-DMN isomer exhibits additional relevant densities localized on both the CH₃ groups. These methylic amplitudes are much more conspicuous than those found for the 1,4-DMN isomer, contributing to the increase of the component by ca. 17%.

Figure 3: Calculated density distributions for the 1,4-DMN (top) and 2,3-DMN (bottom) isomers. Positive and negative are represented by yellow and blue isosurfaces (0.001 a.u.), respectively. CAM-B3LYP/6-31+G* results.

3.3. Static and Dynamic First-Order Hyperpolarizabilities of the DMN Isomers

The calculated static and values of the DMN isomers obtained in vacuum at the CAM-B3LYP/6-31+G* level are collected in Table 3. As for the computed dipole moments [14], due to the mutual disposition of the methyl substituents, the 1,5-DMN and 2,6-DMN isomers are nonpolar compounds. In the present study, besides to the static first-order hyperpolarizabilities we also determined the frequency-dependent properties for the SHG and EOPE NLO phenomena since observed data are nearly always obtained at incident optical fields. In order to minimize resonance enhancements, the dynamic hyperpolarizabilities were evaluated at the value of 1064 nm ( = 0.04282 a.u.), which is sufficiently apart from the observed lowest-energy absorption of DMNs (the experimental values are in the range 274–289 nm) [73]. The dynamic data are included in Table 3. As should be expected, for all the isomers. The dispersion effects enhance the static values by 6–18% and 17–57%, for the EOPE and SHG processes, respectively, the largest increases being found for 1,4-DMN, in agreement with its highest value among the DMN isomers.

Table 3: Calculated , static and dynamic  a.u.) electronic and (10⁻⁵³ C³m³J⁻²) of the dimethylnaphthalene isomersa.

1,2-DMN
0.69
430.3
743.9
808.0
923.4

1,3-DMN
0.54
−412.0
599.6
660.5
782.3

1,4-DMN
0.06
62.5
179.6
211.6
282.2

1,5-DMN
0.00
0.0
0.0
0.0
0.0

1,6-DMN
0.49
−348.8
448.9
474.5
529.0

1,7-DMN
0.62
370.3
628.4
670.1
737.5

1,8-DMN
0.66
0.0
673.3
734.2
859.3

2,3-DMN
0.90
879.8
1157.5
1256.9
1474.9

2,6-DMN
0.00
0.0
0.0
0.0
0.0

2,7-DMN
0.20
0.0
461.7
497.0
541.9

All calculations were carried out at the CAM-B3LYP/6-31+G* level on the B3LYP/6-31G* geometry.

The static and dynamic values concordantly predict the following order:

 1,4-DMN < 1,6-DMN < 2,7-DMN < 1,3-DMN < 1,7-DMN < 1,8-DMN < 1,2-DMN < 2,3-DMN.

Note that the above order is rather different from those obtained for the polarizabilities; a versus relationship cannot be established. In particular, the static and dynamic CAM-B3LYP/6-31+G*  (2,3-DMN) values are calculated to be 5-6 times higher than the corresponding data obtained for the 1,4-DMN isomer. More importantly, we notice that the first-order hyperpolarizabilities of the ,-DMNs are somewhat different from each other for the 2,6-DMN isomer the is zero since it belongs to the symmetry point group, whereas the (2,3-DMN)/(2,7-DMN) ratios are predicted to be 2.5–2.7 by the present static and dynamic computations. Similarly, the values of the ,-DMN isomers are rather different from each other. In fact, on passing from 1,4-DMN to 1,8-DMN, the data increase by a factor of three/four, whereas the (1,5-DMN) value is zero ( symmetry point group). Note that a somewhat different situation occurs for the ,-DMN isomers, which exhibit hyperpolarizabilities much closer to each other.

Interestingly, as for the computed data, and (Table 3) are, respectively, the smallest and greatest values along the series of the DMNs. The order of the values is roughly similar to that of the first-order hyperpolarizabilities, a linear relationship between the and values being established (, ). This result indicates that, as for the data, the values of the DMNs are mainly affected by the mutual disposition of the groups, mesomeric effect enhancements of the hyperpolarizabilities being expected to be marginal.

As for the polarizabilities, we explored the hyperpolarizability differences between the 1,4-DMN and 2,3-DMN isomers through hyperpolarizability density analyses. The density, , is expressed as follows [71, 74]:

By analyzing the main contributing components, for the 1,4-DMN isomer, the largest hyperpolarizability component lies along the -axis. At the CAM-B3LYP/6-31+G*  level, the static (1,4-DMN) value is calculated to be 62.5 × 10⁻⁵³ C³m³J⁻², recovering ca. 35% of the value. Similarly, the component dominates the first-order hyperpolarizability of 2,3-DMN, the ratio being calculated to be 0.76. Thus, we computed the distributions for the component at the CAM-B3LYP/6-31+G*  level using a numerical procedure described in detail in [74]. Figure 4 displays the densities for the 1,4-DMN and 2,3-DMN isomers. As can be seen from the graphical representations, for both the isomers, the most significant amplitude is furnished by the N skeleton, although it is much more expanded in the case of the 2,3-DMN isomer. In addition, different from the 1,4-DMN isomer, for 2,3-DMN, a conspicuous positive contribution is also provided by both the CH₃ groups. On the whole, the above results are in qualitative agreement with the calculated static hyperpolarizabilities, the (2,3-DMN)/(1,4-DMN) and (2,3-DMN)/(1,4-DMN) ratios being predicted to be 14 and 6, respectively.

Figure 4: Calculated density distributions for the 1,4-DMN (top) and 2,3-DMN (bottom) isomers. Positive and negative are represented by yellow and blue isosurfaces (0.001 a.u.), respectively. CAM-B3LYP/6-31+G* results.

4. Conclusions

In this study, we investigated the electronic polarizabilities and first-order hyperpolarizabilities of the series of DMN isomers. The computations were performed in vacuum using the CAM-B3LYP functional and the 6-31+G*  basis set. The response electric properties were obtained in the static and dynamic regimes for the SHG and EOPE NLO phenomena at the value of 1064 nm. The average polarizability varies a little along the series of isomers, whereas both the longitudinal polarizability and anisotropy of polarizability are much more affected by the position of the methyl substituents. In agreement with recent theoretical studies on the polarizabilities [14] and Raman spectra [15], linear relationships are established between the calculated polarizabilities and the experimental biodegradation rates of DMNs, confirming the important role of dispersive and/or inductive interactions for the biodegradative mechanisms of this group of substituted PAH isomers.

The static and frequency-dependent first-order hyperpolarizabilities are strongly dependent on the relative position of the groups. This is especially evident for the ,-DMN and ,-DMN isomers: 2,6-DMN < 2,7-DMN < 2,3-DMN and 1,5-DMN < 1,4-DMN < 1,8-DMN. The present results suggest that some DMN isomers might be distinguished through SHG and EOPE NLO measurements. Density analysis computations are useful for qualitative interpretations of the (hyper)polarizabilities of the DMN isomers.

References

J. C. Arcos and M. F. Argus, “Molecular geometry and carcinogenic activity of aromatic compounds. New perspectives,” Advances in Cancer Research, vol. 11, pp. 305–471, 1969.

D. Hoffmann and E. L. Wynder, “Organic particulate pollutants—chemical analysis and bioassays for carcinogenicity,” in Air Pollution, A. C. Stern, Ed., vol. 2, pp. 361–455, Academic Press, New York, NY, USA, 1977.

H. V. Gelboin and P. O. P. Ts’o, Polycyclic Hydrocarbons and Cancer, Academic Press, New York, NY, USA, 1978.

R. G. Harvey, Polycyclic Aromatic Hydrocarbons: Chemistry and Carcinogenicity, Cambridge University Press, Cambridge, UK, 1991.

B. J. Finlayson-Pitts and J. N. Pitts, Chemistry of the Upper and Lower Atmosphere, Academic Press, San Diego, Calif, USA, 2000.

C. A. Peters, “Coal tar dissolution in water-misclble solvents: experimental evaluation,” Environment Science Technology, vol. 27, no. 13, pp. 2831–2843, 1993.

K. L. Garcia, J. J. Delfino, and D. H. Powell, “Non-regulated organic compounds in Florida sediments,” Water Research, vol. 27, no. 11, pp. 1601–1613, 1993.

L. P. Burkhard and B. R. Sheedy, “Evaluation of screening procedures for bioconcentratable organic chemicals in effluents and sediments,” Environmental Toxicology and Chemistry, vol. 14, no. 4, pp. 697–711, 1995.

R. Atkinson and S. M. Aschmann, “Kinetics of the reactions of naphthalene, 2-methylnaphthalene, and 2,3-dimethylnaphthalene with OH radicals and with O₃ at 295 1 K,” International Journal of Chemical Kinetics, vol. 18, pp. 569–573, 1986.

C. E. Banceu, C. Mihele, D. A. Lane, and N. J. Bunce, “Reactions of methylated naphthalenes with hydroxyl radicals under simulated atmospheric conditions,” Polycyclic Aromatic Compounds, vol. 18, no. 4, pp. 415–425, 2001.

P. T. Phousongphouang and J. Arey, “Rate constants for the gas-phase reactions of a series of alkylnaphthalenes with the OH radical,” Environmental Science and Technology, vol. 36, no. 9, pp. 1947–1952, 2002.

K. H. Wammer and C. A. Peters, “Polycyclic aromatic hydrocarbon biodegradation rates: A Structure-based Study,” Environmental Science and Technology, vol. 39, no. 8, pp. 2571–2578, 2005.

K. H. Wammer and C. A. Peters, “A molecular modeling analysis of polycyclic aromatic hydrocarbon biodegradation by naphthalene dioxygenase,” Environmental Toxicology and Chemistry, vol. 25, no. 4, pp. 912–920, 2006.

V. Librando and A. Alparone, “Electronic polarizability as a predictor of biodegradation rates of dimethylnaphthalenes. An ab initio and density functional theory study,” Environmental Science and Technology, vol. 41, no. 5, pp. 1646–1652, 2007.

A. Alparone and V. Librando, “Raman DFT study of dimethylnaphthalenes: isomer identification and prediction of biodegradation rate coefficients,” Structural Chemistry, vol. 23, pp. 1467–1474, 2012.

V. Librando and A. Alparone, “Structure, vibrational properties and polarizabilities of methylnaphthalene isomers. A quantum-mechanical approach,” Polycyclic Aromatic Compounds, vol. 27, no. 1, pp. 65–94, 2007.

J. S. Salafsky, “Second-harmonic generation as a probe of conformational change in molecules,” Chemical Physics Letters, vol. 381, no. 5-6, pp. 705–709, 2003.

M. Sliwa, K. Nakatani, T. Asahi, P. G. Lacroix, R. B. Pansu, and H. Masuhara, “Polarization and wavelength dependent nonlinear optical properties of a photo-switchable organic crystal,” Chemical Physics Letters, vol. 437, no. 4–6, pp. 212–217, 2007.

M. Sliwa, S. Létard, I. Malfant et al., “Design, synthesis, structural and nonlinear optical properties of photochromic crystals: toward reversible molecular switches,” Chemistry of Materials, vol. 17, no. 18, pp. 4727–4735, 2005.

A. Alparone, “Dipole (hyper)polarizabilities of fluorinated benzenes: an ab initio investigation,” Journal of Fluorine Chemistry, vol. 144, pp. 94–101, 2012.

G. J. M. Velders, J.-M. Gillet, P. J. Becker, and D. Feil, “Electron density analysis of nonlinear optical materials. An ab initio study of different conformations of benzene derivatives,” Journal of Physical Chemistry, vol. 95, no. 22, pp. 8601–8608, 1991.

E. Hendrickx, K. Clays, A. Persoons, C. Dehu, and J. L. Brédas, “The bacteriorhodopsin chromophore retinal and derivatives: an experimental and theoretical investigation of the second-order optical properties,” Journal of the American Chemical Society, vol. 117, no. 12, pp. 3547–3555, 1995.

A. Alparone, A. Millefiori, and S. Millefiori, “Electronic dipole polarizability and hyperpolarizability of formic acid,” Chemical Physics Letters, vol. 409, no. 4-6, pp. 288–294, 2005.

W. Niewodniczański and W. Bartkowiak, “Theoretical study of geometrical and nonlinear optical properties of pyridinum N-phenolate betaine dyes,” Journal of Molecular Modeling, vol. 13, no. 6-7, pp. 793–800, 2007.

Z.-M. Su, H.-L. Xu, Z.-R. Li, S. Muhammad, F. L. Gu, and K. Harigaya, “Knot-isomers of möbius cyclacene: How does the number of knots influence the structure and first hyperpolarizability,” Journal of Physical Chemistry C, vol. 113, no. 34, pp. 15380–15383, 2009.

A. Plaquet, B. Champagne, F. Castet et al., “Theoretical investigation of the dynamic first hyperpolarizability of DHA-VHF molecular switches,” New Journal of Chemistry, vol. 33, no. 6, pp. 1349–1356, 2009.

J. Lipiński and W. Bartkowiak, “Conformation and solvent dependence of the first and second molecular hyperpolarizabilities of charge-transfer chromophores. Quantum-chemical calculations,” Chemical Physics, vol. 245, no. 1-3, pp. 263–276, 1999.

J. Perrenoud-Rinuy, P.-F. Brevet, and H. H. Girault, “Second harmonic generation study of myoglobin and hemoglobin and their protoporphyrin IX chromophore at the water/1,2-dichloroethane interface,” Physical Chemistry Chemical Physics, vol. 4, no. 19, pp. 4774–4781, 2002.

S. A. Mitchell and R. A. McAloney, “Second harmonic optical activity of tryptophan derivatives adsorbed at the air/water interface,” Journal of Physical Chemistry B, vol. 108, no. 3, pp. 1020–1029, 2004.

J. S. Salafsky, “Second-harmonic generation for studying structural motion of biological molecules in real time and space,” Physical Chemistry Chemical Physics, vol. 9, no. 42, pp. 5704–5711, 2007.

R. M. Williams, W. R. Zipfel, and W. W. Webb, “Interpreting second-harmonic generation images of collagen I fibrils,” Biophysical Journal, vol. 88, no. 2, pp. 1377–1386, 2005.

W.-L. Chen, T.-H. Li, P.-J. Su et al., “Second harmonic generation χ tensor microscopy for tissue imaging,” Applied Physics Letters, vol. 94, no. 18, Article ID 183902, 3 pages, 2009.

B. Champagne, E. A. Perpète, S. J. A. Van Gisbergen et al., “Assessment of conventional density functional schemes for computing the polarizabilities and hyperpolarizabilities of conjugated oligomers: an ab initio investigation of polyacetylene chains,” Journal of Chemical Physics, vol. 109, no. 23, pp. 10489–10498, 1998.

M. Kamiya, H. Sekino, T. Tsuneda, and K. Hirao, “Nonlinear optical property calculations by the long-range-corrected coupled-perturbed Kohn-Sham method,” Journal of Chemical Physics, vol. 122, no. 23, Article ID 234111, 10 pages, 2005.

H. Sekino, Y. Maeda, M. Kamiya, and K. Hirao, “Polarizability and second hyperpolarizability evaluation of long molecules by the density functional theory with long-range correction,” Journal of Chemical Physics, vol. 126, no. 1, Article ID 014107, 2007.

T. Yanai, D. P. Tew, and N. C. Handy, “A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP),” Chemical Physics Letters, vol. 393, no. 1–3, pp. 51–57, 2004.

A. Alparone, “Linear and nonlinear optical properties of nucleic acid bases,” Chemical Physics, vol. 410, pp. 90–98, 2013.

D. Jacquemin, E. A. Perpète, G. Scalmani, M. J. Frisch, R. Kobayashi, and C. Adamo, “Assessment of the efficiency of long-range corrected functionals for some properties of large compounds,” Journal of Chemical Physics, vol. 126, no. 14, Article ID 144105, 12 pages, 2007.

P. A. Limacher, K. V. Mikkelsen, and H. P. Lüthi, “On the accurate calculation of polarizabilities and second hyperpolarizabilities of polyacetylene oligomer chains using the CAM-B3LYP density functional,” Journal of Chemical Physics, vol. 130, no. 19, Article ID 194114, 7 pages, 2009.

M. Nakano, T. Minami, H. Fukui et al., “Approximate spin-projected spin-unrestricted density functional theory method: application to the diradical character dependences of the (hyper)polarizabilities in p-quinodimethane models,” Chemical Physics Letters, vol. 501, no. 1-3, pp. 140–145, 2010.

A. Alparone, “Comparative study of CCSD(T) and DFT methods: electronic (hyper)polarizabilities of glycine,” Chemical Physics Letters, vol. 514, no. 1–3, pp. 21–25, 2011.

C.-C. Zhang, H.-L. Xu, Y.-Y. Hu, S.-L. Sun, and Z.-M. Su, “Quantum chemical research on structures, linear and nonlinear optical properties of the Li@n-acenes salt ( = 1, 2, 3, and 4),” Journal of Physical Chemistry A, vol. 115, no. 10, pp. 2035–2040, 2011.

T. Pluta, M. Kolaski, M. Medved’, and S. Budzák, “Dipole moment and polarizability of the low-lying excited states of uracil,” Chemical Physics Letters, vol. 546, pp. 24–29, 2012.

M. Medved, Š. Budzák, and T. Pluta, “Static NLO responses of fluorinated polyacetylene chains evaluated with long-range corrected density functionals,” Chemical Physics Letters, vol. 515, no. 1–3, pp. 78–84, 2011.

A. Alparone, “Evolution of electric dipole (hyper)polarizabilities of b-strand polyglycine single chains: an ab initio and DFT theoretical study,” Journal of Physical Chemistry A, vol. 117, pp. 5184–5194, 2013.

A. Alparone, “The effect of secondary structures on the NLO properties of single chain oligopeptides: a comparison between -strand and α-helix polyglycines,” Physical Chemistry Chemical Physics, vol. 15, pp. 12958–12962, 2013.

A. Alparone, “Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers,” Journal of Molecular Modeling, vol. 19, pp. 3095–3102, 2013.

A. Alparone, “Response electric properties of α-helix polyglycines: a CAM-B3LYP DFT investigation,” Chemical Physics Letters, vol. 563, pp. 88–92, 2013.

M. J. Frisch, G. W. Trucks, H. B. Schlegel, et al., “Gaussian 09, Revision A.02,” Gaussian, Wallingford, Conn, USA, 2009

H. Sekino and R. J. Bartlett, “Frequency dependent nonlinear optical properties of molecules,” The Journal of Chemical Physics, vol. 85, no. 2, pp. 976–989, 1986.

J. E. Rice and N. C. Handy, “The calculation of frequency-dependent hyperpolarizabilities including electron correlation-effects,” International Journal of Quantum Chemistry, vol. 43, pp. 91–118, 1992.

V. Librando, A. Alparone, and Z. Minniti, “Computational study on dipole moment, polarizability and second hyperpolarizability of nitronaphthalenes,” Journal of Molecular Structure, vol. 856, no. 1–3, pp. 105–111, 2008.

M. Torrent-Sucarrat, M. Solà, M. Duran, J. M. Luis, and B. Kirtman, “Basis set and electron correlation effects on ab initio electronic and vibrational nonlinear optical properties of conjugated organic molecules,” Journal of Chemical Physics, vol. 118, no. 2, pp. 711–718, 2003.

A. Alparone and V. Librando, “Physicochemical characterization of environmental mutagens: 3-nitro-6-azabenzo[a]pyrene and its N-oxide derivative,” Monatshefte für Chemie, vol. 143, pp. 1123–1132, 2012.

I. Černušák, V. Kello, and A. J. Sadlej, “Standardized medium-size basis sets for calculations of molecular electric properties: group IIIA,” Collection of Czechoslovak Chemical Communications, vol. 68, no. 2, pp. 211–239, 2003.

V. Keshari, W. M. K. P. Wijekoon, P. N. Prasad, and S. P. Karna, “Hyperpolarizabilities of organic molecules: Ab initio time-dependent coupled perturbed Hartree—Fock—Roothaan studies of basic heterocyclic structures,” Journal of Physical Chemistry, vol. 99, no. 22, pp. 9045–9050, 1995.

J. E. Rice, R. D. Amos, S. M. Colwell, N. C. Handy, and J. Sanz, “Frequency dependent hyperpolarizabilities with application to formaldehyde and methyl fluoride,” The Journal of Chemical Physics, vol. 93, no. 12, pp. 8828–8839, 1990.

H. Sekino and R. J. Bartlett, “Molecular hyperpolarizabilities,” Journal of Chemical Physics, vol. 98, no. 4, pp. 3022–3037, 1993.

F. Sim, S. Chin, M. Dupuis, and J. E. Rice, “Electron correlation effects in hyperpolarizabilities of p-nitroaniline,” Journal of Physical Chemistry, vol. 97, no. 6, pp. 1158–1163, 1993.

G. Maroulis and A. J. Thakkar, “Polarizabillties and hyperpolarizabilities of F2,” The Journal of Chemical Physics, vol. 90, no. 1, pp. 366–370, 1989.

E. F. Archibong and A. J. Thakkar, “Hyperpolarizabilities and polarizabilities of Li-1 and B+: finite-field coupled-cluster calculations,” Chemical Physics Letters, vol. 173, no. 5-6, pp. 579–584, 1990.

E. F. Archibong and A. J. Thakkar, “Static polarizabilities and hyperpolarizabilities, and multipole moments for Cl2 and Br2. Electron correlation and molecular vibration effects,” Chemical Physics Letters, vol. 201, no. 5-6, pp. 485–492, 1993.

G. Maroulis, “Hyperpolarizability of H₂O revisited: accurate estimate of the basis set limit and the size of electron correlation effects,” Chemical Physics Letters, vol. 289, no. 3-4, pp. 403–411, 1998.

D. Xenides and G. Maroulis, “Basis set and electron correlation effects on the first and second static hyperpolarizability of SO₂,” Chemical Physics Letters, vol. 319, no. 5-6, pp. 618–624, 2000.

H. A. Kurtz, J. J. P. Stewart, and K. M. Dieter, “Calculation of the nonlinear optical properties of molecules,” Journal of Computational Chemistry, vol. 11, pp. 82–87, 1990.

D. M. Bishop and G. Maroulis, “Accurate prediction of static polarizabilities and hyperpolarizabilities. A study on FH (X₁Σ⁺),” The Journal of Chemical Physics, vol. 82, no. 5, pp. 2380–2391, 1985.

J. E. Rice and N. C. Handy, “The calculation of frequency-dependent polarizabilities as pseudo-energy derivatives,” The Journal of Chemical Physics, vol. 94, no. 7, pp. 4959–4971, 1991.

J. R. Hammond, K. Kowalski, and W. A. Dejong, “Dynamic polarizabilities of polyaromatic hydrocarbons using coupled-cluster linear response theory,” Journal of Chemical Physics, vol. 127, no. 14, Article ID 144105, 9 pages, 2007.

A. Alparone, “Structural, energetic and response electric properties of cyclic selenium clusters: an ab initio and density functional theory study,” Theoretical Chemistry Accounts, vol. 131, article 1239, pp. 1–14, 2012.

M. F. Vuks, “Determination of optical anisotropy of aromatic molecules from double refraction in crystals,” Optics Spectroscopy, vol. 20, pp. 361–368, 1966.

M. Nakano, I. Shigemoto, S. Yamada, and K. Yamaguchi, “Size-consistent approach and density analysis of hyperpolarizability: second hyperpolarizabilities of polymeric systems with and without defects,” The Journal of Chemical Physics, vol. 103, no. 10, pp. 4175–4191, 1995.

M. Nakano, S. Ohta, K. Tokushima et al., “First and second hyperpolarizabilities of donor-acceptor disubstituted diphenalenyl radical systems,” Chemical Physics Letters, vol. 443, no. 1–3, pp. 95–101, 2007.

N. Mora-Diez, R. J. Boyd, and G. L. Heard, “Effects of alkyl substituents on the excited states of naphthalene: Semiempirical Study,” Journal of Physical Chemistry A, vol. 104, no. 5, pp. 1020–1029, 2000.

S. Yamada, M. Nakano, I. Shigemoto, S. Kiribayashi, and K. Yamaguchi, “Intense electron correlation dependence of the first hyperpolarizabilities β of a nitroxide radical and formaldehyde,” Chemical Physics Letters, vol. 267, no. 5-6, pp. 445–451, 1997.