Realistic calculations of carbon-based disordered systems

Carbon nanotubes rank amongst potential candidates for a new family of nanoscopic devices, in particular for sensing applications. At the same time that defects in carbon nanotubes act as binding sites for foreign species, our current level of control over the fabrication process does not allow one to specifically choose where these binding sites will actually be positioned. In this work we present a theoretical framework for accurately calculating the electronic and transport properties of long disordered carbon nanotubes containing a large number of binding sites randomly distributed along a sample. This method combines the accuracy and functionality of ab initio density functional theory to determine the electronic structure with a recursive Green`s functions method. We apply this methodology on the problem of nitrogen-rich carbon nanotubes, first considering different types of defects and then demonstrating how our simulations can help in the field of sensor design by allowing one to compute the transport properties of realistic nanotube devices containing a large number of randomly distributed binding sites.

Effects of torsional disorder on poly-para-phenylene

We study the effect of thermal disorder on the electronic structure of one-dimensional poly-para-phenylene (PPP). In a real chain the torsion angles between rings are bound to be distributed over a range of values, which depend on temperature, and thus the chain is intrinsically disordered. In this study we simulated this kind of thermally induced off-diagonal disorder through the simple Huckel method. We base our Hamiltonian on ab initio results for the effect of temperature on torsion angles, and the effect of torsion angles on the energy gap. We analyze the electronic structure of 200-monomer-long chains focusing on the density of states, and the associated localization character (measured by the inverse participation ratio). Our results contrast with the usually assumed Gaussian-shaped density of localized states for disordered systems. (C) 2009 Elsevier B.V. All rights reserved.

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.

Spectroscopic and theoretical studies of some N-methoxy-N-methyl-2-[(4 `-substituted) phenylsulfonyl]propanamides

The analysis of the IR carbonyl band of the N-methoxy-N-methyl-2-[(4`-substituted)phenyisulfonyI]propanamides Y-PhSO(2)CH(Me)C(O)N(OMe)Me (Y = OMe 1, Me 2, H 3, Cl 4, NO(2) 5), supported by HF/6-31G(d,p) calculations of 3,, indicated the existence of two gauche conformers (g(1) and g(2)), the g, conformer being the most stable and the least polar one (in gas phase and in solution). Both conformers are present in solution of the non polar solvent (CCl(4)) for 1-5 and in solution of the more polar solvents (CHCl(3)) for 1. 4, 5 and (CH(2)Cl(2)) for 5, while only the g(1) conformer is present in solution of the most polar solvent (CH(3)CN) for 1-5. The g, and g2 conformers correspond to the enantiomeric pairs of diastereomers (diast(1) and diast(2)) whose relative configurations are [C(3)(R)N(R)]/[C(3)(s)N(s)] and [C(3)(R)N(s)]/[C(3)(s)N(R)], respectively. The computed carbonyl frequencies for g(1) (diast(1)) and g(2) (diast(2)) stereoisomers of3 match well the experimental values. The NBO analysis, for 3 shows the important role of the orbital interactions in conformer stabilization and the overall balance of these interactions corroborates that the g, conformer is more stable than the 92 one. The observed abnormal solvent effect on the relative intensities of the carbonyl doublet components is attributed to the molecular crowding in the g2 conformer which hinders its solvation in comparison to the g, conformer (diast(1)). X-ray single crystal analysis performed for 3 shows the existence Of two 92, and g(1b) conformers of diastereomers (diast2, and diast(1b)) whose absolute configurations are [C(3)(R)N(s)] and [C(3)(R)N(R)], respectively. The larger population and. thus, the larger stabilization of the g(2), conformer over the gib form in the crystals may be associated with a larger energy gain deriving from dipole moment coupling in the former conformer along with a series of C-H center dot center dot center dot O electrostatic and hydrogen bond interactions, (C) 2009 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.

Low-lying singlet and triplet electronic states of RhB

The low-lying X-1 Sigma(+), a(3)Delta, A(1)Delta, b(3)Sigma(+), B-1 Pi, c(3)Pi, C-1 Phi, D-1 Sigma(+), E-1 Pi, d(1)Phi, and e(3)Pi electronic states of RhB have been investigated at the ab initio level, using the multistate multiconfigurational second-order perturbation (MS-CASPT2) theory, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the eleven electronic states included in this work, only three (the X-1 Sigma(+), D-1 Sigma(+), and E-1 Pi states) have been investigated experimentally. Potential energy curves, spectroscopic constants, dipole moments, binding energies, and chemical bonding aspects are presented for all electronic states.

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

A theoretical study on the XeF(2) molecule

A relativistic four-component study was performed for the XeF(2) molecule by using the Dirac-Coulomb (DC) Hamiltonian and the relativistic adapted Gaussian basis sets (RAGBSs). The comparison of bond lengths obtained showed that relativistic effects on this property are small (increase of only 0.01 angstrom) while the contribution of electron correlation, obtained at CCSD(T) or CCSD-T levels, is more important (increase of 0.05 angstrom). Electron correlation is also dominant over relativistic effects for dissociation energies. Moreover, the correlation-relativity interaction is shown to be negligible for these properties. The electron affinity, the first ionization potential and the double ionization potential are obtained by means of the Fock-space coupled cluster (FSCC) method, resulting in DC-CCSD-T values of 0.3 eV, 12.5 eV and 32.3 eV, respectively. Vibrational frequencies and some anharmonicity constants were also calculated under the four-component formalism by means of standard perturbation equations. All these molecular properties are, in general, ill satisfactory agreement with available experimental results. Finally, a partition in terms of charge-charge flux-dipole flux (CCFDF) contributions derived by means of the quantum theory of atoms in molecules (QTAIM) in non-relativistic QCISD(FC)/3-21G* calculations was carried out for XeF(2) and KrF(2). This analysis showed that the most remarkable difference between both molecules lies on the charge flux contribution to the asymmetric stretching mode, which is negligible in KrF(2) but important in XeF(2). (c) 2008 Elsevier B.V. All rights reserved.

CASPT2 Study of the Potential Energy Surface of the HSO(2) System

The importance of the HSO(2) system in atmospheric and combustion chemistry has motivated several works dedicated to the study of associated structures and chemical reactions. Nevertheless controversy still exists in connection with the reaction SH + O(2) -> H + SO(2) and also related to the role of the HSOO isomers in the potential energy surface (PES). Here we report high-level ab initio calculation for the electronic ground state of the HSO(2) system. Energetic, geometric, and frequency properties for the major stationary states of the PES are reported at the same level of calculations:,CASPT2/aug-cc-pV(T+d)Z. This study introduces three new stationary points (two saddle points and one minimum). These structures allow the connection of the skewed HSOOs and the HSO(2) minima defining new reaction paths for SH + O(2) -> H + SO(2) and SH + O(2) -> OH + SO. In addition, the location of the HSOO isomers in the reaction pathways have been clarified.

Characterization and spectroscopic analysis of phenol-ethanol hydrogen bonded clusters

Different hydrogen bonded clusters involving phenol and ethanol are studied theoretically using MP2/aug-cc-pVDZ. Nine different 1: 1 clusters are obtained and analyzed according to their stability and spectroscopic properties. Different isomeric forms of ethanol are considered. Attention is also devoted to the spectral shift of the characteristic pi -> pi* transition of phenol. Using TDHF, CIS, CIS(D) and TDB3LYP in aug-cc-pVDZ basis set, all results agree that a red shift is obtained when phenol is the hydrogen donor and a blue shift is obtained in the opposite case. These results are used to rationalize the red shift observed for phenol in liquid ethanol. (C) 2010 Elsevier B.V. All rights reserved.

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

NMR Chemical Shielding and Spin-Spin Coupling Constants of Liquid NH(3): A Systematic Investigation using the Sequential QM/MM Method

The NMR spin coupling parameters, (1)J(N,H) and (2)J(H,H), and the chemical shielding, sigma((15)N), of liquid ammonia are studied from a combined and sequential QM/MM methodology. Monte Carlo simulations are performed to generate statistically uncorrelated configurations that are submitted to density functional theory calculations. Two different Lennard-Jones potentials are used in the liquid simulations. Electronic polarization is included in these two potentials via an iterative procedure with and without geometry relaxation, and the influence on the calculated properties are analyzed. B3LYP/aug-cc-pVTZ-J calculations were used to compute the V(N,H) constants in the interval of -67.8 to -63.9 Hz, depending on the theoretical model used. These can be compared with the experimental results of -61.6 Hz. For the (2)J(H,H) coupling the theoretical results vary between -10.6 to -13.01 Hz. The indirect experimental result derived from partially deuterated liquid is -11.1 Hz. Inclusion of explicit hydrogen bonded molecules gives a small but important contribution. The vapor-to-liquid shifts are also considered. This shift is calculated to be negligible for (1)J(N,H) in agreement with experiment. This is rationalized as a cancellation of the geometry relaxation and pure solvent effects. For the chemical shielding, U(15 N) Calculations at the B3LYP/aug-pcS-3 show that the vapor-to-liquid chemical shift requires the explicit use of solvent molecules. Considering only one ammonia molecule in an electrostatic embedding gives a wrong sign for the chemical shift that is corrected only with the use of explicit additional molecules. The best result calculated for the vapor to liquid chemical shift Delta sigma((15)N) is -25.2 ppm, in good agreement with the experimental value of -22.6 ppm.

Conformational behavior of different possible ways of oligoglycine formation in a solvent-free environment

The possible ways for glycine oligopeptide formation in gas phase, both in the extended P-strand like conformation and folded 2(7)-ribbon like conformations are analyzed using quantum chemical calculations. We focus on the sequential formation of peptide bond through upgradation of the immediate lower order molecule and observe the consequences in other related processes like oligoglycine formation through simultaneous peptide linkage of n glycine monomers and interchange of molecular conformation through peptide linkage. A comparison is made between the structures and binding energies obtained for both conformers. All binding energies are increased by the zero-point energy contribution. The role of electron correlation effects is briefly analyzed. The folded 2(7)-ribbon-like conformations in vacuo are found to be more stable in comparison to the extended structure. (c) 2007 Elsevier B.V. All rights reserved.