Theoretical study and rate constant calculation for the O(3P) + C2H5CN reaction

The complicated microscopic reaction mechanisms of O(³P) with C₂H₅CN on the ground electronic state energy surface have been investigated at the G3(MP2) level of theory based on the geometric parameters optimized at the B3LYP/6-311 + G(d, p) level. Two kinds of H-abstraction and addition–elimination channels are considered, namely methylene-H abstraction, methyl-H abstraction, C-addition/elimination and N-addition/elimination. The kinetics of the title reaction have been studied using the TST and multichannel RRKM methodologies over a wide temperature range of 200–2000 K. The results show that the methylene-H abstraction process is predominant for the whole reaction. With an increase of temperature, H-abstraction from the methyl position channel should be taken into account. The C-addition/elimination process provides a few contributions to the title reaction compared with two kinds of H-abstraction channels over the whole temperature region and the N-addition/elimination channel can be negligible due to the high entrance barrier and unstable products.     

A theoretical investigation of the oxidation states of palladium complexes and their role in the carbonylation reaction

Two different, yet related, topics are discussed: (i) the reduction of palladium (II) in Pd(OAc)₂ complexes reacting with phenyl phosphines and leading to Pd(0) phosphine complexes, and (ii) the carbonylation reaction of allyl chlorides catalyzed by these Pd(0) species. The results show that the overall reduction is an exothermic process that can be accomplished along two different reaction paths, one being clearly favoured over the other. Similarly, three different channels have been determined for the carbonylation reaction that primarily differ in the timing and the way in which the reacting species bind the metal. In the first path (the σ-path), the allyl fragment interacts very weakly with the metal, whereas the CO molecule strongly binds it and reacts with the allyl. The second channel (the π–η³ pathway) is characterized by a π–η³ interaction between the allyl fragment and the palladium, to which the CO molecule binds, before the two units react affording the product. In both cases, two consecutive migrations of the chlorine ‘assist’ the course of the reaction. In the third case (the η² pathway) the allyl fragment initially enters the palladium coordination sphere, and the CO molecule then simultaneously binds it and the phosphorous atom of one phosphine ligand. The first two paths are favoured.     

The accuracy of local MP2 methods for conformational energies

The relative energies of 95 conformers of four peptide models are studied using MP2 and LMP2 methods and correlation consistent basis sets ranging from double-zeta to augmented quintuple-zeta quality. It is found that both methods yield quite similar results, and the differences between MP2 and LMP2 decrease systematically with increasing basis set. Due to reduced intramolecular basis set superposition effects (BSSE), the LMP2 results converge more slowly to the basis set limit for most of these rather small systems. However, for larger peptides, the BSSE has a very large effect on the energy difference between extended and helical structures, leading to a very strong basis set dependence of the canonical MP2 results. It is demonstrated for alanine octapeptides that the basis set error exceeds 30 and 20kJ mol^?1, respectively, if augmented double-zeta and triple-zeta basis sets are used. On the other hand, the LMP2 results are only slightly affected by the basis set size, and, even with augmented double-zeta basis sets, reasonably accurate results are obtained. Furthermore, for the larger systems, the computation times for the LMP2 calculations are shown to be up to one order or magnitude shorter than for canonical MP2 calculations with the same basis set.     

Penalty functions for combining coupled-cluster and perturbation amplitudes in local correlation methods with optimized orbitals

Local active space correlation models based on the coupled-cluster doubles (CCD) model like Generalized Valence Bond Perfect Pairing (GVB-PP) and Imperfect Pairing (IP) are attractive methods for treating electron correlation, because they are computationally inexpensive and can describe strong correlations. However, they suffer from symmetry-breaking (SB) in systems with multiple resonance structures, which arises due to neglected correlations. We investigate the extent to which these problems can be removed by using second-order perturbation theory (PT) for weak correlations coupling three different electron pairs, and (infinite-order) coupled-cluster (CC) theory for stronger correlations involving electrons in only one or two pairs. The resulting Three-Pair Corrected Imperfect Pairing (TIP) method is explored here, and it is shown that to robustly combine CC and PT it is necessary to modify several aspects of the basic method. Most importantly, a penalty function term is introduced to ensure the PT amplitudes remain small. Comparison against CC treatment of the three-pair correlations suggests penalty terms will be beneficial for any hybrid CC/PT method that includes orbital optimization. The TIP method greatly reduces SB in aromatic hydrocarbons and recovers a significantly higher fraction of the valence electron correlation energy than IP.     

Concentration and temperature dependent interaction studies using dielectric and thermodynamic methods on mixtures of anisole with o-toluidine and m-toluidine

Molecular interactions in mixtures of anisole with o-toluidine and anisole and m-toluidine have been studied at three different temperatures using the dielectric method with measurements of the static permittivity and permittivity at optical frequency. From the measured values, the Kirkwood correlation parameter, Bruggeman parameter, excess permittivity and thermodynamic excess free energy were computed for the mixtures. Positive and negative values of excess permittivity were obtained for both mixtures. The excess free energy for the anisole+o-toluidine mixture is positive at all three temperatures, whereas mixed values (positive and negative) are obtained for the anisole+m-toluidine mixture. Alignment of the dipoles in both mixtures was identified by Kirkwood factors. The investigation shows that the interaction between the components changes systematically with concentration and temperature and the change is minimum.     

Geometry optimization in the presence of external forces: a theoretical model for enforced structural changes in molecules

Recently, we proposed a very simple quantum chemical model to simulate the effect of external forces acting on a single molecule [Mol. Phys. 107, 2403 (2009)]. It is based on optimizing the geometry of a molecule with an external force applied to selected pairs of nuclei. In this study we extend this model by considering interactions of external forces not only with the nuclei but also with their electrons, in particular their core electrons, which can be viewed as ‘rigidly’ connected to a nucleus. In the proposed revised model an external force acts on an object which consists of the nucleus of an atom and its 1s core electrons. It is shown in this study that such a model predicts the same conformational (structural) changes in a molecule as our simpler model where the external forces interact with the nuclei only. However, the magnitude of the forces required to cause these changes is now lower and within the range of forces used in real AFM experiments.     

Theoretical studies on the anisotropy of the first hyperpolarizabilities of one- and two-dimensional charge-transfer molecules: role of solvent polarity

The effects of solvent polarity on the geometry and static first hyperpolarizability of one- and two-dimensional charge transfer molecules based on para-nitroaniline (PNA) and 1,3-diamino-4,6-dinitrobenzene (DADB) were studied using a polarizable continuum model (PCM) at the B3LYP/6-31+G(d,p) level of theory. The inclusion of solvent does not change the geometry significantly, but produces a change in the calculated results of the electronic structures and nonlinear optical properties. It is found that solvent polarity has an important influence on the first hyperpolarizability, whereas its influence is negligible on the nonlinearity anisotropy parameter both for PNA and DADB, due to the fact that the off-diagonal component and the diagonal component of the first hyperpolarizability possess similar solvent effect characteristics.     

Effect of curve crossing on absorption and resonance Raman spectra: analytical treatment for delta-function coupling in parabolic potentials

We propose an exactly solvable model for the two-state curve-crossing problem. Our model assumes the coupling to be a delta function. It is used to calculate the effect of curve crossing on the electronic absorption spectrum and the resonance Raman excitation profile.     

A generalized Langevin algorithm for studying permeation across biological ion channels

We present a new leapfrog algorithm for the numerical solution of the generalized Langevin equation (GLE) in the case where the friction kernel is exponentially decaying. Like other leapfrog and Verlet algorithms, our algorithm is second order in velocity and third order in position. It is relatively easy to implement compared with other available algorithms, and would therefore make a good candidate for exploring the effects of finite memory time-scales in situations where modelling the precise functional form of the memory kernel was not important. We have tested this algorithm on a one-dimensional barrier crossing model, and found good asymptotic agreement with limits obtained using Brownian dynamics (BD) simulations, as well as with a theoretical asymptotic limit. We have also used the algorithm to perform a more sophisticated simulation of ion conduction through a KcsA channel. The results are a close match to corresponding results obtained using the Langevin equation, thereby helping to justify the use of Brownian dynamics in KcsA and other similar ion channels.     

Quantum chemical study of the ground-state alcoholic complexation of selected dual luminescent compounds

The ground- and excited-state features of dual luminescent molecules are strongly influenced by the presence of alcoholic additives. Selected ground-state properties of methanol and 1,1,1,3,3,3-hexafluoro-propan-2-ol complexes of 4-aminobenzonitrile, 4-aminopyridine and aniline derivatives were obtained by quantum chemical calculations. The formation enthalpies of the complexes are the most exothermic when the cyano or ring nitrogen interacts with the hydroxyl group of the alcohols. The binding energies are almost doubled when the fluorinated reactant is the hydrogen bond donor. Parallel to the enhancement of the stabilization of the complex, the ground-state dipole moment also increases notably with complex formation. In principle, consideration of this increase is essential in interpretation of the solvatochromatic behaviour of the complexed species.     

苯胺分子S_1态振动光谱及ab initio计算(英文)

利用紫外激光,结合超声速脉冲分子技术和飞行时间质谱仪,在287～296 nm波长范围内实验研究了苯胺分子的共振双光子电离光谱(R2PI).实验中观测到对应S_1←S_0跃迁的0-0带出现在293.86nm(3.4029 cm~(-1))处,并测得对应S_1态的若干振动模和来自S_0态的热带.为了对S_1态的振动模进行标定,分别在HF/6-31+G(d,p)和CIS/6-31+G(d,p)基组水平上对苯胺在S_0态的构型进行优化和在S_1态对其振动频率进行分析计算.另外,利用自然键轨道分析(NBO)讨论了电子激发过程.结果显示,来自氨基氮原子上的一个电子从对应的孤对电子轨道激发到苯环的π反键轨道上,恰好对应的是S_1←S_0跃迁.     

Nonlinear optical and structural properties of M@C N endohedrals (M = Li,Ca and Sc,N = 60 and 70)

We have investigated the static (ω = 0) and frequency-dependent nonlinear optical (NLO) properties of the M@C _N endohedrals (M = Li, Ca, Sc, N = 60 and 70) using the SSH (Su–Schrieffer–Heeger) approximation and sum-over-state (SOS) approach. Also, we study the effects of displacement and alkali, alkaline earth and lanthanide metal atoms and type of cage on the hyperpolarisabilities of the M@C _N endohedrals. The hyperpolarisability magnitudes and spectra are in agreement with experiment and the work of others using the SSH approximation. Our results indicate that the cage-type effect on the NLO spectra of M@C _N endohedrals is dramatic. Also, atom type has little effect on the highest peak value. These relationships between the atom and cage type and hyperpolarisability values may be beneficial to experimentalists when designing new NLO materials with large NLO responses.     

Analysis of multi-configuration density functional theory methods: theory and model application to bond-breaking

We consider the extension of the standard single-determinant Kohn–Sham method to the case of a multi-configuration auxiliary wave function. By applying the rigorous Kohn–Sham method to this case, we construct the proper interacting and auxiliary energy functionals. Following the Hohenberg–Kohn theorem for both energy functionals, we derive the corresponding multi-configuration Kohn–Sham equations, based on a local effective potential. At the end of the analysis we show that, at the ground state, the auxiliary wavefunction must collapse into a single-determinant wave function, equal to the regular KS wavefunction. We also discuss the stability of the wavefunction in multi-configuration density functional theory methods where the auxiliary system is partially interacting, and the remaining (residual) correlation is evaluated as a functional of the density. As an example showing both the challenges and the possibilities, we implement such a procedure for the perfect pairing wavefunction, using a residual correlation functional that is based on the Lee–Yang–Parr functional, and present results for an elementary bond-breaking process.     

A unified description of sum frequency generation,parametric down conversion and two-photon fluorescence

A superoperator non-equilibrium Green's function formalism is presented for computing nonlinear optical processes involving any combination of classical and quantum optical modes. Closed correlation-function expressions based on superoperator time-ordering are derived for the combined effects of causal response and non-causal spontaneous fluctuations. Coherent three wave mixing (sum frequency generation and parametric down conversion) involving one and two quantum optical modes, respectively, are compared with their incoherent counterparts: two-photon-induced fluorescence and two-photon-emitted fluorescence.     

Quantitative analysis of Fermi resonances by harmonic derivatives of perturbation theory corrections

Vibrational perturbation theory has proven to be a highly accurate and efficient method for extending the harmonic approximation in the treatment of polyatomic molecular vibrations. Unfortunately, accidental near-degeneracies of the harmonic vibrational levels can lead to resonance and a breakdown of the perturbation approximation. These resonances can be resolved by the diagonalization of a small effective Hamiltonian derived from either of the usual Rayleigh–Schrödinger or van Vleck perturbation theories. However, the proper choice of states for inclusion in the effective Hamiltonian is crucial to the accuracy of the results, and is not often clearly evident. It is proposed that the analytical partial derivatives of the anharmonic vibrational correction with respect to the various harmonic frequencies, called ‘Harmonic Derivatives’ in this work, can be used as a tool to quantitatively assess the existence and strength of first-order, or Fermi, resonances. These derivatives are shown to concisely and clearly reflect the quality of the perturbation approximation and the effect of its breakdown on the computed vibrational levels.