Search

The thermochemistry of the Si-O-H system has been extensively studied both experimentally and theoretically due to its importance in chemical processes, degradation of silica-protected materials in combustion, and geological processes. In this paper, we review past studies and use quantum mechanical methods to generate a new data set. Molecular geometries were generated with DFT using the B3LYP functional. Energetics were calculated with RCCSD(T) methods extrapolated to the complete basis set (CBS/45) limit. Particular attention was given to the treatment of the vibrational modes. A rigid rotor model was used, corrections for anharmonicity were applied, and the Pitzer-Gwinn treatment of the hindered rotation of the M-OH groups was applied. The generated enthalpies of formation at 298 K are compared to those of experiments and other calculations. Generally, the agreement is good. A set of thermodynamic data (enthalpy of formation at 298 K, entropy at 298 K, and heat capacity polynomial to 3000 K) is presented for SiOH, SiO(OH), Si(OH) 2 , SiO(OH) 2 , Si(OH) 3 , Si(OH) 4 , Si 2 O(OH) 6 , and Si 3 O 2 (OH) 8 . These can be added to any of the common computational thermodynamics packages. The application of these data to high-temperature corrosion and geological problems is discussed.

We employ density functional theory (DFT) to examine reaction mechanisms involving singlet oxygen 1 Δ g ( 1 O 2 ) and 1,2-dimethoxyethane (DME) to probe potential parasitic reactions occurring in Li-O 2 batteries. First, we investigate the attack of 1 O 2 on the ethylene group (-CH 2 -CH 2 -) to form H 2 O 2 and a C-C double bond in a single step. Second, we look at hydroperoxide formation that occurs via a two-step mechanism. We employ an implicit solvent model, Li + coordination, and external electric fields to model the complex electrolyte environment near the cathode of a Li-O 2 battery. The initial barriers for these reactions are decreasing functions of the dielectric constant of the implicit solvent model as well as the strength of the electric field. These initial barriers range between 17 and 26 kcal mol -1 for large dielectric constants and in the presence of electric fields. We discuss the implications of these results on ether-based electrolytes for Li-O 2 batteries.

Hafnium dioxide vaporizes primarily to HfO(g) in a reducing environment. The thermochemistry of HfO(g) is calculated from quantum methods and measured via Knudsen effusion mass spectrometry. For the computations, all-electron and relativistic effective core potential calculations are used. The calculation of an accurate dissociation energy and an entire potential energy curve is reported. These calculations lead to Δ f H°(298) = 63.19 ± 10 kJ/mol, S°(298) = 235.52 J/mol K, and C p (298-2500 K) = (2.741 × 10 -9 )T 3 - (9.853 × 10 -6 )T 2 + (1.295 × 10 -2 )T + 2.761 × 10 -1 J/mol K. Experimentally, HfO(g) is generated from the reaction of Hf(s) and HfO 2 (s) in a specially made Hf Knudsen cell. A third law treatment of the data leads to Δ f H°(298) of 58.4 ± 12.3 kJ/mol, in good agreement with the calculated value.

Thermodynamic quantities are calculated for gaseous hydroxides and oxyhydroxides of Cr, Mn, and La. These would form due to water-vapor-containing environments reacting with Cr-forming alloys or oxide components of potential fuel cell interconnects or anode materials. Structures and vibrational modes for the expected hydroxides and oxyhydroxides are calculated with the B3LYP hybrid functional. Enthalpies of formation from selected reactions for each species are calculated using the CCSD(T)/CBS approach. Results show good agreement with literature estimates, measurements, and calculations. The resultant data is reported as Δ f H °(298), S °(298), and C p ( T ) and put into the database for a free-energy minimizer code. Calculations are presented to show the hydroxide and oxyhydroxide vapor pressures above H 2 O + Cr 2 O 3 , Mn 3 O 4 , and La 2 O 3 , as well as the anode material La 0.75 Sr 0.25 Cr 0.5 Mn 0.5 O 3-δ (LSCM).

Recent experimental and computational evidence indicates that singlet oxygen ( 1 O 2 ) attacks the ethylene group (-CH 2 -CH 2 -) in ethylene carbonate (EC) leading to degradation in Li-ion batteries employing EC as the electrolyte solvent [ J. Phys. Chem. A 2018 , 122 , 8828-8839]. Here, we employ computational quantum chemistry to explore this mechanism in detail for a large set of organic molecules. Benchmark calculations comparing density functional theory to the complete active space second-order perturbation theory and internally contracted multireference configuration interaction indicate that the M11 functional adequately captures trends in the transition-state energies for this mechanism. Based on our results, we recommend that solvents which include the ethylene group should be avoided in Li-ion and Li-O 2 batteries where 1 O 2 is generated unless neighboring functional groups raise the reaction barrier to avoid this decomposition pathway.

Gd 2 O 3 and Yb 2 O 3 are the proposed constituents of advanced coating systems in combustion environments. In such environments, they are exposed to high-temperature water vapor, which would lead to gaseous hydroxide formation. Thermodynamic parameters are reported for YbO n (OH) m and GdO n (OH) m species. We first study the MH, MO, MF, and MCl (M = Yb and Gd) species, where some experimental data exist. Structures and spectroscopic constants were calculated at the B3LYP level. For YbO n (OH) m , the B3LYP approach is used in conjunction with an effective core potential, while for GdO n (OH) m , it was necessary to use all-electron basis sets. Enthalpies of formation were calculated with a BD(T) approach. The enthalpies of formation, entropies, and heat capacities were added to a thermochemical database. Hydroxide and oxyhydroxide vapor pressures are calculated above pure Yb 2 O 3 , Gd 2 O 3 , and Y 2 O 3 in 50% H 2 O/Ar from 1000 to 3000 K. Hydroxide and oxyhydroxide vapor pressures are also calculated for potential coating compositions.

The abstraction of a proton by OH - , O 2 - , and XO 2 - from DME n ···X + , where X is Li, Na, or K, is studied using density functional theory. Both the gas phase and the solution phase are studied. In general, when explicit solvent molecules are added, the difference between the gas-phase and solution results becomes rather small. While the DME n ···X + binding energies differ significantly for various alkali cations, the reaction energies and transition-state energies are far less sensitive to the choice of an alkali cation. XO 2 - has a lower barrier height than OH - , which, in turn, has a lower barrier height than O 2 - . The reaction energies follow the same trends.

The Li + -ligand binding energies are computed for seven ligands and their perfluoro analogs using Density Functional Theory. The bonding is mostly electrostatic in origin. Thus the size of the binding energy tends to correlate with the ligand dipole moment, however, the charge-induced dipole contribution can be sufficiently large to affect the dipolebinding energy correlation. The perfluoro species are significantly less strongly bound than their parents, because the electron withdrawing power of the fluorine reduces the ligand dipole moment.

We have investigated the mid-infrared spectral characteristics of a series of polycyclic aromatic hydrocarbons (PAHs) with straight edges and containing an even or odd number of carbons using density functional theory (DFT). For several even and odd-carbon PAHs, the 8.6/6.2 and 7.6/6.2 intensity ratios computed in emission after the absorption of a 8 eV photon match the observed ratios obtained for three reflection nebulae (RNe), namely NGC 1333, NGC 7023, and NGC 2023. Odd-carbon PAHs are favored, particularly for NGC 1333. Both cations and anions are present with the cations being predominant. Relevant PAHs span sizes ranging from 46 to 103-113 carbons for NGC 7023 and NGC 2023 and from 38 to 127 carbons for NGC 1333 and have symmetries ranging from D 2 h to C s . Our work suggests that even and odd-carbon PAHs with straight edges are viable candidates for the PAH emission seen towards irradiated Photo-Dissociation Regions (PDRs).

Gaseous titanium hydroxide and oxyhydroxide species were studied with quantum chemical methods. The results are used in conjunction with an experimental transpiration study of titanium dioxide (TiO 2 ) in water vapor-containing environments at elevated temperatures to provide a thermodynamic description of the Ti(OH) 4 (g) and TiO(OH) 2 (g) species. The geometry and harmonic vibrational frequencies of these species were computed using the coupled-cluster singles and doubles method with a perturbative correction for connected triple substitutions [CCSD(T)]. For the OH bending and rotation, the B3LYP density functional theory was used to compute corrections to the harmonic approximations. These results were combined to determine the enthalpy of formation. Experimentally, the transpiration method was used with water contents from 0 to 76 mol % in oxygen or argon carrier gases for 20-250 h exposure times at 1473-1673 K. Results indicate that oxygen is not a key contributor to volatilization, and the primary reaction for volatilization in this temperature range is TiO 2 (s) + H 2 O(g) = TiO(OH) 2 (g). Data were analyzed with both the second and third law methods using the thermal functions derived from the theoretical calculations. The third law enthalpy of formation at 298.15 K for TiO(OH) 2 (g) at 298 K was -838.9 ± 6.5 kJ/mol, which compares favorably to the theoretical calculation of -838.7 ± 25 kJ/mol. We recommend the experimentally derived third law enthalpy of formation at 298.15 K for TiO(OH) 2 , the computed entropy of 320.67 J/mol·K, and the computed heat capacity [149.192 + (-0.02539)T + (8.28697 × 10 -6 )T 2 + (-15614.05)/T + (-5.2182 × 10 -11 )/T 2 ] J/mol-K, where T is the temperature in K.

The low-lying singlet and triplet states of MgO have been studied using a SA-CASCF/IC-MRCI approach using the aug-cc-pV5Z basis set. The spectroscopic constants ( r e , w e , and T e ) are in good agreement with the available experimental data. The computed lifetime for the B state is in excellent agreement with two of the three experimental results. The d state lifetime is in good agreement with experiment, while the computed D state lifetime is about twice as long as experiment.

The matrix-isolated, mid-infrared spectra of seven acridine-based polycyclic aromatic nitrogen heterocycles (PANHs) have been measured and compared to their non-nitrogen containing parent molecule. The acridine species investigated include acridine, benz[a]acridine, benz[c]acridine, dibenz[a,j]acridine, dibenz[c,h]acridine, dibenz[a,h]acridine and dibenz[a,c]acridine. The previously reported results for 1 and 2-azabenz[a]anthracenes are included for comparison. The experimentally determined band frequencies and intensities are compared with their B3LYP/6-31G(d) values. The overall agreement between experimental and theoretical values is good and in line with our previous investigations. Shifts, typically to the blue, are noted for the C-H out-of-plane (CH oop ) motions upon insertion of a nitrogen atom. The formation of a bay region upon addition of additional benzene rings to the anthracene/acridine structure splits the solo hydrogen motions into a bay region solo and an external solo hydrogen, with the bay region solo hydrogen coupling to the quartet hydrogen motions and the external solo hydrogen coupling with the duo hydrogen motions resulting in an extreme decrease in intensity for the CH oop solo hydrogen band when the external hydrogen is replaced by a nitrogen atom. The C-C and C-H in-plane region of this acridine series exhibits the characteristic two fold increase in intensity, noted previously for PANHs. The strong ≈1400cm -1 band, which was identified in the previous PANH study, is noted in several molecular species as well as another strong PANH feature between 1480 and 1515cm -1 for several molecules. The presence of these strong bands appear to be primarily responsible for the two-fold increase in the C-H in-plane region's (1100-1600cm -1 ) intensity. The C-H stretching region can be characterized by contributions from the solo (bay or external), duo and quartet hydrogens, similar to what was observed in the dibenzopolyacene compounds., (Published by Elsevier B.V.)

Interactions between pre-cured phenolic polymer chains and a solvent have a significant impact on the structure and properties of the final postcured phenolic resin. Developing an understanding of the nature of these interactions is important and will aid in the selection of the proper solvent that will lead to the desired final product. Here, we investigate the role of the phenolic chain structure and the solvent type on the overall solvation performance of the system through molecular dynamics simulations. Two types of solvents are considered: ethylene glycol (EGL) and H 2 O. In addition, three phenolic chain structures are considered, including two novolac-type chains with either an ortho-ortho (OON) or an ortho-para (OPN) backbone network and a resole-type (RES) chain with an ortho-ortho network. Each system is characterized through a structural analysis of the solvation shell and the hydrogen-bonding environment as well as through a quantification of the solvation free energy along with partitioned interaction energies between specific molecular species. The combination of simulations and the analyses indicate that EGL provides a higher solvation free energy than H 2 O due to more energetically favorable hydrophilic interactions as well as favorable hydrophobic interactions between CH element groups. In addition, the phenolic chain structure significantly affects the solvation performance, with OON having limited intermolecular hydrogen-bond formations, while OPN and RES interact more favorably with the solvent molecules. The results suggest that a resole-type phenolic chain with an ortho-para network should have the best solvation performance in EGL, H 2 O, and other similar solvents.

Ab initio techniques are used to study the interaction of ethylene glycol and water with a phenolic polymer. The water bonds more strongly with the phenolic OH than with the ring. The phenolic OH groups can form hydrogen bonds between themselves. For more than one water molecule, there is a competition between water-water and water-phenolic interactions. Ethylene glycol shows the same effects as those of water, but the potential energy surface is further complicated by CH 2 -phenolic interactions, different conformers of ethylene glycol, and two OH groups on each molecule. Thus, the ethylene glycol-phenolic potential is more complicated than the water-phenolic potential. The results of the ab initio calculations are compared to those obtained using a force field. These calibration studies show that the water system is easier to describe than the ethylene glycol system. The calibration studies confirm the reliability of force fields used in our companion molecular dynamics study of a phenolic polymer in water and ethylene solutions.

Density functional theory (DFT), density functional theory molecular dynamics (DFT-MD), and classical molecular dynamics using polarizable force fields (PFF-MD) are employed to evaluate the influence of Li(+) on the structure, transport, and electrochemical stability of three potential ionic liquid electrolytes: N-methyl-N-butylpyrrolidinium bis(trifluoromethanesulfonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosulfonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]). We characterize the Li(+) solvation shell through DFT computations of [Li(Anion)n]((n-1)-) clusters, DFT-MD simulations of isolated Li(+) in small ionic liquid systems, and PFF-MD simulations with high Li-doping levels in large ionic liquid systems. At low levels of Li-salt doping, highly stable solvation shells having two to three anions are seen in both [pyr14][TFSI] and [pyr13][FSI], whereas solvation shells with four anions dominate in [EMIM][BF4]. At higher levels of doping, we find the formation of complex Li-network structures that increase the frequency of four anion-coordinated solvation shells. A comparison of computational and experimental Raman spectra for a wide range of [Li(Anion)n]((n-1)-) clusters shows that our proposed structures are consistent with experiment. We then compute the ion diffusion coefficients and find measures from small-cell DFT-MD simulations to be the correct order of magnitude, but influenced by small system size and short simulation length. Correcting for these errors with complementary PFF-MD simulations, we find DFT-MD measures to be in close agreement with experiment. Finally, we compute electrochemical windows from DFT computations on isolated ions, interacting cation/anion pairs, and liquid-phase systems with Li-doping. For the molecular-level computations, we generally find the difference between ionization energy and electron affinity from isolated ions and interacting cation/anion pairs to provide upper and lower bounds, respectively, to experiment. In the liquid phase, we find the difference between the lowest unoccupied and highest occupied electronic levels in pure and hybrid functionals to provide lower and upper bounds, respectively, to experiment. Li-doping in the liquid-phase systems results in electrochemical windows little changed from the neat systems.

We employ molecular dynamics (MD) simulation and experiment to investigate the structure, thermodynamics, and transport of N-methyl-N-butylpyrrolidinium bis(trifluoromethylsufonyl)imide ([pyr14][TFSI]), N-methyl-N-propylpyrrolidinium bis(fluorosufonyl)imide ([pyr13][FSI]), and 1-ethyl-3-methylimidazolium boron tetrafluoride ([EMIM][BF4]), as a function of Li-salt mole fraction (0.05 ≤ xLi(+) ≤ 0.33) and temperature (298 K ≤ T ≤ 393 K). Structurally, Li(+) is shown to be solvated by three anion neighbors in [pyr14][TFSI] and four anion neighbors in both [pyr13][FSI] and [EMIM][BF4], and at all levels of xLi(+) we find the presence of lithium aggregates. Pulsed field gradient spin-echo NMR measurements of diffusion and electrochemical impedance spectroscopy measurements of ionic conductivity are made for the neat ionic liquids as well as 0.5 molal solutions of Li-salt in the ionic liquids. Bulk ionic liquid properties (density, diffusion, viscosity, and ionic conductivity) are obtained with MD simulations and show excellent agreement with experiment. While the diffusion exhibits a systematic decrease with increasing xLi(+), the contribution of Li(+) to ionic conductivity increases until reaching a saturation doping level of xLi(+) = 0.10. Comparatively, the Li(+) conductivity of [pyr14][TFSI] is an order of magnitude lower than that of the other liquids, which range between 0.1 and 0.3 mS/cm. Our transport results also demonstrate the necessity of long MD simulation runs (∼200 ns) to converge transport properties at room temperature. The differences in Li(+) transport are reflected in the residence times of Li(+) with the anions (τ(Li/-)), which are revealed to be much larger for [pyr14][TFSI] (up to 100 ns at the highest doping levels) than in either [EMIM][BF4] or [pyr13][FSI]. Finally, to comment on the relative kinetics of Li(+) transport in each liquid, we find that while the net motion of Li(+) with its solvation shell (vehicular) significantly contributes to net diffusion in all liquids, the importance of transport through anion exchange increases at high xLi(+) and in liquids with large anions.

We report experimental spectra in the mid-infrared (IR) and near-IR for a series of dibenzoacenes isolated in Ar matrices. The experiments are supported by Density Functional Theory (DFT) and Time-Dependent DFT (TD-DFT) calculations with both vibrational and electronic transitions studied. For the neutrals, we find good agreement between the experimental and B3LYP and BP86 results for all species studied. The band at about 1440 cm(-1) carries more intensity than in typical PAHs and increases in intensity with the size of the dibenzoacene molecule. For the ions the B3LYP approach fails to yield reasonable IR spectra for most systems and the BP86 approach is used. Electronic transitions dominate the vibrational bands in the mid-IR region for the large dibenzoacene ions. In spite of the very strong electronic transitions, there is still reasonable agreement between theory and experiment for the vibrational band positions. The experimental and theoretical results for the dibenzoacenes are also compared with those for the polyacenes., (Published by Elsevier B.V.)

The Li(+)-BF4(-) and BF4(-)-BF4(-) interactions are studied using second order perturbation theory (MP2) and coupled cluster singles and doubles approach, including the effect of connected triples, CCSD(T). The MP2 and CCSD(T) results are in excellent agreement. Using only the MP2 approach, the interactions of Li(+) with bis(trifluoromethane)sulfonimide anion (TFSI) and Li(+) with bis(fluorosulfonyl)imide anion (FSI) are studied. The results of these high level calculations are compared with density functional theory (DFT) calculations for a variety of functionals and with the APPLE&P force field. The B3LYP approach well reproduces the accurate calculations using both a small and large basis set. The M06 and M06L functionals in the larger basis set are in good agreement with the high level calculations. While the APPLE&P force field does not outperform the best functionals, the APPLE&P results agree better with the accurate results than do some of the functionals tested.

Reaction paths for the loss of CO, H2, and H2O from atomistic models of phenolic resin are determined using the hybrid B3LYP approach. B3LYP energetics are confirmed using CCSD(T). The energetics along the B3LYP paths are also evaluated using the PW91 generalized gradient approximation (GGA), the more approximate self-consistent charge density functional tight binding (SCC-DFTB), and the reactive force field (ReaxFF). Compared with the CCSD(T)/cc-pVTZ level for bond and reaction energies and barrier heights, the B3LYP, PW91, DFTB(mio), DFTB(pbc), and ReaxFF have average absolute errors of 3.8, 5.1, 17.4, 13.2, and 19.6 kcal/mol, respectively. The PW91 is only slightly less accurate than the B3LYP approach, while the more approximate approaches yield somewhat larger errors. The SCC-DFTB paths are in better agreement with B3LYP than are those obtained with ReaxFF.

A systematic comparison of atomistic modeling methods including density functional theory (DFT), the self-consistent charge density-functional tight-binding (SCC-DFTB), and ReaxFF is presented for simulating the initial stages of phenolic polymer pyrolysis. A phenolic polymer system is simulated for several hundred picoseconds within a temperature range of 2500 to 3500 K. The time evolution of major pyrolysis products including small-molecule species and char is examined. Two temperature zones are observed which demark cross-linking versus fragmentation. The dominant chemical products for all methods are similar, but the yields for each product differ. At 3500 K, DFTB overestimates CO production (300-400%) and underestimates free H (~30%) and small C(m)H(n)O molecules (~70%) compared with DFT. At 3500 K, ReaxFF underestimates free H (~60%) and fused carbon rings (~70%) relative to DFT. Heterocyclic oxygen-containing five- and six-membered carbon rings are observed at 2500 K. Formation mechanisms for H2O, CO, and char are discussed. Additional calculations using a semiclassical method for incorporating quantum nuclear energies of molecules were also performed. These results suggest that chemical equilibrium can be affected by quantum nuclear effects at temperatures of 2500 K and below. Pyrolysis reaction mechanisms and energetics are examined in detail in a companion manuscript.

The oxidifluoride molecules, OYF(2) and OLaF(2), are produced via the reactions of laser ablated metal atoms with OF(2) in solid argon. The product structures are characterized using matrix isolation infrared spectroscopy as well as theoretical calculations. Similar to the very recently characterized OScF(2) molecule, OYF(2) is predicted to have a (2)B(2) ground state with C(2v) symmetry while the heavier OLaF(2) has a (2)A″ ground state with near C(2v) symmetry. The unpaired electron is mainly located on the terminal oxygen atom, suggesting radical character for the group 3 OMF(2) molecules. In addition, the closed shell singlet OMF molecules with bent geometries are also observed, and they are found to have triple metal-oxygen bonds with higher stretching frequencies and shorter bond lengths than their OMF(2) counterparts. α-Fluorine transfer from OF(2) to metal centers is predicted to be highly exothermic, which is very favorable for the formation of new OMF(2) and OMF species.

The isolated group 4 metal oxydifluoride molecules OMF(2) (M = Ti, Zr, Hf) with terminal oxo groups are produced specifically on the spontaneous reactions of metal atoms with OF(2) through annealing in solid argon. The product structures and vibrational spectra are characterized using matrix isolation infrared spectroscopy as well as B3LYP density functional and CCSD(T) frequency calculations. OTiF(2) is predicted to have a planar structure while both OZrF(2) and OHfF(2) possess pyramidal structures, all with singlet ground states. Three infrared absorptions are observed for each product molecule, one M-O and two M-F stretching modes, and assignments of these molecules are further supported by the corresponding (18)O shifts. The molecular orbitals of the group 4 OMF(2) molecules show triple bond character for the terminal oxo groups, which are also supported by an NBO analysis. These molecular orbitals include a σ bond (O(2p) + Ti(sd hybrid)), a normal electron pair π bond (O(2p) + Ti(d)), and a dative π bond arising from O lone pair donation to the overlapping Ti d orbital. The M-O bond dissociation energies for OMF(2) are comparable to those in the diatomic oxide molecules. The OTiF intermediate is also observed through two slightly lower frequency bond stretching modes, and its yield is increased in complementary TiO + F(2) experiments. Finally, the formation of group 4 OMF(2) molecules is highly exothermic due to the weak O-F bonds in OF(2) as well as the strong new MO and M-F bonds formed.

The scandium oxydifluoride free radical, OScF(2), is produced by the spontaneous, specific reaction of laser ablated Sc atoms with OF(2) in solid argon and characterized by using matrix infrared spectroscopy and theoretical calculations. The OScF(2) molecule is predicted to have C(2v) symmetry and a (2)B(2) ground state with an unpaired electron located primarily on the terminal oxygen atom, which makes it a scandium difluoride molecule coordinated by a neutral oxygen atom radical in forming the Sc-O single bond. The closed shell singlet OScF molecule with an obtuse bent geometry has a much shorter Sc-O bond of 1.682 Å than that of the OScF(2) radical (1.938 Å) on the basis of B3LYP calculations. The Sc-O bond in OScF consists of two covalent bonds and a dative bond in which the oxygen 2p(π) lone pair donates electron density into an empty Sc 3d orbital thus forming a triple oxo bond. Density functional calculations suggest it is highly exothermic for fluorine transfer from OF(2) to scandium, which favors the formation of the OScF(2) radical species as well as the OScF molecule after fluorine loss., (Copyright © 2012 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

The electronic and geometrical structures of the lowest triplet states of (GaAs) n clusters ( n = 2-16) are studied using density functional theory with generalized gradient approximation (DFT-GGA). It is found that the triplet-state geometries are different from the corresponding singlet-state geometries; for n = 2-8, 10, and 11, the triplets and singlets have different topologies, while the (GaAs) 9, (GaAs) 12, (GaAs) 15, and (GaAs) 16 triplets possess a reduced symmetry, due to Jahn-Teller distortions. Except for GaAs, the singlet states are the ground states. Excitation energies and oscillator strengths are computed for excitations from the ground state to ten singlet states of all (GaAs) n clusters using time-dependent density functional theory. The adiabatic singlet-triplet gap is compared to the vertical gap, and the difference in the eigenvalues of the highest-occupied and lowest-unoccupied molecular orbitals (the HOMO-LUMO gap). While these three values show large oscillations for small n, they approach each other as the cluster size grows. Thus, the HOMO-LUMO gap computed using the DFT-GGA approach presents a rather reliable estimate of the adiabatic singlet-triplet gap.

We have computed the electronic and geometrical structures of thirteen atom manganese clusters in all three charge states, Mn(13) (-), Mn(13) (+), and Mn(13) by using density functional theory with the generalized gradient approximation. Our results for Mn(13) (-) are compared with our anion photoelectron spectrum of Mn(13) (-), published in this paper. Our results for Mn(13) (+) are compared with the previously published photoionization results of Knickelbein [J. Chem. Phys. 106, 9810 (1997)]. There is a good agreement between theoretical and experimental values of ionization and electron attachment energies.

Reactions of laser-ablated group 4 transition-metal atoms with acetylene have been carried out. The ethynyl metal hydrides (HC[triple bond]C-MH) and corresponding pi complexes (M-eta(2)-(C2H2)) are identified in the matrix infrared spectra. The observed M-H and C-M stretching absorptions show that oxidative C-H insertion readily occurs during codeposition and photolysis afterward. The absorptions from the pi complex, on the other hand, are relatively weak in the original deposition spectrum but increase dramatically in the process of annealing. The vinylidene complex, another plausible product, is not identified in this study. The observed spectra and DFT calculations both show that the back-donations from the group 4 metals to the antibonding pi* orbital of C2H2 are extensive such that the group 4 metals form unusually strong pi complexes. Thus, it is the formation of two Ti-C bonds in the group 4 systems than leads to the stronger bonding than that in the group 8 systems. While bonds form, the Ti atom is weakly bound to C2H2, and we still refer to it as a pi complex. Evidence of relativistic effects is also observed in frequency trends for the Ti, Zr, and Hf products.

Electronic and geometrical structures of neutral, negatively, and positively charged (GaAs)n clusters are computed using density functional theory with generalized gradient approximation. All-electron computations are performed on (GaAs)2-(GaAs)9 while effective core potentials (ECPs) are used for (GaAs)9-(GaAs)15. Calibration calculations on GaAs and (GaAs)9 species support the use of the ECP for the larger clusters. The ground-state geometries of (GaAs)n(-) and/or (GaAs)n+ are different from the corresponding neutral ground-state geometry, except for n=7, 9, 12, 14, and 15, where the neutral and ions have similar structures. Beginning with n=6, all atoms are three coordinate, except for (GaAs)10+ and (GaAs)13+. For the larger species, there is a competition between fullerenes built from hexagons and rhombi and geometrical configurations where Ga-Ga and As-As bonds are formed, which results in the formation of pentagons. As expected, the static polarizability varies in the order of anion>neutral>cation, but the values are rather similar for all three charge states. The thermodynamic stability for the loss of GaAs is reported.

Electronic and geometrical structures of iron clusters with associative (FeNO, Fe2NO, Fe3NO, Fe4NO, Fe5NO, and Fe6NO) and dissociative (OFeN, OFe2N, OFe3N, OFe4N, OFe5N, and OFe6N) attachments of NO, as well as the corresponding singly negatively and positively charged ions, are computed using density functional theory with generalized gradient corrections. Both types of isomers are found to be stable and no spontaneous dissociation was observed during the geometry optimizations. The ground states correspond to dissociative attachment of NO for all iron clusters Fe(n), except for Fe and Fe+. All of the OFe(n)N clusters have ferrimagnetic ground states, except for OFe2N, OFe2N-, OFe4N, and OFe4N-, which prefer the ferromagnetic coupling. In the ferrimagnetic states, the excess spin density at one iron atom couples antiferromagnetically to the excess spin densities of all other iron atoms. Relative to the high-spin Fe(n) ground state, the lowest energy ferrimagnetic state quenches the total magnetic moments of iron clusters by 7, which is to be compared with a reduction in the magnetic moment of one in the lowest energy ferromagnetic states. Dissociation of NO on the iron clusters has a pronounced impact on the energetics of reactions; the Fe(n)NO+CO-->Fe(n)N+CO2 channels are exothermic while the OFe6N+CO--> Fe6N+CO2 channels are nearly thermoneutral.

Electronic and geometrical structures of Mn(3)-Mn(10) together with their singly negatively and positively charged ions are computed using density functional theory with generalized gradient approximation. The ground-state spin multiplicities in the neutral series are 16, 21, 4, 9, 6, 5, 2, and 5, for Mn(3)-Mn(10), respectively. Thus, there is a transition from a ferromagnetic ground state to a ferrimagnetic ground state at Mn(5). The energy difference between ferrimagnetic and ferromagnetic states in Mn(n) grows rapidly with increasing n and exceeds 2 eV in Mn(10). The corresponding change from ferro- to ferrimagnetic ground state occurs at Mn(6)(-) and Mn(3)(+) in the anionic and cationic series, respectively. Beginning with Mn(6), the ion spin multiplicities differ from that of the neutral by +/-1 (i.e., they obey the empirical "+/-1 rule"). We found that the energy required to remove an Mn atom is nearly independent of the charge state of an Mn(n) cluster and the number of atoms in the cluster, except for Mn(3). The results of our calculations are in reasonable agreement with experiment, except for the experimental data on the magnetic moments per atom, where, in general, we predict smaller values than the experiment.

Electronic and geometrical structures of Fe4Cn(CO)m (n + m < or = 6) and their singly negatively and positively charged ions are computed using density functional theory with generalized gradient approximation. Isomers with CO bonded directly to the cluster iron atoms and bonded to a carbon atom chemisorbed on the cluster surface are optimized for the Fe4C2CO, Fe4C2(CO)2, Fe4C3CO, and Fe4C4CO series. Optimizations of a large number of differently shaped Fe4C4, Fe4C5, and Fe4C6 clusters are performed to find trends in preferable arrangement of carbon atoms, in particular, to determine the relative energetics of structures with single C atoms versus those with C2 dimers or C3 trimers. The computed total energies are used to estimate the energetics of the Boudouard disproportionation reactions Fe4Cn(CO)m + 00-->Fe4Cn+1(CO)m-1 + CO2.

Rotational potentials are computed for heptamers and nonamers of an amino, nitro-substituted phenylene ethynylene trimer molecule. A herringbone and a parallel-slipped packing arrangement are considered. The effect of tilting the molecules with respect to the surface as well as the effect of the gold support are also taken into account. The herringbone structure with the molecules perpendicular to the surface has a low rotational barrier (2 kcal/mol). Tilting the molecules by 30 degrees increases the rotational barriers significantly (16 kcal/mol). The parallel-slipped structure has rotational barriers of 7 kcal/mol. Including the effect of the gold support increases the rotational barriers for tilted molecules but has little effect on perpendicular molecules.

The ground states of the mixed 3d-metal dimers TiV, TiCr, TiMn, TiFe, TiCo, TiNi, TiCu, TiZn, VCr, VMn, VFe, VCo, VNi, VCu, VZn, CrMn, CrFe, CrCo, CrNi, CrCu, CrZn, MnFe, MnCo, MnNi, MnCu, MnZn, FeCo, FeNi, FeCu, FeZn, CoNi, CoCu, CoZn, NiCu, NiZn, and CuZn along with their singly negatively and positively charged ions are assigned based on the results of computations using density functional theory with generalized gradient approximation for the exchange-correlation functional. Except for TiCo and CrMn, our assignment agrees with experiment. Computed spectroscopic constants (r(e),omega(e),D(o)) are in fair agreement with experiment. The ground-state spin multiplicities of all the ions are found to differ from the spin multiplicities of the corresponding neutral parents by +/-1. Except for TiV, MnFe, and MnCu, the number of unpaired electrons, N, in a neutral ground-state dimer is either N(1)+N(2) or mid R:N(1)-N(2)mid R:, where N(1) and N(2) are the numbers of unpaired 3d electrons in the 3d(n)4s(1) occupation of the constituent atoms. Combining the present and previous results obtained at the same level of theory for homonuclear 3d-metal and ScX (X=Ti-Zn) dimers allows one to construct "periodic" tables of all 3d-metal dimers along with their singly charged ions.

Photoelectron spectroscopy (PES) is combined with density functional theory (DFT) to study the monochromium carbide clusters CrCn- and CrCn (n = 2-8). Well-resolved PES spectra were obtained, yielding structural, electronic, and vibrational information about both the anionic and neutral clusters. Experimental evidence was observed for the coexistence of two isomers for CrC2-, CrC3-, CrC4-, and CrC6-. Sharp and well-resolved PES spectra were observed for CrCn- (n = 4,6,8), whereas broad spectra were observed for CrC5- and CrC7-. Extensive DFT calculations using the generalized gradient approximation were carried out for the ground and low-lying excited states of all the CrCn- and CrCn species, as well as coupled-cluster calculations for CrC2- and CrC2. Theoretical electron affinities and vertical detachment energies were calculated and compared with the experimental data to help the assignment of the ground states and obtain structural information. We found that CrC2- and CrC3- each possess a close-lying cyclic and linear structure, which were both populated experimentally. For the larger CrCn- clusters with n = 4, 6, 8, linear structures are the overwhelming favorite, giving rise to the sharp PES spectral features. CrC7- was found to have a cyclic structure. The broad PES spectra of CrC5- suggested a cyclic structure, whereas the DFT results predicted a linear one.

Figure 18 presents the C-O stretching vibrational frequencies of the first-row transition-metal monocarbonyl cations, neutrals, and anions in solid neon; similar diagrams have been reported for neutral MCO species in solid argon, but three of the early assignments have been changed by recent work and one new assignment added. The laser-ablation method produces mostly neutral atoms with a few percent cations and electrons for capture to make anions; in contrast, thermal evaporation gives only neutral species. Hence, the very recent neon matrix investigations in our laboratory provide carbonyl cations and anions for comparison to neutrals on a level playing field. Several trends are very interesting. First, for all metals, the C-O stretching frequencies follow the order cations > neutrals > anions with large diagnostic 100-200 cm-1 separations, which is consistent with the magnitude of the metal d to CO pi * donation. Second, for a given charge, there is a general increase in C-O stretching vibrational frequencies with increasing metal atomic number, which demonstrates the expected decrease in the metal to CO pi * donation with increasing metal ionization potential. Some of the structure in this plot arises from the extra stability of the filled and half-filled d shell and from the electron pairing that occurs at the middle of the TM row; the plot resembles the "double-humped" graph found for the variation in properties across a row of transition metals. For the anions, the variation with metal atom is the smallest since all of the metals can easily donate charge to the CO ligand. Third, for the early transition-metal Ti, V, and Cr families, the C-O stretching frequencies decrease when going down the family, but the reverse relationship is observed for the late transition-metal Fe, Co, and Ni families. In most of the present discussion, we have referred to neon matrix frequencies; however, the argon matrix frequencies are complementary, and useful information can be obtained from comparison of the two matrix hosts. In most cases, the neon-to-argon red shift for neutral carbonyls is from 11 to 26 cm-1, but a few (CrCO) lie outside of this range. In the case of FeCO and Fe(CO)2, it appears that neon and argon trap different low-lying electronic states. In general, the carbonyl neutrals and anions have similar shifts but carbonyl cations have larger matrix shifts. For example, the FeCO+ fundamental is at 2123.0 cm-1 in neon and 2081.5 cm-1 in argon, a 42.5 cm-1 shift, which is larger than those found for FeCO- (11.7 cm-1) and FeCO (11.7 cm-1). It is unusual for different low-lying electronic states to be trapped in different matrices, but CUO provides another example. The linear singlet state (1047.3, 872.2 cm-1) is trapped in solid neon, and a calculated 1.2 kcal/mol higher triplet state is trapped in solid argon (852.5, 804.3 cm-1) and stabilized by a specific interaction with argon. The bonding trends are well described by theoretical calculations of vibrational frequencies. Table 5 compares the scale factors (observed neon matrix/calculated) for the C-O stretching modes of the monocarbonyl cations, neutrals, and anions of the first-row transition metals observed in a neon matrix using the B3LYP and BP86 density functionals. Most of the calculated carbonyl harmonic stretching frequencies are within 1% of the experimental fundamentals at the BP86 level of theory, while calculations using the B3LYP functional give frequencies that are 3-4% higher as expected for these density functionals and calculations on saturated TM-carbonyls. For second- and third-row carbonyls using the BP86 density functional and the LANL effective core potential in conjunction with the DZ basis set, the agreement between theory and experiment is just as good. For example, the 16 M(CO)1-4 neutral and anion and 2 MCO+ cation (M = Ru, Os) carbonyl frequencies are fit within 1.5%. The 16 species (M = Rh, Ir) are fit within 1%, but the Rh(CO)1-4+ calculations are 2-3% too low and Ir(CO)1-4+ computations are 1-2% too low. In addition to predicting the vibrational frequencies, DFT can be used to calculate different isotopic frequencies, and isotopic frequency ratios can be computed as a measure of the normal vibrational mode in the molecule for an additional diagnostic. For diatomic CO, the 12CO/13CO ratio 1.0225 and C16O/C18O ratio 1.0244 characterize a pure C-O stretching mode. In a series of molecules such as RhCO+, RhCO, and RhCO-, where the metal-CO bonding varies, the Rh-C, C-O vibrational interaction is different and the unique isotopic ratios for the carbonyl vibration are characteristic of that particular molecule. Table 6 summarizes the isotopic ratios observed and calculated for the RhCO+,0,- species. Note that RhCO+ exhibits slightly more carbon-13 and less oxygen-18 involvement in the C-O vibration than CO itself and that this trend increases to RhCO and to RhCO- as the Rh-C bond becomes shorter and stronger. Note also how closely the calculated and observed ratios both follow this trend. In a molecule with two C-O stretching modes, for example, bent Ni(CO)2 exhibits a strong b2 mode at 1978.9 cm-1 and a weak a1 mode at 2089.7 cm-1 in solid neon, and these two modes involve different C and O participations. The symmetric mode shows substantially more C (1.0242) and less O (1.0217) participation than does the antisymmetric mode with C (1.0228) and O (1.0238) involvement, based on the given isotopic frequency ratios, which are nicely matched by DFT calculations (a1 1.0244, 1.0224 and b2 1.0232, 1.0241, respectively). These investigations of vibrational frequencies in unsaturated transition-metal carbonyl cations, neutrals, and anions clearly demonstrate the value of a close working relationship between experiment and theory to identify and characterize new molecular species.

Density functional theory has been employed to calculate the harmonic frequencies and intensities of a range of polycyclic aromatic hydrocarbon (PAH) cations that explore both size and electronic structure effects on the infrared spectroscopic properties of these species. The sample extends the size range of PAH species considered to more than 50 carbon atoms and includes several representatives from each of two heretofore unexplored categories of PAH cations: (1) fully benzenoid PAH cations whose carbon skeleton is composed of an odd number of carbon atoms (C(odd) PAHs); and (2) protonated PAH cations (HPAH+). Unlike the radical electronic structures of the PAH cations that have been the subject of previous theoretical and experimental work, the species in these two classes have a 'closed'-shell electronic configuration. The calculated spectra of circumcoronene, C54H18, in both neutral and (radical) cationic form are also reported and compared with those of the other species. Overall, the C(odd) PAHs spectra are dominated by strong CC stretching modes near 1600 cm(-1) and display spectra that are remarkably insensitive to molecular size. The HPAH+ species evince a more complex spectrum consistent with the added contributions of aliphatic modes and their generally lower symmetry. Finally, for both classes of closed-shell cations, the intensity of the aromatic CH stretching modes is found to increase with molecular size far out of proportion with the number of CH groups, approaching a value more typical of neutral PAHs for the largest species studied.

Thanks to the mid-IR sensitivities of the ISO and IRTS orbiting spectrometers it is now possible to search the diffuse interstellar medium for heretofore inaccessible molecular emission. In view of the recent strong case for the presence of C(7-) (Kirkwood et al. 1998, Tulej et al. 1998),and the fact that carbon chains possess prominent infrared active modes in a very clean portion of the interstellar spectrum, we have analyzed the IRTS spectrum of the diffuse interstellar medium for the infrared signatures of these species. Theoretical and experimental infrared band frequencies and absolute intensities of many different carbon chain species are presented. These include cyanopolyynes, neutral and anionic linear carbon molecules, and neutral and ionized, even-numbered, hydrogenated carbon chains. We show that--as a family--these species have abundances in the diffuse ISM on the order of 10(-10) with respect to hydrogen, values consistent with their abundances in dense molecular clouds. Assuming an average length of 10 C atoms per C-chain implies that roughly a millionth of the cosmically available carbon is in the form of carbon chains and that carbon chains can account for a few percent of the visible to near-IR diffuse interstellar band (DIB) total equivalent width (not DIB number).

Calculations are carried out using density functional theory (DFT) to determine the harmonic frequencies and intensities of 1-methylanthracene, 9-methylanthracene, 9-cyanoanthracene, 2-aminoanthracene, acridine, and their positive ions. The theoretical data are compared with matrix-isolation spectra for these species also reported in this work. The theoretical and experimental frequencies and relative intensities for the neutral species are in generally good agreement, whereas the positive ion spectra are only in qualitative agreement. Relative to anthracene, we find that substitution of a methyl or CN for a hydrogen does not significantly affect the spectrum other than to add the characteristic methyl C-H and C triple bond N stretches near 2900 and 2200 cm-1, respectively. However, addition of NH2 dramatically affects the spectrum of the neutral. Not only are the NH2 modes themselves strong, but this electron-withdrawing group induces sufficient partial charge on the ring to give the neutral molecule spectra characteristics of the anthracene cation. The sum of the absolute intensities is about four times larger for 2-aminoanthracene than those for 9-cyanoanthracene. Substituting nitrogen in the ring at the nine position (acridine) does not greatly alter the spectrum compared with anthracene.

Full configuraton-interaction (FCI) calculations have given an unambiguous standard by which the accuracy of theoretical approaches of incorporating electron correlation into molecular structure calculations can be judged. In addition, improvements in vectorization of programs, computer technology, and algorithms now permit a systematic study of the convergence of the atomic orbital (or so-called one-particle) basis set. These advances are discussed and some examples of the solution of chemical problems by quantum mechanical calculations are given to illustrate the accuracy of current techniques.