Conformations and Vibrational Frequencies of a Precursor of Benzovesamicol Analogues Studied by Density Functional Theories

Bulletin of the Korean Chemical Society.
2014.
Aug,
35(8):
2311-2316

- Received : October 04, 2013
- Accepted : April 17, 2014
- Published : August 20, 2014

Download

PDF

e-PUB

PubReader

PPT

Export by style

Article

Metrics

Cited by

TagCloud

Conformations and vibrational frequencies of the racemic (2
RS
,3
RS
)-5-amino-3-(4-phenylpiperazin-1-yl)-1,2,3,4-tetrahydronaphthalen-2-ol-(I) [(2
RS
,3
RS
)-(I)], a precursor of benzovesamicol analogues, have been carried out using various DFT methods (M06-2X, B3LYP, B3PW91, PBEPBE, LSDA, and B3P86) with basis sets of 6-31G(d), 6-31+G(d,p), 6-311+G(d,p), 6-311++G(d,p), cc-pVTZ, and TZVP. The LSDA/6-31G(d) level of theory shows the best performance in reproducing the X-ray powder structure. However, the PBEPBE/cc-pVTZ level of theory is the best method to predict the vibrational frequencies of (2
RS
,3
RS
)-(I). The potential energy surfaces of racemic pairs (2
RS
,3
RS
)-(I) and -(II) are obtained at the LSDA/6-31G(d) level of theory in the gas phase and in water. The results indicate that (2
RS
,3
RS
)-(I) are more stable by ~0.75 kcal/mol in energy than (2
RS
,3
RS
)-(II) in water, whereas conformer AIIg and BIIg are more stable by ~0.04 kcal/mol than AIg in gas phase. In particular, the hydrogen bond distances between the N of piperazine and the OH of tetrahydronaphthalen become longer in gas, compared with those in the water phase. Vibrational frequencies calculated at the PBEPBE/cc-pVTZ level of theory in the gas phase are larger than those in water, whereas their intensities in the gas phase are weaker than those in water.
RS
,3
RS
)-5-amino-3-(4-phenylpiperazin-1-yl)-1,2,3,4-tetrahydronaphthalen-2-ol-(I) [(2
RS
,3
RS
)-(I)] is a precursor for the preparation of benzovesamicol analogues, which are stereoselective inhibitors of acetylcholine uptake in presynaptic cholinergic vesicles. Benzovesamicol analogues and selected atomic numbers for (2
RS
, 3
RS
)-(I) are shown in
Figure 1
. Radiolabeled benzovesamicol analogues have been extensively used as imaging probes in single-photon emission computedand positron emission tomography in
in vitro
and
in vivo
studies of Alzheimer’s disease.
^{1-8}
Assad
et al
. have recently synthesized (2
RS
,3
RS
)-(I) and analyzed its structure by FT-IR and NMR measurements and X-ray powder diffraction.
^{9}
Precusors of benzovesamicol analogues. (a) (2RS ,3RS )-(I) and (2RS ,3RS )-(II) are racemic compounds. Ψ _{1} and Ψ _{2} are dihedral angles. (b) Selected atomic numbers for (2RS ,3RS )-(I).
DFT has been extensively used to study various aspects of molecules with main group atomic species.
^{10-16}
Approximations in DFT permit the accurate calculation of certain physical quantities. The most widely used approximation in physics is the local-density approximation (LDA). Local spin-density approximation (LSDA) is the straightforward generation of LDA including electron spin. The exchange energy in hybrid methods (
e.g.
, M06-2X, B3LYP, B3PW91, PBEPBE, and B3P86) is combined with the exact energy from Hartree–Fock theory. Adjustable parameters in hybrid functionals are generally fitted to a training set of molecules. Although results obtained from hybrid functionals are usually sufficiently accurate for most applications, a systematic method to improve these functionals has yet to be established and compared with several traditional wave-function-based methods, such as configuration interaction or coupled cluster theory. Thus, calculation errors for the conformational structure and properties of (2
RS
,3
RS
)-(I) should be estimated and compared with other methods or experiments.
^{9}
Applications of different DFT methods and various basis sets to provide accurate conformations and vibrational frequencies of (2
RS
,3
RS
)-(I) have not been analyzed yet. This study has the following two main objectives: (1) to investigate the performance of different DFT methods and effect of different basis sets on predicting the conformational structures and vibrational spectra of (2
RS
,3
RS
)-(I) in gas phase and water solution, and (2) to study the conformational structures and vibrational spectra of (2
RS
,3
RS
)-(II). This study provides theoretical data and insights into the conformational structure preferences of racemic (2
RS
,3
RS
)-(I).
^{17}
All structures were initially optimized at the HF/6-31G(d) starting the X-ray power structure. The structure of (2
RS
,3
RS
)-(I) were reoptimized and its vibrational frequencies were calculated in the gas phase using various DFT methods, including M06-2X,
^{18}
B3LYP,
^{19}
LSDA,
^{20}
B3PW91,
^{21}
PBEPBE,
^{22}
and B3P86,
^{23}
and different basis sets, including 6-31G(d), 6-31+G(d,p), 6-311+G(d,p), 6-311++G(d,p), cc-pVTZ, and cc-pVQZ. In addition, the potential surfaces of (2
RS
,3
RS
)-(I) and -(II) were obtained at the LSDA/6-31G(d) level of theory in the gas phase and in water. The solvation free energies were calculated using the implicit Solvation Model based on Density (SMD) method.
^{24}
The scale factor used for vibrational frequencies was 0.944.
^{25}
Comparison of Structures Optimized with Different Methods and Various Basis Sets for (2RS ,3RS )-(I).
The molecular structure of (2
RS
,3
RS
)-(I) with numerical labels is shown in
Figure 1
. Table S1 (Supporting Information) provides the comparison of RMSDs between calculated coordinates in the gas phase among different DFT methods at the 6-31G(d) basis set and the X-ray data for (2
RS
,3
RS
)-(I). RMSD can obtain a more accurate conformational structure among all the methods tested.
^{26}
RMSD is the measure of the average distance between the atoms (usually the backbone atoms) of a superimposed 3D structure.
The calculated RMSDs are 2.7793 Å for M06-2X, 2.7923 Å for B3LYP, 2.7777 Å for LSDA, 2.7854 Å for B3PW91, 2.7958 Å for PBEPBE, and 2.7849 Å for B3P86. Thus, of all the methods tested, the LSDA method has the least RMSD between the calculated coordinates of the gas phase and experimental set-up (powder).
^{9}
As shown in Table S1, differences exist between the calculated and X-ray parameters because experimental values were obtained from molecules in solid powder state, whereas theoretical values were based on an isolated molecule in gas phase.
The comparison of RMSDs between calculated coordinates in the gas phase and X-ray data using various basis sets at LSDA calculation for (2
RS
,3
RS
)-(I) is shown in Table S2. The RMSDs are as follows: 2.7777 Å for 6-31G(d), 2.7819 Å for 6-31+G(d, p), 2.7850 Å for 6-311+G(d,p), 2.7853 Å for 6-311++G(d,p), 2.7819 Å for cc-PVTZ, and 2.7822 Å for TZVP. Thus, of all basis sets tested, 6-31G(d) has the least RMSD. However, the RMSDs between the calculated coordinates in the gas phase and X-ray data are significantly high. We compare the selected structural parameters of the X-ray and calculated parameters to identify the reason for the significantly high deviations.
Selected structural parameters calculated using the LSDA/6-31G(d) method for (2
RS
,3
RS
)-(I) are shown in
Table 1
. In the comparison of the X-ray-generated structure and the calculated structure, the average error of bond length is approximately 0.61%, the bond angle is approximately 0.37%, and the average absolute dihedral angle is approximately 2.1%. Compared with the dihedral angle (
Ψ
_{1}
= −38.1°) of X-ray data, the calculated angle (
Ψ
_{1}
= 4°) of gas phase is significantly different in particular. Thus, conformational energy barrier stability is dependent on
Ψ
_{1}
, as shown in
Figure 2
. As expected, the energy surface of conformational structures is stabilized with the formation of an intramolecular hydrogen bond at
Ψ
_{2}
= –156° to –162°, as shown in
Table 2
. The dihedral angle (
Ψ
_{1}
) is defined as almost free rotation without steric hindrance. These results indicate that LSDA is consistent with the X-ray data in terms of predicting the conformational structure of (2
RS
,3
RS
)-(I). In addition, the ability of DFT methods to provide reliable molecular structures for the molecule with the main group atomic species is already well-established and used in the scientific field. However, based on RMSD, the LSDA method is selected for structural calculations in this study.
^{a} Ref 9 . ^{b} Intramolecular hydrogen bond distance (O1H---N2). ^{c} See Figure 1 .
Plots of the 2D contour map in the variation in relative energy (ΔE in kcal/mol) with Ψ _{1} and Ψ _{2} for (2RS ,3RS )-(I). The structures of labeled AIg-DIg conformers are labeled to correspond to calculated structures of AIg-DIg.
^{a} Relative energies (ΔE ) in Kcal/mol specify differences in total energy (see Table S3 of Supporting Information). Minus sign represents more stable. ^{b} Intramolecular hydrogen bond between a hydroxyl H atom and an N_{2} atom of the piperazine ring (O1-H1O…N2) (see Figure 2 )
Potential Surfaces of (2RS ,3RS )-(I).
Plots of the two-dimensional (2D) contour map in the variation in relative energy (Δ
E
in kcal/mol) with
y
_{1}
and
y
_{2}
for (2
RS
,3
RS
)-(I) are shown in
Figure 2
. The conformational structures of labeled AIg–DIg conformer are labeled to correspond to the calculated structures of AIg-DIg (
Figure 3
).
Preferred conformations AIg, BIg, CIg, and DIg of (2RS ,3RS )-(I) at the LSDA/6-31G(d) level of theory in the gas phase. Intramolecular hydrogen bonds are represented by a dashed line. The distances of hydrogen bond for AIg, BIg, and DIg are 1.92 Å and CIg is 1.93 Å.
Table 2
lists the molecular dihedral angles and relative energies of benzovesamicol analogues, (2
RS
,3
RS
)-(I) and (2
RS
,3
RS
)-(II), optimized at LSDA/6-31G(d) levels in gas phase and water solution. The conformational structures and relative energies of local minima are almost similar. Although
Ψ
_{1}
of conformers significantly differ in dihedral angles,
Ψ
_{2}
angles resulting from intramolecular hydrogen bonds have similar angles. The conformational structure stability of (2
RS
,3
RS
)-(I) by the relative energies (Δ
E
) in gas phase is calculated. Δ
E
specifies differences in total energy (Table S3), and is expressed as kcal/mol with respect to the stable conformational structure. The lowest relative energy conformation is AIg with intramolecular hydrogen bond between a hydroxyl H atom and an N
_{2}
atom of the piperazine ring (O1-H1O…N2) in
Figure 1
. The hydrogen bond functions in determining the lowest energy conformation for benzovesamicol analogue structures. The distances of the hydrogen bond for (2
RS
,3
RS
)-(I) are calculated as 1.92 and 1.93 Å. All Δ
E
values computed with the LSDA/6-31G(d) method, using either one of the optimized conformers, differ by ~0.39 kcal/mol in (2
RS
,3
RS
)-(I) in gas phase.
The energies of conformers at (2
RS
,3
RS
)-(II) are lower than those at (2
RS
,3
RS
)-(I). Thus, the conformers of (2
RS
,3
RS
)-(II) are more stable than those of (2
RS
,3
RS
)-(I). Figure S1 shows the preferred structures AIIg, BIIg, and CIIg of (2
RS
,3
RS
)-(II) at the LSDA/6-31G(d) level in gas phase. The distances of hydrogen bonds are 1.91 and 1.92 Å. All Δ
E
values computed with the LSDA/6-31G(d) method are shown in
Table 2
. The conformer AIIg and BIIg are more stable than AIg in gas phase.
The conformational structure stabilities of (2
RS
,3
RS
)-(I) and (2
RS
,3
RS
)-(II) by relative energies in water solution are calculated. The lowest relative energy conformer is AIs. Either one of the optimized conformers differ by ~0.75 kcal/mol in water solution. Compared with gas phase, the hydrogen bond distances decrease in water solution. Figure S2 shows the preferred conformational structures AIs, BIs, and CIs of (2
RS
,3
RS
)-(I) at the LSDA/6-31G(d) level in water solution. The distances of hydrogen bonds for AIs, BIs, and CIs are 1.88 Å but 1.85 Å for BIs. (2
RS
,3
RS
)-(II) exhibits the structures of AIIs, BIIs, and CIIs in water. The hydrogen bond distances of AIIs and BIIS are 1.90 Å, but CIIs have hydrogen bond distances of 1.85 Å.
Comparison of Vibrational Frequencies Calculated with Different Methods and Basis Sets for (2RS ,3RS )-(I).
Table 3
shows the calculated results of vibrational frequencies with different DFT methods at 6-31G(d) basis set for (2
RS
,3
RS
)-(I). Table S4 shows a comparison of mean absolute deviation (MAD, cm
^{–1}
) between calculated vibrational frequencies and experimental values using various DFT methods at 6-31G(d) basis set for (2
RS
,3
RS
)-(I). The MAD values are as follows: 193.1 cm
^{–1}
for M06-2X, 133.8 cm
^{–1}
for B3LYP, 80.4 cm
^{–1}
for LSDA, 133.1 cm
^{–1}
for B3PW91, 75.6 cm
^{–1}
for PBEPBE, and 138.6 cm
^{–1}
for B3P86.
^{a} Ref. 9
The simulation spectra of (2
RS
,3
RS
)-(I) are shown in Figures S3 and S4. PBEPBE method has the least MAD between calculated and experimental values, and the deviation error is 2.89%. Thus, PBEPBE is selected for the vibrational frequency calculations. The DFT calculation is desirable for resolving disputes in vibrational assignments, and provides valuable insight for understanding the observed spectral features.
Table 4
shows the calculated results of vibrational frequencies with PBEPBE at different basis set for (2
RS
,3
RS
)-(I). Table S5 shows the comparison of MAD between calculated vibrational frequencies and experimental values using various basis sets for (2
RS
,3
RS
)-(I) by PBEPBE method. The MAD values are as follows: 51.7 cm
^{−1}
for 6-31+G(d,p), 50.7 cm
^{−1}
for 6-311+G(d,p), 50.6 cm
^{−1}
for 6-311++G(d,p), 45.9 cm
^{−1}
for cc-PVTZ, and 122.6 cm
^{−1}
for TZVP. MAD of cc-PVTZ basis set is the lowest at 45.9 cm
^{–1}
, and the deviation error is 1.59%. Thus, PBEPBE/cc-PVTZ is selected for vibrational frequency calculation. The simulation spectra of (2
RS
,3
RS
)-(I) are shown in Figure S5. These results indicate that PBEPBE/cc-PVTZ is consistent with available experimental vibrational frequencies.
^{9}
^{a} Ref. 9
Analysis of Vibrational Frequencies for (2RS ,3RS )-(I) and (2RS ,3RS )-(II).
Table 5
shows the vibrational frequencies of benzovesamicol analogues, (2
RS
,3
RS
)-(I) and (2
RS
,3
RS
)-(II), calculated at PBEPBE/cc-PVTZ level in gas phase and water solution. The frequency of the OH stretching region is 3488 cm
^{–1}
for (2
RS
,3
RS
)-(I) and 3479 cm
^{–1}
for (2
RS
,3
RS
)-(II) in gas phase. In (2
RS
,3
RS
)-(I), the stretching regions of (
CH
=CH, Ar), (
CH_{2}
-NH
_{2}
), (C=C), and (C-N, Piperazine) are 3084, 2853, 1602, and 1140 cm
^{–1}
, respectively. The symmetry stretch values of (NH
_{2}
) and (CH
_{2}
, aliphatic) are 3472 and 2908 cm
^{–1}
, respectively. In (2RS,3RS)-(II), the stretch values of (
CH
=CH, Ar), (
CH_{2}
-NH
_{2}
), (C=C), and (C-N, Piperazine) are 3115, 2869, 1603, and 1140 cm
^{–1}
, respectively. The symmetry stretch values of (NH
_{2}
) and (CH
_{2}
, aliphatic) are 3477 and 2925 cm
^{–1}
, respectively. The frequencies of (2
RS
,3
RS
)-(II) slightly increase than those of (2
RS
,3
RS
)-(I), except for the OH frequency and (C-N, Piperazine). For water solution, the stretch values of (OH), (
CH
=CH, Ar), (
CH_{2}
-NH
_{2}
), (C=C), and (C-N, Piperazine) are 3308, 3036, 2815, 1542, and 1115 cm
^{–1}
in (2
RS
,3
RS
)-(I), respectively. The symmetry stretch values of (NH
_{2}
) and (CH
_{2}
, aliphatic) are 3377 and 2861 cm
^{–1}
, respectively. In (2
RS
,3
RS
)-(II), the stretch values of (OH), (
CH
=CH, Ar), (
CH_{2}
-NH
_{2}
), (C=C), and (C-N, Piperazine) are 3433, 3090, 2890, 1565, and 1125 cm
^{–1}
, respectively. The symmetry stretch values of (NH
_{2}
) and (CH
_{2}
, aliphatic) are 3459 and 2931 cm
^{–1}
, respectively. The frequencies of gas phase are higher than those of water solution because of the polar solvent effect, whereas the intensities of gas phase are weaker than those of water solution in (2
RS
,3
RS
)-(I). The simulation spectra of (2
RS
,3
RS
)-(I) and (2
RS
,3
RS
)-(II) in gas phase and water solution are shown in Figure S6. In (2
RS
,3
RS
)-(II), the frequencies of (CH
_{2}
, aliphatic) and (
CH_{2}
-NH
_{2}
) in gas phase are lower than those of water solution.
^{a} Ref. 9

Benzovesamicol analogues
;
DFT methods
;
Conformations
;
Vibrational frequencies
;
Diagnosis of Alzheimer’s disease

Introduction

Racemic (2
PPT Slide

Lager Image

Computational Methods

All calculations were performed with the Gaussian 09 package.
Results and Discussion

Selected Structure Parameters (Å,o) calculated with LSDA/6-31G(d) method for (2RS,3RS)-(I) and (2RS,3RS)-(II)

PPT Slide

Lager Image

PPT Slide

Lager Image

Molecular dihedral angles, relative energies (Erel, Kcal/mol)a, and hydrogen bond distancebof (2RS,3RS)-(I) and (2RS,3RS)-(II) optimized at the LSDA/6-31G(d) levels in gas phase and water solution

PPT Slide

Lager Image

PPT Slide

Lager Image

Comparison of vibrational frequencies (in nm) calculated with various DFT methods at 6-31G(d) basis set for (2RS,3RS)-(I) in gas phase

PPT Slide

Lager Image

Comparison of vibrational frequencies (cm–1) calculated with PBEPBE at different basis sets for (2RS,3RS)-(I) in gas phase

PPT Slide

Lager Image

Vibrational frequencies (cm–1) of (2RS,3RS)-(I) and (2RS,3RS)-(II) calculated at PBEPBE/cc-PVTZ level in gas phase and water solution

PPT Slide

Lager Image

Conclusion

The main findings of this study are summarized as follows:
- a) The LSDA/6-31G(d) level is consistent with the results from X-ray crystallography in predicting the conformational structure of racemic pairs (2RS,3RS)-(I) and (2RS,3RS)-(II), whereas the PBEPBE/cc-PVTZ level is consistent with the vibrational frequencies.
- b) Compared with calculated structures, bond lengths and bond angles are consistent in X-ray-generated structures, except dihedral angles.
- c) Dihedral angle (ψ1) is highly different between calculated gas phase (4°) and X-ray data (–38.1°). Thus, the dihedral angle (ψ1) is defined as almost free rotation without steric hindrance.
- d) (2RS,3RS)-(I) are more stable than (2RS,3RS)-(II) in water, whereas conformer AIIg and BIIg are more stable than conformer AIg in gas phase.
- e) The hydrogen bond distances become longer in gas, compared with those in water.
- f) The frequencies of gas phase are higher than those of water solution because of the polar solvent effect, whereas the intensities of gas phase are weaker than those of water solution for (2RS,3RS)-(I).

Acknowledgements

Supporting Information Available.Supporting information includes RMSD lists between calculated coordinates in gas phase and experiment using different DFT methods and various basis sets for (2RS,3RS)-(I). Molecular dihedral angles and total energies of (2RS,3RS) and (2RS,3RS)-(II) are shown in Table S3. MAD values between calculated vibrational frequencies and experimental frequencies are given in Tables S4 and S5. Preferred conformational structures and calculated IR spectra are shown in Figures S1–S5. This material is available free of charge.

Alfonso A.
,
Grundahl K.
,
Duerr J. S.
,
Han H. P.
,
Rand J. B.
1993
Science
261
617 -
619
** DOI : 10.1126/science.8342028**

Auld D. S.
,
Kornecook T. J.
,
Bastianetto S.
,
Quirion R.
2002
Prog. Neurobiol.
68
209 -
245
** DOI : 10.1016/S0301-0082(02)00079-5**

Efange S. M. N.
,
Garland E.
,
Staley J. K.
,
Khare A. B.
,
Mash D. C.
1997
Neurobiol. Aging
18
407 -
413
** DOI : 10.1016/S0197-4580(97)00038-9**

Mulholland G. K.
,
Jung Y. W.
1992
J. Labelled Comp. Radiopharm.
31
253 -
259
** DOI : 10.1002/jlcr.2580310402**

Mulholland G. K.
,
Jung Y. W.
,
Wieland D. M.
,
Kilbourn M. R.
,
Kuhl D. E.
1993
J. Labelled Comp. Radiopharm.
33
583 -
591
** DOI : 10.1002/jlcr.2580330704**

Nicolas G.
,
Patrick E.
,
Roger R. F.
,
David J. H.
,
Sylvie C.
,
Lucette G.
,
Peter R.
,
Stefan E.
,
Sylvie B.
,
Michael J. F.
,
Denis G.
,
Michael K.
2007
Synapse
61
962 -
979
** DOI : 10.1002/syn.20450**

Rogers G. A.
,
Parsons S. M.
,
Anderson D. C.
,
Nilsson L. M.
,
Bahr B. A.
,
Kornreich W. D.
,
Kaufman B.
1989
J. Med. Chem.
32
1217 -
1230
** DOI : 10.1021/jm00126a013**

Zea-Ponnce Y.
,
Mavel S.
,
Assaad T.
,
Kruse S. E.
,
Parson S. M.
,
Emond P.
,
Chalon S.
,
Kruse S.
,
Giboureau N.
,
Kassiou M.
,
Guilloteau D.
2005
Bioorg. Med. Chem.
13
745 -
753
** DOI : 10.1016/j.bmc.2004.10.043**

Assad T.
,
Rukiah M.
2011
Acta Cryst.
C67
o378 -
o381

Park J. K.
,
Choe S. J.
2013
J. Porphyrins Phthalocyanines
17
376 -
383
** DOI : 10.1142/S1088424613500582**

Choe S. J.
2012
Bull. Korean Chem. Soc.
33
2861 -
2866
** DOI : 10.5012/bkcs.2012.33.9.2861**

Choe S. J.
,
Shim Y. K.
,
Galindev O. J.
2012
J. Porphyrins Phthalocyanines
16
218 -
226
** DOI : 10.1142/S1088424612004537**

Camacho C.
,
Cheng C.
,
Witek H. A.
,
Lee Y. P.
2010
J. Phys. Chem. A
114
11008 -
11016
** DOI : 10.1021/jp105664k**

Gao H.
,
Wei X.
,
Liu X.
,
Yan T.
2010
J. Phys. Chem. B
114
4056 -
4062
** DOI : 10.1021/jp912005a**

Zakharov A. V.
,
Girichev G. V.
2008
J. Mol. Struct. (Theochem)
851
183 -
196
** DOI : 10.1016/j.theochem.2007.11.008**

Gao H.
,
Yan T.
,
Zhang C.
,
He H.
2008
Mol. Struct. (Theochem)
857
38 -
43
** DOI : 10.1016/j.theochem.2008.02.004**

Frisch M. J.
,
Trucks G. W.
,
Schlegel H. B.
,
Scuseria G. E.
,
Robb M. A.
,
Cheeseman J. R.
,
Scalmani G.
,
Barone V.
,
Mennucci B.
,
Peterson G. A.
,
Nakatsuji H.
,
Caricato M.
,
Hratchian H. P.
,
Izmaylov A. F.
,
Bloino J.
,
Zheng G.
,
Sonnenberg J. L.
,
Hada M.
,
Ehara M.
,
Toyota K.
,
Fukuda R.
,
Hasegawa J.
,
Ishida M.
,
Nakajima T.
,
Hona Y.
,
Kitao O.
,
Nakai H.
,
Vreven T.
,
Montgomery, J. A. Jr.
,
Peralta J. E.
,
Ogliaro F.
,
Bearpark M.
,
Heyd J. J.
,
Brothers E.
,
Kudin K. N.
,
Staroverov V. N.
,
Kobayashi R.
,
Normand J.
,
Raghavachari K.
,
Rendell A.
,
Burant J. C.
,
Lyengar S. S.
,
Tomasi J.
,
Cossi M.
,
Rega N.
,
Millam J. M.
,
Klene M.
,
Knox J. E.
,
Cross J. B.
,
Bakken V.
,
Adamo C.
,
Jaramillo J.
,
Gomperts R.
,
Stratmann R. E.
,
Yazyev O.
,
Austin A. J.
,
Cammi R.
,
Pomelli C.
,
Ochterski J. W.
,
Martin R. L.
,
Morokuma K.
,
Zakrzewski V. G.
,
Voth G. A.
,
Salvador P.
,
Dnnenberg J. J.
,
Dapprich S.
,
Daniels A. D.
,
Farkas O.
,
Foresman J. B.
,
Oritz J. V.
,
Cioslowski J.
,
Fox D. J.
2009
Gaussian Inc.
Wallingford, CT

Zaho Y.
,
Truhlar D. G.
2006
J. Phys. Chem.
110
5121 -
5129
** DOI : 10.1021/jp060231d**

Vosko S. H.
,
Wilk L.
,
Nusair M.
1980
Can. J. Phys.
58
1200 -
1211
** DOI : 10.1139/p80-159**

Perdew J. P.
1986
Phy. Rev. B
33
8822 -
8824
** DOI : 10.1103/PhysRevB.33.8822**

Marenich A. V.
,
Gramer C. J.
,
Truhlar D. G.
2009
J. Phys. Chem. B
113
6378 -
6396
** DOI : 10.1021/jp810292n**

Merrick J. P.
,
Moran D.
,
Radom L.
2007
J. Phy. Chem. A
111
11683 -
11700
** DOI : 10.1021/jp073974n**

Maiorov V. N.
,
Grippen G. M.
1994
J. Mol. Biol.
235
625 -
634
** DOI : 10.1006/jmbi.1994.1017**

Citing 'Conformations and Vibrational Frequencies of a Precursor of Benzovesamicol Analogues Studied by Density Functional Theories
'

@article{ JCGMCS_2014_v35n8_2311}
,title={Conformations and Vibrational Frequencies of a Precursor of Benzovesamicol Analogues Studied by Density Functional Theories}
,volume={8}
, url={http://dx.doi.org/10.5012/bkcs.2014.35.8.2311}, DOI={10.5012/bkcs.2014.35.8.2311}
, number= {8}
, journal={Bulletin of the Korean Chemical Society}
, publisher={Korean Chemical Society}
, author={Park, Jong-Kil
and
Choe, Sang Joon}
, year={2014}
, month={Aug}