MULTIMODE quantum calculations of intramolecular vibrational energies of the water dimer and trimer using ab initio-based potential energy surfaces

We report vibrational configuration interaction calculations of the monomer fundamentals of (H2O)(2), (D2O)(2), (H2O)(3), and (D2O)(3) using the code MULTIMODE and full dimensional ab initio-based global potential energies surfaces (PESs). For the dimer the HBB PES [Huang , J. Chem. Phys 128, 034312 (2008)] is used and for the trimer a new PES, reported here, is used. The salient properties of the new trimer PES are presented and compared to previous single-point calculations and the vibrational energies are compared with experiments. (C) 2008 American Institute of Physics.

Potential energy surface and moleculardynamics of MbNO: existence of an unsuspected FeON minimum

Ligands such as CO, O2, or NO are involved in the biological function of myoglobin. Here we investigate the energetics and dynamics of NO interacting with the Fe(II) heme group in native myoglobin using ab initio and molecular dynamics simulations. At the global minimum of the ab initio potential energy surface (PES), the binding energy of 23.4 kcal/mol and the Fe-NO structure compare well with the experimental results. Interestingly, the PES is found to exhibit two minima: There exists a metastable, linear Fe-O-N minimum in addition to the known, bent Fe-N-O global minimum conformation. Moreover, the T-shaped configuration is found to be a saddle point, in contrast to the corresponding minimum for NO interacting with Fe(III). To use the ab initio results for finite temperature molecular dynamics simulations, an analytical function was fitted to represent the Fe-NO interaction. The simulations show that the secondary minimum is dynamically stable up to 250 K and has a lifetime of several hundred picoseconds at 300 K. The difference in the topology of the heme-NO PES from that assumed previously (one deep, single Fe-NO minimum) suggests that it is important to use the full PES for a quantitative understanding of this system. Why the metastable state has not been observed in the many spectroscopic studies of myoglobin interacting with NO is discussed, and possible approaches to finding it are outlined.

Hydrogen-bond assisted stabilization of the less favored conformation of a tridentate Schiff base ligand in dinuclear nickel(II) complex: an experimental and theoretical study

The Schiff base ligand, HL (2-[1-(3-methylamino-propylimino)-ethyl]-phenol), the 1:1 condensation product of 2-hydroxy acetophenone and N-methyl-1,3-diaminopropane, has been synthesized and characterized by X-ray crystallography as the perchlorate salt [H2L]ClO4 (1). The structure consists of discrete [H2L](+) cations and perchlorate anions. Two dinuclear Ni-II complexes, [Ni2L2(NO2)(2)] (2), [Ni2L2(NO3)(2)] (3) have been synthesized using this ligand and characterized by single crystal X-ray analyses. Complexes 2 and 3 are centrosymmetric dimers in which the Ni-II ions are in distorted fac- and mer-octahedral environments, respectively, bridged by two mu(2)-phenolate ions of deprotonated ligand, L. The plane of the phenyl rings and the Ni2O2 basal plane are nearly coplanar in 2 but almost perpendicular in 3. We have studied and explained this different behavior using high level DFT calculations (RI-BP86/def2-TZVP level of theory). The conformation observed in 3, which is energetically less favorable, is stabilized via intermolecular non-covalent interactions. Under the excitation of ultraviolet light, characteristic fluorescence of compound 1 was observed; by comparison fluorescence intensity decreases in case of compound 3 and completely quenched in compound 2.

Mono-aqua-bridged dinuclear complexes of Cu(II) containing NNO donor Schiff base ligand: Hydrogen-bond-mediated exchange coupling

Two new mono-aqua-bridged dinuclear Cu(II) complexes of tridentate NNO Schiff bases, [Cu-2(mu-H2O)L-2(1)(H2O)(2)](BF4)(2)center dot 2H(2)O (1) and [Cu-2(mu-H2O)L-2(2)(H2O)(2)](BF4)(2)center dot 2H(2)O (2) where HL1 = 2-[1-(2-dimethylamino-ethylimino)-ethyl]-phenol and HL2 =2-[(2-dimethylamino-ethylimino)-methyl]-phenol were synthesized. Both the complexes were characterized by single-crystal X-ray diffraction analyses and variable-temperature magnetic measurements. For both the complexes each Cu(II) ion is in a square-pyramidal environment being bonded to three atoms from the tridentate NNO Schiff base and a terminal H2O molecule in the equatorial plane; a second H2O ligand acts as a bridge between the two Cu(II) centres through the axial positions. Hydrogen bonds between the terminal H2O ligand and the Schiff base of the adjacent centre complete the intra-dimer linkages. Variable-temperature (4-300 K) magnetic susceptibility measurement shows the presence of significant antiferromagnetic coupling for both the complexes (J = -12.2 and -12.5 cm(-1), respectively, for 1 and 2), mediated mainly through the intra-dimer H-bonds.

Properties of singular vectors using convective available potential energy as final time norm

We study the feasibility of using the singular vector technique to create initial condition perturbations for short-range ensemble prediction systems (SREPS) focussing on predictability of severe local storms and in particular deep convection. For this a new final time semi-norm based on the convective available potential energy (CAPE) is introduced. We compare singular vectors using the CAPE-norm with SVs using the more common total energy (TE) norm for a 2-week summer period in 2007, which includes a case of mesoscale extreme rainfall in the south west of Finland. The CAPE singular vectors perturb the CAPE field by increasing the specific humidity and temperature of the parcel and increase the lapse rate above the parcel in the lower troposphere consistent with physical considerations. The CAPE-SVs are situated in the lower troposphere. This in contrast to TE-SVs with short optimization times which predominantly remain in the high troposphere. By examining the time evolution of the CAPE singular values we observe that the convective event in the south west of Finland is clearly associated with high CAPE singular values.

Available potential energy and exergy in stratified fluids

Lorenz���s theory of available p otential energy (APE) remains the main framework for studying the atmospheric and oceanic energy cycles. Because the APE generation rate is the volume integral of a thermodynamic e���ciency times the local diabatic heating/cooling rate, APE theory is often regarded as an extension of the theory of heat engines. Available energetics in classical thermodynamics, however, usually relies on the concept of exergy, and is usually measured relative to a reference state maximising entropy at constant energy, whereas APE���s reference state minimises p otential energy at constant entropy. This review seeks to shed light on the two concepts; it covers local formulations of available energetics, alternative views of the dynamics/thermodynamics coupling, APE theory and the second law, APE production/dissipation, extensions to binary ���uids, mean/eddy decomp ositions, APE in incompressible ���uids, APE and irreversible turbulent mixing, and the role of mechanical forcing on APE production.

Application of hydrogen-bond propensity calculations to an indomethacin���nicotinamide (1:1) co-crystal

The crystal structure of an indomethacin���nicotinamide (1 : 1) cocrystal produced by milling has been determined from laboratory powder X-ray diffraction (PXRD) data. The hydrogen bonding motifs observed in the structure represent one of the most probable of all the possible combinations of donors and acceptors in the constituent molecules.

A unified theory of available potential energy

Traditional derivations of available potential energy, in a variety of contexts, involve combining some form of mass conservation together with energy conservation. This raises the questions of why such constructions are required in the first place, and whether there is some general method of deriving the available potential energy for an arbitrary fluid system. By appealing to the underlying Hamiltonian structure of geophysical fluid dynamics, it becomes clear why energy conservation is not enough, and why other conservation laws such as mass conservation need to be incorporated in order to construct an invariant, known as the pseudoenergy, that is a positive���definite functional of disturbance quantities. The available potential energy is just the non���kinetic part of the pseudoenergy, the construction of which follows a well defined algorithm. Two notable features of the available potential energy defined thereby are first, that it is a locally defined quantity, and second, that it is inherently definable at finite amplitude (though one may of course always take the small���amplitude limit if this is appropriate). The general theory is made concrete by systematic derivations of available potential energy in a number of different contexts. All the well known expressions are recovered, and some new expressions are obtained. The possibility of generalizing the concept of available potential energy to dynamically stable basic flows (as opposed to statically stable basic states) is also discussed.

Available potential energy density for a multicomponent Boussinesq fluid with arbitrary nonlinear equation of state

In this paper, the concept of available potential energy (APE) density is extended to a multicomponent Boussinesq ���uid with a nonlinear equation of state. As shown by previous studies, the APE density is naturally interpreted as the work against buoyancy forces that a parcel needs to perform to move from a notional reference position at which its buoyancy vanishes to its actual position; because buoyancy can be de���ned relative to an arbitrary reference state, so can APE density. The concept of APE density is therefore best viewed as de���ning a class of locally de���ned energy quantities, each tied to a di���erent reference state, rather than as a single energy variable. An important result, for which a new proof is given, is that the volume integrated APE density always exceeds Lorenz���s globally de���ned APE, except when the reference state coincides with Lorenz���s adiabatically re-arranged reference state of minimum potential energy. A parcel reference position is systematically de���ned as a level of neutral buoyancy (LNB): depending on the nature of the ���uid and on how the reference state is de���ned, a parcel may have one, none, or multiple LNB within the ���uid. Multiple LNB are only possible for a multicomponent ���uid whose density depends on pressure. When no LNB exists within the ���uid, a parcel reference position is assigned at the minimum or maximum geopotential height. The class of APE densities thus de���ned admits local and global balance equations, which all exhibit a conversion with kinetic energy, a production term by boundary buoyancy ���uxes, and a dissipation term by internal di���usive e���ects. Di���erent reference states alter the partition between APE production and dissipation, but neither a���ect the net conversion between kinetic energy and APE, nor the di���erence between APE production and dissipation. We argue that the possibility of constructing APE-like budgets based on reference states other than Lorenz���s reference state is more important than has been previously assumed, and we illustrate the feasibility of doing so in the context of an idealised and realistic oceanic example, using as reference states one with constant density and another one de���ned as the horizontal mean density ���eld; in the latter case, the resulting APE density is found to be a reasonable approximation of the APE density constructed from Lorenz���s reference state, while being computationally cheaper.

The ocean's gravitational potential energy budget in a coupled climate model

This study examines, in a uni���ed fashion, the budgets of ocean gravitational potential energy (GPE) and available gravitational potential energy (AGPE) in the control simulation of the coupled atmosphere���ocean general circulation model HadCM3. Only AGPE can be converted into kinetic energy by adiabatic processes. Diapycnal mixing supplies GPE, but not AGPE, whereas the reverse is true of the combined e���ect of surface buoyancy forcing and convection. Mixing and buoyancy forcing, thus, play complementary roles in sustaining the large scale circulation. However, the largest globally integrated source of GPE is resolved advection (+0.57 TW) and the largest sink is through parameterized eddy transports (-0.82 TW). The e���ect of these adiabatic processes on AGPE is identical to their e���ect on GPE, except for perturbations to both budgets due to numerical leakage exacerbated by non-linearities in the equation of state.

Thermobaric control of gravitational potential energy generation by diapycnal mixing in the deep ocean

Sources and sinks of gravitational potential energy (GPE) play a rate-limiting role in the large scale ocean circulation. A key source is turbulent diapycnal mixing, whereby irre- versible mixing across isoneutral surfaces is enhanced by turbulent straining of these surfaces. This has motivated international observational e���orts to map diapycnal mixing in the global ocean. However, in order to accurately relate the GPE supplied to the large scale circulation by diapycnal mixing to the mixing energy source, it is ���rst necessary to determine the ratio, �� , of the GPE generation rate to the available potential energy dissipation rate associated with turbulent mixing. Here, the link between GPE and hydro- static pressure is used to derive the GPE budget for a com- pressible ocean with a nonlinear equation of state. The role of diapycnal mixing is isolated and from this a global cli- matological distribution of �� is calculated. It is shown that, for a given source of mixing energy, typically three times as much GPE is generated if the mixing takes place in bottom waters rather than in the pycnocline. This is due to GPE destruction by cabbelling in the pycnocline, as opposed to thermobaric enhancement of GPE generation by diapycnal mixing in the deep ocean.

Variational calculations of rovibrational states: a precise high-energy potential surface for HCN [and discussion]

We report the results of variational calculations of the rovibrational energy levels of HCN for J = 0, 1 and 2, where we reproduce all the ca. 100 observed vibrational states for all observed isotopic species, with energies up to 18000 cm$^{-1}$, to about $\pm $1 cm$^{-1}$, and the corresponding rotational constants to about $\pm $0.001 cm$^{-1}$. We use a hamiltonian expressed in internal coordinates r$_{1}$, r$_{2}$ and $\theta $, using the exact expression for the kinetic energy operator T obtained by direct transformation from the cartesian representation. The potential energy V is expressed as a polynomial expansion in the Morse coordinates y$_{i}$ for the bond stretches and the interbond angle $\theta $. The basis functions are built as products of appropriately scaled Morse functions in the bond-stretches and Legendre or associated Legendre polynomials of cos $\theta $ in the angle bend, and we evaluate matrix elements by Gauss quadrature. The hamiltonian matripx is factorized using the full rovibrational symmetry, and the basis is contracted to an optimized form; the dimensions of the final hamiltonian matrix vary from 240 $\times $ 240 to 1000 $\times $ 1000.We believe that our calculation is converged to better than 1 cm$^{-1}$ at 18 000 cm$^{-1}$. Our potential surface is expressed in terms of 31 parameters, about half of which have been refined by least squares to optimize the fit to the experimental data. The advantages and disadvantages and the future potential of calculations of this type are discussed.

Analysing the ability to retain sidechain hydrogen-bonds in mutant proteins

Motivation: Hydrogen bonds are one of the most important inter-atomic interactions in biology. Previous experimental, theoretical and bioinformatics analyses have shown that the hydrogen bonding potential of amino acids is generally satisfied and that buried unsatisfied hydrogen-bond-capable residues are destabilizing. When studying mutant proteins, or introducing mutations to residues involved in hydrogen bonding, one needs to know whether a hydrogen bond can be maintained. Our aim, therefore, was to develop a rapid method to evaluate whether a sidechain can form a hydrogen-bond. Results: A novel knowledge-based approach was developed in which the conformations accessible to the residues involved are taken into account. Residues involved in hydrogen bonds in a set of high resolution crystal structures were analyzed and this analysis is then applied to a given protein. The program was applied to assess mutations in the tumour-suppressor protein, p53. This raised the number of distinct mutations identified as disrupting sidechain-sidechain hydrogen bonding from 181 in our previous analysis to 202 in this analysis.

Experimental and theoretical study of the infrared spectra of BrHI- and BrDI-

Gas phase vibrational spectra of BrHI- and BrDI- have been measured from 6 to 17 mum (590-1666 cm-1) using tunable infrared radiation from the free electron laser for infrared experiments in order to characterize the strong hydrogen bond in these species. BrHI-.Ar and BrDI-.Ar complexes were produced and mass selected, and the depletion of their signal due to vibrational predissociation was monitored as a function of photon energy. Additionally, BrHI- and BrDI- were dissociated into HBr (DBr) and I- via resonant infrared multiphoton dissociation. The spectra show numerous transitions, which had not been observed by previous matrix studies. New ab initio calculations of the potential-energy surface and the dipole moment are presented and are used in variational ro-vibrational calculations to assign the spectral features. These calculations highlight the importance of basis set in the simulation of heavy atoms such as iodine. Further, they demonstrate extensive mode mixing between the bend and the H-atom stretch modes in BrHI- and BrDI- due to Fermi resonances. These interactions result in major deviations from simple harmonic estimates of the vibrational energies. As a result of this new analysis, previous matrix-isolation spectra assignments are reevaluated. (C) 2004 American Institute of Physics.

Analytical potentials for triatomic molecules from spectroscopic data V. Application to HOX (X=F, Cl, Br, I)

Analytical potential energy functions are reported for HOX (X=F, Cl, Br, I). The surface for HOF predicts two metastable minima as well as the equilibrium configuration. These correspond to HFO (bent) and OHF (linear). Ab initio calculations performed for the HOF surface confirm these predictions. Comparisons are drawn between the two sets of results, and a vibrational analysis is undertaken for the hydrogen bonded OHF species. For HOCl, one further minimum is predicted, corresponding to HClO (bent), the parameters for which compare favourably with those reported from ab initio studies. In contrast, only the equilibrium configurations are predicted to be stable for HOBr and HOI.