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

Discrete Oligomers and Polymers of Chloride Monohydrate Can Form in Encapsulated Environments: Structures and Infrared Spectra of [Cl4 (H2 O)4 ]4- and {[Cl(H2 O)]- }∞

A discrete tetrachloride tetrahydrate cluster, [Cl₄ (H₂ O)₄ ]^4- , was obtained with a partially-fluorinated triaminocyclopropenium cation, [C₃ (N(CH₂ CF₃ )₂ )(NEt₂ )(NPr₂ )]⁺ . The cluster consists of a [Cl₂ (H₂ O)₂ ]^2- square with each Cl^- coordinated by another H₂ O bridged to another Cl^- . A 1D polymer of chloride monohydrate, {[Cl(H₂ O)]^- }_∞ , was obtained with [C₃ (N(CH₂ CF₃ )₂ )₂ (NBuMe)]⁺ . The tetrameric and polymeric structures were found to be computationally-unstable in the gas phase which indicates that an encapsulated environment is essential for their isolation. DFT and DFTB calculations were carried out on gas-phase [Cl₄ (H₂ O)₄ ]^4- to assist the infrared assignments. Anharmonically-corrected B3LYP transition frequencies were in close agreement with experiment, but DFTB models were only appropriate for qualitative interpretation. Solid-state DFTB calculations allowed the vibrational modes to be assigned. The results found are consistent with "discrete" chloride hydrates.

My Trajectory in Molecular Reaction Dynamics and Spectroscopy

This is the story of a career in theoretical chemistry during a time of dramatic changes in the field due to phenomenal growth in the availability of computational power. It is likewise the story of the highly gifted graduate students and postdoctoral fellows that I was fortunate to mentor throughout my career. It includes reminiscences of the great mentors that I had and of the exciting collaborations with both experimentalists and theorists on which I built much of my research. This is an account of the developments of exciting scientific disciplines in which I was involved: vibrational spectroscopy, molecular reaction mechanisms and dynamics, e.g., in atmospheric chemistry, and the prediction of new, exotic molecules, in particular noble gas molecules. From my very first project to my current work, my career in science has brought me the excitement and fascination of research. What a wonderful pursuit!

Dual Basis Approach for Ab Initio Anharmonic Calculations of Vibrational Spectroscopy: Application to Microsolvated Biomolecules

A dual electronic basis set approach is introduced for more efficient but accurate calculations of the anharmonic vibrational spectra in the framework of the vibrational self-consistent field (VSCF) theory. In this approach, an accurate basis set is used to compute the vibrational spectra at the harmonic level. The results are used to scale the potential surface from a more modest but much more efficient basis set. The scaling is such that at the harmonic level the new, scaled potential agrees with one of the accurate basis sets. The approach is tested in the application of the microsolvated, protected amino acid Ac-Phe-OMe, using the scaled anharmonic hybrid potential in the VSCF and VSCF-PT2 algorithms. The hybrid potential method yields results that are in good accord with the experiment and very close to those obtained in calculations with the high-level, very costly potential from the large basis set. At the same time, the hybrid potential calculations are considerably less expensive. The results of the hybrid calculations are much more accurate than those computed from the potential surface corresponding to the modest basis set. The results are very encouraging for using the hybrid potential method for inexpensive yet sufficiently accurate anharmonic calculations for the spectra of large biomolecules.

A Series of Discrete Dichloride Dihydrates: Characterisation and Symmetry Effects

A series of three discrete dichloride dihydrates [Cl₂ (H₂ O)₂ ]^2- have been isolated with different triaminocyclopropenium (TAC) cations and with different crystallographic symmetries. The cluster exhibits D_2h symmetry with the tris(dimethylamino)cyclopropenium cation [C₃ (NMe₂ )₃ ]⁺ , C_2h symmetry with the fluorinated cation [C₃ (N(CH₂ CF₃ )₂ )(NBu₂ )₂ ]⁺ (containing two 2,2,2-trifluoroethyl substituents) and C_2v symmetry with the more fluorinated [C₃ (N(CH₂ CF₃ )₂ )₂ (NBu₂ )]⁺ cation. The effect of symmetry on the infrared spectra of the dichloride ion-pair clusters, as well as deuterated analogues, has been investigated. The D_2h - and C_2h -symmetric clusters each exhibit two stretching bands in the infrared at 3427 and 3368 cm^-1 for D_2h symmetry and 3444 and 3392 cm^-1 for C_2h symmetry, whereas the C_2v -symmetric cluster exhibits three bands at 3475, 3426 and 3373 cm^-1 . Computational studies were carried out on a [Cl₂ (H₂ O)₂ ]^2- cluster with C_2v symmetry to aid the infrared band assignments.

One water to bind a chloride-chloride ion pair: isolation of discrete [Cl2(H2O)]2- in the solid state

A discrete dichloride ion pair in the form of a monohydrate, [Cl₂(H₂O)]^2-, was isolated using the triaminocyclopropenium cation [C₃(NHex₂)(N(CH₂CF₃)₂)₂]⁺. Although this ion pair is calculated to be unstable in the gas phase, the ionic lattice and weak CH-Cl hydrogen bonds assist the stabilization of the cluster. The D₂O and HDO isotopomers were also prepared and characterized.

A Discrete Dichloride Tetrahydrate Trapped by a Cyclopropenium Cation: Structure and Spectroscopic Properties

A discrete dichloride tetrahydrate cluster, [Cl₂ (H₂ O)₄ ]^2- , was obtained as a salt of the bis(diphenylamino)diethylamino cyclopropenium cation [C₃ (NPh₂ )₂ (NEt₂ )]⁺ and characterized by single-crystal X-ray diffraction and infrared spectroscopy. This chloride-chloride ion-pair cluster consists of a [Cl₂ (H₂ O)₂ ]^2- square with opposite edges bridged by water molecules to give a chair-like structure of the non-hydrogen atoms. The solid-state structure is essentially the same as the calculated gas-phase structure. Infrared spectra were also collected on the deuterium analogue [Cl₂ (D₂ O)₄ ]^2- . Computational studies were carried out on gas-phase [Cl₂ (H₂ O)₄ ]^2- to confirm the infrared band assignments in the solid state. The structure and infrared spectrum are consistent with the discrete nature of the cluster.

Gas phase dynamics, conformational transitions and spectroscopy of charged saccharides: the oxocarbenium ion, protonated anhydrogalactose and protonated methyl galactopyranoside

Protonated intermediates are postulated to be involved in the rate determining step of many sugar reactions. This paper presents a study of protonated sugar species, isolated in the gas phase, using a combination of infrared multiple photon dissociation (IRMPD) spectroscopy, classical ab initio molecular dynamics (AIMD) and quantum mechanical vibrational self-consistent field (VSCF) calculations. It provides a likely identification of the reactive intermediate oxocarbenium ion structure in a d-galactosyl system as well as the saccharide pyrolysis product anhydrogalactose (that suggests oxocarbenium ion stabilization), along with the spectrum of the protonated parent species: methyl d-galactopyranoside-H+. Its vibrational fingerprint indicates intramolecular proton sharing. Classical AIMD simulations for galactosyl oxocarbenium ions, conducted in the temperature range ∼300-350 K (using B3LYP potentials on-the-fly) reveal efficient transitions on the picosecond timescale. Multiple conformers are likely to exist under the experimental conditions and along with static VSCF calculations, they have facilitated the identification of the individual structural motifs of the galactosyl oxocarbenium ion and protonated anhydrogalactose ion conformers that contribute to the observed experimental spectra. These results demonstrate the power of experimental IRMPD spectroscopy combined with dynamics simulations and with computational spectroscopy at the anharmonic level to unravel conformer structures of protonated saccharides, and to provide information on their lifetimes.

Temperature Dependence of Intramolecular Vibrational Bands in Small Water Clusters

Cyclic water clusters are pivotal for understanding atmospheric reactions as well as liquid water, yet the temperature (T) dependence of their dynamics and spectroscopy is poorly studied. The development of highly accurate water potentials, such as MB-pol, partly rectifies this. It remains to account for the quantum nuclear effects (NQE), because quantum nuclear dynamics become increasingly inaccurate at low temperatures. From a practical point of view, we find that NQE can be accounted for simply by subtracting a constant from the frequencies obtained from the velocity autocorrelation functions (VACF) of classical nuclear dynamics, resulting in unprecedented agreement with experiment, mostly within 5 cm^-1. We have performed classical simulations of (H₂O)_n clusters (n = 2-5) from 20 K and up to their melting temperature, calculating both all-atom and partial VACF, thus generating the temperature dependence of the vibrational frequencies (IR and Raman bands). Focusing on the hydrogen-bonded (HBed) OH stretch and HOH bend, we find opposing T dependencies. The HBed OH modes blue shift linearly with T, attributed to ring expansion rather than any specific conformational change. The lowest-frequency Raman concerted mode is predicted to show the largest such shift. In contrast, the HOH bend undergoes a red-shift, with the highest frequency concerted band undergoing the largest red-shift. These results can be explained by a coupled-oscillator model for n hydrogen atoms on a ring, constrained to move either tangentially (stretch) or perpendicularly (bend) to the ring. With increasing temperature and weakening of HBs, the intrinsic force constant increases (stretch) or remains constant (bend), while the nearest-neighbor coupling constant decreases, and this results in the interesting behavior revealed herein. T-dependent Raman studies are required for testing some of these predictions.

Anharmonic effect of the unimolecular dissociation of Glycerol to Glycidol

The unimolecular dissociation rate constants of the dehydration of Glycerol to Glycidol were calculated at the MP2/6–311G(d,p) level using the Rice–Ramsperger–Kassel–Marcus (RRKM) theory. The anharmonic effect of the reactions was examined by comparing the rate constants at temperatures (700–3000[Formula: see text]K) of the canonical case and total energies (25654–53089[Formula: see text]cm[Formula: see text]) of the microcanonical system. The calculations showed that high temperatures are required for the reaction to proceed. As the temperatures and total energies increased, the rate of reactions increased. However, the growth rate of the unimolecular dissociation rate constants was high and slower both in the canonical and microcanonical systems. Comparative analysis showed that the anharmonic effect was most significant for the reaction [Formula: see text] and least significant for the reaction [Formula: see text]. The anharmonic effect became more significant as the temperatures and total energies increased. Compared with the microcanonical situation, the anharmonic effect of the canonical system was more pronounced.

MP4 Study of the Anharmonic Coupling of the Shared Proton Stretching Vibration of the Protonated Water Dimer in Equilibrium and Transition States

The structure and harmonic and anharmonic IR spectra of the protonated water dimer (PWD) were calculated in C₁, C₂, and C_s symmetry at the MP4/acc-pVTZ level of theory. We found that structure and IR spectra are practically identical in C₂ and C₁ symmetry, demonstrating that an equilibrium C₁ configuration of the PWD is not realized. Anharmonic coupling of the shared proton stretching vibration with all other modes in the PWD in C₂ and C_s symmetry was the focus of this investigation. For this purpose, 28 two-dimensional potential energy surfaces (2D PES) were built at the MP4/acc-pVTZ level of theory and the corresponding vibrational Schrödinger equations were solved using the DVR method. Differences in the coupling of the investigated mode with other modes in the C₂ and C_s configurations, along with some factors that determine the red- or blue-shift of the stretching vibration frequency, were analyzed. We obtained a rather reasonable value of the stretching frequency of the bridging proton (1058.4 cm^-1) unperturbed by Fermi resonance. The Fermi resonance between the fundamental vibration ν₇ and the combined vibration ν₂ + ν₆ of the same symmetry was analyzed through anharmonic second-order perturbation theory calculations, as well as by 3D PES constructed using Q₂, Q₆, and Q₇ as normal coordinates. A significant (up to 50%) transfer of intensity from the fundamental vibration to the combined one was found. We have estimated the frequency of the bridging proton stretching vibration in the C_s configuration of the PWD based on calculations of the intrinsic anharmonicity and anharmonic double modes interactions at the MP4/acc-pVTZ level of theory (1261 cm^-1).

Effects of Global and Local Anharmonicities on the Thermodynamic Properties of Sulfuric Acid Monohydrate

We use state-of-the-art electronic structure calculation methods and large basis sets to obtain reliable values for the thermodynamic properties of sulfuric acid monohydrate and study the effects of vibrational anharmonicity on these properties. We distinguish between two forms of vibrational anharmonicity: local anharmonicity, which refers to the anharmonicity of the vibrational modes of a given cluster conformer, and global anharmonicity, which originates from accounting for the presence of different conformers in the first place. In our most accurate approach, we solve the nuclear Schrödinger equation variationally for the intermolecular large-amplitude motions, thus quantum-mechanically accounting for the presence of higher-energy conformers for both reactants and products, while using the standard vibrational perturbational approach for the other vibrational modes. This results in a value of -11.0 kJ/mol for the reaction Gibbs free energy at 298.15 K. When standard vibrational perturbational approaches are employed, the effects of local anharmonicity depend heavily on the choice of the electronic structure calculation basis set. In fact, better results can often be achieved by combining a simple harmonic treatment for the vibrational partition function with a statistical mechanical accounting of global anharmonicity. Thus, we recommend that future studies that intend to include anharmonicity start by accounting for the presence of higher-energy conformers and only then consider whether local anharmonicity calculations are feasible and necessary.

Radiative cooling of H3O+ and its deuterated isotopologues

In conjunction with ab initio potential energy and dipole moment surfaces for the electronic ground state, we have made a theoretical study of the radiative lifetimes for the hydronium ion H3O+ and its deuterated isotopologues. We compute the ro-vibrational energy levels and their associated wavefunctions together with Einstein coefficients for electric dipole transitions. A detailed analysis of the stability of the ro-vibrational states has been carried out and the longest-living states of the hydronium ions have been identified. We report estimated radiative lifetimes and cooling functions for temperatures <200 K. A number of long-living meta-stable states are identified, capable of population trapping.

BNg3F3: the first three noble gas atoms inserted into mono-centric neutral compounds - a theoretical study

Following the study of HXeOXeH and HXeCCXeH, in which two Xe atoms were inserted into H2O and C2H2 theoretically and experimentally, the structures and stability of BNg3F3 (Ng = Ar, Kr and Xe), in which three Ng atoms are inserted into BF3, have been explored theoretically using DFT and ab initio calculations. It is shown that BNg3F3 (Ng = Ar, Kr and Xe) with D3h symmetry are local minima with short B-Ng bond lengths of 1.966, 2.027 and 2.214 Å at the CCSD(T)/aug-cc-pVTZ/LJ18 level, which are close to their covalent limits. Note that although BNg3F3 (Ng = Kr and Xe) are energetically higher than the dissociation products 3Ng + BF3, they are still kinetically stable as metastable species with protecting barriers of 13.38 and 17.99 kcal mol(-1) for BKr3F3 and BXe3F3. Moreover, BKr3F3, the tri-Kr-inserted compound, even has comparable kinetic stability to HXeOXeH and HXeOXeF. In addition, upon the formation of BNg3F3, there is a large amount of charge transferred from B to Ng of at least 0.619 e. The calculated Wiberg Bond Indices (WBI) suggest that B-Ng bonds are naturally singly bonded; the large vibrational frequencies of B-Ng and Ng-F stretching modes and the negative Laplacian electron density of B-Ng bonds confirm further that BNg3F3 are stiff molecules with covalent B-Ng bonds. It should be noted that three Ng atoms inserted into mono-centric neutral molecules have not been reported previously. We hope that the present theoretical study may provide important evidence for the experimental synthesis of BNg3F3.

Calculation of anharmonic effects for the unimolecular dissociation of CH3OOH and its deuterated species CD3OOD using the Rice–Ramsperger–Kassel–Marcus theory

Anharmonic and harmonic rate constants for the unimolecular dissociation of CH3OOH and CD3OOD were calculated using the Rice–Ramsperger–Kassel–Marcus (RRKM) theory at the MP2/6-311++G(3df,3pd) level of theory. The anharmonic effect of the reactions was investigated. Comparison of results for the decompositions of CH3OOH and CD3OOD shows that the direct bond dissociation channel, CH3(D3) O + OH (D), is the most dominant reaction. The anharmonic effect plays an important role in the unimolecular dissociation of both CH3OOH and CD3OOD. For channels CH3(D3) O + OH (D) and CH3(D3) + H (D) O2, the anharmonic effect of the unimolecular dissociation of CD3OOD is more pronounced than that of the unimolecular dissociation of CH3OOH. For channel H2(D2) CO + H2(D2) O, the anharmonic effect of the unimolecular dissociation of CH3OOH is more pronounced than that of the unimolecular dissociation of CD3OOD. The isotope effect is more distinct in the anharmonic oscillator model.

Representing the potential-energy surface of protonated water clusters by high-dimensional neural network potentials

We report a reactive neural network potential for protonated water clusters that accurately represents the density-functional theory potential-energy surface.

Largely reduced grid densities in a vibrational self-consistent field treatment do not significantly impact the resultingwavenumbers

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.

Approximate first-principles anharmonic calculations of polyatomic spectra using MP2 and B3LYP potentials: comparisons with experiment

Anharmonic vibrational spectroscopy calculations using MP2 and B3LYP computed potential surfaces are carried out for a series of molecules, and frequencies and intensities are compared with those from experiment. The vibrational self-consistent field with second-order perturbation correction (VSCF-PT2) is used in computing the spectra. The test calculations have been performed for the molecules HNO3, C2H4, C2H4O, H2SO4, CH3COOH, glycine, and alanine. Both MP2 and B3LYP give results in good accord with experimental frequencies, though, on the whole, MP2 gives very slightly better agreement. A statistical analysis of deviations in frequencies from experiment is carried out that gives interesting insights. The most probable percentage deviation from experimental frequencies is about -2% (to the red of the experiment) for B3LYP and +2% (to the blue of the experiment) for MP2. There is a higher probability for relatively large percentage deviations when B3LYP is used. The calculated intensities are also found to be in good accord with experiment, but the percentage deviations are much larger than those for frequencies. The results show that both MP2 and B3LYP potentials, used in VSCF-PT2 calculations, account well for anharmonic effects in the spectroscopy of molecules of the types considered.

Hybrid MP2/MP4 potential surfaces in VSCF calculations of IR spectra: applications for organic molecules

This study introduces an improved hybrid MP2/MP4 ab initio potential for vibrational spectroscopy calculations which is very accurate, yet without high computational demands. The method uses harmonic vibrational calculations with the MP4(SDQ) potential to construct an improved MP2 potential by coordinate scaling. This improved MP2 potential is used for the anharmonic VSCF calculation. The method was tested spectroscopically for four molecules: butane, acetone, ethylene and glycine. Very good agreement with experiment was found. For most of the systems, the more accurate harmonic treatment considerably improved the MP2 anharmonic results.

O-H anharmonic vibrational motions in Cl(-)···(CH3OH)(1-2) ionic clusters. Combined IRPD experiments and AIMD simulations

The structures of Cl(-)-(Methanol)1,2 clusters have been unraveled combining Infrared Predissociation (IR-PD) experiments and DFT-based molecular dynamics simulations (DFT-MD) at 100 K. The dynamical IR spectra extracted from DFT-MD provide the initial 600 cm(-1) large anharmonic red-shift of the O-H stretch from uncomplexed methanol (3682 cm(-1)) to Cl(-)-(Methanol)1 complex (3085 cm(-1)) as observed in the IR-PD experiment, as well as the subtle supplementary blue- and red-shifts of the O-H stretch in Cl(-)-(Methanol)2 depending on the structure. The anharmonic vibrational calculations remarkably provide the 100 cm(-1) O-H blue-shift when the two methanol molecules are simultaneously organized in the anion first hydration shell (conformer 2A), while they provide the 240 cm(-1) O-H red-shift when the second methanol is in the second hydration shell of Cl(-) (conformer 2B). RRKM calculations have also shown that 2A/2B conformers interconvert on a nanosecond time-scale at the estimated 100 K temperature of the clusters formed by evaporative cooling of argon prior to the IR-PD process.

Anharmonic Rice-Ramsperger-Kassel-Marcus (RRKM) and product branching ratio calculations for the partially deuterated protonated water dimers: dissociation and isomerization

Partially deuterated protonated water dimers, H2O·H(+)·D2O, H2O·D(+)·HDO, and HDO·H(+)·HDO, as important intermediates of isotopic labeled reaction of H3O(+) + D2O, undergo direct dissociation and indirect dissociation, i.e., isomerization before the dissociation. With Rice-Ramsperger-Kassel-Marcus theory and ab initio calculations, we have computed their dissociation and isomerization rate constants separately under the harmonic and anharmonic oscillator models. On the basis of the dissociation and isomerization rate constants, branching ratios of two primary products, [HD2O(+)]∕[H2DO(+)], are predicted under various kinetics models with the harmonic or anharmonic approximation included. The feasible kinetics model accounting for experimental results is shown to include anharmonic effect in describing dissociation, while adopting harmonic approximation for isomerization. Thus, the anharmonic effect is found to play important roles affecting the dissociation reaction, while isomerization rates are shown to be insensitive to whether the anharmonic or harmonic oscillator model is being applied.

Raman spectra of long chain hydrocarbons: anharmonic calculations, experiment and implications for imaging of biomembranes

First-principles anharmonic vibrational calculations are carried out for the Raman spectrum of the C-H stretching bands in dodecane, and for the C-D bands in the deuterated molecule. The calculations use the Vibrational Self-Consistent Field (VSCF) algorithm. The results are compared with liquid-state experiments, after smoothing the isolated-molecule sharp-line computed spectra. Very good agreement between the computed and experimental results is found for the two systems. The combined theoretical and experimental results provide insights into the spectrum, elucidating the roles of symmetric and asymmetric CH(3) and CH(2) hydrogenic stretches. This is expected to be very useful for the interpretation of spectra of long-chain hydrocarbons. The results show that anharmonic effects on the spectrum are large. On the other hand, vibrational degeneracy effects seem to be rather modest at the resolution of the experiments. The degeneracy effects may have more pronounced manifestations in higher-resolution experiments. The results show that first-principles anharmonic vibrational calculations for hydrocarbons are feasible, in good agreement with experiment, opening the way for applications to many similar systems. The results may be useful for the analysis of CARS imaging of lipids, for which dodecane is a representative molecule. It is suggested that first-principles vibrational calculations may be useful also for CARS imaging of other systems.

Effect of magnesium cation on the interfacial properties of aqueous salt solutions

Sodium chloride solutions have been used extensively as a model of seawater in both theoretical and experimental studies of the chemistry of sea salt aerosol. Many groups have found that chloride anions are present at the air-solution interface. This observation has been important for the development of a mechanism for the heterogeneous production of molecular chlorine from chloride in sea salt aerosol. However, while sodium chloride is a major constituent of seawater, it is by no means the only salt present. Seawater contains one Mg(2+) for every eight Na(+). Mg(2+) is naturally occurring in ocean waters from mineral deposits in the Earth's crust and biological sources. Mg(2+) forms a hexahydrate structure, rather than contact ion pairs with chloride anion, and this impacts the ordering of water in solution. In this study, we use molecular dynamics simulations, ab initio calculations, and vibrational sum frequency generation (SFG) spectroscopy to explore the effect of the Mg(2+) cation and its tightly bound solvation shell on the surface propensity of chloride, ion-ion interactions, and water structure of the air-solution interface of concentrated chloride salt solutions. In addition, we provide molecular level details that may be relevant to the heterogeneous reactions of chloride in deliquesced sea salt aerosols. In particular, we show that the presence of the divalent Mg(2+) cation does not modify the surface propensity of chloride compared to Na(+) and hence, its availability to interfacial reaction, although some differences in the behavior of chloride may occur due to specific ion interactions. In this work, we also discuss the SFG free OH band at the surface of salt solutions and conclude that it is often not straightforward to interpret.

Approximate solution of the mode-mode coupling integral: Application to cytosine and its deuterated derivative

Ab initio Hartree-Fock (HF), density functional theory (DFT) and second-order Møller-Plesset (MP2) methods were used to perform harmonic and anharmonic calculations for the biomolecule cytosine and its deuterated derivative. The anharmonic vibrational spectra were computed using the vibrational self-consistent field (VSCF) and correlation-corrected vibrational self-consistent field (CC-VSCF) methods. Calculated anharmonic frequencies have been compared with the argon matrix spectra reported in literature. The results were analyzed with focus on the properties of anharmonic couplings between pair of modes. A simple and easy to use formula for calculation of mode-mode coupling magnitudes has been derived. The key element in present approach is the approximation that only interactions between pairs of normal modes have been taken into account, while interactions of triples or more are neglected. FTIR and Raman spectra of solid state cytosine have been recorded in the regions 400-4000 cm(-1) and 60-4000 cm(-1), respectively. Vibrational analysis and assignments are based on calculated potential energy distribution (PED) values.