Electronic (βe) and vibrational (βv) first-order hyperpolarizabilities of uracil were determined in gas and water solution using the Coulomb-attenuating Density Functional Theory level with the Dunning’s correlation-consistent aug-cc-pVDZ basis set. Frequency-dependent βe values were computed for the Second Harmonic Generation (SHG) and Electric Optical Pockels Effect (EOPE) nonlinear optical phenomena. The Polarized Continuum Model was employed to study the solvent effects on the electronic and vibrational properties. The introduction of solvation contributions increases the βe(static) value by ca. 110%. In comparison, smaller enhancements are found for the βe(EOPE) and βe(SHG) data evaluated at the typical wavelength of 694 nm (by 40–50%). The gas-water hyperpolarizability difference was rationalised through a density analysis study. The magnitudes of the vibrational first-order hyperpolarizabilities are comparable to their electronic counterparts and noticeably increase in solution: βv(EOPE) ~ βe(EOPE) in aqueous phase at λ=694 nm. Analysis of the IR and Raman spectra is useful to elucidate the most important contributing modes to the vibrational first-order hyperpolarizabilities.
1. Introduction
Organic nonlinear optical (NLO) compounds are intensively studied, primarily for their potential use in the design of photonic and optoelectronic devices [1–3]. Biomolecules are attractive NLO materials, having the practical advantage to be already available in nature. Over recent years, DNA-based systems have received great attention for their conductive and NLO applications [4–12]. Nevertheless, characterization of the NLO properties of single nucleic acid bases is still rather incomplete. To the best of our knowledge, experimental response electric properties of the smallest base uracil are not available so far, whereas some theoretical estimates of the electronic polarizabilities (αe) [13–21] and second-order hyperpolarizabilities (γe) [18] were previously reported. However, there is significant interest in exploring the second-order NLO effects, which are important for immediate practical applications. At the microscopic level, the second-order NLO properties are associated with the first-order hyperpolarizability tensor (βijk), which originates from the responses of a molecular system to external electric field strengths Fi, producing an induced dipole moment μi(Fi):
(1)μi(Fi)=μi(0)+αijFj+12!βijkFjFk+13!γijklFjFkFl+⋯.
Recently, pure vibrational contributions to the first-order hyperpolarizability of uracil have been calculated in vacuum through a Lanczos procedure [22], whereas explorations of solvent effects on the electronic and vibrational β values are still lacking to date.
Our current computational study mainly focuses on the electronic (βe) and vibrational (βv) static and dynamic first-order hyperpolarizabilities of uracil. The fundamental role of the vibrational counterparts to the hyperpolarizabilities has been widely documented [23]. The present calculations were performed in gas and water solution under the Polarized Continuum Model (PCM) approximation [24, 25]. There are many indications in the literature showing that calculated first-order hyperpolarizabilities of organic molecules are strongly affected by solvent contributions [26–33]. Solvation effects on the electronic and vibrational α [18, 19] and γ [18] values of uracil have been previously explored by means of PCM Hartree-Fock and DFT computations in carbon tetrachloride, acetonitrile, and water solutions.
2. Computational Methods
The present calculations were performed in the gas phase and water solution (ε=78.3553) with the Gaussian 09 package [34]. The solvent effects were entirely modelled under the PCM approximation as implemented in the Gaussian 09 program. The geometry of uracil was optimized under the planar Cs symmetry using the CAM-B3LYP functional [35] and the polarised and diffuse Dunning’s correlation-consistent aug-cc-pVDZ basis set [36]. The IR and Raman spectra were simulated under the harmonic approximation at the CAM-B3LYP/aug-cc-pVDZ level on the geometries optimized at the same level. The structures are true minima on the potential energy surfaces (no imaginary wavenumbers).
Static βe values were calculated at the CAM-B3LYP/aug-cc-pVDZ//CAM-B3LYP/aug-cc-pVDZ level. We selected the CAM-B3LYP functional and aug-cc-pVDZ basis set considering their satisfactory performances in the prediction of the response electric properties of organic compounds, reproducing adequately first-order hyperpolarizabilities obtained using high-level correlated ab initio methods and larger basis sets [28, 29, 37–45]. The dynamic electronic first-order hyperpolarizabilities [βe(-ωσ;ω1,ω2)] for the Second Harmonic Generation [SHG, βe(-2ω;ω,ω)] and Electric Optical Pockels Effect [EOPE, βe(-ω;ω,0)] NLO phenomena were calculated at the CAM-B3LYP/aug-cc-pVDZ level in the ħω range 0–0.06563 a.u. The highest ħω value corresponds to the wavelength (λ) of 694 nm, which is characteristic of the ruby laser.
Static pure vibrational first-order hyperpolarizabilities were obtained at the CAM-B3LYP/aug-cc-pVDZ//CAM-B3LYP/aug-cc-pVDZ level in vacuum and water solution under the double-harmonic approximation (the used symbols have their standard meaning) [23]:
(2)βijkv=[μα]0,0=∑a3N-6((∂μi∂Qa)0(∂αjk∂Qa)0+(∂μj∂Qa)0(∂αik∂Qa)0=wwiww+(∂μk∂Qa)0(∂αij∂Qa)0)×(ωa2)-1.
By assuming the validity of the infinity frequency approximation [46], the βv(SHG) and βv(EOPE) processes are, respectively,
(3)βijkv(-2ω;ω,ω)ω→∞=0,βijkv(-ω;ω,0)ω→∞=13[μα]0,0.
In this study we report the invariant first-order hyperpolarizabilities (βvec) [47]:
(4)βvec=βx2+βy2+βz2,
where βi (i=x,z) is given by βi=(1/3)∑j=x,y,z(βijj+βjij+βjji).
Atomic units are used throughout the work. Conversion factor to the SI is: 1 a.u. of β(e3a03Eh-2)=3.206361×10-53C3m3J-2.
3. Results and Discussion
Table 1 lists the CAM-B3LYP/aug-cc-pVDZ dipole moments. The largest μ component lies along the z-axis, recovering ca. 96% of the total μ value. The gas phase μ(CAM-B3LYP/aug-cc-pVDZ) of 4.57 D overestimates by 18% the experimental datum obtained by microwave measurements [μ(exp.) = 3.87 D] [48], being in good agreement with the high-level ab initio CCSD(T)/aug-cc-pVDZ estimate of 4.33 D (+5.5%) [19]. The introduction of the solvation contributions increases the μ value by 1.4 D (+37%), in qualitative consistency with the observed increase of 0.26 D when passing from the gas phase [48] to dioxane solution [49].
Dipole moments μ (D) and static electronic first-order hyperpolarizabilities βe (a.u.) of uracila.
Gas
Water
μx
1.21
1.88
μz
4.41
5.96
μ
4.57 (3.87)b
6.25
βxxxe
79.3
183.9
βxyye
19.5
44.0
βxzze
5.7
24.6
βzxxe
−106.5
−240.5
βzyye
−36.5
−56.4
βzzze
78.2
262.3
βxe
104.6
252.5
βze
−64.8
−34.6
βvece
123.0
254.9
aCalculations were carried out at the CAM-B3LYP/aug-cc-pVDZ level on the geometry calculated at the same level.
bReference [48].
Table 1 also includes the static electronic first-order hyperpolarizability tensor components βijje (i=x,z; j=x,y,z) in gaseous and aqueous phases. In gas, βzxxe is in absolute value the predominant component (−106.5 a.u.), whereas in water solution the largest components are βzxxe(-240.5 a.u.) and βzzze(262.3 a.u.). When passing from the gas phase to the water solution, the βxe value increases by about a factor of two, whereas on the contrary |βze| decreases by ca. a factor of two. From the present computations, βxe dominates the first-order hyperpolarizability of both the gaseous and aqueous phases, giving ca. 85% and 99% of the βvece value, respectively.
In order to clarify the solvation effects on the response electric properties, we determined the spatial contributions of electrons to the first-order hyperpolarizabilities by computing density of hyperpolarizability amplitudes, ρjk(2)(r) [50, 51]. The ρjk(2)(r) is defined as derivative of the charge density function ρ(r,F) with respect to applied electric field strengths F (r is the position vector). The ρ(r,F) is usually expanded in powers of F:
(5)ρ(r,F)=ρ(0)(r)+∑jρj(1)(r)Fj+12!∑jρjk(2)(r)FjFk+12!∑jρjkl(3)(r)FjFkFl+⋯,ρjk(2)(r)=∂2ρ(r,F)∂Fj∂Fk|Fj=0,Fk=0,βijke=-12!∫rρjk(2)(r)dr.
For a certain positive-negative ρjk(2)(r) pair, the sign is positive when the direction of the positive to negative density is coincident with the positive direction of the chosen coordinate system (Figure 1), whereas the magnitude is proportional to the distance between the two densities. Following the current calculations, the main contribution to βxe is given by the βxxxe component, recovering ca. 75% (64%) and 73% (72%) of the βxe(βvece) values in gas and water solution, respectively. Therefore, we explored the ρxx(2)(r) densities at the CAM-B3LYP/aug-cc-pVDZ level using the numerical procedure previously illustrated by Yamada and coworkers [51]. The results evaluated at the isosurface of 0.25 a.u. are illustrated in Figure 2. As can be appreciated from the graphical representations, the ρxx(2)(r) distribution in the water solution is almost similar to that predicted in vacuum even if the amplitudes are much more spread out. This result is in some consistency with the calculated static βxxxe values, with the βxxxe(water)/βxxxe(gas) and βvece(water)/βvece(gas) ratios being computed to be 2.3 and 2.1, respectively. Note that the above ratios are somewhat greater than those previously predicted for the average αe(1.3) and γe(1.5) properties by HF/aug-cc-pVDZ computations [18].
Structure of uracil and Cartesian coordinate system. Colours: white (hydrogen), grey (carbon), red (oxygen), and cyan (nitrogen) (colour figure online).
Hyperpolarizability density distributions ρxx(2)(r) of uracil in gas (left) and water solution (right). The yellow and blue surfaces (colour figure online) refer to positive and negative ρxx(2)(r) densities, respectively, computed at the isosurface of 0.25 a.u. CAM-B3LYP/aug-cc-pVDZ results.
Figure 3 displays the frequency-dependent first-order hyperpolarizabilities computed in gaseous and aqueous phases in the 0–0.06563 ħω range for the SHG and EOPE NLO processes. It is important to notice that resonance enhancement effects for the SHG phenomenon are expected to be rather marginal, since the experimental lowest-energy absorption being placed at 5.08 eV (0.1867 a.u.) in vapour [52] and 4.77 eV (0.1753 a.u.) in water solution [53] is sufficiently far from the highest 2 ħω value of 0.1307 a.u. Not surprisingly, in the gas phase βvece(-2ω;ω;ω)>βvece(-ω;ω;0)>βvece(0;0;0). The dispersion effects evaluated at the ħω=0.06563 a.u. increase the static values by 13.7% for the EOPE and by 61.5% for the SHG process. On the other hand, in aqueous phase the dispersion effects are significantly reduced, mainly due to incomplete responses of polar solvents as modelled by the PCM treatment [31–33, 54]. As can be appreciated from Figure 3, in water solution βvece(-2ω;ω;ω)<βvece(-ω;ω;0) for the ħω values between 0.01 and 0.04 a.u., whereas βvece(-2ω;ω;ω)~βvece(-ω;ω;0) at ħω~ 0.045 a.u. It is worth noting that at the ħω value of 0.06563 a.u. the dispersion effect is negative for the EOPE phenomenon, with the βvece(-ω;ω;0)(water) value being decreased by ca. 16.0% with respect to the static datum. On the other hand, in the case of the SHG process in water solution, the dispersion effect at ħω=0.06563 a.u. is still positive as for the gas phase, even if it is noticeably inferior (+15.5%). As a consequence, although the static and dynamic electronic first-order hyperpolarizabilities in water solution are greater than the corresponding data in gas (compare the curves in Figure 3), the dispersion effects reduce the βvece(water)/βvece(gas) ratios, which are predicted to be 2.1, 1.53 and 1.48, respectively, for the static, EOPE, and SHG processes at ħω=0.06563 a.u.
Frequency-dependent electronic first-order hyperpolarizability of uracil in gas and water solution. CAM-B3LYP/aug-cc-pVDZ results. The reported data refer to the βvece(-ωσ;ω1,ω2) values obtained at ħω=0.06563 a.u.
Beside the electronic first-order hyperpolarizability, we explored the solvation effects on the pure vibrational counterpart for the EOPE phenomenon. In a recent theoretical study, Christiansen and coworkers have determined the pure vibrational first-order hyperpolarizabilities of uracil in gas using VCI computations and the Lanczos algorithm [22]. However, their reported data refer to the static and SHG process and are not directly comparable to our results. Figure 4 shows the CAM-B3LYP/aug-cc-pVDZ βvecv(-ω;ω;0) data in gaseous and aqueous phases over the 0–4000 cm−1 wavenumbers range. The largest contributions originate from the spectral region between 1500 and 2000 cm−1. Table 2 summarizes the main vibrational contributions to the βvecv(-ω;ω;0) values, also including the vibrational wavenumbers, IR intensities (IIR), and Raman activities (ARaman). The highest-energy region (wavenumbers >3000 cm−1), entirely characterized by C–H and N–H stretching vibrations, furnishes only modest contributions to the βvecv(-ω;ω;0) data, principally owing to the high-wavenumber values. On the other hand, the low-energy modes (wavenumbers <1000 cm−1) produce small and moderate βvecv(-ω;ω;0) contributions, since the vibrational transitions are rather weak in the IR [IIR∝(∂μi/∂Qa)2] and Raman [ARaman∝ (∂αi/∂Qa)2] spectra. As for the pure vibrational polarizabilities [19], in both the gas phase and water solution, the largest βvecv(-ω;ω;0) values originate from the C=O stretching vibrations with the nonnegligible contribution of the in-plane N-H bending deformation (modes ν25 and ν26). A graphical representation of the atomic displacement vectors involved in these vibrational modes is displayed in Figure 5. These transitions located at 1802 cm−1 (ν25) and 1828 cm−1 (ν26) by the present calculations in vacuum exhibit the strongest absorption peaks in the IR spectrum (Figure 6), with the IIR values of 902 and 607 km/mol, respectively. Note that the ν25 and ν26 modes are also active in the Raman spectra (ARaman= 58 and 29 Å4/amu, resp.). As a result, the βvecv(-ω;ω;0) values originated by the ν25 and ν26 transitions contribute, respectively, to ca. 89% and 50% of the total βvecv(-ω;ω;0) datum. In water solution (Figure 6), the wavenumbers of the ν25 and ν26 vibrations are downward shifted, respectively, by ca. 80 cm−1 (−4.5%) and 60 cm−1 (−3.2%) with respect to the gas phase, with the IIR values being concomitantly increased by ca. 131% and 45%, respectively. In addition, when passing from the gaseous to the aqueous phase, the ARaman(ν25) and ARaman(ν26) values enhance by ca. a factor of two and four, respectively. Therefore as for the βe values, the solvent contributions are expected to play a crucial role also for the vibrational first-order hyperpolarizabilities of uracil, increasing the βvecv(-ω;ω;0) values of the ν25 and ν26 modes by 69.4 a.u. (+120%) and 38.5 a.u. (+119%). As can be appreciated by the data reported in Table 2, for both the phases other relevant βvecv(-ω;ω;0) contributions are given by the ring stretching modes as well as by the in-plane ring bending deformations, owing to their relatively low wavenumbers and moderate IIR and ARaman values. On the whole, the introduction of solvent contributions increases the total βvecv(-ω;ω;0) by ca. 130 a.u., with the βvecv(-ω;ω;0)(water)/βvecv(-ω;ω;0)(gas) ratio being predicted to be ca. three. Finally, it is worth noting that, at ħω= 0.06563 a.u. on going from the gas phase to water solution, the βvecv(-ω;ω;0)/βvece(-ω;ω;0) ratio is almost doubled, increasing from 0.46 to 0.91.
Selected vibrational contributions to the first-order hyperpolarizabilities of uracila.
Mode no.
Wavenumbers (cm−1)
IIR (km/mol)
ARaman (Å4/amu)
Descriptionb
βvece (−ω; ω, 0) (a.u.)c
Gas
ν4
411
20
1
τring
11.5
ν5
524
22
2
δring
13.8
ν11
783
4
22
δring
10.6
ν23
1523
127
11
νring + δN-H
12.5
ν24
1711
56
30
νring + δC-H
10.4
ν25
1802
902
58
νC=O + δN-H
57.8
ν26
1828
607
29
νC=O + δN-H
32.4
Total
65.0 (139.8)
d
Water
ν5
530
45
4
δring
29.5
ν23
1527
236
42
νring + δN-H
35.6
ν24
1699
159
75
νring + δC-H
29.6
ν25
1721
2086
115
νC=O + δN-H
127.2
ν26
1769
877
116
νC=O + δN-H
70.9
Total
193.9 (214.1)
d
aCalculations were carried out at the CAM-B3LYP/aug-cc-pVDZ level on the geometry calculated at the same level. The contributions with percentage ≥15% of the total βvecv (−ω; ω, 0) value were considered.
bv: stretching, δ: in-plane bending, τ: torsion.
cThe value in parentheses refers to the CAM-B3LYP/aug-cc-pVDZ βvece (−ω; ω, 0) value at ħω = 0.06563 a.u.
dThe βvecv(-ω;ω,0)/βvece(-ω;ω,0) ratios are 0.46 and 0.91 in gas and water solution, respectively.
Contribution of each normal mode to the vibrational first-order hyperpolarizability of uracil in gas and water solution. CAM-B3LYP/aug-cc-pVDZ results.
Atom vector displacements of the ν C=O + δN-H modes ν25 and ν26.
IR and Raman spectra of uracil in gas and water solution in the 1000–1900 cm−1 wavenumbers range. Lorentz line shapes with a full width at half maximum of 10 cm−1 were used. CAM-B3LYP/aug-cc-pVDZ results.
4. Conclusions
We have examined the effects of solvation on the static and frequency-dependent electronic and vibrational first-order hyperpolarizabilities of uracil. The properties were modeled in vacuum as well as in water solution using the PCM approach. The calculations were carried out using the long-range corrected CAM-B3LYP functional with the Dunning’s correlation-consistent aug-cc-pVDZ basis set. The introduction of solvent contributions significantly increases both the electronic and vibrational first-order hyperpolarizabilities. However, the dispersion effects on the electronic hyperpolarizabilities for the EOPE and SHG NLO phenomena are noticeably reduced when passing from the gas phase to the water solution. The magnitudes of the vibrational properties are comparable to the electronic counterparts, with the βvecv/βvece ratio increasing with the solvation and βvecv(water)~βvece(water) for the EOPE process at the characteristic wavelength of 694 nm. The most relevant contributing modes to the βv values principally involve the very intense infrared C=O stretching + N-H in-plane bending deformation vibrations.
Conflict of Interests
The author declares that there is no conflict of interests regarding the publication of this paper.
ChemlaD. S.ZyssJ.19871-2San Diego, Calif, USAAcademic PressPrasadP. N.WilliamsD. J.1991New York, NY, USAJohn Wiley & SonsNalwaH. S.MiyataS.1997Boca Raton, Fla, USACRC PressWangL.YoshidaJ.OgataN.SasakiS.KajiyamaT.Self-assembled supramolecular films derived from marine deoxyribonucleic acid (DNA)-cationic surfactant complexes: large-scale preparation and optical and thermal properties2001134127312812-s2.0-003483852310.1021/cm000869gYanH.ZhangX.ShenZ.SeemanN. C.A robust DNA mechanical device controlled by hybridization topology2002415686762652-s2.0-003701198410.1038/415062aPorathD.CunibertiG.FeliceR. D.SchusterG. B.Charge transport in DNA based devices2004237Heidelberg, GermanySpringer183227EndresR. G.CoxD. L.SinghR. R. P.Colloquium: the quest for high-conductance DNA20047611952142-s2.0-184276302710.1103/RevModPhys.76.195GroteJ. G.HagenJ. A.ZettsJ. S.NelsonR. L.DiggsD. E.StoneM. O.YaneyP. P.HeckmanE.ZhangC.SteierW. H.JenA. K.-Y.DaltonL. R.OgataN.CurleyM. J.ClarsonS. J.HopkinsF. K.Investigation of polymers and marine-derived DNA in optoelectronics200410825858485912-s2.0-314266185310.1021/jp038056dCohenH.NoguesC.NaamanR.PorathD.Direct measurement of electrical transport through single DNA molecules of complex sequence20051023311589115932-s2.0-2384455026510.1073/pnas.0505272102GuerraC. F.Van Der WijstT.BickelhauptF. M.Nanoswitches based on DNA base pairs: why adenin-thymine is less suitable than guanine-cytosine200679197119792-s2.0-3374894858510.1002/cphc.200600216StecklA. J.DNA-a new material for photonics?200711352-s2.0-3424748837910.1038/nphoton.2006.56CohenH.SapirT.BorovokN.MolotskyT.Di FeliceR.KotlyarA. B.PorathD.Polarizability of G4-DNA observed by electrostatic force microscopy measurements2007749819862-s2.0-3424816863310.1021/nl070013bCamposC. T.JorgeF. E.Basis set convergence of electric properties in HF and DFT calculations of nucleic acid bases200910922852932-s2.0-5814921552710.1002/qua.21792LabelloN. P.FerreiraA. M.TcurtzH. A.An augmented effective core potential basis set for the calculation of molecular polarizabilities20052614146414712-s2.0-2774451883310.1002/jcc.20282JasienP. G.FitzgeraldG.Molecular dipole moments and polarizabilities from local density functional calculations: application to DNA base pairs1990934255425602-s2.0-0000715825JohnsonR. C.PowerT. D.HoltJ. S.ImmarapornB.MonatJ. E.SissokoA. A.YanikM. M.ZagorodnyA. V.CybulskiS. M.Electron-correlated calculations of electric properties of nucleic acid bases19961004818875188812-s2.0-33751061034BaschH.GarmerD. R.JasienP. G.KraussM.StevensW. J.Electrical properties of nucleic acid bases198916365145222-s2.0-0010894596AlparoneA.Linear and nonlinear optical properties of nucleic acid bases2013410909810.1016/j.chemphys.2012.11.005MillefioriS.AlparoneA.Tautomerism and polarizability in uracil: coupled cluster and density-functional theory study20043031-227362-s2.0-304274518610.1016/j.chemphys.2004.05.002PlutaT.KolaskiM.Medved’M.BudzákS.Dipole moment and polarizability of the low-lying excited states of uracil20125462429HaleyT. P.GraybillE. R.CybulskiS. M.Ab initio calculations of dispersion coefficients for nucleic acid base pairs2006124202-s2.0-3454792634910.1063/1.2197832204301HansenM. B.SeidlerP.GyrffyW.ChristiansenO.A Lanczos-chain driven approach for calculating damped vibrational configuration interaction response functions2010133112-s2.0-7795696204410.1063/1.3467889114102BishopD. M.Molecular vibration and nonlinear optics19981041402-s2.0-54649084336TomasiJ.MennucciB.CammiR.Quantum mechanical continuum solvation models20051058299930942-s2.0-2404447340110.1021/cr9904009ScalmaniG.FrischM. J.Continuous surface charge polarizable continuum models of solvation. I. General formalism2010132112-s2.0-7794968389210.1063/1.3359469114110BartkowiakW.MisiaszekT.Solvent effect on static vibrational and electronic contribution of first-order hyperpolarizability of π-conjugated push-pull molecules: quantum-chemical calculations200026133533572-s2.0-003466975810.1016/S0301-0104(00)00262-7CammiR.MennucciB.TomasiJ.Solvent effects on linear and nonlinear optical properties of donor- acceptor polyenes: investigation of electronic and vibrational components in terms of structure and charge distribution changes199812034883488472-s2.0-003247544810.1021/ja980823cAlparoneA.Theoretical study of the electronic (hyper)polarizabilities of amino acids in gaseous and aqueous phases20119761881902-s2.0-84555169190AlparoneA.Corrigendum to “Theoretical study of the electronic (hyper)polarizabilities of amino acids in gaseous and aqueous phases” [Comput. Theor. Chem. 976 (2011) 188–190]20129801442-s2.0-8455516919010.1016/j.comptc.2011.11.033AlparoneA.MillefioriS.Gas and solution phase electronic and vibrational (hyper)polarizabilities in the series formaldehyde, formamide and urea: CCSD(T) and DFT theoretical study20054164–62822882-s2.0-2844445423510.1016/j.cplett.2005.09.088CammiR.CossiM.MennucciB.TomasiJ.Analytical Hartree-Fock calculation of the dynamical polarizabilities α, β, and γ of molecules in solution19961052310556105642-s2.0-0000987405CammiR.CossiM.MennucciB.TomasiJ.Solvent effects on static and dynamic polarizability and hyperpolarizabilities of acetonitrile1997436-4375675752-s2.0-003157433310.1016/S0022-2860(97)00209-3MennucciB.CammiR.CossiM.TomasiJ.Solvent and vibrational effects on molecular electric properties. Static and dynamic polarizability and hyperpolarizabilities of urea in water19984261–31911982-s2.0-0010018614FrischM. J.TrucksG. W.SchlegelH. B.ScuseriaG. E.RobbM. A.CheesemanJ. R.ScalmaniG.BaroneV.MennucciB.PeterssonG. A.NakatsujiH.CaricatoM.LiX.HratchianH. P.IzmaylovA. F.BloinoJ.ZhengG.SonnenbergJ. L.HadaM.EharaM.ToyotaK.FukudaR.HasegawaJ.IshidaM.NakajimaT.HondaY.KitaoO.NakaiH.VrevenT.MontgomeryJ. J. A.PeraltaJ. E.OgliaroF.BearparkM.HeydJ. J.BrothersE.KudinK. N.StaroverovV. N.KobayashiR.NormandJ.RaghavachariK.RendellA.BurantJ. C.IyengarS. S.TomasiJ.CossiM.RegaN.MillamJ. M.KleneM.KnoxJ. E.CrossJ. B.BakkenV.AdamoC.JaramilloJ.GompertsR.StratmannR. E.YazyevO.AustinA. J.CammiR.PomelliC.OchterskiJ. W.MartinR. L.MorokumaK.ZakrzewskiV. G.VothG. A.SalvadorP.DannenbergJ. J.DapprichS.DanielsA. D.FarkasO.ForesmanJ. B.OrtizJ. V.CioslowskiJ.FoxD. J.2009Wallingford, Conn, USAGaussianYanaiT.TewD. P.HandyN. C.A new hybrid exchange-correlation functional using the Coulomb-attenuating method (CAM-B3LYP)20043931–351572-s2.0-314277129710.1016/j.cplett.2004.06.011WoonD. E.DunningT. H.Gaussian basis sets for use in correlated molecular calculations. IV. Calculation of static electrical response properties19941004297529882-s2.0-26444474547AlparoneA.Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers2013193095310210.1007/s00894-013-1838-4AlparoneA.The effect of secondary structures on the NLO properties of single chain oligopeptides: a comparison between β-strand and α-helix polyglycines2013151295812962MedvedM.BudzákS.PlutaT.Static NLO responses of fluorinated polyacetylene chains evaluated with long-range corrected density functionals20115151–378842-s2.0-8005390305610.1016/j.cplett.2011.08.083JacqueminD.PerpèteE. A.MedvedM.ScalmaniG.FrischM. J.KobayashiR.AdamoC.First hyperpolarizability of polymethineimine with long-range corrected functionals2007126192-s2.0-3424905903410.1063/1.2741246191108UrdanetaJ.BermúdezY.ArrietaF.RosalesM.CubillánN.HernándezJ.CastellanoO.SoscúnH.Theoretical study in gas phase of linear and nonlinear optical properties of the ortho-, meta- and para-nitrophenol isomers2010126127372-s2.0-7795149930410.1007/s00214-009-0668-zAlparoneA.Comparative study of CCSD(T) and DFT methods: electronic (hyper)polarizabilities of glycine20115141–321252-s2.0-8005292279610.1016/j.cplett.2011.08.010ZhangC.-C.XuH.-L.HuY.-Y.SunS.-L.SuZ.-M.Quantum chemical research on structures, linear and nonlinear optical properties of the Li@n-acenes salt (n=1, 2, 3, and 4)201111510203520402-s2.0-7995258186010.1021/jp110412nKrawczykP.DFT study of linear and nonlinear optical properties of donor-acceptor substituted stilbenes, azobenzenes and benzilideneanilines20101646596682-s2.0-7795460434010.1007/s00894-009-0623-xMillefioriS.AlparoneA.Theoretical determination of the vibrational and electronic (hyper)polarizabilities of C4H4X (X=O, S, Se, Te) heterocycles2000211249525012-s2.0-003421416710.1039/b001020kBishopD. M.HasanM.KirtmanB.A simple method for determining approximate static and dynamic vibrational hyperpolarizabilities199510310415741592-s2.0-0000032412KeshariV.WijekoonW. M. K. P.PrasadP. N.KarnaS. P.Hyperpolarizabilities of organic molecules: Ab initio time-dependent coupled perturbed Hartree-Fock-Roothaan studies of basic heterocyclic structures19959922904590502-s2.0-33751155296BrownR. D.GodfreyP. D.McNaughtonD.PierlotA. P.Microwave spectrum of uracil19881107232923302-s2.0-0023923033KulakowskaI.GellerM.LesyngB.Barrier to rotation and conformation of the NR2 group in cytosine and its derivatives. Part II. Experimental and theoretical dipole moments of methylated cytosines197540744204292-s2.0-0016786835NakanoM.ShigemotoI.YamadaS.YamaguchiK.Size-consistent approach and density analysis of hyperpolarizability: second hyperpolarizabilities of polymeric systems with and without defects199510310417541912-s2.0-33845765038YamadaS.NakanoM.ShigemotoI.KiribayashiS.YamaguchiK.Intense electron correlation dependence of the first hyperpolarizabilities β of a nitroxide radical and formaldehyde19972675-64454512-s2.0-0031588868ClarkL. B.PeschelG. G.TinocoI.Jr.Vapor spectra and heats of vaporization of some purine and pyrimidine bases19656910361536182-s2.0-0001456998VoetD.GratzerW. B.CoxR. A.DotyP.Absorption spectra of nucleotides, polynucleotides, and nucleic acids in the far ultraviolet19631319320810.1002/bip.360010302CammiR.FredianiL.MennucciB.RuudK.Multiconfigurational self-consistent field linear response for the polarizable continuum model: theory and application to ground and excited-state polarizabilities of para-nitroaniline in solution200311912581858272-s2.0-014202662210.1063/1.1603728