The search for a potentially new molecular species, SiAs: A theoretical contribution

Electronic states of a new molecular species, SiAs, correlating with the three lowest dissociation channels are characterized at a high-level of theory using the CASSCF/MRCI approach along with quintuple-xi quality basis sets. This characterization includes potential energy curves, vibrational energy levels, spectroscopic parameters, dipole and transition dipole moment functions, transition probabilities, and radiative lifetimes. For the ground state (X-2 Pi), an assessment of spin-orbit effects and the interaction with the close-lying A(2)Sigma(+) state is also reported. Similarities and differences with other isovalent species such as SiP and CAs are also discussed.

On the Mechanisms of Triplet Excited State Population in 8-Azaadenine

The photophysics of 8-azaadenine (8-AA) has been studied with the CASPT2//CASSCF protocol and ANO-L double-zeta basis sets. Stationary equilibrium structures, surface crossings, minimum energy paths, and linear interpolations have been used to study possible mechanisms to populate the lowest triplet state, T-1 (3)(pi pi*), capable of sensitizing molecular oxygen. Our results show that two main mechanisms can occur after photoexcitation to the S-2 (1)(pi pi*) state. The first one is through the S-2/S-1 conical intersection (((1)pi pi*/(1)n pi*)(Cl)), leading to the S-1 ((1)n pi*) state minimum, (S-1 ((1)n pi*))(min), where a singlet-triplet crossing, ((1)n pi*/(3)pi pi*)(STC), is accessible. The second one starts with the ((1)pi pi*/(3)n pi*)(STC) at the (S-2((1)pi pi*))(min), from which the system can evolve to the (T-2 ((3)n pi*))(min), with subsequent population of the T-1 excited electronic state, due to the ((3)n pi*/(3)pi pi*)(Cl) conical intersection.

A theoretical contribution characterizing a potentially new molecular species: MgAs

A new molecular species, MgAs, is investigated theoretically for the first time at the CASSCF/MRCI level using quintuple-zeta quality basis sets. Potential energy curves for the lowest-lying electronic states are presented as well as the associated spectroscopic constants. Dipole and transition moment functions for selected states complement this characterization. Estimates of transition probabilities and radiative life-times for the most important transitions are also reported. The effect of spin-orbit interactions is clearly reflected on the potential energy curves. Comparisons with BeAs, BeN, and BeP are made where pertinent. (C) 2011 Elsevier B.V. All rights reserved.

Theoretical study of the XP3 (X = Al, B, Ga) clusters

The lowest singlet and triplet states of AlP3, GaP3 and BP3 molecules with C-s, C-2v and C-3v symmetries were characterized using the B3LYP functional and the aug-cc-pVTZ and aug-cc-pVQZ correlated consistent basis sets. Geometrical parameters and vibrational frequencies were calculated and compared to existent experimental and theoretical data. Relative energies were obtained with single point CCSD(T) calculations using the aug-cc-pVTZ, aug-cc-pVQZ and aug-cc-pV5Z basis sets, and then extrapolating to the complete basis set (CBS) limit. (C) 2011 Elsevier B.V. All rights reserved.

Theoretical study of the absorption and nonradiative deactivation of 1-nitronaphthalene in the low-lying singlet and triplet excited states including methanol and ethanol solvent effects

The photophysics of the 1-nitronaphthalene molecular system, after the absorption transition to the first singlet excited state, is theoretically studied for investigating the ultrafast multiplicity change to the triplet manifold. The consecutive transient absorption spectra experimentally observed in this molecular system are also studied. To identify the electronic states involved in the nonradiative decay, the minimum energy path of the first singlet excited state is obtained using the complete active space self-consistent field//configurational second-order perturbation approach. A near degeneracy region was found between the first singlet and the second triplet excited states with large spin-orbit coupling between them. The intersystem crossing rate was also evaluated. To support the proposed deactivation model the transient absorption spectra observed in the experiments were also considered. For this, computer simulations using sequential quantum mechanic-molecular mechanic methodology was used to consider the solvent effect in the ground and excited states for proper comparison with the experimental results. The absorption transitions from the second triplet excited state in the relaxed geometry permit to describe the transient absorption band experimentally observed around 200 fs after the absorption transition. This indicates that the T-2 electronic state is populated through the intersystem crossing presented here. The two transient absorption bands experimentally observed between 2 and 45 ps after the absorption transition are described here as the T-1 -> T-3 and T-1 -> T-5 transitions, supporting that the intermediate triplet state (T-2) decays by internal conversion to T-1. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4738757]

The radical SeCl: A theoretical contribution to the characterization of its low-lying electronic states

All doublet and quartet electronic states correlating with the first dissociation channel of SeCl and some Rydberg states are investigated theoretically at the CASSCF/MRCI level of theory using extended basis sets, including the contribution of spin-orbit effects. The similarity of the potential energy curves with those of SeF suggests that spectroscopic constants for the ground (X (2)Pi) and the first excited quartet (a(4)Sigma) of SeCl could also be determined via an emission resulting from the reaction of selenium with atomic chlorine. The coupling constant of the ground state at R-e is estimated as -1610 cm (1). The potential energy curves calculated and the derived spectroscopic constants do not support the interpretation and assignment of the scarce transitions recorded experimentally as due to (2)Pi-(2)Pi emissions. That the few observed lines might arise from transitions from the state b(4)Sigma(-)(1/2) to a very high vibrational level of the state a(4)Sigma(-)(1/2) is an open possibility, however, the number of vibrational states and the calculated Delta G(1/2) differ significantly from the reported ones. (C) 2012 Elsevier B. V. All rights reserved.

Ionization of chlorophyll-c(2) in liquid methanol

The ionization of chlorophyll-c(2) in liquid methanol was investigated by a sequential quantum mechanical/Monte Carlo approach. Focus was placed on the determination of the first ionization energy of chlorophyll-c(2). The results show that the first vertical ionization energy (IE) is red-shifted by 0.47 +/- 0.24 eV relative to the gas-phase value. The red-shift of the chlorophyll-c(2) IE in the liquid phase can be explained by Mg center dot center dot center dot OH hydrogen bonding and long-ranged electrostatic interactions in solution. The ionization threshold for chlorophyll-c2 in liquid methanol is close to 6 eV. (C) 2012 Elsevier B.V. All rights reserved.

Cyclic Voltammetry and Computational Chemistry Studies on the Evaluation of the Redox Behavior of Parabens and other Analogues

Parabens are antimicrobial preservatives widely used in pharmaceutical, cosmetic and food industries. The alkyl chain connected to the ester group defines some important physicochemical characteristics of these compounds, including the partition coefficient and redox properties. The voltammetric and computational analyses were carried out in order to evaluate the redox behavior of these compounds and other phenolic analogues. A strong correlation between chemical substituents inductive effects of parabens with redox potentials was observed. Using cyclic voltammetry and glassy carbon working electrode, only one irreversible anodic peak was observed around 0.8 V for methylparaben (MP), ethylparaben (EP), propylparaben (PP), butylparaben (BP), benzylparaben (BzP) and p-substituted phenolic analogues. The electrodonating inductive effect of alkyl groups was demonstrated by the anodic oxidation potential shift to lower values as the carbon number increases and, therefore the parabens (and other phenolic analogues) oxidation processes to the quinonoidic forms showed great dependence on the substituent pattern.

Coupled-Cluster-Berechnungen von Parametern der Kernspin-Resonanz-Spektroskopie

Coupled-Cluster-Berechnungen von Parametern derKernspin-Resonanz-Spektroskopie Dissertationsschrift von Alexander A.Auer, Mainz 2002 Im Rahmen einer Studie der Berechnung von 13C-Verschiebungenwerdendie Einfluesse von Elektronenkorrelation, Basissatz,Gleichgewichtsgeometrie sowie Schwingungs- und Rotationseffekten separat betrachtet.Dabei zeigt sich, dass dieCoupled-Cluster-Singles-Doubles-Methode mitstoerungstheoretischer Behandlung der Dreifachanregungen(CCSD(T)) mit entsprechend grossen Basissaetzen bei Beruecksichtigung derNullpunktsschwingungseffekte Ergebnisse mit ca. 1 ppm Abweichung zum Experiment liefert. Eine Analyse der Elektronenkorrelationseffekte beiCoupled-Cluster- (CC-) Berechnungen von indirekten Spin-Spin-Kopplungskonstanten zeigt, dassCC-Methoden mit Hartree-Fock-Orbitalrelaxation zur Berechnung derKopplungskonstanten ungeeignet sind. Eine Loesung ist die Verwendung unrelaxierter CC-Methoden,in denendie HF-Orbitalrelaxation aus der Berechnung der gestoertenWellenfunktion ausgeschlossen wird. Full-Configuration-Interaction-Berechnungen fuer Borhydridzeigen,dass auf CC-Singles-Doubles-Niveau (CCSD) 94% und aufCC-Singles-Doubles-Triples-Niveau (CCSDT) 99% der Korrelationseffekte beschrieben werden. Weiterhin istdie Beruecksichtigung der Nullpunktsschwingung sowie die Wahl eines ausreichend grossen Basissatzes wichtig. Auf Grundlage der vorangegangenen Studien werden im letztenTeil zwei Beispiele zur Anwendung hochgenauer Berechnungen vonNMR-Parametern vorgestellt.Im Rahmen einer Studie der Spin-Spin-Kopplungskonstanten vonCyclopentan wird eine Karplus-Beziehungzwischen den Kopplungskonstanten und der Konformation desMolekuels aufgestellt, desweiteren werden die NMR-Parameter von Methylidinphosphanuntersucht.

Hochgenaue Berechnungen von Energien und molekularen Eigenschaften mit Hilfe der Coupled-Cluster-Theorie

Coupled-Cluster-Theorie (CC) ist in der heutigen Quantenchemie eine der erfolgreichsten Methoden zur genauen Beschreibung von Molekülen. Die in dieser Arbeit vorgestellten Ergebnisse zeigen, daß neben den Berechnungen von Energien eine Reihe von Eigenschaften wie Strukturparameter, Schwingungsfrequenzen und Rotations-Schwingungs-Parameter kleiner und mittelgrofler Moleküle zuverlässig und präzise vorhergesagt werden können. Im ersten Teil der Arbeit wird mit dem Spin-adaptierten Coupled-Cluster-Ansatz (SA-CC) ein neuer Weg zur Verbesserung der Beschreibung von offenschaligen Systemen vorgestellt. Dabei werden zur Bestimmung der unbekannten Wellenfunktionsparameter zusätzlich die CC-Spingleichungen gelöst. Durch dieses Vorgehen wird gewährleistet, daß die erhaltene Wellenfunktion eine Spineigenfunktion ist. Die durchgeführte Implementierung des Spin-adaptierten CC-Ansatzes unter Berücksichtigung von Einfach- und Zweifachanregungen (CCSD) für high-spin Triplett-Systeme wird ausführlich erläutert. Im zweiten Teil werden CC-Additionsschemata vorgestellt, die auf der Annahme der Additivität von Elektronenkorrelations- und Basissatzeffekten basieren. Die etablierte Vorgehensweise, verschiedene Beiträge zur Energie mit an den Rechenaufwand angepassten Basissätzen separat zu berechnen und aufzusummieren, wird hier auf Gradienten und Kraftkonstanten übertragen. Für eine Beschreibung von Bindungslängen und harmonischen Schwingungsfrequenzen mit experimenteller Genauigkeit ist die Berücksichtigung von Innerschalenkorrelationseffekten sowie Dreifach- und Vierfachanregungen im Clusteroperator der Wellenfunktion nötig. Die Basissatzkonvergenz wird dabei zusätzlich mit Extrapolationsmethoden beschleunigt. Die quantitative Vorhersage der Bindungslängen von 17 kleinen Molekülen, aufgebaut aus Atomen der ersten Langperiode, ist so mit einer Genauigkeit von wenigen Hundertstel Pikometern möglich. Für die Schwingungsfrequenzen dieser Moleküle weist das CC-Additionsschema basierend auf den berechneten Kraftkonstanten im Vergleich zu experimentellen Ergebnissen einen mittleren absoluten Fehler von 3.5 cm-1 und eine Standardabweichung von 2.2 cm-1 auf. Darüber hinaus werden zur Unterstützung von experimentellen Untersuchungen berechnete spektroskopische Daten einiger größerer Moleküle vorgelegt. Die in dieser Arbeit vorgestellten Untersuchungen zur Isomerisierung von Dihalogensulfanen XSSX (X= F, Cl) oder die Berechnung von Struktur- und Rotations-Schwingungs-Parametern für die Moleküle CHCl2F und CHClF2 zeigen, daß bereits störungstheoretische CCSD(T)-Näherungsmethoden qualitativ gute Vorhersagen experimenteller Resultate liefern. Desweiteren werden Diskrepanzen von experimentellen und berechneten Bindungsabständen bei den Molekülen Borhydrid- und Carbenylium durch die Berücksichtigung des elektronischen Beitrages zum Trägheitsmoment beseitigt.

Large-scale coupled-cluster calculations

Coupled-cluster theory provides one of the most successful concepts in electronic-structure theory. This work covers the parallelization of coupled-cluster energies, gradients, and second derivatives and its application to selected large-scale chemical problems, beside the more practical aspects such as the publication and support of the quantum-chemistry package ACES II MAB and the design and development of a computational environment optimized for coupled-cluster calculations. The main objective of this thesis was to extend the range of applicability of coupled-cluster models to larger molecular systems and their properties and therefore to bring large-scale coupled-cluster calculations into day-to-day routine of computational chemistry. A straightforward strategy for the parallelization of CCSD and CCSD(T) energies, gradients, and second derivatives has been outlined and implemented for closed-shell and open-shell references. Starting from the highly efficient serial implementation of the ACES II MAB computer code an adaptation for affordable workstation clusters has been obtained by parallelizing the most time-consuming steps of the algorithms. Benchmark calculations for systems with up to 1300 basis functions and the presented applications show that the resulting algorithm for energies, gradients and second derivatives at the CCSD and CCSD(T) level of theory exhibits good scaling with the number of processors and substantially extends the range of applicability. Within the framework of the ’High accuracy Extrapolated Ab initio Thermochemistry’ (HEAT) protocols effects of increased basis-set size and higher excitations in the coupled- cluster expansion were investigated. The HEAT scheme was generalized for molecules containing second-row atoms in the case of vinyl chloride. This allowed the different experimental reported values to be discriminated. In the case of the benzene molecule it was shown that even for molecules of this size chemical accuracy can be achieved. Near-quantitative agreement with experiment (about 2 ppm deviation) for the prediction of fluorine-19 nuclear magnetic shielding constants can be achieved by employing the CCSD(T) model together with large basis sets at accurate equilibrium geometries if vibrational averaging and temperature corrections via second-order vibrational perturbation theory are considered. Applying a very similar level of theory for the calculation of the carbon-13 NMR chemical shifts of benzene resulted in quantitative agreement with experimental gas-phase data. The NMR chemical shift study for the bridgehead 1-adamantyl cation at the CCSD(T) level resolved earlier discrepancies of lower-level theoretical treatment. The equilibrium structure of diacetylene has been determined based on the combination of experimental rotational constants of thirteen isotopic species and zero-point vibrational corrections calculated at various quantum-chemical levels. These empirical equilibrium structures agree to within 0.1 pm irrespective of the theoretical level employed. High-level quantum-chemical calculations on the hyperfine structure parameters of the cyanopolyynes were found to be in excellent agreement with experiment. Finally, the theoretically most accurate determination of the molecular equilibrium structure of ferrocene to date is presented.

Accurate Predictions of Water Cluster Formation, (H2O)n=2-10

An efficient mixed molecular dynamics/quantum mechanics model has been applied to the water cluster system. The use of the MP2 method and correlation consistent basis sets, with appropriate correction for BSSE, allows for the accurate calculation of electronic and free energies for the formation of clusters of 2−10 water molecules. This approach reveals new low energy conformers for (H2O)n=7,9,10. The water heptamer conformers comprise five different structural motifs ranging from a three-dimensional prism to a quasi-planar book structure. A prism-like structure is favored energetically at low temperatures, but a chair-like structure is the global Gibbs free energy minimum past 200 K. The water nonamers exhibit less complexity with all the low energy structures shaped like a prism. The decamer has 30 conformers that are within 2 kcal/mol of the Gibbs free energy minimum structure at 298 K. These structures are categorized into four conformer classes, and a pentagonal prism is the most stable structure from 0 to 320 K. Results can be used as benchmark values for empirical water models and density functionals, and the method can be applied to larger water clusters.

Semiempirical Study of the Bergman Reaction: Towards a Computationally Efficient and Accurate Method for Modeling Enediyne Anticancer Antibiotics

Over the past 7 years, the enediyne anticancer antibiotics have been widely studied due to their DNA cleaving ability. The focus of these antibiotics, represented by kedarcidin chromophore, neocarzinostatin chromophore, calicheamicin, esperamicin A, and dynemicin A, is on the enediyne moiety contained within each of these antibiotics. In its inactive form, the moiety is benign to its environment. Upon suitable activation, the system undergoes a Bergman cycloaromatization proceeding through a 1,4-dehydrobenzene diradical intermediate. It is this diradical intermediate that is thought to cleave double-stranded dna through hydrogen atom abstraction. Semiempirical, semiempiricalci, Hartree–Fock ab initio, and mp2 electron correlation methods have been used to investigate the inactive hex-3-ene-1,5-diyne reactant, the 1,4-dehydrobenzene diradical, and a transition state structure of the Bergman reaction. Geometries calculated with different basis sets and by semiempirical methods have been used for single-point calculations using electron correlation methods. These results are compared with the best experimental and theoretical results reported in the literature. Implications of these results for computational studies of the enediyne anticancer antibiotics are discussed.

Structures of Cage, Prism, and Book Isomers of Water Hexamer from Broadband Rotational Spectroscopy

Theory predicts the water hexamer to be the smallest water cluster with a three-dimensional hydrogen-bonding network as its minimum energy structure. There are several possible low-energy isomers, and calculations with different methods and basis sets assign them different relative stabilities. Previous experimental work has provided evidence for the cage, book, and cyclic isomers, but no experiment has identified multiple coexisting structures. Here, we report that broadband rotational spectroscopy in a pulsed supersonic expansion unambiguously identifies all three isomers; we determined their oxygen framework structures by means of oxygen-18–substituted water (H218O). Relative isomer populations at different expansion conditions establish that the cage isomer is the minimum energy structure. Rotational spectra consistent with predicted heptamer and nonamer structures have also been identified.

Atmospheric Implications for Formation of Clusters of Ammonium and 1−10 Water Molecules

A mixed molecular dynamics/quantum mechanics model has been applied to the ammonium/water clustering system. The use of the high level MP2 calculation method and correlated basis sets, such as aug-cc-pVDZ and aug-cc-pVTZ, lends confidence in the accuracy of the extrapolated energies. These calculations provide electronic and free energies for the formation of clusters of ammonium and 1−10 water molecules at two different temperatures. Structures and thermodynamic values are in good agreement with previous experimental and theoretical results. The estimated concentration of these clusters in the troposphere was calculated using atmospheric amounts of ammonium and water. Results show the favorability of forming these clusters and implications for ion-induced nucleation in the atmosphere.