Next Article in Journal
Insecticidal Activities of Bark, Leaf and Seed Extracts of Zanthoxylum heitzii against the African Malaria Vector Anopheles gambiae
Next Article in Special Issue
Effects of High Hydrostatic Pressure on Escherichia coli Ultrastructure, Membrane Integrity and Molecular Composition as Assessed by FTIR Spectroscopy and Microscopic Imaging Techniques
Previous Article in Journal
A Sensitive A3B Porphyrin Nanomaterial for CO2 Detection
Previous Article in Special Issue
Comparative Study of the Structural and Vibroelectronic Properties of Porphyrin and Its Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Largely Reduced Grid Densities in a Vibrational Self-Consistent Field Treatment Do Not Significantly Impact the ResultingWavenumbers

1
Institute of Analytical Chemistry and Radiochemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, Innsbruck 6020, Austria
2
Institute of General, Inorganic and Theoretical Chemistry, Center for Chemistry and Biomedicine, University of Innsbruck, Innrain 80-82, Innsbruck 6020, Austria
*
Author to whom correspondence should be addressed.
Molecules 2014, 19(12), 21253-21275; https://doi.org/10.3390/molecules191221253
Submission received: 20 October 2014 / Revised: 25 November 2014 / Accepted: 9 December 2014 / Published: 17 December 2014
(This article belongs to the Special Issue Advances of Vibrational Spectroscopic Technologies in Life Sciences)

Abstract

:
Especially for larger molecules relevant to life sciences, vibrational self-consistent field (VSCF) calculations can become unmanageably demanding even when only first and second order potential coupling terms are considered. This paper investigates to what extent the grid density of the VSCF’s underlying potential energy surface can be reduced without sacrificing accuracy of the resulting wavenumbers. Including single-mode and pair contributions, a reduction to eight points per mode did not introduce a significant deviation but improved the computational efficiency by a factor of four. A mean unsigned deviation of 1.3% from the experiment could be maintained for the fifteen molecules under investigation and the approach was found to be applicable to rigid, semi-rigid and soft vibrational problems likewise. Deprotonated phosphoserine, stabilized by two intramolecular hydrogen bonds, was investigated as an exemplary application.

Graphical Abstract

1. Introduction

The continuously increasing availability of computational resources as well as the development of accurate and efficient quantum chemical approaches have made computational vibrational spectroscopy an indispensable field complementing experimental techniques. Nowadays, almost every quantum chemical software package enables the analysis of second-order properties of the energy (e.g., IR and Raman absorptions, IR intensities, electric dipole polarizabilities, nuclear magnetic resonance chemical shifts, spin-spin coupling constants). Arguably, the prediction of IR and Raman absorptions is among the most important applications considering the prevalence of these technique. While most fundamental IR absorptions can be assigned quite satisfyingly on an empirical basis or via normal coordinate analysis, especially the near-infrared region is cluttered by overtones and combination excitations that are cumbersome to assign.
With regard to computational approaches, the harmonic oscillator approximation (HOA) is the most fundamental technique for obtaining vibrational spectroscopic data. However, due to the rigorous assumption introduced, in that the bond potential exhibits a harmonic shape, significant deviation from experiment is frequently observed. The most simple solution accounting for the lack of anharmonicities is to introduce empirical scaling factors that are multiplied with the harmonically approximated absorptions. Scott and Radom derived scaling factors for a vast number of ab initio and semi-empirical methods and a large number of basis sets [1,2]. It has to be stressed, however, that an empirical scaling may not be applied on a system-independent basis, even though the scaling factors have been derived for a rather large set of molecules. Noteworthy, a harmonic bond potential cannot be assumed a proper basis for considerations towards excitations involving more than one quantum of energy.
Especially due to these fundamental deficiencies, further corrective techniques have been developed that account for anharmonicities in an explicit manner. The two most prominent approaches are the vibrational self-consistent field (VSCF) method [3,4,5,6,7,8] and the vibrational second order perturbative ansatz (VPT2) [9,10]. VPT2 is a technique that relies on the computation of higher-order derivatives of the energy. Gaussian [11] is probably the most prominent commercial software package incorporating a VPT2 algorithm which relies on third and semi-diagonal fourth derivatives of the energy with respect to the nuclear coordinates. This method is applied regularly [12,13] but it may only be employed safely to vibrational problems where the harmonic part of the potential is dominant. Moreover, the computationally demanding generation of third and fourth derivatives limits the approach to rather small systems. VSCF on the other hand utilizes a separability ansatz which enables the anharmonicity of each mode to be accounted for by screening the respective potential in a point-wise manner. Mode interactions may be evaluated by computing d dimensional grids and the maximum possible dimensionality is determined by the number of vibrational degrees of freedom (N ) featured by the molecule under investigation. In a conventional VSCF calculation, each mode’s potential is characterized by 16 grid points [14]. In the past, scientific papers have been reported that employed grids of higher (r = 32) [15] and lower resolutions (r = 8) [16,17,18]. Roy et al. [19] recently presented data based on a variety of grid densities ranging from 8 to 16 points. Since they reported overall deviations summed over a set of reference molecules, one could conclude that this rather diverse choice of grid resolutions does not impose a significant error on the wavenumbers. However, to date there still seem to exist ambiguities with regard to a generally applicable grid density that is computationally feasible. Hence, and considering the massive computational effort involved in high resolution VSCF calculations, it seemed promising to conduct a systematic study addressing these issues. 15 reference molecules are investigated in detail and an application to the lowest energy conformer of phosphoserine is presented. At this point it seems worth mentioning that while many quantum chemical software packages such as GAMESS [20], NWChem [21] and MOLPRO [22] incorporate VSCF implementations, GAMESS is the program employed in this work due to the fact that it is freely available and compatible with the most prominent computer operating systems.

2. Methods

In this section, a brief overview describing the procedures underlying a VSCF evaluation is given. Particularities of the involved techniques are presented while reference to the original literature is given for further details.

2.1. Energy Minimization

A prerequisite for every computational spectroscopic analysis is a proper structural ensemble. For the majority of cases, the user would want to obtain absorption data corresponding to an equilibrium structure and thus, an energetically favorable reference geometry is required. This implies that at least a local energy minimum (or the global energy minimum) is obtained and this is reflected by a nearly zero gradient of the energy g(q)∆q with respect to the nuclear coordinates q. While the program’s default criteria suffice for energetical and structural considerations in many cases, a VSCF calculation requires a thoroughly minimized geometry. This condition is realized by setting up rigorous cutoff tolerances for g(q)∆q. We considered a molecular geometry properly minimized when each gradient contribution is smaller than 0.000001 Eh· Bohr−1. Every calculation reported herein was performed at the Møller-Plesset level of theory, accounting for electron correlation effects via a perturbative ansatz [23,24,25]. The frozen-core approximation [26] was employed, explicitly correlating all but the core electrons. Dunning’s correlation-consistent polarized basis set of triple-ζ quality has been utilized since it can be considered a routinely employed set of functions that is known to deliver results of acceptable accuracy [27]. The appropriate point group was imposed onto the molecular geometry where applicable in order to make use of the efficient symmetry optimized SCF algorithm in GAMESS.

2.2. Second Derivatives of the Energy

A potential V(q) may be expressed as a Taylor series
V ( q ) = V 0 ( q ) + g ( q ) Δ q + 1 2 Δ q T Δ q +
and for the case that the reference geometry resembles an energy minimum and assuming a harmonic potential shape, all but the first and third terms vanish. This simplified approach enables the evaluation of the elements of ℍr,s by analytical or numerical analysis:
r , s = 2 V ( q ) q r q s
A transformation into mass-weighted coordinates and diagonalization of the Hessian matrix ℍr,s yields the eigenvectors corresponding to the normal modes. At this point, computational spectroscopic data at the HOA level is available. It seems worth mentioning that the further description of extended vibrational analysis is based on mass-weighted normal coordinates or other non-redundant coordinates for the sake of clarity, whereas the expansion of the potential energy surfaces was realized in rectilinear coordinates (Section 3.1). VSCF calculations require a Hessian corresponding to the equilibrium geometry, as will be discussed in the next paragraph. The Hessian matrices in this work have been obtained semi-numerically involving two displacements of ±0.01 Bohr about the reference geometry for each atom.

2.3. The VSCF Routine

Even though VSCF accounts for anharmonicities explicitly, it still requires ℍr,s as a reference state. This is owed to the fact, that VSCF assumes a molecule’s vibrational wave function Ψ(Q1,..,N) to be separable into single mode wave functions ψ i ( n ) :
Ψ ( Q 1 , , N ) = i = 1 N ψ i ( n ) ( Q i )
While this product ansatz can be implemented rather efficiently, it inherently limits the accuracy of the VSCF technique since it represents decoupled vibrational modes. There are, besides the methods mentioned in Section 2.4, techniques available which are able to circumvent this limitation [28,29]. Introducing the variational principle [30,31], the single mode VSCF equation [6,32] is formulated as:
[ 1 2 2 Q i 2 + V i ( n ) ¯ ( Q i ) ] ψ i ( n ) ( Q i ) = ε i ( n ) ψ i ( n ) ( Q i )
where V i ( n ) ¯ ( Q i ) is mode Qi’s effective potential. Equation (4) neglects parts of the full Watson Hamiltonian for non-linear molecules [33] which accounts for N vibrational degrees of freedom:
Ĥ = 1 2 i = 1 N 2 Q i 2 + V ( Q 1 , , N ) + 1 2 α β π ˆ α μ α β π ˆ β 1 8 α μ α α
The second and third part of the Hamiltonian represent the generalized inverse of the effective moment of inertia µαβ and the vibrational angular momentum operator π ˆ α . Given that the full Watson Hamiltonian is cumbersome to implement computationally, most implementations resort to neglecting the last two terms of Ĥ. The error introduced is small especially for fundamental excitations as has been shown by Bowman et al. [34].
The VSCF equations are solved in an iterative manner until self-consistency is achieved. Illustratively, VSCF may be seen as an approach where each vibration is influenced by the mean field of the surrounding modes. However, especially when certain mode interactions exhibit unique properties, this mean field ansatz does no longer properly describe the vibrational problem sufficiently well. A number of post-VSCF techniques accounting for this deficiency have been presented and they are discussed in the subsequent section.
Since the VSCF technique involves a grid-based screening of bond potentials and their interactions, a discussion of the so-called hierarchical expansion is justified. For a system exhibiting N normal modes, the VSCF expansion is denoted as:
V ( Q 1 , , N ) = i = 1 N V i d i a g ( Q i ) + i < j N V i , j p a i r s ( Q i , Q j ) + i < j < k N V i , j , k t r i p l e s ( Q i , Q j , Q k ) +
A complete interaction scheme is realized by imposing N -dimensionality onto Equation (6) but obviously, such a treatment scales unfavorably with the studied system size. Noteworthy, most applications of VSCF theory truncate the expansion largely by including only i = 1 N V i d i a g ( Q i ) and i < j N V i , j p a i r s ( Q i , Q j ) contributions [35,36,37,38]. This simplification is valid for a large number of applications and the inclusion of higher-order interaction terms is hardly justifiable for other than the smallest molecules or for cases, where convergence is only achieved when including such terms. Programs such as MULTIMODE [34,39], Molpro [22] and MIDAScpp enable coupling orders of d > 3 yielding IR bands of exceptional accuracy [7,40]. In this article, we confine our discussion to i = 1 N V i d i a g ( Q i ) and i < j N V i , j p a i r s ( Q i , Q j ) contributions since this study aims at larger molecules where an incorporation of higher coupling terms would be prohibitive. The number of grid points due (Np) is calculated via Equation (7) .
N p = r × N + N × ( N 1 ) 2 × r 2
Assuming r = 16 and taking glycine as an example, the pair-wise approximation would require 71,040 points to be computed while a VSCF calculation involving also three-mode interactions already gives rise to 8,361,344 energy evaluations, a nearly 120 fold increase in computational burden. Herein, we will examine the impact of the grid density on the quality of the absorption data by conducting VSCF calculations involving between 6 and 16 grid points per mode. For the displacements underlying the PES scan, a symmetric grid range of [−4 ω i 0.5 , +4 ω i 0.5 ], with ωi being the harmonic frequency of the i-th normal mode, was chosen. The grid points within these boundaries have been set-up in an equidistant manner which implies that for even grid densities (i.e., 6, 8, 10, 12, 14 and 16 points), equilibrium is located between the two innermost grid points. For odd grid densities (i.e., 7, 9, 11, 13 and 15 points), equilibrium is described by one distinct grid point. For every resolution r < 16, the grid points are interpolated to the original resolution of r = 16 by a polynomial fit.

2.4. Extensions of the VSCF Approach

As mentioned earlier, the VSCF assumption does not hold for many applications. Therefore, a perturbative approach has been suggested, accounting for the error introduced by the mean-field VSCF ansatz [15,41,42]. A prerequisite is that the difference between the "true" energy and the mean-field VSCF energy is small. An expansion known from electronic structure theory [23] is then formulated for the energy and wave function of a state n that is truncated at second order:
E n E n 0 + λ E n 1 + λ 2 E n 2
Ψ n Ψ n 0 + λ Ψ n 1 + λ 2 Ψ n 2
with H = H n 0 + λ∆V and an insertion of En and Ψn in H Ψn = EnΨn, the energy contribution at second order may be formulated as:
E n 2 = m n Ψ n 0 | Δ V | Ψ m 0 Ψ m 0 | Δ V | Ψ n 0 E n 0 E m 0
Herein, Ψ n 0 and Ψ m 0 are the unperturbed vibrational product wave functions and E n 0 and E m 0 are the unperturbed VSCF energies. This technique, which is commonly abbreviated as second order perturbation theory augmented (PT2)-VSCF or correlation-corrected (CC)-VSCF, has proven useful for many applications. However, since obtaining correlation-corrected VSCF energies requires significant computational resources, the basic PT2-VSCF approach is limited to small systems. Due to the inherent need for an efficient solution of Equation (10), Gerber’s work group developed a solution based on the assumption of orthogonal vibrational single mode wave functions (Equation (3)) which leads to the annihilation of diagonal elements in Equation (10) [43]. Tests involving Glyn peptides showed that the runtime improvement is greater for larger molecules: for monomeric Gly, a speedup factor close to 6 was observed while the correlation-corrected wavenumbers of tetraglycine have been evaluated more than 16 times faster as with the conventional PT2-VSCF ansatz. This acceleration technique, and also the fact that this accelerated PT2-VSCF technique is readily implemented in GAMESS [20], are probably the main reasons for PT2-VSCF being a routinely employed VSCF correction. It has to be stressed, however, that the evaluation of the VSCF equations with its further corrections may not be confused with the preceding and highly demanding evaluation of the potential energy grid (vide supra).
Problems may arise when degenerate vibrational states are present. PT2-VSCF can fail here due to a close to zero denominator in Equation (4) which can lead to a largely overestimated perturbative correction. Hence, the degenerate PT2-VSCF method has been developed [44] but to date it has been implemented in GAMESS exclusively for degeneracies arising from fundamental excitations. A more generally available simple solution available in the GAMESS program code is that contributions involving denominators falling below a critical value are excluded from the treatment. Most applications reporting PT2-VSCF derived data rely on this simplification and therefore, we will also confine our calculations to this simplified technique.
Besides perturbation theory augmented VSCF, significant effort is put into post-VSCF techniques involving configuration interaction [45,46,47,48], coupled-cluster [49,50,51] and multi-configurational SCF theory [52,53]. While such methods yield highly accurate data, they are to date only applicable to small vibrational problems with less than 20 atoms.
Due to the popularity of the PT2-VSCF method and the fact that results of good quality at manageable computational cost are available also for larger molecules, all data presented in this paper is corrected exclusively with this perturbative ansatz.

3. Results and Discussion

The fifteen molecules under investigation gave rise to 176 distinct vibrational degrees of freedom, considering all fundamental stretching, deformation and torsional vibrations. PT2-VSCF can yield questionable or sometimes even divergent results for very low-lying and floppy torsions, which is due to the fact that the PES expansion is conventionally carried out in Cartesian coordinates [18,19,54]. The error induced through an anharmonic correction of such a vibration can exceed the boundaries of accuracy known from the HOA [54]. Hence, the usual procedure is to either treat such vibrations harmonically or to describe the underlying displacements in internal coordinates. A substitution for physically more meaningful internal coordinates was proposed by Njegic and Gordon [54] and could be shown to yield good results for formamide and thioformamide [55] as well as for H2O2 but by introducing a new expansion technique for the kinetic energy operator [56]. Nonetheless, setting up internal coordinates that properly describe the displacements underlying a VSCF treatment is by no means a trivial task. The GAMESS code is able to identify each normal mode’s contributions to particular internal coordinates [57], but the user still has to input a balanced description of each vibrational degree of freedom which can become unmanageably difficult for larger systems. Hence, and especially when large molecules are investigated, the majority of users of VSCF theory resort to a PES expansion in Cartesian coordinates and the contributions of critical torsions are omitted when self-consistency is not achieved.
For glycine (C2H5NO2), three normal modes (i.e., the N-Cα-Ccarb-O torsion, the NH2 group torsion and δCcarbO2,oop) had to be excluded from the VSCF treatment since they are known to lead to divergence during a perturbation theory corrected VSCF evaluation [19]. Similarly, for methanol (CH3OH) the CO axis torsion has been omitted from the PT2-VSCF treatment due to an inadequate description of this particularly floppy torsional mode within the VSCF framework [58]. Dimethylether also exhibits two floppy torsions involving the C-O axes that have been excluded likewise. For ethane, one normal mode near 290 cm−1 [59,60,61,62] was omitted due to its floppy character but the other six missing modes did arise from degenerate states in νCH3,as, δCH3,as and ρCH3. Experimentally, these degeneracies cannot be distinguished and since the PT2-VSCF derived values did not exhibit significant numerical discrepancies, each pair of degeneracies is presented as a single mean value for the sake of visibility.

3.1. Performance of the PT2-VSCF Approach

The computationally obtained absorption data are compared to experimental data in Table 1. The column headers indicate the employed number of grid points during the VSCF evaluations. As a measure of quality, the mean absolute percentage error µ for each molecule and each grid density is calculated according to Equation (11):
μ = 100 N i = 1 N | E x p t i V S C F i E x p t i |
The obtained values for µ are generally smaller than 2%, which is in good agreement with the recent work by Roy et al. [19] who concluded that MP2/cc-pVTZ based PT2-VSCF evaluations exhibit a mean unsigned error of under 2%. MP2 is a topical ab initio method that accounts for electron correlation effects to a certain extent and it seems as if for most molecules, this method indeed delivers satisfactory results. Acetonitrile is an exception due to its CN triple bond. However, the particularities of this molecule are discussed elsewhere [63] and it was found that state-of-the-art quantum chemical methods [64,65,66,67,68,69] are required for a proper description of νCN. Importantly, the more or less ubiquitous stretching motions arising from a methyl group are not described in a reliable manner and this has been recently ascribed to the nature of the MP2 method [19]. Conversely, it was shown that when higher order coupling terms (i.e., i < j < k N V i , j , k t r i p l e s ( Q i , Q j , Q k ) ) are included in a DPT2-VSCF [44] treatment, MP2 indeed seems to be a viable ab initio method of choice [63]. VSCF data relies on the HOA and hence makes use of the rigid rotor approximation [70] which inherently is co-responsible for discrepancies between experiment and theory. Considering that still a large number of approximations are involved even in state-of-the-art VSCF calculations and their underlying ab initio energy evaluations, it must be concluded that computationally derived spectroscopic data are prone to a certain degree of fortuitous error compensation. However, both the data presented by Roy et al. [19] and our results indicate that triple-ζ based MP2 calculations in a VSCF treatment rather reliably deliver data with no more than 2% of unsigned error.
Table 1. Computed and experimental absorption data and the calculated values for µ in %. n.a. means “not applicable” and n.o. means “not observed”.
Table 1. Computed and experimental absorption data and the calculated values for µ in %. n.a. means “not applicable” and n.o. means “not observed”.
H2O (C2v)Sym.678910111213141516Expt. [59,71,72]
νOH2 ,asB2287539973731377237663766376537653765376537653756
νOH2 ,sA1304137853657369136843683368336823682368236823652
δOH2A1121616881570159015861586158615861586158615861595
μ 21.35.30.80.60.60.60.50.50.50.50.5n.a.
CO2 (Dh)Sym.678910111213141516Expt. [59,73,74]
νCO2,as g + 179125472359239123862386238523852385238523852349
νCO2,s u + 99314081255131913161315131513151315131513151285
δCO2μ487699645654653653653653653653653667
δCO2μ487699645654653653653653653653653667
μ 25.16.92.32.02.02.02.02.02.02.02.0n.a.
CH2O (C2v)Sym.678910111213141516Expt. [59,75,76,77]
νCH2,asB1218430002809283728332833283328322832283228322843
νCH2,sA1225029052801282328182818281828172817281728172782
νCOA1131218551719174517411740174017401740174017401746
δCH2A1115716141502152115181518151815181518151815181500
ρCH2B194513401240125712541254125412541254125412541249
ωCH2B289112501162117611741173117311731173117311731167
μ 23.06.40.80.80.70.70.70.70.70.70.7n.a.
C2H2 (Dh)Sym.678910111213141516Expt. [59,60,78,79,80]
νC2H2,s g + 269934903375339033843380338433843384338433843374
νC2H2,as u + 252334053199323332283227322732273227322732273289
νCC g + 146920461912193719321932193219321932193219321974
δC2H2,su557792733746743743743743743743743730
δC2H2,su557792733746743743743743743743743730
δC2H2,asg443638589602598598598598598598598612
δC2H2,asg443638589602598598598598598598598612
μ 24.55.22.01.71.81.81.81.81.81.81.8n.a.
HCOOH (Cs)Sym.678910111213141516Expt. [59,81,82,83,84]
νOHA'284136823530355935503557355235523552355235523570
νCHA'224631342959298329782977297729772977297729772943
νC=OA'133719101764179017861786178617851785178517851770
δCHA'104314801367138413811381138113811381138113811387
δOHA'98213671265128012771277127712771277127712771387
νC-OA'84811481077108910861086108610861086108610861105
δCHA"78611141029104210401040104010401040104010401033
δOH,oopA"515631588623609608609609609609609638
δOCOA'469659611619617617617617617617617625
μ 22.86.02.11.41.71.71.71.71.71.71.7n.a.
CH4 (Td)Sym.678910111213141516Expt. [59,60,85,86]
νCH3,asF2233332103000303530293029302830283028302830283019
νCH3,asF2229931772990302030153015301430143014301430143019
νCH3,asF2233032092999303330283027302730273027302730273019
νCH4,sA1258829682948294329452944294329432943294329432917
δCH4,asE116616411523154215391539153915391539153915391534
δCH4,asE116616411523154215391539153915391539153915391534
δCH3,sF299313901293130913061306130613061306130613061306
δCH3,sF299213901293130913061306130613061306130613061306
δCH3,sF299313901294130913071307130713061306130613061306
µ 22.35.90.90.40.30.30.30.30.30.30.3n.a.
CH3Cl (C3v)Sym.678910111213141516Expt. [59,87,88,89,90]
νCH3,asE233832033002303530303029302930293029302930293039
νCH3,asE230631922974301930163015301530153015301530143039
νCH3,sA1233930742964298929872986298629852985298529852937
δCH3,asE110615551445146314601460146014601460146014601452
δCH3,asE110615561444146214601459145914591459145914591452
δCH3,sA1103614641357137513721372137113711371137113711355
ρCH3E77310991017103110281028102810281028102810281017
ρCH3E77610991020103310311031103110311031103110311017
νCClA1585802749760758758758758758758758732
µ 23.07.00.91.41.21.21.21.21.21.21.2n.a.
CH3OH (Cs)Sym.678910111213141516Expt. [59,60,91,92,93]
νOHA’281239003646369436873686368636863686368636863681
νCH3,asA’230131982984303530303029302830283028302830283000
νCH3,asA”223030912881291229082908290729072907290729072960
νCH3,sA’236329512849295729472942294129412941294129412844
δCH3,asA’113015871475149314901490149014901490149014901477
δCH3,asA”111615741459147814751475147514751474147414741477
δCH3,sA’110415631447146614631463146314631463146314631455
δOHA’102014581342136213591359135813581358135813581345
ρCH3A”84311781069108310811081108010801080108010801060
νCOA’78910601015102810261026102510251025102510251033
µ 22.96.40.91.31.21.21.21.21.11.11.1n.a.
CH3CN (C3v)Sym.678910111213141516Expt. [59,94,95,96,97,98,99,100]
νCH3,asE228131532940297229672967296729662966296629663009
νCH3,asE232231722976299429902989298929882988298829883009
νCH3,sA1232230592951297029682967296729672967296729672954
νCNA1165623082152217821742174217321732173217321732267
δCH3,asE109915471438145614531453145314531453145314531454
δCH3,asE109915501437145614531453145314531453145314531454
δCH3 ,sA1104714871378139713941394139413931393139313931389
ρCH3E78811211037105210491049104910491049104910491041
ρCH3E78911221038105310501050105010501050105010501041
νCCA1692983923932930930930930930930930920
δCCNE270395361367366366366366366366366361
δCCNE268395359365364364364364364364364361
µ 24.46.41.11.11.01.01.01.01.01.01.0n.a.
C2H4(D2h)Sym.678910111213141516Expt. [59,60,101,102]
νCH2,asB2u239933583122316131553154315431543154315431543106
νCH2,asB1g237933283095313331273126312631263126312631263103
νCH2,sAg229731523047305530523051305130513051305130513026
νCH2,sB3u231832533023306030553054305430543054305430542989
νCCAg127317501628164816441644164416441644164416441623
δCH2B3u109515541439145814551455145514551455145514551444
δCH2Ag105114321345136113591358135813581358135813581342
ρCH2B1g92213181217123512321232123212321232123212321236
τCH2Au79511291046106010581058105810581058105810581023
ωCH2B1u7301036960973971971971971971971971949
ωCH2B2g7121013939951949949949949949949949943
ρCH2B2u619893824836834834834834834834834826
µ 23.47.70.81.51.41.41.31.31.31.31.3n.a.
C2H4O (C2v)Sym.678910111213141516Expt. [59,60,103,104]
νCH2,asB2237233133083312031143114311331133113311331133065
νCH2,asA223613296306731043098309830973097309730973097n.o.
νCH2,sA1237830972939304730643058305730563056305630563018
νCH2,sB1229532142988302430193018301830183018301830183006
δCH2,sA1115316341500151815151515151515151515151515151498
δCH2,asB1112115901472149214891489148914891489148914891472
νCCA1107013281277128812861285128512851285128512851270
ρCH2A28761253115711741171117111711171117111711171n.o.
τCH2B287112441149116611631163116311631163116311631142
ωCH2B185912291135115211491149114911491149114911491151
ωCH2A188212031128114111391138113811381138113811381148
τCH2A27831117103210471044104410441044104410441044n.o.
νOC2,sA1669946877889887887887887887887887877
νOC2,asB1624891822833832831831831831831831872
ρCH2B2615885816829827827827827827827827821
µ 23.36.51.21.41.41.31.31.31.31.31.3n.a.
CH3NH2(Cs)Sym.678910111213141516Expt. [59,105,106,107]
νNH2,asA”262935173318334833433343334333433343334333433411
νNH2,sA’253934313362335233473347334633463346334633463349
νCH3,asA”228331102919294729432942294229422942294229422985
νCH3,asA’228330792907296029572956295629562955295529552963
νCH3,sA’220130102833288929112906290429042904290429042816
δNH2A’122317331596161016081608160816081608160816081642
δCH3,asA”113216011483150315001494150015001500150015001481
δCH3,asA’111915811467148714841483148314831483148314831463
δCH3,sA’108315341421144014371437143614361436143614361450
ρCH3A”9981420131413311328132813281328132813281328n.o.
ρCH3A’89112371152116911671167116711671167116711671144
νCNA’84511091045105910571057105710571057105710571050
τNH2A”7251052967981979979979979979979979n.o.
ωNH2A’637909840852851850850850850850850816
CN axis torsionA”245316345337307313316313313314314304
µ 22.85.92.42.31.61.71.81.71.71.71.7n.a.
C2H6(D3d)Sym.678910111213141516Expt. [59,60,61,62,108,109]
νCH3,asEu230932212999303730313030303030303030303030302985
νCH3,asEg229031852967301130073006300630063006300630062969
νCH3,sA1g247030062950298130042996299529952994299429942954
νCH3,sA2u225031722943298229762975297529752975297529752896
δCH3,asEg112315961474149414911491149114911491149114911468
δCH3,asEu112315931475149614931493149214921492149214921469
δCH3,sA1g105915011393141114081408140814081408140814081388
δCH3,sA2u104214821372139213891389138813881388138813881379
ρCH3Eg90912941197121512121212121212121212121212121190
νCCA1g8181051100310101008100810081008100810081008995
ρCH3Eu620899828842839839839839839839839822
µ 22.37.50.61.71.61.61.61.61.61.61.6n.a.
CH3OCH3(C2v)Sym.678910111213141516Expt. [59,60,110,111]
νCH3,asA1233931422945301630133012301230113011301130112996
νCH3,asB1231331452949297629762976297529752975297529752996
νCH3,asA2224231262906294229362936293629362936293629352952
νCH3,asB2223631172899293429282928292829282928292829282925
νCH3,sA1220629492825286928522856285528552855285528542817
νCH3,sB1220030802861289728912891289128912891289128912817
δCH3,asA1114215971488150615031503150315031503150315031464
δCH3,asB2112215871475149414921491149114911491149114911464
δCH3,asB1111815891475149414911491149114911491149114911464
δCH3,asA2111115791466148514821482148214821482148214821464
δCH3,sA1111415741460147914761475147514751475147514751452
δCH3,sB1108215411426144514421442144214421442144214421452
ρCH3A194413491248126512631262126212621262126212621244
ρCH3B189712751181119711941194119411941194119411941227
ρCH3B289312731180119411951190118911881188118811881179
ρCH3A287112451153116811661166116611661166116611661150
νC2O,asB183912071113112911271127112611261126112611261102
νC2O,sA1780988941950949949949948948948948928
δC2OA1333452423429428428428428428428428418
µ 23.07.21.11.71.61.51.51.51.51.51.5n.a.
C2H5NO2(Cs)Sym.678910111213141516Expt. [112,113]
νOHA’277337373524355735523551355135503550355035503560
νNH2,asA”260935563337337033653364336433643364336433643410
νNH2,sA’26293480334633573352335133513351335133503350n.o.
νCH2,asA”22683101290729372932293229322931293129312931n.o.
νCH2,sA’226530592982295729532952295229522952295229522958
νC=OA’134619211777180217981798179817981798179817981779
δNH2A’122817501602162016171617161716171617161716171630
δCH2A’108515321421144014371437143714371437143714371429
νC−OA’104314821372139113881388138813871387138713871373
τCH2,NH2A”10301466135513741371137113711371137013701370n.o.
ωCH2A’9681382127512931290129012901290129012901290n.o.
δCCN,oopA”8811262116411811178117811781178117811781178n.o.
νCNA’90512451156117211691169116911691169116911691136
δCO2A’83211861097111311101110111011101110111011101101
δCNH2A’7161002922938935935934934934934934907
δC=O,oopA”685990910924922922922922922922922883
νC −CA’658863815823822822822822822822822801
δC=O,ipA’486678629637636636636636636636636619
C-O axis torsionA”369529475485483483483483483483483500
O-C2N axis shearA’365495461467466466466466466466466463
O=C2N axis shearA’194286262266266266266266266266266n.o.
µ 22.87.51.51.81.71.71.71.71.71.71.7n.a.

3.2. Performance at Reduced Grid Densities

As a VSCF evaluation with 16 grid points per mode is computationally highly demanding, a reduction of the grid density is desirable. From the data in Table 1 it becomes evident that for most cases, a reduction to 10 grid points yields identical results as the corresponding calculation with 16 grid points, indicating that the computational demand in a VSCF treatment of i = 1 N V i d i a g ( Q i ) and i < j N V i , j p a i r s ( Q i , Q j ) contributions can be reduced safely by more than 50%. It has to be kept in mind, however, that for every set of reduced densities, the grid is interpolated to the original resolution and the displacement boundaries are not affected by the reduction. Hence, as long as the points on the potential energy surface are chosen properly, the accuracy of the VSCF results is not affected. Noteworthy, the default routine in GAMESS chooses the displacement boundaries in dependence to the underlying normal mode’s wavenumber but the user is nonetheless able to alter these settings if desired.
Table 2. Obtained values for µ in %, their arithmetic means and the corresponding t-values according to Equation (12). n.a. means “not applicable”.
Table 2. Obtained values for µ in %, their arithmetic means and the corresponding t-values according to Equation (12). n.a. means “not applicable”.
Number of Grid Points
Molecule678910111213141516
H2O21.35.30.80.60.60.60.50.50.50.50.5
CO225.16.92.32.02.02.02.02.02.02.02.0
CH2O23.06.40.80.80.70.70.70.70.70.70.7
C2H224.55.22.01.71.81.81.81.81.81.81.8
HCOOH22.86.02.11.41.71.71.71.71.71.71.7
CH422.35.90.90.40.30.30.30.30.30.30.3
CH3Cl23.07.00.91.41.21.21.21.21.21.21.2
CH3OH22.96.40.91.31.21.21.21.21.11.11.1
CH3CN24.46.41.11.11.01.01.01.01.01.01.0
C2H423.47.70.81.51.41.41.31.31.31.31.3
C2H4O23.36.51.21.41.41.31.31.31.31.31.3
CH3NH222.85.92.42.31.61.71.81.81.71.71.7
C2H622.37.50.61.71.61.61.61.61.61.61.6
CH3OCH323.07.21.11.71.61.51.51.51.51.51.5
C2H5NO222.87.51.51.81.71.71.71.71.71.71.7
Mean µ23.16.51.31.41.31.31.31.31.31.31.3
t-Value102.9625.590.092.430.810.781.711.340.662.39n.a.
Reviewing the values for µ in Table 2 permits the conclusion that major errors are starting to appear at grid densities <8 points. The GAMESS manual states that the VSCF code is “thought to give accuracy to 50 cm−1 for the larger fundamentals” when MP2 with a triple-ζ basis set is employed and hence, the boundaries of accuracy are wider than the error inflicted by a reduction of the grid density to 8 points. In order to provide statistical evidence as to what extent a reduction of the grid density is viable, paired t-tests [114] of the obtained values for µ have been conducted with the following equation:
t = | i = 1 χ μ r e f μ c a l c χ | s d × χ
In Equation (12), µref is the reference value at a grid density of 16 points. µcalc is the mean absolute percentage deviation at the examined grid density. χ is the number of investigated molecules and the t-value is calculated with sd as the corrected sample standard deviation of the respective grid density. For a two-tailed problem, critical values for t can be specified according to Table 3.
Table 3. Critical boundaries according to Student’s t-distribution for a two-tailed problem.
Table 3. Critical boundaries according to Student’s t-distribution for a two-tailed problem.
t-Value Difference
t ≤ 2.145Insignificant
2.145 < t ≤ 2.977Probable
2.977 < t ≤ 4.140Significant
4.140 < tHighly significant
Except for the data obtained at grid densities of 6 and 7 points, no significant deviations from the reference data are observed, indicating that even a reduction to 8 grid points is a very viable option. The computational effort involved in a pair-approximated VSCF evaluation (Equation (7)) may thus be reduced by a factor of nearly 4. It has to be stressed, however, that the obtained t-value for r = 8 (i.e., 0.09) may not be confused as to indicate the objectively best agreement with the reference values since the statistical analysis is based on unsigned error values. Moreover, the anharmonic correction seems to be slightly over-estimated at r = 8, yielding more red-shifted absorptions, as can be learned from Table 1. This over-correction is responsible for a slightly improved agreement with experiment for most of the molecules, especially since it is known that MP2 tendentially yields blue-shifted IR absorptions [19]. Overall, it may be concluded that VSCF treatments involving at least 8 grid points reproduce the experimental values within the same boundaries of accuracy as the reference values with a mean absolute percentage deviation from experiment of ~1.3% (Table 2).
In their recent paper, Roy et al., reported that soft and semi-rigid molecules have to be treated at higher grid densities than rigid molecules and they mention ethylene oxide as a particular example where the VSCF equations converge at no less than 12 grid points. Our results contradict this statement since our VSCF calculations of ethylene oxide did indeed converge for grid densities ≥6 points and a density of 8 points yielded the same quality of results as observed for the other molecules. Given the fact that the identical code, the same ab initio method and the same basis set were employed, this outcome is rather puzzling. This disagreement may be ascribed to the use of symmetry during our VSCF calculations and the underlying reference geometries having been minimized to the most stringent criterion available. Overall, none of our VSCF evaluations exhibited convergence problems apart from the few excluded torsional degrees of freedom that have been mentioned earlier.

3.3. Application to Deprotonated Phosphoserine

Phosphorylation of proteins is considered a major signal transduction mechanism, mainly occurring at the OH terminus of the amino acids serine, threonine and tyrosine. As a major constituent of biofluids, phosphoserine was subject of a gas-chromatographic investigation of human urine samples [115]. In autopsied Alzheimer’s disease brain tissue, l-phosphoserine was found in elevated concentrations [116] and only very recently, plasma phosphoserine levels were found to be upregulated in sepsis patients [117]. Deprotonated phosphoserine, abbreviated as [pSer-H], was investigated via infrared multiple photon dissociation (IRMPD) spectroscopy and hence, experimental data of a few major gas phase IR absorptions of the molecule is available [118]. The authors identified the lowest energy conformer of [pSer-H] as the one exhibiting hydrogen bonds between the carboxylic OH and the phosphate O as well as between the phosphate OH and the amino N. This conformer was chosen as a benchmark for the 8 point VSCF analysis due to the fact that it represents rigid and weak structural motifs likewise. Figure 1 shows the minimum geometry obtained at the MP2/cc-pVTZ level of theory which served as a reference state for the VSCF evaluation. Table 4 lists the PT2-VSCF computed absorptions of [pSer-H], the respective intensities and the available experimental data [118]. Due to the cage-like structure, many normal modes couple among each other and hence, especially lower vibrations cannot be ascribed to distinct normal modes. In such cases, the main vibrational contributions are given in Table 4. For illustrative purposes, the calculated IR spectrum of [pSer-H] is shown in Figure 2. The intensities have been computed using harmonic dipole derivatives at the MP2/cc-pVTZ level of theory and band broadening was introduced via a Lorentzian function using a band width at half height of 20 cm−1.
Figure 1. The lowest energy conformer of [pSer-H] exhibits a hydrogen bond between the carboxylic OH and the phosphate O and between the phosphate OH and the amino N.
Figure 1. The lowest energy conformer of [pSer-H] exhibits a hydrogen bond between the carboxylic OH and the phosphate O and between the phosphate OH and the amino N.
Molecules 19 21253 g001
Figure 2. The computed IR spectrum of [pSer-H]. For band assignments, refer to Table 4.
Figure 2. The computed IR spectrum of [pSer-H]. For band assignments, refer to Table 4.
Molecules 19 21253 g002
Table 4. The main absorptions of [pSer-H] and the experimentally obtained [118] values.
Table 4. The main absorptions of [pSer-H] and the experimentally obtained [118] values.
Band NumberPT2-VSCF (cm−1)Intensity (km mol−1) IRMPD (cm−1) [118]Main Contributions
13430225n.o.νOH
2334611n.o.νNH2,as
33139125n.o.νNH2,s
4294521n.o.νCαH
5291547n.o.νCβH2,s
6288835n.o.νCβH2,as
717343591728νCcarb =O
81565391610δNH2
915122441461δCcarbOH,ip
144911419δCβH2
13804n.o.τNH2 and δNCαH
13719n.o.ωCβH2
133768n.o.νCcarb −O and τCβH2
1013154841291νP =O
129221n.o.ωCβH2 and δCcarbCαH and νCcarb −O
11125124n.o.νCαCβ and δCβCαH and νCα −N and τNH2
12306n.o.νCcarbCα and δCcarbCαH and τNH2
12117660n.o.δCcarbOH,oop
13115751n.o.νCβO
141113991108δPOH
15108474n.o.νNCαCcarb,as and νCβO
1610601411052ωNH2
1710171211028νP − O H −bonded
18957113n.o.νNCαCβ,as
1991517n.o.νNCαCβ,s and νNCαCcarb,as and νCcarbCαCβ,as
2083098836νP −OH
82324812δCcarbCαCβ and δCcarbO2,oop
7877n.o.δCcarbO2 and δCcarbCαCβ
21736120738νP −OCβ and δCcarbO2
2269239n.o.δCαNCβCcarb,umbrella and νP − OCbeta
51341n.o.δPOH,oop
Considering the 11 available experimental values, µ was computed for the corresponding VSCF data as 1.38%. Since we did calculate the IR absorptions of [pSer-H] exclusively using 8 grid points per mode, no definitive conclusion can be drawn whether higher grid densities would enable an improved accuracy for this larger molecule. A re-evaluation of the potential energy surface at grid densities up to 16 points would be unmanageably expensive but nonetheless, the data presented fits well into the error boundaries as observed for the 15 molecules discussed earlier.
The band assignment in Table 4 clearly shows that empirical considerations are only applicable to the large fundamentals. For a proper description of low lying vibrations that are determined by a number of torsions, computational approaches are an important option aiding the spectroscopist. In their experimental work, Scuderi et al., have employed computational techniques during their band assignment of [pSer-H] as well, but they resorted to scaled HOA absorption data at the B3LYP/6-311+G(d,p) level of theory.

4. Conclusions

VSCF theory has grown an important field in computational spectroscopy and this is owed to the increasingly available computational resources and to efficient and easily applicable algorithms. With this systematic study it has been demonstrated that a largely reduced grid density in a VSCF evaluation does not only spare considerable resources but also does not significantly affect the resulting absorption data. It was found that the convergence of VSCF equations is not impaired by such a reduced-effort technique even when highly problematic densities of <8 grid points are employed. While further investigations are required with regard to the role of symmetry in a VSCF treatment, it was found that reduced grid densities may be safely applied to a wide set of molecular entities. Application to [pSer-H] showed that a 8 point VSCF calculation yields accurate data for a molecule exhibiting weak interactions. The routine presented herein should, in conjunction with other recent developments [35], prove valuable for larger molecules relevant to life sciences where a conventional VSCF treatment with 16 grid points per mode is not feasible.

Author Contributions

O.M.D.L., B.M.R., G.K.B. and C.W.H. designed research; O.M.D.L. performed research; O.M.D.L. analyzed data; O.M.D.L. wrote the paper. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Scott, A.P.; Radom, L. Harmonic Vibrational Frequencies: An Evaluation of Hartree-Fock, Møller-Plesset, Quadratic Configuration Interaction, Density Functional Theory, and Semiempirical Scale Factors. J. Phys. Chem. 1996, 100, 16502–16513. [Google Scholar]
  2. Sinha, P.; Boesch, S.E.; Gu, C.; Wheeler, R.A.; Wilson, A.K. Harmonic Vibrational Frequencies: Scaling Factors for HF, B3LYP, and MP2 Methods in Combination with Correlation Consistent Basis Sets. J. Phys. Chem. A 2004, 108, 9213–9217. [Google Scholar] [CrossRef]
  3. Bowman, J.M. Self-Consistent Field Energies and Wavefunctions for Coupled Oscillators. J. Chem. Phys. 1978, 68, 608–610. [Google Scholar] [CrossRef]
  4. Cohen, M.; Greita, S.; McEarchran, R. Approximate and Exact Quantum Mechanical Energies and Eigenfunctions for a System of Coupled Oscillators. Chem. Phys. Lett. 1979, 60, 445–450. [Google Scholar] [CrossRef]
  5. Gerber, R.; Ratner, M. A Semiclassical Self-Consistent Field (SC SCF) Approximation for Eigenvalues of Coupled-Vibration Systems. Chem. Phys. Lett. 1979, 68, 195–198. [Google Scholar]
  6. Bowman, J.M. The Self-Consistent-Field Approach to Polyatomic Vibrations. Acc. Chem. Res. 1986, 19, 202–208. [Google Scholar] [CrossRef]
  7. Carter, S.; Culik, S.J.; Bowman, J.M. Vibrational Self-Consistent Field Method for Many-Mode Systems: A New Approach and Application to the Vibrations of CO Adsorbed on Cu(100). J. Chem. Phys. 1997, 107, 10458–10469. [Google Scholar] [CrossRef]
  8. Carney, G.D.; Sprandel, L.L.; Kern, C.W. Variational Approaches to Vibration-Rotation Spectroscopy for Polyatomic Molecules. In Advances in Chemical Physics; Prigogine, I., Rice, S.A., Eds.; John Wiley & Sons,Inc.: Hoboken, NJ, USA, 2007; Volume 37, pp. 305–379. [Google Scholar]
  9. Barone, V. Vibrational Zero-Point Energies and Thermodynamic Functions Beyond the Harmonic Approximation. J. Chem. Phys. 2004, 120, 3059–3065. [Google Scholar] [CrossRef] [PubMed]
  10. Barone, V. Anharmonic Vibrational Properties by a Fully Automated Second-order Perturbative Approach. J. Chem. Phys. 2005, 122. [Google Scholar] [CrossRef]
  11. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuseria, G.E.; Robb, M.A.; Cheeseman, J.R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 09, Revision D.01; Gaussian Inc.: Wallingford, CT, USA, 2013.
  12. Barone, V.; Biczysko, M.; Bloino, J. Fully Anharmonic IR and Raman Spectra of Medium-Size Molecular Systems: Accuracy and Interpretation. Phys. Chem. Chem. Phys. 2014, 16, 1759–1787. [Google Scholar] [CrossRef] [PubMed]
  13. Biczysko, M.; Bloino, J.; Carnimeo, I.; Panek, P.; Barone, V. Fully Ab Initio IR Spectra for Complex Molecular Systems from Perturbative Vibrational Approaches: Glycine as a Test Case. J. Mol. Struct. 2012, 1009, 74–82. [Google Scholar] [CrossRef]
  14. Chaban, G.M.; Jung, J.O.; Gerber, R.B. Ab Initio Calculation of Anharmonic Vibrational States of Polyatomic Systems: Electronic Structure Combined with Vibrational Self-Consistent Field. J. Chem. Phys. 1999, 111, 1823–1829. [Google Scholar] [CrossRef]
  15. Jung, J.O.; Gerber, R.B. Vibrational Wave Functions and Spectroscopy of (H2O)n, n = 2, 3, 4, 5: Vibrational Self-Consistent Field with Correlation Corrections. J. Chem. Phys. 1996, 105, 10332–10348. [Google Scholar] [CrossRef]
  16. Chaban, G.M.; Jung, J.O.; Gerber, R.B. Anharmonic Vibrational Spectroscopy of Hydrogen-Bonded Systems Directly Computed from Ab Initio Potential Surfaces: (H2O)n, n = 2, 3; Cl(H2O)n, n = 1, 2; H+(H2O)n, n = 1, 2; H2O-CH3OH. J. Phys. Chem. A 2000, 104, 2772–2779. [Google Scholar] [CrossRef]
  17. Chaban, G.M.; Jung, J.O.; Gerber, R.B. Anharmonic Vibrational Spectroscopy of Glycine: Testing of Ab Initio and Empirical Potentials. J. Phys. Chem. A 2000, 104, 10035–10044. [Google Scholar] [CrossRef]
  18. Miller, Y.; Chaban, G.; Gerber, R. Theoretical Study of Anharmonic Vibrational Spectra of HNO3, HNO3-H2O, HNO4: Fundamental, Overtone and Combination Excitations. Chem. Phys. 2005, 313, 213–224. [Google Scholar] [CrossRef]
  19. Roy, T.K.; Carrington, T.; Gerber, R.B. Approximate First-Principles Anharmonic Calculations of Polyatomic Spectra Using MP2 and B3LYP Potentials: Comparisons with Experiment. J. Phys. Chem. A 2014, 118, 6730–6739. [Google Scholar] [CrossRef] [PubMed]
  20. Schmidt, M.W.; Baldridge, K.K.; Boatz, J.A.; Elbert, S.T.; Gordon, M.S.; Jensen, J.H.; Koseki, S.; Matsunaga, N.; Nguyen, K.A.; Su, S.; et al. General Atomic and Molecular Electronic Structure System. J. Comput. Chem. 1993, 14, 1347–1363. [Google Scholar] [CrossRef]
  21. Valiev, M.; Bylaska, E.; Govind, N.; Kowalski, K.; Straatsma, T.; Dam, H.V.; Wang, D.; Nieplocha, J.; Apra, E.; Windus, T.; et al. NWChem: A Comprehensive and Scalable Open-Source Solution for Large Scale Molecular Simulations. Comput. Phys. Commun. 2010, 181, 1477–1489. [Google Scholar] [CrossRef]
  22. Werner, H.; Knowles, P.J.; Knizia, G.; Manby, F.R.; Schütz, M.; Celani, P.; Korona, T.; Lindh, R.; Mitrushenkov, A.; Rauhut, G.; et al. MOLPRO, Version 2012.1, A Package of Ab Initio Programs. 2012. Available online: http://www.molpro.net/ (accessed on 15 December 2014).
  23. Møller, C.; Plesset, M.S. Note on an Approximation Treatment for Many-Electron Systems. Phys. Rev. 1934, 46, 618–622. [Google Scholar] [CrossRef]
  24. Fletcher, G.D.; Schmidt, M.W.; Gordon, M.S. Developments in Parallel Electronic Structure Theory. In Advances in Chemical Physics; Prigogine, I., Rice, S.A., Eds.; John Wiley & Sons, Inc.: Hoboken, NJ, USA, 2007; Volume 110, pp. 267–294. [Google Scholar]
  25. Aikens, C.M.; Gordon, M.S. Parallel Unrestricted MP2 Analytic Gradients Using the Distributed Data Interface. J. Phys. Chem. A 2004, 108, 3103–3110. [Google Scholar] [CrossRef]
  26. Sachs, E.S.; Hinze, J.; Sabelli, N.H. Frozen Core Approximation, a Pseudopotential Method Tested on Six States of NaH. J. Chem. Phys. 1975, 62, 3393–3398. [Google Scholar] [CrossRef]
  27. Dunning, T.H. Gaussian Basis Sets for Use in Correlated Molecular Calculations. I. The Atoms Boron Through Neon and Hydrogen. J. Chem. Phys. 1989, 90, 1007–1023. [Google Scholar] [CrossRef]
  28. Bowman, J.M.; Carrington, T.; Meyer, H.D. Variational Quantum Approaches for Computing Vibrational Energies of Polyatomic Molecules. Mol. Phys. 2008, 106, 2145–2182. [Google Scholar] [CrossRef]
  29. Chan, M.; Carrington, T.; Manzhos, S. Anharmonic Vibrations of the Carboxyl Group in Acetic Acid on TiO2: Implications For Adsorption Mode Assignment in Dye-Sensitized Solar Cells. Phys. Chem. Chem. Phys. 2013, 15, 10028–10034. [Google Scholar] [CrossRef] [PubMed]
  30. Sakurai, J.J.; Tuan, S.F. Modern Quantum Mechanics; Addison-Wesley: Boston, MA, USA, 1985; Volume 1. [Google Scholar]
  31. Griffiths, D.J.; Harris, E.G. Introduction to Quantum Mechanics; Prentice Hall: Upper Saddle River, NJ, USA, 1995; Volume 2. [Google Scholar]
  32. Ratner, M.A.; Gerber, R.B. Excited Vibrational States of Polyatomic Molecules: The Semiclassical Self-Consistent Field Approach. J. Phys. Chem. 1986, 90, 20–30. [Google Scholar] [CrossRef]
  33. Watson, J.K. Simplification of the Molecular Vibration-Rotation Hamiltonian. Mol. Phys. 1968, 15, 479–490. [Google Scholar] [CrossRef]
  34. Bowman, J.M.; Carter, S.; Huang, X. MULTIMODE: A Code to Calculate Rovibrational Energies of Polyatomic Molecules. Int. Rev. Phys. Chem. 2003, 22, 533–549. [Google Scholar] [CrossRef]
  35. Roy, T.K.; Gerber, R.B. Vibrational Self-Consistent Field Calculations for Spectroscopy of Biological Molecules: New Algorithmic Developments and Applications. Phys. Chem. Chem. Phys. 2013, 15, 9468–9492. [Google Scholar] [CrossRef] [PubMed]
  36. Gregurick, S.K.; Fredj, E.; Elber, R.; Gerber, R.B. Vibrational Spectroscopy of Peptides and Peptide-Water Complexes: Anharmonic Coupled-Mode Calculations. J. Phys. Chem. B 1997, 101, 8595–8606. [Google Scholar] [CrossRef]
  37. Gregurick, S.K.; Liu, J.H.Y.; Brant, D.A.; Gerber, R.B. Anharmonic Vibrational Self-Consistent Field Calculations as an Approach to Improving Force Fields for Monosaccharides. J. Phys. Chem. B 1999, 103, 3476–3488. [Google Scholar] [CrossRef]
  38. Kowal, A.T. First-principles Calculation of Geometry and Anharmonic Vibrational Spectra of Thioformamide and Thioformamide-d2. J. Chem. Phys. 2006, 124. [Google Scholar] [CrossRef]
  39. Carter, S.; Bowman, J.M.; Handy, N.C. Extensions and Tests of MULTIMODE: A Code to Obtain Accurate Vibration/Rotation Energies of Many-Mode Molecules. Theor. Chem. Acc. 1998, 100, 191–198. [Google Scholar]
  40. Rheinecker, J.; Bowman, J.M. The Calculated Infrared Spectrum of Cl H2O Using a New Full Dimensional Ab Initio Potential Surface and Dipole Moment Surface. J. Chem. Phys. 2006, 125. [Google Scholar] [CrossRef]
  41. Norris, L.S.; Ratner, M.A.; Roitberg, A.E.; Gerber, R.B. Møller-Plesset Perturbation Theory Applied to Vibrational Problems. J. Chem. Phys. 1996, 105, 11261–11267. [Google Scholar] [CrossRef]
  42. Pele, L.; Gerber, R.B. On the Number of Significant Mode-Mode Anharmonic Couplings in Vibrational Calculations: Correlation-Corrected Vibrational Self-Consistent Field Treatment of Di-, Tri-, and Tetrapeptides. J. Chem. Phys. 2008, 128. [Google Scholar] [CrossRef]
  43. Pele, L.; Brauer, B.; Gerber, R. Acceleration of Correlation-Corrected Vibrational Self-Consistent Field Calculation Times for Large Polyatomic Molecules. Theor. Chem. Acc. 2007, 117, 69–72. [Google Scholar] [CrossRef]
  44. Matsunaga, N.; Chaban, G.M.; Gerber, R.B. Degenerate Perturbation Theory Corrections for the Vibrational Self-Consistent Field Approximation: Method and Applications. J. Chem. Phys. 2002, 117, 3541–3547. [Google Scholar]
  45. Bowman, J.M.; Christoffel, K.; Tobin, F. Application of SCF-SI Theory to Vibrational Motion in Polyatomic Molecules. J. Phys. Chem. 1979, 83, 905–912. [Google Scholar] [CrossRef]
  46. Christoffel, K.M.; Bowman, J.M. Investigations of Self-Consistent Field, SCF CI and Virtual Stateconfiguration Interaction Vibrational Energies for a Model Three-Mode System. Chem. Phys. Lett. 1982, 85, 220–224. [Google Scholar] [CrossRef]
  47. Christiansen, O.; Luis, J.M. Beyond Vibrational Self-Consistent-Field Methods: Benchmark Calculations for the Fundamental Vibrations of Ethylene. Int. J. Quantum Chem. 2005, 104, 667–680. [Google Scholar] [CrossRef]
  48. Neff, M.; Rauhut, G. Toward Large Scale Vibrational Configuration Interaction Calculations. J. Chem. Phys. 2009, 131. [Google Scholar] [CrossRef] [PubMed]
  49. Kaulfuss, U.B.; Altenbokum, M. Anharmonic Oscillator as a Test of the Coupled-Cluster Method. Phys. Rev. D 1986, 33, 3658–3664. [Google Scholar] [CrossRef]
  50. Nagalakshmi, V.; Lakshminarayana, V.; Sumithra, G.; Prasad, M.D. Coupled Cluster Description of Anharmonic Molecular Vibrations. Application to O3 and SO2. Chem. Phys. Lett. 1994, 217, 279–282. [Google Scholar] [CrossRef]
  51. Thomsen, B.; Yagi, K.; Christiansen, O. Optimized Coordinates in Vibrational Coupled Cluster Calculations. J. Chem. Phys. 2014, 140. [Google Scholar] [CrossRef]
  52. Culot, F.; Liévin, J. A Multiconfigurational SCF Computational Method for the Resolution of the Vibrational Schrödinger Equation in Polyatomic Molecules. Theor. Chim. Acta 1994, 89, 227–250. [Google Scholar] [CrossRef]
  53. Meyer, H.D. Studying Molecular Quantum Dynamics With the Multiconfiguration Time-Dependent Hartree Method. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2012, 2, 351–374. [Google Scholar] [CrossRef]
  54. Njegic, B.; Gordon, M.S. Exploring the Effect of Anharmonicity of Molecular Vibrations on Thermodynamic Properties. J. Chem. Phys. 2006, 125. [Google Scholar] [CrossRef]
  55. Bounouar, M.; Scheurer, C. The Impact of Approximate VSCF Schemes and Curvilinear Coordinates on the Anharmonic Vibrational Frequencies of Formamide and Thioformamide. Chem. Phys. 2008, 347, 194–207. [Google Scholar]
  56. Strobusch, D.; Scheurer, C. Hierarchical Expansion of the Kinetic Energy Operator in Curvilinear Coordinates for the Vibrational Self-Consistent Field Method. J. Chem. Phys. 2011, 135. [Google Scholar] [CrossRef]
  57. Boatz, J.A.; Gordon, M.S. Decomposition of Normal-Coordinate Vibrational Frequencies. J. Phys. Chem. 1989, 93, 1819–1826. [Google Scholar] [CrossRef]
  58. Yagi, K.; Hirao, K.; Taketsugu, T.; Schmidt, M.W.; Gordon, M.S. Ab Initio Vibrational State Calculations with a Quartic Force Field: Applications to H2CO, C2H4, CH3OH, CH3CCH and C6H6. J. Chem. Phys. 2004, 121, 1383–1389. [Google Scholar] [CrossRef] [PubMed]
  59. Shimanouchi, T. Tables of Molecular Vibrational Frequencies Consolidated Volume I. Natl. Bur. Stand. 1972, 1, 1–164. [Google Scholar]
  60. Herzberg, G. Infrared and Raman Spectra of Polyatomic Molecules; Van Nostrand: New York, NY, USA, 1949. [Google Scholar]
  61. Smith, L.G. The Infra-Red Spectrum of C2H6. J. Chem. Phys. 1949, 17, 139–167. [Google Scholar] [CrossRef]
  62. Weiss, S.; Leroi, G.E. Direct Observation of the Infrared Torsional Spectrum of C2H6, CH3CD3, and C2D6. J. Chem. Phys. 1968, 48, 962–967. [Google Scholar] [CrossRef]
  63. Lutz, O.M.; Rode, B.M.; Bonn, G.K.; Huck, C.W. The Impact of Highly Correlated Potential Energy Surfaces on the Anharmonically Corrected IR Spectrum of Acetonitrile. Spectrochim. Acta Part A 2014, 131, 545–555. [Google Scholar] [CrossRef]
  64. Bartlett, R.J.; Purvis, G.D. Many-Body Perturbation Theory, Coupled-Pair Many-Electron Theory, and the Importance of Quadruple Excitations for the Correlation Problem. Int. J. Quantum Chem. 1978, 14, 561–581. [Google Scholar] [CrossRef]
  65. Pople, J.A.; Krishnan, R.; Schlegel, H.B.; Binkley, J.S. Electron Correlation Theories and Their Application to the Study of Simple Reaction Potential Surfaces. Int. J. Quantum Chem. 1978, 14, 545–560. [Google Scholar] [CrossRef]
  66. Gauss, J.; Stanton, J.F. Perturbative Treatment of Triple Excitations in Coupled-cluster Calculations of Nuclear Magnetic Shielding Constants. J. Chem. Phys. 1996, 104, 2574–2583. [Google Scholar] [CrossRef]
  67. Pople, J.A.; Head-Gordon, M.; Raghavachari, K. Quadratic Configuration Interaction. A General Technique for Determining Electron Correlation Energies. J. Chem. Phys. 1987, 87, 5968–5975. [Google Scholar] [CrossRef]
  68. Piecuch, P.; Włoch, M. Renormalized Coupled-Cluster Methods Exploiting Left Eigenstates of the Similarity-Transformed Hamiltonian. J. Chem. Phys. 2005, 123. [Google Scholar] [CrossRef]
  69. Piecuch, P.; Włoch, M.; Gour, J.R.; Kinal, A. Single-Reference, Size-Extensive, Non-Iterative Coupled-Cluster Approaches to Bond Breaking and Biradicals. Chem. Phys. Lett. 2006, 418, 467–474. [Google Scholar] [CrossRef]
  70. Dennison, D.M. The Infrared Spectra of Polyatomic Molecules Part I. Rev. Mod. Phys. 1931, 3, 280–345. [Google Scholar] [CrossRef]
  71. Barker, E.F.; Sleator, W.W. The Infrared Spectrum of Heavy Water. J. Chem. Phys. 1935, 3, 660–663. [Google Scholar] [CrossRef]
  72. Benedict, W.S.; Gailar, N.; Plyler, E.K. Rotation-Vibration Spectra of Deuterated Water Vapor. J. Chem. Phys. 1956, 24, 1139–1165. [Google Scholar] [CrossRef]
  73. Amat, G.; Pimbert, M. On Fermi Resonance in Carbon Dioxide. J. Mol. Spectrosc. 1965, 16, 278–290. [Google Scholar] [CrossRef]
  74. Suzuki, I. General Anharmonic Force Constants of Carbon Dioxide. J. Mol. Spectrosc. 1968, 25, 479–500. [Google Scholar] [CrossRef]
  75. Davidson, D.W.; Stoicheff, B.P.; Bernstein, H.J. The Infrared and Raman Spectra of Formaldehyde-d1 Vapor. J. Chem. Phys. 1954, 22, 289–294. [Google Scholar] [CrossRef]
  76. Blau, H.H.; Nielsen, H.H. The Infrared Absorption Spectrum of Formaldehyde Vapor. J. Mol. Spectrosc. 1957, 1, 124–132. [Google Scholar] [CrossRef]
  77. Nakagawa, T.; Kashiwagi, H.; Kurihara, H.; Morino, Y. Vibration-Rotation Spectra of Formaldehyde: Band Contour Analysis of the ν2 and ν3 Fundamentals. J. Mol. Spectrosc. 1969, 31, 436–450. [Google Scholar] [CrossRef]
  78. Lafferty, W.J.; Thibault, R.J. High Resolution Infrared Spectra of C 2 12 H2, C12C13H2, and C 2 13 H2. J. Mol. Spectrosc. 1964, 14, 79–96. [Google Scholar] [CrossRef]
  79. Scott, J.; Rao, K. Infrared Absorption Bands of Acetylene: Part I. Analysis of the Bands at 13.7µ. J. Mol. Spectrosc. 1965, 16, 15–23. [Google Scholar] [CrossRef]
  80. Scott, J.; Rao, K.N. Infrared Absorption Bands of Acetylene: Part III. Analysis of the Bands at 3µ. J. Mol. Spectrosc. 1966, 20, 438–460. [Google Scholar] [CrossRef]
  81. Williams, V.Z. Infra-Red Spectra of Monomeric Formic Acid and Its Deuterated Forms. I. High Frequency Region. J. Chem. Phys. 1947, 15, 232–242. [Google Scholar] [CrossRef]
  82. Wilmshurst, J.K. Vibrational Assignment in Monomeric Formic Acid. J. Chem. Phys. 1956, 25, 478–480. [Google Scholar] [CrossRef]
  83. Millikan, R.C.; Pitzer, K.S. Infrared Spectra and Vibrational Assignment of Monomeric Formic Acid. J. Chem. Phys. 1957, 27, 1305–1308. [Google Scholar] [CrossRef]
  84. Miyazawa, T.; Pitzer, K.S. Internal Rotation and Infrared Spectra of Formic Acid Monomer and Normal Coordinate Treatment of Out-Of-Plane Vibrations of Monomer, Dimer, and Polymer. J. Chem. Phys. 1959, 30, 1076–1086. [Google Scholar] [CrossRef]
  85. Allen, H.C.; Plyler, E.K. ν3 Band of Methane. J. Chem. Phys. 1957, 26, 972–973. [Google Scholar] [CrossRef]
  86. Herranz, J.; Morcillo, J.; Gómez, A. The ν2 Infrared Band of CH4 and CD4. J. Mol. Spectrosc. 1966, 19, 266–282. [Google Scholar] [CrossRef]
  87. King, W.T.; Mills, I.M.; Crawford, B. Normal Coordinates in the Methyl Halides. J. Chem. Phys. 1957, 27, 455–457. [Google Scholar] [CrossRef]
  88. Jones, E.; Popplewell, R.; Thompson, H. Vibration—Rotation Bands of Methyl Chloride. Spectrochim. Acta 1966, 22, 669–680. [Google Scholar]
  89. Maki, A.G.; Thibault, R. Analysis of Some Perturbations in the ν4 Band of Methyl Chloride. J. Chem. Phys. 1968, 48, 2163–2167. [Google Scholar] [CrossRef]
  90. Morillon-Chapey, M.; Graner, G. Fine Structure in the ν1 Band of CH3Cl near 2970 cm−1. J. Mol. Spectrosc. 1969, 31, 155–169. [Google Scholar] [CrossRef]
  91. Tanaka, C.; Kuratani, K.; Mizushima, S.I. In-Plane Normal Vibrations of Methanol. Spectrochim. Acta 1957, 9, 265–269. [Google Scholar] [CrossRef]
  92. Van Thiel, M.; Becker, E.D.; Pimentel, G.C. Infrared Studies of Hydrogen Bonding of Methanol by the Matrix Isolation Technique. J. Chem. Phys. 1957, 27, 95–99. [Google Scholar] [CrossRef]
  93. Falk, M.; Whalley, E. Infrared Spectra of Methanol and Deuterated Methanols in Gas, Liquid, and Solid Phases. J. Chem. Phys. 1961, 34, 1554–1568. [Google Scholar] [CrossRef]
  94. Nakagawa, I.; Shimanouchi, T. Rotation-Vibration Spectra and Rotational, Coriolis Coupling, Centrifugal Distortion and Potential Constants of Methyl Cyanide. Spectrochim. Acta 1962, 18, 513–539. [Google Scholar] [CrossRef]
  95. Tolonen, A.; Koivusaari, M.; Paso, R.; Schroderus, J.; Alanko, S.; Anttila, R. The Infrared Spectrum of Methyl Cyanide Between 850 and 1150 cm−1: Analysis of the ν4, ν7, and 3ν18 Bands with Resonances. J. Mol. Spectrosc. 1993, 160, 554–565. [Google Scholar] [CrossRef]
  96. Venkateswarlu, P. The Rotation-Vibration Spectrum of Methyl Cyanide in the Region 1.6 µ–20 µ. J. Chem. Phys. 1951, 19, 293–298. [Google Scholar] [CrossRef]
  97. Venkateswarlu, P. The Infrared Spectrum of Methyl Cyanide. J. Chem. Phys. 1952, 20, 923–923. [Google Scholar] [CrossRef]
  98. Thompson, H.W.; Williams, R.L. The Infra-red Spectra of Methyl Cyanide and Methyl Isocyanide. Trans. Faraday Soc. 1952, 48, 502–513. [Google Scholar]
  99. Parker, F.; Nielsen, A.; Fletcher, W. The Infrared Absorption Spectrum of Methyl Cyanide Vapor. J. Mol. Spectrosc. 1957, 1, 107–123. [Google Scholar] [CrossRef]
  100. Huet, T. The ν1 and ν5 Fundamental Bands of Methyl Cyanide. J. Mol. Struct. 2000, 517, 127–131. [Google Scholar] [CrossRef]
  101. Arnett, R.L.; Crawford, B.L. The Vibrational Frequencies of Ethylene. J. Chem. Phys. 1950, 18, 118–126. [Google Scholar] [CrossRef]
  102. Allen, H.C.; Plyler, E.K. The Structure of Ethylene from Infrared Spectra. J. Am. Chem. Soc. 1958, 80, 2673–2676. [Google Scholar]
  103. Lord, R.C.; Nolin, B. Vibrational Spectra of Ethylene Oxide and Ethylene Oxide-d4. J. Chem. Phys. 1956, 24, 656–658. [Google Scholar]
  104. Potts, W. The Fundamental Vibration Frequencies of Ethylene Oxide and Ethylene Imine. Spectrochim. Acta 1965, 21, 511–527. [Google Scholar] [CrossRef]
  105. Kirby-Smith, J.S.; Bonner, L.G. The Raman Spectra of Gaseous Substances 1. Apparatus and the Spectrum of Methylamine. J. Chem. Phys. 1939, 7, 880–883. [Google Scholar] [CrossRef]
  106. Gray, A.P.; Lord, R.C. Rotation-Vibration Spectra of Methyl Amine and Its Deuterium Derivatives. J. Chem. Phys. 1957, 26, 690–705. [Google Scholar] [CrossRef]
  107. Tsuboi, M.; Hirakawa, A.Y.; Ino, T.; Sasaki, T.; Tamagake, K. Amino Wagging and Inversion in Methylamines. J. Chem. Phys. 1964, 41, 2721–2734. [Google Scholar] [CrossRef]
  108. Cole, A.; Lafferty, W.; Thibault, R. Rotational Fine Structure of the Perpendicular Band, ν7, of Ethane. J. Mol. Spectrosc. 1969, 29, 365–374. [Google Scholar] [CrossRef]
  109. Nakagawa, I.; Shimanouchi, T. Rotation-Vibration Spectra and Rotational, Coriolis Coupling and Potential Constants of Ethane, Ethane-d6 and Ethane-1, 1, 1-d3. J. Mol. Spectrosc. 1971, 39, 255–274. [Google Scholar] [CrossRef]
  110. Taylor, R.C.; Vidale, G.L. Raman Spectrum and Vibrational Assignments of Gaseous Dimethyl Ether. J. Chem. Phys. 1957, 26, 122–123. [Google Scholar] [CrossRef]
  111. Fateley, W.; Miller, F.A. Torsional Frequencies in the Far Infrared-II: Molecules With Two or Three Methyl Rotors. Spectrochim. Acta 1962, 18, 977–993. [Google Scholar]
  112. Stepanian, S.G.; Reva, I.D.; Radchenko, E.D.; Rosado, M.T.S.; Duarte, M.L.T.S.; Fausto, R.; Adamowicz, L. Matrix-Isolation Infrared and Theoretical Studies of the Glycine Conformers. J. Phys. Chem. A 1998, 102, 1041–1054. [Google Scholar] [CrossRef]
  113. Brauer, B.; Chaban, G.M.; Gerber, R.B. Spectroscopically-Tested, Improved, Semi-Empirical Potentials for Biological Molecules: Calculations for Glycine, Alanine and Proline. Phys. Chem. Chem. Phys. 2004, 6, 2543–2556. [Google Scholar] [CrossRef]
  114. Goulden, C. Methods of Statistical Analysis; Wiley: New York, NY, USA, 2007; pp. 50–55. [Google Scholar]
  115. Kataoka, H.; Nakai, K.; Katagiri, Y.; Makita, M. Analysis of Free and Bound O-Phosphoamino Acids in Urine by Gas Chromatography with Flame Photometric Detection. Biomed. Chromatogr. 1993, 7, 184–188. [Google Scholar] [CrossRef] [PubMed]
  116. Mason, R.; Trumbore, M.W.; Pettegrew, J.W. Membrane Interactions of a Phosphomonoester Elevated Early in Alzheimer’s Disease. Neurobiol. Aging 1995, 16, 531–539. [Google Scholar] [CrossRef] [PubMed]
  117. Chiarla, C.; Giovannini, I.; Siegel, J.H. High Phosphoserine in Sepsis: Panel of Clinical and Plasma Amino Acid Correlations. SpringerPlus 2014, 3. [Google Scholar] [CrossRef] [PubMed]
  118. Scuderi, D.; Correia, C.F.; Balaj, O.P.; Ohanessian, G.; Lemaire, J.; Maitre, P. Structural Characterization by IRMPD Spectroscopy and DFT Calculations of Deprotonated Phosphorylated Amino Acids in the Gas Phase. ChemPhysChem 2009, 10, 1630–1641. [Google Scholar] [CrossRef] [PubMed]
  • Sample Availability: Not available.

Share and Cite

MDPI and ACS Style

Lutz, O.M.D.; Rode, B.M.; Bonn, G.K.; Huck, C.W. Largely Reduced Grid Densities in a Vibrational Self-Consistent Field Treatment Do Not Significantly Impact the ResultingWavenumbers. Molecules 2014, 19, 21253-21275. https://doi.org/10.3390/molecules191221253

AMA Style

Lutz OMD, Rode BM, Bonn GK, Huck CW. Largely Reduced Grid Densities in a Vibrational Self-Consistent Field Treatment Do Not Significantly Impact the ResultingWavenumbers. Molecules. 2014; 19(12):21253-21275. https://doi.org/10.3390/molecules191221253

Chicago/Turabian Style

Lutz, Oliver M. D., Bernd M. Rode, Günther K. Bonn, and Christian W. Huck. 2014. "Largely Reduced Grid Densities in a Vibrational Self-Consistent Field Treatment Do Not Significantly Impact the ResultingWavenumbers" Molecules 19, no. 12: 21253-21275. https://doi.org/10.3390/molecules191221253

APA Style

Lutz, O. M. D., Rode, B. M., Bonn, G. K., & Huck, C. W. (2014). Largely Reduced Grid Densities in a Vibrational Self-Consistent Field Treatment Do Not Significantly Impact the ResultingWavenumbers. Molecules, 19(12), 21253-21275. https://doi.org/10.3390/molecules191221253

Article Metrics

Back to TopTop