Molecular dynamics simulation of liquid trimethylphosphine

Structural and dynamical properties of liquid trimethylphosphine (TMP), (CH(3))(3)P, as a function of temperature is investigated by molecular dynamics (MD) simulations. The force field used in the MD simulations, which has been proposed from molecular mechanics and quantum chemistry calculations, is able to reproduce the experimental density of liquid TMP at room temperature. Equilibrium structure is investigated by the usual radial distribution function, g(r), and also in the reciprocal space by the static structure factor, S(k). On the basis of center of mass distances, liquid TMP behaves like a simple liquid of almost spherical particles, but orientational correlation due to dipole-dipole interactions is revealed at short-range distances. Single particle and collective dynamics are investigated by several time correlation functions. At high temperatures, diffusion and reorientation occur at the same time range as relaxation of the liquid structure. Decoupling of these dynamic properties starts below ca. 220 K, when rattling dynamics of a given TMP molecules due to the cage effect of neighbouring molecules becomes important. (C) 2011 American Institute of Physics. [doi: 10.1063/1.3624408]

Structure, stability, depolarized light scattering, and vibrational spectra of fullerenols from all-electron density-functional-theory calculations

We have investigated the stability, electronic properties, Rayleigh (elastic), and Raman (inelastic) depolarization ratios, infrared and Raman absorption vibrational spectra of fullerenols [C(60)(OH)(n)] with different degrees of hydroxylation by using all-electron density-functional-theory (DFT) methods. Stable arrangements of these molecules were found by means of full geometry optimizations using Becke's three-parameter exchange functional with the Lee, Yang, and Parr correlation functional. This DFT level has been combined with the 6-31G(d,p) Gaussian-type basis set, as a compromise between accuracy and capability to treat highly hydroxylated fullerenes, e.g., C(60)(OH)(36). Thus, the molecular properties of fullerenols were systematically analyzed for structures with n=1, 2, 3, 4, 8, 10, 16, 18, 24, 32, and 36. From the electronic structure analysis of these molecules, we have evidenced an important effect related to the weak chemical reactivity of a possible C(60)(OH)(24) isomer. To investigate Raman scattering and the vibrational spectra of the different fullerenols, frequency calculations are carried out within the harmonic approximation. In this case a systematic study is only performed for n=1-4, 8, 10, 16, 18, and 24. Our results give good agreements with the expected changes in the spectral absorptions due to the hydroxylation of fullerenes.

HFE gene mutations in patients with primary iron overload: Is there a significant improvement in molecular diagnosis yield with HFE sequencing?

Rare HFE variants have been shown to be associated with hereditary hemochromatosis (HH), an iron overload disease. The low frequency of the HFE p.C282Y mutation in HH-affected Brazilian patients may suggest that other HFE-related mutations may also be implicated in the pathogenesis of HH in this population. The main aim was to screen for new HFE mutations in Brazilian individuals with primary iron overload and to investigate their relationship with HH. Fifty Brazilian patients with primary iron overload (transferrin saturation >50% in females and 60% in males) were selected. Subsequent bidirectional sequencing for each HFE exon was performed. The effect of HFE mutations on protein structure were analyzed by molecular dynamics simulation and free binding energy calculations. p.C282Y in homozygosis or in heterozygosis with p.H63D were the most frequent genotypic combinations associated with HH in our sample population (present in 17 individuals, 34%). Thirty-six (72.0%) out of the 50 individuals presented at least one HFE mutation. The most frequent genotype associated with HH was the homozygous p.C282Y mutation (n = 11, 22.0%). One novel mutation (p.V256I) was indentified in heterozygosis with the p.H63D mutation. In silico modeling analysis of protein behavior indicated that the p.V256I mutation does not reduce the binding affinity between HFE and beta 2-microglobulin ((beta 2M) in the same way the p.C282Y mutation does compared with the native HFE protein. In conclusion, screening of HFE through direct sequencing, as compared to p.C282Y/p.H63D genotyping, was not able to increase the molecular diagnosis yield of HH. The novel p.V256I mutation could not be implicated in the molecular basis of the HH phenotype, although its role cannot be completely excluded in HH-phenotype development. Our molecular modeling analysis can help in the analysis of novel, previously undescribed, HFE mutations. (C) 2010 Elsevier Inc. All rights reserved.

Calculations of vibrational frequencies, Raman activities and degrees of depolarization for complexes involving water, methanol and ethanol

Raman activities and degrees of depolarization are reported for 14 complexes involving methanol, ethanol and water using the MP2/aug-cc-pVDZ model. For ethanol both trans and gauche isomers are considered. The red-shifts of the OH stretching and the blue shifts of the bending tau(CO-OH) mode were analyzed for the proton-donor molecules upon hydrogen bond. The shift of the nu(CO) stretching mode of the alcohol molecules are also analyzed and found to be specific giving characterization of the amphoteric relation, being positive for the proton-acceptor and negative for the proton-donor molecule. (c) 2008 Elsevier B.V. All rights reserved.

Ab Initio Calculations of Structural Evolution and Conductance of Benzene-1,4-dithiol on Gold Leads

By performing at) initio density functional theory (DFT) calculations and electronic transport simulations based on the OFT nonequilibrium Green`s functions method we investigate how the conformational changes of a benzene-1,4-dithiol molecule bonded to gold affect the molecular transport as the electrodes are separated from each other. In particular we consider the full evolution of the stretching process until the Junction breaking point and compare results obtained with a standard semilocal exchange and correlation functional to those computed with a self-interaction corrected method. We conclude that the inclusion of self-interaction corrections is fundamental for describing both the molecule conductance and its stability against conformational fluctuations.

Degenerate Two-Photon Absorption in All-Trans Retinal: Nonlinear Spectrum and Theoretical Calculations

In this work we investigate the degenerate two-photon absorption spectrum of all-trans retinal ill ethanol employing the Z-scan technique with femtosecond pulses, The two-photon absorption (2PA) spectrum presents a monotonous increase as the excitation wavelength approaches the one-photon absorption band and it peak at 790 nm. We attribute the 2PA hand to the mixing of states (1)B(u)+-like and vertical bar S(1)>, which are strongly allowed by one- and two-photon, respectively. We modeled the 2PA spectrum by using the sum-over-states approach and obtained spectroscopic parameters of the electronic transitions to vertical bar S >, vertical bar S(2)> (""(1)Bu(+)""), vertical bar S(3)>, and vertical bar S(4)> singlet-excited states. The results were compared with theoretical predictions of one- and two-photon transition calculations using the response Functions formalism within the density functional theory framework with the aid of the CAM-B3LYP functional.

Combined Crystallographic and Solution Molecular Dynamics Study of Allosteric Effects in Ester and Ketone p-tert-Butylcalix[4]arene Derivatives and Their Complexes with Acetonitrile, Cd(II), and Pb(II)

We describe here a procedure to bridge the gap in the field of calixarene physicochemistry between solid-state atomic-resolution structural information and the liquid-state low-resolution thermodynamics and spectroscopic data. We use MD simulations to study the kinetics and energetics involved in the complexation of lower rim calix[4]arene derivatives (L), containing bidentate ester (1) and ketone (2) pendant groups, with acetonitrile molecule (MeCN) and Cd2+ and Pb2+ ions (M2+) in acetonitrile solution. On one hand, we found that the prior inclusion of MeCN into the calix to form a L(MeCN) adduct has only a weak effect in preorganizing the hydrophilic cavity toward metal ion binding. On the other hand, the strong ion-hydrophilic cavity interaction produces a wide open calix which enhances the binding of one MeCN molecule (allosteric effect) to stabilize the whole (M2+)1(MeCN) bifunctional complex. We reach two major conclusions: (i) the MD results for the (M2+)1(MeCN) binding are in close agreement with the ""endo"", fully encapsulated, metal complex found by X-ray diffraction and in vacuo MD calculations, and (ii) the MD structure for the more flexible 2 ligand, however, differs from the also endo solid-state molecule. In fact, it shows strong solvation effects at the calixarene lower bore by competing MeCN molecules that share the metal coordination sphere with the four C=O oxygens of an ""exo"" (M2+)2(MeCN) complex.

Vibrational properties, crystal X-ray diffraction structure and quantum chemical calculations on a divalent sulfur substituted phthalimide: 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio]

Structural and conformational properties of 1H-Isoindole-1,3(2H)-dione, 2-[(methoxycarbonyl)thio] (S-phthalimido O-methyl thiocarbonate) are analyzed using a combined approach including X-ray diffraction, vibrational spectra and theoretical calculation methods. The vibrational properties have been studied by infrared and Raman spectroscopies along with quantum chemical calculations (B3LYP and B3PW91 functional in connection with the 6-311++G** and aug-cc-pVDZ basis sets). The crystal structure was determined by X-ray diffraction methods. The substance crystallizes in the monoclinic P2(1)/c space group with a = 6.795(1), b = 5.109(1), c = 30.011(3) angstrom, beta = 90.310(3)degrees and Z = 4 molecules per unit cell. The conformation adopted by the N-S-C=O group is syn (C=O double bond in synperiplanar orientation with respect to the N-S single bond). The experimental molecular structure is well reproduced by the MP2/aug-cc-pVDZ method. (C) 2010 Elsevier B.V. All rights reserved.

Experimental and theoretical investigation of molecular structure and conformation of the 4-isopropylthioxanthone

In this study, the molecular structure and conformational analyses of the 4-isopropylthioxanthone (4-ITX) are reported according to experimental and theoretical results. The compound crystallizes in the centrosymmetric P (1) over bar space group with only one molecule in the asymmetric unit, presenting the most stable conformation, in which the three fused-rings adopt a planar geometry, and the isopropyl group assumes a torsional angle with less sterical hindrance. The structural and conformational analyses were performed using theoretical calculations such as Hartree-Fock (HF), DFT method in combination with 6-311G(d,p) and 6-31++G(d,p) and the results were compared with infrared spectroscopy (FT-IR) and X-ray diffraction (XRD). The supramolecular assembly of 4-ITX is kept by non-classical C-H center dot center dot center dot O hydrogen bonds and weak interactions such as pi-pi stacking. 4-ITX was also studied by (1)H and (13)C NMR spectroscopy. UV-Vis absorption spectroscopic properties of the 4-ITX showed the long-wavelength maximum shifts towards high energy when the solvent polarity increases. (C) 2011 Elsevier B.V. All rights reserved.

Predicting and characterizing new molecular species at a high-level theoretical approach: HSeCl, HClSe, and SeCl

New molecular species HSeCl, HClSe, and SeCl were investigated at a high level of theory, CCSD(T), with a series of correlation consistent basis sets with extrapolation to the CBS limit. Account has been taken for valence-only and core-valence correlation effects, and of anharmonic effects on the vibrational frequencies. HSeCl is 43.25 kcal mol (1) more stable than HClSe. A barrier (Delta G(#)) of 47.20 kcal mol (1) separates these species. Internuclear distances are generally overestimated by 0.008 angstrom in the valence-only correlation calculations. Inclusion of anharmonicity leads to much improved vibrational frequencies. For SeCl, we estimate Delta H(f) (0 K) = 23.96 and Delta H(f) (298.15 K) = 24.64 kcal mol (1); for HSeCl, we had 4.20 and 4.97 kcal mol (1), respectively. (C) 2010 Elsevier B.V. All rights reserved.

Identification of species formed after pyridine adsorption on iron, cobalt, nickel and silver electrodes by SERS and theoretical calculations

The adsorption of pyridine (py) on Fe, Co, Ni and Ag electrodes was studied using surface-enhanced Raman scattering (SERS) to gain insight into the nature of the adsorbed species. The wavenumber values and relative intensities of the SERS bands were compared to the normal Raman spectrum of the chemically prepared transition metal complexes. Raman spectra of model clusters M(4)(py) (four metal atoms bonded to one py moiety) and M(4)(alpha-pyridil) where M = Ag, Fe, Co or Ni were calculated by density functional theory (DFT) and used to interpret the experimental SERS results. The similarity of the calculated M(4)(py) spectra with the experimental SERS spectra confirm the molecular adsorption of py on the surface of the metallic electrodes. All these results exclude the formation of adsorbed alpha-pyridil species, as suggested previously. Copyright (C) 2009 John Wiley & Sons, Ltd.

Exploring new molecular species on the (1)[H, Se, F] singlet potential energy surface: Energetics, structures, IR spectra, and heats of formation

Structural, energetic, and vibrational properties of new molecular species, HSeF and HFSe, the associated transition state, and dissociation fragments are investigated using a state-of-the-art theoretical approach, CCSD(T)/CBS. HSeF is a normal covalently bonded molecule 38.98 kcal mol (1) more stable than the complex HF-Se, which shows an unusual structure with a central fluorine atom and a bond angle of 101.8 degrees.A barrier (Delta G(#)) of 49.01 kcal mol (1) separates the two species. Vibrational frequencies are also quite distinct. Heats of formation are evaluated for the diatomic fragments and HSeF. Final Delta(f)H values depend on the experimental accuracy of those of Se(g) and H(2)Se. (c) 2009 Elsevier B.V. All rights reserved.

The intramolecular charge transfer in a donor-pi-acceptor dianion probed by resonance Raman spectroscopy and quantum chemical calculations

The dideprotonation of 4-(4-nitrophenylazo)resorcinol generates an anionic species with substantial electronic pi delocalization. As compared to the parent neutral species, the anionic first excited electronic transition, characterized as an intramolecular charge transfer (ICT) from the CO(-) groups to the NO(2) moiety, shows a drastic red shift of ca. 200 nm in the lambda(max) in the UV-vis spectrum, leading to one of the lowest ICT energies observed (lambda(max) = 630 nm in dimethyl sulfoxide (DMSO)) in this class of push-pull molecular systems. Concomitantly, a threefold increase in the molar absorptivity (epsilon(max)) in comparison to the neutral species is observed. The resonance Raman enhancement profiles reveal that in the neutral species the chromophore involves several modes, as nu(C-N), nu(N=N), nu(C=C) and nu(s)(NO(2)), whereas in the dianion, there is a selective enhancement of the NO(2) vibrational modes. The quantum chemical calculations of the electronic transitions and vibrational wavenumbers led to a consistent analysis of the enhancement patterns observed in the resonance Raman spectra. Copyright (C) 2009 John Wiley & Sons, Ltd.

Predicting new molecular species of potential interest to atmospheric chemistry: The isomers HSI and HIS

Structural, vibrational, and energetic properties of new molecular species, HSI and HIS are investigated for the first time using a state-of-the-art theoretical approach. These molecules can be easily differentiated by their geometric parameters and vibrational spectra. HSI is much more stable, and a direct unimolecular isomerization is very unlikely. Kinetics estimates predict that only at low temperatures there is a possibility of isolating HIS. For HS-I, we estimate a bond dissociation energy of 46.25 kcal/mol, and a heat of formation at 298.15 K of 12.84 kcal/mol. For the H(2)S + I(2) -> HSI + HI reaction enthalpy, we found 8.40 kcal/ mol. (C) 2008 Elsevier B.V. All rights reserved.

Azo-hydrazone tautomerism in protonated aminoazobenzenes: Resonance Raman spectroscopy and quantum-chemical calculations

The protonation effect on the vibrational and electronic spectra of 4-aminoazobenzene and 4-(dimethylamino)azobenzene was investigated by resonance Raman spectroscopy, and the results were discussed on the basis of quantum-chemical calculations. Although this class of molecular systems has been investigated in the past concerning the azo-hydrazone tautomerism, the present work is the first to use CASSCF/CASPT2 calculations to unveil the structure of both tautomers as well the nature of the molecular orbitals involved in chromophoric moieties responsible for the resonance Raman enhancement patterns. More specifically both the resonance Raman and theoretical results show clearly that in the neutral species, the charge transfer transition involves mainly the azo moiety, whereas in the protonated forms there is a great difference, depending on the tautomer. In fact, for the azo tautomer the transition is similar to that observed in the corresponding neutral species, whereas in the hydrazone tautomer such a transition is much more delocalized due to the contribution of the quinoid structure. The characterization of protonated species and the understanding of the tautomerization mechanism are crucial for controlling molecular properties depending on the polarity and pH of the medium.