Search

The doublet potential energy surfaces involved in the decomposition of the nitromethane radical anion (CH(3)NO(2) (-)) have been studied by using the multistate extension of the multiconfigurational second-order perturbation method (MS-CASPT2) in conjunction with large atomic natural orbital-type basis sets. A very low energy barrier is found for the decomposition reaction: CH(3)NO(2) (-)--[CH(3)NO(2)](-)--CH(3)+NO(2) (-). No evidence has been obtained on the existence of an isomerization channel leading to the initial formation of the methylnitrite anion (CH(3)ONO(-)) which, in a subsequent reaction, would yield nitric oxide (NO). In contrast, it is suggested that NO is formed through the bimolecular reaction: CH(3)+NO(2) (-)--[CH(3)O-N-O](-)--CH(3)O(-)+NO. In particular, the CASSCF/MS-CASPT2 results indicate that the methylnitrite radical anion CH(3)ONO(-) does not represent a minimum energy structure, as concluded by using density functional theory (DFT) methodologies. The inverse symmetry breaking effect present in DFT is demonstrated to be responsible for such erroneous prediction.

Triplet energy transfer (TET) from aromatic donors to 1,3,5,7-cyclooctatetraene (COT) is an extreme case of "nonvertical" behavior, where the transfer rate for low-energy donors is considerably faster than that predicted for a thermally activated (Arrhenius) process. To explain the anomalous TET of COT and other molecules, a new theoretical model based on transition state theory for nonadiabatic processes is proposed here, which makes use of the adiabatic potential energy surfaces (PES) of reactants and products, as computed from high-level quantum mechanical methods, and a nonadiabatic transfer rate constant. It is shown that the rate of transfer depends on a geometrical distortion parameter gamma=(2g(2)/kappa(1))(1/2) in which g stands for the norm of the energy gradient in the PES of the acceptor triplet state and kappa(1) is a combination of vibrational force constants of the ground-state acceptor in the gradient direction. The application of the model to existing experimental data for the triplet energy transfer reaction to COT from a series of pi,pi(*) triplet donors, provides a detailed interpretation of the parameters that determine the transfer rate constant. In addition, the model shows that the observed decrease of the acceptor electronic excitation energy is due to thermal activation of C=C bond stretchings and C-C bond torsions, which collectively change the ground-state COT bent conformation (D(2d)) toward a planar triplet state (D(8h)).

International audience; We present here the results of time-resolved absorption and emission experiments for 4-(dimethylamino)benzonitrile in solution, which suggest that the fluorescent intramolecular charge transfer (ICT) state may differ from the twisted ICT (TICT) state observed in transient absorption.

Some linear and nonlinear optical (NLO) properties of Ni(SCH)4 and several of its derivatives have been computed by employing a series of basis sets and a hierarchy of methods (e.g., HF, DFT, coupled cluster, and multiconfigurational techniques). The electronic structure of Ni(SCH)4 has been also analyzed by using CASSCF/CASPT2, ab initio valence bond, and DFT methods. In particular we discuss how the diradicaloid character (DC) of Ni(SCH)4 significantly affects its NLO properties. The quasidegeneracy of the two lowest-energy singlet states 1 mathg and 1 math1u, the clear DC nature of the former, and the very large number of low-lying states enhance the NLO properties values. These particular features are used to interpret the NLO properties of Ni(SCH)4. The DC of the considered derivatives has been estimated and correlated with the NLO properties. CASVB computations have shown that the structures with Ni(II) are the dominant ones, while those with Ni(0) and Ni(IV) have negligible weight. The weights of the four diradical structures were discussed in connection with the weight of the structures, where all the electrons are paired. Comparative discussion of the properties of Ni(SCH)4 with those of tetrathia fulvalene demonstrates the very large effect of Ni on the properties of the Ni-dithiolene derivatives. A similar remarkable effect on the NLO properties is produced by one or two methyl or C3S groups. The considered Ni-dithiolene derivatives have exceptionally large NLO properties. This feature in connection with their other physical properties makes them ideal candidates for photonic applications. Luis.Serrano@uv.es

Using Cholesky decomposition and density fitting to approximate the electron repulsion integrals, an implementation of the complete active space self-consistent field (CASSCF) method suitable for large-scale applications is presented. Sample calculations on benzene, diaquo-tetra- μ -acetato-dicopper(II), and diuraniumendofullerene demonstrate that the Cholesky and density fitting approximations allow larger basis sets and larger systems to be treated at the CASSCF level of theory with controllable accuracy. While strict error control is an inherent property of the Cholesky approximation, errors arising from the density fitting approach are managed by using a recently proposed class of auxiliary basis sets constructed from Cholesky decomposition of the atomic electron repulsion integrals. © 2008 American Institute of Physics.

Vertical excitation energies as well as related absolute photoabsorption oscillator strength data are very scarce in the literature for methane. In this study, we have characterized the three existing series of low-lying Rydberg states of CH4 by computing coupled cluster linear response (CCLR) vertical excitation energies together with oscillator strengths in the molecular-adapted quantum defect orbital formalism from a distorted Cs geometry selected on the basis of outer valence green function calculations. The present work provides a wide range of data of excitation energies and absolute oscillator strengths which correspond to the Rydberg series converging to the three lower ionization potential values of the distorted methane molecule, in energy regions for which experimentally measured data appear to be unavailable. Jose.V.Pitarch@uv.es Alfredo.Sanchez@uv.es Jose.Sanchez@uv.es

Psoralen photophysics has been studied on quantum chemistry grounds using the multiconfigurational second-order perturbation method CASPT2. Absorption and emission spectra of the system have been rationalized by computing the energies and properties of the low-lying singlet and triplet excited states. The S1 ππ* state has been determined to be responsible of the lowest absorption and fluorescence bands and to initially carry the population in the photophysical processes related to the phototherapeutic properties of psoralen derivatives. The low-lying T1 ππ* state is, on the other hand, protagonist of the phosphorescence, and its prevalent role in the reactivity of psoralen is suggested to be related to the elongation of the pyrone ring C3–C4 bond, where the spin density is distributed on both carbon atoms. Analysis of energy gaps and spin-orbit coupling elements indicates that the efficient photophysical process leading to the population of the lowest triplet state does not take place at the Franck-Condon region but along the S1 relaxation path. Juan.Serrano@uv.es Luis.Serrano@uv.es Manuela.Merchan@uv.es

Magnetic susceptibility and nuclear magnetic shielding at the nuclei of bis-heteropentalenes formed by two furan units ([2,3-b], [3,2-b], [3,4-b], and [3,4-c] isomers) have been computed by several approximated techniques and a large Gaussian basis set to achieve near Hartree–Fock estimates. Ab initio models of the ring currents induced by a magnetic field normal to the molecular plane were obtained for the three isomeric systems of higher symmetry, showing that the π electrons give rise to intense diamagnetic circulation. The π currents are responsible for enhanced magnetic anisotropy and strong out-of-plane proton deshielding. The theoretical findings are used to build up a “diatropicity matrix” for two fused five-membered heterocyclic systems. Inmaculada.Garcia-Cuesta@uv.es ; Alfredo.Sanchez@uv.es

The intermolecular potential-energy surface pertaining to the interaction between benzene and N2 is investigated theoretically and experimentally. Accurate intermolecular interaction energies are evaluated for the benzene–N2 van der Waals complex using the coupled cluster singles and doubles including connected triples [CCSD(T)] method and the aug-cc-pVDZ basis set extended with a set of 3s3p2d1f1g midbond functions. After fitting the energies to an analytic function, the intermolecular Schrödinger equation is solved to yield energies, rotational constants, and Raman-scattering coefficients for the lowest intermolecular levels of several benzene–N2 isotopomers. Experimentally, intermolecular Raman spectra of jet-cooled h6- and d6-benzene–N2 measured at 0.03 cm−1 resolution by mass-selective, ionization-detected stimulated Raman spectroscopies are reported. Seven intermolecular bands are assigned for each isotopomer, including transitions involving intermolecular bending and stretching vibrations and internal rotation about the benzene C6 axis. These Raman data, together with measured rotational constants and binding energies obtained by other groups on benzene–N2, agree well with the theoretical results. Such agreement points to the promise of the quantum chemical methodology employed herein in future investigations of larger van der Waals complexes. tbonpe@uv.es

The vertical spectrum of tetracyanoethylene was studied using coupled cluster theory. It was found that the lowest singlet-singlet transition, which corresponds to the excitation from the highest occupied molecular orbital (HOMO) to the lowest unoccupied molecular orbital (LUMO) excitation, occurs at 5.16 eV in the gas phase and is lowered approximately 0.1 eV due to solvent effects in acetonitrile. A parallel study on the ethene spectrum showed the quality of the basis sets and methods used, by placing the V state 7.92 eV above the ground state and giving an energy for the 0-0 transition of 5.42 eV to be compared with the experimental value of 5.50 eV.

A study of the intermolecular potential-energy surface (IPS) and the intermolecular states of the perprotonated and perdeuterated benzene–He complex is reported. From a fit to ab initio data computed within the coupled cluster singles and doubles including connected triples model for 280 interaction geometries, an analytic IPS including two- to four-body atom–atom terms is obtained. This IPS, and two other Lennard-Jones atom–atom surfaces from the literature, are each employed in dynamically exact (within the rigid-monomer approximation) calculations of J = 0 intermolecular states of the isotopomers. Rotational constants and Raman-scattering coefficients for intermolecular vibrational transitions are also calculated for each of the three surfaces. The calculated results are compared with experimental results reported herein pertaining to intermolecular Raman spectra of benzene–He. The calculated rotational constants are compared with experimental values from the literature. The fitted IPS of this work leads to calculated observables that match the experimental results very well. The IPSs from the literature are not as successful, specifically in regard to the intermolecular Raman spectra. Thomas.Bondo@uv.es

The dissociative electron transfer reaction CH3Cl+e−→CH3•+Cl− in aqueous solution is studied by using a QM/MM method. In this work the quantum subsystem (a methylchloride molecule plus an electron) is described using density functional theory while the solvent (300 water molecules) is described using the TIP3P classical potential. By means of molecular dynamics simulations and the thermodynamic integration technique we obtained the potential of mean force (PMF) for the carbon–chlorine bond dissociation of the neutral and radical anion species. Combining these two free energy curves we found a quadratic dependence of the activation free energy on the reaction free energy in agreement with Marcus’ relationship, originally developed for electron transfer processes not involving bond breaking. We also investigated dynamical aspects by means of 60 dissociative trajectories started with the addition of an extra electron to different configurations of a methylchloride molecule in solution. The PMF shows the existence of a very flat region, in which the system is trapped during some finite time if the quantum subsystem quickly losses its excess kinetic energy transferring it to the solvent molecules. One of the most important factors determining the effectiveness of this energy transfer seems to be the existence of close contacts (hydrogen bonds) between the solute and the solvent. Estanislao.Silla@uv.es Ignacio.Tunon@uv.es

The lowest-energy band of the electronic spectrum of pyrrole has been studied with vibrational resolution by using multiconfigurational second-order perturbation theory (CASPT2) and its multistate extension (MS–CASPT2) in conjunction with large atomic natural orbital-type basis sets including Rydberg functions. The obtained results provide a consistent picture of the recorded spectrum in the energy region 5.5–6.5 eV and confirm that the bulk of the intensity of the band arises from a ππ∗ intravalence transition, in contradiction to recent theoretical claims. Computed band origins for the 3s,3p Rydberg electronic transitions are in agreement with the available experimental data, although new assignments are suggested. As illustrated in the paper, the proper treatment of the valence–Rydberg mixing is particularly challenging for ab initio methodologies and can be seen as the main source of deviation among the recent theoretical results as regards the position of the low-lying valence excited states of pyrrole. serrano@uv.es ; merchan@uv.es

Glycine and alanine aminoacids chemistry in solution is explored using a hybrid three parameters density functional (B3PW91) together with a continuum model. Geometries, energies, and vibrational spectra of glycine and alanine zwitterions are studied at the B3PW91/6-31+G∗∗ level and the results compared with those obtained at the HF and MP2/6-31+G∗∗ levels. Solvents effects are incorporated by means of an ellipsoidal cavity model with a multipolar expansion (up to sixth order) of the solute’s electrostatic potential. Our results confirm the validity of the B3PW91 functional for studying aminoacid chemistry in solution. Taking into account the more favorable scaling behavior of density functional techniques with respect to correlated ab initio methods these studies could be extended to larger systems. silla@uv.es ; tunon@uv.es

A theoretical study of the vibrational structure of the math 2A′ ground and math 2A′ excited states of the formyl radical, HCO, and its deuterated form, DCO, has been performed. The potential energy surfaces have been computed by means of a multiconfigurational perturbative method, CASPT2. The computed geometries and the harmonic and anharmonic frequencies are successfully compared to the available experimental information. The vibrational intensities of the transition math 2A′↔math 2A′ have been computed both for absorption and emission. The results lead to accurate determinations of several structural parameters and some reassignments of the vibrational transitions of the so-called hydrocarbon flame bands of the formyl radical. serrano@uv.es

Ab initio results for the electronic spectrum of all-trans-retinal and its truncated model 3-methyl-all-trans (10-s-cis)-2,4,6,8,10-undecapentaen-1-al are presented. The study includes geometry determination of the ground state. Vertical excitation energies have been computed using multiconfigurational second-order perturbation theory through the CASPT2 formalism. The lowest singlet excited state in gas phase is predicted to be of nπ∗ character. The lowest triplet state corresponds, however, to a ππ∗ state. The most intense feature of the spectrum is due to the strongly dipole-allowed ππ∗ transition, in accordance with the observed maximum in the one-photon spectra. The vertical excitation energies of the Bu- and Ag-like states are found close, the latter ≈ 1 eV higher than the maximum in the two-photon spectra. Solvent effects and nonvertical nature of the observed maximum in the two-photon spectra are invoked in rationalizing the deviation with respect to the best present estimate for the Ag-like state. In addition, qualitative aspects of the one-bond photoisomerization about the C11 = C12 double bond of retinal are considered. The overall isomerization picture from 11-cis into all-trans-retinal, as taking place mainly along the triplet manifold, agrees with experimental evidence. merchan@uv.es ; gonzaler@uv.es

Limited previous experience with the mean value total dressing (MVTD) method had shown that MVTD energies for closed shell systems are generally better than CCSD(T) ones compared to FCI. The method, previously published as total dressing 2′(td-2′), is based on the single reference intermediate Hamiltonian theory. It is not a CC method but deals in a great part with the same physical effects that CC methods that incorporate amplitudes of triples such as CCSDT or its CCSDT-1n approaches. A number of test calculations comparing to diverse CC methods, as well as FCI and experiment when available, have been performed. The tests concern equilibrium energies in NH3 and CH2, equilibrium energies and distances in some diatomics (BF, NO+, CN+, C2, BeO), different bond breaking situations (H2O, BH, HF, SiH2) and spectroscopic properties of different bonding conditions (Li2, LiNa, LiBe+, NeH+, and O3). The results are in general closer to the full CCSDT ones in the equilibrium regions and close to CCSDT-1 along most dissociation curves. A few exceptions to this rule are analyzed with the help of an approach to MVTD that does not take into account the effects of linked quadriexcitations. Such analysis suggests the interest of improving the treatment of effects of linked triples in the MVTD model. The separate contributions of linked and unlinked triples and quadruples are also analyzed for some of the above diatomics representing different behaviors of bond breaking. The interest of such analysis is illustrated in the NeH+ molecule. The MVTD results show, in general, a high quality, provided that the nature of the correlation problem does not become largely multiconfigurational, as occurs in multiple bond dissociation or in the asymmetric stretching of ozone. garciain@uv.es ; sanchezm@uv.es ; sanchez@uv.es