Electron correlation effects in the 4f14 shell

This paper serves a twofold purpose. First, Löwdin's inner projection in both nonperturbative and perturbative forms is applied to the quartic anharmonic oscillator. Inner projection with perturbation theory yields rational approximations to Brillouin–Wigner-type perturbation expansions. These lower bounds are compared with [N ? 1, N] Padé approximants to the Rayleigh–Schrödinger perturbation series for this problem. These Padés are also expressible as the even convergents, w_2N, of a Stieltjes-type continued fraction. The latter representation has certain advantages with respect to its Padé counterpart. Inner projection without perturbation theory provides significantly better results than the perturbative version. The application of inner projection techniques to a perturbed hydrogen atom is not straightforward. The usual problems associated with the continuum spectrum of hydrogen are present. By means of a nonunitary “tilting” transformation associated with the Lie group SO(4, 2), these problems may be bypassed. In the SO(4, 2)-reformulated eigenvalue problem, a reinterpretation of the basic variables, as developed by Silverstone and Moats, yields a new Hamiltonian that permits direct use of the inner projection method. This method has been applied to the ground state of the hydrogen atom in a magnetic field, using both four- and eight-dimensional basis manifolds. This represents the first application of inner projection to this problem.     

Aufbau principle and singlet-triplet gap in spherical Hooke atoms

Singlet and triplet spin state energies for three-dimensional Hooke atoms, that is, electrons in a quadratic confinement, with even number of electrons (2, 4, 6, 8, 10) is discussed using Full-CI and CASSCF type wavefunctions with a variety of basis sets and considering perturbative corrections up to second order. The effect of the screening of the electron–electron interaction is also discussed by using a Yukawa-type potential with different values of the Yukawa screening parameter (λ_ee = 0.2, 0.4, 0.6, 0.8, 1.0). Our results show that the singlet state is the ground state for two and eight electron Hooke atoms, whereas the triplet is the ground spin state for 4-, 6-, and 10-electron systems. This suggests the following Aufbau structure 1s < 1p < 1d with singlet ground spin states for systems in which the generation of the triplet implies an inter-shell one-electron promotion, and triplet ground states in cases when there is a partial filling of electrons of a given shell. It is also observed that the screening of electron–electron interactions has a sizable quantitative effect on the relative energies of both spin states, specially in the case of two- and eight-electron systems, favoring the singlet state over the triplet. However, the screening of the electron–electron interaction does not provoke a change in the nature of the ground spin state of these systems. By analyzing the different components of the energy, we have gained a deeper understanding of the effects of the kinetic, confinement and electron–electron interaction components of the energy.     

Improved Hylleraas calculations for ground state energies of lithium ISO–electronic sequence

The ²S ground state of lithium iso–electronic sequence is calculated by the use of Hylleraas-type wave functions. A 92 term one-spin wave function was used for lithium atom calculations. The energy obtained was ?7.478031 a.u. as compared with the previous best value of ?7.478025 a.u. calculated by Larsson. In addition, improved energies for Z = 4 to 8 were calculated by the use of 60 term wave functions. This work thus provides the lowest ab initio ground state energies for lithium sequence to date.     

Ab initio quasirelativistic calculations on electronic transitions in ICl by the multireference many‐body perturbation theory

A quasirelativistic perturbative method of ab initio calculations on ground and excited molecular electronic states and transition properties within the relativistic effective core potential approximation is presented and discussed. The method is based on the construction of a state‐selective many‐electron effective Hamiltonian in the model space spanned by an appropriate set of Slater determinants by means of the second‐order many‐body multireference perturbation theory. The neglect of effective spin–orbit interactions outside of the model space allows the exploitation of relatively high nonrelativistic symmetry during the evaluation of perturbative corrections and therefore dramatic reduction of the cost of computations without any contraction of the model‐space functions. One‐electron transition properties are evaluated via the perturbative construction of spin‐free transition density matrices. Illustrative calculations on the X0⁺ ? A1, B0⁺, and (ii)1 transitions in the ICl molecule are reported. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002     

Intruder state avoidance multireference Møller-Plesset perturbation theory

A new perturbation approach is proposed that enhances the low‐order, perturbative convergence by modifying the zeroth‐order Hamiltonian in a manner that enlarges any small‐energy denominators that may otherwise appear in the perturbative expansion. This intruder state avoidance (ISA) method can be used in conjunction with any perturbative approach, but is most applicable to cases where small energy denominators arise from orthogonal‐space states—so‐called intruder states—that should, under normal circumstances, make a negligible contribution to the target state of interests. This ISA method is used with multireference Møller–Plesset (MRMP) perturbation theory on potential energy curves that are otherwise plagued by singularities when treated with (conventional) MRMP; calculation are performed on the 1³Σ state of O₂; and the 2¹Δ, 3¹Δ, 2³Δ, and 3³Δ states of AgH. This approach is also applied to other calculations where MRMP is influenced by intruder states; calculations are performed on the ³Π_u state of N₂, the ³Π state of CO, and the 2¹A′ state of formamide. A number of calculations are also performed to illustrate that this approach has little or no effect on MRMP when intruder states are not present in perturbative calculations; vertical excitation energies are computed for the low‐lying states of N₂, C₂, CO, formamide, and benzene; the adiabatic ¹A₁–³B₁ energy separation in CH₂, and the spectroscopic parameters of O₂ are also calculated. Vertical excitation energies are also performed on the Q and B bands states of free‐base, chlorin, and zinc–chlorin porphyrin, where somewhat larger couplings exists, and—as anticipated—a larger deviation is found between MRMP and ISA‐MRMP. © 2002 Wiley Periodicals, Inc. J Comput Chem 10: 957–965, 2002     

Nearly Degenerate Isomers of C(BH)2: Cumulene,Carbene, or Carbone?

The ground electronic state of C(BH)₂ exhibits both a linear minimum and a peculiar angle‐deformation isomer with a central B‐C‐B angle near 90°. Definitive computations on these species and the intervening transition state have been executed by means of coupled‐cluster theory including single and double excitations (CCSD), perturbative triples (CCSD(T)), and full triples with perturbative quadruples (CCSDT(Q)), in concert with series of correlation‐consistent basis sets (cc‐pVXZ, X=D, T, Q, 5, 6; cc‐pCVXZ, X=T, Q). Final energies were pinpointed by focal‐point analyses (FPA) targeting the complete basis‐set limit of CCSDT(Q) theory with auxiliary core correlation, relativistic, and non‐Born–Oppenheimer corrections. Isomerization of the linear species to the bent form has a minuscule FPA reaction energy of 0.02 kcal mol^?1 and a corresponding barrier of only 1.89 kcal mol^?1. Quantum tunneling computations reveal interconversion of the two isomers on a timescale much less than 1 s even at 0 K. Highly accurate CCSD(T)/cc‐pVTZ and composite c～CCSDT(Q)/cc‐pCVQZ anharmonic vibrational frequencies confirm matrix‐isolation infrared bands previously assigned to linear C(BH)₂ and provide excellent predictions for the heretofore unobserved bent isomer. Chemical bonding in the C(BH)₂ species was exhaustively investigated by the atoms‐in‐molecules (AIM) approach, molecular orbital plots, various population analyses, local mode vibrations and force constants, unified reaction valley analysis (URVA), and other methods. Linear C(BH)₂ is a cumulene, whereas bent C(BH)₂ is best characterized as a carbene with little carbone character. Weak B–B attraction is clearly present in the unusual bent isomer, but its strength is insufficient to form a CB₂ ring with a genuine boron–boron bond and attendant AIM bond path.     

Enzymatic preparation of optically active silicon-containing amino acids and their application

Optically active 3-trimethyl silylalanine (TMS-Ala) was prepared by hydrolysis of N-acetyl-dl-TMS-Ala catalyzed by acylase I (aminoacylase; N-acylamino-acid amidohydrolase, EC3.5.1.14). Acylase I from porcine kidney (PKA) was found to be more effective than that from Aspergillus melleus in the preparation of l-TMS-Ala. Under the optimized conditions, optically pure l-TMS-Ala (>99% enantiomeric excess, ee) was obtained with a 72% yield. Furthermore, a highly optically pure d-TMS-Ala (96% ee) could also be obtained with a 76% yield by chemical hydrolysis of the residual substrate. Enzymatic synthesis of peptides containing TMS-Ala was also attempted in ethyl acetate. Benzyloxycarbonyl (Z)-l-TMS-Ala served as the substrate for thermolysin, whereas l-TMS-Ala-OMe was inactive as the amino component. In the case of inhibitory activity of dipeptides toward thermolysin, l-Leu-(l-TMS-Ala) was found to be a more potent inhibitor than l-Leu-l-Leu, which is known to be one of the most effective inhibitors of thermolysin among the dipeptides consisting of natural aminoacids.     

Theoretical Study on the Formation of H‐ and O‐Atoms,HONO, OH,NO, and NO2 from the Lowest Lying Singlet and Triplet States in Ortho‐Nitrophenol Photolysis

The photolysis of nitrophenols was proposed as a source of reactive radicals and NOx compounds in polluted air. The S₀ singlet ground state and T₁ first excited triplet state of nitrophenol were investigated to assess the energy dependence of the photofragmentation product distribution as a function of the reaction conditions, based on quantum chemical calculations at the G3SX//M06–2X/aug‐cc‐pVTZ level of theory combined with RRKM master equation calculations. On both potential energy surfaces, we find rapid isomerization with the aci‐nitrophenol isomer, as well as pathways forming NO, NO₂, OH, HONO, and H‐, and O‐atoms, extending earlier studies on the T₁ state and in agreement with available work on other nitroaromatics. We find that accessing the lowest photofragmentation channel from the S₀ ground state requires only 268 kJ/mol of activation energy, but at a pressure of 1 atm collisional energy loss dominates such that significant fragmentation only occurs at internal energies exceeding 550 kJ/mol, making this surface unimportant for atmospheric photolysis. Intersystem crossing to the T₁ triplet state leads more readily to fragmentation, with dissociation occurring at energies of ～450 kJ/mol above the singlet ground state even at 1 atm. The main product is found to be OH + nitrosophenoxy, followed by formation of hydroxyphenoxy + NO and phenyloxyl + HONO. The predictions are compared against available experimental data.     

Relativistic effects in HgHe and HgXe CCSD(T) ground state potential curves. Low‐density viscosity simulations of Hg:Xe mixture

The comparison of coupled cluster with single and double excitations and with perturbative correction of triple excitations [CCSD(T)] ground state potential curves of mercury with rare gases (RG): HgHe and HgXe, at several levels of theory is presented. The scalar relativistic (REL) effects and spin‐orbit coupling effects in the ground state potential curves of these weakly bounded dimers are considered. The CCSD(T) ground state potential curves at the level of the Dirac‐Coulomb Hamiltonian (DCH) are compared with CCSD(T) curves at the level of 4‐component spin‐free modified DCH, the scalar 2nd order Douglas‐Kroll‐Hess (DKH2) and the nonrelativistic (NR‐LL) (Lévy‐Leblond) Hamiltonian. In addition, London‐Drude formula and SCF interaction energy curves are employed in the analysis of different contributions of REL effects in dissociation energies of HgRG and Hg₂ dimers. Moreover, the large anharmonicity of the HgHe ground state potential curve is highlighted. The computationally less demanding scalar DKH2 Hamiltonian is employed to calculate the HgXe, Hg₂, and Xe₂ all electron CCSD(T) ground state potential curves in highly augmented quadruple zeta basis sets. These potential curves are used to simulate the shear viscosity of mercury, xenon, and mercury‐xenon (Hg:Xe) mixture. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2011     

Excited state properties,fluorescence energies,and lifetimes of a poly(fluorene‐phenylene), based on TD‐DFT investigation

The structural and electronic properties of fluorene‐phenylene copolymer (FP)_n, n = 1–4 were studied by means of quantum chemical calculations based on density functional theory (DFT) and time dependent density functional theory (TD‐DFT) using B3LYP functional. Geometry optimizations of these oligomers were performed for the ground state and the lowest singlet excited state. It was found that (FP)_n is nonplanar in its ground state while the electronic excitations lead to planarity in its S₁ state. Absorption and fluorescence energies were calculated using TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods. Vertical excitation energies and fluorescence energies were obtained by extrapolating these values to infinite chain length, resulting in extrapolated values for vertical excitation energy of 2.89 and 2.87 eV, respectively. The S₁ ← S₀ electronic excitation is characterized as a highest occupied molecular orbital to lowest unoccupied molecular orbital transition and is distinguishing in terms of oscillator strength. Fluorescence energies of (FP)_n calculated from TD‐B3LYP/SVP and TD‐B3LYP/SVP+ methods are 2.27 and 2.26 eV, respectively. Radiative lifetimes are predicted to be 0.55 and 0.51 ns for TD‐B3LYP/SVP and TD‐B3LYP/SVP+ calculations, respectively. These fundamental information are valuable data in designing and making of promising materials for LED materials. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2010     

Time‐dependent density functional calculations on the electronic spectra of the neutral nickel complex [Ni(LISQ)2] (LISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1−)) and its monoanion and dication

The electronic spectrum of the neutral nickel complex [Ni(L^ISQ)₂] (L^ISQ = 3,5‐di‐tert‐butyl‐o‐diiminobenzosemiquinonate(1^?)) and the spectra of its anion and dication have been calculated by means of time‐dependent density functional theory. The electronic ground state of the neutral complex exhibits an open shell singlet diradical character. The mandatory multireference problem for this electronic ground state has been treated approximately by using the unrestricted and spin symmetry broken Kohn‐Sham Slater determinant as the wave function for the noninteracting reference system in the time‐dependent density functional calculations. A reasonable agreement with observed transition energies and band intensities has been achieved. This holds also for the long wavelength transitions that are shown to be of charge transfer type. The charge distributions in the electronic ground state and the corresponding low lying excited states, however, are rather similar. Thus, the known failure of standard time‐dependent density functional theory to describe improperly long range charge transfer transitions is absent in this work. The applied computational scheme might be adequate for calculating electronic spectra of transition metal complexes with noninnocent ligands. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009     

Six-dimensional tunneling dynamics of internal rotation in hydrogen peroxide molecule and its isotopomers

The six-dimensional torsion-vibration Hamiltonian of the H₂O₂ molecule and its H/D- and ¹⁸O/¹⁶O-isotopomers is derived. The Hamiltonian includes the kinetic energy operator, which depends on the tunneling coordinate, and the potential energy surface represented as a quartic polynomial with respect to the small-amplitude transverse coordinates. Parameters of the Hamiltonian were obtained from DFT calculations of the equilibrium geometries, eigenvectors, and eigenfrequencies of normal vibrations at the stationary points corresponding to the ground state and both the cis- and trans-transition states, carried out with the B3LYP density functional and 6-311+G(2d,p) basis set. The quantum dynamics problem is solved using the perturbative instanton approach generalized for the excited states situated above the barrier top. Vibration-tunneling spectra are calculated for the ground state and low-lying excited states with energies below 2000 cm^–1. Strong kinematic and squeezed potential couplings between the large-amplitude torsional motion and bending modes are shown to be responsible for the vibration-assisted tunneling and for the dependence of tunneling splittings on the quantum numbers of small-amplitude transverse vibrations. Mode-specific isotope effects are predicted.     

Ab initio characterization of size dependence of electronic spectra for linear anionic carbon clusters C(n) (-) (n = 4-17)

In this article, we determine the ground‐state equilibrium geometries of the linear anionic carbon clusters C (n = 4–17) by means of the density functional theory B3LYP, CAM‐B3LYP, and coupled cluster CCSD(T) calculations, as well as their electronic spectra obtained by the multireference second‐order perturbation theory CASPT2 method. These studies indicate that these linear anions possess doublet ²∏_g or ²∏_u ground state, and the even‐numbered clusters are generally acetylenic, whereas the odd‐numbered ones are essentially cumulenic. The energy differences, electron affinities, and incremental binding energies of C chains all exhibit a notable tread of parity alternation, with n‐even chains being more stable than n‐odd ones. In addition, the predicted vertical excitation energies from the ground state to four low‐lying excited states are in reasonably good agreement with the available experimental observations, and the calculations for the higher excited electronic transitions can provide accurate information for the experimentalists and spectroscopists. Interestingly, the absorption wavelengths of the 1²∏_u/g ← X²∏_g/u transitions of the n‐even clusters show a nonlinear trend of exponential growth, whereas those of the n‐odd counterparts are found to obey a linear relationship as a function of the chain size, as shown experimentally. Moreover, the absorption wavelengths of the transitions to the higher excited states of C series have the similar linear size dependence as well. © 2011 Wiley Periodicals, Inc. J Comput Chem, 2011     

The Evolution of Bonding and Thermodynamic Properties of Boron‐Doped Small Carbon Clusters: An Ab Initio Study

Theoretical studies on BC_n (n=1–6) clusters are carried out using density functional theory, Møller–Plesset second‐order perturbation theory (MP2), coupled‐cluster calculations including up to triple excitations (CCSD(T)), and higher‐level approaches. All possible isomers depending on the positions of the boron atom are generated and the lowest‐energy isomers are determined for doublet and quartet electronic states. The three potential evolution paths of the clusters are determined as a function of their size. The energetic and electronic consequences for the increased size of structures differ significantly, which leads to representatives of the ground electronic state from different structural groups. The ab initio calculated thermal functions allow enhancements to the available atomization energies and improve the agreement between the calculated and experimental heat content.     

Low-temperature heat capacity of heptacopper(II) complex [Cu7(μ3-Cl)2(μ3-OH)6-(d-pen-disulfide)3]

A single crystal calorimetry of a heptacopper(II) complex of [Cu₇(μ₃-Cl)₂(μ₃-OH)₆-(d-pen-disulfide)₃] which has a double-cubane structure supported by d-penicillaminedisulfides has been performed at low-temperature region below 8 K. This compound is a metal complex which contains seven Cu(II)s in a cluster unit. These Cu(II)s are magnetically coupled each other by strong intra-complex interactions. The heat capacities under magnetic fields exhibit Schottky type anomalies explained by the Zeeman splitting of the doublet ground state of the complex. The g-value of the ground state is evaluated as 1.86 from the systematic analysis of the Schottky peak under magnetic fields. The first excited state of the cluster seems to be separated at least by several Kelvins, which is consistent with the theoretical calculations and magnetic susceptibility results.     

Conical Intersections and Low‐Lying Electronic States of Tetrafluoroethylene

The low‐lying electronic states of tetrafluoroethylene (C₂F₄) are characterized theoretically for the first time using equation‐of‐motion coupled cluster theory (EOM‐CCSD), and complete active space self‐consistent field (CASSCF) and second‐order perturbation theory (CASPT2). Computations are performed for vertical excitation energies, equilibrium geometries, minimum‐energy conical intersections, and potential energy curves along three geometric coordinates: 1) twisting of the F?C?C?F dihedral angle, 2) pyramidalization of the CF₂ group, and 3) migration of a fluorine atom resulting in an ethylidene‐like (CF₃CF) structure. The results suggest two relaxation pathways from the Rydberg‐3s excited electronic state to the ground state. These relaxation pathways are discussed in conjunction with the femtosecond photoionization spectroscopy results of Trushin et al. [ChemPhysChem­ 2004 , 5, 1389].     

Multireference configuration interaction calculations for the F(2P)+HCl-->HF+Cl(2P) reaction: a correlation scaled ground state (1 2A') potential energy surface

This paper presents a new ground state (1 (2)A(')) electronic potential energy surface for the F((2)P)+HCl-->HF+Cl((2)P) reaction. The ab initio calculations are done at the multireference configuration interaction+Davidson correction (MRCI+Q) level of theory by complete basis set extrapolation of the aug-cc-pVnZ (n=2,3,4) energies. Due to low-lying charge transfer states in the transition state region, the molecular orbitals are obtained by six-state dynamically weighted multichannel self-consistent field methods. Additional perturbative refinement of the energies is achieved by implementing simple one-parameter correlation energy scaling to reproduce the experimental exothermicity (DeltaE=-33.06 kcalmol) for the reaction. Ab initio points are fitted to an analytical function based on sum of two- and three-body contributions, yielding a rms deviation of <0.3 kcalmol for all geometries below 10 kcalmol above the barrier. Of particular relevance to nonadiabatic dynamics, the calculations show significant multireference character in the transition state region, which is located 3.8 kcalmol with respect to F+HCl reactants and features a strongly bent F-H-Cl transition state geometry (theta approximately 123.5 degrees ). Finally, the surface also exhibits two conical intersection seams that are energetically accessible at low collision energies. These seams arise naturally from allowed crossings in the C(infinityv) linear configuration that become avoided in C(s) bent configurations of both the reactant and product, and should be a hallmark of all X-H-Y atom transfer reaction dynamics between ((2)P) halogen atoms.     

The lowest electronic states of the benzoyl ion. A molecular orbital study with configuration interaction

A configuration interaction method of molecular orbital theory shown to be accurate for the calculation of excited state energies in several aromatic systems was applied to the problem of excited states of the benzoyl ion. The excited singlet and triplet states of the benzoyl ion lie at least 3.8 eV and 2.6 eV, respectively, above the ground state. These results are not in agreement with a postulated state 20 kcal above the ground state. On the other hand, the charge distribution in excited states does agree with that postulated for the 20 kcal state. The p-hydroxy and p-cyano substituents do not greatly influence the charge distribution between the ring and the carbonyl group in either the ground or lowest excited states.     

Ab initio description of photoabsorption and electron transfer in a doubly-linked porphyrin-fullerene dyad

Structure, photoabsorption and excited states of two representative conformations obtained from molecular dynamics (MD) simulations of a doubly-linked porphyrin-fullerene dyad DHD6ee are studied by using both DFT and wavefunction based methods. Charge transfer from the donor (porphyrin) to the acceptor (fullerene) and the relaxation of the excited state are of special interest. The results obtained with LDA, GGA, and hybrid functionals (SVWN, PBE, and B3LYP, respectively) are analyzed with emphasis on the performance of used functionals as well as from the point of view of their comparison with wavefunction based methods (CCS, CIS(D), and CC2). Characteristics of the MD structures are retained in DFT optimization. The relative orientation of porphyrin and fullerene is significantly influencing the MO energies, the charge transfer (CT) in the ground state of the dyad and the excitation of ground state CT complex (g-CTC). At the same time, the excitation to the locally excited state of porphyrin is only little influenced by the orientation or cc distance. TD-DFT underestimates the excitation energy of the CT state, however for some cases (with relatively short donor-acceptor separations), the use of a hybrid functional like B3LYP alleviates the problem. Wavefunction based methods and CC2 in particular appear to overestimate the CT excitation energies but the inclusion of proper solvation models can significantly improve the results.