- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Articles in Press ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents
ISRN Organic Chemistry
Volume 2012 (2012), Article ID 274980, 7 pages
Formation of 5-Hydroxy-3-methoxy-1,4-naphthoquinone and 8-Hydroxy-4-methoxy-1,2-naphthoquinone from Juglone
1Department of Biochemistry and Food Chemistry, University of Turku, Vatselankatu 2, 20014 Turku, Finland
2Department of Chemistry, University of Turku, Vatselankatu 2, 20014 Turku, Finland
Received 25 August 2012; Accepted 12 September 2012
Academic Editors: I. Alonso, N. Farfan, and F. Felluga
Copyright © 2012 Bastian Blauenburg et al. 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.
From the treatment of 5-hydroxy-1,4-naphthoquinone (juglone) with acetic anhydride and H2SO4 followed subsequently by treatment with methanolic HCl, 5-hydroxy-3-methoxy-1,4-naphthoquinone (3-methoxy juglone) and 8-hydroxy-4-methoxy-1,2-naphthoquinone were obtained as products rather than the anticipated product 2,5-dihydroxy-1,4-naphthoquinone (2-hydroxy juglone). The reaction and the identification of the products are discussed in terms of NMR and DFT calculations.
As part of our ongoing investigation  into the C-ribosylation of (R)-prealnumycin (1) by AlnA, an enzyme produced by Streptomyces albus (Scheme 1), it was decided to test various 1,4-naphthoquinone substrates against AlnA in order to identify the structural motifs necessary for, or inhibitory to, reaction.
Thus, we considered 2,5-dihydroxy-1,4-naphthoquinone (2, 2-hydroxy juglone, Figure 1) as a test substrate since it could potentially reveal aspects of the charge required to be present on C-3 for the reaction to proceed. Additionally, the facile preparation of 2,5-dihydroxy-7-methyl-1,4-naphthoquinone (2-hydroxy-7-methyl juglone, 3) from 5-hydroxy-7-methyl-1,4-naphthoquinone (7-methyl juglone, 4) had been reported  by treatment first with acetic anhydride and H2SO4 and then subsequently with methanolic HCl. After treating 5-hydroxy-1,4-naphthoquinone (5, juglone) as prescribed, we obtained not 2-hydroxy juglone (2) as anticipated, but rather 5-hydroxy-3-methoxy-1,4-naphthoquinone (6, 3-methoxy juglone) and 8-hydroxy-4-methoxy-1,2-naphthoquinone (7). Mahapatra et al. , in turn, seem to have taken their preparation from Lillie and Musgrave  who performed precisely the same reaction. Interestingly, whilst the former workers only reported the one product, namely, 2-hydroxy-7-methyl juglone (3), the latter pair reported the formation of both 3 (as the major product) and its regioisomer 3,5-dihydroxy-7-methyl-1,4-naphthoquinone (3-hydroxy-7-methyl juglone, 8). Mahapatra et al.  only reported 1H NMR chemical shifts () for their compound, and indeed only one OH signal—which they assigned as HO-5—with no carbon or correlation spectra reported. As best we can ascertain, it appears that Mahapatra et al.  based their structural assignment on the work of Lillie and Musgrave  and the of HO-5 (vide infra). This structural assignment therefore carries with it some concern, not because it is a labile proton per se, but rather because of its value in relation to the reported value . The correctness of this structural assignment is scrutinized in this work.
Both 3-methoxy juglone (6) and 8-hydroxy-4-methoxy-1,2-naphthoquinone (7) are known compounds, 6 in particular has been studied well synthetically [4, 5], is a widely spread naturally occurring compound present, for example, in a number of species of the genus Juglans [6–9], and possesses significant biological activity [7–9]. Compound 7 is unknown outside the laboratory [4, 10, 11], and overall there are only limited reports on it, though it too has potential biological activity of significance . Herein we report on the synthesis of these two compounds and discuss the reaction and identification of the products in terms of NMR and DFT calculations.
2. Results and Discussion
Our original intent was to follow analogously the preparation of Mahapatra et al.  who reported the synthesis of what they considered to be 2-hydroxy-7-methyl juglone (3) from 7-methyl juglone (4) and apply the same conditions to juglone (5) in order to obtain 2-hydroxy juglone (2). After performing the synthetic preparation, we isolated two major components in comparable amounts and readily identified them as 3-methoxy juglone (6) and 8-hydroxy-4-methoxy-1,2-naphthoquinone (7). Whilst fully expecting the first eluting compound from column chromatography over silica gel to be 2-hydroxy juglone (2) rather than 3-methoxy juglone (6), it was immediately evident that a methoxy group was present in the molecule by inspection of the 1H NMR spectrum (see Table 1), and furthermore, its general location was evident from the long-range coupling and an NOE to a vinyl-type proton. Although the methoxy group was then expected to be attached at C-2, the correlation of the vinyl-type proton (H-2) to C-8a in the HMBC spectrum was categorical in placing it at C-3 in addition to an NOE between the methoxy group protons and HO-5. Following these pivotal connectivities, the rest of the NMR analysis was consistent with the assigned structure. Turning our intention to the second isolate 7 eluting from the column, it exhibited some strongly divergent NMR spectral differences to 6 (e.g., , , etc.), and the reason quickly became apparent as this isomer of 6 was not a 1,4-quinoid structure, but rather a 1,2-quinoid structure. Notable correlations identifying this compound as 8-hydroxy-4-methoxy-1,2-naphthoquinone (7) were H-2 to C-8a and H-8 to C-1 (conspicuous by its of 168.15 ppm) in the HMBC spectrum and the NOE between the methoxy group and H-8.
Though the obtainment of the methoxy derivatives 3-methoxy juglone (6) and 8-hydroxy-4-methoxy-1,2-naphthoquinone (7) was unexpected, the presence of methoxy groups may be rationalized simply by methylation of the freed hydroxyls after hydrolysis (Scheme 2). The surprising occurrence of 7 can be accounted for by either of the following: if 4,8-dihydroxy-1,2-naphthoquinone (9) is not produced directly during the hydrolytic deacetylation concomitant with spontaneous oxidation of the tetraacetylated intermediate 10, then it arises through prototropic tautomerism from 3,5-dihydroxy-1,4-naphthoquinone (3-hydroxy juglone, 11) with which it is in equilibrium with in solution . This has been demonstrated by obtaining 7 as part of a mixture together with 6 when 2-hydroxy juglone (2) is methylated . Given the results of our reaction together with those of Mahapatra et al.  and Lillie and Musgrave , it seems that despite the deceptive simplicity of the reaction, it can be highly variable in terms of regiospecificity since sometimes attack at C-2 dominates , sometimes attack occurs exclusively at C-3 (this work), or sometimes attack seemingly occurs exclusively at C-2  (vide infra).
To first address the question of regiospecificity, DFT calculations revealed that 2-hydroxy juglone (2) is 1.40 kcal mol−1 more stable than 3-hydroxy juglone (11), similar results were obtained for 2-hydroxy-7-methyl juglone (3, 0.00 kcal mol−1) and 3-hydroxy-7-methyl juglone (8, 1.24 kcal mol−1) in consideration of the work of Mahapatra et al.  and Lillie and Musgrave , and furthermore, similar results were also obtained for 2-methoxy juglone (12, 0.00 kcal mol−1) and 3-methoxy juglone (6, 0.66 kcal mol−1) in consideration of the methoxy derivatives isolated in this study. Thus, in each case, the 2-substituted derivative was calculated to be the thermodynamic product, which is consistent with the results of Mahapatra et al.  and Lillie and Musgrave  but does not account for the observations herein. To consider if the C-2 position was more susceptible to nucleophilic attack, the atomic charges on C-2 and C-3 were examined: in juglone (5), the Mulliken charge on C-2, , was calculated to be −0.182 cf. at −0.176. These differences did not differ appreciably for 7-methyl juglone (4) where charges of −0.184 and −0.175, respectively, were found for and . These charge differences are considered inconsequential to the course of the reaction, either between 2- and 3-substituted products or for differences observed between 7-methyl juglone (4) and juglone (5).
Kinetically however, it would appear that, based on the energetic preference by 5.24 kcal mol−1 for the intermediate in the case of attack by acetic anhydride on protonated juglone (5H) at the C-3 position over the C-2 position (Scheme 3), 3-hydroxy juglone (11) would be the preferred product after oxidative hydrolysis. Interestingly, attack by acetate to yield the enolate ion provides the opposite effect, and then some with attack at the C-2 position preferred over the C-3 position by 9.39 kcal mol−1; that is, 2-hydroxy juglone (2) would be expected as the preferred end-product. Since the reaction was conducted under acid-catalyzed conditions and is likely to be irreversible, it is concluded that 3-hydroxy juglone (11) should be the dominant product, essentially in concert with observations in this work since 3-substituted products (6 and 7) were obtained. Kinetically one can anticipate similar results for 7-methyl juglone (4) as for the reaction of juglone (5). However it is clear that there must be a fine balance between the competing pathways, and even slight perturbation of the conditions could lead to very different results being obtained.
To gain further insight into these structures, DFT calculations of their 1H and 13C NMR chemical shifts ( and , resp.) were performed (see Table 2) as these have proven to be of considerable assistance in structural analysis [14–16]. Interestingly, there is nothing in the calculated of the nonlabile protons in 2-hydroxy juglone (2) and 3-hydroxy juglone (11) to categorically distinguish between the two structures. The same holds for the calculated of the nonlabile protons in 2-hydroxy-7-methyl juglone (3) and 3-hydroxy-7-methyl juglone (8) and for the methoxy derivatives 2-methoxy juglone (12) and 3-methoxy juglone (6), with the latter pair even more similar. Thus, to reliably distinguish between the two compounds in any of these pairs based on the of the nonlabile protons, one must be in possession of both compounds. The exception is for the exchangeable HO-5 protons where significant differences are clearly evident. Usually though, labile protons are considered unreliable for assignment purposes based on experimental however, the experimental of OH protons can, if certain conditions are met, namely conditions conducive to slow exchange of the protons, be treated as reliable for assignment purposes . With respect to calculated , although the chemical shifts of labile protons are generally problematic to calculate in absolute terms, they are normally reliably found in their correct chemical shift order. When exchangeable protons are strongly hydrogen bonded though, they can be accurately calculated, and this is clearly evident in the case of 2 and 11 (experimental values of 12.28 ppm and 11.04 ppm, resp. ) and for 6. Thus, HO-5 is the only proton with a decisive indication of isomer identification in these pairs of compounds.
Mahapatra et al.  only observed one OH signal, assigned as the HO-5 proton, in which case the other OH signal must therefore be broad due to exchange despite the possibility of it too being intramolecularly hydrogen bonded. Therefore, under such conditions the reliability of the observed OH proton must be treated with caution. However, one inference can be made: the reported  value for , 11.69 ppm, is intermediate between values reported  for 2-hydroxy-7-methyl juglone (3), 12.23 ppm, and 3-hydroxy-7-methyl juglone (8), 10.99 ppm. If a mixture of 3 and 8 was present, then the observed intermediate value may be a consequence of intermolecular proton exchange. Of note, exchange can be appreciable on the NMR timescale at the observational frequency of 200 MHz used by Mahapatra et al. . A mixture of 3 and 8 would also exhibit a depressed m.p., and indeed a depressed m.p. of 10°C below that of the lower melting point isomer 3 was in fact reported  by them. Since there are insufficient differences in the calculated of the nonlabile protons to decisively affect isomer identification of the pure compounds with only one isomer on hand, one cannot conclude from the reported  which compound(s) was present in the sample. Indeed, for a mixture, differences might only be seen with decent resolution which could easily be lost by even slightly poor field homogeneity and which would be problematic in any case at 200 MHz. Mahapatra et al.  only described H-6 and H-8 as singlets, whilst we observed multiplets for the analogous protons in 7-methyl juglone (4). Furthermore, they did not do more to determine which compound they had or even if they had a mixture of compounds (e.g., IR and HRMS are insufficient for this purpose in the current case). Therefore, although we cannot unequivocally confirm that Mahapatra et al.  have the wrong structure(s) for the compound they claim as 2-hydroxy-7-methyl juglone, we must conclude that they had, most probably, a mixture of 3 and 8.
The DFT-calculated 13C NMR chemical shifts are presented in Table 3 wherein it can be seen that for each isomeric pair, four carbons are able to provide substantial differentiation for assignment purposes. There is clearly no problem in identifying our 3-methoxy juglone (6) as the correct structure based on comparison of the calculated for 2-methoxy juglone (12) and 6, in particular, , , and . Even assignment is not necessary as , , and are all well dispersed from one another so one only needs to look for chemical shift differences () in the appropriate region. Of special note though is the between and , which has been calculated as 11.02 ppm for 12 whilst for 6 it is 2.65 ppm the latter comparing rather well to the observed value of 0.99 ppm. At the very least, 13C NMR acquisition, if not affording isomer identification (correlation spectra notwithstanding) outright, readily alludes to the presence of isomers in a sample and should have been the method of choice in the study of Mahapatra et al. .
The treatment of 5-hydroxy-1,4-naphthoquinone (5, juglone) with acetic anhydride under acidic conditions can lead to variable results as the course of the reaction is very dependent on the exact conditions applied. Herein the treatment of juglone (5) with acetic anhydride and H2SO4 and then subsequent treatment with methanolic HCl led to the isomeric compounds 5-hydroxy-3-methoxy-1,4-naphthoquinone (6, 3-methoxy juglone) and 8-hydroxy-4-methoxy-1,2-naphthoquinone (7) being obtained instead of the desired 2-hydroxy juglone (2). To correctly identify the products of a reaction in which regiospecificity may be variable or lacking altogether, it is important that appropriate analytical methods be applied, for example, 13C NMR, correlation spectra, or molecular modeling, to properly identify isomers or at least ascertain the presence of other isomers in the sample. In the case of comparison to the literature data, care needs to be taken that conditions are compliant with the measurements being undertaken to ensure correct interpretation of the results. Certainly, reliance on more than just one value for identification should be exercised for the dangers of not doing so are self-evident.
NMR spectra were acquired using a Bruker Avance NMR spectrometer equipped with either a 5 mm inverse or a 5 mm normal configuration probe, both with z-gradient capability, at a field strength of 11.75 T operating at 500 and 125 MHz for 1H and 13C nuclei, respectively, at 25°C with samples contained in CDCl3. The chemical shifts of 1H and 13C nuclei are reported relative to TMS incorporated as an internal standard ( ppm for both 1H and 13C). General NMR experimental details for 1D 1H, 13C, and DEPT and standard gradient-selected 2D DQF-COSY, HSQC, and HMBC spectra have been previously described [17–20] for routine structural determinations and chemical shift assignments. For 1D selective COSY (set on a value of 7 Hz for vicinal coupling whilst for long-range , values varied in the range 0.25–1 Hz) and NOESY [21, 22] (mixing times: 0.5, 0.75 s) experiments, selective excitation was effected using a 180° 50 ms Gaussian-shaped pulse via excitation sculpting . Spin analysis was performed using the Perch iteration software [12, 13] for the extraction of and .
4.2. Molecular Modeling
DFT quantum chemical calculations were performed using Gaussian 09  (version A.01) and analyzed using GaussView (version 3.07). Geometry optimization of the structures in the gas phase was performed using the M06-2X hybrid metadensity functional [25, 26] with the 6-31G(d) basis set in tandem with vibrational analysis and thermochemistry calculations at the same level of theory. Vibrational analyses, invoking the keyword freq = noraman, were conducted to confirm that optimized structures were true minima on the potential energy surface by not providing imaginary frequencies and to obtain the thermodynamic contributions at 298.15 K and 1 atm wherein frequencies were left unscaled. If necessitated by the presence of imaginary frequencies, structures were re-optimized (together also with vibrational analyses and thermochemistry calculations) using tight convergence criteria by invoking the keywords opt = tight and int = ultrafine.
Absolute chemical shieldings were calculated for geometry-optimized structures in the gas phase using the GIAO method  and the B3LYP functional [28, 29] with the cc-pVTZ basis set. Chemical shifts were determined by subtracting calculated shieldings from the calculated shieldings of the reference compound TMS for which both and ppm. The chemical shifts of 1H and 13C nuclei were calibrated following literature methodology  using the following relationships:
4.3. Synthetic Preparation of 3-Methoxy Juglone (6) and 8-Hydroxy-4-methoxy-1,2-naphthoquinone (7)
Following the method of Mahapatra et al. , a mixture of juglone (0.6 mmol), acetic acid anhydride (3 mL), and concentrated H2SO4 (0.1 mL) was kept overnight and then poured onto ice. The reaction mixture was extracted with chloroform, dried with MgSO4, and then filtered. The residue obtained after removal of the solvent in vacuo was crystallized from chloroform-hexane and the solid material collected by filtration, taken up in methanolic HCl (2 M), and then refluxed for 1 hour. After extraction of the compound using diethyl ether, column chromatography over silica gel (CHCl3–CH3OH, 95 : 5) resulted in early fractions containing 6 followed by later fractions containing 7. Analysis of the reaction products was accomplished by HPLC (SunFire C18 3.5 μm, 2.1 × 150 mm) equipped with a UV-Vis DAD and using a gradient of 15–100% acetonitrile in water with 0.1% formic acid whereby 7 was found to elute before 6. Both compounds were subjected to semipreparative HPLC (SunFire Prep C18 5 μm, 10 × 250 mm) for final purification using the same gradient of 15–100% acetonitrile in water with 0.1% formic acid. UV-Vis: for 6, 281, 411 nm; for 7, 286, 422 nm. For 1H and 13C NMR data of 6 and 7, see Table 1.
The Academy of Finland (Grant nos. 136060 and 251486) is thanked for financial support and the CSC-IT Center for Science Ltd. is thanked for providing computational resources.
- T. Oja, K. D. Klika, L. Appassamy et al., “Biosynthetic pathway toward carbohydrate-like moieties of alnumycins contains unusual steps for C-C bond formation and cleavage,” Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 16, pp. 6024–6029, 2012.
- A. Mahapatra, S. P. N. Mativandlela, B. Binneman et al., “Activity of 7-methyljuglone derivatives against Mycobacterium tuberculosis and as subversive substrates for mycothiol disulfide reductase,” Bioorganic and Medicinal Chemistry, vol. 15, no. 24, pp. 7638–7646, 2007.
- T. J. Lillie and O. C. Musgrave, “Ebenaceae extractives. Part 7. Use of hydroxy-proton shifts of juglone derivatives in structure elucidation,” Journal of the Chemical Society, Perkin Transactions 1, no. 4, pp. 355–359, 1977.
- H. Laatsch, “Methylation of 3,-dihydroxylated 2,-binaphtho-1,4;,-quinones,” Zeitschrift für Naturforschung B, vol. 45, no. 3, pp. 393–400, 1990.
- G. Wurm and U. Geres, “1,4-Naphthoquinones, XVII: studies on the synthesis of 3-hydroxy-juglone derivatives,” Archiv der Pharmazie, vol. 322, no. 3, pp. 155–157, 1989.
- H.-Y. Yu, X. Li, F.-Y. Meng, H.-F. Pi, P. Zhang, and H.-L. Ruan, “Naphthoquinones from the root barks of Juglans cathayensis Dode,” Journal of Asian Natural Products Research, vol. 13, no. 7, pp. 581–587, 2011.
- H. Ruan and H. Yu, “Application of two naphthoquinone compounds in preparing the medical preparations for treating neoplasm,” Chinese Patent, CN 102552226 A 20120711, 2012.
- W. Tang and F. Chen, “Preparation of antitumor agent from Juglans regia,” Chinese Patent, CN 101185671 A 20080528, 2008.
- M.-l. Sun, Z.-q. Song, and G.-z. Fang, “Antifungal activities and active components of alcoholic extracts from leaves of Juglans mandshurica Maxim,” Linchan Huaxue Yu Gongye, vol. 27, pp. 81–84, 2007.
- K. Ditrich, G. Hamprecht, E. Ammermann, G. Lorenz, and H. Laatsch, “Preparation of 1,2-naphthoquinones as fungicides,” German Patent, DE 3926747 A1 19910214, 1991.
- G. Wurm and U. Geres, “1,4-Naphthoquinones, X: reactions of the methyl ethers of 1,5-dihydroxynaphthalene with singlet oxygen,” Archiv der Pharmazie, vol. 318, no. 10, pp. 931–937, 1985.
- Peak Research NMR Software, Perch Solutions, Kuopio, Finland, 2003, http://www.perchsolutions.com.
- R. Laatikainen, M. Niemitz, U. Weber, J. Sundeun, T. Hassinen, and J. Vepsäläinen, “General strategies for total-lineshape-type spectral analysis of NMR spectra using integral-transform iterator,” Journal of Magnetic Resonance, Series A, vol. 120, no. 1, pp. 1–10, 1996.
- S. Böhm, J. Tomaščiková, J. Imrich et al., “Computational study to assign structure, tautomerism, E/Z and s-cis/s-trans isomerism, π-delocalization, partial aromaticity, and the ring size of 1,3-thiazolidin-4-ones and 1,3-thiazin-4-ones formed from thiosemicarbazides,” Journal of Molecular Structure, vol. 916, no. 1–3, pp. 105–118, 2009.
- P. Tähtinen, A. Bagno, K. D. Klika, and K. Pihlaja, “Modeling NMR parameters by DFT methods as an aid to the conformational analysis of cis-fused 7a(8a)-methyl octa(hexa)hydrocyclopenta[d][1,3]oxazines and [3,1]benzoxazines,” Journal of the American Chemical Society, vol. 125, no. 15, pp. 4609–4618, 2003.
- K. Pihlaja, P. Tähtinen, K. D. Klika, T. Jokela, A. Salakka, and K. Wähälä, “Experimental and DFT 1H NMR study of conformational equilibria in trans-4′,7-dihydroxyisoflavan--ol and trans-isoflavan-4-ol,” Journal of Organic Chemistry, vol. 68, no. 18, pp. 6864–6869, 2003.
- J. Mäki, P. Tähtinen, L. Kronberg, and K. D. Klika, “Restricted rotation/tautomeric equilibrium and determination of the site and extent of protonation in bi-imidazole nucleosides by multinuclear NMR and GIAO-DFT calculations,” Journal of Physical Organic Chemistry, vol. 18, no. 3, pp. 240–249, 2005.
- P. Virta, A. Koch, M. U. Roslund et al., “Synthesis, characterisation and theoretical calculations of 2,6-diaminopurine etheno derivatives,” Organic and Biomolecular Chemistry, vol. 3, no. 16, pp. 2924–2929, 2005.
- K. D. Klika, J. Bernát, J. Imrich, I. Chomča, R. Silanpää, and K. Pihlaja, “Unexpected formation of a spiro acridine and fused ring system from the reaction between an N-acridinylmethyl-substituted thiourea and bromoacetonitrile under basic conditions,” Journal of Organic Chemistry, vol. 66, no. 12, pp. 4416–4418, 2001.
- E. Balentová, J. Imrich, J. Bernát et al., “Stereochemistry, tautomerism, and reactions of acridinyl thiosemicarbazides in the synthesis of 1,3-thiazolidines,” Journal of Heterocyclic Chemistry, vol. 43, no. 3, pp. 645–656, 2006.
- J. Stonehouse, P. Adell, J. Keeler, and A. J. Shaka, “Ultrahigh-quality NOE spectra,” Journal of the American Chemical Society, vol. 116, no. 13, pp. 6037–6038, 1994.
- K. Stott, J. Stonehouse, J. Keeler, T.-L. Hwang, and A. J. Shaka, “Excitation sculpting in high-resolution nuclear magnetic resonance spectroscopy: application to selective NOE experiments,” Journal of the American Chemical Society, vol. 117, no. 14, pp. 4199–4200, 1995.
- T.-L. Hwang and A.-J. Shaka, “Water suppression that works. Excitation sculpting using arbitrary waveforms and pulsed field gradients,” Journal of Magnetic Resonance, Series A, vol. 112, no. 2, pp. 275–279, 1995.
- M. J. Frisch, G. W. Trucks, H. B. Schlegel et al., Gaussian09, Revision A.01, Gaussian, Wallingford, Con, USA, 2009.
- Y. Zhao and D. G. Truhlar, “The M06 suite of density functionals for main group thermochemistry, thermochemical kinetics, noncovalent interactions, excited states, and transition elements: two new functionals and systematic testing of four M06-class functionals and 12 other functionals,” Theoretical Chemistry Accounts, vol. 120, no. 1–3, pp. 215–241, 2008.
- Y. Zhao and D. G. Truhlar, “Density functionals with broad applicability in chemistry,” Accounts of Chemical Research, vol. 41, no. 2, pp. 157–167, 2008.
- K. Wolinski, J. F. Hinton, and P. Pulay, “Efficient implementation of the gauge-independent atomic orbital method for NMR chemical shift calculations,” Journal of the American Chemical Society, vol. 112, no. 23, pp. 8251–8260, 1990.
- A. D. Becke, “Density-functional thermochemistry. III. The role of exact exchange,” Journal of Chemical Physics, vol. 98, no. 7, pp. 5648–5652, 1993.
- C. Lee, W. Yang, and R. G. Parr, “Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density,” Physical Review B, vol. 37, no. 2, pp. 785–789, 1988.