Thermodynamic Stabillity of Hydrogen-Bonded Systems in Polar and Nonpolar Environments

The thermodynamic properties of a selected set of benchmark hydrogen-bonded systems (acetic acid dimer and the complexes of acetic acid with acetamide and methanol) was studied with the goal of obtaining detailed information on solvent effects on the hydrogen-bonded interactions using water, chloroform, and n-heptane as representatives for a wide range in the dielectric constant. Solvent effects were investigated using both explicit and implicit solvation models. For the explicit description of the solvent, molecular dynamics and Monte Carlo simulations in the isothermal isobaric (NpT) ensemble combined with the free energy perturbation technique were performed to determine solvation free energies. Within the implicit solvation approach, the polarizable continuum model and the conductor-like screening model were applied. Combination of gas phase results with the results obtained from the different solvation models through an appropriate thermodynamic cycle allows estimation of complexation free energies, enthalpies, and the respective entropic contributions in solution. Owing to the strong solvation effects of water the cyclic acetic acid dimer is not stable in aqueous solution. In less polar solvents the double hydrogen bond structure of the acetic acid dimer remains stable. This finding is in agreement with previous theoretical and experimental results. A similar trend as for the acetic acid dimer is also observed for the acetamide complex. The methanol complex was found to be thermodynamically unstable in gas phase as well as in any of the three solvents. (C) 2010 Wiley Periodicals, Inc. J Comput Chem 31: 2046-2055, 2010

Electronic properties of liquid hydrogen fluoride: A sequential quantum mechanical/Born-Oppenheimer molecular dynamics approach

The electronic properties of liquid hydrogen fluoride (HF) were investigated by carrying out sequential quantum mechanics/Born-Oppenheimer molecular dynamics. The structure of the liquid is in good agreement with recent experimental information. Emphasis was placed on the analysis of polarisation effects, dynamic polarisability and electronic excitations in liquid HF. Our results indicate an increase in liquid phase of the dipole moment (similar to 0.5 D) and isotropic polarisability (5%) relative to their gas-phase values. Our best estimate for the first vertical excitation energy in liquid HF indicates a blue-shift of 0.4 +/- 0.2 eV relative to that of the gas-phase monomer (10.4 eV). (C) 2010 Elsevier B.V. All rights reserved.

Many-body energy decomposition of hydrogen-bonded glycine clusters in gas-phase

A detailed analysis of the many-body contribution to the interaction energies of the gas-phase hydrogen-bonded glycine clusters, (Gly)(N), N = 1-4 is presented. The energetics of the hydrogen-bonded dimer, trimer and tetramer complexes have been analyzed using density-functional theory. The magnitude of the two-through four-body energy terms have been calculated and compared. The relaxation energy and the two-body energy terms are the principal contributors to the total binding energy. Four-body contribution is negligible. However, the three-body contribution is found to be sizable and the formation of the cyclic glycine trimer presents geometric strains that make it less favorable. (C) 2010 Elsevier B.V. All rights reserved.

A model for the stabilization of atomic hydrogen centers in borate glasses

A previously proposed model describing the trapping site of the interstitial atomic hydrogen in borate glasses is analyzed. In this model the atomic hydrogen is stabilized at the centers of oxygen polygons belonging to B-O ring structures in the glass network by van der Waals forces. The previously reported atomic hydrogen isothermal decay experimental data are discussed in the light of this microscopic model. A coupled differential equation system of the observed decay kinetics was solved numerically using the Runge Kutta method. The experimental untrapping activation energy of 0.7 x 10(-19) J is in good agreement with the calculated results of dispersion interaction between the stabilized atomic hydrogen and the neighboring oxygen atoms at the vertices of hexagonal ring structures. (C) 2009 Elsevier B.V. All rights reserved.

Pressure-induced electrical and structural anomalies in Pb(1-x)Ca(x)TiO(3) thin films grown at various oxygen pressures by chemical solution route

Lead calcium titanate (Pb(1-x)Ca(x)TiO(3) or PCT) thin films have been thermally treated under different oxygen pressures, 10, 40 and 80 bar, by using the so-called chemical solution deposition method. The structural, morphological, dielectric and ferroelectric properties were characterized by x-ray diffraction, FT-infrared and Raman spectroscopy, atomic force microscopy and polarization-electric-field hysteresis loop measurements. By annealing at a controlled pressure of around 10 and 40 bar, well-crystallized PCT thin films were successfully prepared. For the sample submitted to 80 bar, the x-ray diffraction, Fourier transformed-infrared and Raman data indicated deviation from the tetragonal symmetry. The most interesting feature in the Raman spectra is the occurrence of intense vibrational modes at frequencies of around 747 and 820 cm(-1), whose presence depends strongly on the amount of the pyrochlore phase. In addition, the Raman spectrum indicates the presence of symmetry-breaking disorder, which would be expected for an amorphous (disorder) and mixed pyrochlore-perovskite phase. During the high-pressure annealing process, the crystallinity and the grain size of the annealed film decreased. This process effectively suppressed both the dielectric and ferroelectric behaviour. Ferroelectric hysteresis loop measurements performed on these PCT films exhibited a clear decrease in the remanent polarization with increasing oxygen pressure.

Hydrogen adsorption on boron doped graphene: an ab initio study

(i) The electronic and structural properties of boron doped graphene sheets, and (ii) the chemisorption processes of hydrogen adatoms on the boron doped graphene sheets have been examined by ab initio total energy calculations. In (i) we find that the structural deformations are very localized around the boron substitutional sites, and in accordance with previous studies (Endo et al 2001 J. Appl. Phys. 90 5670) there is an increase of the electronic density of states near the Fermi level. Our simulated scanning tunneling microscope (STM) images, for occupied states, indicate the formation of bright (triangular) spots lying on the substitutional boron (center) and nearest-neighbor carbon (edge) sites. Those STM images are attributed to the increase of the density of states within an energy interval of 0.5 eV below the Fermi level. For a boron concentration of similar to 2.4%, we find that two boron atoms lying on the opposite sites of the same hexagonal ring (B1-B2 configuration) represents the energetically most stable configuration, which is in contrast with previous theoretical findings. Having determined the energetically most stable configuration for substitutional boron atoms on graphene sheets, we next considered the hydrogen adsorption process as a function of the boron concentration, (ii). Our calculated binding energies indicate that the C-H bonds are strengthened near boron substitutional sites. Indeed, the binding energy of hydrogen adatoms forming a dimer-like structure on the boron doped B1-B2 graphene sheet is higher than the binding energy of an isolated H(2) molecule. Since the formation of the H dimer-like structure may represent the initial stage of the hydrogen clustering process on graphene sheets, we can infer that the formation of H clusters is quite likely not only on clean graphene sheets, which is in consonance with previous studies (Hornekaer et al 2006 Phys. Rev. Lett. 97 186102), but also on B1-B2 boron doped graphene sheets. However, for a low concentration of boron atoms, the formation of H dimer structures is not expected to occur near a single substitutional boron site. That is, the formation (or not) of H clusters on graphene sheets can be tuned by the concentration of substitutional boron atoms.

2-(4,5,6,7-tetrafluorobenzimidazol-2-yl)-4,4,5,5-tetramethyl-4,5-dihydro-1H-imidazole-3-oxide-1-oxyl, a hydrogen-bonded organic quasi-1D ferromagnet

The title radical (F4BlmNN) is a stable nitronylnitroxide that forms hydrogen-bonded NH center dot center dot center dot ON chains in the solid state. The chains assemble the F4BlmNN molecules to form stacked contacts between the radical groups, in a geometry that is expected to exhibit ferromagnetic (FM) exchange based on spin polarization (SP) models. The experimental magnetic susceptibility of F4BlmNN confirms the expectation, showing 1-D Heisenberg chain FM exchange behavior over 1.8-300 K with an intrachain exchange constant Of J(chain)/k = +22 K. At lower temperatures, ac magnetic susceptibility and variable field heat capacity measurements show that F4BlmNN acts as a quasi-1-D ferromagnet. The dominant ferromagnetic exchange interaction is attributable to overlap between spin orbitals of molecules within the hydrogen-bonded chains, consistent with the SP model expectations. The chains appear to be antiferromagnetically exchange coupled, giving cusps in the ac susceptibility and zero field heat capacity at lower temperatures. The results indicate that the sample orders magnetically at about 0.7 K. The magnetic heat capacity ordering cusp shifts to lower temperatures as external magnetic field increases, consistent with forming a bulk antiferromagnetic phase below a Neel temperature of T-N(0) = 0.72 K, with a critical field of H-c approximate to 1800 Oe. The interchain exchange is estimated to be zJ/k congruent to (-)0.1 K.

Electrical transport in disordered and ordered magnetic domains under pressures and magnetic fields

Magnetic M( T, H, P) and electrical transport.( T, H, P) measurements in a strong spin-lattice-charge coupled La(0.7)Ca(0.3)MnO(3) system have been conducted. The application of H and P leads to the formation of different magnetic domain structures in the vicinity and below the polaronic-to-ferromagnetic transition temperature. The charge mobility is more sensitive to the variation of the spatial wave function overlap between Mn(3+) eg and O(2-) 2p orbitals due to the applied compacting pressure rather than the relative spin orientation between neighbouring Mn ions when the magnetic field is applied. In spite of the presence of different magnetic domain structures due to the sample history, the effect of magnetic field and pressure is less pronounced at lower temperatures on electrical transport properties.

Dirac equation in noncommutative space for hydrogen atom

We consider the energy levels of a hydrogen-like atom in the framework of theta-modified, due to space noncommutativity, Dirac equation with Coulomb field. It is shown that on the noncommutative (NC) space the degeneracy of the levels 2S(1/2), 2P(1/2) and 2P(3/2) is lifted completely, such that new transition channels are allowed. (C) 2009 Elsevier B.V. All rights reserved.

Large Changes of Static Electric Properties Induced by Hydrogen Bonding: An ab Initio Study of Linear HCN Oligomers

We report the partitioning of the interaction-induced static electronic dipole (hyper)polarizabilities for linear hydrogen cyanide complexes into contributions arising from various interaction energy terms. We analyzed the nonadditivities of the studied properties and used these data to predict the electric properties of an infinite chain. The interaction-induced static electric dipole properties and their nonadditivities were analyzed using an approach based on numerical differentiation of the interaction energy components estimated in an external electric field. These were obtained using the hybrid variational-perturbational interaction energy decomposition scheme, augmented with coupled-cluster calculations, with singles, doubles, and noniterative triples. Our results indicate that the interaction-induced dipole moments and polarizabilities are primarily electrostatic in nature; however, the composition of the interaction hyperpolarizabilities is much more complex. The overlap effects substantially quench the contributions due to electrostatic interactions, and therefore, the major components are due to the induction and exchange induction terms, as well as the intramolecular electron-correlation corrections. A particularly intriguing observation is that the interaction first hyperpolarizability in the studied systems not only is much larger than the corresponding sum of monomer properties, but also has the opposite sign. We show that this effect can be viewed as a direct consequence of hydrogen-bonding interactions that lead to a decrease of the hyperpolarizability of the proton acceptor and an increase of the hyperpolarizability of the proton donor. In the case of the first hyperpolarizability, we also observed the largest nonadditivity of interaction properties (nearly 17%) which further enhances the effects of pairwise interactions.

N-H...S=C hydrogen bond in O-alkyl N-methoxycarbonyl thiocarbamates, ROC(S)N(H)C(O)O)CH(3) (R = CH(3)(-), CH(3)CH(2)-)

Pure O-methyl N-methoxycarbonyl thiocarbamate CH(3)OC(S)N(H)C(O)OCH(3) (I) and O-ethyl N-methoxycarbonyl thiocarbamate, CH(3)CH(2)OC(S)N(H)C(O)OCH(3) (II), are quantitatively prepared by the addition reaction between the CH(3)OC(O)NCS and the corresponding alcohols. The compounds are characterized by multinuclear ((1)H and (13)C) and bi-dimensional ((13)C HSQC) NMR, GC-MS and FTIR spectroscopy techniques. Structural and conformational properties are analyzed using a combined approach involving crystallographic data, vibration spectra and theoretical calculations. The low-temperature (150 K) crystal structure of II was determined by X-ray diffraction methods. The substance crystallizes in the monoclinic space group P2(1)/n with a = 4.088(1)angstrom. b = 22.346(1)angstrom, c = 8.284(1)angstrom, beta = 100.687(3)degrees and Z = 4 molecules per unit cell. The conformation adopted by the thiocarbamate group -OC(S)N(H)- is syn (C=S double bond in synperiplanar orientation with respect to the N-H single bond), while the methoxycarbonyl C=O double bond is in antiperiplanar orientation with respect to the N-H bond. The non-H atoms in II are essentially coplanar and the molecules are arranged in the crystal lattice as centro-symmetric dimeric units held by N-H center dot center dot center dot S=C hydrogen bonds Id(N center dot center dot center dot S) = 3.387(1)angstrom, <(N-H center dot center dot center dot S) = 166.4(2)degrees]. Furthermore, the effect of the it electronic resonance in the structural and vibrational properties is also discussed. (C) 2009 Elsevier Ltd. All rights reserved.

Structural effects and hydrogen bonds on N,N `-di(methoxycarbonylsulfenyl)urea, [CH(3)OC(O)SNH](2)CO, studied by experimental and theoretical methods

Pure N,N`-di(methoxycarbonylsulfenyl)urea, [CH(3)OC(O)SNH](2)CO, is quantitatively prepared by the hydrolysis reaction of CH(3)OC(O)SNCO and characterized by (1)H NMR, GC-MS and FTIR spectroscopy techniques. Structural and conformational properties are analyzed using a combined approach with data obtained from X-ray diffraction, vibrational spectra and theoretical calculation methods. The IR and Raman spectra for normal and deuterated species are reported. The crystal structure of [CH(3)OC(O)SNH](2)CO was determined by X-ray diffraction methods. The substance crystallizes in the orthorhombic P2(1)2(1)2 space group with a = 9.524(2), b = 12.003(1), c = 4.481 (1) angstrom, and Z = 2 moieties in the unit cell. The molecule is sited on a twofold crystallographic axis (C(2)) parallel to c and shows the anti-anti conformation (S-N single bonds antiperiplanar with respect to the opposite C-N single bonds in sulfenyl-urea-sic group). Neighboring molecules are arranged in a chain motif that extends along the C(2)-axis and is held by bifurcated NH center dot center dot center dot O center dot center dot center dot HN intermolecular bonds. A local planar symmetry is observed in the crystal for the central -SN(H)C(O)N(H)S- skeleton. Experimental and calculated data allow to trace this structural feature to the occurrence of N-H center dot center dot center dot O=C hydrogen bonding interactions. Calculated vibrational and structural properties are in good agreement with the experimentally determined features. (C) 2008 Elsevier B.V. All rights reserved.

Analysis of Agonist and Antagonist Effects on Thyroid Hormone Receptor Conformation by Hydrogen/Deuterium Exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T(3)) or antagonist (NH(3)). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, C-terminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T(3), but not NH(3), increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T(3) but not NH(3.) We present data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011)

The Determination of Total Hydrogen, Carbon, Nitrogen and Sulfur in Silicates, Silicate Rocks, Soils and Sediments

Determinations of the volatile elements carbon, hydrogen, sulfur and nitrogen in many geological RM, performed with the LECO CHN and SC analysers, are presented. The method allowed the determination of S in concentrations from a few % m/m to 0.001% m/m or less, of C from % m/m to 0.01% m/m and of H from % m/m to 0.004% m/m. Accuracy was usually better than the XRF method (for S). All obtained values passed the Sutarno-Steger test, which establishes that vertical bar(mean(analysed) - mean(certified))vertical bar/ S(certified) < 2, for the cases with an appropriate number of determinations (n > 10 for each element). It was possible to perform routine determination of C, H and S with the instrumentation, coupled with the determination of major and minor elements in geological materials. Determination of nitrogen could also be performed on an exploratory basis, with improvements in the method dependent on the future availability of more reference materials with reliable composition of this element.

Discontinuous Galerkin approximation of two-phase flows in heterogeneous porous media with discontinuous capillary pressures

We design and investigate a sequential discontinuous Galerkin method to approximate two-phase immiscible incompressible flows in heterogeneous porous media with discontinuous capillary pressures. The nonlinear interface conditions are enforced weakly through an adequate design of the penalties on interelement jumps of the pressure and the saturation. An accurate reconstruction of the total velocity is considered in the Raviart-Thomas(-Nedelec) finite element spaces, together with diffusivity-dependent weighted averages to cope with degeneracies in the saturation equation and with media heterogeneities. The proposed method is assessed on one-dimensional test cases exhibiting rough solutions, degeneracies, and capillary barriers. Stable and accurate solutions are obtained without limiters. (C) 2010 Elsevier B.V. All rights reserved.