Theoretical chemistry: current applications to photochemistry and thermochemistry

A historical perspective is given contrasting challenges and advances in theoretical chemistry at the time the first issue of Theoretical Chemistry Accounts appeared in 1962 and the progress achieved since then as expressed in current state-of-the-art applications in photochemistry and thermochemistry.

Computation of (3)J(HH) coupling constants with a combination of density functional theory and semiempirical calculations. Application to complex molecules

This work evaluates the efficiency of economic levels of theory for the prediction of (3)J(HH) spin-spin coupling constants, to be used when robust electronic structure methods are prohibitive. To that purpose, DFT methods like mPW1PW91. B3LYP and PBEPBE were used to obtain coupling constants for a test set whose coupling constants are well known. Satisfactory results were obtained in most of cases, with the mPW1PW91/6-31G(d,p)//B3LYP/6-31G(d,p) leading the set. In a second step. B3LYP was replaced by the semiempirical methods PM6 and RM1 in the geometry optimizations. Coupling constants calculated with these latter structures were at least as good as the ones obtained by pure DFT methods. This is a promising result, because some of the main objectives of computational chemistry - low computational cost and time, allied to high performance and precision - were attained together. (C) 2012 Elsevier B.V. All rights reserved.

Structure and electronic properties of hydrated mesityl oxide: a sequential quantum mechanics/molecular mechanics approach

The hydration of mesityl oxide (MOx) was investigated through a sequential quantum mechanics/molecular mechanics approach. Emphasis was placed on the analysis of the role played by water in the MOx syn-anti equilibrium and the electronic absorption spectrum. Results for the structure of the MOx-water solution, free energy of solvation and polarization effects are also reported. Our main conclusion was that in gas-phase and in low-polarity solvents, the MOx exists dominantly in syn-form and in aqueous solution in anti-form. This conclusion was supported by Gibbs free energy calculations in gas phase and in-water by quantum mechanical calculations with polarizable continuum model and thermodynamic perturbation theory in Monte Carlo simulations using a polarized MOx model. The consideration of the in-water polarization of the MOx is very important to correctly describe the solute-solvent electrostatic interaction. Our best estimate for the shift of the pi-pi* transition energy of MOx, when it changes from gas-phase to water solvent, shows a red-shift of -2,520 +/- 90 cm(-1), which is only 110 cm(-1) (0.014 eV) below the experimental extrapolation of -2,410 +/- 90 cm(-1). This red-shift of around -2,500 cm(-1) can be divided in two distinct and opposite contributions. One contribution is related to the syn -> anti conformational change leading to a blue-shift of similar to 1,700 cm(-1). Other contribution is the solvent effect on the electronic structure of the MOx leading to a red-shift of around -4,200 cm(-1). Additionally, this red-shift caused by the solvent effect on the electronic structure can by composed by approximately 60 % due to the electrostatic bulk effect, 10 % due to the explicit inclusion of the hydrogen-bonded water molecules and 30 % due to the explicit inclusion of the nearest water molecules.

A simple analysis of the influence of the solvent-induced electronic polarization on the N-15 magnetic shielding of pyridine in water

Electronic polarization induced by the interaction of a reference molecule with a liquid environment is expected to affect the magnetic shielding constants. Understanding this effect using realistic theoretical models is important for proper use of nuclear magnetic resonance in molecular characterization. In this work, we consider the pyridine molecule in water as a model system to briefly investigate this aspect. Thus, Monte Carlo simulations and quantum mechanics calculations based on the B3LYP/6-311++G (d,p) are used to analyze different aspects of the solvent effects on the N-15 magnetic shielding constant of pyridine in water. This includes in special the geometry relaxation and the electronic polarization of the solute by the solvent. The polarization effect is found to be very important, but, as expected for pyridine, the geometry relaxation contribution is essentially negligible. Using an average electrostatic model of the solvent, the magnetic shielding constant is calculated as -58.7 ppm, in good agreement with the experimental value of -56.3 ppm. The explicit inclusion of hydrogen-bonded water molecules embedded in the electrostatic field of the remaining solvent molecules gives the value of -61.8 ppm.

SCIENTIFIC LANGUAGE RHETORIC: FROM THE THEORETICAL BASIS TO THE PRODUCTION OF DIDACTIC MATERIAL FOR UNDERGRADUATE CHEMISTRY TEACHING

This work is part of a study that focused on analyzing the contributions of didactic activities related to scientific language rhetoric characteristics aimed at developing students' abilities to identify such characteristics in chemistry scientific texts and critical reading of those texts. In this study, we present the theoretical basis adopted to determine the scientific discourse characteristics and for the production of the didactic material used in those activities. Latour, Coracini and Campanario studies on persuasive rhetorical strategies present in scientific articles aided the production of such material.

Correlation Between Molecular Conformation, Packing, and Dynamics in Oligofluorenes: A Theoretical/Experimental Study

Fluorene-based systems have shown great potential as components in organic electronics and optoelectronics (organic photovoltaics, OPVs, organic light emitting diodes, OLEDs, and organic transistors, OTFTs). These systems have drawn attention primarily because they exhibit strong blue emission associated with relatively good thermal stability. It is well-known that the electronic properties of polymers are directly related to the molecular conformations and chain packing of polymers. Here, we used three oligofluorenes (trimer, pentamer, and heptamer) as model systems to theoretically investigate the conformational properties of fluorene molecules, starting with the identification of preferred conformations. The hybrid exchange correlation functional, OPBE, and ZINDO/S-CI showed that each oligomer exhibits a tendency to adopt a specific chain arrangement, which could be distinguished by comparing their UV/vis electronic absorption and C-13 NMR spectra. This feature was used to identify the preferred conformation of the oligomer chains in chloroform-cast films by comparing experimental and theoretical UV/vis and C-13 NMR spectra. Moreover, the oligomer chain packing and dynamics in the films were studied by DSC and several solid state NMR techniques, which indicated that the phase behavior of the films may be influenced by the tendency that each oligomeric chain has to adopt a given conformation.

Electronic stopping cross sections for protons in Al2O3: an experimental and theoretical study

The electronic stopping cross section (SCS) of Al2O3 for proton beams is studied both experimentally and theoretically. The measurements are made for proton energies from 40 keV up to 1 MeV, which cover the maximum stopping region, using two experimental methods, the transmission technique at low energies (similar to 40-175 keV) and the Rutherford backscattering at high energies (approximate to 190-1000 keV). These new data reveal an increment of 16% in the SCS around the maximum stopping with respect to older measurements. The theoretical study includes electronic stopping power calculations based on the dielectric formalism and on the transport cross section (TCS) model to describe the electron excitations of Al2O3. The non-linear TCS calculations of the SCS for valence electrons together with the generalized oscillator strengths (GOS) model for the core electrons compare well with the experimental data in the whole range of energies considered.

Theoretical models for the antitrypanosomal activity of thiosemicarbazone derivatives

Thiosemicarbazones are cruzain inhibitors which have been identified as potential antitrypanosomal agents. In this work, several molecular properties were calculated at the density functional theory (DFT)/B3LYP/6-311G* level for a set of 44 thiosemicarbazones. Unsupervised and supervised pattern recognition techniques (hierarchical cluster analysis, principal component analysis, kth-nearest neighbors, and soft independent modeling by class analogy) were used to obtain structureactivity relationship models, which are able to classify unknown compounds according to their activities. The chemometric analyses performed here revealed that 12 descriptors can be considered responsible for the discrimination between high and low activity compounds. Classification models were validated with an external test set, showing that predictive classifications were achieved with the selected variable set. The results obtained here are in good agreement with previous findings from the literature, suggesting that our models can be useful on further investigations on the molecular determinants for the antichagasic activity. (C) 2012 Wiley Periodicals, Inc.

Molecular chemistry and the missing mass problem in planetary nebulae

Context. Detections of molecular lines, mainly from H-2 and CO, reveal molecular material in planetary nebulae. Observations of a variety of molecules suggest that the molecular composition in these objects differs from that found in interstellar clouds or in circumstellar envelopes. The success of the models, which are mostly devoted to explain molecular densities in specific planetary nebulae, is still partial however. Aims. The present study aims at identifying the influence of stellar and nebular properties on the molecular composition of planetary nebulae by means of chemical models. A comparison of theoretical results with those derived from the observations may provide clues to the conditions that favor the presence of a particular molecule. Methods. A self-consistent photoionization numerical code was adapted to simulate cold molecular regions beyond the ionized zone. The code was used to obtain a grid of models and the resulting column densities are compared with those inferred from observations. Results. Our models show that the inclusion of an incident flux of X-rays is required to explain the molecular composition derived for planetary nebulae. We also obtain a more accurate relation for the N(CO)/N(H-2) ratio in these objects. Molecular masses obtained by previous works in the literature were then recalculated, showing that these masses can be underestimated by up to three orders of magnitude. We conclude that the problem of the missing mass in planetary nebulae can be solved by a more accurate calculation of the molecular mass.

Interfacial concentrations of chloride and bromide in zwitterionic micelles with opposite dipoles: Experimental determination by chemical trapping and a theoretical description

Interfacial concentrations of chloride and bromide ions, with Li+, Na+, K+, Rb+, Cs+, trimethylammonium (TMA(+)), Ca2+, and Mg2+ as counterions, were determined by chemical trapping in micelles formed by two zwitterionic surfactants, namely N-hexadecyl-N,N-dimethyl-3-ammonio-1-propanesulfonate (HPS) and hexadecylphosphorylcholine (HDPC) micelles. Appropriate standard curves for the chemical trapping method were obtained by measuring the product yields of chloride and bromide salts with 2,4,6-trimethyl-benzenediazonium (BF4) in the presence of low molecular analogs (N,N,N-trimethyl-propane sulfonate and methyl-phosphorylcholine) of the employed surfactants. The experimentally determined values for the local Br- (Cl-) concentrations were modeled by fully integrated non-linear Poisson Boltzmann equations. The best fits to all experimental data were obtained by considering that ions at the interface are not fixed at an adsorption site but are free to move in the interfacial plane. In addition, the calculation of ion distribution allowed the estimation of the degree of ion coverage by using standard chemical potential differences accounting for ion specificity. (C) 2012 Elsevier Inc. All rights reserved.

Experimental and Theoretical Study on the One- and Two-Photon Absorption Properties of Novel Organic Molecules Based on Phenylacetylene and Azoaromatic Moieties

This Article reports a combined experimental and theoretical analysis on the one and two-photon absorption properties of a novel class of organic molecules with a pi-conjugated backbone based on phenylacetylene (JCM874, FD43, and FD48) and azoaromatic (YB3p2S) moieties. Linear optical properties show that the phenylacetylene-based compounds exhibit strong molar absorptivity in the UV and high fluorescence quantum yield with lifetimes of approximately 2.0 ns, while the azoaromatic-compound has a strong absorption in the visible region with very low fluorescence quantum yield. The two-photon absorption was investigated employing nonlinear optical techniques and quantum chemical calculations based on the response functions formalism within the density functional theory framework. The experimental data revealed well-defined 2PA spectra with reasonable cross-section values in the visible and IR. Along the nonlinear spectra we observed two 2PA allowed bands, as well as the resonance enhancement effect due to the presence of one intermediate one-photon allowed state. Quantum chemical calculations revealed that the 2PA allowed bands correspond to transitions to states that are also one-photon allowed, indicating the relaxation of the electric-dipole selection rules. Moreover, using the theoretical results, we were able to interpret the experimental trends of the 2PA spectra. Finally, using a few-energy-level diagram, within the sum-over-essential states approach, we observed strong qualitative and quantitative correlation between experimental and theoretical results.

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.

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.

Conformational preferences for some 3-(4 '-substituted phenylsulfonyl)-1-methyl-2-piperidones through spectroscopic and theoretical studies

The analysis of the infrared (IR) carbonyl band of some 3-(4'-substituted phenylsulfonyl)-1-methyl-2-piperidones 1-5 bearing as substituents: OMe 1, Me 2, H 3, Cl 4 and NO2 5, supported by B3LY13/6-31G(d,p) calculations along with NBO analysis (for 1, 3 and 5) and X-ray diffraction (for 5), indicated the existence of three stable conformations i.e. quasi-axial (q-ax), syn-clinal (s-cl) and quasi-equatorial (q-eq). In the gas phase, the q-ax conformer is calculated as the most stable (ca. 88%) and the least polar, the s-cl conformer is less stable (ca. 12%) but more polar, and the q-eq conformer is the least stable (ca. 1%) and the most polar of the three conformers evaluated. The sum of the most important orbital interactions from NBO analysis and the trend of the electrostatic interactions accounts for the relative populations as well as for the v(CO) frequencies of the q-ax. s-cl and q-eq conformers calculated in the gas phase. The unique IR v(CO) band in CCl4 may be ascribed to the most stable q-ax conformer. The more intense (60%) high frequency doublet component in CHCl3 may be assigned to the summing up of the least stable q-eq and the less stable s-cl conformers, as their frequencies are almost coincident. The occurrence of only a single v(CO) band in both CH2Cl2 and CH3CN supports the fact that the v(CO) band of the two more polar conformers appear as a single band. Additional support to this rationalization is given by the single point PCM method, which showed a progressive increase of the q-eq + s-cl/q-ax population ratio going from the gas phase to CCl4, to CHCl3, to CH2Cl2 and to CN3CN. X-ray single crystal analysis of 5 indicates that this compound displays a quasi-axial geometry with respect to the [O=C-CH-S] moiety, and that the 2-piperidone ring assumes a slightly distorted half-chair conformation. In the crystal packing, molecules of 5 are arranged into supramolecular layers linked through C-H center dot center dot center dot O interactions along with it pi center dot center dot center dot pi interactions between adjacent benzene rings. (C) 2012 Elsevier B.V. All rights reserved.