Spin contamination in XOO radicals X = F, Cl, Br, HO: How is the investigation of the HOOO radical affected

We have analysed the effect of spin contamination in the wavefunction of HOOO. At least, two solutions can be found for the HF wavefunction. One, lower in energy, presents a high spin contamination and gives qualitatively incorrect structural parameters. On the other hand, the less contaminated HF reference gives structural parameters that are in better agreement with experiment, and positive spin densities on all atoms. Some of the problems described during previous investigations of HOOO can now be traced to problems in the HF reference. For the first time we report a CCSD(T) estimation of the structure of HOOO cis employing a HF reference with small spin contamination. (C) 2008 Elsevier B.V. All rights reserved.

Thiothionyl bromide (S=SBr2): The elusive isomer of dibromodisulfane (BrSSBr)

CCSD(T) with a series of correlation consistent basis up to quadruple-zeta is used to investigate the structures, vibrational spectra, relative stability, heats of formation, and barrier to isomerization of S=SBr2 and BrSSBr. It represents the most accurate and detailed characterization of these molecules to date. We show that the frequency mode at 302 cm(-1), detected in various studies and assigned to impurities by some authors, and to the anti-symmetric SBr stretch in BrSSBr by others, thus in fact corresponds to the anti-symmetric SBr stretch in the elusive S=SBr2 species; it thus corroborates and complements an earlier partial IR spectra study attributable to S=SBr2. Including corrections for relativistic and core-valence correlation effects, we also predict 26.33 (12.74) kcal/mol for Delta H-f (298.15 K) of S=SBr2 (BrSSBr). For the S=SBr2 -> BrSSBr reaction, our best estimates for the Gibbs free energy and the enthalpy of the reaction at 298.15 K are -13.71 and -13.44 kcal/mol, respectively. For a value of Delta G(#) equal to 23.52 kcal/mol, we estimate a TST rate constant, at 298.15 K, of 3.57 x 10(-5) s(-1). (c) 2007 Elsevier B.V. All rights reserved.

Spectroscopic and theoretical studies of some N,N-diethyl-2-[(4 `-substituted)phenylsulfonyl]acetamides

The analysis of the IR carbonyl band of the N,N-diethyl-2-[(4`-substituted)phenylsulfonyl]acetamides Et(2)NC(O)CH(2)S(O)(2)-C(6)H(4)-Y (Y = OMe 1, Me 2,1-13, Cl 4, Br 5, NO(2) 6) supported by B3LYP/6-31G(d,p) calculations for 3, indicated the existence of three pairs (anti and syn) of cis (c) and gauche (g(1) and g(2)) conformers in the gas phase, being the gauche conformers significantly more stable than the cis ones. The anti geometry is more stable than the syn one, for each pair of cis and gauche conformers. The summing up of the orbital (NBO analysis) and electrostatic interactions justifies quite well the populations and the v(CO) frequencies of the anti and syn pairs of c, g(1) and g(2) conformers. The IR higher carbonyl frequency component whose population is ca. 10%, in CCl(4), may be ascribed to the least stable and most polar cis conformer pair (in the gas phase) and the lower frequency component whose population is ca. 90%, to the summing up of the populations of the two most stable and least polar gauche conformer pairs (g(1) and g(2)) (in the gas phase). The reversal of the cis(c)/gauche (g(1) + g(2)) population ratio observed in chloroform ca. 60% (cis)/40% (gauche) and the occurrence of the most polar cis(c) conformer only, in acetonitrile, strongly suggests the coalescence of the two gauche components in a unique carbonyl band in solution. A further support to this rationalization is given by the single point PCM solvation model performed by HF/6-31G(d,p) method, which showed a progressive increase of the c/(g(1) + g(2)) ratio going from gas to CCl(4), to CHCl(3) and to CH(3)CN. X-ray single crystal analysis of 4 indicates that this compound assumes, in the solid state, the syn-clinal (gauche) conformation with respect to the [O=C-CH(2)-S] moiety, and the most stable anti geometry relative to the [C(O)N(CH(2)CH(3))(2)] fragment. In order to obtain larger energy gain from the crystal packing the molecules of 4 are linked in centrosymmetric dimers through two C-H center dot center dot center dot O interactions (C-H([O-Ph])center dot center dot center dot O([SO2])) forming a step ladder. (C) 2011 Elsevier B.V. All rights reserved.

Spectroscopic and theoretical studies of some 2-ethylsulfinyl-(4 `-substituted)-phenylthioacetates

The analysis of the IR nu(co) bands of the 2-ethylsulfinyl-(4`-substituted)-phenylthioacetates 4`-Y-C(6)H(4)SC(O)CH(2)S(O)Et (Y = NO(2) 1, Cl 2, Br 3, H 4, Me 5, OMe 6) supported by B3LY/6-31G(d,p) calculations along with the NBO analysis for 1.4 and 6 and X-ray analysis for 3, indicated the existence of four gauche (q-g-syn, g(3)-syn. g(1)-atin and q-g(2)-syn) conformers for 1-6 The calculations reproduce quite well the experimental results, i e the computed q-g-syn and g3-syn conformers correspond in the IR spectrum (in solution), to the nu(co) doublet higher frequency component of larger intensity, while the computed grant, conformer correspond to the nu(co) doublet lower frequency component (in solution) NBO analysis showed that the n(s) -> pi(center dot)(c1=o2), no(co) -> sigma(c1-s3), no(co) -> sigma(c1-c4) orbital interactions are the main factors which stabilize the q-g-syn, g(3)-syn, g(1)-anti and q-g(2)-syn conformers for 1, 4 and 6 The no(co) -> sigma(c1-s3) interaction which stabilizes the q-g-syn, g(3)-syn and q-g(2)-syn conformers into a larger extent than the granti conformer, is responsible for the larger tto frequencies of the former conformers relative to the latter one. The q-g-syn, g(3)-syn and q-g(2)-syn conformers are further stabilized sigma(c4-s5) -> pi(co)center dot (strong). pi(co)/sigma(c1-c4,) no(co) -> sigma(c6-H17[Et]) (weak) and pi(co)/sigma(c4-c5) pi(co) (strong) orbital interactions. The q-g-syn conformer is also stabilized by sigma(c4-s5) -> pi(center dot)(co) (strong), pi(co)/sigma(c4-c5).no(co) -> sigma(c6-H17[Et]), pi(C9=C11[ph]) -> sigma(c4-H6x-CH2]) (weak). no((SO)) -> sigma(C11-H23[ph]) (medium) pi(co)/sigma(c4-c5)(strong) orbital interactions. The q-g-syn conformei is further stabilized by the n(S5) O((C))(8-) S((SO))(8+) attractive Coulornbic interaction while the q-g(2)-syn conformer is destabilized by the n55 0,8c-0) repulsive Coulombic interaction. This analysis indicates the following conformer stabilization order. q-g-syn, g(3)-syn > g(1)-anti >> q-g(2)-syn X-ray single crystal analysis of 3 indicates that it assumes in the solid a distorted q-g(2)-syn geometry which is stabilized through almost the same orbital and Coulombic interaction which takes place for the q-g(2)-syn conformer, in the gas, along with dipole moment coupling and a series intermolecular C-HO0 interactions. (C) 2010 Elsevier B V All rights reserved

Oxindole-Schiff base copper(II) complexes interactions with human serum albumin: Spectroscopic, oxidative damage, and computational studies

CD and EPR were used to characterize interactions of oxindole-Schiff base copper(II) complexes with human serum albumin (HSA). These imine ligands form very stable complexes with copper, and can efficiently compete for this metal ion towards the specific N-terminal binding site of the protein, consisting of the amino acid sequence Asp-Ala-His. Relative stability constants for the corresponding complexes were estimated from CD data, using the protein as competitive ligand, with values of log K(CuL) in the range 15.7-18.1, very close to that of [Cu(HSA)] itself, with log K(CuHSA) 16.2. Some of the complexes are also able to interfere in the a-helix structure of the protein, while others seem not to affect it. EPR spectra corroborate those results, indicating at least two different metal species in solution, depending on the imine ligand. Oxidative damage to the protein after incubation with these copper(II) complexes, particularly in the presence of hydrogen peroxide, was monitored by carbonyl groups formation, and was observed to be more severe when conformational features of the protein were modified. Complementary EPR spin-trapping data indicated significant formation of hydroxyl and carbon centered radicals, consistent with an oxidative mechanism. Theoretical calculations at density functional theory (DFT) level were employed to evaluate Cu(II)-L binding energies, L -> Cu(II) donation, and Cu(II) -> L back-donation, by considering the Schiff bases and the N-terminal site of HSA as ligands. These results complement previous studies on cytotoxicity, nuclease and pro-apoptotic properties of this kind of copper(II) complexes, providing additional information about their possibilities of transport and disposition in blood plasma. (C) 2009 Elsevier Inc. All rights reserved.

A quantum chemical study on the formation of phosphorus mononitride

The chemical mechanism of the (1)PN formation was successfully studied by using the CCSD(T)/6-311++G(3df,3pd) level of theory. The (1)NH(3) + (3)PH and (4)P + NH(3) reaction paths are not energetically favorable to form the (1)PN molecule. However, the (3)NH + (3)PH, (4)N + (3)PH(3), (4)N + (3)PH, (4)P + (3)NH, and (4)P + (2)NH(2) reaction paths to form the (1)PN molecule are only energetically favorable by taking place through specific transition states to form the (1)PN molecule. The NH(3) + (3)PH, (4)N + (1)PH(3), NH(3) + (4)P, and (4)N + (2)PH(2) reactions are spin-forbidden and the probability of hopping for these reactions was estimated to be 0 by the Landau-Zener theory. This is the first detailed study on the chemical mechanism for the (1)PN formation. (C) 2009 Elsevier B.V. All rights reserved.

Interleukin-22 forms dimers that are recognized by two interleukin-22R1 receptor chains

Interleukin-22 (IL-22) is a class 2 cytokine whose primary structure is similar to that of interleukin 10 (IL-10) and interferon-gamma (IFN-gamma). IL-22 induction during acute phase immune response indicates its involvement in mechanisms of inflammation. Structurally different from IL-10 and a number of other members of IL-10 family, which form intertwined inseparable V-shaped dimers of two identical polypeptide chains, a single polypeptide chain of IL-22 folds on itself in a relatively globular structure. Here we present evidence, based on native gel electrophoresis, glutaraldehyde cross-linking, dynamic light scattering, and small angle x-ray scattering experiments, that human IL-22 forms dimers and tetramers in solution under protein concentrations assessable by these experiments. Unexpectedly, low-resolution molecular shape of IL-22 dimers is strikingly similar to that of IL-10 and other intertwined cytokine dimeric forms. Furthermore, we determine an ab initio molecular shape of the IL-22/IL-22R1 complex which reveals the V-shaped IL-22 dimer interacting with two cognate IL-22R1 molecules. Based on this collective evidence, we argue that dimerization might be a common mechanism of all class 2 cytokines for the molecular recognition with their respective membrane receptor. We also speculate that the IL-22 tetramer formation could represent a way to store the cytokine in nonactive form at high concentrations that could be readily converted into functionally active monomers and dimers upon interaction with the cognate cellular receptors.

Formation of atomic carbon chains from graphene nanoribbons

The formation of one-dimensional carbon chains from graphene nanoribbons is investigated using ab initio molecular dynamics. We show under what conditions it is possible to obtain a linear atomic chain via pulling of the graphene nanoribbons. The presence of dimers composed of two-coordinated carbon atoms at the edge of the ribbons is necessary for the formation of the linear chains, otherwise there is simply the full rupture of the structure. The presence of Stone-Wales defects close to these dimers may lead to the formation of longer chains. The local atomic configuration of the suspended atoms indicates the formation of single and triple bonds, which is a characteristic of polyynes.

Adiabatic intramolecular movements for water systems

An effective treatment of the intramolecular degrees of freedom is presented for water, where these modes are decoupled from the intermolecular ones, ""adiabatically"" allowing these coordinates to be positioned at their local minimum of the potential energy surface. We perform ab initio Monte Carlo simulations with the configurational energies obtained via density functional theory. We study a water dimer as a prototype system, and even in this simple case the intramolecular relaxations are very important to properly describe properties such as the dipole moment. We show that rigid simulations do not correctly sample the phase space, resulting in an average dipole moment smaller than the one obtained with the adiabatic model, which is closer to the experimental result. (c) 2008 American Institute of Physics.

First-principles study of the adsorption of atomic and molecular hydrogen on BC(2)N nanotubes

The adsorption of atomic and molecular hydrogen on armchair and zigzag boron carbonitride nanotubes is investigated within the ab initio density functional theory. The adsorption of atomic H on the BC(2)N nanotubes presents properties which are promising for nanoelectronic applications. Depending on the adsorption site for the H, the Fermi energy moves toward the bottom of the conduction band or toward the top of the valence band, leading the system to exhibit donor or acceptor characteristics, respectively. The H(2) molecules are physisorbed on the BC(2)N surface for both chiralities. The binding energies for the H(2) molecules are slightly dependent on the adsorption site, and they are near to the range to work as a hydrogen storage medium.

Thermodynamic self-consistency issues related to the Cluster Variation Method: The case of the BCC Cr-Fe (Chromium-Iron) system

The Cluster Variation Method (CVM), introduced over 50 years ago by Prof. Dr. Ryoichi Kikuchi, is applied to the thermodynamic modeling of the BCC Cr-Fe system in the irregular tetrahedron approximation, using experimental thermochemical data as initial input for accessing the model parameters. The results are checked against independent data on the low-temperature miscibility gap, using increasingly accurate thermodynamic models, first by the inclusion of the magnetic degrees of freedom of iron and then also by the inclusion of the magnetic degrees of freedom of chromium. It is shown that a reasonably accurate description of the phase diagram at the iron-rich side (i.e. the miscibility gap borders and the Curie line) is obtained, but only at expense of the agreement with the above mentioned thermochemical data. Reasons for these inconsistencies are discussed, especially with regard to the need of introducing vibrational degrees of freedom in the CVM model. (C) 2008 Elsevier Ltd. All rights reserved.

Assignment of the absolute configuration at the sulfur atom of thioridazine metabolites by the analysis of their chiroptical properties: The case of thioridazine 2-sulfoxide

Thioridazine (THD) is a commonly prescribed phenotiazine neuroleptic drug, which is extensively biotransformed in the organism producing as main metabolites sulfoxides and a sulfone by sulfur oxidation Significant differences have been observed in the activity of the THD enantiomers as well as for its main metabolites, and enantioselectivity phenomena have been proved in the metabolic pathway. Here the assignment of the absolute configuration at the sulfur atom of enantiomeric THD-2-sulfoxide (THD-2-SO) has been carried out by circular dichroism (CD) spectroscopy The stereoisomers were separated by HPLC on Chiralpak AS column, recording the CD spectra for the two collected enantiomeric fractions The theoretical electronic CD spectrum has been obtained by the TDDFT/B3LYP/6-31G*. as Boltzmann averaging of the contributions calculated for the most stable conformations of the drug The comparison of the simulated and experimental spectra allowed the absolute configuration at the sulfur atom of the four THD-2-SO stereoisomers to be assigned The developed method should be useful for a reliable correlation between stereochemistry and activity and/or toxicity

Evaluation of the enthalpy of formation, proton affinity, and gas-phase basicity of gamma-butyrolactone and 2-pyrrolidinone by isodesmic reactions

The knowledge of thermochemical parameters such as the enthalpy of formation, gas-phase basicity, and proton affinity may be the key to understanding molecular reactivity. The obtention of these thermochemical parameters by theoretical chemical models may be advantageous when experimental measurements are difficult to accomplish. The development of ab initio composite models represents a major advance in the obtention of these thermochemical parameters,. but these methods do not always lead to accurate values. Aiming at achieving a comparison between the ab initio models and the hybrid models based on the density functional theory (DFT), we have studied gamma-butyrolactone and 2-pyrrolidinone with a goal of obtaining high-quality thermochemical parameters using the composite chemical models G2, G2MP2, MP2, G3, CBS-Q, CBS-4, and CBS-QB3; the DFT methods B3LYP, B3P86, PW91PW91, mPW1PW, and B98; and the basis sets 6-31G(d), 6-31+G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p), 6-311G(d), 6-311+G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p), aug-cc-pVDZ, and aug-cc-pVTZ. Values obtained for the enthalpies of formation, proton affinity, and gas-phase basicity of the two target molecules were compared to the experimental data reported in the literature. The best results were achieved with the use of DFT models, and the B3LYP method led to the most accurate data.

Excited state electronic polarization and reappraisal of the n <- pi* emission of acetone in water

Electronic polarization of the acetone molecule in the excited n -> pi* state is considered and its influence on the solvent shift in the emission spectrum is analyzed. Using an iterative procedure the electronic polarizations of both the ground and the excited states are included and compared with previous results obtained with Car-Parrinello dynamics. Analysis of the emission transition obtained using CIS(D)/aug-cc-pVDZ on statistically uncorrelated solute-solvent structures, composed of acetone and twelve explicit water molecules embedded in the electrostatic field of remaining 263 water molecules, corroborates that the solvent effect is mild, calculated here between 80 and 380 cm (1). (c) 2010 Published by Elsevier B.V.

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