Ab initio investigation of the n --> pi* transitions in thiocarbonyl dyes

The n --> pi* transitions in more than 100 thiocarbonyl dyes have been calculated with an ab initio procedure relying on the combination of time-dependent density functional theory (TD-DFT) for evaluating excited states and the polarizable continuum model (PCM) for modeling the bulk solvent effects on both the geometrical and electronic structures. Two hybrid functionals (B3LYP and PBE0) and several basis sets, some including f polarization functions, have been used. B3LYP provides the most accurate raw estimates, but once simple linear regression is performed, both functionals give similar results with a small advantage for PBE0. By use of the latter, the mean absolute deviation with respect to experiment is limited to 0.06 eV whereas less than 20% of the estimates differ from absorption data by more than 0.10 eV. To assess the validity limits of our model, compounds containing multiple C=S chromophores have also been considered.     

Coupled cluster calculation of the n --> pi* electronic transition of acetone in aqueous solution

The combined linear response coupled cluster/molecular mechanics (CC/MM) scheme including mutual polarization effects in the coupling Hamiltonian is applied together with supermolecular CC methods to the study of the gas-to-aqueous solution blue shift of the n --> pi* excitation energy in acetone. The aug-cc-pVDZ basis set is found to be adequate for the calculation of this excitation energy. In the condensed phase, the shift in the excitation energy is obtained by statistical averaging over 800 solute-solvent configurations extracted from a molecular dynamics simulation. We find the shift to be around 1100-1200 cm(-1) depending on the specific model used to describe solvent polarization. The importance of including explicit polarization in both the molecular dynamics simulation as well as the CC/MM calculations is emphasized. Furthermore, the significant dependence of the excitation energy on the CO bond length of acetone is discussed.     

On the influence of nuclear fluctuations on calculated NMR shieldings of benzene and ethylene: a Feynman path integral–ab initio investigation*

The absolute magnetic shieldings of benzene and ethylene have been theoretically studied under the conditions of thermal equilibrium, i.e., under explicit consideration of the nuclear degrees of freedom. For this purpose we have combined the Feynman path integral quantum Monte Carlo (PIMC) formalism with the gauge‐including atomic orbital (GIAO) approach in the Hartree–Fock (HF) approximation. The HF operator has been employed to derive the NMR parameters of the two hydrocarbons via an ensemble averaging over large sets of molecular configurations that are populated in thermal equilibrium. The nuclear fluctuations are responsible for a deshielding of the nuclei relative to the shieldings at the vibrationless minimum of the potential energy surface (PES). The influence of the nuclear degrees of freedom is largest for the isotropic part of the ¹³C shielding tensor. The theoretical results can be explained on the basis of simple geometrical considerations. The bond lengths in thermal equilibrium are larger than the bond lengths at the minimum of the PES. This length enhancement is the prerequisite for a deshielding of the nuclei in thermal equilibrium. The vibrational corrections of the nuclear magnetic resonance (NMR) parameters of benzene and ethylene are quantum driven; classical thermal degrees of freedom of the nuclei are of minor importance. Conceptual problems of theoretical studies of NMR parameters on the basis of a single molecular geometry are emphasized. The influence of the spatial uncertainty of the nuclei becomes decisive in molecules with light atoms. It is pointed out that the combination of the PIMC formalism with electronic Hamiltonians of state‐of‐the‐art quality renders possible accurate determinations of NMR parameters. © 2002 John Wiley & Sons, Inc. Int J Quantum Chem 86: 280–296, 2002     

An ab initio study of solvent polarity and hydrogen bonding effects on the nitrogen NMR shieldings of N,N‐dimethylacetamidine

Density functional theory combined with the polarizable continuum model (PCM) and continuous set of gauge transformations method is applied to investigate the effects of solvent polarity on the nitrogen NMR shieldings of N, N‐dimethylacetamidine. Hydrogen bonding effects on shielding are likewise calculated using a supermolecule approach, where the imino group of the solute is hydrogen bonded with solvent. Theoretical results are compared with published experimental data. The PCM shielding calculations utilizing PCM‐optimized solute geometries yield results comparable to those obtained with the supermolecule approach. Geometry optimization of the solute appears to be more important in PCM shielding calculations than in the supermolecule approach. The large solvent shifts observed in water can only be reproduced when the N·H distance used in the calculation indicates full proton transfer from water to the imino nitrogen of the solute. Copyright © 2002 John Wiley & Sons, Ltd.     

Quantum wavepacket ab initio molecular dynamics: an approach for computing dynamically averaged vibrational spectra including critical nuclear quantum effects

We have introduced a computational methodology to study vibrational spectroscopy in clusters inclusive of critical nuclear quantum effects. This approach is based on the recently developed quantum wavepacket ab initio molecular dynamics method that combines quantum wavepacket dynamics with ab initio molecular dynamics. The computational efficiency of the dynamical procedure is drastically improved (by several orders of magnitude) through the utilization of wavelet-based techniques combined with the previously introduced time-dependent deterministic sampling procedure measure to achieve stable, picosecond length, quantum-classical dynamics of electrons and nuclei in clusters. The dynamical information is employed to construct a novel cumulative flux/velocity correlation function, where the wavepacket flux from the quantized particle is combined with classical nuclear velocities to obtain the vibrational density of states. The approach is demonstrated by computing the vibrational density of states of [Cl-H-Cl]-, inclusive of critical quantum nuclear effects, and our results are in good agreement with experiment. A general hierarchical procedure is also provided, based on electronic structure harmonic frequencies, classical ab initio molecular dynamics, computation of nuclear quantum-mechanical eigenstates, and employing quantum wavepacket ab initio dynamics to understand vibrational spectroscopy in hydrogen-bonded clusters that display large degrees of anharmonicities.     

Solvent effects on the UV (n --> pi*) and NMR (13C and 17O) spectra of acetone in aqueous solution. An integrated car-parrinello and DFT/PCM approach

The accurate reproduction of ultraviolet and nuclear magnetic resonance spectra of acetone in aqueous solution is used as a test of an integrated computational tool rooted in the density functional theory, the polarizable continuum model, and the Car-Parrinello dynamics. The analysis and solution of several conceptual and practical issues results in a robust and effective approach, which also can be used by nonspecialists and provides a general and powerful complement to experimental techniques.     

Vibrational dynamics of the hydrogen bond in H(2)S-HF: Fourier-transform-infrared spectra and ab initio theory

The rotationally resolved infrared spectrum of the hydrogen bonded complex H(2)S-HF and of its isotopomer D(2)S-DF in the HF/DF stretching range have been observed in a supersonic jet Fourier-transform infrared (FTIR) experiment and indicate a predissociation lifetime of 130 ps for H(2)S-HF. Complementary spectra taken at a temperature of 190 K in a cell without resolved rotational structure indicate the presence of strong anharmonic couplings between low frequency intermolecular modes and the HF donor stretch mode previously observed in other complexes with heavier acceptor molecules without rotational fine structure. The anharmonic analysis of the hot band progressions and of the rotational data confirm the coupling mechanism. The coupling constants and the absolute frequency of the hydrogen bonded stretch mode are in excellent agreement with theoretical predictions based on adiabatic variational calculations on potential surfaces computed at MP2 and CCSD(T) level. Complementary calculations with a perturbational approach further confirm the coupling model.     

On the charge distribution in ethanes and disilanes and correlations with equilibrium bond lengths; an ab initio study

The equilibrium electronic wave-functions for a series of fluoro- and chloro-ethanes and disilanes of general formula M₂H_6−nX_n, (M=C, Si; X=F, Cl), were analysed by the most commonly used methods for electron distribution, using the Mulliken and Löwdin populations, natural atomic orbital (NAO) populations and atoms in molecules (AIM) electron densities. Although the numerical values for local atomic charges vary greatly, all the methods correlate, but in markedly differing ways. The Mulliken charges seem the most selective in relation to systematic change of substituents in the current type of molecular structure. A number of examples occur where the AIM charges at C, Si centres are effectively identical in different molecules, where some differences might have been anticipated. These are often distinguished by Mulliken populations. The fluoroethanes exemplify this, since a plot of the AIM charges (for example on either the F or H centres) against the Mulliken charges for all members of the series, shows three nearly parallel lines, corresponding to those centres with 0, 1 or 2 fluorine atoms on the centre under study. The bond critical points at which the AIM charges are determined seem to be counter to intuition in some cases. This is a density rather than atomic orbital size issue however. The Mulliken and NAO charges seem more reasonable than those from the AIM method. There is an unexpected correlation of the local bond dipoles from the Mulliken analyses, with the calculated equilibrium bond lengths. These correlations lead to bond length values for the non-polarised bonds MX, which agree with data based on covalent radii for some bonds.     

Two-way effects between hydrogen bond and intramolecular resonance effect: an ab initio study on complexes of formamide and its derivatives with water

Ab initio calculations up to MP2/aug-cc-pVTZ//MP2/cc-pVTZ level, including natural charge population and natural resonance theory analyses, have been carried out to study the two-way effects between hydrogen bond (H-bond) and the intramolecular resonance effect by using the H-bonded complexes of formamide ( FAO) and its derivatives ( FAXs, X represents the heavy atoms in the substituent groups, CH 2, NH, SiH 2, PH, and S) with water as models. Unlike NH 3 and NH 2CH 3 which prefer being H-bond acceptors ( HA) to form H-bond with water, the amino groups in the six monomers, because of the resonance effect, prefer being H-bond donors ( HD) rather HA. Six monomers can all form HD complexes with water, and only two ( FAC and FASi) with the weakest resonance effect are able to form HA complexes with water. The HD H-bond and resonance effect enhance each other (positive two-way effects) whereas the HA H-bond and resonance effect weaken each other (negative two-way effects). The H-bond energies in the six HD complexes are nearly linearly correlated with the weights of the dipolar resonance in Pauling's model and the N-C bond lengths; the correlation coefficients are 0.91 and 0.93, respectively. The positive two-way effects also happens in FAO-water complex, in which the FAO CO group serves as HA ( HA co ). Interestingly, when the HD and HA co H-bonds are present in FAO H-bond complex simultaneously, the enhancements are much more significant, and the energies of the two types of H-bonds are much larger than those when only one type of H-bond is present, reflecting the cooperative effects. By using the knowledge to the two-way effects, we computationally designed a molecule ( FAO- BH 3 ) to increase H-bond energy. Because of the oxygen lone pair donation to the empty pi orbital of BH 3, FAO- BH 3 has a much stronger resonance effect than FAO. As a result, the H-bond energy (-5.55 kcal/mol) in HD H 2O ... FAO- BH 3 complex is much greater than the -3.30 kcal/mol in the HD H 2O...FAO complex. The two-way effects can be rationalized as follows: the resonance effect leads to intramolecular charge shifts in the monomers which facilitate or prevent the charge donation or acceptation of their H-bond partners. Therefore, the H-bonds are strengthened or weakened. In reverse, the charge donations or acceptations of their H-bond partners facilitate or prevent the intramolecular charge shifts in the monomer moieties, which enhance or weaken the resonance effect. The understanding to the two-way effects may be helpful in drug design and refinement by modulating the H-bond strength and in building empirical H-bond models to study large biological molecules. The study supports Pauling's resonance model.     

Structure and dynamics of phosphate ion in aqueous solution: an ab initio QMCF MD study

A simulation of phosphate in aqueous solution was carried out employing the new QMCF MD approach which offers the possibility to investigate composite systems with the accuracy of a QMMM method but without the time consuming creation of solute-solvent potential functions. The data of the simulations give a clear picture of the hydration shells of the phosphate anion. The first shell consists of 13 water molecules and each oxygen of the phosphate forms in average three hydrogens bonds to different solvent molecules. Several structural parameters such as radial distribution functions and coordination number distributions allow to fully characterize the embedding of the highly charged phosphate ion in the solvent water. The dynamics of the hydration structure of phosphate are described by mean residence times of the solvent molecules in the first hydration shell and the water exchange rate.     

Gas phase and water solution conformational analysis of the herbicide diuron (DCMU): an ab initio study

In the present work, the conformational equilibrium for the herbicide diuron (DCMU) has been investigated using high level ab initio calculations. The solvent effect was included through two different continuum models: (1) the real cavity IPCM method and (2) the standard dipole Onsager model SCRF. The effect due to solute-solvent hydrogen-bond interactions was analyzed considering a hybrid discreet-continuum model. At the Hartree-Fock level, the gas phase results showed that only the trans forms (A and B) are present in the equilibrium mixture, with the relative concentrations found to be 33% (A) and 67% (B) (HF/6-311+G**//6-31G**). When the electronic correlation effect is included (MP2/6-31G*//HF/6-31G*), a relative stabilization of the cis forms was observed, with the conformational distribution calculated as 38% (A), 50% (B), 6% (C) and 6% (D). The trans conformations were found to be completely planar, which has been considered to be a prerequisite for the herbicide binding. In water solution, the trans conformation A should be the most abundant conformer, the IPCM and SCRF values being ca. 100% and ca. 85% respectively. The IPCM calculations with the isodensity level set to 0.0005 present a conformational distribution close to that obtained from the hybrid model [92% (A) and 8% (B)], which has been considered our best solvent approach. Regarding the biological action of urea-type herbicides, the results presented here are important, because some QSAR studies have suggested that the partition coefficient is related to the herbicide activity, so the conformational equilibrium may play a role in the biological action. Received: 23 February 1998 / Accepted: 28 May 1998 / Published online: 19 August 1998     

Electronic spectra of adenine and 2-aminopurine: an ab initio study of energy level diagrams of different tautomers in gas phase and aqueous solution

Ground and lowest two singlet excited state geometries of four tautomeric forms (N9H, N7H, N3H and N1H) of each of adenine and 2-aminopurine (2AP) were optimized using an ab initio approach employing a mixed basis set (6-311 + G* on the nitrogen atom of the amino group and 4-31G basis set on the other atoms). Excited states were generated employing configuration interaction involving single electron excitations (CIS). Subsequently, the different species were solvated in water employing the self-consistent reaction field (SCRF) approach along with the corresponding gas phase optimized geometries. Thus the observed absorption and fluorescence spectra of adenine and 2AP have been explained successfully. It is concluded that both the N9H and N7H forms of 2AP would contribute to absorption and fluorescence spectra. Further, the fluorescence of 2AP would be absorbed by its cation in which both the N9 and N7 atoms are protonated, the fluorescence of which can have an anti-Stokes component. Among the different tautomers of adenine, the N9H form would be present dominantly in the ground state in aqueous solutions but the N7H form would be produced by energy transfer and subsequent fluorescence. The N3H form of adenine appears to be responsible for the observed absorption near 300 nm by its solutions intermittently exposed to ultraviolet radiation. The rings of the different species related to 2AP and adenine remain almost planar in the pi-pi* and n-pi* singlet excited states as in the ground state. The pyramidal character of the amino group is usually less in the pi-pi* excited states than that in the corresponding ground or n-pi* excited states. Molecular electrostatic potential (MEP) maps of the molecules provide useful clues regarding phototautomerism.     

The effects of fluorine and chlorine substitution on bond lengths in ethanes and disilanes: comparisons of ab initio and experimental information

Ab initio and some density functional theory calculations of bond lengths in fluoro- and chloro-ethanes and disilanes are reported with a precision of ±0.0001 Å under strictly comparable conditions. The resulting changes in MH and MX (M=C, Si; X=F, Cl) bond length are analysed for the effects of halogens substituted in geminal (), or vicinal (gauche or trans) positions. The shortening effect of halogen on an MH bond is markedly reduced or even reversed by the introduction of electron correlation at the MP2 or B3LYP level. MX bonds are little affected. gauche halogen consistently shortens both MH and MX bonds, while trans halogen has no effect on an MH bond but a small and variable effect on the MX bond.

The reality of these calculated changes in bond length is tested in two ways. MH bond lengths are plotted against experimental values of the isolated stretching frequencies ν^isMH, which themselves correlate well with experimental r₀ bond lengths. Agreement on the resulting substituent effects is generally good for the gauche and trans effects of halogen but variable for effects. Unobserved ν^isMH values are predicted from computed bond lengths in fluoroethanes, chloroethanes and chlorodisilanes.

Calculated MX and MM bond lengths are compared with experimental values, notably those from electron diffraction studies amongst the ethanes. Most calculations underestimate the changes found experimentally in CF and CCl bond lengths. CC bond length changes are underestimated in fluoroethanes and overestimated in the chloro-compounds.

The ‘offset’ value (r_e(calc)−r_e(true)) for a CH or SiH bond calculated with a given basis set and level of theory in most cases varies markedly throughout the series of compounds. The same is true for CF, CCl, CC and SiSi bonds if the corresponding offset values for the r_a lengths are constant.

The need is stressed for extended experimental work on many of the compounds, especially the disilanes. It is recommended that structures should be refined with ab initio derived constraints on the bond lengths involved and differences between spectroscopic and diffraction-based geometries reconciled through the calculation of r_z structures.     

Hydrogen bonding and covalent effects in binding of cisplatin to purine bases: ab initio and atoms in molecules studies

Ab initio and density functional calculations are employed to investigate the role of hydrogen bonding in the binding of cisplatin to the purine bases guanine and adenine. Through the use of the theory of atoms in molecules (AIM), it is shown that hydrogen bonds are ubiquitous in such systems, with N-H...N and N-H...Cl interactions present in addition to the expected N-H...O. This in turn means that the known stability of cisplatin-guanine complexes cannot be ascribed solely to hydrogen bonding and allows decomposition of total binding energy into contributions from covalent and hydrogen bonds. To do so, a new method for predicting hydrogen bond energies from bond critical point properties is proposed, employing partial least-squares analysis to remove the family dependence of simple models. Still more hydrogen bond motifs are found in bifunctional complexes of the general type purine-[Pt(NH(3))(2)](2+)-purine, including purine...purine contacts, though again the energetics of these are insufficient to explain the observed trends in stability. Finally, the effect of platination on the pairing of guanine with cytosine is studied in a similar manner, revealing large redistributions of hydrogen bonding but surprisingly small overall changes in pairing energy.     

Electronic spectra of linear isoelectronic clusters C2n+1S and C2n+1Cl+ (n = 0-4): an ab initio study

Structures and stabilities of linear carbon chains C2n+1S and C2n+1Cl+ (n=0-4) in their ground states have been investigated by the CCSD and B3LYP approaches. The CASSCF calculations have been used to determine geometries of selected excited states of both isoelectronic series. Linear C2n+1S cluster has a cumulenic carbon framework, whereas its isoelectronic C2n+1Cl+ has a dominant character of acetylenic structure in the vicinity of terminal Cl. The vertical excitation energies of low-lying excited states have been calculated by the CASPT2 method. Calculations show that the excitation energies have nonlinear size dependence. The 2(1)Sigma+<--X1Sigma+ transition energy in C2n+1S has a limit of 1.78 eV, as the chain size is long enough. The predicted vertical excitation energies for relatively strong 1(1)Pi<--X1Sigma+ and 2(1)Sigma+<--X1Sigma+ transitions are in reasonable agreement with available experimental values. The spin-orbit effect on the spin-forbidden transition in both series is generally small, and the enhancement of the spin-forbidden transition by spin-orbit coupling exhibits geometrical and electronic structural dependence.     

An ab initio investigation on (CO2)n and CO2(Ar)m clusters: geometries and IR spectra

An ab initio investigation on CO(2) homoclusters is done at MPWB1K6-31++G(2d) level of theory. Electrostatic guidelines are found to be useful for generating initial structures of (CO(2))(n) clusters. The ab initio minimum energy geometries of (CO(2))(n) with n=2-8 are T shaped, cyclic, trigonal pyramidal, tetragonal pyramidal, tetragonal bipyramidal, pentagonal bipyramidal, and pentagonal bipyramid with one CO(2) molecule attached to it. A test calculation on (CO(2))(20) cluster is also reported. The geometric parameters of the energetically most favored (CO(2))(n) clusters match quite well their experimental counterparts (wherever available) as well as those derived from molecular dynamics studies. The effect of clustering is quantified through the asymmetric C-O stretching frequency shift relative to the single CO(2) molecule. (CO(2))(n) clusters show an increasing blueshift from 1.8 to 9.6 cm(-1) on increasing number of CO(2) molecules from n=2 to 8. The energetics and geometries of CO(2)(Ar)(m) clusters have also been explored at the same level of theory. The geometries for m=1-6 show a predominant T type of the argon-CO(2) molecule interaction. Higher clusters with m=7-12 show that the argon atoms cluster around the oxygen atom after the saturation of the central carbon atom. The CO(2)(Ar)(m) clusters exhibit an increasing redshift in the C-O asymmetric stretch relative to CO(2) molecule of 0.7-5.6 cm(-1) with increasing number of argon atoms through m=1-8.     

Synthesis, structural characterization and ab initio calculation of dipyridyltetraazathiapentalene: a highly conjugative polycyclic molecule with hypervalent N–S–N bond

A highly conjugative polyheterocyclic compound, tetraazathiapentalene fused with pyridine rings, was synthesized by reacting 2-aminopyridine with carbon disulfide. The single crystal X-ray determination reveals that the molecule crystallizes in monoclinic space group C2/c, with the following unit cell dimensions: a=11.062(2), b=9.030(1), c=20.898(5) Å, β=102.98(1)°, V=2034.00(3) Å³, Z=8, and that a hypervalent N–S–N bond exists in the molecule. Ab initio calculations predict its IR and ¹H NMR spectra that are coincident with the experimental ones and reveal the bonding nature of the hypervalent N–S–N bond and the electronic structure of the molecule.     

Simulation of the photodeactivation of formamide in the nO-pi* and pi-pi* states: an ab initio on-the-fly surface-hopping dynamics study

The short-time photodynamics (1 ps) of formamide in its low-lying singlet excited n(O)-pi(*) and pi-pi(*) states have been investigated by the direct trajectory surface-hopping method based on multiconfigurational ab initio calculations. The simulations showed that in both states, the primary deactivation process is C-N bond dissociation. In the ground state, the energy is transferred to (a) translational motion of the HCO and NH(2) fragments, (b) additional C-H dissociation from the vibrationally hot HCO fragment, or (c) formation of NH(3) and CO. In addition to the C-N dissociation pathway, C-O bond fission is found to be an additional primary deactivation path in the pi-pi(*) dynamics. From fractional occupations of trajectories, lifetimes of formamide were estimated: tau(S(1))=441 fs and tau(S(2))=66 fs.     

Hydrogen bonding and solvent effects on the lowest 1(n, pi*) excitations of triazines in water

Our method for estimating solvent effects on electronic spectra in media with strong solute-solvent interactions is applied here to calculate the absorption and fluorescence solvatochromatic shifts of dilute triazines in water. First, the ab initio CASSCF method is used to estimate the gas-phase electronic excitation properties and state charge distributions; second, Monte Carlo simulations are performed to elucidate liquid structures around the ground and excited state solute; finally, the solvent shift is evaluated based on the gas-phase charge distributions and the explicit solvent structures. For the dilute triazine solutions, simulations predict one linear (different) hydrogen bond attached to each nitrogen atom. Upon the first (1)(n, pi*)electronic excitation one hydrogen bond is completely broken. For the absorption and fluorescence spectra, our calculations demonstrated that the specific solvent-solute interaction, in any electronic state, plays a critical role in the determination of solvent shifts.