How does trimethylamine N-oxide counteract the denaturing activity of urea?

Trimethylamine N-oxide, TMAO, stabilizes globular proteins and is able to counteract the denaturing activity of urea. The mechanism of this counteraction has remained elusive up to now. A rationalization is proposed grounded on the same theoretical model used to clarify the origin of cold denaturation, and the denaturing activity of GdmCl versus the stabilizing one of Gdm(2)SO(4) [G. Graziano, Phys. Chem. Chem. Phys., 2010, 12, 14245-14252; G. Graziano, Phys. Chem. Chem. Phys., 2011, 13, 12008-12014]. The fundamental quantities are: (a) the difference in the solvent-excluded volume on passing from the N-state to the D-state, calculated in water and in aqueous osmolyte solution; (b) the difference in energetic attractions of the N-state and the D-state with the surrounding solvent molecules, calculated in water and in aqueous osmolyte solution. In aqueous 8 M urea + 4 M TMAO solution, the first quantity is so large and positive to counteract the second one that is large and negative due to preferential binding of urea molecules to the protein surface. This happens because aqueous 8 M urea + 4 M TMAO solution has a volume packing density markedly larger than that of water, rendering the cavity creation process much more costly. The volume packing density increase reflects the strength of the attractions of water molecules with both urea and TMAO molecules. This mechanism readily explains why TMAO counteraction is operative even though urea molecules are preferentially located on the protein surface.     

Contrasting the denaturing effect of guanidinium chloride with the stabilizing effect of guanidinium sulfate

Guanidinium chloride, GdmCl, is a strong denaturing agent of globular proteins, whereas guanidinium sulfate, Gdm(2)SO(4), is a stabilizing agent of globular proteins. The stabilizing activity of Gdm(2)SO(4) is unexpected because the denaturant capability of GdmCl is due to direct interactions of Gdm(+) ions with protein surface groups. It is shown that the statistical thermodynamic approach devised to explain the molecular origin of cold denaturation [G. Graziano, Phys. Chem. Chem. Phys., 2010, 12, 14245-14252] can provide a rationalization of the different behaviour of GdmCl and Gdm(2)SO(4) towards globular proteins. The fundamental quantity is the reversible work to create in the aqueous solution a cavity suitable to host the D-state and a cavity suitable to host the N-state. In aqueous GdmCl solutions, this contribution is not large enough to overwhelm the conformational entropy gain upon unfolding and the direct attractions between Gdm(+) ions and protein surface groups; in aqueous Gdm(2)SO(4) solutions, it is so large that it overwhelms the two destabilizing contributions. Sulfate ions, due to their high charge density, interact strongly with water molecules producing a number density increase, that, in turn, renders the cavity creation process very costly, reversing the denaturing power of Gdm(+) ions and stabilizing the N-state of globular proteins.     

Solvent effects on the local structure of p-nitroaniline in supercritical water and supercritical alcohols

Raman spectra of p-nitroaniline in supercritical water and supercritical alcohols were measured, and the effects of solvents on the NO 2 and NH 2 stretching modes were investigated. The intensity and frequency of the NO 2 stretching mode significantly changed as a function of the solvent density and temperature. The frequency of the NO 2 stretching mode correlated with the absorption peak energy of the S 1<--S 0 transition. On the other hand, the vibrational frequency of the NH 2 stretching mode did not correlate with the absorption peak shift, although it had a large frequency shift as a function of the density. The correlation between the NO 2 frequency and absorption peak energy suggested that the solvent effects of supercritical water and supercritical alcohols were similar to those for nonpolar solvents. The density functional calculation using the polarizable continuum model and p-nitroaniline-water clusters qualitatively reproduced the density dependence of the NO 2 stretching mode as well as the solvent polarity dependence. Detailed vibrational analysis revealed that the coupling between the NO 2 and C-NH 2 vibrational motions at the harmonic level has an important effect on the intensity and frequency shift of the NO 2 stretching mode. The frequency shift of the NH 2 stretching mode correlated with the degree of hydrogen bonding between the solvent molecules estimated from NMR measurements [Hoffmann M. M.; Conradi, M. S. J. Phys. Chem. B. 1998, 102, 263]. The existence of intermolecular hydrogen bonding around the NH 2 group was demonstrated even at low-density conditions.     

Effect of nitrogen adsorption on the mid-infrared spectrum of water clusters

Experimental Fourier-transform infrared spectra and DFT calculated infrared spectra are compared to investigate the effect of adsorbed nitrogen on the OH-stretch band complex of water clusters. Using a collisional cooling experiment, pure as well as partially and completely N(2)-covered water clusters consisting of 20-200 water molecules have been generated in thermal equilibrium in the aerosol phase within the temperature range of 5-80 K. Computational IR-spectra simulations have been performed for discrete pure and N(2)-covered water clusters including 10, 15, 20, and 30 water molecules. The adsorbed N(2) molecules especially affect the three-coordinated water molecules at the cluster surface which could be observed as a blue shift of the companion O-H band at 2900 cm(-1) and a red shift of the dangling O-H band at 3700 cm(-1) by about 20 cm(-1) in both cases. The most striking effect of the N(2) adsorbate is an intensity increase of the dangling O-H band by a factor of 3-5. Furthermore, the onset temperature of nitrogen adsorption at the water cluster surface was experimentally found to be roughly 30 K for cluster sizes of about 100 water molecules. Experimental and computational results are in good agreement. The presented results are based on and support the work of V. Buch, J. P. Devlin, and co-workers (e.g., J. Phys. Chem. B, 1997; J. Phys. Chem. A, 2003; Int. Rev. Phys. Chem., 2004).     

Water-carbon dioxide mixtures at high temperatures and pressures: local order in the water rich phase investigated by vibrational spectroscopy

Raman scattering combined with near- and mid-infrared absorption spectroscopies was used to investigate the evolution of the local order in the water rich phase of water-CO(2) mixtures under isobaric heating (T=40-360 degrees C, P=250 bars). The quantitative analysis of the spectra shows that tetramers and larger oligomers are the main constituents of water at moderate temperatures below 80 degrees C. As the temperature increases, the dimer and trimer concentrations considerably increase at the expense of larger oligomers. Finally, water dimers are predominant at the highest temperature investigated close to the temperature of total miscibility of the mixture (T=366 degrees C, P=250 bars). This result is consistent with our previous investigation [R. Oparin T. Tassaing, Y. Danten, and M. Besnard, J. Chem. Phys. 120, 10691 (2004)] on water dissolved in the CO(2) rich phase where we found that close to the temperature of total miscibility water also exists mainly under dimeric form. The current study combined with that mentioned above provides a model investigation of the evolution of the state of aggregation of water molecules in binary mixture involving a hydrophobic solvent in a wide range of temperature.     

Molecular dynamics simulations of surface-initiated melting of nitromethane

The melting of nitromethane initiated at solid-vacuum interfaces has been investigated using molecular dynamics nvt simulations with a realistic force field [D. C. Sorescu et al., J. Phys. Chem. B 104, 8406 (2000)]. The calculated melting point (251+/-5 K) is in good agreement with experiment (244.73 K) and values obtained previously (approximately 255.5 and 266.5+/-8 K) using other simulation methods [P. M. Agrawal et al., J. Chem. Phys. 119, 9617 (2003)]. Analyses of the molecular orientations and diffusion during the simulations as functions of the distance from the exposed surfaces show that the melting is a direct crystal-to-liquid transition, in which the molecules first gain rotational freedom, then mobility. There is a slight dependence of the melting temperature on the exposed crystallographic face.     

The importance of polarizability in the modeling of solubility: quantifying the effect of solute polarizability on the solubility of small nonpolar solutes in popular models of water

In recent work by Paschek [J. Chem. Phys. 120, 6674 (2004)] and others [see H. Docherty et al., J. Chem. Phys. 125, 074510 (2006) for a review] it has been suggested that, when coupled to a simple Lennard-Jones model for various small nonpolar solute molecules, the most common models of water (e.g., SPC/E and TIP4P) fail to reproduce quantitatively the solubility of small nonpolar solute molecules in water due in part to failing to account for polarization of the solute molecule. Given the importance of such systems as test-case prototype models of the solubility of proteins and biomolecules, in this work, we investigate the impact of using a polarizable solute model with the SPC/E, TIP3P, TIP4P, TIP4P-Ew, and TIP4P/2005 rigid water models. Specifically we consider Ne, Ar, Kr, Xe, and methane as solutes. In all cases we observe that the use of a polarizable solute improves agreement between experiment and simulations, with the best agreement seen for the largest solutes, Kr, CH(4), and Xe and the modern reparametrizations of the TIP4P model, i.e., the TIP4P-Ew and TIP4P/2005 models.     

Heat capacity effects associated with the hydrophobic hydration and interaction of simple solutes: a detailed structural and energetical analysis based on molecular dynamics simulations

We examine the SPCE [H. J. C. Berendsen et al., J. Chem. Phys. 91, 6269 (1987)] and TIP5P [M. W. Mahoney and W. L. Jorgensen, J. Chem. Phys 112, 8910 (2000)] water models using a temperature series of molecular-dynamics simulations in order to study heat-capacity effects associated with the hydrophobic hydration and interaction of xenon particles. The temperature interval between 275 and 375 K along the 0.1-MPa isobar is studied. For all investigated models and state points we calculate the excess chemical potential for xenon employing the Widom particle insertion technique. The solvation enthalpy and excess heat capacity is obtained from the temperature dependence of the chemical potentials and, alternatively, directly by Ewald summation, as well as a reaction field based method. All three methods provide consistent results. In addition, the reaction field technique allows a separation of the solvation enthalpy into solute/solvent and solvent/solvent parts. We find that the solvent/solvent contribution to the excess heat capacity is dominating, being about one order of magnitude larger than the solute/solvent part. This observation is attributed to the enlarged heat capacity of the water molecules in the hydration shell. A detailed spatial analysis of the heat capacity of the water molecules around a pair of xenon particles at different separations reveals that even more enhanced heat capacity of the water located in the bisector plane between two adjacent xenon atoms is responsible for the maximum of the heat capacity found for the desolvation barrier distance, recently reported by Shimizu and Chan [J. Am. Chem. Soc. 123, 2083 (2001)]. The about 60% enlarged heat capacity of water in the concave part of the joint xenon-xenon hydration shell is the result of a counterplay of strengthened hydrogen bonds and an enhanced breaking of hydrogen bonds with increasing temperature. Differences between the two models with respect to the heat capacity in the xenon-xenon contact state are attributed to the different water model bulk heat capacities, and to the different spatial extension of the structure effect introduced by the hydrophobic particles. Similarities between the different states of water in the joint xenon-xenon hydration shell and the properties of stretched water are discussed.     

Monte Carlo simulations of critical cluster sizes and nucleation rates of water

We have calculated the critical cluster sizes and homogeneous nucleation rates of water at temperatures and vapor densities corresponding to experiments by Wolk and Strey [J. Phys. Chem B 105, 11683 (2001)]. The calculations have been done with an expanded version of a Monte Carlo method originally developed by Vehkamaki and Ford [J. Chem. Phys. 112, 4193 (2000)]. Their method calculates the statistical growth and decay probabilities of molecular clusters. We have derived a connection between these probabilities and kinetic condensation and evaporation rates, and introduce a new way for the calculation of the work of formation of clusters. Three different interaction potential models of water have been used in the simulations. These include the unpolarizable SPC/E [J. Phys. Chem. 91, 6269 (1987)] and TIP4P [J. Chem. Phys. 79, 926 (1983)] models and a polarizable model by Guillot and Guissani [J. Chem. Phys. 114, 6720 (2001)]. We show that TIP4P produces critical cluster sizes and a temperature and vapor density dependence for the nucleation rate that agree well with the experimental data, although the magnitude of nucleation rate is constantly overestimated by a factor of 2 x 10(4). Guissani and Guillot's model is somewhat less successful, but both the TIP4P and Guillot and Guissani models are able to reproduce a much better experimental temperature dependency of the nucleation rate than the classical nucleation theory. Using SPC/E results in dramatically too small critical clusters and high nucleation rates. The water models give different average binding energies for clusters. We show that stronger binding between cluster molecules suppresses the decay probability of a cluster, while the growth probability is not affected. This explains the differences in results from different water models.     

Parallel algorithm for the computation of the Hartree-Fock exchange matrix: gas phase and periodic parallel ONX

In this paper we present an efficient parallelization of the ONX algorithm for linear computation of the Hartree-Fock exchange matrix [J. Chem. Phys. 106, 9708 (1997)]. The method used is based on the equal time (ET) partitioning recently introduced [J. Chem. Phys. 118, 9128 (2003)] and [J. Chem. Phys. 121, 6608 (2004)]. ET exploits the slow variation of the density matrix between self-consistent-field iterations to achieve load balance. The method is presented and some benchmark calculations are discussed for gas phase and periodic systems with up to 128 processors. The current parallel ONX code is able to deliver up to 77% overall efficiency for a cluster of 50 water molecules on 128 processors (2.56 processors per heavy atom) and up to 87% for a box of 64 water molecules (two processors per heavy atom) with periodic boundary conditions.     

Effects of nonpolar solutes on the thermodynamic response functions of aqueous mixtures

We investigate the effect of adding nonpolar solutes at atmospheric pressure on water's temperature of maximum density, isothermal compressibility, and isobaric heat capacity, using a statistical mechanical model of water solutions [H. S. Ashbaugh, T. M. Truskett, and P. G. Debenedetti, J. Chem. Phys. 116, 2907 (2002)]. We find that the temperature of maximum density increases with solute hydrophobicity, as characterized by its size, and decreases with its van der Waals attractive parameter a, in agreement with experiment. We predict similar trends for the addition of solutes on the isothermal compressibility and isobaric heat capacity: solute hydrophobicity causes an upward shift in water's anomalies, whereas dispersive interactions as measured by the solute's van der Waals attractive parameter shift the anomalies to lower temperatures. The locus along which the competing contributions of solute size sigma and interaction strength a to the shift in water's response functions balance each other obeys the scaling relationship sigma6 approximately a.     

Electronic structure and spectrum of UO2 2+ and UO2Cl4 2-

A theoretical study is presented of the electronic spectra of the UO(2) (2+) and UO(2)Cl(4) (2-) ions, based on multiconfigurational perturbation theory (CASSCF/CASPT2), combined with a recently developed method to treat spin-orbit coupling [P.-A. Malmqvist et al., Chem. Phys. Lett. 357, 230 (2002); B. O. Roos and P.-A. Malmqvist, Phys. Chem. Chem. Phys. 6, 2919 (2004)]. The results are compared to the experimental spectroscopic data obtained for uranyl ions in Cs(2)UO(2)Cl(4) crystals from Denning [Struct. Bonding (Berlin) 79, 215 (1992)] and to previous theoretical calculations performed using a combined configuration-interaction spin-orbit treatment [Z. Zhang and R. M. Pitzer, J. Phys. Chem. A 103, 6880 (1999); S. Matsika and R. M. Pitzer, J. Phys. Chem. A. 105, 637 (2001)]. As opposed to the latter results, the calculations performed in this work point to a significant effect of the weakly bound equatorial chlorine ligands on the excitation energies.     

A rationale for the contrasting activity (towards globular proteins) of tert-butyl alcohol and trimethylamine N-oxide

tert-Butyl alcohol, TBA, and trimethylamine N-oxide, TMAO, even though they are isosteric molecules, have contrasting activity towards globular proteins: the former destabilizes the native state, whereas the latter stabilizes it. TBA addition to water causes a density decrease and it tends to form self-aggregates; TMAO addition to water causes a density increase and it does not show any tendency to form self-aggregates. By inserting such experimental information in the framework of the statistical thermodynamic model devised to rationalize the conformational stability of globular proteins [G. Graziano, Phys. Chem. Chem. Phys., 2010, 12, 14245-14252], it emerges that: (a) TBA is a destabilizing agent because its addition to water causes a significant decrease in the magnitude of the contribution due to the solvent-excluded volume decrease associated with protein folding; (b) TMAO is a stabilizing agent because its addition to water causes a significant increase in the magnitude of the contribution due to the solvent-excluded volume decrease associated with protein folding.     

Thermodynamic anomalies in a lattice model of water: solvation properties

We investigate a lattice-fluid model of water, defined on a three-dimensional body-centered-cubic lattice. Model molecules possess a tetrahedral symmetry, with four equivalent bonding arms. The model is similar to the one proposed by Roberts and Debenedetti [J. Chem. Phys. 105, 658 (1996)], simplified by removing distinction between "donors" and "acceptors." We focus on the solvation properties, mainly as far as an ideally inert (hydrophobic) solute is concerned. As in our previous analysis, devoted to neat water [J. Chem. Phys. 121, 11856 (2004)], we make use of a generalized first-order approximation on a tetrahedral cluster. We show that the model exhibits quite a coherent picture of water thermodynamics, reproducing qualitatively several anomalous properties observed both in pure water and in solutions of hydrophobic solutes. As far as supercooled liquid water is concerned, the model is consistent with the second critical-point scenario.     

Molecular dynamics simulations of melting and the glass transition of nitromethane

Molecular dynamics simulations have been used to investigate the thermodynamic melting point of the crystalline nitromethane, the melting mechanism of superheated crystalline nitromethane, and the physical properties of crystalline and glassy nitromethane. The maximum superheating and glass transition temperatures of nitromethane are calculated to be 316 and 160 K, respectively, for heating and cooling rates of 8.9 x 10(9) Ks. Using the hysteresis method [Luo et al., J. Chem. Phys. 120, 11640 (2004)] and by taking the glass transition temperature as the supercooling temperature, we calculate a value of 251.1 K for the thermodynamic melting point, which is in excellent agreement with the two-phase result [Agrawal et al., J. Chem. Phys. 119, 9617 (2003)] of 255.5 K and measured value of 244.73 K. In the melting process, the nitromethane molecules begin to rotate about their lattice positions in the crystal, followed by translational freedom of the molecules. A nucleation mechanism for the melting is illustrated by the distribution of the local translational order parameter. The critical values of the Lindemann index for the C and N atoms immediately prior to melting (the Lindemann criterion) are found to be around 0.155 at 1 atm. The intramolecular motions and molecular structure of nitromethane undergo no abrupt changes upon melting, indicating that the intramolecular degrees of freedom have little effect on the melting. The thermal expansion coefficient and bulk modulus are predicted to be about two or three times larger in crystalline nitromethane than in glassy nitromethane. The vibrational density of states is almost identical in both phases.     

Polarizable and flexible model for ethanol

A polarizable, flexible model for ethanol is obtained based on an extensive series of B3LYP/6-311++G(d,p) calculations and molecular dynamics simulations. The ethanol model includes electric-field dependence in both the atomic charges and the intramolecular degrees of freedom. Field-dependent intramolecular potentials have been attempted only once previously, for OH and HH stretches in water [P. Cicu et al., J. Chem. Phys. 112, 8267 (2000)]. The torsional potential involving the hydrogen-bonding hydrogen in ethanol is found to be particularly field sensitive. The methodology for developing field-dependent potentials can be readily generalized to other molecules and is discussed in detail. Molecular dynamics simulations of bulk ethanol are performed and the results are assessed based on comparisons with the self-diffusion coefficient [N. Karger et al., J. Chem. Phys. 93, 3437 (1990)], dielectric constant [J. T. Kindt and C. A. Schmuttenmaer, J. Phys. Chem. 100, 10373 (1996)], enthalpy of vaporization [R. C. Wilhoit and B. J. Zwolinski, J. Phys. Chem. Ref. Data, Suppl. 2, 2 (1973)], and experimental interatomic distributions [C. J. Benmore and Y. L. Loh, J. Chem. Phys. 112, 5877 (2000)]. The simultaneous variation of the atomic charges and the intramolecular potentials requires modified equations of motion and a multiple time step algorithm has been implemented to solve these equations. The article concludes with a discussion of the bulk structure and properties with an emphasis on the hydrogen bonding network.     

Molecular nitrogen-N2 properties: the intermolecular potential and the equation of state

Quantum mechanical (QM) high precision calculations were used to determine N(2)-N(2) intermolecular interaction potential. Using QM numerical data the anisotropic potential energy surface was obtained for all orientations of the pair of the nitrogen molecules in the rotation invariant form. The new N(2)-N(2) potential is in reasonably good agreement with the scaled potential obtained by van der Avoird et al. using the results of Hartree-Fock calculations [J. Chem. Phys. 84, 1629 (1986)]. The molecular dynamics (MD) of the N(2) molecules has been used to determine nitrogen equation of state. The classical motion of N(2) molecules was integrated in rigid rotor approximation, i.e., it accounted only translational and rotational degrees of freedom. Fincham [Mol. Simul. 11, 79 (1993)] algorithm was shown to be superior in terms of precision and energy stability to other algorithms, including Singer [Mol. Phys. 33, 1757 (1977)], fifth order predictor-corrector, or Runge-Kutta, and was therefore used in the MD modeling of the nitrogen pressure [S. Krukowski and P. Strak, J. Chem. Phys. 124, 134501 (2006)]. Nitrogen equation of state at pressures up to 30 GPa (300 kbars) and temperatures from the room temperature to 2000 K was obtained using MD simulation results. Results of MD simulations are in very good agreement (the error below 1%) with the experimental data on nitrogen equation of state at pressures below 1 GPa (10 kbars) for temperatures below 1800 K [R. T. Jacobsen et al., J. Phys. Chem. Ref. Data 15, 735 (1986)]. For higher temperatures, the deviation is slightly larger, about 2.5% which still is a very good agreement. The slightly larger difference may be attributed to the vibrational motion not accounted explicitly by rigid rotor approximation, which may be especially important at high temperatures. These results allow to obtain reliable equation of state of nitrogen for pressures up to 30 GPa (300 kbars), i.e., close to molecular nitrogen stability limit, determined by Nellis et al. [Phys. Rev. Lett. 53, 1661 (1984)].     

Existence of a size-dependent diffusivity maximum for uncharged solutes in water and its implications

Recent studies suggest that there exists a size-dependent diffusivity maximum in binary mixtures interacting via Lennard-Jones potential when the size of one of the two components is varied (Ghorai, P. K.; Yashonath, S. J. Phys. Chem., 2005, 109, 5824). We discuss in the present paper the importance of the existence of a size-dependent maximum for an uncharged solute in liquid or amorphous solid water and its relation to the ionic conductivity maximum in water. We report molecular dynamics investigations into the size dependence of the self-diffusivity, D, of the uncharged solutes in water at low temperatures (30 K) with immobile as well as mobile water. We find that a maximum in self-diffusivity exists as a function of the size of solute diffusing within water at low temperatures but not at high temperatures. This is due to the relatively weak interactions between the solute and the water compared to the kinetic energy at room temperature. Previously, we have shown that a similar maximum exists for guests sorbed in zeolites and is known as the levitation effect (LE). Thus, it appears that the existence of a size-dependent maximum is universal and extends from zeolites to simple liquids to solvents of polyatomic species. We examine the implications of this for the size-dependent maximum in ionic conductivity in polar solvents known for over a hundred years. These results support the view that the size-dependent maximum seen for ions in water has its origin in the LE (see Ghorai, P. Kr.; Yashonath, S.; Lynden-Bell, R. M. J. Phys. Chem. 2005, 109, 8120).     

Modeling collective behavior of molecules in nanoscale direct deposition processes

We present a theoretical model describing the collective behavior of molecules in nanoscale direct deposition processes such as dip-pen nanolithography. We show that strong intermolecular interactions combined with nonuniform substrate-molecule interactions can produce various shapes of molecular patterns including fractal-like structures. Computer simulations reveal circular and starlike patterns at low and intermediate densities of preferentially attractive surface sites, respectively. At large density of such surface sites, the molecules form a two-dimensional invasion percolation cluster. Previous experimental results showing anisotropic patterns of various chemical and biological molecules correspond to the starlike regime [P. Manandhar et al., Phys. Rev. Lett. 90, 115505 (2003); J.-H. Lim and C. A. Mirkin, Adv. Mater. (Weinheim, Ger.) 14, 1474 (2002); D. L. Wilson et al., Proc. Natl. Acad. Sci. U.S.A. 98, 13660 (2001); M. Su et al., Appl. Phys. Lett. 84, 4200 (2004); R. McKendry et al., Nano Lett. 2, 713 (2002); H. Zhou et al., Appl. Surf. Sci. 236, 18 (2004); G. Agarwal et al., J. Am. Chem. Soc. 125, 580 (2003)].     

Complexes of type C6H7(+)·L (L = N2 and CO2) studied by explicitly correlated coupled cluster theory

Complexes of the benzenium ion (C(6)H(7)(+)) with N(2) or CO(2) have been studied by explicitly correlated coupled cluster theory at the CCSD(T)-F12x (x = a, b) level [T. B. Adler et al., J. Chem. Phys. 127, 221106 (2007)] and the double-hybrid density functional B2PLYP-D [T. Schwabe and S. Grimme, Phys. Chem. Chem. Phys. 9, 3397 (2007)]. Improved harmonic vibrational wavenumbers for C(6)H(7)(+) have been obtained by CCSD(T?)-F12a calculations with the VTZ-F12 basis set. Combining them with previous B2PLYP-D anharmonic contributions we arrive at anharmonic wavenumbers which are in excellent agreement with recent experimental data from p-H(2) matrix isolation IR spectroscopy [M. Bahou et al., J. Chem. Phys. 136, 154304 (2012)]. The energetically most favourable conformer of C(6)H(7)(+)·N(2) shows a π-bonded structure similar to C(6)H(7)(+)·Rg (Rg = Ne, Ar) [P. Botschwina and R. Oswald, J. Phys. Chem. A 115, 13664 (2011)] with D(e) ≈ 870 cm(-1). For C(6)H(7)(+)·CO(2), a slightly lower energy is calculated for a conformer with the CO(2) ligand lying in the ring-plane of the C(6)H(7)(+) moiety (D(e) ≈ 1508 cm(-1)). It may be discriminated from other conformers through a strong band predicted at 1218 cm(-1), red-shifted by 21 cm(-1) from the corresponding band of free C(6)H(7)(+).