Combining genetic algorithm and simulated annealing: A molecular geometry optimization study

We introduce a new hybrid approach to determine the ground state geometry of molecular systems. Firstly, we compared the ability of genetic algorithm (GA) and simulated annealing (SA) to find the lowest energy geometry of silicon clusters with six and 10 atoms. This comparison showed that GA exhibits fast initial convergence, but its performance deteriorates as it approaches the desired global extreme. Interestingly, SA showed a complementary convergence pattern, in addition to high accuracy. Our new procedure combines selected features from GA and SA to achieve weak dependence on initial parameters, parallel search strategy, fast convergence and high accuracy. This hybrid algorithm outperforms GA and SA by one order of magnitude for small silicon clusters (Si6 and Si10). Next, we applied the hybrid method to study the geometry of a 20-atom silicon cluster. It was able to find an original geometry, apparently lower in energy than those previously described in literature. In principle, our procedure can be applied successfully to any molecular system. © 1998 Elsevier Science B.V.

Experimental and theoretical studies on molecular geometry and polymorphism in crystalline drug substances

DUE TO COPYRIGHT RESTRICTIONS ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY AND INFORMATION SERVICES WITH PRIOR ARRANGEMENT

Uso de modelagem molecular no estudo dos conceitos de nucleofilicidade e basicidade

Molecular modeling enables the students to visualize the abstract relationships underlying theoretical concepts that explain experimental data on the molecular and atomic levels. With this aim we used the free software "Arguslab 4.0.1" (semi-empirical method) to study the reaction of 1-chloropropane with ethoxide in solution, known to lead to methyl propyl ether, through the S N2 mechanism, and propene, through the E2 mechanism. This tool allows users to calculate some properties (i. e. heat formation or electric charges) and to produce 3D images (molecular geometry, electrostatic potential surface, etc.) that render the comprehension of the factors underlying the reaction's progress, which are related to the structure of the reagents, and the process kinetic clearer and easier to understand by the students

1-(3-cyanophenyl)-3-(2-furoyl)thiourea

The title compound, C13H9N3O2S, was synthesized from furoyl isothiocyanate and 3-aminobenzonitrile in dry acetone. The thiourea group is in the thioamide form. The thiourea fragment makes dihedral angles of 3.91 (16) and 37.83 (12)degrees with the ketofuran group and the benzene ring, respectively. The molecular geometry is stabilized by N-H center dot center dot center dot O hydrogen bonds. In the crystal structure, centrosymmetrically related molecules are linked by two intermolecular N-H center dot center dot center dot S hydrogen bonds to form dimers.

Synthesis, characterization and reactivity of new cobalt, palladium and nickel complexes bearing α-diimines and P, O ligands

This work describes the synthesis and characterization of a series of new α-diimine and P,O, β-keto and acetamide phosphines ligands, and their complexation to Ni(II), Co(II),Co(III) and Pd(II) to obtain a series of new compounds aiming to study their structural characteristics and to test their catalytic activity. All the compounds synthesized were characterized by the usual spectroscopic and spectrometric techniques: Elemental Analysis, MALDI-TOF-MS spectrometry, IR, UV-vis, 1H, 13C and 31P NMR spectroscopies. Some of the paramagnetic compounds were also characterized by EPR. For the majority of the compounds it was possible to solve their solid state structure by single crystal X-ray diffraction. Tests for olefin polymerization were performed in order to determine the catalytic activity of the Co(II) complexes. Chapter I presents a brief introduction to homogenous catalysis, highlighting the reactions catalyzed by the type of compounds described in this thesis, namely olefin polymerization and oligomerization and reactions catalyzed by the complexes bearing α-diimines and P,O type ligands. Chapter II is dedicated to the description of the synthesis of new α-diimines cobalt (II) complexes, of general formula [CoX2(α-diimine)], where X = Cl or I and the α-diimines are bis(aryl)acenaphthenequinonediimine) (Ar-BIAN) and 1,4-diaryl-2,3-dimethyl-1,4-diaza-1,3-butadiene (Ar-DAB). Structures solved by single crystal X-ray diffraction were obtained for all the described complexes. For some of the compounds, X-band EPR measurements were performed on polycrystalline samples, showing a high-spin Co(II) (S = 3/2) ion, in a distorted axial environment. EPR single crystal experiments on two of the compounds allowed us to determine the g tensor orientation in the molecular structure. In Chapter III we continue with the synthesis and characterization of more cobalt (II)complexes bearing α-diimines of general formula [CoX2(α-diimine)], with X = Cl or I and α-diimines are bis(aryl)acenaphthenequinonediimine) (Ar-BIAN) and 1,4-diaryl-2,3-dimethyl- 1,4-diaza-1,3-butadiene (Ar-DAB). The structures of three of the new compounds synthesized were determined by single crystal X-ray diffraction. A NMR paramagnetic characterization of all the compounds described is presented. Ethylene polymerization tests were done to determine the catalytic activity of several of the Co(II) complexes described in Chapter II and III and their results are shown. In Chapter IV a new rigid bidentate ligand, bis(1-naphthylimino)acenaphthene, and its complexes with Zn(II) and Pd(II), were synthesized. Both the ligand and its complexes show syn and anti isomers. Structures of the ligand and the anti isomer of the Pd(II) complex were solved by single crystal X-ray diffraction. All the compounds were characterized by elemental analysis, MALDI-TOF-MS spectrometry, and by IR, UV-vis, 1H, 13C, 1H-1H COSY, 1H-13C HSQC, 1H-13C HSQC-TOCSY and 1H-1H NOESY NMR when necessary. DFT studies showed that both conformers of [PdCl2(BIAN)] are isoenergetics and can be obtain experimentally. However, we can predict that the isomerization process is not available in square-planar complex, but is possible for the free ligand. The molecular geometry is very similar in both isomers, and only different orientations for naphthyl groups can be expected. Chapter V describes the synthesis of new P, O type ligands, β-keto phosphine, R2PCH2C(O)Ph, and acetamide phosphine R2PNHC(O)Me, as well as a series of new cobalt(III) complexes namely [(η5-C5H5)CoI2{Ph2PCH2C(O)Ph}], and [(η5- C5H5)CoI2{Ph2PNHC(O)Me}]. Treating these Co(III) compounds with an excess of Et3N, resulted in complexes η2-phosphinoenolate [(η5-C5H5)CoI{Ph2PCH…C(…O)Ph}] and η2- acetamide phosphine [(η5-C5H5)CoI{Ph2PN…C(…O)Me}]. Nickel (II) complexes were also obtained: cis-[Ni(Ph2PN…C(…O)Me)2] and cis-[Ni((i-Pr)2PN…C(…O)Me)2]. Their geometry and isomerism were discussed. Seven structures of the compounds described in this chapter were determined by single crystal X-ray diffraction. The general conclusions of this work can be found in Chapter VI.

Estudo da geometria da uréia por métodos ab initio e simulação computacional de líquidos

A study was carried out on the urea geometries using ab initio calculation and Monte Carlo computational simulation of liquids. The ab initio calculated results showed that urea has a non-planar conformation in the gas phase in which the hydrogen atoms are out of the plane formed by the heavy atoms. Free energies associated to the rotation of the amino groups of urea in water were obtained using the Monte Carlo method in which the thermodynamic perturbation theory is implemented. The magnitude of the free energy obtained from this simulation did not permit us to conclude that urea is non-planar in water.

A Theoretical investigation of a potential high energy density compound 3,6,7,8-tetranitro-3,6,7,8-tetraaza-tricyclo[3.1.1.1(2,4)]octane

The B3LYP/6-31G (d) density functional theory (DFT) method was used to study molecular geometry, electronic structure, infrared spectrum (IR) and thermodynamic properties. Heat of formation (HOF) and calculated density were estimated to evaluate detonation properties using Kamlet-Jacobs equations. Thermal stability of 3,6,7,8-tetranitro-3,6,7,8-tetraaza-tricyclo [3.1.1.1(2,4)]octane (TTTO) was investigated by calculating bond dissociation energy (BDE) at the unrestricted B3LYP/6-31G(d) level. Results showed the N-NO2 bond is a trigger bond during the thermolysis initiation process. The crystal structure obtained by molecular mechanics (MM) methods belongs to P2(1)/C space group, with cell parameters a = 8.239 Å, b = 8.079 Å, c = 16.860 Å, Z = 4 and r = 1.922 g cm-3. Both detonation velocity of 9.79 km s-1 and detonation pressure of 44.22 GPa performed similarly to CL-20. According to the quantitative standards of energetics and stability, TTTO essentially satisfies this requirement as a high energy density compound (HEDC).

Synthesis and characterization of solid 2-methoxycinnamylidenepyruvic acid

The 2-methoxycinnamylidenepyruvic acid (2-MeO-HCP) was synthesized and characterized for nuclear magnetic resonance (¹H and 13C NMR), mass spectrometry (MS), Infrared spectroscopy (FTIR) and differential scanning calorimetry (DSC). The application of DSC for purity determination is well documented in literature and is used in the analysis of pure organic compounds. The molecular geometry and vibrational frequencies of 2-MeO-HCP have been calculated.

Organisation moléculaire dirigée par le groupe CONH2 en 2D et 3D

Notre étude a pour objet la conception, la synthèse ainsi que l’étude structurale d’architectures supramoléculaires obtenues par auto-assemblage, en se basant sur les concepts de la tectonique moléculaire. Cette branche de la chimie supramoléculaire s’occupe de la conception et la synthèse de molécules organiques appelées tectons, du grec tectos qui signifie constructeur. Le tecton est souvent constitué de sites de reconnaissance branchés sur un squelette bien choisi. Les sites de reconnaissance orientés par la géométrie du squelette peuvent participer dans des interactions intermoléculaires qui sont suffisamment fortes et directionnelles pour guider la topologie du cristal résultant. La stratégie envisagée utilise des processus d'auto-assemblage engageant des interactions réversibles entre les tectons. L’auto-assemblage dirigé par de fortes interactions intermoléculaires directionnelles est largement utilisé pour fabriquer des matériaux dont les composants doivent être positionnés en trois dimensions (3D) d'une manière prévisible. Cette stratégie peut également être utilisée pour contrôler l’association moléculaire en deux dimensions (2D), ce qui permet la construction de monocouches organisées et prédéterminées sur différents types des surfaces, tels que le graphite.Notre travail a mis l’accent sur le comportement de la fonction amide comme fonction de reconnaissance qui est un analogue du groupement carboxyle déjà utilisé dans plusieurs études précédentes. Nous avons étudié le comportement d’une série de composés contenant un noyau plat conçu pour faciliter l'adsorption sur le graphite et modifiés par l'ajout de groupes amide pour favoriser la formation de liaisons hydrogène entre les molécules ainsi adsorbées. La capacité de ces composés à former de monocouches organisées à l’échelle moléculaire en 2D a été examinée par microscopie à effet tunnel, etleur organisation en 3D a également été étudiée par cristallographie aux rayons X. Dans notre étude, nous avons systématiquement modifié la géométrie moléculaire et d'autres paramètres afin d'examiner leurs effets sur l'organisation moléculaire. Nos résultats suggèrent que les analyses structurales combinées en 2D et 3D constituent un important atout dans l'effort pour comprendre les interactions entre les molécules adsorbées et l’effet de l’interaction avec la surface du substrat.

Implementation and application of basis set superposition error-correction schemes to the theoretical modeling of weak intermolecular interactions

This thesis deals with the so-called Basis Set Superposition Error (BSSE) from both a methodological and a practical point of view. The purpose of the present thesis is twofold: (a) to contribute step ahead in the correct characterization of weakly bound complexes and, (b) to shed light the understanding of the actual implications of the basis set extension effects in the ab intio calculations and contribute to the BSSE debate. The existing BSSE-correction procedures are deeply analyzed, compared, validated and, if necessary, improved. A new interpretation of the counterpoise (CP) method is used in order to define counterpoise-corrected descriptions of the molecular complexes. This novel point of view allows for a study of the BSSE-effects not only in the interaction energy but also on the potential energy surface and, in general, in any property derived from the molecular energy and its derivatives A program has been developed for the calculation of CP-corrected geometry optimizations and vibrational frequencies, also using several counterpoise schemes for the case of molecular clusters. The method has also been implemented in Gaussian98 revA10 package. The Chemical Hamiltonian Approach (CHA) methodology has been also implemented at the RHF and UHF levels of theory for an arbitrary number interacting systems using an algorithm based on block-diagonal matrices. Along with the methodological development, the effects of the BSSE on the properties of molecular complexes have been discussed in detail. The CP and CHA methodologies are used for the determination of BSSE-corrected molecular complexes properties related to the Potential Energy Surfaces and molecular wavefunction, respectively. First, the behaviour of both BSSE-correction schemes are systematically compared at different levels of theory and basis sets for a number of hydrogen-bonded complexes. The Complete Basis Set (CBS) limit of both uncorrected and CP-corrected molecular properties like stabilization energies and intermolecular distances has also been determined, showing the capital importance of the BSSE correction. Several controversial topics of the BSSE correction are addressed as well. The application of the counterpoise method is applied to internal rotational barriers. The importance of the nuclear relaxation term is also pointed out. The viability of the CP method for dealing with charged complexes and the BSSE effects on the double-well PES blue-shifted hydrogen bonds is also studied in detail. In the case of the molecular clusters the effect of high-order BSSE effects introduced with the hierarchical counterpoise scheme is also determined. The effect of the BSSE on the electron density-related properties is also addressed. The first-order electron density obtained with the CHA/F and CHA/DFT methodologies was used to assess, both graphically and numerically, the redistribution of the charge density upon BSSE-correction. Several tools like the Atoms in Molecules topologycal analysis, density difference maps, Quantum Molecular Similarity, and Chemical Energy Component Analysis were used to deeply analyze, for the first time, the BSSE effects on the electron density of several hydrogen bonded complexes of increasing size. The indirect effect of the BSSE on intermolecular perturbation theory results is also pointed out It is shown that for a BSSE-free SAPT study of hydrogen fluoride clusters, the use of a counterpoise-corrected PES is essential in order to determine the proper molecular geometry to perform the SAPT analysis.

Coriolis perturbations in the infrared spectrum of allene

High resolution infrared spectra of the ν9 and ν10 perpendicular fundamentals of the allene molecule are reported, in which the J structure in the sub-bands has been partially resolved. Analysis of the latter shows that the vibrational origin ν9 = 999 cm−1, some 35 cm−1 below previous assignments. The pronounced asymmetry in the intensity distribution of the rotational structure which this assignment implies is shown to be expected theoretically, due to the Coriolis perturbations involved, and it is interpreted in terms of the sign and magnitude of the ratio of the dipole moment derivatives in the two fundamentals. The results of this analysis are shown to be in good agreement with observations on allene-1.1-d2, where similar intensity perturbations are observed, and with an independent analysis of the ν8 band of allene-h4. The A rotational constant of allene-h4 is found to have the value 4.82 ± 0.01 cm−1, and for the molecular geometry we obtain r(CH) = 1.084 A, r(CC) = 1.308 A, and HCH = 118.4°. A partial analysis of the rotational structure of the hot bands (ν9 + ν11 − ν11) and (ν10 + ν11 − ν11) is presented; these provide an example of a strong Coriolis interaction between nearly degenerate A1A2 and B1B2 pairs of vibrational levels. Some localized rotational perturbations in the ν9 and ν10 fundamentals are also noted, and their possible interpretations are discussed.

The general harmonic force field of methyl fluoride

The complete general harmonic force field of methyl flouride was recalculated using the most recent literature frequency, Coriolis ζ, and centrifugal distortion data for 12CH3F, 13CH3F, 12CD3F, 12CHD2F and 12CH2DF. The anharmonic corrections applied to the observed frequency data and the adopted molecular geometry are considered to be more realistic than those used hitherto. There is excellent overall agreement between the fitted force constants and the highest quality ab initio force field currently available.

Salt and ionic cocrystalline forms of amides: protonation of carbamazepine in aqueous media

The products of reactions of the pharmaceutical amide carbamazepine (CBZ) with strong acids under aqueous conditions were investigated by both powder and single crystal X-ray diffraction. Despite previous claims to the contrary, it was found that salt forms with CBZ protonated at the amide O atom could be isolated from reactions with both HCl and HBr. These forms include the newly identified hydrate phase [CBZ(H)][Cl]·H O. Reactions with other mineral acids (HI and HBF ) gave ionic cocrystalline (ICC) forms (CBZ· [acridinium][I ]·2.5I and CBZ·[H O ] [BF ] ·H O) as well as the salt form CBZ·[CBZ(H)][BF ]·0.5H O. Reaction 2 4 3 2 5 2 0.25 4 0.25 2 4 2 of CBZ with a series of sulfonic acids also gave salt forms, namely, [CBZ(H)][O SC H ], [CBZ(H)][O SC H (OH)]· 3 6 5 3 6 4 0.5H O, [CBZ(H)] [O SCH CH SO ], and [CBZ(H)][O SC H (OH) (COOH)]·H O. CBZ and protonated CBZ(H) 2 2 3 2 2 3 3 6 3 2 moieties can be differentiated in the solid state both by changes to molecular geometry and by differing packing preferences

Crystal structure of garciniaphenone and evidences on the relationship between keto-enol tautomerism and configuration

Garciniaphenone (=rel-(1R,5R,7R)-3-benzoyl-4-hydroxy-8,8-dimethyl-1,7-bis(3-methylbut-2-en-1-yl)bicyclo[3.3.1]non-3-ene-2,9-dione; 1). a novel natural product, was isolated from a hexane extract of Garcinia brasiliensis fruits. The crystal structure of 1 as well as the selected geometrical and Configurational features were compared with those of known related polyprenylated benzophenones. Garciniaphenone is the first representative of polyprenylated benzophenones without a prenyl substituent at C(5). Notably, the absence of a 5-prenyl substituent has an impact on the molecular geometry. The tautomeric form of 1 in the solid state was readily established by a residual-electronic-density map generated by means of a difference Fourier analysis, and there is an entirely delocalized six-membered chelate ring encompassing the keto-enol moiety. The configuration at C(7) was used to rationalize the nature of the keto-enol tautomeric form within 1. The intermolecular array in the network is maintained by nonclassical intermolecular H-bonds.

Crystalline structure of mangiferin, a C-glycosyl-substituted 9H-xanthen-9-one isolated from the stem bark of Mangifera indica

The crystalline structure of mangiferin (= 2-beta-D-glucopyranosyl-1,3,6,7-tetrahydroxy-9H-xanthen-9-one; 1), a biologically active xanthenone C-glycoside, isolated from the stem bark of Mangifera indica (Anacardiaceae), was unambiguously determined by single-crystal X-ray diffraction (XRD). The crystal structure is summarized as follows: triclinic, P1, a = 7.6575(5), b = 11.2094(8), c = 11.8749(8) angstrom, alpha = 79.967(5), beta = 87.988(4), gamma = 72.164(4)degrees, V = 955.3(1) angstrom(3), and Z = 2. The structure also shows two molecules in the asymmetric unit cell and five crystallization H2O molecules. The packing is stabilized by several intermolecular H-bonds involving either the two symmetry-independent mangiferin molecules 1a and 1b, or the H2O ones.